example2.m

rang doppler imaging and motion compensation中的源代码

%
% EXAMPLE2.M
%
% This  is an example program to generate a simulated point
% imagery with translational and rotational motion. The
% generated signature is then motion compensated and its 
% image is displayed in spatial domain.
%
% Radar parameter must be changed for each application.
% Translational motion parameter estimates must be provided
% at motioncomp function in a low vector form.
%
% bandw: Radar bandwidth.
% carf: Carrier frequency.
% prf: Pulse repetion frequency.
% nn: Number of pulses in a burst.
% mm: Number of bursts.
% position[x position, y position]: Position of point in the target.
% motionp[velocity, acceleration, jerk, angular velocity, angular acceleration]:
% provides translational and rotational motion of the target.
% dbb: Signal to noise ratio.
%

   clear all;
   close; 
   bandw=2.7650e+8;   
   carf=3.0e+9;
   prf=5.40e+3;
   nn=256;
   mm=64;
   cc=299792500;
   Pfa=1e-3;
   tint=1/prf;
   amplen=17;
   thre=0.6;
   sizemat=9;
   delf=bandw/(nn-1);
   delran=cc/(2*bandw);
   radarp=[bandw;carf;prf;mm;nn;cc];

   position=[-11.0 0;-9.0 0];
   motionp=[-20;5;0;0.02;1e-35];
   sposition=[0;0];
   numit=91;
   iter=100;

   bandfreq=tranfreqm(mm,nn,carf,bandw);
   ifour=4*pi/cc*(-i);
   bc=bandfreq.*ifour;
   tt=timegen(prf,mm,nn);
   bt=tt(:,1);
   mmnn=mm*nn; 
   dbb=0:0;

   posc=[];
   sda=[];
   sdb=[];
   raovv=[];
   for db = dbb
   acv=[];
   apos=[];


   Umnnoise=signalgen(radarp,position,motionp,sposition,db,'none');
   [origin,d7data]=loaddata(radarp,'simuldat',Umnnoise);
   adata=dispdft(d7data,delran,1,0,'none');

   bandfreq=tranfreqm(mm,nn,carf,bandw);
   ifour=4*pi/cc*(-i);
   bc=bandfreq.*ifour;
   tt=timegen(prf,mm,nn);
   bt=tt(:,1);
  
      motmat=motioncomp(origin,[0 0],tt,bc);
      mq=1*mm;
      nq=1*nn;
      de=[delran, 0.7871];
      resul1=rdisplay('graydisp',motmat,'imagery',[mq,nq],'rect',de,0);
      xlabel('SLANT RANGE (m)')
      ylabel('CROSS RANGE (m)')
  end