防灾数学分享 http://blog.sciencenet.cn/u/fzmath 防灾科技学院数学教研室

博文

复变函数问题的MATLAB求解之二——留数的计算

已有 6619 次阅读 2018-5-28 17:37 |系统分类:教学心得

 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://blog.sciencenet.cn/blog-292361-1116161.html

上一篇:达芬方程的分叉图
下一篇:2018高考数学真题解析与点评Word全打包
收藏 IP: 124.238.130.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-30 13:03

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部