【正文】 IFFTsubplot(221)G=20*log10(abs(s_ut)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)的方位角39。)ylabel(39。成像范圍39。)title(39。(圖 一)原始信號39。)subplot(222)G=20*log10(abs(src_fut)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(fu,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)的方位角39。)ylabel(39。成像范圍39。)title(39。(圖 二)距離多普勒域頻譜39。)subplot(223)G=20*log10(abs(s2rc_fut)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(fu,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)的方位角39。)ylabel(39。成像范圍39。)title(39。(圖 三)RMC后的RD域頻譜39。)subplot(224)G=20*log10(abs(s2rcac_ut)+1e6)。gm=max(max(G))。gn=gm60。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)的方位角39。)ylabel(39。成像范圍39。)title(39。(圖 四) 目標(biāo)圖象39。)%%CS算法產(chǎn)生Stripmap SAR的回波thetaT=0。 %T平臺波束斜視角 thetaT=thetaT*pi/180。 %radthetaR=0。 %R平臺波束斜視角 thetaR=thetaR*pi/180。c=3e8。 %光速 a=。 %載頻lambda=c/a。 %波長%%測繪帶區(qū)域X0=200。 %方位向[X0,X0] Rtc=3000。 Rrc=3000。 Rc=(Rtc+Rrc)/2。R0=150。 %距離向[RcR0,Rc+R0]%%距離向(成像范圍),r/t domainTr=。 % (200m)Br=150e6。 %LFM信號帶寬 150MHzKr=Br/Tr。 %調(diào)頻斜率Nr=1024。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域序列%%方位向(成像目標(biāo)方位角,Cross成像范圍),x/u domainv=100。 %SAR 平臺速度Lsar=300。 %合成孔徑長度Na=512。 x=linspace(X0,X0,Na)。%u域序列u=x/v。 du=2*X0/v/Na。fu=linspace(1/2/du,1/2/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調(diào)頻中心頻率fdr=ftdr+frdr。 %Doppler調(diào)頻斜率%％目標(biāo)位置Ntar=3。 %目標(biāo)個數(shù)Ptar=[Rrc,0,1 %距離向坐標(biāo),方位向坐標(biāo),sigma Rrc+50,50,1 Rrc+50,50,1]。%%產(chǎn)生回波s_ut=zeros(Nr,Na)。U=ones(Nr,1)*u。 %擴充為矩陣T=t39。*ones(1,Na)。for i=1:1:Ntar rn=Ptar(i,1)。xn=Ptar(i,2)。sigma=Ptar(i,3)。 rtn=rn+RtcRrc。 RT=sqrt(rtn^2+(rtn*tan(thetaT)+xnv*U).^2)。 RR=sqrt(rn^2+(rn*tan(thetaT)+xnv*U).^2)。 R=RT+RR。 DT=TR/c。 phase=pi*Kr*DT.^22*pi/lambda*R。 s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)Tr/2).*(abs(v*Uxn)Lsar/2)。end。%方位向ffts_kt=fftshift(fft(fftshift(s_ut).39。)).39。%CS變換kc=4*pi/lambda。kc=kc*ones(1,Na)。kx=fu/v。p_kx0=sqrt(kc.^2kx.^2)。 %相位項泰勒展開的系數(shù)函數(shù)p_kx1=2*kc/c/p_kx0。p_kx2=2.*kx.^2/c^2./p_kx0.^3。C_kx=(c*p_kx1/2+1)。Ks_r=12*Kr*Rc.*p_kx2。Ks_kx_r=Kr/pi./Ks_r。r0=Rc。s2_ut=exp(j*pi*C_kx.*ones(Nr,1)*Ks_kx_r.*(t39。*ones(1,Na)2*r0*(1+C_kx)/c).^2)。 %設(shè)計的線性調(diào)頻信號S_cs=s_kt.*s2_ut。%距離向fftS_kw=fftshift(fft(fftshift(S_cs)))。%距離向匹配濾波w=2*pi*f。rmc_r=exp(j.*w*2*C_kx*r0/c).*exp(j.*w.^2/4/pi/Kr/(1+C_kx))。rmc_r=rmc_r39。*ones(1,Na)。S_rmc=S_kw.*rmc_r。%距離向ifftS_kt=fftshift(ifft(fftshift(S_rmc)))。d_kxr=4*pi/c^2*Kr*C_kx*(1+C_kx).*(Rcr0).^2。 %CS變換帶來的相位誤差S_kt=S_kt.*exp(j*d_kxr)。 %消除相位誤差%方位向匹配濾波FU=ones(Nr,1)*fu。H_kx=exp(j*pi/fdr*(FUfdc).^2)。 %方位向壓縮因子I_ut=S_kt.*H_kx。I_ut=fftshift(ifft(fftshift(I_ut.39。))).39。subplot(221)G=20*log10(abs(s_ut)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)方位角39。) ylabel(39。成像范圍39。)title(39。(圖 一)原始信號39。)subplot(222)G=20*log10(abs(S_rmc)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)方位角39。) ylabel(39。成像范圍39。)title(39。(圖 二) 距離向匹配濾波后頻譜39。)subplot(223)G=20*log10(abs(S_kt)+1e6)。gm=max(max(G))。gn=gm40。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)方位角39。) ylabel(39。成像范圍39。)title(39。(圖 三) 消除相位誤差后頻譜39。)subplot(224)G=20*log10(abs(I_ut)+1e6)。gm=max(max(G))。gn=gm60。 %顯示動態(tài)范圍40dBG=255/(gmgn)*(Ggn).*(Ggn)。imagesc(x,rRc,G),colormap(gray)grid on,axis tight,xlabel(39。成像目標(biāo)方位角39。) ylabel(39。成像范圍39。) title(39。(圖 四) 目標(biāo)圖象39。)29 / 30