STAP仿真程序（全自由度空时自适应处理）

%opt2d.m: 全自由度空时自适应处理
%--------------------------------------------------------------------------
%start         : 2004.11.04  AM 10:21   yunrisheng
%Latest change : 
%--------------------------------------------------------------------------

clear
tic
%杂波仿真参数
N = 12;                        % 阵元个数
M = 10;                        % 相干脉冲数
CNR = 30;                      % 杂噪比
beta = 1;                      % 杂波折叠系数（beta = 2*v*T/d）
sita_a = -90:.9:90.;         % 杂波单元个数               
sita = sita_a*pi/180;
[NN N_bin] = size（sita）;
%目标参数
sita_t = -25;                  % 目标DOA
omiga_t = 0.4;                 % 目标Doppler
SNR = 0;                       % 信噪比

%空间导向矢量和时间导向矢量
%空间频率和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）;


%目标空时信号
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）;

%计算杂波协方差矩阵
R = zeros（M*NM*N）;                     
S = zeros（M*NN_bin）;                   
ksai = 10^（CNR/10）*（randn（1N_bin）+j*randn（1N_bin））/sqrt（2）;               %服从正态分布的随机幅值，方差为1
for ii = 1:N_bin
    S（:ii） = kron（aN（:ii）bN（:ii））;  
    R = R + ksai（ii）.*（S（:ii）*S（:ii）‘）;       
end

%干扰协方差矩阵，杂噪比为30dB
R = R +eye（M*N）;