exa6_10.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 33 行

% exa060206_pwelch.m, for example 6.2.6 
% to estimate the PSD of x(n) by Welch average;

clear;
Fs=1000;% 采样频率
% 产生含有噪声的序列
var=sqrt(1/exp(1.0));
n=0:1/Fs:1;
N=length(n);
e=var*randn(1,N);
w0=100*pi;
w1=50*pi;
xn=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e;

% 绘制信号波形
subplot(211)
plot(n,abs(xn))
xlabel('n')
title('x(n)=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e(n)')

% 用 Welch 平均估计序列的功率谱；
nfft=1024;
window=hamming(100);
noverlap=20;
range='half';
[Pxx,f] = pwelch(xn,window,noverlap,nfft,Fs,range);
plot_Pxx=10*log10(Pxx);
subplot(212)
plot(f,plot_Pxx)
title('Welch 法平均功率谱估计')