我用逃逸时间算法画的Julia集的彩色分形图

liliangbiao

详细请见孙博文的《分形算法与程序设计-Visual C++实现》。他是用Vc++实现的，我用Matlab实现的，但是我的图像和他得到的彩色不太一样，好像他的好些！

咕噜噜

能不能贴出程序来分享一下

无水1324

看了一下你的两个帖子，基本上都是图片欣赏类的，建议发到图片版？

liliangbiao

不是欣赏啊！我这里是讨论啊！因为我说了孙博文先生画的要比我做的好些，我现在向大家请教如何在Matlab里面实现他那种色彩特别逼真的分形图，我在Matlab里面调试了好久，也没有实现，不知道大家能不能帮这个忙？

咕噜噜

那你把源程序发上来啊，一者大家可能没时间找书再自己编程序，二者毕竟研究方向都不相同，不可能有那么充足的时间专门研究一下这个，还请见谅

liliangbiao

请浏览我在研学论坛 混沌与分形上边的发表的帖子吧！

咕噜噜

:@) 不好意思，我很少进研学
而且这里大多数会员可能也没有注册过那里，发到这里共享不是很好吗？为什么一定要去那里？

liliangbiao

呵呵 无所谓了，这些天发帖子都累的很了，算了吧，要是真正的想学一样东西，就不怕累吧！
所以我觉得去研学没有什么错误的。互相学习一下！

[ 本帖最后由 mjhzhjg 于 2007-8-26 15:46 编辑 ]

咕噜噜

:@) 去研学当然没错误，只是我觉得你发到这里让大家更方便共享而已
大家来这里当然想学东西，可是也要为方便大家考虑不是吗

hohoo

到研学论坛看了 没有找到你的程序呀

liliangbiao

有的！

hohoo

我在研学论坛上找了一下 只找到liliangbiao这个Mandelbrot Set集的分形程序 自己没有试 大家看看
%这是画的Mandelbrot Set集的分形程序(彩色图像)！经过试验和调试是没有错误的！
Nmax = 50; scale = 0.005;
xmin = -2.4; xmax = 1.2;
ymin = -1.5; ymax = 1.5;
% Generate X and Y coordinates and Z complex values
[x,y]=meshgrid(xmin:scale:xmax, ymin:scale:ymax);
z = x+1i*y;
% Generate w accumulation matrix and k counting matrix
w = zeros(size(z));
k = zeros(size(z));
% Start off with the first step ...
N = 0;
% While N is less than Nmax and any k's are left as 0
while N<Nmax & ~all(k()
% Square w, add z
w = w.^2+z;
% Increment iteration count
N = N+1;
% Any k locations for which abs(w)>4 at this iteration and no
% previous iteration get assigned the value of N
k(~k & abs(w)>4) = N;
end
% If any k's are equal to 0 (i.e. the corresponding w's never blew up) set
% them to the final iteration number
k(k==0) = Nmax;
% Open a new figure
figure
% Display the matrix as a surface
s=pcolor(x,y,k);
% If you truly want the Mandelbrot curve in B&W, comment the above line and
% uncomment these two
% s = pcolor(x, y, mod(k, 2));
% colormap([0 0 0;1 1 1])
% Turn off the edges of the surface (because the cells are so small, the
% edges would drown out any useful information if we left them black)
set(s,'edgecolor','none')
% Adjust axis limits, ticks, and tick labels
axis([xmin xmax -ymax ymax])
fontsize=15;
set(gca,'xtick',[xmin:0.4:xmax],'FontSize',fontsize)
set(gca,'ytick',[-ymax:0.5:ymax],'FontSize',fontsize)
xlabel('Re z','FontSize',fontsize)
ylabel('Im z','FontSize',fontsize)

yzsldj

这里有更简单的Mandelbrot Set集的分形程序：

x=linspace(-2,0.5,300);
y=linspace(-1.25,1.25,300);      
[X,Y]=meshgrid(x,y);
n_iter=200; c=0; a=2;           
W=mandelbrot(X,Y,c,a,n_iter);                     
pcolor(X,Y,W), shading flat;
axis('square'); colormap prism(256);


function W=mandelbrot(X,Y,c,a,n_iter)           
C=X+i*Y;
W=zeros(size(X));
Z=c;
for k=1:n_iter,
     Z=Z.^2+C;                                 
     i0=find(abs(Z)>a);
   W(i0)=k;
  Z(i0)=NaN;
end
i0=find(W==0); W(i0)=NaN;

zhangnan3509

这么做好像不怎么厚道吧

huang880116

有没有画julia集的程序啊？在matlab上实现