longteng4.m

Welch法：Welch法对Bartlett法进行了两方面的修正：一是选择适当的窗函数w(n)

%Welch法对Bartlett法进行了两方面的修正：
%一是选择适当的窗函数w(n)，并在周期图计算前直接加进去，
%加窗的优点是无论什么样的窗函数均可使谱估计非负。
%二是在分段时,可使各段之间有重叠，这样会使方差减小
% author:闫朝星 2007.12.18
% version 0.1
%---------------------------------
Fs=600;
n=0:1/Fs:1;
xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+randn(size(n));
nfft=512;
window=boxcar(100);%矩形窗
window1=hamming(100);%海明窗
window2=blackman(100);%blackman窗
noverlap=20;%数据无重叠
range='half';%频率间隔为[0 Fs/2],计算一半的频率
[Pxx,f]=pwelch(xn,window,noverlap,nfft,Fs,range);
[Pxx1,f]=pwelch(xn,window1,noverlap,nfft,Fs,range);
[Pxx2,f]=pwelch(xn,window2,noverlap,nfft,Fs,range);
plot_Pxx=10*log10(Pxx);
plot_Pxx1=10*log10(Pxx1);
plot_Pxx2=10*log10(Pxx2);
figure(1)
plot(f,plot_Pxx);
figure(2)
plot(f,plot_Pxx1);
figure(3)
plot(f,plot_Pxx2);