Mandelbrot闆嗗拰Julia闆嗙殑鍒嗗舰鍥句箣matlab瀹炵幇 - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

Mandelbrot

闆嗗拰

Julia

闆嗙殑鍒嗗舰鍥句箣

matlab

瀹炵幇

 

 

鍩轰簬閫冮�告椂闂寸畻娉?/p>

 

1.

 Mandelbrot

闆?/p>

 

function

 Mandelbrot(res,iter,xc,yc,xoom)

 

%Mandelbrot

 

% res

鏄�鐩�鏍囧垎杈ㄧ巼锛?/p>

iter

鏄�寰�鐜�娆℃暟锛岋�?/p>

xc,yc

锛夋槸鍥惧儚涓�蹇冿�?/p>

xoom

鏄�鏀惧ぇ鍊嶆暟

 

x0=xc-2/xoom;x1=xc+2/xoom;

 

y0=yc-2/xoom;y1=yc+2/xoom;

 

x=linspace(x0,x1,res);

 

y=linspace(y0,y1,res);

 

[xx,yy]=meshgrid(x,y);

 

z=xx+yy*1i;             

 

C=z;

 

N=zeros(res,res); 

%

鍒濆�嬪�?/p>

N

锛屾渶缁堟牴鎹?/p>

N

锛屽�瑰悇鐐硅繘琛屾煋鑹�

 

tic                  

%

 

鏄剧ず

tic

鍜?/p>

toc

闂寸殑绋嬪簭杩愯�屾椂闂�

  

 

for

 k=1:iter     

 

   z=z.^2+C;        

%

瀵圭┖闂翠笂姣忕偣閮借繘琛岃凯浠?/p>

 

   N(abs(z)>4)=k;  

%

閫冮�稿崐寰勪负

4

锛岃�烘煇鐐归�冮�革紝璁板綍閫冮�告椂闂?/p>

k,

鏈�閫冮�稿垯鏃堕棿涓?/p>

0

 

   z(abs(z)>4)=0;

 

   C(abs(z)>4)=0;

 

end

 

imshow(N,[]);   

 

toc

 

end 

 

>>Mandelbrot(512,100,0,0,1)