几秒钟内的 Numpy FFT

Question

我想了解 python 中的 fft 是如何工作的（或一般情况下）。
当我有一个记录了几秒钟的信号时，我只能显示一秒钟的 FFT。
有没有办法获取所有数据并将其平均？我以前在 LabView 中做过这个。有人可以帮忙吗？

这是我正在使用的示例代码。

from numpy import linspace, sin, pi, log10, average, arange
from matplotlib.pyplot import plot, show, figure
from numpy.fft import fft, fftfreq

N = 1000 * 60 * 4
dt = 1 / 1000
x = linspace(0, N*dt, N)
freq = linspace(10, 200, N)
sinsweep = sin(x*freq)
mavg = [average(sinsweep[i*60:(i+1)*60]) for i in range(int(N/60))]
plot(freq, sinsweep, '.')
plot(linspace(10, 200, int(N/60)), mavg, '.')

f = figure()
t = arange(60)
sp = fft(mavg, n=60)
freq = fftfreq(t.shape[-1])
plot(sp.imag)
show()

Answer 1

我对您的代码进行了一些修改，以获得漂亮的光谱。首先我增加点数来验证香农准则。还有一些技巧
提高速度。

from numpy import linspace, sin, pi, log10, average, arange
from matplotlib.pyplot import plot, show, figure
from numpy.fft import fft, fftfreq
close()
N = 10000 * 60 * 4
dt = 1 / 10000
t = arange(0, N*dt, dt)
freq = linspace(10, 200, N)
sinsweep = sin(t*freq)
mavg = sinsweep.reshape(-1,60).mean(1)
tm=t[::60]
figure('signal')
plot(tm, mavg)    
sp = fft(mavg)
freq = fftfreq(tm.size,dt*60)
valid=freq>0
figure('spectrum')
plot(freq[valid],abs(sp[valid]))
show()

为

这是你所期望的吗？

几秒钟内的 Numpy FFT

Numpy FFT over few seconds

python

numpy

fft