数学实验之数值仿真..docxVIP

《数学实验》报告 学 院：电子与信息学院 专业班级：信息工程2班 学 号： 201130301153姓 名：黄锌 实验名称：迭代与分形 实验日期：2013年4月16日星期二实验三迭代与分形实验目的与要求了解分形几何的基本情况；了解通过迭代方式产生分形图的方法；了解matlab软件中简单的程序结构；掌握matlab软件中plot, fill等函数的基本用法；实验任务对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。自己构造生成元（要有创意），按照图形迭代的方式产生分形图，用计算机编制程序绘制出它的图形，并计算其分形维数。设计过程任务1：对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。分析：要想实现该效果，必定要用到迭代。首先可以知道等边三角形本质上就是连在一起的三条等长线段，将对一条线段的动作循环执行3次即可。而对于每一条线段这就是一条Koch.(1)代码程序如下：functionxuehua(k)for j=0:2if 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)];endn=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; endn=4*n;clear pp=[r;q2];clear rendplot(p(:,1),p(:,2))hold on;axis equalend(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:15if 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=[2.5,-2.5;7.5,-2.5];elseif j==9; p=[7.5,-2.5;7.5,-7.5]elseif j==10; p=[7.5,-7.5;2.5,-7.5];elseif j==11; p=[2.5,-7.5;2.5,-2.5]elseif j==12; p=[5,-2.5;2.5,-5];elseif j==13; p=[2.5,-5;5,-7.5]elseif j==14; p=[5,-7.5;7.5,-5];else j==15; p=[7.5,-5;5,-2.5]endn=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; endn=4*n;clear pp=[r;q2];clear rendplot(p(:,1),p(:,2),r)hold on;axis equale