微分方程組matlab的解析解數(shù)值解-閱讀頁(yè)
2025-01-14 11:38本頁(yè)面
【正文】 y(2)=y(1)*y(3)。取 t0=0， tf=12，輸入命令： [T,Y]=ode45(39。,[0 12],[0 1 1])。39。*39。+39。 dy(1)=y(2)。 2. 取 x0=0， xf=，建立主程序 ： x0=0， xf= [x,y]=ode15s(39。,[x0 xf],[0 0])。) hold on y=0::2。) 結(jié)論 : 導(dǎo)彈大致在（ 1， ）處擊中乙艦To Matlab(ff6)令 y1=y,y2=y1’，將方程（ 3）化為一階微分方程組。 dy(1)=5*(1y(1))/sqrt((1y(1))^2+(ty(2))^2)。 取 t0=0， tf=2，建立主程序 ： [t,y]=ode45(39。,[0 2],[0 0])。 plot(1,Y,39。)， hold on plot(y(:,1),y(:,2),39。) To Matlab(chase2)5. 結(jié)果見圖 1導(dǎo)彈大致在（ 1， ）處擊中乙艦，與前面的結(jié)論一致 .圖 1 圖 2返 回 在 ，按二分法逐步修改 tf，即分別取 tf=1，,…, 直到 tf=，得圖 2.結(jié)論：時(shí)刻 t=，導(dǎo)彈在（ 1， ）處擊中乙艦。 dy(1)=20*(10+20*cos(t)y(1))/sqrt ((10+20*cos(t)y(1))^2+(20+15*sin(t)y(2))^2)。取 t0=0， tf=10，建立主程序 ： t0=0。 [t,y]=ode45(39。,[t0 tf],[0 0])。 X=10+20*cos(T)。 plot(X,Y,39。) hold on plot(y(:,1),y(:,2),39。) 在 ，不斷修改 tf的值 ,分別取 tf=5, , ,…, 至 ,狗剛好追上慢跑者 .To Matlab(chase3)建立 m文件 : function dy=eq4(t,y) dy=zeros(2,1)。 dy(2)=5*(20+15*sin(t)y(2))/sqrt ((10+20*cos(t) y(1))^2+(20+15*sin(t)y(2))^2)。tf=10。eq439。 T=0::2*pi。 Y=20+15*sin(T)。39。*39。 dx(1)=x(1)*(*x(2))。其次，建立主程序 ： [t,x]=ode45(39。,[0 15],[25 2])。39。*39。 可以猜測(cè)： x1（ t）與 x2（ t）都是周期
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1