%function?Farfield_DOA?（resignal）
clc;
clear;
close?all;
Freq?=700;????????????????????%Frequency
snr?=40;??????????????????????%Source?to?noise?ratio
M?=12;??????????????????????????%The?number?of?microphones
Cspeed?=?340;????????????????????%?The?speed?of?sound
Lambda?=?Cspeed/340;?????????????%Wavelength
d?=?Lambda/2;?????????????????????%Array?elements?spacing
Dirsources?=?[2035];?????????????%Direction?of?sound?sources
Numsources?=?length（Dirsources）;??????%Number?of?sound?sources
?
for?q?=?1:Numsources
????a（:q）=exp（-1i?*?2?*?pi?*?Freq?*?d?*（0:M-1）‘*sin（Dirsources（q）*pi/180）/Cspeed）;?????%
end
?
Energy?=?[12];????????????????%Sources?energy
y?=?a*Energy‘;?????????????????%Received?signal?of?array
Y?=?awgn（ysnr‘measured‘）;?????%Adding?noise
?
grid?=?0.1;???????????????????????%The?value?of?grid?node??
dirmin?=?-90;???????????????????%The?min?and?max?DOA?value?of?the?source?plane
dirmax?=?90;
Dirmicrogrid?=?dirmin:grid:dirmax;???????%Direction?between?microphone?and?grid?node
Numgrid?=?length（Dirmicrogrid）;??????????%Number?of?grid
?
A?=?zeros（MNumgrid）;????????????????????%Measurement?matrix?（1/sqrt（M））*e
for?j?=?1:Numgrid
????A（:j）=exp（-1i?*?2?*?pi?*?Freq?*?d?*?（0:M-1）‘?*?sin（Dirmicrogrid（j）?*?pi?/?180）?/?Cspeed）;
end
tic
resignal?=?CS_OMP（YA）?;???????%Recover?using?OMP
resignal?=?real（resignal‘）;
toc
?
nonsource?=?nonzeros（resignal）
?
figure（1）;
?
plot（Dirmicrogridresignal‘k*-‘）?????????????????????????????????%??重建信號
hold?on
plot（DirsourcesEnergy‘ks‘）
title（[‘SNR?=?‘num2str（snr）]）
xlabel（‘兩個信源方向角‘）
ylabel（‘兩個信源幅度值‘）
dirction?=?（find（resignal）?-1）?*grid?-90

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-10-31?06:56??CS_DOA\
?????文件????????1621??2017-10-31?06:50??CS_DOA\CS_DOA.m
?????目錄???????????0??2017-10-31?06:57??__MACOSX\
?????目錄???????????0??2017-10-31?06:57??__MACOSX\CS_DOA\
?????文件?????????120??2017-10-31?06:50??__MACOSX\CS_DOA\._CS_DOA.m
?????文件?????????497??2017-10-31?06:48??CS_DOA\CS_OMP.m
?????文件?????????120??2017-10-31?06:48??__MACOSX\CS_DOA\._CS_OMP.m