scipy.signal.
residuez#
- scipy.signal.residuez(b, a, tol=0.001, rtype='avg')[源代码]#
计算 b(z) / a(z) 的部分分式展开。
如果 M 是分子 b 的阶数,N 是分母 a 的阶数
b(z) b[0] + b[1] z**(-1) + ... + b[M] z**(-M) H(z) = ------ = ------------------------------------------ a(z) a[0] + a[1] z**(-1) + ... + a[N] z**(-N)
则部分分式展开 H(z) 定义为
r[0] r[-1] = --------------- + ... + ---------------- + k[0] + k[1]z**(-1) ... (1-p[0]z**(-1)) (1-p[-1]z**(-1))
如果存在任何重复的根(距离小于 tol),则部分分式展开具有如下形式的项
r[i] r[i+1] r[i+n-1] -------------- + ------------------ + ... + ------------------ (1-p[i]z**(-1)) (1-p[i]z**(-1))**2 (1-p[i]z**(-1))**n
此函数用于 z 的负幂的多项式,例如 DSP 中的数字滤波器。对于正幂,请使用
residue
。有关该算法的详细信息,请参阅
residue
的注释。- 参数:
- barray_like
分子多项式系数。
- aarray_like
分母多项式系数。
- tolfloat,可选
两个根被认为相等时,它们之间的距离容差。默认值为 1e-3。有关详细信息,请参阅
unique_roots
。- rtype{‘avg’, ‘min’, ‘max’},可选
用于计算表示一组相同根的根的方法。默认值为 ‘avg’。有关详细信息,请参阅
unique_roots
。
- 返回:
- rndarray
对应于极点的残差。对于重复的极点,残差的顺序与按幂分数升序排列相对应。
- pndarray
按大小升序排列的极点。
- kndarray
直接多项式项的系数。
另请参阅