【Matlab代码】Koch曲线

步骤：
1，将基本的图形即线段用数组表示出来。由于是线段，我们只需要列出端点坐标即可；
2，将第一步得到的数组做如下处理：首先将线段缩小到原来的三分之一。然后将得到的线段复制四份：将第一段放在最左边，将第二段逆时针旋转60°，左端点与第一段的右端点相接；将第三段顺时针旋转60°，左端点与第二段的右端点相接；最后将第四段左端点与第三段右端点相接。这样我们得到了最初屏幕上的第五个图。

function Koch(k）
p=[0 10];
for m=1:k
q1=p/3;
q2=10/3+(q1.*exp(i*pi/3));
q3=(10/3+10/3*exp(li*pi/3))+(q1.*exp(-li*pi/3));
q4=20/3+p/3;
p=[q1 q2 q3 q4];
end
figure
plot(p)
axis equal

function Koch(k）
p=[0 10];
for m=1:k
q1=p/3;
q2=10/3+(q1.*exp(i*pi/3));
q3=(10/3+10/3*exp(li*pi/3))+(q1.*exp(-li*pi/3));
q4=20/3+p/3;
p=[q1 q2 q3 q4];
end
figure
a=real(p);
b=image(p);
fill(a,b,'o')
axis equal