a=6;
b=3;
t=0:0.01:1;
x=a*cos(degtorad(t*360));
y=b*sin(degtorad(t*(2*360)));
v=[x;y];
x_center = 1;
y_center = 3;
center = repmat([x_center; y_center], 1, length(x));
set(gca,'color',[0,0,0]);
axis([-8,10,-6,10])
for i=-pi+1:0.1:pi
theta = i;
R = [cos(theta) -sin(theta); sin(theta) cos(theta)];
vo = R*(v - center) + center;
x_rotated = vo(1,:);
y_rotated = vo(2,:);
hold on
plot(x_rotated, y_rotated,'g');
pause(0.1)
end
set(gca,'color',[0,0,0]);
%%%偏差部分
% dif = 0:10:360;
% r = 2;
% dif_X = r*cos(degtorad(dif));
% dif_Y = r*sin(degtorad(dif));
% for i=1:length(dif)
% hold on
% plot(x+dif_X(i),y+dif_Y(i))
% end