题目

matlab用四阶龙格库塔法解微分方程组；
dy(1)=-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));
dy(2)=(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);
以上述计算结果再算一下表达式；
x=3000*sin(70*2*pi/360)-y(1)*sin(y(2));
z=3000*cos(70*2*pi/360)-12*t-y(1)*cos(y(2));
（每一组x,z对应一组y1,y2）并画出x,z的图像.

提问时间：2020-12-16

答案

function dy=test(t,y)
dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);];
[t,y]=ode45('test',[0.01,1],[1,1])
plot(t,y(:,1),t,y(:,2));
x=3000*sin(70*2*pi/360)-y(1).*sin(y(2));
z=3000*cos(70*2*pi/360)-12*t-y(1).*cos(y(2));
plot(t,x,t,z)
最后的这图像不好.你可再看看.方法是没问题的.
你可到MATLAB中文论坛去请教关于MATLAB的问题,那儿高手多的很.