合成孔徑雷達點目標仿真MATLAB程序

1、合成孔徑雷達點目標仿真matlab程序合成孔徑雷達成像系統(tǒng)點目標仿真源程序：clcclose allc=3e8; %光速fc=1e9; %載波頻率lambda=c/fc; %波長%成像區(qū)域xmin=0;xmax=50;yc=10000;y0=500;%sar基本參數v=100; %雷達平臺速度h=0; %雷達平臺高度r0=sqrt(yc2+h2);d=4; %天線孔徑長度lsar=lambda*r0/d; %合成孔徑長度tsar=lsar/v; %合成孔徑時間ka=-2*v2/lambda/r0;%線性調頻率ba=abs(ka*tsar);prf=2*ba; %脈沖重復頻率prt=1/prf;

2、ds=prt; %脈沖重復周期 nslow=ceil(xmax-xmin+lsar)/v/ds);%脈沖數nslow=2nextpow2(nslow); %量化為2的指數sn=linspace(xmin-lsar/2)/v,(xmax+lsar/2)/v,nslow); %創(chuàng)建時間向量prt=(xmax-xmin+lsar)/v/nslow; %更新 prf=1/prt; % 更新脈沖重復頻率fa=linspace*prf,*prf,nslow);tr=5e-6; %脈沖寬度br=30e6; %調頻信號帶寬kr=br/tr; %調頻率fsr=2*br; %快時間域取樣頻率dt=1/fsr; %

3、快時間域取樣間隔rmin=sqrt(yc-y0)2+h2);rmax=sqrt(yc+y0)2+h2+(lsar/2)2);nfast=ceil(2*(rmax-rmin)/c/dt+tr/dt);nfast=2nextpow2(nfast); tm=linspace(2*rmin/c,2*rmax/c+tr,nfast);dt=(2*rmax/c+tr-2*rmin/c)/nfast; %更新fsr=1/dt;fr=linspace*fsr,*fsr,nfast);dy=c/2/br; %距離分辨率dx=d/2; %方位分辨率ntarget=3; %目標數目ptarget=xmin,yc,1

4、 %目標位置 xmin,yc+10*dy,1 xmin+20*dx,yc+50*dy,1; k=ntarget; %目標數目n=nslow; %慢時間采樣數m=nfast; %快時間采樣數t=ptarget; %目標位置%合成孔徑回波仿真srnm=zeros(n,m);for k=1:1:k sigma=t(k,3); dslow=sn*v-t(k,1); r=sqrt(dslow.2+t(k,2)2+h2); tau=2*r/c; dfast=ones(n,1)*tm-tau'*ones(1,m); phase=pi*kr*dfast.2-(4*pi/lambda)*(r'*

5、ones(1,m); srnm=srnm+sigma*exp(j*phase).*(0<dfast&dfast<tr).*(abs(dslow)<lsar/2)'*ones(1,m);end%距離壓縮tr=tm-2*rmin/c;refr=exp(j*pi*kr*tr.2).*(0<tr&tr<tr);%距離壓縮參考函數f_refr=fft(refr);sr=zeros(n,m);for k2=1:1:m temp1=fft(srnm(k2,:); fsrnm=temp1.*conj(f_refr); sr(k2,:)=ifft(fsrnm

6、);end %方位壓縮ta=sn-(xmin+xmax)/2/v;refa=exp(j*pi*ka*ta.2).*(abs(ta)<tsar/2); %方位壓縮參考函數f_refa=fft(refa);sa=zeros(n,m);for k3=1:1:m temp2=fft(sr(:,k3); f_sa=temp2.*conj(f_refa.'); sa(:,k3)=fftshift(ifft(f_sa);end row=tm*c/2;col=sn*v;%繪圖%回波雷達數據figure(1) subplot(211) imagesc(abs(srnm);title('sa

7、r data') subplot(212) imagesc(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);contour(row,col,ga,*a,a,'b');grid onx1=sqrt(h2+min(abs(ptarget(:,2).2)-5*dy;x2=sqrt(h2+max(abs(ptarget(:,2).2)+5*dy;y1=min(ptarget(:,1)-5*dx; y2=max(ptarget(:,1)+5*dx; axis(x1,x2,y1,y2)%三維成像圖 figure(5) mesh(ga(400:600),(200:500);axis tight xl