《matlab生成分形树》由会员分享，可在线阅读，更多相关《matlab生成分形树（17页珍藏版）》请在金锄头文库上搜索。
1、分形图形的生成元生成,由生成元产生的分形是一种规则分形，是数学家按照一定规则构造出来的，相当于物理学的模型。构造特点：最终图形是按照一定规则R对初始F0不断修改得到的。,分形图形的生成元生成,分形树,生成方法：,选取一条直线段F0，将该线段三等分，在等分点上各画一条长度为原线段长度三分之一的线段，并与原线段成固定夹角，得到图形F1。,分形树,生成方法：,再将F1的每一段都按上述方法修改，直至无穷，则最后得到的极限图象，称为分形树。,分形树,MATLAB程序：,p=0 0;0 10; n=2; %初始节点为2 plot(p(:,1),p(:,2), k); Hold on；%图形保持 A=cos
2、(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,:);,分形树,MATLAB程序：,p1=p(j,:);p2=p(j+1,:);,p1,p2,分形树,MATLAB程序：,d=(p2-p1)/3; r(i,:)=p1;i=i+1; %矩阵r的第一行存储点p1 q1=p1+d;,q1,分形树,MATLAB程序：,r(i,:)=q1;i=i+1; %矩阵r的第二行存储点q1 r(i,:)=q1;i=i+1; %矩阵r的第三行存储点q1 q2=q1+d*A;,q2,分形树,MATL
3、AB程序：,r(i,:)=q2;i=i+1; %矩阵r的第四行存储点q2 r(i,:)=q1;i=i+1; %矩阵r的第五行存储点q1 q3=p1+2*d;,q3,分形树,MATLAB程序：,r(i,:)=q3;i=i+1; %矩阵r的第六行存储点q3 r(i,:)=q3;i=i+1; %矩阵r的第七行存储点q3 q4=q3+d*A;,q4,分形树,MATLAB程序：,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
4、),xy(:,2),k); %联结q1 q2 Hold on；%图形保持 xy=q3;q4; plot(xy(:,1),xy(:,2),k); %联结q3 q4 Hold on；%图形保持 end,分形树,MATLAB程序：,p=r; %更新节点矩阵 n=5*n; %更新节点数 end axis equal axis off,分形树,levy曲线,主型：,如图所示，levy曲线的主型是以单位长线段为斜边的等腰直角三角形去掉斜边后留下的两条直角边且向上凸起。,levy曲线,MATLAB程序：,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;,levy曲线,MATLAB程序：,r(i,:)=q1;i=i+1; r(i,:)=q1;i=i+1; r(i,:)=p2;i=i+1; end n=n*2;p=r; end for m=1:2:n plot(p(m,1) p(m+1,1),p(m,2) p(m+1,2),k); hold on; end axis equal,levy曲线,
