exa6_21.m

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

%exa060310_MVSE.m, for example 6.3.10
%to test MVSE.m;

clear all;
%信号序列产生；
N=1000;
n=1:N;
w0=1.4*pi;
var=sqrt(1/exp(1.0));
e=var*randn(1,N);
xn=exp(j*pi*n-j*pi)+exp(j*w0*n-j*0.7*pi)+e;

% 绘制信号波形
subplot(211)
plot(n,xn)
xlabel('数据采样点n')
ylabel('幅度')
title('x(n)=exp(j*pi*n-j*pi)+exp(j*w0*n-j*0.7*pi)+e(n)')

% 估计自相关函数；
m=-500:500;
[r,lag]=xcorr(xn,500,'biased');

R=[r(501) r(502) r(503) r(504);
   r(500) r(501) r(502) r(503);
   r(499) r(500) r(501) r(502);
   r(498) r(499) r(500) r(501)];

[V,D]=eig(R);
V3=[V(1,3),V(2,3),V(3,3),V(4,3)].';
V3=[V(1,4),V(2,4),V(3,4),V(4,4)].';

p=0:3;
wm=[0:0.002*pi:2*pi];
B=[(exp(-j)).^(wm'*p)];
A=B.';

%最小方差功率谱估计
z=A'*inv(R)*A;
Z=diag(z');
pmv=1./Z;
subplot(212)
plot(wm/pi,pmv);
title('基于最小方差的功率谱估计')
ylabel('功率谱幅度(dB)')
xlabel('角度频率,w/pi')
grid on;