【复变函数—留数：有理函数的部分分式展开MATLAB实现】

概念说明

1.有理函数是一类能求出原函数的重要函数类，其依据就是它有部分分式展开。

2.对有理函数的展开方法一般使用待定系数法，即通过解一组线形方程组来确定系数。

3.在MATLAB的行向量中，num和 den分别表示有理分式的分子和分母的系数。

4.MATLAB中有理函数部分分式展开语句命令为:

[r, p, k]=residue (num, den)

其中r是部分分式的系数数组，也是留数数组，p是极点数组，k是余项。

实例说明

给出一个函数，如下，求其部分分式展开：

那么，我们可以直接利用该命令:

[r, p, k]=residue (num, den)

就可以求出该函数的部分分式展开。

但是注意：

这个函数的分母不是多项式形式，这里是3次幂可以算出；
如果次数高了，怎么求其分母系数？
那么，我们可以用如下MATLAB程序：

clear all;clc
%%有理函数的部分分式展开，分母展开为多项式的系数
a=input('输入:a=');%比如这个函数a=[1 1]
n=input('输入：n=');%比如这个函数n=1（因为做了conv运算，这里n=幂次-2。）

b=conv(a,a);
for i=1:n
    b=conv(b,a);
end
b

那么得到的b，即为分母展开后的行向量，即是den向量。

所以，再调用：

num=[1 4 6];
den=b;
[r, p, k]=residue (num, den)

运行之后，即可得到r，p，k3个数组，即r是部分分式的系数数组，也是留数数组，p是极点数组，k是余项。

如图：

K这里是空数组，说明该函数展开后，余项为0。

回顾知识点

有理函数的部分分式展开：

代码语句	解释说明
[r, p, k]=residue (num, den)	r是部分分式的系数数组，也是留数数组，p是极点数组，k是余项
conv(b,a)	卷积/多项式相乘