分形幾何中一些經(jīng)典圖形的matlab畫(huà)法(已修改)
2025-04-28 23:40 本頁(yè)面
【正文】 分形幾何中一些經(jīng)典圖形的Matlab畫(huà)法（1）function koch(a1,b1,a2,b2,n)%koch(0,0,9,0,3)%a1,b1,a2,b2為初始線段兩端點(diǎn)坐標(biāo),n為迭代次數(shù)a1=0。b1=0。a2=9。b2=0。n=3。%第i1次迭代時(shí)由各條線段產(chǎn)生的新四條線段的五點(diǎn)橫、縱坐標(biāo)存儲(chǔ)在數(shù)組A、B中[A,B]=sub_koch1(a1,b1,a2,b2)。for i=1:n for j=1:length(A)/5。 w=sub_koch2(A(1+5*(j1):5*j),B(1+5*(j1):5*j))。 for k=1:4 [AA(5*4*(j1)+5*(k1)+1:5*4*(j1)+5*(k1)+5),BB(5*4*(j1)+5*(k1)+1:5*4*(j1)+5*(k1)+5)]=sub_koch1(w(k,1),w(k,2),w(k,3),w(k,4))。 end end A=AA。 B=BB。endplot(A,B)hold onaxis equal%由以(ax,ay),(bx,by)為端點(diǎn)的線段生成新的中間三點(diǎn)坐標(biāo)并把這五點(diǎn)橫、縱坐標(biāo)依次分別存%儲(chǔ)在數(shù)組A,B中function [A,B]=sub_koch1(ax,ay,bx,by)cx=ax+(bxax)/3。cy=ay+(byay)/3。ex=bx(bxax)/3。ey=by(byay)/3。L=sqrt((excx).^2+(eycy).^2)。alpha=atan((eycy)./(excx))。if (excx)0 alpha=alpha+pi。enddx=cx+cos(alpha+pi/3)*L。dy=cy+sin(alpha+pi/3)*L。A=[ax,cx,dx,ex,bx]。B=[ay,cy,dy,ey,by]。%把由函數(shù)sub_koch1生成的五點(diǎn)橫、縱坐標(biāo)A,B順次劃分為四組,分別對(duì)應(yīng)四條折線段中%每條線段兩端點(diǎn)的坐標(biāo),并依次分別存儲(chǔ)在4*4階矩陣k中,k中第i(i=1,2,3,4)行數(shù)字代表第%i條線段兩端點(diǎn)的坐標(biāo)function w=sub_koch2(A,B)a11=A(1)。b11=B(1)。a12=A(2)。b12=B(2)。a21=A(2)。b21=B(2)。a22=A(3)。b22=B(3)。a31=A(3)。b31=B(3)。a32=A(4)。b32=B(4)。a41=A(4)。b41=B(4)。a42=A(5)。b42=B(5)。w=[a11,b11,a12,b12。a21,b21,a22,b22。a31,b31,a32,b32。a41,b41,a42,b42]。圖1 Von Koch曲線（2）Levy function levy(n)% levy(16),n為levy曲線迭代次數(shù)%x1,y1,x2,y2為初始線段兩端點(diǎn)坐標(biāo),nn為迭代次數(shù)n=16。x1=0。y1=0。x2=1。y2=0。%第i1次迭代時(shí)由各條線段產(chǎn)生的新兩條線段的三端點(diǎn)橫、縱坐標(biāo)存儲(chǔ)在數(shù)組X、Y中[X,Y]=levy1(x1,y1,x2,y2)。for i=1:n for j=1:length(X)/3
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1