circle.m

来自「用Matlab编写的导弹比例导引运动、圆形运动、三维曲线运动模型」· M 代码 · 共 21 行

clear
r=200;
g=pi/4;
w=pi/100;
T=0.1;
m=-1;
W=[0.5;1;1.5;2;2.5];  
for j=1:5
      m=m*(-1);
      w=m*W(j);        
    for i=1:200 
    x(i)=r*cos(g+w*(i-1)*T);
    y(i)=r*sin(g+w*(i-1)*T);
    a=pi-w*i*T;
    b=w*T;
    if a<b
         break;
    end  
  end
end
plot(x,y)