合成孔徑雷達(dá)點目標(biāo)仿真MATLAB程序.docx

《合成孔徑雷達(dá)點目標(biāo)仿真MATLAB程序.docx》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、合成孔徑雷達(dá)成像系統(tǒng)點目標(biāo)仿真源程序：clccloseallC=3e8;%光速Fc=1e9;%載波頻率lambda=C/Fc;%波長%成像區(qū)域Xmin=0;Xmax=50;Yc=10000;Y0=500;%SAR基本參數(shù)V=100;%雷達(dá)平臺速度H=0;%雷達(dá)平臺高度R0=sqrt(Yc^2+H^2);D=4;%天線孔徑長度Lsar=lambda*R0/D;%合成孔徑長度Tsar=Lsar/V;%合成孔徑時間Ka=-2*V^2/lambda/R0;%線性調(diào)頻率Ba=abs(Ka*Tsar);PR

2、F=2*Ba;%脈沖重復(fù)頻率PRT=1/PRF;ds=PRT;%脈沖重復(fù)周期Nslow=ceil((Xmax-Xmin+Lsar)/V/ds);%脈沖數(shù)Nslow=2^nextpow2(Nslow);%量化為2的指數(shù)sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%創(chuàng)建時間向量PRT=(Xmax-Xmin+Lsar)/V/Nslow;%更新PRF=1/PRT;%更新脈沖重復(fù)頻率fa=linspace(-0.5*PRF,0.5*PRF,Nslo

3、w);Tr=5e-6;%脈沖寬度Br=30e6;%調(diào)頻信號帶寬Kr=Br/Tr;%調(diào)頻率Fsr=2*Br;%快時間域取樣頻率dt=1/Fsr;%快時間域取樣間隔Rmin=sqrt((Yc-Y0)^2+H^2);Rmax=sqrt((Yc+Y0)^2+H^2+(Lsar/2)^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*Rm

4、ax/C+Tr-2*Rmin/C)/Nfast;%更新Fsr=1/dt;fr=linspace(-0.5*Fsr,0.5*Fsr,Nfast);DY=C/2/Br;%距離分辨率DX=D/2;%方位分辨率Ntarget=3;%目標(biāo)數(shù)目Ptarget=[Xmin,Yc,1%目標(biāo)位置Xmin,Yc+10*DY,1Xmin+20*DX,Yc+50*DY,1];K=Ntarget;%目標(biāo)數(shù)目N=Nslow;%慢時間采樣數(shù)M=Nfast;%快時間采樣數(shù)T=Ptarget;%目標(biāo)位置%合成孔徑回波仿真Srnm

5、=zeros(N,M);fork=1:1:Ksigma=T(k,3);Dslow=sn*V-T(k,1);R=sqrt(Dslow.^2+T(k,2)^2+H^2);tau=2*R/C;Dfast=ones(N,1)*tm-tau'*ones(1,M);phase=pi*Kr*Dfast.^2-(4*pi/lambda)*(R'*ones(1,M));Srnm=Srnm+sigma*exp(j*phase).*(0

6、nes(1,M));end%距離壓縮tr=tm-2*Rmin/C;Refr=exp(j*pi*Kr*tr.^2).*(0

7、bs(ta)

8、c(angle(Srnm))%灰度圖figure(2)colormap(gray)imagesc(row,col,255-abs(Sr));title('距離壓縮'),xlabel('距離向'),ylabel('方位向');figure(3)colormap(gray)imagesc(row,col,255-abs(Sa));title('方位壓縮'),xlabel('距離向'),ylabel('方位向');%輪廓圖figure(4)Ga=abs(Sa);a=max(max(Ga));conto