《用matlab求留數(shù)》由會(huì)員分享，可在線閱讀，更多相關(guān)《用matlab求留數(shù)(3頁(yè)珍藏版)》請(qǐng)?jiān)诮痄z頭文庫(kù)上搜索。

1、收稿日期: 2006) 05- 29作者簡(jiǎn)介: 賈新民(1956) ), 男, 四川鄰水縣人, 副教授, 研究方向: 計(jì)算機(jī)應(yīng)用1用 MATLAB求 留 數(shù)賈新民( 昌吉學(xué)院計(jì)算機(jī)工程系 新疆 昌吉 831100)摘 要: 文章給出了用 MATLAB求復(fù)變函數(shù)的留數(shù)的方法。關(guān)鍵詞:MATLAB; 復(fù)變函數(shù); 高階極點(diǎn); 留數(shù)中圖分類(lèi)號(hào): O174文獻(xiàn)標(biāo)識(shí)碼: A文章編號(hào): 1671- 6469(2006)03- 0101- 03Matlab 是集數(shù)值計(jì)算、 符號(hào)運(yùn)算及圖形處理 等功能為一體的功能強(qiáng)大的科學(xué)計(jì)算平臺(tái)。Mat -lab 在高等數(shù)學(xué)、 力學(xué)、 機(jī)械、 電工、 信號(hào)處理、 自動(dòng)控制等領(lǐng)。

2、域都有廣泛的應(yīng)用。本文介紹利用 Matalb 求復(fù)函數(shù)留數(shù)的方法。1. 留數(shù)及其求法:留數(shù)又稱(chēng)為殘數(shù)。復(fù)變函數(shù) f( z) 在 p 點(diǎn)的留數(shù)定義為 f( z) 在 p 點(diǎn)羅朗展開(kāi)中( z- p)- 1的系數(shù)c- 1。f( z) 在 p點(diǎn)的留數(shù)記為 res f ( p) 留數(shù)的求法: 設(shè) z0是函數(shù) f( z) 的 m( m 1) 階極點(diǎn), 則Res f( p) =1 ( m- 1) !limzypdm- 1 dzm- 1 ( z- p)mf( z) 。( 1)下面通過(guò)例子說(shuō)明幾種特殊情形下用Matlab求留數(shù)的方法。 2. f( z) 是有理分式的情形Matalab 的信號(hào)處理工具箱中的函數(shù) 。

3、residue的作用是求有理分式展開(kāi)為部分分式時(shí)的展開(kāi)系數(shù), 巧妙地應(yīng)用該函數(shù), 可以求出有理分式的留數(shù) ( 數(shù)值解) 。使用語(yǔ)法為 r, p = residue( 分子多項(xiàng)式, 分母多項(xiàng)式)其中 r 是部分分式的系數(shù)數(shù)組( 留數(shù)數(shù)組) , p 是極點(diǎn)數(shù)組。使用時(shí)需要注意, 當(dāng)函數(shù)有重極點(diǎn)時(shí), 對(duì)同一個(gè)極點(diǎn) p, 存在幾個(gè)展開(kāi)系數(shù) r, 這幾個(gè)r 中只有與相同極點(diǎn)中第一個(gè)對(duì)應(yīng)的 r 是( z-p)- 1的系數(shù)即與極點(diǎn) p 對(duì)應(yīng)的留數(shù), 其余的不是留數(shù)。例1 求函數(shù) f( z) =z+ 3 z2+ 3z+ 2的極點(diǎn)和留數(shù)。 解: matlab 程序及結(jié)果為a= 1, 3 ; b= 1, 3, 2。

4、 ; % a、 b 分別 是 f( z) 分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)矩陣 r, p = residue( a, b) %r 和 p分別是 f( z)的留數(shù)數(shù)組和極點(diǎn)數(shù)組r =- 12p =- 2- 1從結(jié)果看出: res f( - 2) = - 1, res f( - 1) = - 2。例2 求函數(shù) f( z) =z15 ( z2+ 1)2( z4+ 2)3在各極點(diǎn)的留數(shù)。解: 所給函數(shù)一共有 7 個(gè)極點(diǎn): 以及 p= ? i, p=42e(P+ 2KP) i/4, ( k= 0, 1, 2, 3) 以及 p= 。用residue 函數(shù)時(shí), 首先要將 f( z) 分母展開(kāi)為多項(xiàng)式, 用手工展。

