%空時自適應處理降維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