Matlab的多项式留数与极点的计算

Matlab的多项式留数与极点的计算

以下面的多项式为例：

运算代码：

clc
clear
close

syms p
% 定义多项式
Zp=(5*p^5+71*p^3+70*p)/(2*p^6+35*p^4+117*p^2+36);
% 提取分子与分母
[I,D]=numden(Zp);
I=double(coeffs(I,p,"All"));%分子
D=double(coeffs(D,p,"All"));%分母
% 计算极点与留数
[r,p] = residue(I,D);

disp(['该多项式的留数为：',num2str(r')])
disp(['该多项式的极点为：',num2str(p')])
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

运算结果（e-16为浮点计算误差，可以忽略）：

标准结果：