music.m

来自「阵列天线的DOA估计MUSIC算法」· M 代码 · 共 14 行

%music算法以自相关矩阵分析为基础，将系统自相关矩阵分为两个子空间：信号子空间和噪声子空间。
%为求功率谱估计，music算法计算信号空间和噪声空间的
%特征值向量函数，使得在正弦信号频率处功率谱现峰值，而其他地方函数最小。
例子：用music法估计信号x(t)=sin(2pif1t)+2sin(2pif2t)+w(t)的功率谱。其中f1=50Hz，
f2=120Hz,w(t)为白噪声，采样频率Fs=1000Hz.
程序如下：
N=1000;
n=0:N-1;
Fs=1024;
t=n/Fs;
xn=sin(2*pi*100*t)+2*sin(2*pi*200*t)＋randn(1,N);
subplot(211);
pmusic(xn,7)