分形樹__Matlab

1、%這是一個(gè)生成樹的主函數(shù)，它的輸入分別為每叉樹枝的縮短比、樹枝的偏角、生長次數(shù)%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%小提示：若用做函數(shù)，請(qǐng)將虛線框內(nèi)語句刪去。functionf=tree(w,dtheata,NN)%虛線框%cle

2、ar;clc;clf;w=0.8;dtheata二pi/6;NN=8;%建議生長次數(shù)NN不要超過10%虛線框%n=2NN;%從主枝算起，共需生成2NN個(gè)樹枝forNNK=1:nx1=0;y1=0;r1=1;theata1=pi/2;dataway=ten2twoN(NNK,NN);%把每一個(gè)樹枝的編號(hào)轉(zhuǎn)化為一個(gè)NN位的二進(jìn)制數(shù)forNNL=1:NNifdataway(NNL)=0 x2,y2,r2,theata2=antmoveleft(xl,yl,rl,theatal,w,dtheata);%若路徑數(shù)組上對(duì)應(yīng)的數(shù)字為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);%否則，數(shù)字為1，向右生長x1=x2;y1=y2;r1=r2;theata1=theata2;holdon%pause(eps)endendendholdoff%這是一個(gè)十進(jìn)制轉(zhuǎn)換為二進(jìn)制的函數(shù)，適用于正整數(shù)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%該函數(shù)的功能是將十進(jìn)制轉(zhuǎn)化為指定位皿NN)的二進(jìn)制數(shù)functiondataway=ten2twoN(NNK,NN);dataway=ten2two(NNK);a,b=size(dataway);Ncc=NN-b;%需要補(bǔ)充的0的個(gè)數(shù)forcc=1:Nccdataway=0,dataway;end%這是一個(gè)顯示平面中點(diǎn)的運(yùn)動(dòng)的函數(shù)，它的輸入為起始點(diǎn)的坐標(biāo)、%將要運(yùn)動(dòng)的方位角、將要運(yùn)動(dòng)的長度functionx2,y2=antmove(x1,y1,r1,theata)x2=x1+r1*cos(theata);y2=y1+r1*sin(theata);plot(x

5、1,x2,y1,y2)%這是一個(gè)顯示平面中點(diǎn)的運(yùn)動(dòng)的函數(shù)，它的輸入為起始點(diǎn)的坐標(biāo)、%將要運(yùn)動(dòng)的方位角、將要運(yùn)動(dòng)的長度%參見函數(shù)antmove%與antmove不同的是，該函數(shù)的返回值中多了兩個(gè)變量%這兩個(gè)變量r2,theata2為下一步點(diǎn)的移動(dòng)準(zhǔn)備了數(shù)據(jù)%theata角增大，表示點(diǎn)向左移動(dòng)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;%這是一個(gè)顯示平面中點(diǎn)的運(yùn)動(dòng)的函數(shù)，它的輸入為起始點(diǎn)的坐標(biāo)、%將要運(yùn)動(dòng)的方位角、將要運(yùn)動(dòng)的長度%參見函數(shù)antmove%與antmove不同的是，該函數(shù)的返回值中多了兩個(gè)變量%這兩個(gè)變量r2,theata2為下一步點(diǎn)的移動(dòng)準(zhǔn)備了數(shù)據(jù)%theata角減小，表示點(diǎn)向右移動(dòng)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文件中再運(yùn)行！%注意：把這些程序全部保存在名為tree的M文件中再運(yùn)行！%注意：把這些程序全部保存在名為tr