合成孔径雷达点目标仿真MATLAB程序

1、合成孔径雷达点目标仿真matlab程序合成孔径雷达成像系统点目标仿真源程序：clcclose allc=3e8; %光速fc=1e9; %载波频率lambda=c/fc; %波长%成像区域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); %创建时间向量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