Koch曲线、Sierpinski三角形、Cantor集的计算机实现.docVIP

Koch曲线、Sierpinski三角形、Cantor集的计算机实现 Koch曲线 function koch(level) xl=zeros(level,1); xr= xl; xr(level)=1; yl=xl; yr=yl; r=sqrt(1/3^2-1/6^2); hold on, subkoch(xl,xr,yl,yr,level,r); hold off; axis equal; axis tight; axis off; %------------------------------------------------------------ function subkoch(xl,xr,yl,yr,level,r) % draws a line at the lowest iteration level if level==1 plot([xl(1) xr(1)],[-yl(1) -yr(1)],b-) ,return end % branch in the next level level = level-1; % left subbranch xl(level) = xl(level+1); yl(level) = yl(level+1); xr(level) = 1/3*xr(level+1)+2/3*xl(level+1); yr(level) = 1/3*yr(level+1)+2/3*yl(level+1); subkoch(xl,xr,yl,yr,level,r); % middle left subbranch xl(level) = xr(level); yl(level) = yr(level); xr(level) = .5*xr(level+1)+.5*xl(level+1)-r*(yl(level+1)-yr(level+1)); yr(level) = .5*yr(level+1)+.5*yl(level+1)+r*(xl(level+1)-xr(level+1)); subkoch(xl,xr,yl,yr,level,r); % middle right subbranch xl(level) = xr(level); yl(level) = yr(level); xr(level) = 2/3*xr(level+1)+1/3*xl(level+1); yr(level) = 2/3*yr(level+1)+1/3*yl(level+1); subkoch(xl,xr,yl,yr,level,r); % right subbranch xl(level) = xr(level); yl(level) = yr(level); xr(level) = xr(level+1); yr(level) = yr(level+1); subkoch(xl,xr,yl,yr,level,r); koch(4); 2. Sierpinski三角形 function mysierpinski(x,y,L,n) %x,y为三角形中心点坐标，L为三角形边长，n为递归深度 hold on if n==1 x1=x-L/2; y1=y-L*tan(pi/6)/2; %计算三角形顶点的坐标 x2=x+L/2; y2=y-L*tan(pi/6)/2; x3=x; y3=y+L*tan(pi/6)/2; plot([x1;x2],[y1;y2]); plot([x2;x3],[y2;y3]); plot([x3;x1],[y3;y1]) %画三角形的边 else x01=x-L/4; y01=y-L*tan(pi/6)/4; %计算小三角形中心的坐标 x02=x+L/4; y02=y-L*tan(pi/6)/4; x03=x; y03=y+L*tan(pi/6)/4; mysierpinski(x01,y01,L/2,n-1) %递归调用 mysierpinski(x02,y02,L/2,n-1) mysierpinski(x03,y03,L/2,n-1) end mysierpinski(10,10,3,4); 3.Cantor集 function cantor3(n) if nargin==0; n=6; end hold on; [x,L]=ccp(0,0.5,n,0); axis off; function [x,L]=ccp(x,L,n,y); plot([x-L,x+L],[y,y],r,linewidth,2); if n0; [x,L]=ccp(x-L/1.5,L/3,n-1,y-1); L=L*3;x=x+L/1.5; [x,L]=ccp(x+L/1.5,L