MUSIC算法源码

clear all
close all
derad=pi/180;%角度变弧度
radeg=180/pi;
twpi=2*pi;
kelm=8;%阵元数
dd=0.5;%阵元间距
d=0:dd:（kelm-1）*dd;
iwave=3;%信源数
theta=[10 30 60];%入射角度
snr=6;%信噪比
n=500;%采样数
A=exp（-j*twpi*d.‘*sin（theta*derad））;%方向矢量
S=randn（iwaven）;%信源信号
X=A*S;%接收信号
X1=awgn（Xsnr‘measured‘）;%添加噪声
Rxx=X1*X1‘/n;%计算协方差矩阵
[EVD]=eig（Rxx）;%特征值分解
EVA=diag（D）‘;
[EVAI]=sort（EVA）;%特征值从小到大排序
EVA=fliplr（EVA）;%左右翻转，从大到小排序
EV=fliplr（EV（:I））;%对应特征矢量