双通道dpca程序

%%产生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*Lsa