iterative_z25cd.m

来自「分形计算的几个程序原码」· M 代码 · 共 22 行

%映射:z--z^2.5+c 的动态(映射后的分形图是动态显示的)的MATLAB程序 
function iterative_z25cd; 
% z--->z^2.5+c; 射后的分形图是动态显示的 
N=256; 
[x,y]=meshgrid(linspace(-2,2,N));
c=x+i*y;
z=zeros(N); 
N=zeros(N);
figure;
set(gcf,'DoubleBuffer','on'); 
h=imshow(N);
for k=1:80; 
    z=z.^2.5+c; 
    r=abs(z);
    z(r>4)=0;
    c(r>4)=0;
    N(r>4)=k; 
    set(h,'CData',N); 
    pause(0.5); 
end