迭代与分形.docVIP

《数学实验》报告 学 院： 电子与信息学院 专业班级： 信息工程2班 学 号： 201130301153 姓 名： 黄 锌 实验名称： 迭代与分形 实验日期： 实验三 迭代与分形 实验目的与要求 了解分形几何的基本情况； 了解通过迭代方式产生分形图的方法； 了解matlab软件中简单的程序结构； 掌握matlab软件中plot, fill等函数的基本用法； 实验任务 对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。 自己构造生成元（要有创意），按照图形迭代的方式产生分形图，用计算机编制程序绘制出它的图形，并计算其分形维数。 设计过程 任务1：对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。 分析：要想实现该效果，必定要用到迭代。首先可以知道等边三角形本质上就是连在一起的三条等长线段，将对一条线段的动作循环执行3次即可。而对于每一条线段这就是一条Koch. (1)代码程序如下： function xuehua(k) for j=0:2 if j==0; p=[0,0;10,0]; elseif j==1; p=[5,-5*sqrt(3);0,0]; else j==2; p=[10,0;5,-5*sqrt(3)]; end n=1; A=[cos(pi/3),-sin(pi/3);sin(pi/3),cos(pi/3)]; for s=1:k j=0; for i=1:n q1=p(i,:); q2=p(i+1,:); d=(q2-q1)/3; j=j+1;r(j,:)=q1; j=j+1;r(j,:)=q1+d; j=j+1;r(j,:)=q1+d+d*A; j=j+1;r(j,:)=q1+2*d; end n=4*n; clear p p=[r;q2]; clear r end plot(p(:,1),p(:,2)) hold on; axis equal end (2)运行结果如下： 当k=0时， 当k=1时， 当k=2时， 当k=3时， 当k=4时， 当k=5时， (3) Koch雪花的面积 设三角形的边长为a 当k=0时， S= 当k=1时， S=+3 当k=2时， S=+3+12 当k=3时， S=+3+12+48 当k=4时， S=+3+12+48+192 ………. 当k=m时， S=+3+12+48+192+….+3 = = ） = 根据上述推导可以确定： (4) Koch雪花分形维数。 相似形个数：m=12， 边长放大倍数：c=3， 任务2：对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。 代码程序如下： function xuehua1(k) for j=0:15 if j==0; p=[0,0;10,0]; elseif j==1; p=[0,-10;0,0]; elseif j==2; p=[10,0;10,-10]; elseif j==3; p=[10,-10;0,-10]; elseif j==4; p=[0,-5;5,0]; elseif j==5; p=[5,0;10,-5] elseif j==6; p=[10,-5;5,-10]; elseif j==7; p=[5,-10;0,-5] elseif j==8; p=