几个分形的matlab实现.pdfVIP

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