实验一计算复变函数极限微分积分留数泰勒级数展开式.pptxVIP

实验一计算复变函数极限微分积分留数泰勒级数展开式会计学第1页/共17页MATLAB实现内容1、MATLAB求复变函数极限2、MATLAB求复变函数微分3、MATLAB求复变函数积分4、MATLAB求复变函数在孤立奇点的留数5、MATLAB求复变函数的泰勒级数展开式第2页/共17页MATLAB基本命令1、MATLAB求复变函数极限用函数limit求复变函数极限【Matlab源程序】 syms z f=;limit(f,z,z0)返回极限结果例1 求 在z=0的极限 解 syms z;f=z*exp(z)/(sin(z))limit(f,z,0)ans = 1例2 设求第3页/共17页解 【Matlab源程序】syms z f=sin(z)/z;limit(f,z,0)ans=1limit(f,z,1+i). .ans=1/2*sin(1)*cosh(1)-1/2*i*sin(1)*cosh(1)+1/2*i*cos(1)*sinh(1)+1/2*cos(1)*sinh(1)例3设第4页/共17页2、 MATLAB求复变函数微分用函数diff求复变函数极限【Matlab源程序】 syms z f=（）; diff(f,z)返回微分结果解 syms z f=exp(z)/((1+z)*(sin(z)));diff(f) ans = exp(z)/(1+z)/sin(z)-exp(z)/(1+z)^2/sin(z)-exp(z)/(1+z)/sin(z)^2*cos(z)例 4 求积分第5页/共17页3、 MATLAB求复变函数积分(1)用函数int求解非闭合路径的积分.【Matlab源程序】syms z a bf=int(f,z,a,b)返回积分结果解syms z x1=int(cosh(3*z),z,pi/6*i,0)x2=int((z-1)*exp(-z),z,0,i)结果为：x1 = -1/3*ix2 = -i/exp(i)例5 计算积分的值．第6页/共17页(2) 用函数int 求解闭合路径的积分.解 【Matlab源程序】 syms t z z=2*cos(t)+i*2*sin(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若只输出6位有效数值，使用语句 vpa(inc,6)结果为 ans =.238258e-5+.313254e-4*i第7页/共17页4、 MATLAB求复变函数在孤立奇点的留数（1）f(z)=p(z)/q(z)；p(z)、q(z)都是按降幂排列的 多项式 用函数residue求f(z)=p(z)/q(z)在孤立奇点的留数 【Matlab源程序】 [R,P,K]= residue (B,A) 返回留数，极点 说明：向量B为f(z)的分子系数； 向量A为f(z)的分母系数； 向量R为留数； 向量P为极点位置； 向量k为直接项: 例6 求函数在奇点处的留数．例7 计算积分的值，其中C是正向圆周第8页/共17页解 [R,P,K]= residue([1,0,1],[1,1])结果为：R= 2P = -1K = 1 -1．解 先求被积函数的留数 [R,P,K]= residue ([1,0],[1,0,0,0,-1]) 结果为： 可见在圆周第9页/共17页R = 0.2500 0.2500 -0.2500 + 0.0000i -0.2500 - 0.0000iP = -1.0000 1.0000 0.0000 + 1.0000i 0.0000 - 1.0000iK = []内有四个极点,所以积分值等于S=2*pi*i*sum(R)结果为S =0故原积分第10页/共17页（2）如果已知函数奇点z0的重数为m,则可用下面的MATLAB语句求出相应的留数．R=limit(F*(z-z0),z,z0) %单奇点R=limit(diff(F*(z-z0)^m,z,m-1)/prod(1:m-1);z,z0) % m重奇点例8 求函数在孤立奇点处的留数．解分析原函数可知：是三重奇点，是单奇点，因此可以直接使用下面的第11页/共17页MATLAB语句分别求出这两个奇点的留数．syms zf=sin(z+pi/3)*exp(-2*z)/(z^3*(z-1))R=limit(diff(f*z^3,z,2)/prod(1:2),z,0)结果为：R = -1/4*3^(1/2)+1/2?;limit(f*(z-1),z,1)ans = 1/2*exp(-2)*sin(1)+1/2*exp(-2)*cos