Matlab在复变函数中应用数学实验（一）华中科技大学数学系二○○一年十月MATLAB在复变函数中的应用复变函数的运算是实变函数运算的一种延伸，但由于其自身的一些特殊的性质而显得不同，特别是当它引进了“留数”的概念，且在引入了Taylor级数展开Laplace 变换和Fourier变换之后而使其显得更为重要了。

使用MATLAB来进行复变函数的各种运算；介绍留数的概念及MAT–LAB的实现；介绍在复变函数中有重要应用的Taylor展开（Laurent展开Laplace变换和Fourier变换）。

1 复数和复矩阵的生成在MATLAB中，复数单位为)1ji，其值在工作空间中都显示为=sq rt=(-0+。

.1i00001.1 复数的生成复数可由iz+=。

a=语句生成，也可简写成biaz*+b另一种生成复数的语句是)exp(ithetar=，也可简写成)=，z*exp(theta*irz*其中theta为复数辐角的弧度值，r为复数的模。

1.2 创建复矩阵创建复矩阵的方法有两种。

（1）如同一般的矩阵一样以前面介绍的几种方式输入矩阵例如：)]iA**ii=+3[i*-+*,),235336exp(23,exp(9（2）可将实、虚矩阵分开创建，再写成和的形式例如：)2,3(re=；randim=；)2,3(randim i re com *+=]5466.07271.05681.02897.07027.05341.08385.03420.03704.03412.03093.06602.0[i ii i ii com ++++++=注意 实、虚矩阵应大小相同。

2 复数的运算1．复数的实部和虚部复数的实部和虚部的提取可由函数real 和imag 实现。

调用形式 )(x real返回复数x 的实部)(x imag返回复数x 的虚部2．共轭复数复数的共轭可由函数conj 实现。

调用形式)(x conj返回复数x 的共轭复数3．复数的模和辐角复数的模和辐角的求解由功能函数abs 和angle 实现。

调用形式 )(x abs 复数x 的模)(x angle复数x 的辐角例：求下列复数的实部与虚部、共轭复数、模与辐角 （1）i231+ （2）i i i --131 （3）ii i 2)52)(43(-+（4）i i i +-2184由MATLAB 输入如下：]21^48^,2/)52()43(),1/(3/1),23/(1[i i i i i i i i i i a +*--*=--+==ai iii0000.30000.10000.135000.35000.25000.11538.02308.0-----)(a real %实部=ans0.23081.5000 –3.50001.0000)(a imag%虚部=ans–0.1538–2.5000 –13.0000 –3.0000)(a conj%共轭复数=ans0.2308+0.1538i 1.5000+2.5000i–3.5000+13.0000i1.0000+3.0000i)(a abs%模=ans0.27742.9155 13.4629 3.1623)(a angle%辐角=ans–0.5880–1.0304–1.8228-1.24904．复数的乘除法复数的乘除法运算由“/”和“*”实现。

例 复数的乘除法演示。

)3/exp(4i pi x *==xi 4641.30000.2-)5/exp(3i pi y *= =yi 7634.14271.2-)5/exp(31i pi y **==1yi 7634.14271.2+y x /=ansi 5423.02181.1-1/y x=ans I 3260.11394.0-由此例可见，i 5/)( 相当于)5/()(i * ，和i *5/)( 不相等。

5．复数的平方根复灵敏的平方根运算由函数sprt 实现。

调用形式)(x sprt返回复数x 的平方根值6．复数的幂运算复数的幂运算的形式为n x ^，结果返回复数x 的n 次幂。

例 求下列各式的值 )6/1()^1(-=ans0.8660+0.5000 i7．复数的指数和对数运算复数的指数和对数运算分别由函数exp和log实现。

调用形式exp(x返回复数x的以e为底的指数值))log(x返回复数x的以e为底的对数值例求下列式的值（参见参考资料【4】P.68.2–15）。

log(i-)ans=0-.15708i+-log(i)34ans=6094.1+.2i21438．复数的三角函数运算复数的三角函数运算函数参见下面的复数三角函数复数三角函数表9． 复数方程求根复数方程求根或实方程的复数根求解也由函数solve 实现。

见下面的例子. 例 求方程083=+x 所有的根（参见参考资料【4】P.32.1–16）。

)083^('=+'x solve=ans [ –2])]2/1(^31[*-i)]2/1(^31[*+i3 留数留数定义：设a 是)(z f 的孤立奇点，C 是a 的充分小看邻域内一条把a 点包含在其内部的闭路，积分⎰C dz z f i)(21π称为)(z f 在a 点的留数或残数，记作]),([Re a z f s 。

在MATLAB 中，可由函数residue 实现。

residue 留数函数（部分分式展开）),(],,[A B residue K P R =函数返回留数，极点和2个多项式比值)(/)(s A s B 的部分分式展开的直接项。

)()()()2()2()1()1()()(s K n P s n R P s R P s R s A s B +-++-+-= 如果没有重根，则向量B 和A 为分子、分母以s 降幂排列的多项式系数，留数返回为向量R 、极点在向量P 的位置，直接项返回到向量K 。

极点的数目)()(1)(P len g th R len g th A len g th n ==-=。

如果)()(A length B length <，则直接项系数为空；否则1)()()(+-=A length B length K length 。

如果存在M 重极点即有)1()(-+==m j P j P 则展开项包括以下形式mj P s m j R j P s j R j P s j R ))(()1())(()1()()(2--+++-++- ),,(],[K P R residue A B = 有3个输入变量和2个输出变量，函数转换部分因式展开还为系数为B 和A 的多项式比的形式。

注意：数值上讲，分式多项式的部分因式展开实际上代表了一类病态问题。

如果分母多项式)(S A 是一个近似有重根的多项式，则在数值上的一点微小变化，包括舍入误差都可能造成极点和留数结果上的巨大变化。

因此使用状态空间和零点—极点表述的方法是可取的。

例 求如下函数的奇点处的留数。

zz z 212-+ 在MATLAB 实现如下])0,2,1[],1,1([],,[-=residue k p r=r 1.5000–0.5000=p2=k[ ]所以可得5.0]0),([Re;5.1]2),([Re-==zfszfs。

例计算下面的积分⎰-C dzzz1 4其中C为正向圆周2||=z。

（参见参考资料【4】P.158.例2）解：先求被积函数的留数])1,0,0,0,1[],0,1([],,[-=residuekpr=r0.25000.2500–0.2500–0.0000 i–0.250+0.0000 i=p–1.00001.00000.0000+1.0000 i0.0000–1.0000 i=k[ ]可见在圆周2||=z内有四个极点，所以积分值等于)25.025.025.025.0(2=--+**pi。

4 Taylor级数展开Taylor级数开展在复变函数中有很重要的地位，如分析复变函数的解析性等。

函数)(xf在0xx=点的Taylor级数开展为+-''+-'+-+=!3/)0)(0(!2/)0)(0()0)(0(0)(3^2^x x x f x x x f x x x f x x f在MATLAB 中可由函数taylor 来实现。

taylor 泰勒级数展开 )(f taylor 返回f 函数的五次幂多项式近似。

此功能函数可有3个附加参数。

),(n f taylor 返回1-n 次幂多项式。

),(a f taylor 返回a 点附近的幂多项式近似。

),(x r taylor使用独立变量代替函数)(f findsym 。

例 求下列函数在指定点的泰勒开展式（参见参考资料【4】P.143.12）。

（1）10,/12-=z z （2）4/0,pi z tgz =；MATLAB 实现为： )1,2^/1(-x taylor5)^1(64)^1(53)^1(42)^1(323+*++*++*++*+*+x x x x x)4/),(tan(pi x taylor=ans *+*-*+*-*+*-*+3/103)^4/1(3/82)^4/1(22/121pi x pi x pi x5)^4/1(15/644)^4/1(pi x pi x *-*+*-例 再看下面的展开式 )10,/)(sin(x x taylor=ans8^362880/16^5040/14^120/12^6/11x x x x *+*-*+*-展开式说明0=x 是此函数的伪奇点！这里的taylor 展开式运算实质上是符号运算，因此在MATLAB 中执行此命令前应先定义符号变量syms,，否则MATLAB将给出出错信息！xz5 Laplace变换及其逆变换1．Laplace变换L=返回以默认独立变量T对符号函数F的Laplace变换。

函数)laplace(F返回默认为s的函数。

如果)F=，则Laplace函数返回tF(s的函数)LL=。

其中定义L为对t的积分(tsFtL*=。

s*-(t,0i n f)exp(),)(int()(tFl a p l a c e等价于, L=以t代替s的Laplace变换。

),(t)FlaplaceFttx-=。

xL**,0))(infexp(),int(()laplaceL=以z代替s的Laplace变换（相对于w的积分）。

wF)(z,,wz)int((wL*(=。

*-zF(z,))w,exp(),,0Flaplace等价于inf)例如：syms a s t w x)5^laplace(xans=120s/6^paplace*a(exp(s))ans=t-1a/()wlaplace*x(sin(t),)ans=)2^2^/(w t w +),),((t w w x cons laplace *=ans)2^2^/(x t t + )),2/3(^(t sym x laplace=ans)2/5(^/)2/1(^4/3t pi *))))((((''x F sym diff laplace=ans)0(),),((F s s x x F laplace -*2．Laplace 逆变换)(L ilaplace F = 返回以默认独立变量s 的数量符号L 的Laplace 变换，默认返回t 的函数。