exa6_18.m

matlab应用实例。如：FIR、IIR数字滤波器设计等。

%exa060307_pmcov.m, for example 6.3.7
%to test pmcov.m;

clear all;
Fs=1000; % 采样频率
%信号序列产生；
n=0:1/Fs:.3;
w0=200*pi;
w1=400*pi;
xn=cos(w0*n)+sin(w1*n)+randn(size(n));

% 绘制信号波形
subplot(211)
plot(n,xn)
xlabel('时间(s)')
ylabel('幅度')
title('x(n)=cos(w0*n)+sin(w1*n)+e(n)')
ymax_xn=max(xn)+0.2;
ymin_xn=min(xn)-0.2;
axis([0 0.3 ymin_xn ymax_xn])

% 使用改进协方差方法估计序列功率谱；
p=floor(length(xn)/3)+1;
nfft=1024;
[xpsd,f]=pmcov(xn,p,nfft,Fs,'half');

% 绘制功率谱估计
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(212)
plot(f,xpsd);
title('基于改进协方差方法的功率谱估计')
ylabel('功率谱密度(dB)')
xlabel('频率(Hz)')
grid on;
ymin = min(xpsd)-2;
ymax = max(xpsd)+2;
axis([0 Fs/2 ymin ymax])