分形幾何中一些經(jīng)典圖形的matlab畫法-免費(fèi)閱讀
【正文】 markersize39。 end end endx=a*x+b*y+e。d=M2(5)。e=M1(3)。% r為[0,1]區(qū)間內(nèi)產(chǎn)生的n維隨機(jī)數(shù)組r=rand(1,n)。w1=。 w(1,t)=i。,1) 一維元胞自動(dòng)機(jī)畫Sierpinski三角形178。s(i,j+1)==1) s(i+1,j)=1。t=1。.39。 y=2*y。w=1。,)圖5 Julia 集(6) function sierpinski(a,b,c,d,n,m,r)%sierpinski(0,0,1,1,12,200,200)%(a,b),(c,d)收斂區(qū)域左上角和右下角坐標(biāo),m為分辨率% n為逃逸時(shí)間,需要反復(fù)試探,r逃逸半徑a=0。 y=r*sin(alpha)。 if wx0 alpha=atan(wy/wx)。y=image(z0)。,39。f=M3(6)。c=M2(4)。b=M1(2)。y=0。圖3 分形樹（4）function sierpinski_ifs(n,w1,w2,w3)%sierpinski_ifs(10000,1/3,1/3,1/3)%w1,w2,w3出現(xiàn)頻率n=10000。b1=Y(1)。y3=y2+L*2/3*sin(a+b)。 w2=W(:,5:8)。hold onW=tree2(X,Y)。圖2 Levy 曲線（3） function tree(n,a,b)% tree(8,pi/8,pi/8),n為分形樹迭代次數(shù)%a,b為分枝與豎直方向夾角%x1,y1,x2,y2為初始線段兩端點(diǎn)坐標(biāo),nn為迭代次數(shù)n=8。a12=X(2)。 end X=XX。x1=0。b41=B(4)。b21=B(2)。dy=cy+sin(alpha+pi/3)*L。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。b2=0。n=3。cy=ay+(byay)/3。A=[ax,cx,dx,ex,bx]。a22=A(3)。a42=A(5)。y1=0。 Y=YY。b12=Y(2)。a=pi/8。w1=W(:,1:4)。 w=[w1。x4=x2+L*2/3*cos(ab)。a2=X(2)。w1=1/3。% r為[0,1]區(qū)間內(nèi)產(chǎn)生的n維隨機(jī)數(shù)組r=rand(1,n)。e=M1(3)。d=M2(5)。 end end endx=a*x+b*y+e。markersize39。x=1。 end if wx0 alpha=pi+atan(wy/wx)。 B(1,k)=x。b=0。for i=1:m x0=a+(ca)*(i1)/m。 else x=2*x。,39。w=zer
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號-1