王福昌
复变函数问题的MATLAB求解之二——留数的计算
2018-5-28 17:37
阅读:6612

 1、复函数导数

    复函数导数与实数函数的求导法则基本相同

例1、 已知复变函数f(z)=\frac{z^2+3z+4}{(z-1)^5}, 求 f^{(3)}(-i\sqrt{5}).

>> f=(z^2+3*z+4)/(z-1)^5;f3 =diff(f,z,3);d3 = subs(f3,z,-sqrt(-5))

 d3 = 

- 30/(1 + 5^(1/2)*1i)^6 + (90*(- 3 + 5^(1/2)*2i))/(1 + 5^(1/2)*1i)^7 + (210*(1 + 5^(1/2)*3i))/(1 + 5^(1/2)*1i)^8

>> simplify(d3)

ans =

3565/4374 - (5^(1/2)*650i)/2187

2、留数的计算

法则I:a为f(z)的一阶极点, Res[f(z),a]=lim_{z\rightarrow a}(z-a)f(z)

法则II:a为P(z)/Q(z)的一阶极点, Res[P(z)/Q(z),a]=P(a)/Q'(a)

法则III:a为f(z)的m阶极点,Res[f(z),a]=1/(m-1)!lim_{z\rightarrow a}d^{m-1}/dz^{m-1} [(z-a)^m f(z)]

例2 求复变函数 f(z)=\frac{1}{z^3(z-1)}sin z exp(-2z) 的留数

解:0为三阶极点,1为简单极点,分别用法则III和法则I可得

>> syms z; f = sin(z+pi/3)*exp(-2*z)/(z^3*(z-1));r1 = limit(diff(z^3*f,z,2)/prod(1:2),z,0)

 r1 =

 1/2 - 3^(1/2)/4


 >> r1 = limit((z-1)*f,z,1)

 r1 =

 exp(-2)*sin(pi/3 + 1)


转载本文请联系原作者获取授权,同时请注明本文来自王福昌科学网博客。

链接地址:https://wap.sciencenet.cn/blog-292361-1116161.html?mobile=1

收藏

分享到:

当前推荐数:0
推荐到博客首页
网友评论0 条评论
确定删除指定的回复吗?
确定删除本博文吗?