【问题标题】：Matlab filter with Python lfilter ValueError object of too small depthMatlab过滤器与深度太小的Python lfilter ValueError对象
【发布时间】：2019-07-18 11:55:54
【问题描述】：

我想在 Python 中过滤一个信号，灵感来自于 Matlab 代码。 Matlab有函数filter，应该和scipy.signal.lfilter类似（来自问题：Matlab filter() with SciPy lfilter()）。但是，我仍然收到ValueError: object of too small depth for desired array。

Matlab 代码（以 Octave 执行）：

% Matlab
x = [1.0485e-04  -2.4193e-04  -3.0078e-04  1.5750e-03  -1.9698e-03  1.3902e-04  2.7568e-03  -3.8059e-03  2.0123e-03  3.3257e-03]
xfilt = filter(1, [1 -0.992217938], x);
disp(xfilt);

% output
1.0485e-04  -1.3790e-04  -4.3760e-04   1.1408e-03  -8.3788e-04  -6.9233e-04   2.0699e-03  -1.7522e-03   2.7378e-04   3.5974e-03

Python：

# Python
from scipy.signal import lfilter

x = np.array([1.0485e-04, -2.4193e-04, -3.0078e-04, 1.5750e-03, -1.9698e-03, 1.3902e-04, 2.7568e-03, -3.8059e-03, 2.0123e-03, 3.3257e-03])
lfilter(1, np.array([1, -0.992217938]), x, axis=0)

导致错误的原因：

ValueError                                Traceback (most recent call last)
<ipython-input-87-d5c23d362b45> in <module>
      1 x = np.array([1.0485e-04, -2.4193e-04, -3.0078e-04, 1.5750e-03, -1.9698e-03, 1.3902e-04, 2.7568e-03, -3.8059e-03, 2.0123e-03, 3.3257e-03])
----> 2 print(lfilter(1, np.array([1, -0.992217938]), x, axis=0))

~/anaconda3/envs/*env*/lib/python3.6/site-packages/scipy/signal/signaltools.py in lfilter(b, a, x, axis, zi)
   1378     else:
   1379         if zi is None:
-> 1380             return sigtools._linear_filter(b, a, x, axis)
   1381         else:
   1382             return sigtools._linear_filter(b, a, x, axis, zi)

ValueError: object of too small depth for desired array

系统

Python：3.6.8
Scipy：1.2.0

试过

基于问题“Matlab filter not compatible with Python lfilter”，我尝试将axis=0 添加到lfilter，但仍然收到ValueError。

问题

如何在 Python 中做 Matlab 代码？

【问题讨论】：

标签： python matlab scipy filtering

【解决方案1】：

SciPy 的 lfilter 期望 b 参数是一维数组（或“类数组”，例如列表），而不是标量。例如，

lfilter([1], np.array([1, -0.992217938]), x, axis=0)

【讨论】：

谢谢。 Python 中的输出与 Matlab 输出相同。