分形幾何中一些經(jīng)典圖形的matlab畫法-文庫(kù)吧資料
2025-04-22 23:40本頁(yè)面
【正文】 if wx0 alpha=pi+atan(wy/wx)。 wy=ycy。% A為產(chǎn)生的服從標(biāo)準(zhǔn)正態(tài)分布的n維隨機(jī)數(shù)組A=randn(1,n)。 B=zeros(2,n)。x=1。% z^2+c=z0,x=real(z0)。cx=。cy=image(c)。markersize39。.39。k=k+1。B(1,k)=x。 end end endx=a*x+b*y+e。d=M3(5)。e=M3(3)。 else if r(i)w1+w2+w3 a=M3(1)。d=M2(5)。e=M2(3)。 else if r(i)w1+w2 a=M2(1)。d=M1(5)。e=M1(3)。for i=1:n if r(i)w1 a=M1(1)。% 當(dāng)1/3=r(i)2/3時(shí),進(jìn)行M2對(duì)應(yīng)的壓縮映射。k=1。% r為[0,1]區(qū)間內(nèi)產(chǎn)生的n維隨機(jī)數(shù)組r=rand(1,n)。x=0。M2=[ 0 0 0]。w3=1/3。w1=1/3。w=[a1,b1,a2,b2,a3,b3,a4,b4]。a4=X(3)。a3=X(1)。a2=X(2)。%把由函數(shù)tree1生成的X,Y順次劃分為兩組,分別對(duì)應(yīng)兩分枝兩個(gè)端點(diǎn)的坐標(biāo),并存儲(chǔ)在一維%數(shù)組w中function w=tree2(X,Y)a1=X(1)。plot(a,b)axis equalhold onX=[x2,x3,x4]。a=[x3,x2,x4]。x4=x2+L*2/3*cos(ab)。 endendx3=x2+L*2/3*cos(a+b)。 else if (x2x1)0 a=pi+atan((y2y1)/(x2x1))。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)。 w=[w1。 end w1=W(:,1:4)。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為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。w1=W(:,1:4)。plot([x1,x2],[y1,y2])hold on[X,Y]=tree1(x1,y1,x2,y2,a,b)。x2=0。x1=0。a=pi/8。a21,b21,a22,b22]。b22=Y(3)。b21=Y(2)。b12=Y(2)。b11=Y(1)。Y=[y1,y3,y2]。y3=1/2*(x1+x2+y1+y2)。 Y=YY。 [XX(3*2*(j1)+3+1:3*2*(j1)+3+3),YY(3*2*(j1)+3+1:3*2*(j
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1