用Mathematica7绘制分形图Mandelbrot集合的程序问题,急!!

用Mathematica7绘制分形图Mandelbrot集合的程序问题,急!!,第1张

cx,cy的意思就是参数Z的实部和虚部的值。这个程序中令cx=1,cy=0,则迭代方程还是Zn+1 =Zn^k+C,n趋向正无穷时就有z=z^n+c了

这程序设c的模小于2,迭代次数小于50次

用了个Block函数将z作为局部分量和全局变量区分清楚

通过画密度函数画出c的图像

实际上,这个程序可简化成

M1[x_,y_]:=

Block[{z,k=0},

z=x+y*I

While[(Abs[z]<2.0)&&(k<50),++k

z=z^2+(x+y*I)]

Return[k]

]

M2[px_,py_,p_]:=

Block[{t},

t=DensityPlot[M1[xx,yy],{xx,-1.5,0.5},{yy,-1.2,1.2},p,

Mesh->False,ColorFunction->Hue]

Return[t]

]

Mandelbrot=M2[x,y,{PlotPoints->120,PlotLabel->"Mandelbrot集"}]

你看看行不行

这次我们用循环实现了一个Mandelbrot图形的计算。在《xialulee大战pythonchallenge—...下面是程序运行的结果: 本回答由网友推荐 评论 xiangjuan314 采纳率:94% ...


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/11084022.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-13
下一篇 2023-05-13

发表评论

登录后才能评论

评论列表(0条)

保存