校正线性距离走动的距离多普勒算法.doc

%校正线性距离走动的距离多普勒算法%注：波束射线有一定的斜视角，距离弯曲的影响仍然可以忽略，但距离走动的影响必须加以考虑%%%%%%%%%%%%%%%%%%%%%%clear;clc;close all;%%========================================================C=3e8; %电磁波传播速度Fc=1e9; %载频lambda=C/Fc; %波长 Xmin=0; %目标方位区域[Xmin,Xmax]，地面沿航线方向Xmax=50; Yc=10000; %条带中心，地面垂直航线距离Y0=500; %条带宽度区域[Yc-Y0,Yc+Y0] V=100; %雷达平台速度H=0; %雷达平台高度Rb=sqrt(Yc^2+H^2);%条带中心到航线的垂直距离 xsj=30/180*pi; %天线的斜视角,弧度表示D=4; %天线方位向有效长度Lsar=Rb*(tan(xsj+1/2*(lambda/D))-tan(xsj-1/2*(lambda/D))); %合成孔径长度Tsar=Lsar/V; %合成孔径时间 R0=Rb*sec(xsj);Ka=-2*V^2*cos(xsj).^2/lambda/R0; %多普勒调制频率Ba=abs(Ka*Tsar); %多普勒带宽PRF=2*Ba; PRT=1/PRF; ds=PRT; %慢时间域的抽样间隔Nslow=ceil((Xmax-Xmin+Lsar)/V/ds); %慢时间域的抽样数Nslow=2^nextpow2(Nslow); sn=linspace((Xmin-Rb*(tan(xsj+1/2*(lambda/D))-tan(xsj)))/V,(Xmax+Rb*(tan(xsj)-tan(xsj-1/2*(lambda/D))))/V,Nslow);PRT=(Xmax-Xmin+Lsar)/V/Nslow; %更新xmin=Xmin-Rb*(tan(xsj+1/2*(lambda/D))-tan(xsj));xmax=Xmax+Rb*(tan(xsj)-tan(xsj-1/2*(lambda/D)));L1=-Rb*(tan(xsj+1/2*(lambda/D))-tan(xsj));L2=Rb*(tan(xsj)-tan(xsj-1/2*(lambda/D)));PRF=1/PRT;ds=PRT; Tr=5e-6; %发射脉冲宽度Br=30e6; %发射脉冲带宽Kr=Br/Tr; Fsr=2*Br; %快时间域的抽样频率dt=1/Fsr; if xsj=1/2*(lambda/D) Rmin=sqrt((Yc-Y0)^2+H^2);else Rmin=sqrt(((Yc-Y0)^2+H^2)*(1+tan(xsj-1/2*(lambda/D)).^2));endRmax=sqrt(((Yc+Y0)^2+H^2)*(1+tan(xsj+1/2*(lambda/D)).^2));Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt);%快时间域抽样数Nfast=2^nextpow2(Nfast); tm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast);dt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast; %更新Fsr=1/d