电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

分形图之matlab实现VIP免费

分形图之matlab实现_第1页
1/7
分形图之matlab实现_第2页
2/7
分形图之matlab实现_第3页
3/7
分形图之matlab 实现 1. Mandelbrot集 function Mandelbrot(res,iter,xc,yc,xoom) %Mandelbrot % res是目标分辨率,iter是循环次数,(xc,yc)是图像中心,xoom是放大倍数 x0=xc-2/xoom;x1=xc+2/xoom; y0=yc-2/xoom;y1=yc+2/xoom; x=linspace(x0,x1,res); y=linspace(y0,y1,res); [xx,yy]=meshgrid(x,y); z=xx+yy*1i; C=z; N=zeros(res,res); %初始化N,最终根据N,对各点进行染色 tic % 显示tic和toc间的程序运行时间 for k=1:iter z=z.^2+C; %对空间上每点都进行迭代 N(abs(z)>4)=k; %逃逸半径为4,诺某点逃逸,记录逃逸时间k,未逃逸则时间为0 z(abs(z)>4)=0; C(abs(z)>4)=0; end imshow(N,[]); toc end >>Mandelbrot(512,100,0,0,1) > > Mandelbrot(512,128,-1.478,0,300) 2.Ju lia 集 function Julia(c,res,iter,xc,yc,xoom) %Julia集 %c为参数, res是目标分辨率,iter是循环次数,(xc,yc)是图像中心,xoom是放大倍数 x0=xc-2/xoom;x1=xc+2/xoom; y0=yc-2/xoom;y1=yc+2/xoom; x=linspace(x0,x1,res); y=linspace(y0,y1,res); [xx,yy]=meshgrid(x,y); z=xx+yy*1i; N=zeros(res,res); C=c*ones(res,res); for k=1:iter z=z.^2+C; N(abs(z)>2)=k; C(abs(z)>2)=0; z(abs(z)>2)=0; end colormap jet; image(x,y,N); axis square; end >>Ju lia(i,512,200,0,0,1) >> Ju lia(i,512,200,0,0,2000) >>Ju lia(0.1+0.7i,512,200,0,0,1) >>Ju lia(-0.8-0.21i,512,200,0,0,1)

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

分形图之matlab实现

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部