几个分形的matlab实现资料

1、几个分形的matlab实现摘要：给出几个分形的实例，并用 matlab编程实现方便更好的理解分形，欣赏其带来的 数学美感关键字：Koch曲线 实验 图像一、问题描述：从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成山丘形图形如下在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边 代替，再次形成新的图形如此迭代，形成Koch分形曲线。二、算法分析：考虑由直线段（2个点）产生第一个图形（5个点）的过程。图1中，设R和P5分别为原始直线段的两个端点，现需要在直线段的中间依次插入三个点P2，P3，P4。显然F2位于线段三分之一处，P4位于线段三分之二

2、处，F3点的位置可看成是由 卩4点以卩2点为轴心，0逆时针旋转60而得。旋转由正交矩阵A-nCOS()3nsin()-si n(-)cosL)3实现。算法根据初始数据（P和P5点的坐标），产生图1中5个结点的坐标。结点的坐标数组形成一个5 2矩阵，矩阵的第一行为 R的坐标，第二行为 P2的坐标 ，第五行为P5的坐标。矩阵的第一列元素分别为5个结点的x坐标，第二列元素分别为 5个结点的y坐标。进一步考虑Koch曲线形成过程中结点数目的变化规律。设第k次迭代产生的结点数为m，第k1次迭代产生的结点数为nk1，则nk和nk1中间的递推关系为n1=4nk- 3。三、实验程序及注释：p=0 0;10 0

3、;%两初始两个点的坐标，第一列为x坐标，第二列为y坐标n=2;%曲结点数A=cos(pi/3) -si n( pi/3);s in (pi/3) cos(pi/3);%旋转矩阵for k=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;%用向量方法计算迭代后处于4k+

4、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(0 10 0 10)四、实验数据记录：由第三部分的程序，可得到如下的Koch分形曲线：图2五、注记：1参照实验方法，可绘制如下生成元的Koch分形曲线:图3此时，旋转矩阵为:%两初始两个点的坐标，第一列为x坐标，第二列为y坐标%为结点数%旋转矩阵程序和曲线如下:p=0 0;10 0;n=2;A=0 -1;1

5、0;for k=1:4d=diff(p)/3;%diff计算相邻两个点的坐标之差，得到相邻两点确定的向量%则d就计算岀每个向量长度的三分之一，与题中将线段三等分对应m=5* n-4;%迭代公式q=p(1:n-1,:);%以原点为起点，前n-1个点的坐标为终点形成向量p(2:5:m,:)=q+d; p(3:5:m,:)=q+d+d*A;p(4:5:m,:)=q+2*d+d*A; p(5:5:m,:)=q+2*d;n=m;end涮向量方法计算迭代后处于 涮向量方法计算迭代后处于 涮向量方法计算迭代后处于 涮向量方法计算迭代后处于 淞代后新的结点数目5k+2位置上的点的坐标5k+3位置上的点的坐标5

6、k+4位置上的点的坐标 5k位置上的点的坐标plot(p(:,1),p(:,2)淤岀每相邻两个点的连线p(6:5:m,:)=p(2:n,:);%迭代后处于5k+1位置上的点的坐标为迭代前的相应坐标axis(0 10 0 10)图4由于中间三分之一部分是一个正方形时，有很多连接的部分。所以我们将高度压缩到原来的0.7倍，即中间部分为一个长与宽之比为1:0.7的矩形时，得到程序和曲线如下：p=0 0;10 0;%两初始两个点的坐标，第一列为x坐标，第二列为y坐标n=2;%n结点数A=0 -1;1 0;%旋转矩阵for k=1:4d=diff(p)/3;%diff计算相邻两个点的坐标之差，得到相邻两

7、点确定的向量%则d就计算岀每个向量长度的三分之一，与题中将线段三等分对应m=5* n-4;%迭代公式q=p(1：n-1,:);%以原点为起点，前n-1个点的坐标为终点形成向量p(6:5:m,:)=p(2: n,:); p(2:5:m,:)=q+d; p(3:5:m,:)=q+d+0.7*d*A; p(4:5:m,:)=q+2*d+0.7*d*A;%迭代后处于5k+1位置上的点的坐标为迭代前的相应坐标 涮向量方法计算迭代后处于5k+2位置上的点的坐标涮向量方法计算迭代后处于5k+3位置上的点的坐标涮向量方法计算迭代后处于5k+4位置上的点的坐标涮向量方法计算迭代后处于淞代后新的结点数目5k位置上

8、的点的坐标p(5:5:m,:)=q+2*d;n=m;endplot(p(:,1),p(:,2)淤岀每相邻两个点的连线axis(0 10 0 10)图52 参照实验方法，我们由四边形的四个初始点出发，对于四边形的每条边，生成元如 下：可得到火焰般的图形。程序和曲线如下：p=0 10;10 0;0 -10;-10 0;0 10;%两四边形四个顶点的坐标，其中第五个点与第一个点重合，以便于绘图%第一列为x坐标，第二列为y坐标n=5; %曲结点数A=cos(-pi/3) -sin(-pi/3);sin(-pi/3) cos(-pi/3);%旋转矩阵，顺时针旋转 60 度for k=1:5d=diff(

9、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;endplot(p(:,1),p(:,2)axis(-10 10 -10 10)JU L|LIJ|i-1D-6-6.-4024 B 010图73 .参照实验方法，由下列的生成元，绘制Koch分形曲线:图8分析：为了绘图方便，我们将结点数处理一下，把第一次迭代产生的六个点看成十个点，即图中有五条线段(1 2,34,5 6,7 8,9 10)，我们将每条线段的每个端点看成新的两

10、个结点，这样我们就可以很方便地用plot绘图了。程序和曲线如下：p=0 0;10 10;%f为初始两个点的坐标，第一列为x坐标，第二列为y坐标n=2;%曲结点数A=cos(pi/3) -si n( pi/3);si n( 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:4d=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,:);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: