dispdft.m

rang doppler imaging and motion compensation中的源代码

   function [fdata,rangei] = dispdft(data,delran,bnum,sub,dy)
   %
   % Display or save the image of a burst in spatial domain.
   % delran = range resolution given cc/(2*bandw).
   % bnum = burst index number to be displayed.
   % sub = subplot number. if one plot is drawn, then sub = 0.
   % dy = 'disp', 'dispdeps' or 'none'
   %

   [m n] = size(data);
   leng=length(data);
   if m ~= leng
      temp=m;
      m = n;
      n = temp ;
      data = data.';
   end

  lda=length(data);
  pdc=(-lda/2:(lda/2-1))*delran;
  fdata=fft(data);

  bleng=length(bnum);
  if bleng ~= 1
  II=[];
  for it=1:n 
  bdata=abs(fftshift(fdata(:,it)));
  [Y,I]=max(bdata);
  II=[II,I];
  end
  rangei=II(bnum).*delran-lda/2*delran;
  else
  bdata=abs(fftshift(fdata(:,bnum)));
  cdata=normalize(bdata);
  end

  if strcmp(dy,'disp') == 1 | strcmp(dy,'dispdeps') == 1
  if sub ~= 0
  subplot(1,2,sub)
  plot(pdc,cdata)
  if sub == 1
  title('(a)')
  else
  title('(b)')
  end
  xlabel('Range (m)')
  ylabel('Amplitude')
  else
  plot(pdc,bdata)
  end
  if (sub == 2) & (strcmp(dy,'dispdeps')) == 1
  print fouriabs -deps
  end
  end