Question

我正在尝试为Matlab频率响应函数实现Python等价物

[h,f] = freqz(b, 1, 512, 12.5)

在here中描述。我目前的尝试

f, h = scipy.signal.freqz(b, 1)

没有给出预期的结果。尝试参数worN和whole（请参阅here）似乎无法解决问题。

应该怎么做？

修改

Matlab示例：

>> [h, f] = freqz(1:5, 1, 512, 12.5)

h =

  15.0000 + 0.0000i
  14.9976 - 0.2454i
  14.9902 - 0.4907i
  14.9780 - 0.7358i
  14.9609 - 0.9806i
  14.9389 - 1.2250i
  ...

f =

         0
    0.0122
    0.0244
    0.0366
    0.0488
    0.0610
    ...

Python示例：

>>> f, h = scipy.signal.freqz(range(1,6), 1)
>>> h
array([ 15.00000000 +0.j        ,  14.99755288 -0.24541945j,
        14.99021268 -0.49073403j,  14.97798292 -0.73583892j,
        14.96086947 -0.98062944j,  14.93888050 -1.22500102j,
        ...])
>>> f
array([ 0.        ,  0.00613592,  0.01227185,  0.01840777,  0.02454369,
        0.03067962,  ...])

换句话说，Scipy函数为h提供了良好的值，但f的值不匹配。

Answer 1

在两种语言freqz中，第一个参数需要分子系数b，而不是像你写的那样a。应该是

freqz(b, a, ...)

看起来您正在尝试查找FIR滤波器的响应，其中只有分子系数且a始终为1。

Answer 2

signal.signal.freqz没有以Hz为单位返回频率的选项，因此您之后必须自己缩放频率。

相当于Matlab的

[h, f] = freqz(b, a, n, fs)

使用freqz中的scipy.signal：

w, h = freqz(b, a, worN=n)
f = fs * w / (2*np.pi)

例如，

In [15]: import numpy as np

In [16]: from scipy.signal import freqz

In [17]: w, h = freqz(range(1,6), 1, worN=512)

In [18]: h[:6]
Out[18]: 
array([ 15.00000000+0.j        ,  14.99755288-0.24541945j,
        14.99021268-0.49073403j,  14.97798292-0.73583892j,
        14.96086947-0.98062944j,  14.93888050-1.22500102j])

In [19]: w[:6]
Out[19]: 
array([ 0.        ,  0.00613592,  0.01227185,  0.01840777,  0.02454369,
        0.03067962])

In [20]: f = 12.5*w/(2*np.pi)

In [21]: f[:6]
Out[21]: 
array([ 0.        ,  0.01220703,  0.02441406,  0.03662109,  0.04882812,
        0.06103516])

Python中的matlab freqz函数

2 个答案: