分形几何中一些经典图形的Matlab画法.docxVIP

分形几何中一些经典图形的 Ma lab 画法 ———————————————————————————————— 作者 : ———————————————————————————————— 日期 : 分形几何中一些经典图形的Ｍ atlab 画法 （ 1）Ｋ och  曲线程序  koｃ h.m function ｋ oｃh(a１ ,b1,a2,b２， n) %ｋ ocｈ（ 0,0,９,0,3） %a1，b1,a２ ,b２为初始线段两端点坐标 a1=0;ｂ１ =0; ａ 2＝ 9;ｂ 2=0;ｎ＝３ ;  ,ｎ为迭代次数 %第ｉ－ 1 次迭代时由各条线段产生的新四条线段的五点横、纵坐标存储在数组 [ Ａ ,B］ =sｕｂ＿ koch1（ a１ ,b1,ａ 2,b2)；  A、B 中 ｆｏ r  i=1:n foｒ j=1:le ｎ gth(A)/5; w＝ sub＿ koch2( Ａ (1+5*( ｊ -１）： 5＊ j),Ｂ (1+ ５＊（ j －１）：５＊ｊ )); for k=1: ４ [AA （ 5*4＊ (j- １）＋ 5*(k-1)+ １： 5*4 ＊ (j-1)+5*( ｋ -1)+5),BB(5 ＊ 4* （ j－ 1)＋５ *(k － 1)+1 ：５＊ 4*( ｊ -1)+5 ＊ (k－ 1） +５） ] ＝ sub_koch1（ w(k,1 ） ,ｗ (k，２ )，ｗ (ｋ ,３ ),ｗ k,4)); eｎ d end A=ＡA； B＝BB; nd plｏｔ (A,B) hoｌｄ on aｘis ｅ qｕ al %由以 (aｘ ,aｙ ),(bx,b ｙ )为端点的线段生成新的中间三点坐标并把这五点横、纵坐标依次分别存％储在数组 A,B 中 unｃ tion [A,B]=s ｕb_koch １ (ax,ay,ｂ x,by) cx=ax+(bx- ａ x)／ 3; cy=ay+ （ｂ y-ay)/3; x=bx － (ｂ x-ax)/3; ey＝ bｙ -(by- ａ y)/3 ； =ｓ qrt((ex-cx). ＾ 2+(ey-cy).^2 ） ; aｌpha=atan（ (ｅｙ－ cy). ／(eｘ－ cx)); i ｆ (eｘ -cx) ０ ａｌ pha=alｐha＋ pi; end dx=cx+co ｓ (alpｈａ +ｐｉ／ 3）* Ｌ； dy=ｃ y+sin （ aｌ pha+ｐ i/3)* Ｌ ; A=[a ｘ ,cx,dx,ex,bx]; B＝ [ay,ｃ y,dy,ey,bｙ］ ; %把由函数 suｂ _koc ｈ１生成的五点横、纵坐标 A,B 顺次划分为四组，分别对应四条折线段中 %每条线段两端点的坐标，并依次分别存储在  4*4  阶矩阵  k 中 ,k 中第  i(i= １， 2,3,4）行数字 代表第 %i 条线段两端点的坐标 ｆｕ ncｔ iｏｎ w= ｓ ub_ｋ och2(Ａ ,B) a11=A（１ )； b11=B（ 1)； ａ1２＝Ａ（ 2)； b１ 2＝B(2); a21=A( ２ )； b2１ =B(2 ）； 22=A （ 3);b22=B(3); a31＝ A （3);b ３1=B （ 3); a3２＝ A（ 4);b ３２ =B(4 ）； a4１ =Ａ (4);b ４1=B(4) ； a42＝ A( ５ );b4２ =B(5); w=[ ａ１１ ,b11,a1２ ,ｂ 1２； a21，b21，a２ 2,ｂ２ 2;a３ 1，b3１， a３２ ,b3２ ; ａ41,ｂ 41，a42, ｂ  4 ２  ］  ; 图 1 Ｖ on Ｋｏ ch 曲线 （ 2） Levy  曲线程序  l ｅｖ y.m function levy （ n) % levy(16),n 为 lｅ vｙ曲线迭代次数 %x1, ｙ 1,x2,y2 为初始线段两端点坐标  ,ｎ n 为迭代次数 n＝１ 6; x1=0;y1=0; x2＝ 1； y２ =0; %第ｉ -１次迭代时由各条线段产生的新两条线段的三端点横、纵坐标存储在数组Ｘ、 Y 中 [X ， Y］ =lev ｙ 1(ｘ１ ,y1,x2,y2) ； f ｏ r i= １ :n for  ｊ＝ 1:ｌ enｇtｈ（Ｘ )／ 3 ｗ =ｌ eｖ y2(X （ 1+3＊ (j － 1):３＊ j),Ｙ (1+3*(j-1):3*j) ） ; ［ X Ｘ (3＊ 2* （ｊ -1） +1:3*2 ＊ (j － 1)+ ３ ),Ｙ Y(3*2*(j-1  ） +１ :3* ２ *(j-1  ）+3) ］ =le ｖ y1(w(1, １ ),w(1,2) ，ｗ (1,3),w(1,4)) ； [XX(3*2*(j-