matlab 画x a的直线方程式,用matlab以动画的方式绘制出摆线 x=a(t-sint) y=a(1-cost) （a自己赋值）的渐屈线...

clear;

clc;

close;

a=1;

syms t

x=a*(t-sin(t));

y=a*(1-cos(t));

ezplot(x,y,[0,2*pi]),grid on;hold on;

dy=diff(y)/diff(x);

dyy=diff(dy)/diff(x);

xx=x-(1+dy^2)*dy/dyy;%渐屈线的坐标

yy=y+(1+dy^2)/dyy;

M=50;

t=0;

xxx=subs(xx);

yyy=subs(yy);

H1=plot(xxx,yyy,'r');hold on;grid on;axis([0,7,-2.5,2.5]);

x1=subs(x);

y1=subs(y);

H2=plot([x1,xxx],[y1,yyy],'k--');

H3=plot(x1,y1,'ko');

H4=plot(xxx,yyy,'ro');

tt=linspace(0,2*pi,M);

for i=1:M

pause(0.2);

t=tt(1:i);

xxx=subs(xx);

yyy=subs(yy);

x1=subs(x);

y1=subs(y);

set(H1,'xdata',xxx,'ydata',yyy);

<