机载雷达STAP降维3DT算法

%空时自适应处理降维3DT算法
clc;clear all;close all;
load clutter_matrix.mat;

[NKL]=size（clutter_matrix）;
N=16;
K=10;
CNR=60;
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