yycure3.m

来自「这是matlab编程的一个例子」· M 代码 · 共 17 行

x=0:0.0645:0.4515;
x1=0:0.0645:0.1935;
x2=0.1935:0.0645:0.4515;
y1=100*0.02*sin(10.47*x1)+100*0.52*x1.*sin(10.47*x1);
y2=100*0.02*sin(10.47*x2)+10*sin(10.47*x2)-100*0.388*x2.*sin(10.47*x2);
plot(x1,y1,'b*-',x2,y2,'b*-')
theAxes=axis;
fmat=moviein(30);
for j=1:30
    y1=100*0.02*sin(10.47*x1-j*2*pi/30)+100*0.52*x1.*sin(10.47*x1-j*2*pi/30);
    y2=100*0.02*sin(10.47*x2-j*2*pi/30)+17.51*sin(10.47*x2-j*2*pi/30)-100*0.388*x2.*sin(10.47*x2-j*2*pi/30);
    z=0;
    plot(x1,y1,'r*-',x2,y2,'r*-',x,z,'b*-',(x+0.0215),z,'b--',(x+0.043),z,'b--')
    axis([0 0.5 -12 12])
    fmat(:,j)=getframe;
end
movie(fmat,30);