%opt2d.m:?全自由度空時(shí)自適應(yīng)處理
%--------------------------------------------------------------------------
%start?????????:?2004.11.04??AM?10:21???yunrisheng
%Latest?change?:?
%--------------------------------------------------------------------------

clear
tic
%雜波仿真參數(shù)
N?=?12;????????????????????????%?陣元個(gè)數(shù)
M?=?10;????????????????????????%?相干脈沖數(shù)
CNR?=?30;??????????????????????%?雜噪比
beta?=?1;??????????????????????%?雜波折疊系數(shù)（beta?=?2*v*T/d）
sita_a?=?-90:.9:90.;?????????%?雜波單元個(gè)數(shù)???????????????
sita?=?sita_a*pi/180;
[NN?N_bin]?=?size（sita）;
%目標(biāo)參數(shù)
sita_t?=?-25;??????????????????%?目標(biāo)DOA
omiga_t?=?0.4;?????????????????%?目標(biāo)Doppler
SNR?=?0;???????????????????????%?信噪比

%空間導(dǎo)向矢量和時(shí)間導(dǎo)向矢量
%空間頻率和Dopple頻率滿足?omiga_d?=?beta?*?omiga_s
omiga_s?=?pi*sin（sita）;??????
omiga_d?=?beta*omiga_s;???????

aN?=?zeros（NN_bin）;
bN?=?zeros（MN_bin）;

aN?=?exp（-j*[0:N-1]‘*omiga_s）./sqrt（N）;
bN?=?exp（-j*[0:M-1]‘*omiga_d）./sqrt（M）;


%目標(biāo)空時(shí)信號(hào)
aN_t?=?zeros（N1）;
bN_t?=?zeros（M1）;

aN_t?=?exp（-j*pi*[0:N-1]‘*sin（sita_t*pi/180））/sqrt（N）;
bN_t?=?exp（-j*pi*[0:M-1]‘*omiga_t）/sqrt（M）;

S_t?=?zeros（M*N1）;
S_t?=?kron（aN_tbN_t）;

%計(jì)算雜波協(xié)方差矩陣
R?=?zeros（M*NM*N）;?????????????????????
S?=?zeros（M*NN_bin）;???????????????????
ksai?=?10^（CNR/10）*（randn（1N_bin）+j*randn（1N_bin））/sqrt（2）;???????????????%服從正態(tài)分布的隨機(jī)幅值，方差為1
for?ii?=?1:N_bin
????S（:ii）?=?kron（aN（:ii）bN（:ii））;??
????R?=?R?+?ksai（ii）.*（S（:ii）*S（:ii）‘）;???????
end

%干擾協(xié)方差矩陣，雜噪比為30dB
R?=?R?+eye（M*N）;??