原文<a href="http://math.ucr.edu/home/baez/roots/" target="_blank">http://math.ucr.edu/home/baez/roots/</a><br/>我试着画了一个龙的图案，matlab程序是<br/>clc;<br/>clear all;<br/>z=0.372-0.547i; <br/>k=0;<br/>fz=0;<br/>for a=1:10^4<br/> k=k+1;<br/> for b=1:50 %一般50次迭代就收敛了<br/> if mod(b,round(rand(1)))<1 %这里的转换规则是我自己瞎写的，不同的规则对应龙中不同的点<br/> fz=1+z*fz;<br/> else<br/> fz=1-z*fz;<br/> end<br/> end<br/> dragon(k)=fz;<br/>end<br/>plot(dragon,'r.');<br/>axis equal;<br/>很白痴，但还好，跟大哥画的一样。<br/>把z改成0.8+0.1i就是羽毛形状了。<br/>绝对最终形状的就是初值z，大家可以试试不同的z得到什么结果
↧