example1.asv

rang doppler imaging and motion compensation中的源代码

%
% EXAMPLE1.M
%
% This  is an example program to generate a simulated point
% imagery with translational and rotational motion. Given
% generated signature, The target motion is estimated by iterative
% method of motion compensation using weighted least squares method.
% The motion estimates is found via optimization for fast global
% minimum search.
%
% 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;     % bandw: Radar bandwidth.
   carf=3.0e+9;         % carf: Carrier frequency.     
   prf=5.40e+3;         % prf: Pulse repetion frequency.
   nn=256;              % nn: Number of pulses in a burst.
   mm=64;               % mm: Number of bursts.
   cc=299792500;
   Pfa=1e-3;            % Pfa: probability of false alarm
   tint=1/prf;          % tint: 脉冲时间间隔
   amplen=17;
   thre=0.6;
   sizemat=9;
   delf=bandw/(nn-1);   % delf: 频率步长
   delran=cc/(2*bandw); %距离分辨率
   radarp=[bandw;carf;prf;mm;nn;cc]; %% 雷达参数

   position=[-10.5 0;-9.5 0];  % position[x position, y position]: Position of point in the target.
   motionp=[-20;5;0;1e-35;1e-35]; % motionp[velocity, acceleration, jerk, angular velocity, angular acceleration]:
   sposition=[0;0]; %%散射点坐标
   numit=91;
   iter=100;

   bandfreq=tranfreqm(mm,nn,carf,bandw);    %% 计算每个脉冲的频率,,生成mm*nn距阵
   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');
   [d7q,indexc]=cleaning3(d7data,origin,mm,nn,Pfa);
   SNR=10^(db/10)*ones(1,64);
   RSNR=SNR./indexc;
   ASNR=mean(RSNR);
   md7q=sum(sum(abs(d7q)))/mmnn;
   [smt, emt]=CIAproc(d7q,amplen,thre,sizemat,RSNR);
   
   bd7q=[];
   for ita = 1:mm
   for itb=1:sizemat
   cmt=smt(ita,itb);
   dmt=emt(ita,itb);
   if  cmt ~= 0
   ad7q=sum(abs(d7q(ita,cmt:dmt)));
   bd7q=[bd7q;ad7q];
   end
   end
   end
   
   mmt=(emt-smt)';
   rmt=nonzeros(mmt);

   [cv,ca]=cramer2(radarp,ASNR,'dimension',rmt,bd7q,md7q);
   vv=0;
   aa=0;
   cr=sqrt(cv^2+ca^2)

   mf=fmine5('none','minimize',d7q,tt,bc,vv,aa,smt,emt,tint,sizemat,bt,cr)
   end