计算数字滤波器的频率响应。
给定数字滤波器的M-order分子b和N-order分母a,计算其频率响应:
jw -jw -jwM
jw B(e ) b[0] + b[1]e + ... + b[M]e
H(e ) = ------ = -----------------------------------
jw -jw -jwN
A(e ) a[0] + a[1]e + ... + a[N]e
参数:
b:array_like线性滤波器的分子。如果b的维数大于1,则假定系数存储在第一维中,并且b.shape[1:],a.shape[1:],并且频率阵列的形状必须与广播兼容。
a:array_like线性滤波器的分母。如果b的维数大于1,则假定系数存储在第一维中,并且b.shape[1:],a.shape[1:],并且频率阵列的形状必须与广播兼容。
worN:{None, int, array_like}, 可选参数如果是单个整数,则以那么多的频率进行计算(默认值为N = 512)。这是一种方便的替代方法:
np.linspace(0, fs if whole else fs/2, N, endpoint=False)
使用快速进行FFT计算的数字可以加快计算速度(请参见注释)。
如果为数组,则以给定的频率计算响应。这些单位与fs相同。
whole:bool, 可选参数通常,频率是从0到奈奎斯特频率fs /2(unit-circle的upper-half)计算得出的。如果整为真,则计算从0到fs的频