资料试验微分方程数值解.pptxVIP

会计学;例 用向前欧拉公式解方程;;y=adveula(ode121,0,1,1,10) y,= ,,,,0,,,,,,,,,,,,,1.0000 ,,,,0.1000,,,,1.0959 ,,,,0.2000,,,,1.1841 ,,,,0.3000,,,,1.2662 ,,,,0.4000,,,,1.3434 ,,,,0.5000,,,,1.4164 ,,,,0.6000,,,,1.4860 ,,,,0.7000,,,,1.5525 ,,,,0.8000,,,,1.6165 ,,,,0.9000,,,,1.6782 ,,,,1.0000,,,,1.7379;改进的欧拉公式可以从几何上加以解释如下：;欧拉方法可以推广到解高阶方程和方程组，如：;,注意到;,,,,,,,,满足l1+,l2=1,,,2a,l2=1,的公式;3,,四阶龙格-库塔方法的MATLAB实现;§2 应用问题求解;设走私船在t时刻的位置为x1(t),y1(t),则;,例2、食饵-捕食者模型（弱肉强食）;-d;ts=0:0.1:15; x0=[25,2]; [t,x]=ode45(shier,ts,x0) result=[t,x] plot(t,x),grid,gtext(x1(t)),gtext(x2(t)), figure; plot(x(:,1),x(:,2)),grid,... ,,,xlabel(x1),ylabel(x2);食饵-捕食者模型的数值解;从前面所得的方程组初值问题;上机内容;4,,,火箭发射;第一阶段，燃料燃烧推动阶段，为变加速运动，初值问题为;function,dx=odehj1(t,x) global,M,r,N,k g=9.8; dx=zeros(2,1); dx(1)=x(2); dx(2)=(N-(M-r*t)*g-k*x(2)^2)/(M-r*t);;clear global,M,Q,r,N,k,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%声明全局变量 M=1400;Q=1080;N=32000;r=18;k=0.4;T=Q/r;x0=[0,0]; [t1,x]=ode45(odehj1,[0:T],x0);,,,,,,,,,,,,,,,,,,,,,,,%求第一阶段解 h1=x(end,1);v1=x(end,2);,,,,,,,,,%引擎关闭时火箭高度和速度 y0=[h1,v1];,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%第二阶段初值 [t2,y]=ode45(odehj2,[T:T+,20],y0);,,,,,,,,,,,,,,,,,,,,,,%求第二阶段解 s=max(find(y(:,2)=0));,,,,,,,,,,,,,,,,,,,,%速度为零时的时间分量 hm=y(s,1);,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%最大高度 t=[t1;t2(2:end)];h=[x(:,1);y(2:end,1)];,,,,,,,,,,,%高度与时间 v=[x(:,2);y(2:end,2)];,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%速度函数;L=length(h);dt=(T+12)/(L-1);,,,,,,,,,,,,,,,,,,%用三点公式求加速度 ac(1)=(-3*v(1)+4*v(2)-v(3))/2/dt;,,,,,,,,,,,,,,,,,,,,,,,%起始点加速度 for,k=2:L-1,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, ,,,ac(k)=(v(k+1)-v(k-1))/2/dt;,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%中间点加速度 end, ac(L)=(v(L-2)-4*v(L-1)+3*v(L))/2/dt;,,,,,,,,,,,,,,,,%最后点加速度 a1=ac(length(t1)-1);,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%关闭引擎时加速度 am=ac(length(t1)+s-1);,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,%最高??加速度 result1=[T,h1,v1,a1],,,,,,,,,,,,,,,,%引擎关闭时高度,速度和加速度 result2=[T+t2(s),hm,am],,,,,,,,,,%达到最高点时间高度和加速度 subplot(1,3,1),plot(t,h,b-),axis([0,80,0,14