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) – 直接多项式项的系数。

警告

此函数可能会同步设备。