8阵元均匀圆阵MUSIC DOA算法的matlab程序代码

%%--------------------8阵元圆阵DOA MUSIC算法---------------------%%
%%             开始时间：2012.7.19                               %%
%%             修改时间：2012.7.20                               %%
%%--------------------------------------------------------------%%

clear;
clc;
%采样点数
N = 100;
%空域天线个数
M = 8;
%信源个数
K = 3;

%入射波的方向
Direction1 = 100 * pi/180; 
Direction2 = 5 * pi/180;
Direction3 = 10 * pi/180;
Direction = [Direction1;Direction2;Direction3];

%入射波的信噪比
SNR1 = 30;
SNR2 = 40;
SNR3 = 35;

%入射波的幅值，这里假设噪声的功率为1
A1 = 10^（SNR1/20）;
A2 = 10^（SNR2/20）;
A3 = 10^（SNR3/20）;

%-----------------假设入射信号为正弦波--------------%%
%入射波频率
f1 = 0.1;
f2 = 0.2;
f3 = 0.3;

%入射波的随机相位这里假设其服从0带2pi平均分布
phase1 = 2*pi*rand;
phase2 = 2*pi*rand;
phase3 = 2*pi*rand;

%入射信号表达式
%复指数形式
signal1 = A1 * exp（j*（2*pi*f1*（0:N-1）+phase1））;
signal2 = A2 * exp（j*（2*pi*f2*（0:N-1）+phase2））;
signal3 = A3 * exp（j*（2*pi*f3*（0:N-1）+phase3））;
signal = [ signal1; signal2; signal3];
%%-------------------------------------------%%

%%-----------------构造方向向量---------------%%
Radius_MUS = 0.8;    %八角阵元半径

d = 0.7654 * Radius_MUS   %阵元间距

c = 3 * 10^8 ; % 光速
f1_radio = 2.4 * 10^9;  %2.4GHz频率
f2_radio = 5.8 * 10^9;  %5.8GHz频率
f600M_radio = 6 * 10^8;
f800M_radio = 8 * 10^8;

wavelength1 = c/f1_radio; 
wavelength2 = c/f2_radio;
wavelength600M = c/f600M_radio; 
wavelength800M = c/f800M_radio;

wavelengt