MATLAB程式应用.DOC

MATLAB程式應用 碎形(Fractals) 班級：二技機電三甲 姓名：謝岳廷 學號指導老師：趙春棠 老師 碎行簡介： ◎何謂碎形？? ? 所謂碎形指的是：一個幾何物件，在越來越細微的尺度上，不斷重複自我模仿，以固定的方式重複製造細節。如果以高倍顯微鏡放大局部，看來和全體依舊相同。其外形不像直線、弧線和古典歐氏幾何的表面那麼的平滑，而是相當破碎、不規則、且不論巨觀或微觀都不連續，而是跳動的驚人的擬似? ? 古典數學家認為碎形圖形(當然當時並沒有碎形這名稱)面目可憎，稱之為「病態的形狀」，根本不屑一顧。然而，曼德布洛特對其情有獨鍾，他在擔任IBM研究員期間，嘗試用電腦產生碎形圖形。經由合適的方程式及演算法則，當曼氏的處女作人工海岸線出現在眾人面前時，大家不禁大吃一驚那不就是紐西蘭嗎？於是人們不得不同意，碎形確實是自然界中的一份子。 曼德布洛特何許人也？? ? 面對自然的複雜性，我們不禁要尋思：這複雜性是隨機的意外？還是有跡可循的？? ? 法國數學家曼德布洛特(Benoit Mandelbrot, 1924- )獨創了數學的一個全新支派碎形論(fractals)。碎形研究的開展，告訴我們任何奇怪的形狀也有其意義；這些凹摺和糾結，絕非古典幾何形狀的瑕疵，它們常是開啟物質本質的鑰匙。數學家們預言：碎形論將成為未來數學的新主流。 關於Manderbrot Set的性質 程式碼 %Mandelbrot.mclear all;close all;Xlower=-2;Xupper=2;Ylower=-2;Yupper=2;dpi=400;precision=50;Re=ones(dpi,1)*linspace(Xlower,Xupper,dpi);Im=(linspace(Yupper,Ylower,dpi))*ones(1,dpi);A=Re+i*Im;hold on;for p=1:dpi??? for q=1:dpi?????? Z=A(p,q);?????? counter=0;?????? while(counterprecision)????????? if abs(Z)2???????????? break????????? end????????? Z=Z^2+A(p,q);????????? counter=counter+1;?????? end?????? if counter==precision????????? plot(real(A(p,q)),imag(A(p,q)),b)?????? else????????? switch(mod(counter,4))????????? case 0???????????? plot(real(A(p,q)),imag(A(p,q)),c)????????? case 1???????????? plot(real(A(p,q)),imag(A(p,q)),g)????????? case 2???????????? plot(real(A(p,q)),imag(A(p,q)),y)????????? case 3???????????? plot(real(A(p,q)),imag(A(p,q)),m)????????? end?????? end??? endendaxis equal;axis ([Xlower Xupper Ylower Yupper]);hold off; 程式說明： 程式碼中的變數的介紹： Xlower,Xupper,Ylower,Yupper－圖形X軸與Y軸的上界與下界。 dpi－解析度，將圖形範圍的長與寬各分割成dpi*dpi的點陣列進行計算。就我的PC而言，此數值設定在400以下較恰當，1分鐘內就可以跑出圖形。過多的dpi會記憶體不足，而且跑了老半天沒有結果。 Re,Im,A－矩陣A就是記錄了點陣列的dpi*dpi的複數矩陣，先分別算出實部Re與虛部Im再合併而得。precision,counter－precision就是迭代次數限制，計算時從矩陣A逐一取出點 進行迭代，以counter紀錄該點被迭代的次數。迭代出現模長大於2者判定發散，則counter必小於precision。達到迭代次數限制未發散即判定收斂，則counter等於precision。precision建議是20以上，太小的話圖形近似度太低，用50。進一步來說，圖形的區域較大時，precision取太大沒有意義，因為細微的區域本來就不清楚﹔圖形區域很小時，precision應該取很大，否則圖形的正確性不夠。將會呈現這個概念。plot(real(A(p,q