%Made by huxuec
clear
R=2;v=0.5;t=60;%半徑,速度,時間
a=-pi/2;%m點的起始位置
w=v/R;%角速度
tt=0:0.1:t;
[nn,n]=size(tt)
b=(0:20)*pi/10
cx=v.*tt;
cy=ones(1,n).*R; %圓心坐標
rx=[]
ry=[]
for i=1:n
rx(i,:)=cx(i)+R.*cos(b);
ry(i,:)=cy(i)+R.*sin(b); %圓上點的坐標
end
for i=1:n
mx(i)=cx(i)+R.*cos(a-w.*tt(i));
my(i)=cy(i)+R.*sin(a-w.*tt(i)); %m坐標
%---------圓圈----------------m軌跡---------------m點.--------
plot(rx(i,:),ry(i,:),'b',mx(1:i),my(1:i),'r',mx(i),my(i),'.')
axis([-R v*t+R 0 (v*t+2*R)]) %規定視窗大小
drawnow
end