Koch分形曲线分析和总结.docx

Koch分形曲线

分形原理

这是一类复杂的平面曲线，可用算法描述。从一条直线段开始，将线段中间三分之一部分用等边三角形的两条边代替，形成具有5个结点的图形(图1)；在新的图形中，又将图中每一直线段中间的三分之一部分都用一等边三角形的两条边代替，再次形成新的图形(图2)，这时，图形中共有17个结点。

图1 第一次迭代 图2 第二次迭代

这种迭代继续进行下去可以形成Koch分形曲线。在迭代过程中，图形中的点将越来越多，而曲线最终显示细节的多少将取决于迭代次数和显示系统的分辩率。

算法分析

算法分析：考虑由直线段(2个点)产生第一个图形(5个点)的过程。设P和P

1 5

分别为原始直线段的两个端点。现在需要在直线段的中间依次插入三个点

P,P,P 产生第一次迭代的图形(图1)。显然，P位于P点右端直线段的三分之

2 3 4 2 1

一处，P

位于P点右端直线段的三分之二处；而P

点的位置可以看成是由P点

4 1 3 4

绕P旋转60度(逆时针方向)而得到的，故可以处理为向量PP

经正交变换而得

2 24

到向量PP。算法如下：

23

P

2

?P?(P

1 5

?P)/3；

1

P

4

?P?2(P

1 5

?P)/3；

1

(3) P?P?(P?P)?AT；

3 2 4 2

在(3)中，A为正交矩阵：

? ? ??

?cos ? sin?

A?? 3 3?

?sin? cos??

?? 3 3??

算法根据初始数据(P和P点的坐标)，产生图1中5个结点的坐标。结点的

1 5

坐标数组形成一个 5×2 矩阵，矩阵的第一行为 P

1

的坐标，第二行为P的坐

2

标，……，第五行为P

5

的坐标。矩阵的第一列元素分别为5个结点的X坐标，

第二列元素分别为5个结点的Y坐标。

进一步考虑Koch曲线形成过程中结点数目的变化规律。设第k次迭代产生

结点数为n

k

，第k+1次迭代产生结点数为n

k?1

，则n

k

和n

k?1

之间的递推关系式为

n

k?1

?4n

k

3。

MATLAB实现

p=[0 0;10 0];n=2;

A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)];for k=1:5

d=diff(p)/3;m=4*n-3;

q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);p(2:4:m,:)=q+d;

p(3:4:m,:)=q+d+d*A;

p(4:4:m,:)=q+2*d;

n=m;

endplot(p(:,1),p(:,2),k)axisequal

axisoff

程序运行后，可得如下分形曲线图形：

图3第五次迭代图形

图1所示的简单图形被称为是koch曲线的生成元。下面给出其它几种生成元及对应曲线：

图4树形生成元及对应曲线

图5第二种树形生成元及对应曲线

矩形生成元

第一种生成方式

第二种生成方式

图6矩形生成元及两种对应曲线

附源程序

第一种树形源程序：

p=[00;1010];n=2;

line(p(:,1),p(:,2),Color,k);

A=[cos(pi/6)-sin(pi/6);sin(pi/6)cos(pi/6)];fork=1:4

i=1;

forj=1:2:n

p1=p(j,:);p2=p(j+1,:);

endp=r;

d=(p2-p1)/3;r(i,:)=p1;i=i+1;q1=p1+d;r(i,:)=q1;i=i+1;

r(i,:)=q1;i=i+1;q2=q1+d*A;r(i,:)=q2;i=i+1;

r(i,:)=q1;i=i+1;q3=p1+2*d;r(i,:)=q3;i=i+1;

r(i,:)=q3;i=i+1;q4=q3+d*A;r(i,:)=q4;i=i+1;

r(i,:)=q3;i=i+1;

r(i,:)=p2;i=i+1;xy=[q1;q2];

line(xy(:,1),xy(:,2),Color,k);xy=[q3;q4];line(xy(:,1),xy(:,2),Color,k);

n=5*n;

end

axisequalaxisoff

第二种树形源程序

clear

p=[00;015];

b=pi/8;e=1.5;

a=[cos(b)-sin(b);sin(b)cos(b)];

n=2;

fork=1:4

i=1;

forj=1:2:n

p1=p(j,:);

p2=p(j+1