DPCA 程序编写

%%%%%%%%%%%%%%%%% 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=