微分方程數(shù)值解實(shí)驗(yàn)-文庫吧
2025-04-01 23:19 本頁面
【正文】 38/3*z3。 %eval_r(abz)。 ksx=10*(y2x2)。 ksy=x2*(28z2)y2。ksz=x2*y28/3*z2。 ktx=10*(y1x1 )。 kty=x1*(28z1)y1。ktz=x1*y18/3*z1。 x=x3+h/12*(23*kfx16*ksx+5*ktx)。y=y3+h/12*(23*kfy16*ksy+5*kty)。z=z3+h/12*(23*kfz16*ksz+5*ktz)。end[t,y]=ode45(Rk4,[0,30],[12,2,9])suptitle(39。RungeKutta4階法39。) %總標(biāo)題subplot(2,2,1)。plot(t,y(:,1))。grid on。legend(39。x關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。x39。,39。FontSize39。,14)。subplot(2,2,2)。plot(t,y(:,2))。grid on。legend(39。y關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。y39。,39。FontSize39。,14)。subplot(2,2,3)plot(t,y(:,3))。grid on。legend(39。z關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。z39。,39。FontSize39。,14)。subplot(2,2,4)plot3(y(:,1),y(:,2),y(:,3))。grid on。legend(39。x,y,z的空間關(guān)系圖39。,1)。xlabel(39。x39。,39。FontSize39。,14)。ylabel(39。y39。,39。FontSize39。,14)。zlabel(39。y39。,39。FontSize39。,14)。view(40,60)。 %鎖定同樣的視圖，便于比較a=1:1:3000。b=1:1:3000。c=1:1:3000。t=1:1:3000。x1=5。 y1=5。 z1=10。 x2=。 y2=。 z2=。x3=。 y3=。 z3=。for i=1:1:3000 [x,y,z]=adams4(x1,y1,z1,x2,y2,z2,x3,y3,z3,)。a(i)=x。b(i)=y。c(i)=z。x1=x2。y1=y2。z1=z2。x2=x3。y2=y3。z2=z3。x3=x。y3=y。z3=z。fprintf(39。x=%f39。39。y=%f39。39。z=%f\n39。,x,y,z)。%顯示迭代值end suptitle(39。Adams4階外插法39。) %總標(biāo)題subplot(2,2,1)。plot(t,a)。grid on。 legend(39。x關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。x39。,39。FontSize39。,14)。subplot(2,2,2)。plot(t,b)。grid on。 legend(39。y關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。y39。,39。FontSize39。,14)。subplot(2,2,3)。plot(t,c)。grid on。 legend(39。z關(guān)于t 的變化關(guān)系圖39。,1)。xlabel(39。t39。,39。FontSize39。,14)。ylabel(39。z39。,39。FontSize39。,14)。subplot(2,2,4)plot3(a,b,c)。grid on。 legend(39。x,y,z的空間關(guān)系圖39。,1)。xlabel(
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1