%空時降維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）））‘