matlab 实现分形图形-科赫曲线绘制

科赫曲线如下
https://upload.wikimedia.org/wikipedia/commons/f/fd/Von_Koch_curve.gif


matlab代码如下

x=[-1,0,1];

y=[0,sqrt(3),0];

c=complex(x,y);

c(4)=c(1);

plot(real(c), imag(c));

axis([-1.5 1.5 -1 2]);

axis square;

pause(1);

for i=1:8

    n=size(c);

    n=n(2);

    n=n-1;

    c1=complex(zeros(1,4*n),zeros(1,4*n));

    for j=1:n

        c1(4*j-3)=c(j);

        c1(4*j-2)=(2*c(j)+c(j+1))/3;

        c1(4*j)=(c(j)+2*c(j+1))/3;

        c1(4*j-1)=c1(4*j-2)+(c1(4*j)-c1(4*j-2))*complex(0.5,sqrt(3)/2);

    end

    c=c1;

    c(4*n+1)=c(1);

    plot(real(c), imag(c));

    axis([-1.5 1.5 -1 2]);

    axis square;

    pause(1);

end