%%--------------------8陣元圓陣DOA?MUSIC算法---------------------%%
%%?????????????開始時間：2012.7.19???????????????????????????????%%
%%?????????????修改時間：2012.7.20???????????????????????????????%%
%%--------------------------------------------------------------%%

clear;
clc;
%采樣點數(shù)
N?=?100;
%空域天線個數(shù)
M?=?8;
%信源個數(shù)
K?=?3;

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

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

%入射波的幅值，這里假設(shè)噪聲的功率為1
A1?=?10^（SNR1/20）;
A2?=?10^（SNR2/20）;
A3?=?10^（SNR3/20）;

%-----------------假設(shè)入射信號為正弦波--------------%%
%入射波頻率
f1?=?0.1;
f2?=?0.2;
f3?=?0.3;

%入射波的隨機相位這里假設(shè)其服從0帶2pi平均分布
phase1?=?2*pi*rand;
phase2?=?2*pi*rand;
phase3?=?2*pi*rand;

%入射信號表達式
%復(fù)指數(shù)形式
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];
%%-------------------------------------------%%

%%-----------------構(gòu)造方向向量---------------%%
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