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')])
运算结果(e-16为浮点计算误差,可以忽略):
标准结果: