Question

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的图像。

Answer 1

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的问题，那儿高手多的很。