分形树__Matlab

1、%这是一个生成树的主函数，它的输入分别为每叉树枝的缩短比、树枝的偏角、生长次数%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%小提示：若用做函数，请将虚线框内语句删去。functionf=tree(w,dtheata,NN)%虚线框%cle

2、ar;clc;clf;w=0.8;dtheata二pi/6;NN=8;%建议生长次数NN不要超过10%虚线框%n=2NN;%从主枝算起，共需生成2NN个树枝forNNK=1:nx1=0;y1=0;r1=1;theata1=pi/2;dataway=ten2twoN(NNK,NN);%把每一个树枝的编号转化为一个NN位的二进制数forNNL=1:NNifdataway(NNL)=0 x2,y2,r2,theata2=antmoveleft(xl,yl,rl,theatal,w,dtheata);%若路径数组上对应的数字为0，则向左生长x1=x2;y1=y2;r1=r2;theata1=theat

3、a2;holdon%pause(eps)elsex2,y2,r2,theata2=antmoveright(xl,yl,rl,theatal,w,dtheata);%否则，数字为1，向右生长x1=x2;y1=y2;r1=r2;theata1=theata2;holdon%pause(eps)endendendholdoff%这是一个十进制转换为二进制的函数，适用于正整数functionyushu=ten2two(x)yushu=;shang=2;ifx=1yushu=1;elsewhileshang1shang=fix(x/2);yushu=mod(x,2),yushu;x=shang;end

4、yushu=1,yushu;end%该函数的功能是将十进制转化为指定位皿NN)的二进制数functiondataway=ten2twoN(NNK,NN);dataway=ten2two(NNK);a,b=size(dataway);Ncc=NN-b;%需要补充的0的个数forcc=1:Nccdataway=0,dataway;end%这是一个显示平面中点的运动的函数，它的输入为起始点的坐标、%将要运动的方位角、将要运动的长度functionx2,y2=antmove(x1,y1,r1,theata)x2=x1+r1*cos(theata);y2=y1+r1*sin(theata);plot(x

5、1,x2,y1,y2)%这是一个显示平面中点的运动的函数，它的输入为起始点的坐标、%将要运动的方位角、将要运动的长度%参见函数antmove%与antmove不同的是，该函数的返回值中多了两个变量%这两个变量r2,theata2为下一步点的移动准备了数据%theata角增大，表示点向左移动functionx2,y2,r2,theata2=antmoveleft(x1,y1,r1,theata1,w,dtheata)x2=x1+r1*cos(theata1);y2=y1+r1*sin(theata1);plot(x1,x2,y1,y2);r2=r1*w;theata2=theata1+dthea

6、ta;%这是一个显示平面中点的运动的函数，它的输入为起始点的坐标、%将要运动的方位角、将要运动的长度%参见函数antmove%与antmove不同的是，该函数的返回值中多了两个变量%这两个变量r2,theata2为下一步点的移动准备了数据%theata角减小，表示点向右移动functionx2,y2,r2,theata2=antmoveright(x1,y1,r1,theata1,w,dtheata)x2=x1+r1*cos(theata1);y2=y1+r1*sin(theata1);plot(x1,x2,y1,y2);r2=r1*w;theata2=theata1-dtheata;%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tree的M文件中再运行！%注意：把这些程序全部保存在名为tr