数学实验matlabKh雪花.docVIP

《数学实验》报告3分形实例 电气二班陆展辉201430222325 (51) 问题描述 对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形阁称为Koch 雪花。编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。 2、 对一条横向线段，先将其等分成4段，然后再将第二段向上移，将第三段向下 移，再将第四段的相邻端点连接起来，迭代一次后变成阁3-21.继续迭代得到的 分形图，称为Minkowski香肠。编制程序绘制出它的图形，并计算它的分形维数。 图3-21 M i nkows k i香肠一次迭代 问题分析与实验过程 实验过程.? 1.仿照Koch曲线代码对三角形的每条边进行Koch曲线化, 函数的输入参数有三角形的边长R和迭代次数k,输出Koch 雪花图形以及雪花所围面积 (1)代码如下： function xuehua(k) % k为迭代次数 forj=0:2 %依次对3条边进行Koch曲线运算 ifj==O; p=[0,0;10,0]; elseif j==l; p=[5,-5*sqrt(3);O,O]; else j==2; p=[10,0;5,-5*$qrt(3)]; end n=l； %存放线段的数量，初始值为1 A=[cos(pi/3),-sin(pi/3);sin(pi/3),cos(pi/3)]; % 用于计算新的结点 for s=l:k j二0; % j为行数 for i=l:n ql=p(i，:)； ％ S前线段的起点坐标 q2=p(i+l，:)； ％ A前线段的终点坐标 d=(q2-ql)/3; j=j+l;r(j,:)=ql; %原起点存入 r j=j+l;r(j，:)=ql+d; %新 1 点存入 r j=j+l;r(j,:)=ql+d+d*A,; %新 2 点存入 r j=j+l;r(j，:)=ql+2*d; %新 3 点存入 r end n=4*n; %全部线段迭代一次后，线段数量乘4 clear p %清空p，注意：最后一个终点q2不在r中 p=lr;q2J; 条边的全部结点 clear r end plot(p(:,l),p(:,2)) % 连接各个结点 hold on; axis equal end 不同k对应不同的图像如下: k=l 么2k=l 时 S= 么2 k=l 时 S= 4 + 12 k=2 时 S= 4 + 12 +27 (2) Koch雪花面积推导如下所示: k=0 时 S= 4 TOC \o 1-5 \h \z V3 , 73 , 73 2 4^3 , 厂* 厂- 广 厂- k=3 时 S=4 +12 + 27 + 243 nk=n n k=n 时 S= 4 + 12 + 每一次迭加，所产生的新三角形的边长变为上一次的3,数量为上一次的4倍. 73 9 73 ——厂細 —— S= 4 + 4 1、2 (3* (i) +……+3* 4-) * (士)2 ) _ _ n 1 V|r2 ￡[3*4(/-°*(-)2] =4 + 4 “=i 3 曲线总面积无穷大。 (3)综上所述可得Koch雪花的分形维数为: 根据迭代的规律得到： 相似形个数：m=6边长放大倍数：c=3, d = In m + In c. = In 6 + In 3 = i 631 2、绘制Minkowski香肠 (i)编辑实现题目迭代的函数 在Matlab中，编制一个函数来绘制Minkowski香肠的阁形。具体代码如下: function Minkowski (k) % 显示迭代k次后的Minkowski曲线图 p=[0, 0; 10, 0] ; %存放结点坐标，每行一个点，初始值为两结点的坐 n=l ；A=[0, 1;-1 n=l ； A=[0, 1;-1，0]; for s=l:k j=0; for i=l:n; %实现迭代过程，计算所有的结点的坐标 ql=P (i，：)； %目前线段的起点坐标 q2=p (i+1，：)； %目前线段的终点坐标 d=(q2-ql)/4； J-j+1 (j，：)=ql; % 原起点存入r j=j+l;r (j，：)=ql+d; % 新1 点存入r j=j+l;r (j, :)=ql+d+d*A; % 新2点存入r j=j+l；r (j, :)=ql+2*d+d*A; % 新3点存入r j=j+l;r (j, ：)=ql+2*d+d*A，； % 新4点存入r j=j+l；r (j, :)=ql+3*d+d*A’ ； % 新5点存入r j=j+l;r (j，：)=ql+3*d; % 新6点存入r end %原终点作为下条线段的起点，在迭代下条线段时存入 n=n*7； %全部线段迭代一次后，线段数量乘7 clear p %清空p，注意：最后一