matlab与复变函数.docVIP

matlab与复变函数 复数和复矩阵的生成 在matlab中，虚数单位为i=j=sqrt(-1),其值在工作空间中都显示为0 + 1.0000i 1．复数的生成 复数可由z=a+b*i语句生成，也可简写成z=a+bi 其中a为复数的实部，b为复数的虚部。 复数可由z=r*exp(i*theta)生成，也可简写成z=r*exp(theta i) 其中theta为复数辐角的弧度值，r为复数的模。 2．创建复矩阵 同一般的矩阵一样直接输入其元素 【例1】 A=[3+6*i -2-8*i;5*exp(i*7) 9-5i] 可用实部、虚部分别创建矩阵，再写成和的形式 【例2】 re=rand(3,2); im=randn(3,2); com=re+im*i com = 0.9501 - 0.4326i 0.4860 + 0.2877i 0.2311 - 1.6656i 0.8913 - 1.1465i 0.6068 + 0.1253i 0.7621 + 1.1909i 函数complex( )可以生成复数、复矩阵 【例3】 complex(3,8) ans = 3.0000 + 8.0000i 复数的运算 real(z) 返回复数z的实部 imag(z) 返回复数z的虚部 conj(z) 返回复数z的共轭复数 abs(z) 返回复数z的模（绝对值） angle(z) 返回复数z的辐角 sqrt(z) 返回复数z的平方根值 exp(z) 返回复数z的以e为底的指数值 log(z) 返回复数z的以e为底的对数值 复数的三角函数和反三角函数与实函数符号相同 以上参数z也可以为复矩阵。 【例4】 z=3+6*i; re=real(z); im=imag(z); ab=abs(z); an=angle(z); sq=sqrt(z); ex=exp(z); lo=log(z); si=sin(z); ac=acos(z); 【例5】求方程z^3=8=0所有的根。 solve(z^3+8=0) ans = [ -2] [ 1-i*3^(1/2)] [ 1+i*3^(1/2)] 复变函数的极限、导数、积分与级数 求复变函数的极限、导数、积分与级数等的命令与实函数相同，不再重复，仅给出范例。 【例6】求极限。 syms z f=z/(1+z); limit(f,z,1+5*i) ans = 27/29+5/29*i 【例7】求在的导数。 syms z f=log(1+sin(z)); df=diff(f,z) df = cos(z)/(1+sin(z)) vdf=subs(df,z,i/2) vdf = 0.8868 - 0.4621i 【例8】计算。 syms z f=z*cos(z); intf=int(f,z,0,i) intf = cosh(1)-sinh(1)-1 vpa(intf,8) ans = -.6321206 【例9】计算积分为正向圆周：。 syms t z; z=2*exp(i*t); f=1/((z+i)^10*(z-1)*(z-3)); inc=int(f*diff(z),t,0,2*pi) inc = 779i*pi+237/312500000*pi 【例10】求幂级数的收敛半径。 syms n; fn=exp(i*pi/n); R=abs(limit(fn^(1/n),n,inf) R=abs(limit(fn^(1/n),n,inf)) R = 1 【例11】求函数在点的泰勒展开式。 syms z taylor(tan(z),i,2) ans = i*tanh(1)+(1-tanh(1)^2)*(z-i) 有理分式函数留数的计算 [R,P,K]=residue(B,A) 返回留数、极点和两个多项式比值B(s)/A(s)的部分分式展开的直接项 向量B和A为分子、分母以s降幂排列的多项式系数，向量R是返回的留数，向量P是返回的极点，向量K由B(s)/A(s)的商的多项式系数组成，如果length(B)length(A)，则K为空向量，否则，length(K)=length(B)-length(A)+1。 如果存在m级极点即有，则展开式包括以下形式： [B,A]=residue(R,P,K) 根据已知的奇点P、奇点的留数R和K来计算有理分式函数的系数B和A 当复变函数的分子和分母部分不是多项式时，则residue( )命令不 能解决问题，此时主要根据留数计算规则进行处理。 【例12】计算函数的留数。 B=[1 3 0 2];A=[1 6 -1]