Matlab在复变函数中应用.doc

Matlab在复变函数中应用数学实验（一）西安交通大学理学院 二八年十一月MATLAB在复变函数中的应用复变函数的运算是实变函数运算的一种延伸，但由于其自身的一些特殊的性质而显得不同，特别是当它引进了“留数”的概念，且在引入了Taylor级数展开Laplace变换和Fourier变换之后而使其显得更为重要了。使用MATLAB来进行复变函数的各种运算；介绍留数的概念及MATLAB的实现；介绍在复变函数中有重要应用的Taylor展开和Laurent展开；利用Matlab 实现Laplace变换和Fourier变换。1 复数和复矩阵的生成在MATLAB中，复数单位为，其值在工作空间中都显示为。1.1 复数的生成复数可由语句生成，也可简写成。另一种生成复数的语句是，也可简写成，其中theta为复数辐角的弧度值，r为复数的模。1.2 创建复矩阵创建复矩阵的方法有两种。（1）如同一般的矩阵一样以前面介绍的几种方式输入矩阵例如：（2）可将实、虚矩阵分开创建，再写成和的形式例如：；注意 实、虚矩阵应大小相同。2 复数的运算1复数的实部和虚部复数的实部和虚部的提取可由函数real和imag实现。调用形式返回复数的实部返回复数的虚部2共轭复数复数的共轭可由函数conj实现。调用形式返回复数的共轭复数3复数的模和辐角复数的模和辐角的求解由功能函数abs和angle实现。调用形式复数的模复数的辐角例：求下列复数的实部与虚部、共轭复数、模与辐角（1）（2）（3）（4）由MATLAB输入如下： %实部0.23081.50003.50001.0000 %虚部0.15382.500013.00003.0000 %共轭复数0.2308+0.1538i 1.5000+2.5000i 3.5000+13.0000i1.0000+3.0000i%模0.27742.915513.46293.1623%辐角0.58801.03041.8228-1.24904复数的乘除法复数的乘除法运算由“/”和“”实现。例 复数的乘除法演示。由此例可见，相当于，和不相等。5复数的平方根复灵敏的平方根运算由函数sqrt实现。调用形式sqrt(x)返回复数的平方根值6复数的幂运算复数的幂运算的形式为，结果返回复数的次幂。例 求下列各式的值0.8660+0.5000 i7复数的指数和对数运算复数的指数和对数运算分别由函数exp和log实现。调用形式返回复数x的以e为底的指数值返回复数x的以e为底的对数值例 求下列式的值（参见参考资料【4】P.68.215）。8复数的三角函数运算复数的三角函数运算函数参见下面的复数三角函数复数三角函数表函数名函 数 功 能函数名函 数 功 能返回复数的正弦函数值返回复数的反正弦值返回复数的余弦函数值返回复数的反余弦值返回复数的正切函数值返回复数的反正切值返回复数的余切函数值返回复数的反余切值返回复数的正割函数值返回复数的反正割值返回复数的余割函数值返回复数的反余割值返回复数的双曲正弦值返回复数的双曲余切值返回复数的双曲余弦值返回复数的双曲正割值返回复数的双曲正切值返回复数的双曲余割值9 复数方程求根复数方程求根或实方程的复数根求解也由函数solve实现。见下面的例子.例 求方程所有的根（参见参考资料【4】P.32.116）。 23 留数留数定义：设a是的孤立奇点，C是a的充分小看邻域内一条把a点包含在其内部的闭路，积分称为在a点的留数或残数，记作。在MATLAB中，可由函数residue实现。residue 留数函数（部分分式展开）函数返回留数，极点和2个多项式比值的部分分式展开的直接项。如果没有重根，则向量B和A为分子、分母以s降幂排列的多项式系数，留数返回为向量R、极点在向量P的位置，直接项返回到向量K。极点的数目。如果，则直接项系数为空；否则。如果存在M重极点即有则展开项包括以下形式有3个输入变量和2个输出变量，函数转换部分因式展开还为系数为B和A的多项式比的形式。注意：数值上讲，分式多项式的部分因式展开实际上代表了一类病态问题。如果分母多项式是一个近似有重根的多项式，则在数值上的一点微小变化，包括舍入误差都可能造成极点和留数结果上的巨大变化。因此使用状态空间和零点极点表述的方法是可取的。例 求如下函数的奇点处的留数。在MATLAB实现如下1.50000.500020 所以可得。例 计算下面的积分其中C为正向圆周。（参见参考资料【4】P.158.例2）解：先求被积函数的留数0.25000.25000.25000.0000 i0.250+0.0000 i1.00001.00000.0000+1.0000 i0.00001.0000 i 可见在圆周内有四个极点，所以积分值等于。4 Taylor级数展开Taylor级数开展在复变函数中有很重要的地位，如分析复变函数的解析性等。函数在点的Taylor级数开展为在MATLAB中可由函数taylor来实现。taylor泰勒级数展开返回函数的五次幂多项式近似。此功能函数可有3个附加参数。返回次幂多项式。返回点附近的幂多项式近似。使用独立变量代替函数。例 求下列函数在指定点的泰勒开展式（参见参考资料【4】P.143.12）。（1）（2）；MATLAB实现为：例 再看下面的展开式展开式说明是此函数的伪奇点！这里的展开式运算实质上是符号运算，因此在MATLAB中执行此命令前应先定义符号变量，否则MATLAB将给出出错信息！共形映射成像特点复变函数中的共形映射可以解释为从一个平面到另一个平面的单叶保角变换。其最大的特点是能够将一个不规则的区域变换到一个易于我们研究的区域。我们知道，在流体力学中，复数的引入使许多问题的研究变得简单，例如流量、环流，流量、复速度的表示和研究都得到一定的简化。而在研究一个物体表面流体的作用时，如果物体表面是不规则的，例如机翼剖面，那么研究难度将会增大很多。同样在电场，电势研究方面，对于不规则物体的电场研究也是一个难点。而保形变换能够将不规则形状变换到一个规则的形状，例如圆形或者矩形，那么在这个规则的形状下研究，问题将得到很大的简化。由于共形映射具有诸如保角等优良的性质，它已经不单单只为解决数学问题提供帮助，在更为广阔的方面都有所应用。它的各种变换形式更是在不同领域的学科中发挥着重要作用。比较著名的如默比乌斯变换，施瓦茨-克里斯托费尔变换等等。在弹性力学、静电学、热流与流体力学等学科中，共形映射都能采取一种化繁为简的方法解决很多重要而复杂的实际问题。复变函数的理论发展和实际应用日趋成熟，复变函数研究的中心对象解析函数也在理论上形成了较为完整的理论体系，然而解析函数的可视化发展却甚为缓慢。Visual Complex Analysis一书将理论上难以理解的复变函数理论通过二维复平面表现出来，在共形映射的成像方面也介绍甚少。而在各种关于复变函数理论的书籍中，共形映射部分也仅仅是理论上的给出基本共形映射，诸如线性变换，指数变换等变换的区域描述。各种常见共形映射的图形，诸如：线性函数，指数函数，幂函数，儒科夫斯基函数等，通过比较不容单页性区域上的成像，加深这些函数的成像原理和过程。最后给出了一个用Maltab做的复变函数作图工具，能够实现常见区域上的共形映射的二维、三维的动态或者静态图形。线性函数 一般具有形式 的变化成为线性变化。简记为.【注】条件ad-bc0是必要的，否则导致L(z)恒为常数。例3.1 求将上半平面变到的变换，并做出该共形映射的图像，满足： , 。解：从到单位圆的变换为：从单位圆到到圆的函数为复合得到：，所以 所求变化为：实现上述变换的程序以及结果如下：程序3-1：x=-20:0.1:20; %产生上半平面的格点y=0.1:0.1:20;xx,yy=meshgrid(x,y);z=xx+yy*i;w=5*i*(z-i)./(z+i)+1/7*pi %函数subplot(2,2,1);cplxmap(z,w); %自带函数作图subplot(2,2,2);surf(real(z),real(w),imag(w),imag(z); %图二colorbar(vert);title(三维图形);subplot(2,2,3);plot(real(z),imag(z); %原上半平面格点title(原上半平面格点);w=5*i*(z-i)./(z+i)+1/7*pi;subplot(2,2,4);plot(real(w),imag(w); %映射后的图形title(映射后的图形);运行结果：如图3-1图3-1为了详细得到映射过程，即具体从原来的上半平面中的哪部分到新图形的哪部分，我们可以设点原来上半平面的区域范围来作图：程序3-2：x=0:0.1:20; %产生上半平面的格点y=0.1:0.1:20;xx,yy=meshgrid(x,y);z=xx+yy*i;w=5*i*(z-i)./(z+i)+1/7*pi %函数subplot(2,2,1);plot(real(z),imag(z);title(第一象限图);subplot(2,2,2);plot(real(w),imag(w); %图二axis equal;title(作图映射区域);x=-20:0.1:0; %产生上半平面的格点y=0.1:0.1:20;xx,yy=meshgrid(x,y);z=xx+yy*i;w=5*i*(z-i)./(z+i)+1/7*pi %函数subplot(2,2,3);plot(real(z),imag(z); %原上半平面格点title(第二象限);subplot(2,2,4);plot(real(w),imag(w); %映射后的图形axis equal;title(映射后的图形); 运行结果：如图3-2图3-2幂函数：讨论幂函数 ，其中n1且为自然数, 除了z=0及z=外，它处处具有不为0的导数，因而这些点都是保角的。且由于幂函数将角形区域保形变换成角形区域 所以它的单页性区域为：0 arg(z) a ()。例3.2 将 变到平面 .解：用函数程序3-4：a=cplxgrid(40);ll=a.(1/360);z=ll.(36);w=z.7;subplot(3,2,1);cplxmap(z,w);title(w=z7的三维图形);colorbar(vert);subplot(3,2,2);surf(real(z),real(w),imag(w),imag(z);title(w=z7的三维图形);colorbar(vert);a=cplxgrid(40);ll=a.(1/360);z=ll.(36);w=z.3w2=z.10;subplot(3,2,3);cplxmap(z,w);title(w=z3的三维图形);colorbar(vert);subplot(3,2,4);surf(real(z),imag(z),real(w2),imag(w2);title(w=z5的三维图形);colorbar(vert);subplot(3,2,5);plot(real(z),imag(z);title(z平面);axis equal;subplot(3,2,6);plot(real(w),imag(w);title(w平面);%平衡坐标轴axis equal;运行结果：如图3-4图3-43.3 根式函数作为幂函数的逆变换将z平面上的角形区域保形变换成角形区域例3.3 做出和的图形，并进行对比。解：程序3-5：z=cplxgrid(40); w=z.(1/3); w2=z.(1/10);subplot(2,2,1); cplxmap(z,w); % z.(1/3)的三维图像title(w=z.(1/3); subplot(2,2,2);cplxmap(z,w2); %z.(1/10)的三维图形title(w=z.(1/6); subplot(2,2,3);plot(real(z),imag(z); %z的平面图形title(z); subplot(2,2,4);plot(real(w),imag(w); %z.(1/3)的平面图形title(w=z.(1/3);运行结果：如图3-5图3-5指数函数指数函数在任意有限点均有，因此它在z平面是保角的。单页区域是带形区域：，可将它变到：例3.5 做出 的三维和二维图形，单页性区域为： （1） 三维图形：程序3-8：x=-20:0.1:20;y=0:0.1:13*pi/7;xx,yy=meshgrid(x,y);z=xx+yy*i;w=exp(z);subplot(2,1,1);cplxmap(z,w);%surf(real(w),imag(w),real(z),imag(z);colorbar(vert);view(10,3)subplot(2,1,2);cplxmap(w,z);运行结果：如图3-8图3-8(2)二维图形程序3-9：x=-20:0.1:20;y=0:0.1:13*pi/7;xx,yy=meshgrid(x,y);z=xx+yy*i;w=exp(z);subplot(2,1,1);plot(real(z),imag(z);subplot(2,1,2);plot(real(w),imag(w);运行结果：如图3-9图3-93.5 对数函数作为指数函数的逆变换将w平面上的角形区域保形变换成z平面上的带形区域g：，单页性区域为：的扇形区域，保形变换为：。例3.6 做出函数 在如下区域的图形； (1)； (2) ；(3)； (4) .解：程序3-10：a=cplxgrid(40);for i=0:3 z=exp(i*pi/2)*a.(i+1)/4); w=log(z); subplot(2,2,i+1); cplxmap(z,w);end运行结果：如图3-10图3-10程序3-11：a=cplxgrid(40);for k=0:3 z=exp(k+1)*pi/4*i)*a.(k+1)/4); w=log(z); subplot(4,2,2*k+1); plot(real(z),imag(z); axis equal; subplot(4,2,2*k+2); plot(real(w),imag(w);end运行结果：如图3-11图3-11儒科夫斯基变换具有形式的变换，成为儒科夫斯基变换。不是一般性，取a=1，得到。如果令，即则将圆变换为椭圆其中。单页性区域可记为：（1）单位圆盘：；（2）单位圆盘外部：；（3）上半平面：；（4）下半平面：。例3.10对于圆周内部和圆周外部所围成的区域用儒科夫斯基变换。程序3-17：for x=-1.9:0.01:1.1 if x1 y=(-sqrt(1.52-(x+0.4)2)+0.9):0.01:(s