exm060222_1.m

好多m文件

%平面上的四个点和它们构成的方块
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