我想在python中过滤信号,灵感来自MATLAB代码。MATLAB具有filter
的函数,该函数应该类似于scipy.signal.lfilter
(来自scipy lfilter()的问题:MATLAB Filter())。但是,我仍然得到ValueError: object of too small depth for desired array
。
MATLAB代码(在八度执行):
% 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过滤器与Python Lfilter不兼容"的问题,我尝试将axis=0
添加到lfilter
,但我仍然有ValueError。
问题
如何在Python中执行MATLAB代码?
Scipy的lfilter
期望b
参数是1-D数组(或"类似阵列",例如列表),而不是标量。例如,
lfilter([1], np.array([1, -0.992217938]), x, axis=0)