作业帮用Mathematica7绘制分形图Mandelbrot集合的程序问题,这是可行的程序Fx1[x_,y_,cx_,cy_,n_] :=Block[{z,ct = 0},z = x + y*I;While[(Abs[z] < 2.0) && (ct < 50),++ct;z = z^n + (cx + cy*I)*(x + y*I)];Return[ct];]Ht1[cx_,cy_,n_,pu_List,po_List
来源：学生作业帮助网 编辑：作业帮 时间：2024/11/05 22:00:34
![用Mathematica7绘制分形图Mandelbrot集合的程序问题,这是可行的程序Fx1[x_,y_,cx_,cy_,n_] :=Block[{z,ct = 0},z = x + y*I;While[(Abs[z] < 2.0) && (ct < 50),++ct;z = z^n + (cx + cy*I)*(x + y*I)];Return[ct];]Ht1[cx_,cy_,n_,pu_List,po_List](/uploads/image/z/5673016-64-6.jpg?t=%E7%94%A8Mathematica7%E7%BB%98%E5%88%B6%E5%88%86%E5%BD%A2%E5%9B%BEMandelbrot%E9%9B%86%E5%90%88%E7%9A%84%E7%A8%8B%E5%BA%8F%E9%97%AE%E9%A2%98%2C%E8%BF%99%E6%98%AF%E5%8F%AF%E8%A1%8C%E7%9A%84%E7%A8%8B%E5%BA%8FFx1%5Bx_%2Cy_%2Ccx_%2Ccy_%2Cn_%5D+%3A%3DBlock%5B%7Bz%2Cct+%3D+0%7D%2Cz+%3D+x+%2B+y%2AI%3BWhile%5B%28Abs%5Bz%5D+%3C+2.0%29+%26%26+%28ct+%3C+50%29%2C%2B%2Bct%3Bz+%3D+z%5En+%2B+%28cx+%2B+cy%2AI%29%2A%28x+%2B+y%2AI%29%5D%3BReturn%5Bct%5D%3B%5DHt1%5Bcx_%2Ccy_%2Cn_%2Cpu_List%2Cpo_List)
用Mathematica7绘制分形图Mandelbrot集合的程序问题,这是可行的程序Fx1[x_,y_,cx_,cy_,n_] :=Block[{z,ct = 0},z = x + y*I;While[(Abs[z] < 2.0) && (ct < 50),++ct;z = z^n + (cx + cy*I)*(x + y*I)];Return[ct];]Ht1[cx_,cy_,n_,pu_List,po_List
用Mathematica7绘制分形图Mandelbrot集合的程序问题,
这是可行的程序
Fx1[x_,y_,cx_,cy_,n_] :=
Block[{z,ct = 0},
z = x + y*I;
While[(Abs[z] < 2.0) && (ct < 50),++ct;
z = z^n + (cx + cy*I)*(x + y*I)];
Return[ct];
]
Ht1[cx_,cy_,n_,pu_List,po_List,pl_List] :=
Block[{kok},
kok = DensityPlot[
Fx1[xx,yy,cx,cy,n],{xx,pu[[2]],pu[[3]]},{yy,po[[2]],
po[[3]]},pl,
Mesh -> False,ColorFunction -> Hue];
Return[kok]
]
M1 = Ht1[1,0,2,{x,-1.5,0.5},{y,-1.2,1.2},{PlotPoints -> 120,
PlotLabel -> "X"}]
另外这个程序的思路
1.选定参数Z=cx+cy*I,取方程为Zn+1 =Zn^k+C*Z进行迭代.
2.Z的模小于2,迭代次数不超过50.
3.对于Z在平面上表示时,设xx为x的初始迭代坐标,yy为y的初始迭代坐标坐标表示.
4.所取的变量范围为（cx,cy）,分别为x与y的范围,迭代次数为n.
我对上面的有些疑问,首先M集合的迭代方程应该是Zn+1 =Zn^k+C,思路和程序中为什么用Z=Z^n+c*Z迭代还能对呢?
然后思路4,所取的变量范围为（cx,复数C应该是任意变化的,为什么有范围呢?如果不是变化范围的意思,那么cx,最后赋初值时cx=1,cy=0代表着什么?
用Mathematica7绘制分形图Mandelbrot集合的程序问题,这是可行的程序Fx1[x_,y_,cx_,cy_,n_] :=Block[{z,ct = 0},z = x + y*I;While[(Abs[z] < 2.0) && (ct < 50),++ct;z = z^n + (cx + cy*I)*(x + y*I)];Return[ct];]Ht1[cx_,cy_,n_,pu_List,po_List
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]Hue];
Return[t]
]
Mandelbrot=M2[x,y,{PlotPoints->120,PlotLabel->"Mandelbrot集"}]
你看看行不行