复平面的迭代.pptVIP

参考书：《分形算法与程序设计》 复平面函数的迭代算法及作图 （1） 当初值︱z︱＜1 时，F(z)→0； （2） 当初值︱z︱＞1 时, F(z)→∞； 结论：实轴上有两个不动点，即0和∞，称为“平庸吸引子”。 （3） 当初值︱z︱= 1 时，F(z)在圆周∣Z∣=1上变化，呈现出不稳定的性态。 （1）.定义一个半径充分大的圆 确定以W为中心，R为半径的圆，R为正整数且足够大，使圆包含 W，再定义一个发散区域： V={(x,y)∈R2;(x2+y2)1/2R 当轨道进入V区域，轨道发散。 （ 2）. 确定最大迭代次数 （ 3）. 用不同的颜色绘制逃逸和不逃逸的点，也可以在逃逸的点中用确定该点是逃逸点的迭代数来确定绘制的颜色，这样可以观察到点逃逸的快慢。 Mandelbrot集 牛顿分形 Newton奠定了经典力学、光学和微积分学的基础。 但是除了创造这些自然科学的基础学科外，他还建立了一些方法，这些方法虽然比不上整个学科那么有名，但已被证明直到今天还是非常有价值的。 例如，牛顿建议用一个逼近方法求解一个方程的根。猜测一个初始点，然后使用函数的一阶导数，用切线逐渐逼近方程的根。如方程 Z^6 + 1 = 0有六个根，用牛顿的方法“猜测”复平面上各点最后趋向方程的那一个根，你就可以得到一个怪异的分形图形。和平常的Julia分形一样，你能永远放大下去，并有自相似性。 牛顿分形图形中的颜色显示每个答案的种类及性质，即迭代到目的地花费的时间， Newton分形 牛顿迭代法作图 对于复平面上面的任一点c=a+b*i，i为虚数单位，我们都可以得到一个复平面上的二次函数f(z)=z2+c，z为复变量。 用这个f(z)进行迭代，迭代初值为0，迭代次数为无穷次，如果最终的迭代所得到的值是一个有限值，那么点c就属于是Mandelbrot集,这些点c称为收敛点；若不能得到有限值，则点c不属于Mandelbrot集，称为发散点。 不用0作为初始迭代点，而用复平面上面的其它任意一点作为初始迭代点，再根据这些点的最终迭代结果给每个初始迭代点着色，我们就可以得到Julia集。 具体地说，Mandelbrot记录的是整个区域上的c值情况，而Julia集是取一固定的c 值后，观察复平面上每一点(x，y)在迭代中的表现，并把结果记录下来。 作图步凑： 1) 选择参数a or c、迭代次数L、收敛值M； 2) 选择象素点(x,y)； 3) 设置迭代计数i，初值为0； 4) 用公式Z=Za+C计算Z的新值； 5) 递增迭代次数i，即i=i+1； 6) 重复步骤5、6，一直到终结条件(i＜L或｜Z｜＜M)得到满足； 7) 在屏幕上显示此象素点(x,y)； 影响复迭代分形图的若干因素分析 Logistic方程的变形Z＝Z2+C对应着两类著名的分形集-Mandelbrot集(C平面分形)和Julia集(Z平面分形)。将Z=Z2+C中的指数2换成可变的参数α，即得到复迭代变换Z=Za+C ，其对应的分形图也分为两种情况：C平面分形和Z平面分形。这是两类不同的分形图，而它们又是有联系的。　　生成这两类分形图的一个通用过程如下： 1) 选择指数α、迭代次数L、收敛值M；2) 选择一象素点(x,y)；3) 根据所希望的分形形式选择C和Z(即生成C平面分形，还是Z平面分形)；4) 设置迭代计数i，初值为0；5) 用公式Z=Zα+C计算Z的新值； 6) 递增迭代次数i，即i=i+1；7) 重复步骤5、6，一直到终结条件(i＜L或｜Z｜＜M)得到满足；　　 8) 在屏幕上显示此象素点(x,y)；9) 重复步骤2～8，对所希望的所有象素点进行处理。　　 分析上述对应于复迭代变换Z=Zα+C的分形图生成过程，我们发现，参数C、指数α、迭代次数L、收敛值M的取值对最后生成的图形都是有影响的。 一.参数C的取值对Z平面分形的影响 研究此问题时假定α=2,则此时复迭代变换函数为 Z=Z2+C 现在考虑C=0和C≠0两种情况。 (1)　C=0　　复迭代变换简化为Z=Z2，根据初始点Z0的不同，有以下三种情况：　　.｜Z0｜＜1则｛Zn｝很快趋向于零，得到的图形只是一个点，称此点为吸附点(attractor)。　　.｜Z0｜＞1，则｛Zn｝→∞，Z0经过变换后无限发散，得不到图形。　　.｜Z0｜=1，则｛Zn｝永远保持在单位圆｜Z｜=1上，此时得到的图形是圆心在原点、半径为1的圆。 （2）C≠0　　如图1所示，不同的参数C，对应着形态各异的分形图象。而对于确定的参数C，其对应分形图的局部放大图与整体图像极其相似，因而具有严格的自似性。 不同参数C对应的Z平面分形图及其放大图，如下： 指