几个分形的matlab实现

精品文档 1欢迎下载 几个分形的 matlab 实现 摘要 给出几个分形的实例 并用 matlab 编程实现方便更好的理解分形 欣赏其带来 的数学美感 关键字 Koch 曲线 实验 图像 一 问题描述 从一条直线段开始 将线段中间的三分之一部分用一个等边三角形的两边代替 形成 山丘形图形如下 图1 在新的图形中 又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条 边代替 再次形成新的图形如此迭代 形成Koch分形曲线 二 算法分析 考虑由直线段 2个点 产生第一个图形 5个点 的过程 图1中 设和分别为 1 P 5 P 原始直线段的两个端点 现需要在直线段的中间依次插入三个点 显然位 2 P 3 P 4 P 2 P 于线段三分之一处 位于线段三分之二处 点的位置可看成是由点以点为轴心 4 P 3 P 4 P 2 P 逆时针旋转60 0而得 旋转由正交矩阵 3 cos 3 sin 3 sin 3 cos A 实现 算法根据初始数据 和点的坐标 产生图1中5个结点的坐标 结点的坐标数组 1 P 5 P 形成一个矩阵 矩阵的第一行为的坐标 第二行为的坐标 第五行为的25 1 P 2 P 5 P 坐标 矩阵的第一列元素分别为5个结点的坐标 第二列元素分别为5个结点的坐标 xy 进一步考虑Koch曲线形成过程中结点数目的变化规律 设第次迭代产生的结点数为k 第次迭代产生的结点数为 则和中间的递推关系为 k n1 k 1 k n k n 1 k n34 1 kk nn 精品文档 2欢迎下载 三 实验程序及注释 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就计算出每个向量长度的三分之一 与题中将线段三等分对应 m 4 n 3 迭代公式 q p 1 n 1 以原点为起点 前n 1个点的坐标为终点形成向量 p 5 4 m p 2 n 迭代后处于4k 1位置上的点的坐标为迭代前的相应坐标 p 2 4 m q d 用向量方法计算迭代后处于4k 2位置上的点的坐标 p 3 4 m q d d A 用向量方法计算迭代后处于4k 3位置上的点的坐标 p 4 4 m q 2 d 用向量方法计算迭代后处于4k位置上的点的坐标 n m 迭代后新的结点数目 end plot p 1 p 2 绘出每相邻两个点的连线 axis 0 10 0 10 四 实验数据记录 由第三部分的程序 可得到如下的Koch分形曲线 图2 五 注记 参照实验方法 可绘制如下生成元的Koch 分形曲线 精品文档 3欢迎下载 图3 此时 旋转矩阵为 01 10 2 cos 2 sin 2 sin 2 cos A 程序和曲线如下 p 0 0 10 0 P为初始两个点的坐标 第一列为x坐标 第二列为y坐标 n 2 n为结点数 A 0 1 1 0 旋转矩阵 for k 1 4 d diff p 3 diff计算相邻两个点的坐标之差 得到相邻两点确定的向量 则d就计算出每个向量长度的三分之一 与题中将线段三等分对应 m 5 n 4 迭代公式 q p 1 n 1 以原点为起点 前n 1个点的坐标为终点形成向量 p 6 5 m p 2 n 迭代后处于5k 1位置上的点的坐标为迭代前的相应坐标 p 2 5 m q d 用向量方法计算迭代后处于5k 2位置上的点的坐标 p 3 5 m q d d A 用向量方法计算迭代后处于5k 3位置上的点的坐标 p 4 5 m q 2 d d A 用向量方法计算迭代后处于5k 4位置上的点的坐标 p 5 5 m q 2 d 用向量方法计算迭代后处于5k位置上的点的坐标 n m 迭代后新的结点数目 end plot p 1 p 2 绘出每相邻两个点的连线 axis 0 10 0 10 图4 由于中间三分之一部分是一个正方形时 有很多连接的部分 所以我们将高度压缩到原来 的0 7倍 即中间部分为一个长与宽之比为1 0 7的矩形时 得到程序和曲线如下 p 0 0 10 0 P为初始两个点的坐标 第一列为x坐标 第二列为y坐标 精品文档 4欢迎下载 n 2 n为结点数 A 0 1 1 0 旋转矩阵 for k 1 4 d diff p 3 diff计算相邻两个点的坐标之差 得到相邻两点确定的向量 则d就计算出每个向量长度的三分之一 与题中将线段三等分对应 m 5 n 4 迭代公式 q p 1 n 1 以原点为起点 前n 1个点的坐标为终点形成向量 p 6 5 m p 2 n 迭代后处于5k 1位置上的点的坐标为迭代前的相应坐标 p 2 5 m q d 用向量方法计算迭代后处于5k 2位置上的点的坐标 p 3 5 m q d 0 7 d A 用向量方法计算迭代后处于5k 3位置上的点的坐标 p 4 5 m q 2 d 0 7 d A 用向量方法计算迭代后处于5k 4位置上的点的坐标 p 5 5 m q 2 d 用向量方法计算迭代后处于5k位置上的点的坐标 n m 迭代后新的结点数目 end plot p 1 p 2 绘出每相邻两个点的连线 axis 0 10 0 10 图 参照实验方法 我们由四边形的四个初始点出发 对于四边形的每条边 生成元如 下 精品文档 5欢迎下载 图 6 可得到火焰般的图形 程序和曲线如下 p 0 10 10 0 0 10 10 0 0 10 P为四边形四个顶点的坐标 其中第五个点与第一个点重合 以便于绘图 第一列为x坐标 第二列为y坐标 n 5 n为结点数 A cos pi 3 sin pi 3 sin pi 3 cos pi 3 旋转矩阵 顺时针旋转60度 for k 1 5 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 p 3 4 m q 2 d d A p 4 4 m q 2 d n m end plot p 1 p 2 axis 10 10 10 10 精品文档 6欢迎下载 图 7 参照实验方法 由下列的生成元 绘制 Koch 分形曲线 图 8 分析 为了绘图方便 我们将结点数处理一下 把第一次迭代产生的六个点看成十个 点 即图中有五条线段 我们将每条线段 的每个端点看成新的两个结点 这样我们就可以很方便地用 plot 绘图了 程序和曲线如下 p 0 0 10 10 P为初始两个点的坐标 第一列为x坐标 第二列为y坐标 n 2 n为结点数 A cos pi 3 sin pi 3 sin pi 3 cos pi 3 B cos pi 3 sin pi 3 sin pi 3 cos pi 3 旋转矩阵A对应于第一次逆时针旋转60度 旋转矩阵B对应于第二次顺时针旋转60度 for k 1 4 d diff p 3 d1 d 1 2 n 取每条线段对应的向量 m 5 n 迭代公式 q1 p 1 2 n 1 p 10 10 m p 2 2 n 精品文档 7欢迎下载 p 1 10 m p 1 2 n 迭代后处于10k与10k 1位置上的点的坐标为迭代前的相应坐标 p 2 10 m q1 d1 用向量方法计算迭代后处于10k 2 10k 3 10k 5位置上的点的坐标 都相同 p 3 10 m p 2 10 m p 4 10 m q1 d1 d1 A 用向量方法计算迭代后处于10k 4位置上的点的坐标 p 5 10 m p 2 10 m p 6 10 m q1 2 d1 用向量方法计算迭代后处于10k 6 10k 7 10k 9位置上的点的坐标 都相同 p 7 10 m p 6 10 m p 8 10 m q1 2 d1 d1 B p 9 10 m p 6 10 m n m 迭代后