科赫雪花曲線的MATLAB編程實(shí)現(xiàn)

1、 科赫雪花曲線的MATLAB編程實(shí)現(xiàn)2.1 經(jīng)一次迭代的科赫曲線MATLAB實(shí)現(xiàn)程序如下：x1=1 2 2.5 3 4; y1=0 0 0 0 0; h1=plot(x1,y1,'linewidth',2,'erasemode','xor'); axis equal axis off for g=linspace(0,1,40)*sin(pi/3); y1(3)=g; set(h1,'ydata',y1); drawnow; end 一次迭代所得科赫曲線如圖一： 圖一：2.2 經(jīng)二次迭代的科赫曲線MATLAB 實(shí)現(xiàn)程序如下：x2=

2、x1(1); y2=y1(1); for k=2:length(x1); t=linspace(x1(k-1),x1(k),4) ; tt=t(2),mean(t),t(3:4); x2=x2,tt; t=linspace(y1(k-1),y1(k),4); tt=t(2),mean(t),t(3:4); y2=y2,tt; end A=angle(y2(4:4:end)-y2(2:4:end)*i+(x2(4:4:end)-x2(2:4:end); for g=linspace(0,1,40)*sin(pi/3)/3; y2(3:4:end)=(y2(4:4:end)+y2(2:4:end)

3、/2+imag(g*exp(i*(A+pi/2); x2(3:4:end)=(x2(4:4:end)+x2(2:4:end)/2+real(g*exp(i*(A+pi/2) ; set(h1,'ydata',y2,'xdata',x2); drawnow; end 二次迭代后所得科赫曲線如圖二： 圖二2.3 經(jīng)三次迭代的科赫曲線MATLAB 實(shí)現(xiàn)程序如下x3=x2(1); y3=y2(1); for k=2:length(x2); t=linspace(x2(k-1),x2(k),4); tt=t(2),mean(t),t(3:4); x3=x3,tt; t=l

4、inspace(y2(k-1),y2(k),4); tt=t(2),mean(t),t(3:4); y3=y3,tt; end A=angle(y3(4:4:end)-y3(2:4:end)*i+(x3(4:4:end)-x3(2:4:end); for g=linspace(0,1,40)*sin(pi/3)/9; y3(3:4:end)=(y3(4:4:end)+y3(2:4:end)/2+imag(g*exp(i*(A+pi/2); x3(3:4:end)=(x3(4:4:end)+x3(2:4:end)/2+real(g*exp(i*(A+pi/2); set(h1,'ydat

5、a',y3,'xdata',x3); drawnow; end 三次迭代后所得科赫曲線如圖三： 圖三2.4經(jīng)四次迭代的科赫曲線MATLAB 實(shí)現(xiàn)程序如下x4=x3(1); y4=y3(1); for k=2:length(x3); t=linspace(x3(k-1),x3(k),4); tt=t(2),mean(t),t(3:4); x4=x4,tt; t=linspace(y3(k-1),y3(k),4); tt=t(2),mean(t),t(3:4); y4=y4,tt; end A=angle(y4(4:4:end)-y4(2:4:end)*i+(x4(4:4:

6、end)-x4(2:4:end); for g=linspace(0,1,40)*sin(pi/3)/27; y4(3:4:end)=(y4(4:4:end)+y4(2:4:end)/2+imag(g*exp(i*(A+pi/2); x4(3:4:end)=(x4(4:4:end)+x4(2:4:end)/2+real(g*exp(i*(A+pi/2); set(h1,'ydata',y4,'xdata',x4); drawnow; end四次迭代后所得科赫曲線如圖四： 圖四2.5 經(jīng)五次迭代的科赫曲線MATLAB 實(shí)現(xiàn)程序如下x5=x4(1);y5=y4(1)

7、;for k=2:length(x4); t=linspace(x4(k-1),x4(k),4); tt=t(2),mean(t),t(3:4); x5=x5,tt; t=linspace(y4(k-1),y4(k),4); tt=t(2),mean(t),t(3:4); y5=y5,tt; end A=angle(y5(4:4:end)-y5(2:4:end)*i+(x5(4:4:end)-x5(2:4:end); for g=linspace(0,1,40)*sin(pi/3)/81; y5(3:4:end)=(y5(4:4:end)+y5(2:4:end)/2+imag(g*exp(i*(A+pi/2); x5(3:4:end)=(x5(4:4:end)+x5