几个分形的matlab实现.docxVIP

03 0 3 3 3 3 几个分形的 matlab 实现 摘要：给出几个分形的实例，并用 matlab 编程实现方便更好的理解分形，欣赏其带来的 数学美感 关键字： Koch 曲线 实验 图像 一、问题描述： 从一条直线段开始， 将线段中间的三分之一部分用一个等边三角形的两边代替， 形成山 丘形图形如下 图 1 在新的图形中， 又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边 代替，再次形成新的图形如此迭代，形成 Koch 分形曲线。 二、算法分析： 考虑由直线段（ 2个点）产生第一个图形（ 5个点）的过程。图 1中，设 P1 和 P5 分别为 原始直线段的两个端点， 现需要在直线段的中间依次插入三个点 P2， P3， P4。 显然 P2 位于 线段三分之一处， P4 位于线段三分之二处， 逆时针旋转 60 而得。旋转由正交矩阵 P3 点的位置可看成是由 P4 点以 P2 点为轴心， 实现。 算法根据初始数据（ P1和 cos( ) sin( ) A sin( ) cos( ) P5 点的坐标），产生图 1中5个结点的坐标。结点的坐标数组 形成一个 5 2 矩阵，矩阵的第一行为 P1 的坐标，第二行为 P2 的坐标……，第五行为 P5 的 坐标。矩阵的第一列元素分别为 5个结点的 x 坐标，第二列元素分别为 5个结点的 y 坐标。 进一步考虑 Koch 曲线形成过程中结点数目的变化规律。设第 k次迭代产生的结点数为 nk ，第 k 1次迭代产生的结点数为 nk 1 ，则 nk和 nk 1 中间的递推关系为 nk 1 4nk 3。 三、实验程序及注释： 1 / 7 p=[0 p=[0 0;10 0]; n=2; %P为初始两个点的坐标 %n为结点数 A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)]; for k=1:4 , 第一列为 x坐标 , 第二列为 y坐标 %旋转矩阵 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; %diff 计算相邻两个点的坐标之差 , 得到相邻两点确定的向量 %则d就计算出每个向量长度的三分之一 , 与题中将线段三等分对应 %迭代公式 %以原点为起点 , 前n-1 个点的坐标为终点形成向量 %迭代后处于 4k+1位置上的点的坐标为迭代前的相应坐标 %用向量方法计算迭代后处于 4k+2位置上的点的坐标 p(3:4:m,:)=q+d+d*A; p(4:4:m,:)=q+2*d; n=m; end %用向量方法计算迭代后处于 %用向量方法计算迭代后处于 %迭代后新的结点数目 4k+3位置上的点的坐标 4k位置上的点的坐标 plot(p(:,1),p(:,2)) %绘出每相邻两个点的连线 axis([0 10 0 10]) 四、实验数据记录： 由第三部分的程序，可得到如下的 Koch分形曲线： 图2 五、注记： １．参照实验方法，可绘制如下生成元的 Koch 分形曲线： 图3 2 / 7 此时，旋转矩阵为： 程序和曲线如下： p=[0 0;10 0]; n=2; A=[0 -1;1 0]; for k=1:4 d=diff(p)/3; m=5*n-4; q=p(1:n-1,:); p(6:5:m,:)=p(2:n,:); A 2 22 2cos( ) 2 2 2 2 sin( ) cos( ) 0 1 1 0 %P为初始两个点的坐标 , 第一列为 x坐标 , 第二列为 y坐标 %n为结点数 %旋转矩阵 %diff 计算相邻两个点的坐标之差 , 得到相邻两点确定的向量 %则d就计算出每个向量长度的三分之一 , 与题中将线段三等分对应 %迭代公式 %以原点为起