cupyx.scipy.signal.residuez#
- cupyx.scipy.signal.residuez(b, a, tol=0.001, rtype='avg')[源代码]#
计算 b(z) / a(z) 的部分分式展开。
如果分子 b 的次数为 M,分母 a 的次数为 N
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 的注意事项。
- 参数:
b (array_like) – 分子多项式系数。
a (array_like) – 分母多项式系数。
tol (float, optional) – 用于确定两个根在它们之间距离上是否相等的容差。默认为 1e-3。更多详情请参见 unique_roots。
rtype ({'avg', 'min', 'max'}, optional) – 用于计算代表一组相同根的根的方法。默认为 'avg'。更多详情请参见 unique_roots。
- 返回:
r (ndarray) – 与极点对应的留数。对于重复极点,留数按幂分数升序排列。
p (ndarray) – 按幅值升序排列的极点。
k (ndarray) – 直接多项式项的系数。
警告
此函数可能会同步设备。