%%%%%%%%%%%%%%%%%?DPCA方法檢測動目標?（c）?鄧彬?2006.4?%%%%%%%%%%%%%%%%%%%%%%%%

%%產生Stripmap?SAR的回波
clear?all
c=3e8;%光速
%%測繪帶區域
X0=204.8;%方位向[-X0X0]，條帶長度的“一半”
Rc=10000;%條帶中心參考距離
R0=480;%距離向[Rc-R0Rc+R0]條帶寬度的“一半”
lambda=0.02;%波長
fc=c/lambda;%載頻

%%距離向（Range）r/t?domain
Tr=40e-6;%LFM信號脈寬
Br=10e6;?%LFM信號帶寬（=Kr*Tr）Hz
Kr=Br/Tr;?%調頻斜率
Fsr=40e6;%徑向快時間采樣率，應>Br
?????????%還應盡量使下面的Nr＝2若干次冪，以利于FFT。
Nr=ceil（Fsr*4*R0/c）;
r=Rc+linspace（-R0R0Nr）;
t=2*r/c;%t域序列
dt=R0*4/c/Nr;%采樣周期
f=linspace（-1/2/dt1/2/dtNr）;%f域序列

%%方位向（AzimuthCross-Range）x/u?domain
Va=200;%SAR?平臺速度
Lsar=150;%合成孔徑長度
PRF=1000;%Hz應>多普勒帶寬Ba（k*Tsar）還應盡量使下面的Na＝2若干次冪，以利于FFT
Ta=1/PRF;
D=Va*Ta;%兩孔徑間距
Na=ceil（PRF*2*X0/Va）;
x=linspace（-X0X0Na）;%u域序列
u=x/Va;
du=2*X0/Va/Na;
fu=linspace（-1/2/du1/2/duNa）;%fu域序列
fdc=0;%Doppler調頻中心頻率
fdr=-2*（Va）^2/lambda/Rc;%Doppler調頻斜率
Ba=abs（fdr*Lsar/Va）;%計算多普勒帶寬

%目標位置
Ptar=[Rc0001?%距離向坐標方位向坐標距離向速度方位向速度，sigma?????????????
??????Rc301001];
Ntar=size（Ptar1）;%目標個數
%%%%%%%%%%%%%%%?1、產生距壓、徙校后的回波?%%%%%%%%%%%
s_ut=zeros（NrNa）;
s_ut2=zeros（NrNa）;
U=ones（Nr1）*u;%擴充為矩陣
T=t‘*ones（1Na）;
for?i=1:1:Ntar
????rn=Ptar（i1）;xn=Ptar（i2）;
????vrn=Ptar（i3）;vcn=Ptar（i4）;
????sigma=Ptar（i5）;
????
????R=sqrt（（（Va-vcn）*（U-xn/Va））.^2+（rn+vrn*（U-xn/Va））.^2）;%位于x＝0處的動目標波程，王玲《兵工學報》vr正方向改為遠離。??
????phase=-4*pi*R/lambda;??
????tmp=zeros（size（T））;tmp（128:）=ones（1size（T2））;%?距離向上本應是乘Sinc函數，但可用delta函數近似，即tmptmp在128處，故目標距離均應設在Rc處
????s_ut=s_ut+sigma*exp（j*phase）.*（abs（Va*U-xn）????
????
????R2=sqrt（（（Va-vcn）*（U-xn/Va）+Va*Ta-D-vrn*Ta）.^2+（rn+vrn*（U-xn/Va）+vrn*Ta）.^2?）;
????phase2=-4*pi*R2/lambda;?
????s_ut2=s_ut2+sigma*exp（j*phase2）.*（abs（Va*U-xn）end;

%%%%%%%%%%%%%%%?第一路成像?%%%%%%%%%%%%%
FU=ones（Nr1）*fu;
p0_2fu=