%平面上的四個點和它們構成的方塊
p1=[-0.5,0,1]';p2=[-0.5,1,1]';p3=[0.5,1,1]';p4=[0.5,0,1]';
Sq=[p1,p2,p3,p4,p1];
%平移投影：沿x軸移動0.5 , 沿y軸移動1 。
dx=0.5;dy=1;T=[1,0,dx;0,1,dy;0,0,1];
%旋轉投影：逆時針旋轉30度。
th=pi/6;R=[cos(th),-sin(th),0;sin(th),cos(th),0;0,0,1];
%刻度投影：x方向放大到2倍，y方向放大到3倍。
alpha=2;beta=3;S=[alpha,0,0;0,beta,0;0,0,1];
E=eye(3,3);%為編程方便，設計一個單位陣。在它投影下，任何形狀都不變。	<10>
TRS={E,T,R,S};				%用元胞數組存放四個變換矩陣				<11>
ss={'r^','rd','rp','rh'};	%用元胞數組存放數據點形設置值				<12>
tt={'Original Square','Translation','Rotation','Scaling'};		
								%用元胞數組存放四張子圖的圖名				<13>
for i=1:4
	W=TRS{i}*Sq;							%進行仿射變換
	subplot(2,2,i)
	for k=1:4
		plot(W(1,k),W(2,k),ss{k});		%用不同點形標志四個頂點		<19>
		axis([-3,3,-1,5]),axis equal
		hold on							%使以后圖形畫在當前子圖上		<21>
	end
	plot(W(1,:),W(2,:))					%連接四個頂點使圖形封閉		<23>
	grid on								%畫過坐標格線
	title(tt{i})							%給子圖題寫圖名
	hold off								%使以后圖形不再畫在當前子圖上
end