分形图形迭代生成与性质分析.pdfVIP

实验四迭代问题二——利用迭代法制作

分形图形

班级：统计0502学号：1303050206：欧钊锋

实验目的：

⚫以迭代的方法介绍分形的基本特性以及用函数作出分形图形的方法，使我们对其有一定

的了解

实验理论：

先利用函数形的制作分形图生成元，然后运用迭代法和line函数生成分形图形

实验内容：

1.用计算机绘出koch曲线，sierpinski三角形及花草树木的图形，你能否自己构造一些生成

元，并由此绘出相应的分形图行？任取分形图形的一个局部并将它放大，它同原来的分形图

形有什么关系？

2.从一个正三角形出发，用koch曲线的生成元做迭代得到的极限图形称为koch雪花曲线

（1）试计算雪花曲线的边长及面积，它们是否有限？你如何解释所得出的结论？

（2）雪花曲线是否光滑（即每一点是否有曲线存在？）

（3）其它的一些分形是否具有类似的性质？

3假设生成元由M条线段构成，每段长度为原线段长度的1/C，定义该分形的维数为

d=log(n)/log(C)。想，这样定义分形的维数有什么道理？

程序设计：

Koch雪花图形

functionouzhaoxuehua(x1,y1,x2,y2,n)

ifn1

ouzhaoxuehua(x1,y1,(2*x1+x2)/3,(2*y1+y2)/3,n-1)

ouzhaoxuehua((2*x1+x2)/3,(2*y1+y2)/3,(x1+x2)/2+(y1-y2)*sqrt(3)/6,(y1+y2)/2-(x1-x2)*sqrt(3)/

6,n-1)

ouzhaoxuehua((x1+x2)/2+(y1-y2)*sqrt(3)/6,(y1+y2)/2-(x1-x2)*sqrt(3)/6,(x1+2*x2)/3,(y1+2*y2)/

3,n-1)

ouzhaoxuehua((x1+2*x2)/3,(y1+2*y2)/3,x2,y2,n-1)

else

x=[x1,(2*x1+x2)/3,(x1+x2)/2+(y1-y2)*sqrt(3)/6,(x1+2*x2)/3,x2]

y=[y1,(2*y1+y2)/3,(y1+y2)/2-(x1-x2)*sqrt(3)/6,(y1+2*y2)/3,y2]

line(x,y)

end

axisequal

functionouzhaohuacaoshumu(x1,y1,x2,y2,n)

ifn1

ouzhaohuacaoshumu(x1,y1,(2*x1+x2)/3,(2*y1+y2)/3,n-1)

ouzhaohuacaoshumu((2*x1+x2)/3,(2*y1+y2)/3,(x1+x2)/2+(y1-y2)/(4*sin(pi/4)),(y1+y2)/2-(x1-x2

)/(4*sin(pi/4)),n-1)

ouzhaohuacaoshumu((2*x1+x2)/3,(2*y1+y2)/3,(x1+2*x2)/3,(y1+2*y2)/3,n-1)

ouzhaohuacaoshumu((x1+2*x2)/3,(y1+2*y2)/3,(x1+5*x2)/6-(y1-

y2)/(4*sin(pi/4)),(y1+5*y2)/6+(x1-x2)/(4*sin(pi/4)),n-1)

ouzhaohuacaoshumu((x1+2*x2)/3,(y1+2*y2)/3,x2,y2,n-1)

else

x=[x1,(2*x1+x2)/3,(x1+x2)/2+(y1-y2)/(4*sin(pi/4)),(2*x1+x2)/3,(x1+2*x2)/3,(x1+5*x2)/6-(y1-y

2)/(4*sin(pi/4)),(x1+2*x2)/3,x2]

y=[y1,(2*y1+y2)/3,(y1+y2)/2-(x1-x2)/(4*sin(pi/4)),