与窗口FFT振幅校正

180(30天)

显示旧的评论

CduP 2020年4月9日

1
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/516655-fft-with-windowing-amplitude-correction

评论道: 艾略特马丁森2021年5月17日

我想计算一个压力信号的FFT与正确的振幅,正如我想从FFT计算产生的声压级。

基本上,我当我能正确地计算FFT振幅不使用一个窗口,但我想用一个汉宁窗,我不知道如何规范化的窗口。在下面的代码中,这篇文章 //www.tianjin-qmedu.com/matlabcentral/answers/33653-psd-estimation-fft-vs-welch, FFT是由窗口规范化权力通过FFT /(窗口‘*窗口)。然而,FFT长度除以Nx从来没有规范化的信号。但是,正确的PSD是发现。我能理解,结果从pwelch之前也应该扩展可用于获得正确的声压水平。

                         Fs = 1024;
                        
                         t = (0:1 / Fs: 1 - 1 / Fs)。”;
                        
                         x =罪(2 *π* t * 200);
                        
                         Nx =长度(x);
                        
                         %窗口数据
                        
                         w =汉宁(Nx);
                        
                         xw = x。* w;
                        
                         %计算能力
                        
                         nfft = Nx;
                        
                         X = fft (xw nfft);
                        
                         mx = abs (X) ^ 2;
                        
                         %正常化的窗口。乘以2 (DC &尼奎斯特除外)
                        
                         %计算片面的频谱。除以Fs计算
                        
                         %谱密度。
                        
                         mx = mx / (w * w);
                        
                         NumUniquePts = nfft / 2 + 1;
                        
                         mx = mx (1: NumUniquePts);
                        
                         mx (2: end-1) = mx (2: end-1) * 2;
                        
                         Pxx1 = mx / Fs;
                        
                         Fx1 = (0: NumUniquePts-1) * Fs / nfft;
                        
                         [Pxx2, Fx2] = pwelch (x, w, 0 nfft Fs);
                        
                         情节(Fx1 10 * log10 (Pxx1) Fx2, 10 * log10 (Pxx2),“:”);
                        
                         传奇(“通过FFT PSD”,“通过pwelch PSD”)