5、開(kāi)是比較費(fèi)事的。matlab 的信號(hào)處理工具箱中求卷積的函數(shù) conv 可以用來(lái)進(jìn) 行多項(xiàng)式乘積。用Matlab 求解該問(wèn)題的程序和結(jié)101昌吉學(xué)院學(xué)報(bào) 2006 年第 3 期果如下:a= 1, 0, 1 ;a1= conv( a, a) ; a2= 1, 0, 0, 0, 2 ;a3= conv( a2, a2) ;a4= conv( a3, a2) ;a5= conv( a1, a4) % a5 是分母多項(xiàng)式 系數(shù)矩陣a5 =Columns 1 through 13 1 0 2 0 7 0 12 0 18 0 24 0 20Columns 14 through 170 16 0 8b= 1。

6、, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; %分子多項(xiàng)式系數(shù)矩陣r, p = residue( b, a5) % 調(diào)用 residue 函數(shù)求極點(diǎn)和留數(shù)r = 0. 2037 +0. 1637i- 0. 0207 +0. 0976i- 0. 0139 +0. 0049i 0. 2037 - 0. 1637i- 0. 0207 - 0. 0976i- 0. 0139 - 0. 0049i0. 2037 - 0. 1637i 0. 0207 + 0. 0976i- 0. 0139 - 0. 0049i0. 2037 + 0. 1637i 0. 0。

7、207 - 0. 0976i- 0. 0139 + 0. 0049i0. 0926 + 0. 0000i- 0. 0000 + 0. 0093i 0. 0926 - 0. 0000i- 0. 0000 - 0. 0093ip =0. 8409 + 0. 8409i0. 8409 + 0. 8409i0. 8409 + 0. 8409i0. 8409 - 0. 8409i0. 8409 - 0. 8409i0. 8409 - 0. 8409i- 0. 8409 + 0. 8409i- 0. 8409 + 0. 8409i- 0. 8409 + 0. 8409i- 0. 8409 - 0. 840。

8、9i- 0. 8409 - 0. 8409i- 0. 8409 - 0. 8409i - 0. 0000 + 1. 0000i- 0. 0000 + 1. 0000i- 0. 0000 - 1. 0000i- 0. 0000 - 1. 0000i為了方便, 將最后結(jié)果重新排列為p = r=0. 8409 + 0. 8409i 0. 2037 + 0. 1637i0. 8409 + 0. 8409i - 0. 0207 + 0. 0976i0. 8409 + 0. 8409i - 0. 0139 + 0. 0049i0. 8409 - 0. 8409i 0. 2037 - 0. 1637i0. 。

9、8409 - 0. 8409i - 0. 0207 - 0. 0976i0. 8409 - 0. 8409i - 0. 0139 - 0. 0049i- 0. 8409 + 0. 8409i 0. 2037 - 0. 1637i- 0. 8409 + 0. 8409i 0. 0207 + 0. 0976i- 0. 8409 + 0. 8409i - 0. 0139 - 0. 0049i- 0. 8409 - 0. 8409i 0. 2037 + 0. 1637i- 0. 8409 - 0. 8409i 0. 0207 - 0. 0976i- 0. 8409 - 0. 8409i - 0. 01。

10、39 + 0. 0049i- 0. 0000 + 1. 0000i 0. 0926 + 0. 0000i- 0. 0000 + 1. 0000i - 0. 0000 + 0. 0093i- 0. 0000 - 1. 0000i 0. 0926 - 0. 0000i- 0. 0000 - 1. 0000i - 0. 0000 - 0. 0093i以上結(jié)果中左邊一列是函數(shù)的極點(diǎn)數(shù)組( n重極點(diǎn)算 n 個(gè)極點(diǎn)) , 右邊一列是部分分式的展開(kāi)系數(shù)。由前所述, 只有與重極點(diǎn) p 的第一行對(duì)應(yīng)的 r 值( 即標(biāo)有下劃線的值) 才是留數(shù)。即res f(42eiP/4) = 0. 2037+ 0. 1637i。

11、,res f(42e- iP/4) = 0. 2037- 0. 1637i,res f(42ei3P/4) = 0. 2037- 0. 1637i,res f(42e- i3P/4) = 0. 2037+ 0. 1637i,res f( i) = 0. 0926 res f( - i) = 0. 0926。f( ) 可以根據(jù) f( z) 在所有極點(diǎn)的留數(shù)和等于 0的定理求出為res f( ) = - 26i= 1res f( pi) = - 1。 3. 利用Matlab 的符號(hào)運(yùn)算功能求留數(shù)的數(shù) 102昌吉學(xué)院學(xué)報(bào) 2006 年第 3 期學(xué)表達(dá)式用residue 函數(shù)求留數(shù)雖然方便快捷, 但有。

