%?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??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????24576??2012-12-23?14:32??光盤說明.doc

?????文件????2738917??2013-06-09?09:39??陣列信號處理的理論與應用（第2版）中仿真圖.pdf

?????文件???????1122??2010-08-16?09:41??MATALB?程序\1.MUSIC算法MATLAB程序\doa_music.m

?????文件???????1233??2012-12-02?13:20??MATALB?程序\10.L型陣下基于2D-MUSIC的二維DOA估計算法\main_2d_doa_music.m

?????文件????????177??2012-12-02?13:35??MATALB?程序\11.一階前向預測的DOA估計算法MATLAB程序\array_manifold.m

?????文件???????1663??2012-12-02?13:39??MATALB?程序\11.一階前向預測的DOA估計算法MATLAB程序\lp_doaestimation.m

?????文件????????605??2012-12-02?13:34??MATALB?程序\11.一階前向預測的DOA估計算法MATLAB程序\receive_signal.m

?????文件????????294??2012-12-02?13:31??MATALB?程序\11.一階前向預測的DOA估計算法MATLAB程序\steering_vector.m

?????文件???????2131??2012-12-02?13:47??MATALB?程序\12.?寬帶信號DOA估計的非相干信號子空間（ISM）MATLAB程序\ISM_code.m

?????文件???????2138??2012-12-02?14:06??MATALB?程序\13.信源數估計MDL算法MATLAB程序\mdl_sourcenumber.m

?????文件????????959??2012-12-02?14:20??MATALB?程序\14.均勻線陣DOA估計Unitary?-ESPRIT算法MATLAB程序\test_unitary_esprit.m

?????文件???????2649??2012-12-02?14:20??MATALB?程序\14.均勻線陣DOA估計Unitary?-ESPRIT算法MATLAB程序\unitary_esprit.m

?????文件???????1673??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\AssignAgain.m

?????文件????????546??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\ColToRow.m

?????文件????????526??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\ConjQVector.m

?????文件????????693??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\FVeToQVe.m

?????文件????????502??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\hpc.m

?????文件????????981??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\IsoMatrix.m

?????文件???????2978??2012-12-03?10:00??MATALB?程序\15?四元數MUSIC的MATLAB程序\main.m

?????文件????????525??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\PolarSource.m

?????文件????????813??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\QuatMuti.m

?????文件????????796??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\SteerVector.m

?????文件????????822??2012-12-03?09:58??MATALB?程序\15?四元數MUSIC的MATLAB程序\VectorMulti.m

?????文件???????1438??2012-12-10?10:19??MATALB?程序\16.LCMV波束形成算法MATLAB程序\lcmv_beamform.m

?????文件???????1950??2012-12-10?10:15??MATALB?程序\17.LMS自適應波束形成MATLAB程序\Lms_beamform.m

?????文件????????659??2012-06-19?14:18??MATALB?程序\18.陣列信號處理中DOA估計演示平臺MATLAB程序\doa_music.m

?????文件???????2130??2012-10-14?10:13??MATALB?程序\18.陣列信號處理中DOA估計演示平臺MATLAB程序\frame.fig

?????文件???????3501??2012-06-19?14:18??MATALB?程序\18.陣列信號處理中DOA估計演示平臺MATLAB程序\frame.m

?????文件?????????97??2012-06-19?14:18??MATALB?程序\18.陣列信號處理中DOA估計演示平臺MATLAB程序\khatri_rao.m

?????文件???????4861??2012-10-14?10:13??MATALB?程序\18.陣列信號處理中DOA估計演示平臺MATLAB程序\onedimensional.fig

............此處省略48個文件信息