分形算法程序.pdfVIP

Julia 集的算法程序 c=-0.4+0.5i; [x,y]=meshgrid(linspace(-2,2,200)); z=x+y*i; N=zeros(size(x)); C=c*ones(size(x)); for k=1:50; z=z.^2+C; N(abs(z)2)=k; C(abs(z)2)=0; z(abs(z)2)=0; end imshow(N,[]) M 集算法程序 function liusir(pmin,pmax,qmin,qmax) % 逃逸时间算法 % Examples: % liusir(-2,2,-2,2) dp=(pmax-pmin)/200; dq=(qmax-qmin)/200; [Xp,Yp]=meshgrid(pmin:dp:pmax,qmin:dq:qmax); X=zeros(size(Xp)); Y=X; Z=ones(size(X)); N=Z; for n=1:300; x1=X.*X-Y.*Y+Xp; y1=2*X.*Y+Yp; r=x1.^2+y1.^2; Z(r16)=0; N(r16)=n; X=x1.*Z; Y=y1.*Z; end figure; surf(Xp,Yp,round(N/20)) shading interp colormap(flipud(jet)) view(2) grid off 另一种算法程序 p=[-2.25 .75];q=[-1.5 1.5];a=200; dp=(p(2)-p(1))/a; dq=(q(2)-q(1))/a; [X,Y]=meshgrid(p(1):dp:p(2),q(1):dq:q(2)); P0=X+i*Y;x=zeros(size(P0)); for j=1:300 x=x.*x+P0; modx=abs(x); x(modx100)=0; P0(modx100)=0; end plot(P0(P0~=0),.) sierpinski 毯的程序 function Sierpinski_rug(n) % 生成Sierpinski 地毯 % Examples: % Sierpinski_rug(7) global A_d global Rm global Lm % axis off; Lm=+R-L-R+; Rm=-L+R+L-; x0=0; y0=0; anglez0=0; A_d=pi/3; ch=A_Iterate(n); run_order(x0,y0,anglez0,ch); % clear; %***************************************************************************** * function ch=A_Iterate(n); if n==0; ch=L; elseif n0; ch=a_Ite(A_Iterate(n-1)); end %***************************************************************************** function cha=a_Ite(ch) global Rm global Lm cha=[]; for s=1:length(ch); if ch(s)==R; cha=strcat(cha,Rm); elseif ch(s)==L; cha=strcat(cha,Lm); elseif ch(s)==+ | -; cha=strcat(cha,ch(s)); end end %***************************************************************************** * function run_order(x0,y0,anglez0,ch) global A_d for s=1:length(ch); if ch(s)==+; anglez0=anglez0+A_d; % + elseif ch(s)==-; anglez0=anglez0-A_d; % - elseif ch(s)==R | L; x=x0+cos(anglez0