资源描述：
《几个分形的matlab实现》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、几个分形的matlab实现摘要：给出几个分形的实例，并用matlab编程实现方便更好的理解分形，欣赏其带来的数学美感关键字：Koch曲线实验图像一、问题描述：从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成山丘形图形如下图1在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边代替，再次形成新的图形如此迭代，形成Koch分形曲线。二、算法分析：考虑由直线段（2个点）产生第一个图形（5个点）的过程。图1中，设和分别为原始直线段的两个端点，现需要在直线段的中间依次插入三个点，，
2、。显然位于线段三分之一处，位于线段三分之二处，点的位置可看成是由点以点为轴心，逆时针旋转600而得。旋转由正交矩阵实现。算法根据初始数据（和点的坐标），产生图1中5个结点的坐标。结点的坐标数组形成一个矩阵，矩阵的第一行为的坐标，第二行为的坐标……，第五行为的坐标。矩阵的第一列元素分别为5个结点的坐标，第二列元素分别为5个结点的坐标。进一步考虑Koch曲线形成过程中结点数目的变化规律。设第次迭代产生的结点数为，第次迭代产生的结点数为，则和中间的递推关系为。三、实验程序及注释：p=[00;100];%P为初始两个点的坐标,第
3、一列为x坐标,第二列为y坐标n=2;%n为结点数A=[cos(pi/3)-sin(pi/3);sin(pi/3)cos(pi/3)];%旋转矩阵fork=1:4d=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;
4、%用向量方法计算迭代后处于4k+2位置上的点的坐标p(3:4:m,:)=q+d+d*A';%用向量方法计算迭代后处于4k+3位置上的点的坐标p(4:4:m,:)=q+2*d;%用向量方法计算迭代后处于4k位置上的点的坐标n=m;%迭代后新的结点数目endplot(p(:,1),p(:,2))%绘出每相邻两个点的连线axis([010010])四、实验数据记录：由第三部分的程序，可得到如下的Koch分形曲线：图2五、注记：１．参照实验方法，可绘制如下生成元的Koch分形曲线：图3此时，旋转矩阵为：程序和曲线如下：p=[00
5、;100];%P为初始两个点的坐标,第一列为x坐标,第二列为y坐标n=2;%n为结点数A=[0-1;10];%旋转矩阵fork=1:4d=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;%用向量方法计算迭代后处于
6、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;%迭代后新的结点数目endplot(p(:,1),p(:,2))%绘出每相邻两个点的连线axis([010010])图4由于中间三分之一部分是一个正方形时，有很多连接的部分。所以我们将高度压缩到原来的0.7倍，即中间部分为一个长与
7、宽之比为1:0.7的矩形时，得到程序和曲线如下：p=[00;100];%P为初始两个点的坐标,第一列为x坐标,第二列为y坐标n=2;%n为结点数A=[0-1;10];%旋转矩阵fork=1:4d=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位置上的点的坐标为迭代前的相应坐
8、标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;%
