clc;
clear?all;
close?all;
%%?ture?DOD?and?DOA
DOD?=?[-40??-20??20???20];
DOA?=?[?30??-50??40??-50];
%?DOD?=?[-40??30?20];
%?DOA?=?[?30??30?50];
DOD_DOA?=?DOD?+?1i?*?DOA;???????????%?the?ture?location
%%?parameters
SNR?=?0;
Mt?=?3;?????????????????????????????%?transmitting?array
Mr?=?3;?????????????????????????????%?recieving?array
K?=?length（DOA）;????????????????????%?mumber?of?targets
Num_of_pulse?=?100;?????????????????%?the?number?of?independent?pulses
T?=?512;????????????????????????????%?the?code?length?per?pulse
%%???randomly?generating?the?reflection?coefficient?of?each?target??
Beta?=?zeros（Num_of_pulseK）;?????????????????????%?the?RCS?of?targets?in?different?pulses
std_beta?=?sqrt（0.5?+?0.5?*?rand（1K））;???????????%?the?standard?deviation?of?RCS
for?k?=?1?:?K
????Beta（:k）?=?std_beta（k）?*?（randn（Num_of_pulse1）?+?1i?*?randn（Num_of_pulse1））?/?sqrt（2）;????%?the?RCS?of?the?targets?（maybe?complex?valued）
end
%%?the?transmitted?waveforms
H?=?hadamard（T）;
S?=?（1?+?1i）/sqrt（2）?*?H（1:Mt:）;?????????????????????????%?orthogonal?transmiting?signals
%%??the?recieved?data?and?matched?filtering

theta?=?DOD?*?pi?/?180;???????????????????????????????????%?converting?to?radians
A_theta?=?exp（1i?*?pi?*?[0:1:Mt-1]‘?*?sin（theta））;???????%?half?wavelength?is?assumed
phi?=?DOA?*?pi?/?180;?????????????????????????????????????%?converting?to?radians
B_phi?=?exp（1i?*?pi?*?[0:1:Mr-1]‘?*?sin（phi））;???????????%?half?wavelength?is?assumed

R?=?zeros（MrMtNum_of_pulse）;?
for?n?=?1?:?Num_of_pulse?
????beta?=?Beta（n:）;
????Xn?=?B_phi?*?diag（beta）?*?conj（A_theta‘）?*?S;

????%?---------------?Receiver?Noise?Mode----------------
%??????sigma?=?sqrt（mean（var（Xn.‘）））?*?10^（-SNR/20）;?????????????????????????%?std?of?the?noise
%??????Noise?=?sigma?*?（randn（size（Xn））?+?1i?*?randn（size（Xn）））?/?sqrt（2）;???%?sensor?noise
%??????Yn?=?Xn?+?Noise;????????????
%??????R（::n）?=?Yn?*?S‘/T;?????????????????????????????????????????????????%?matched?d