第六章复平面上的迭代-Read.docVIP

第四章 复平面上的迭代 4.1 复数四则运算与迭代 迪万内在《混沌动力系统引论》中说：“复解析函数动力学是一个比较专门的课题。……20 世纪20年代，这一领域在数学家法图和朱丽亚的指导下曾经极度繁荣，……但在70年代的新繁荣时期以前有一段休眠期。后来，主要由曼德布罗特的美丽动人的计算机作图和道阿弟、哈伯德和沙利文(D.Sullivan)等人引人入胜的数学工作，又一次把人们的注意力引向复平面丰富多彩的初等映射的动力学行为。” 对于本书作者和读者而言，“解析函数”、“复变函数”等概念都可暂且不管，知道在复平面上对简单的多项式进行迭代能够生成美妙的分形图形，就足够了。在计算机十分普及以前，要做这种练习几乎是不可能的。当年朱丽亚和法图是靠抽象的数学思维，艰难地研究复迭代的。进入80年代，计算机产业大发展，正巧分形几何学又到来了，一切条件都已成熟。曼德布罗特在干柴上点了一把火，复迭代复兴了，这把火一直烧到现在，而且火势仍不减当年。 复数迭代比实迭代(见第八章)稍稍麻烦一些，进行复迭代需要知道复数的四则运算。我们先复习一下高中学过的内容。 任意复数可以写成三种形式：1)代数形式：z=x+yi；2)三角形式：z=r (cosθ+isinθ)；3)指数形式：z=re^(iθ)。 其中i=SQRT(-1)是虚数单位，x=Re(z)是复数的实部，y=Im(z )是复数的虚部,r=｜z｜=｜x+yi｜=SQRT(x^2+y ^2)叫做复数的模，θ叫做复数的辐角。在复平面上，实部一般用横坐标表示， 虚部一般用纵坐标表示。 设z_1，z_2，z_3都是复数，复数基本运算如下： 恒等式：e^(iθ)=cosθ+isinθ, 加法：z_1+z_2=(x_1+x_2)+(y_1+ y_2)i, 减法：z_1-z_2=(x_1-x_2)+(y_1- y_2)i， 乘法：z_1·z_2=(x_1x_2- y_1y_2)+(x_1y_2+x_2y_1)i， 除法：z_1/z_2=(x_1x_2+y_1y_2)/ (x^2_2+y^2_2)+(x_2y_1-x_1y_2) / (x^2_2+y^2_2)i， 余弦：cos(x + iy)= cosxcoshy-isinx sinhy, 正弦：sin(x + iy)= sinxcoshy+icosx sinhy, 指数：e^(x + iy)= e^xcosy+ie^x siny, 其中 双曲余弦：coshx= [e^x + e^(-x)]/2, 双曲正弦：sinhx= [e^x - e^(-x)]/2. 在复迭代程序中，这些基本运算经常用到，最好将它们都编制成函数(过程)，存到一个单元中去，使用时可以把它们当成“标准函数(过程)”调用。设单元为COMPLEX.TPU， COMPLEX.PAS的内容如下： unit complex;｛PASCAL单元：复数四则运算｝ interface procedure add (x1, y1, x2, y2 : real; var x3, y3 : real); {复数加法} procedure sub (x1, y1, x2, y2 : real; var x3, y3 : real); {复数减法} procedure mult(x1, y1, x2, y2 : real; var x3, y3 : real); {复数乘法} procedure cdiv(x1, y1, x2, y2 : real; var x3, y3 : real); {复数除法} function cosh(x : real) : real; {双曲余弦函数} function sinh(x : real) : real; {双曲正弦函数} procedure csin(x,y : real; var x1, y1 :real);{复数正弦} procedure ccos(x,y : real; var x1, y1 :real);{复数余弦} procedure cexp(x,y : real; var x1, y1 :real);{复数指数} implementation procedure add (x1, y1, x2, y2 : real; var x3, y3 : real);｛复数加法｝ {计算z_3 = z_1 + z_2， 其中:z_1 = x_1 + iy_1,z_2 = x_2 + iy_2,z_3 = x_3 + iy_3} begin x3 := x1+x2;y3 := y1+y2 end; procedure sub (x1, y1, x2, y2 : real;