%?DOA?estimation?by?MUSIC?
%?Developed?by?xiaofei?zhang?（南京航空航天大學?電子工程系?張小飛）
%?EMAIL:zhangxiaofei@nuaa.edu.cn
clear?all
close?all
derad?=?pi/180;????????%?deg?->?rad
radeg?=?180/pi;
twpi?=?2*pi;
kelm?=?8;???????????????%?陣列數量
dd?=?0.5;???????????????%?space?
d=0:dd:（kelm-1）*dd;?????%?
iwave?=?3;??????????????%?number?of?DOA
theta?=?[20?50?60];?????%?角度
snr?=?10;???????????????%?input?SNR?（dB）
n?=?500;?????????????????%?
A=exp（-j*twpi*d.‘*sin（theta*derad））;%%%%?direction?matrix
S=randn（iwaven）;
X=A*S;
X1=awgn（Xsnr‘measured‘）;%添加噪聲
Rxx=X1*X1‘/n;%協方差
InvS=inv（Rxx）;?%%%%求逆
[EVD]=eig（Rxx）;%%%%?求矩陣Rxx的全部特征值，構成對角陣D，并求特征向量構成EV的列向量
EVA=diag（D）‘;%生成對角矩陣
[EVAI]=sort（EVA）;%則I的每一列是與EVA中列向量的元素相對應的置換向量？什么是置換矩陣
EVA=fliplr（EVA）;%左右翻轉矩陣
EV=fliplr（EV（:I））;

%?MUSIC
for?iang?=?1:361
????????angle（iang）=（iang-181）/2;
????????phim=derad*angle（iang）;
????????a=exp（-j*twpi*d*sin（phim））.‘;
????????L=iwave;????
????????En=EV（:L+1:kelm）;%如何形成噪聲子空間
????????SP（iang）=（a‘*a）/（a‘*En*En‘*a）;%為什么分母不為一
??
end
???
%?
SP=abs（SP）;
SPmax=max（SP）;
SP=10*log10（SP/SPmax）;%?
h=plot（angleSP）;%？h
set（h‘Linewidth‘2）%set的作用
xlabel（‘angle?（degree）‘）
ylabel（‘magnitude?（dB）‘）
axis（[-90?90?-60?0]）%規定坐標的范圍
set（gca?‘XTick‘[-90:30:90]）%規定坐標軸的間隔
grid?on??%網格線

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\
?????目錄???????????0??2018-03-15?09:36??陣列信號處理的理論和應用-MATLAB?程序\1.MUSIC算法MATLAB程序\
?????文件??????????22??2008-01-02?09:28??陣列信號處理的理論和應用-MATLAB?程序\1.MUSIC算法MATLAB程序\Untitled.m
?????文件????????1389??2017-11-30?20:25??陣列信號處理的理論和應用-MATLAB?程序\1.MUSIC算法MATLAB程序\doa_music1.m
?????文件????????1297??2018-03-15?09:55??陣列信號處理的理論和應用-MATLAB?程序\1.MUSIC算法MATLAB程序\doa_music_org.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\2.ESPRIT?算法MATLAB程序\
?????文件?????????675??2017-12-03?16:53??陣列信號處理的理論和應用-MATLAB?程序\2.ESPRIT?算法MATLAB程序\main.m
?????文件????????1632??2017-12-03?16:56??陣列信號處理的理論和應用-MATLAB?程序\2.ESPRIT?算法MATLAB程序\tls_esprit.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\3.Root-MUSIC算法MATLAB程序\
?????文件?????????911??2017-12-03?16:58??陣列信號處理的理論和應用-MATLAB?程序\3.Root-MUSIC算法MATLAB程序\root_music.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\
?????文件????????1312??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\Unitary_esprit.m
?????文件?????????200??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\qq.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\5.空間平滑MUSIC算法的MATLAB程序\
?????文件????????1373??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\5.空間平滑MUSIC算法的MATLAB程序\main.m
?????文件?????????240??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\5.空間平滑MUSIC算法的MATLAB程序\mssp.m
?????文件?????????174??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\5.空間平滑MUSIC算法的MATLAB程序\ssp.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\
?????文件????????1627??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\genspacetimechan.m
?????文件????????4668??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\jade.m
?????文件????????1040??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\main.m
?????文件?????????464??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\qtrans.m
?????文件?????????498??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\raisedcos_filter.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\7.傳播算子DOA估計算法MATLAB程序\
?????文件????????1513??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\7.傳播算子DOA估計算法MATLAB程序\PM.m
?????目錄???????????0??2018-03-14?15:24??陣列信號處理的理論和應用-MATLAB?程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\
?????文件?????????152??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\R_hankel.m
?????文件????????1477??2017-10-23?15:35??陣列信號處理的理論和應用-MATLAB?程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\matrix_pencil_L.m
?????文件?????????685??2010-10-25?11:20??陣列信號處理的理論和應用-MATLAB?程序\光盤說明.txt