www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著，每章的MATLAB源代码程序 > 第21章/代码21-1.txt

    
p=[0  0;10  0];
n=2;
A=[cos(pi/3)  -sin(pi/3);sin(pi/3)  cos(pi/3)];
for  k=1:5
    d=diff(p)/3;m=4*n-3;
    q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);  
        p(2:4:m,:)=q+d;                 
    p(3:4:m,:)=q+d+d*A'; 
    p(4:4:m,:)=q+2*d;
        n=m; 
end
plot(p(:,1),p(:,2),'k')                    % 画图
axis equal                                 % 设置显示比例
axis off                                   % 不显示坐标轴