www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著,每章的MATLAB源代码程序 > 第21章/代码21-3.txt
clear % 清理工作空间 p=[0 0;0 15]; b=pi/8;e=1.5; a=[cos(b) -sin(b);sin(b) cos(b)]; n=2; for k=1:4 i=1; for j=1:2:n p1=p(j,:); p2=p(j+1,:); d=(p2-p1)/3; new(i,:)=p1;i=i+1; q1=p1+d; new(i,:)=q1;i=i+1; new(i,:)=q1;i=i+1; new(i,:)=p1+d+e*d*a';i=i+1; new(i,:)=q1;i=i+1; new(i,:)=p1+d+e*d*a;i=i+1; new(i,:)=p1+d;i=i+1; q2=p1+2*d; new(i,:)=q2;i=i+1; new(i,:)=q2;i=i+1; new(i,:)=p1+2*d+d*a;i=i+1; new(i,:)=q2;i=i+1; new(i,:)=p1+2*d+d*a';i=i+1; new(i,:)=p1+2*d;i=i+1; new(i,:)=p2;i=i+1; end n=n*7;p=new; end for m=1:2:n line([p(m,1) p(m+1,1)],[p(m,2) p(m+1,2)],'Color','k'); end axis equal % 设置显示比例 axis off % 不显示坐标轴