%?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?=?[10?30?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）;%%%%?
EVA=diag（D）‘;
[EVAI]=sort（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）;
set（h‘Linewidth‘2）
xlabel（‘angle?（degree）‘）
ylabel（‘magnitude?（dB）‘）
axis（[-90?90?-60?0]）
set（gca?‘XTick‘[-90:30:90]）
grid?on??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-06-24?23:51??陣列信號處理的理論和應用\
?????目錄???????????0??2013-09-07?11:13??陣列信號處理的理論和應用\程序\
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\1.MUSIC算法MATLAB程序\
?????文件????????1122??2010-10-25?11:20??陣列信號處理的理論和應用\程序\1.MUSIC算法MATLAB程序\doa_music.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\2.ESPRIT?算法MATLAB程序\
?????文件?????????674??2010-10-25?11:20??陣列信號處理的理論和應用\程序\2.ESPRIT?算法MATLAB程序\main.m
?????文件????????1584??2010-10-25?11:20??陣列信號處理的理論和應用\程序\2.ESPRIT?算法MATLAB程序\tls_esprit.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\3.Root-MUSIC算法MATLAB程序\
?????文件?????????908??2010-10-25?11:20??陣列信號處理的理論和應用\程序\3.Root-MUSIC算法MATLAB程序\root_music.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\
?????文件?????????200??2010-10-25?11:20??陣列信號處理的理論和應用\程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\qq.m
?????文件????????1312??2010-10-25?11:20??陣列信號處理的理論和應用\程序\4.面陣中二維角度估計?Unitary?-ESPRIT算法MATLAB程序\Unitary_esprit.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\5.空間平滑MUSIC算法的MATLAB程序\
?????文件????????1373??2010-10-25?11:20??陣列信號處理的理論和應用\程序\5.空間平滑MUSIC算法的MATLAB程序\main.m
?????文件?????????240??2010-10-25?11:20??陣列信號處理的理論和應用\程序\5.空間平滑MUSIC算法的MATLAB程序\mssp.m
?????文件?????????174??2010-10-25?11:20??陣列信號處理的理論和應用\程序\5.空間平滑MUSIC算法的MATLAB程序\ssp.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\
?????文件????????1627??2010-10-25?11:20??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\genspacetimechan.m
?????文件????????4668??2010-10-25?11:20??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\jade.m
?????文件????????1040??2010-10-25?11:20??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\main.m
?????文件?????????464??2010-10-25?11:20??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\qtrans.m
?????文件?????????498??2010-10-25?11:20??陣列信號處理的理論和應用\程序\6.角度和時延聯合估計（JADE）算法MATLAB程序\raisedcos_filter.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\7.傳播算子DOA估計算法MATLAB程序\
?????文件????????1513??2010-10-25?11:20??陣列信號處理的理論和應用\程序\7.傳播算子DOA估計算法MATLAB程序\PM.m
?????目錄???????????0??2013-09-07?11:04??陣列信號處理的理論和應用\程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\
?????文件????????1477??2010-10-25?11:20??陣列信號處理的理論和應用\程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\matrix_pencil_L.m
?????文件?????????152??2010-10-25?11:20??陣列信號處理的理論和應用\程序\8.基于增廣矩陣束的L型陣列的二維DOA估計MATLAB程序\R_hankel.m
?????文件?????????685??2010-10-25?11:20??陣列信號處理的理論和應用\程序\光盤說明.txt
?????文件????21786907??2019-06-24?23:51??陣列信號處理的理論和應用\陣列信號處理的理論和應用.pdf