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