matlab绘制peano（皮亚诺）曲线和koch分形曲线

飞霜墨染褚

matlab绘制peano（皮亚诺）曲线和koch（科赫曲线，雪花曲线）分形曲线

koch曲线
matlab plot函数绘制koch曲线程序，程序还是比较简单的，这里只绘制出了雪花的三分之一

1. function koch_curve(number)%number代表koch的阶数，范围为大于等于2  
   
2. figure  
   
3. set(gcf,'position',[0,0,1920,1080]);%设置窗口分辨率，[0,0]和[1920,1080]分别为窗口左上角和右下角坐标可根据自己的屏幕分辨率调整，注释掉这句则使用matlab默认窗口分辨率  
   
4. n=2;  
   
5. koch1=[0,0;1,0];  
   
6.     for i=1:number  
   
7.         koch2=zeros(4*n-3,2);  
   
8.         k=2;  
   
9.         for j=2:n  
   
10.             koch2(k,:)=[(koch1(j-1,1)*2+koch1(j,1))/3,(koch1(j-1,2)*2+koch1(j,2))/3];  
    
11.             koch2(k+1,:)=[(koch1(j-1,1)+koch1(j,1)+sqrt(3)*(koch1(j-1,2)-koch1(j,2))/3)/2,(koch1(j-1,2)+koch1(j,2)-sqrt(3)*(koch1(j-1,1)-koch1(j,1))/3)/2];  
    
12.             koch2(k+2,:)=[(koch1(j,1)*2+koch1(j-1,1))/3,(koch1(j,2)*2+koch1(j-1,2))/3];  
    
13.             koch2(k+3,:)=koch1(j,:);  
    
14.             k=k+4;  
    
15.         end  
    
16.         n=4*n-3;  
    
17.         x=koch2(:,1);  
    
18.         y=koch2(:,2);  
    
19.         plot(x,y)  
    
20.         axis equal   
    
21.         koch1=koch2;  
    
22.         pause(1);  
    
23.     end   
    
24. end  

复制代码

matlab 绘制皮亚诺曲线

1. function peano_curve(n)  
   
2.     peano_old=[0,0;0,1;0.5,1;0.5,0;1,0;1,1];  
   
3.     x=peano_old(:,1);  
   
4.     y=peano_old(:,2);  
   
5.     figure  
   
6.     set(gcf,'position',[0,0,1920,1080]);  
   
7.     plot(x,y)  
   
8.     axis equal  
   
9.     for i=1:n-1  
   
10.         p1=[peano_old(:,1),2+1/(3^i-1)-peano_old(:,2)];  
    
11.         p1=p1(length(p1):-1:1,:);  
    
12.         p2=[p1(:,1),4+3/(3^i-1)-p1(:,2)];  
    
13.         p2=p2(length(p2):-1:1,:);  
    
14.         peano_new=[peano_old;p1;p2];  
    
15.         p1=[2+1/(3^i-1)-peano_new(:,1),peano_new(:,2)];  
    
16.         p1=p1(length(p1):-1:1,:);  
    
17.         p2=[4+3/(3^i-1)-p1(:,1),p1(:,2)];  
    
18.         p2=p2(length(p2):-1:1,:);  
    
19.         peano_new=[peano_new;p1;p2];  
    
20.         peano_old=peano_new/(3+2/(3^i-1));  
    
21.         x=peano_old(:,1);  
    
22.         y=peano_old(:,2);  
    
23.         plot(x,y)  
    
24.         axis equal  
    
25.         pause(1)  
    
26.     end  
    
27. end

复制代码