【发布时间】:2017-01-11 08:43:31
【问题描述】:
我正在尝试使用我的代码计算周期图
from scipy import signal
import numpy as np
import matplotlib.pyplot as plt
x = [line.rstrip('\n') for line in open('27000.dat')]
x = np.array(x)
fs=64
f, Pxx_den = signal.periodogram(x, fs)
plt.semilogy(f, Pxx_den)
plt.xlabel('frequency [Hz]')
plt.ylabel('PSD [V**2/Hz]')
plt.show()
但我得到了
File "m3.py", line 9, in <module>
f, Pxx_den = signal.periodogram(x, fs)
File "/home/milenko/miniconda2/lib/python2.7/site-packages/scipy/signal/spectral.py", line 141, in periodogram
scaling, axis)
File "/home/milenko/miniconda2/lib/python2.7/site-packages/scipy/signal/spectral.py", line 273, in welch
return_onesided, scaling, axis)
File "/home/milenko/miniconda2/lib/python2.7/site-packages/scipy/signal/spectral.py", line 391, in csd
mode='psd')
File "/home/milenko/miniconda2/lib/python2.7/site-packages/scipy/signal/spectral.py", line 824, in _spectral_helper
scale = 1.0 / (fs * (win*win).sum())
TypeError: ufunc 'multiply' did not contain a loop with signature matching types dtype('S64') dtype('S64') dtype('S64')
为什么? 我在 Ubuntu 16.04,Python 2.7.11 上。
【问题讨论】:
标签: python