分形幾何中一些經(jīng)典圖形的matlab畫法-文庫(kù)吧
2025-04-01 23:40 本頁(yè)面
【正文】 w=levy2(X(1+3*(j1):3*j),Y(1+3*(j1):3*j))。 [XX(3*2*(j1)+1:3*2*(j1)+3),YY(3*2*(j1)+1:3*2*(j1)+3)]=levy1(w(1,1),w(1,2),w(1,3),w(1,4))。 [XX(3*2*(j1)+3+1:3*2*(j1)+3+3),YY(3*2*(j1)+3+1:3*2*(j1)+3+3)]=levy1(w(2,1),w(2,2),w(2,3),w(2,4))。 end X=XX。 Y=YY。endplot(X,Y)hold onaxis equal%由以(x1,y1),(x2,y2)為端點(diǎn)的線段生成新的中間點(diǎn)坐標(biāo)并把(x1,y1),(x2,y2)連同新點(diǎn)橫、縱坐%標(biāo)依次分別存儲(chǔ)在數(shù)組X,Y中function [X,Y]=levy1(x1,y1,x2,y2)x3=1/2*(x1+x2+y1y2)。y3=1/2*(x1+x2+y1+y2)。X=[x1,x3,x2]。Y=[y1,y3,y2]。%把由函數(shù)levy1生成的三點(diǎn)橫、縱坐標(biāo)X,Y順次劃分為兩組,分別對(duì)應(yīng)兩條折線段中每條線%段兩端點(diǎn)的坐標(biāo),并依次分別存儲(chǔ)在2*4階矩陣w中,w中第i(i=1,2)行數(shù)字代表第i條線段%兩端點(diǎn)的坐標(biāo)function w=levy2(X,Y)a11=X(1)。b11=Y(1)。a12=X(2)。b12=Y(2)。a21=X(2)。b21=Y(2)。a22=X(3)。b22=Y(3)。w=[a11,b11,a12,b12。a21,b21,a22,b22]。圖2 Levy 曲線（3） function tree(n,a,b)% tree(8,pi/8,pi/8),n為分形樹(shù)迭代次數(shù)%a,b為分枝與豎直方向夾角%x1,y1,x2,y2為初始線段兩端點(diǎn)坐標(biāo),nn為迭代次數(shù)n=8。a=pi/8。b=pi/8。x1=0。y1=0。x2=0。y2=1。plot([x1,x2],[y1,y2])hold on[X,Y]=tree1(x1,y1,x2,y2,a,b)。hold onW=tree2(X,Y)。w1=W(:,1:4)。w2=W(:,5:8)。% w為2^k*4維矩陣,存儲(chǔ)第k次迭代產(chǎn)生的分枝兩端點(diǎn)的坐標(biāo),% w的第i(i=1,2,…,2^k)行數(shù)字對(duì)應(yīng)第i個(gè)分枝兩端點(diǎn)的坐標(biāo)w=[w1。w2]。for k=1:n for i=1:2^k [X,Y]=tree1(w(i,1),w(i,2),w(i,3),w(i,4),a,b)。 W(i,:)=tree2(X,Y)。 end w1=W(:,1:4)。 w2=W(:,5:8)。 w=[w1。w2]。end%由每個(gè)分枝兩端點(diǎn)坐標(biāo)(x1,y1),(x2,y2)產(chǎn)生兩新點(diǎn)的坐標(biāo)(x3,y3),(x4,y4),畫兩分枝圖形,并把%(x2,y2)連同新點(diǎn)橫、縱坐標(biāo)分別存儲(chǔ)在數(shù)組X,Y中function [X,Y]=tree1(x1,y1,x2,y2,a,b)L=sqrt((x2x1)^2+(y2y1)^2)。if (x2x1)==0 a=pi/2。 else if
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1