thorn_fractal.m

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

function Thorn_Fractal; 
% x=x./cos(y)+c_x 
% y=y./cos(x)+c_y; 
% Author's email:[e-mail] 
% reference: 
% http://astronomy.swin.edu.au/~pbourke/fractals/thorn/ 
warning off; 
L=2; 
th=3*L; 
[x,y]=meshgrid(linspace(-L,L,300)); 
N=100; 
c_x=1*ones(size(x)); 
c_y=0.2*ones(size(x)); 
D=N*ones(size(x)); 
for k=1:N; 
   x1=x./cos(y)+c_x; 
   y1=y./cos(x)+c_y; 
   D(abs(x1+i*y1)>th)=k; 
   c_x(abs(x1+i*y1)>th)=0; 
   c_y(abs(x1+i*y1)>th)=0; 
   x=x1;y=y1; 
   x(abs(x1+i*y1)>th)=0; 
   y(abs(x1+i*y1)>th)=0; 
end 
figure; 
imagesc(D);