空时自适应处理3DT算法仿真

%空时降维3dt方法
%2012-2-27

clc;clear all;close all;
load clutter_matrix.mat;

[NKL]=size（clutter_matrix）;
N=16;
K=10;
CNR=40;
Rc=clutter_matrix*clutter_matrix‘/L;
noise=max（max（Rc））/（10^（CNR/10））*eye（N*K）;
Rx=Rc+noise;
anoise=max（max（noise））;                          %噪声功率
Qs=eye（N）;                                       %空域降维矩阵
psi0=pi/2;

% 全维STAP
fd=-1:1/50:1;
inv_Rx=inv（Rx）;
for i=1:length（fd）;
    Ss=exp（j*pi*（0:N-1）‘*cos（psi0））;            %目标方向确定时，Ss固定，但目标doppler频率未知，故每一个fd有一个最优权矢量wopt
    St=exp（j*pi*（0:K-1）‘*fd（i））;
    S=kron（StSs）;
    wopt=inv_Rx*S/（S‘*inv_Rx*S）;
    IF（i）=abs（wopt‘*S）^2*（10^（CNR/10）+1）*anoise/（wopt‘*Rx*wopt）;
end
figure
plot（fd10*log10（abs（IF）））
xlabel（‘2f_d/f_r‘）;ylabel（‘IF/dB‘）;
grid on


% index=0;
% for fd=-1:1/50:1
%     index=index+1;
%     Qt=[exp（j*pi*（0:K-1）.‘*（fd-1/K））exp（j*pi*（0:K-1）.‘*fd）exp（j*pi*（0:K-1）.‘*（fd+1/K））];%时域降维矩阵
%     Q=kron（QtQs）;
%     Ry=Q‘*Rx*Q;
%     inv_Ry=inv（Ry）;
%     Ss=exp（j*pi*（0:N-1）.‘*cos（psi0））;
%     St=exp（j*pi*（0:K-1）.‘*fd）;
%     g1=（exp（j*pi*（0:K-1）.‘*（fd-1/K）））‘