曼德尔布罗分形—碗形程序分享

scott198510

1. function bowl(arg1,arg2)
   
2. %
   
3. %~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   
4. % 本程序采用Mandelbrot Fractals （曼德尔布罗分形）生成三维碗形状
   
5. % 关于分形几何的细节可参阅：曼德尔布罗经典著作《大自然的分形几何》
   
6. % 输入参数优选：argu1>=40;argu2>=300;
   
7. arg1=40 ;
   
8. arg2=300;
   
9. dist1=-0.25;
   
10. dist2=0;
    
11. l=2.5;
    
12. xaxis=linspace(dist1-sin(l)*10,dist1+sin(l)*10,arg2);
    
13. yaxis=linspace(dist1-sin(l)*10,dist1+sin(l)*10,arg2);
    
14. [xtrans,ytrans]=meshgrid(xaxis,yaxis);
    
15. t=zeros(arg2);
    
16. c=xtrans+i*ytrans;
    
17. for k=1:arg1;
    
18. t=t.^2+log(c+101);
    
19. tem=-abs(t);
    
20. a = ( tem <-.5);
    
21. b = 1-a;
    
22. u = abs(c.*c);
    
23. d = b.*sin(-0.2*u);
    
24. end
    
25. e = d-1;
    
26. colormap(jet);
    
27. mesh(0.01*a+e);
    
28. set(gcf,'color','w');

复制代码

% 曼德尔布罗分形生成三维碗形状

scott198510

1. function density(arg1,arg2)
   
2. %
   
3. %~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   
4. %     本程序采用Mandelbrot Fractals （曼德尔布罗分形）表示密度的演化模型
   
5. %      输入参数优选：argu1>=40;argu2>=300;
   
6. 
   
7. arg1=41;
   
8. arg2=301;
   
9. dist1=-0.25;
   
10. dist2=0;
    
11. l=2.5;
    
12. xaxis=linspace(dist1-sin(l)*10,dist1+sin(l)*10,arg2);
    
13. yaxis=linspace(dist1-sin(l)*10,dist1+sin(l)*10,arg2);
    
14. [xtrans,ytrans]=meshgrid(xaxis,yaxis);
    
15. t=zeros(arg2);
    
16. c=xtrans+i*ytrans;
    
17. for k=1:arg1;
    
18. t=t.^+c*(c+101);
    
19. tem=-abs(t);
    
20. a = ( tem <-.5);
    
21. b = 1-a;
    
22. u = abs(c.*c);
    
23. d = b.*exp(-0.2*sin(u));
    
24. end
    
25. e = d-1;
    
26. colormap(jet);
    
27. mesh(0.01*a+e);
    
28. set(gcf,'color','w');

复制代码

再来一个，密度演化模型

电脑克星LK

chyanog 发表于 2013-10-5 16:42
bowl的运行结果不对啊，出来的是一个平面

我运行的碗程序挺好的，结果跟楼主给的图一样。

552006599

好牛，佩服，佩服