cramer2.m

rang doppler imaging and motion compensation中的源代码

   function [delv,dela] = cramer2(radarp,snr,dim,L,bd,mumu)
   % 
   % Calculates Cramer-Rao lower bound of the velocity
   % and acceleration estimates. It provides the precision
   % achievable in the motion estimates so that it limites
   % the number of iteration in motion estimation.
   %
   
    bandw=radarp(1);
    carf=radarp(2);
    prf=radarp(3);
    mm=radarp(4);
    nn=radarp(5);
    cc=radarp(6);
    delran=cc/(2*bandw)
    lsnr=sqrt(2*snr)
    T=mm*nn/prf;  

    if strcmp(dim,'twodimension') == 1
    delv=sqrt(18*delran^2./(mm^3*nn*T^2*pi^2*lsnr));
    dela=sqrt(270*delran^2./(mm^3*nn*T^4*pi^2*lsnr));
    elseif strcmp(dim,'onedimension') == 1
    delv=sqrt(18*delran^2./(mm*nn*T^2*pi^2*lsnr));
    dela=sqrt(270*delran^2./(mm*nn*T^4*pi^2*lsnr));
    else
    sizel=length(L)-1;
    k1=(lsnr*pi*T)/(6*nn*delran);
    delv=sqrt(1/sum((k1*bd).^2));
    k2=(lsnr*pi*T^2)/(sqrt(2160)*nn*delran);
    dela=sqrt(1/sum((k2*bd).^2));
    end