尝试使用 scipy.signal.cwt 时出现奇怪的错误:
我有一些列表 c，我想像这样进行连续小波变换:
scipy.signal.cwt(np.array(c), scipy.signal.morlet, np.arange(.01,.1,.01))
我得到一个奇怪的错误:
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-55-5af5e14b96cd> in <module>()
----> 1 sig.cwt(a, sig.morlet, np.arange(.01,.1,.01))
/usr/local/lib/python2.7/site-packages/scipy/signal/wavelets.pyc in cwt(data, wavelet, widths)
359 wavelet_data = wavelet(min(10 * width, len(data)), width)
360 output[ind, :] = convolve(data, wavelet_data,
--> 361 mode='same')
362 return output
/usr/local/lib/python2.7/site-packages/scipy/signal/signaltools.pyc in convolve(in1, in2, mode)
270
271 if np.iscomplexobj(kernel):
--> 272 return correlate(volume, kernel[slice_obj].conj(), mode)
273 else:
274 return correlate(volume, kernel[slice_obj], mode)
/usr/local/lib/python2.7/site-packages/scipy/signal/signaltools.pyc in correlate(in1, in2, mode)
129 in1zpadded = np.zeros(ps, in1.dtype)
130 sc = [slice(0, i) for i in in1.shape]
--> 131 in1zpadded[sc] = in1.copy()
132
133 if mode == 'full':
ValueError: could not broadcast input array from shape (66467) into shape (66466)
是什么导致了这个错误？
最佳答案
scipy.signal.cwt 的第三个参数是宽度，它必须大于 1，所以将代码更改为:
scipy.signal.cwt(np.array(c), scipy.signal.morlet, np.arange(.01,.1,.01) * len(c))
关于python - scipy.signal.cwt 在 correlate() 中出现值错误，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/23232012/