在不使用内置FFT函数的情况下，使用Python查找功率谱

Question

我想使用以下公式使用Python查找功率谱：

功率谱

我正在尝试绘制与此有关的数组：

P = []
for k in range(0,int(N/2)):
    P.append((2/N)*(sum(x[k]*np.cos(2*np.pi*nu*t[k]))**2+(sum(x[k]*np.sin(2*np.pi*nu*t[k])))**2))

其中nu是频率：

nu = []
for j in range(0, int(N/2), 2):
nu.append(j/T)

x和t来自数据集，而N只是所有x的总和。

Python不断告诉我：“不能将序列乘以'float类型的非整数'。在这里，我在做什么错？ 我确定这与将错误类型的数据相互相乘有关。 但我不确定如何更改此设置。

我知道有一个内置的FFT函数，但是我认为让该函数正常工作对我很有帮助。

Answer 1

nu是一个列表，因此代码段2*np.pi*nu*t[k]给您错误。尝试类似的方法：

for i in nu:
    P.append((2/N)*(sum(x*np.cos(2*np.pi*nu[i]*t))**2+(sum(x*np.sin(2*np.pi*nu[i]*t)))**2))

这样，您求和的数量为(2/N)*(sum(x*np.cos(2*np.pi*nu[i]*t))**2+(sum(x*np.sin(2*np.pi*nu[i]*t)))**2)对于每个i (2/N)*(sum(x*np.cos(2*np.pi*nu[i]*t))**2+(sum(x*np.sin(2*np.pi*nu[i]*t)))**2) 。