12、兩個(gè)缺陷: 第一, 要求函數(shù) f( z) 的形式必須是有理分 式; 第二, 只能求數(shù)值解, 不能求解析解。當(dāng)函數(shù) f( z) 的形式不是有理分式時(shí), 求函數(shù)在極點(diǎn) p 處的留數(shù)的方法只能先判斷極點(diǎn)的重?cái)?shù), 然后利用Matlab 的符號(hào)運(yùn)算功能根據(jù)公式( 1) 進(jìn)行計(jì)算, 可以求出留數(shù)的解析形式。例3 已知 f( z) =1- e2x z4, 求 res f( 0) 。解: z= 0 是 f( z) 的 3 階極點(diǎn), 將 m= 3 代入求留數(shù)的公式( 1) , 輸入以下Matlab 語(yǔ)句:limit( diff( ( . z3* ( 1- exp( 2* z) ) / z 4. ) , 2) ,。

13、. z. , 0) / 2 % 用單撇號(hào)括起來(lái)表明是符號(hào)運(yùn)算運(yùn)算結(jié)果為ans=- 4/ 3即 res f( 0) = - 4/ 3。例4 設(shè) f( z) =ez- 1 sin3z, 求 res f( kP ) , k= 1, 2,3,解: 當(dāng) kX 0 時(shí), z= kP是 f( z) 的 3 階極點(diǎn), 在 Matlab 的命令窗口輸入語(yǔ)句:limit( diff( ( -( z- pi) 3* ( exp( z) - 1) / ( sin( z) )3. ) , 2) , . z. , pi) / 2運(yùn)算結(jié)果為 ans=- exp( pi) + 1/ 2即 res f( P) = - eP+1。

14、 2;當(dāng)k= 2 時(shí), 將上述語(yǔ)句中 pi 改為 2* pi, 得limit( diff( ( -( z- 2* pi) 3* ( exp( z) - 1) / ( sin( z) ) 3. ) , 2) , . z. , 2* pi) / 2ans= exp( 2* pi) - 1/ 2即 res f( 2P) = e2P-1 2。同樣可以求出 k= 3P時(shí)的留數(shù)為res f( 3P) = - e3P+1 2。歸納以上結(jié)果有res f( kP) = ( - 1)kekP-( - 1)k2。f( z) 的形式是有理分式時(shí), 也可以用符號(hào)運(yùn)算方法求留數(shù)。例5 設(shè) f( z) =z15 ( z2+。

15、 1)2( z4+ 2)3, 求 res f( i) 。解: z= i 是 f( z) d 的二階極點(diǎn), 在 Matlab 的命 令窗口, 輸入語(yǔ)句:limit( diff( ( . ( z- i) 2* z15/ ( ( z2+ 1) 2* ( z4+ 2) 3) . ) , 1) , . z. , i) 回車(chē)后可得結(jié)果為ans=5/ 54即res f( i) =5 54( = 0. 092592592) , 例 2 數(shù)值 計(jì)算是保留了 4 位小數(shù)的結(jié)果。4. 簡(jiǎn)單的評(píng)述上面介紹的兩種求復(fù)變函數(shù)的留數(shù)的方法各有長(zhǎng)處。看起來(lái)似乎用符號(hào)運(yùn)算方法既可以求一 般函數(shù)的留數(shù), 也可以求有理分式的留數(shù), 而且得到的結(jié)果是解析的, 應(yīng)該更有效。但事實(shí)上不一定。仍以例 2 的函數(shù)為例, p= 21/ 4eiP/4是函數(shù)的 3 階極點(diǎn), 根據(jù)例 3 和例 4, 用語(yǔ)句 limit( diff( (. ( z-2( 1/ 4) * exp( i* pi/ 4) ) 3* z15/ ( ( z2+ 1) 2* ( z4+ 2) 3) . ) , 2) , . z. , 2( 1/ 4) * exp( i* pi) ) / 2 就可以求出函數(shù)在該點(diǎn)的留數(shù)。輸入上述語(yǔ)句后, 經(jīng) 過(guò)一段時(shí)間的等待后, 結(jié)果卻是好幾百位的一個(gè)數(shù), 可見(jiàn), 至少在這種情況下, 數(shù)值解法雖然不。

總結(jié)

以上是生活随笔為你收集整理的matlab 求留数,用matlab求留数的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。