MATLAB生成分形树.ppt

分形图形的生成元生成 分形图形的生成元生成 分形树 生成方法： 分形树 生成方法： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 MATLAB程序： 分形树 levy曲线 主型： levy曲线 MATLAB程序： levy曲线 MATLAB程序： levy曲线 * 由生成元产生的分形是一种规则分形，是数学家按照一定规则构造出来的，相当于物理学的模型。构造特点：最终图形是按照一定规则R对初始F0不断修改得到的。 选取一条直线段F0，将该线段三等分，在等分点上各画一条长度为原线段长度三分之一的线段，并与原线段成固定夹角，得到图形F1。 再将F1的每一段都按上述方法修改，直至无穷，则最后得到的极限图象，称为分形树。 p=[0 0;0 10]; n=2; %初始节点为2 plot(p(:,1),p(:,2), k); Hold on；%图形保持 A=[cos(pi/6) -sin(pi/6);sin(pi/6) cos(pi/6)]; for k=1:4 i=1; for j=1:2:n p1=p(j,:);p2=p(j+1,:); p1=p(j,:);p2=p(j+1,:); p1 p2 d=(p2-p1)/3; r(i,:)=p1;i=i+1; %矩阵r的第一行存储点p1 q1=p1+d; p1 p2 q1 r(i,:)=q1;i=i+1; %矩阵r的第二行存储点q1 r(i,:)=q1;i=i+1; %矩阵r的第三行存储点q1 q2=q1+d*A; p1 p2 q1 q2 r(i,:)=q2;i=i+1; %矩阵r的第四行存储点q2 r(i,:)=q1;i=i+1; %矩阵r的第五行存储点q1 q3=p1+2*d; p1 p2 q1 q2 q3 r(i,:)=q3;i=i+1; %矩阵r的第六行存储点q3 r(i,:)=q3;i=i+1; %矩阵r的第七行存储点q3 q4=q3+d*A; p1 p2 q1 q2 q3 q4 r(i,:)=q4;i=i+1; %矩阵r的第八行存储点q4 r(i,:)=q3;i=i+1; %矩阵r的第九行存储点q3 r(i,:)=p2;i=i+1; %矩阵r的第十行存储点p2 xy=[q1;q2]; plot(xy(:,1),xy(:,2),‘k’); %联结q1 q2 Hold on；%图形保持 xy=[q3;q4]; plot(xy(:,1),xy(:,2),k); %联结q3 q4 Hold on；%图形保持 end p=r; %更新节点矩阵 n=5*n; %更新节点数 end axis equal axis off 如图所示，levy曲线的主型是以单位长线段为斜边的等腰直角三角形去掉斜边后留下的两条直角边且向上凸起。 p=[0 0;10 0];n=2; A=[cos(pi/4) -sin(pi/4);sin(pi/4) cos(pi/4)]; for k=1:10 i=1; for j=1:2:n p1=p(j,:);p2=p(j+1,:); d=(p2-p1)*sqrt(2)/2; r(i,:)=p1;i=i+1; q1=p1+d*A; r(i,:)=q1;i=i+1; r(i,:)=q1;i=i+1; r(i,:)=p2;i=i+1; end n=n*2;p