rd算法

《rd算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、%%产生StripmapSAR的回波clearallthetaT=0;%T平台波束斜视角thetaT=thetaT*pi/180;%radthetaR=0;%R平台波束斜视角thetaR=thetaR*pi/180;c=3e8;%光速fc=1.5e9;%载频lambda=c/fc;%波长%%测绘带区域X0=200;%方位向[-X0,X0]Rtc=3000;Rrc=3000;Rc=(Rtc+Rrc)/2;R0=150;%距离向[Rc-R0,Rc+R0]%%距离向(Range),r/tdomainTr=1.5e-6;%LFM信号脉宽1.5us(200m)

2、Br=150e6;%LFM信号带宽150MHzKr=Br/Tr;%调频斜率Nr=512;r=Rc+linspace(-R0,R0,Nr);t=2*r/c;%t域序列dt=R0*4/c/Nr;%采样周期f=linspace(-1/2/dt,1/2/dt,Nr);%f域序列%%方位向(Azimuth,Cross-Range),x/udomainv=100;%SAR平台速度Lsar=300;%合成孔径长度Na=1024;x=linspace(-X0,X0,Na);%u域序列u=x/v;du=2*X0/v/Na;fu=linspace(-1/2/du,1/2

3、/du,Na);%fu域序列ftdc=v*sin(thetaT);ftdr=-(v*cos(thetaT))^2/lambda/Rtc;frdc=v*sin(thetaR);frdr=-(v*cos(thetaR))^2/lambda/Rrc;fdc=ftdc+frdc;%Doppler调频中心频率fdr=ftdr+frdr;%Doppler调频斜率%％目标位置Ntar=3;%目标个数Ptar=[Rrc,0,1%距离向坐标,方位向坐标,sigmaRrc+50,-50,1Rrc+50,50,1];%%产生回波s_ut=zeros(Nr,Na);U=on

4、es(Nr,1)*u;%扩充为矩阵T=t'*ones(1,Na);fori=1:1:Ntarrn=Ptar(i,1);xn=Ptar(i,2);sigma=Ptar(i,3);rtn=rn+Rtc-Rrc;RT=sqrt(rtn^2+(rtn*tan(thetaT)+xn-v*U).^2);RR=sqrt(rn^2+(rn*tan(thetaT)+xn-v*U).^2);R=RT+RR;DT=T-R/c;phase=pi*Kr*DT.^2-2*pi/lambda*R;s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)

5、2).*(abs(v*U-xn)

6、tshift(fft(fftshift(src_ut).')).';%距离多普勒域%%二次距离压缩,距离迁移校正原理仿真src_fuf=fftshift(fft(fftshift(src_uf).')).';%距离压缩后的二维频谱F=f'*ones(1,Na);%扩充为矩阵FU=ones(Nr,1)*fu;p0_2f=exp(j*pi/fc^2/fdr*(FU.*F).^2+j*pi*fdc^2/fc/fdr*F-j*pi/fc/fdr*FU.^2.*F);s2rc_fuf=src_fuf.*p0_2f;s2rc_fut=fftshift(ifft(

7、fftshift(s2rc_fuf)));%距离多普勒域%%方位压缩p0_2fu=exp(j*pi/fdr*(FU-fdc).^2);%方位向压缩因子s2rcac_fut=s2rc_fut.*p0_2fu;%方位压缩s2rcac_fuf=fftshift(fft(fftshift(s2rcac_fut)));%距离方位压缩后的二维频谱s2rcac_ut=fftshift(ifft(fftshift(s2rcac_fut).')).';%方位向IFFTsubplot(221)G=20*log10(abs(s_ut)+1e-6);gm=max(max(G

8、));gn=gm-40;%显示动态范围40dBG=255/(gm-gn)*(G-gn).*(G>gn);im