exa6_22.m

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

%exa060401_MUSIC.m, for example 6.4.1
%to test MUSIC.m;

clear all;
%信号序列产生；
randn('state',1); 
N=1000; %采样点数
n=0:N-1;  
en=randn(1,N);
xn=exp(j*pi/2*n)+2*exp(j*pi/4*n)+exp(j*pi/3*n)+en;  

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

%用修正协方差方法估计相关矩阵
X=corrmtx(xn,7,'mod'); 

%用MUSIC算法估计序列的功率谱
subplot(212)
pmusic(X,2);
title('基于MUSIC算法的功率谱估计')
ylabel('功率谱幅度(dB)')
xlabel('归一化频率(rad/sample)')
grid on;