Question

我正在尝试使用我的代码计算周期图

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。

Answer 1

与this question一样，模糊的错误代码来自于对字符串数据而不是浮点数据运行numpy函数。致rstrip returns a string的电话。我不知道您的数据文件是什么样的，但在x = x.astype(np.float)之后添加x = np.array(x)可以解决您的问题。

TypeError：ufunc'splaly'不包含带签名匹配类型的循环dtype（'S64'）dtype（'S64'）dtype（'S64'）

1 个答案: