代码21-3.txt

matlab宝典书籍源码 全部代码 希望大家喜欢

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                                                  % 不显示坐标轴