exa6_9.m

新版Matlab 7辅助信号处理技术与应用

% exa060204_psd_Bartlett.m,  for example 6.2.5 
% to test psd_Bartlett.m  ;

clear all;
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(311)
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)')

%参数设置
nfft=1024;
window=boxcar(1001);
noverlap=0;
p=0.9;
%计算序列的PSD
[Pxx,Pxxc]=psd(xn,nfft,Fs,window,noverlap,p);
%绘制图形
index=0:round(nfft/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
plot_Pxxc=10*log10(Pxxc(index+1));
%功率谱的对数曲线
subplot(312)
plot(k,plot_Pxx)
ylabel('功率谱的对数曲线')
%置信区间中的功率谱曲线
subplot(313)
plot(k,[plot_Pxx plot_Pxx-plot_Pxxc plot_Pxx+plot_Pxxc])
xlabel('置信区间中的功率谱曲线')