正确的摆幅fft

Question

我有一个简单的鼻窦信号，但有些杂音。 如果我尝试一种简单的FFT算法，则会得到信号（23）和频率（2）的幅度。 如果我尝试使用带有Hanning窗口的Matlab进行安装，则得到的频率正确，但是幅度是错误的。 如何获得真实振幅？ 这是我正在使用的代码：

time = 0:0.01:50;
frequency = 2;
amplitude = 23;
y = amplitude * sin (2 * pi * frequency * time);
y= y + 6 * randn(size(time));
y = y - mean(y);
N = length(y);
Fs = 100;
NFFT = 2^nextpow2(N);
M = 4396;
w = hanning(M);
[Pyy,Fy] = pwelch(y, w,[],M,Fs);
plot(Fy,Pyy);

Answer 1

恕我直言，由于频谱泄漏 ，您不能这样做。 如果看到正在使用的Hanning窗口，则它是低通滤波器（居中并归一化时）。 它将降低主频率上的功率，并将引入其他频率上的功率，即平滑。 但是，它是在频域中具有无限持续时间的非周期性信号，无法计算。

figure; 
plot(w);
title(['Hanning window with ', num2str(M), 'points']);

您可以在这里阅读更多内容。