%% 旋转
x=[0 1 1 0 0];
y=[0 0 1 1 0];
plot(x,y)
h=patch(x,y,'b');
h.FaceAlpha=0.5;
hold on
axis([-2 2 -2 2])
axis equal
for i=1:40 %旋转40次
M=[cos(pi/40*i) -sin(pi/40*i);
sin(pi/40*i) cos(pi/40*i)];
T=M*[x;y];
plot(T(1,:),T(2,:))
h2=patch(T(1,:),T(2,:),'r');
h2.FaceAlpha=0.1;
pause(0.2)
end
%% X轴压缩
x=[0 1 1 0 0];
y=[0 0 1 1 0];
plot(x,y)
h=patch(x,y,'b');
h.FaceAlpha=0.5;
hold on
axis([-2 2 -2 2])
axis equal
M=[0.2 0;0 1];
T=M*[x;y];
plot(T(1,:),T(2,:))
h2=patch(T(1,:),T(2,:),'r');
h2.FaceAlpha=0.1;