ycoarse.m

请不要用做非法用途

g=9.8;
PI=3.1415926;
k=0.2*g;
t=0.01;
dphi=PI/180;
 x=[70*dphi/60
    80*dphi/60
    10*dphi];
xx=[x(1)
   x(2)];
for i=1:1000   
tt=[cos(xx(2))*cos(x(3))-sin(xx(1))*sin(xx(2))*sin(x(3)) cos(xx(2))*sin(x(3))+sin(xx(1))*sin(xx(2))*cos(x(3))
   -cos(xx(1))*sin(x(3)) cos(xx(1))*cos(x(3))];
dx=tt*(-k)*xx;
xx=dx*t+xx;
data(i)=xx(1)*57.3*60;
data1(i)=xx(2)*57.3*60;
end;
m=1:1000;
n=m/10;
plot(n,data(m),n,data1(m));
ylabel('');
xlabel('t/s');