CFAR_MTI_rxh2.rar_CFAR检测matlab_MATLAB cfar_cfar matlab_mti cfar

function [coordinate_tar]=CFAR_MTI_rxh2(dif,threshold,eliminate,k,nrn_unit,nrn_protect,nan_unit,nan_protect)%k为多普勒中心的通道号 %eliminate为杂波所覆盖区域一般的多普勒通道数 %恒虚警处理%%%%%% 9,2,4,1 [nrn,nan]=size(dif);%%7407*130 if k-eliminate<=0; dif(:,k-eliminate+nan:nan)=0; dif(:,1:k+eliminate)=0; elseif k+eliminate>nan; dif(:,k-eliminate:nan)=0; dif(:,1:k+eliminate-nan)=0; else dif(:,k-eliminate:k+eliminate)=0; end f3=abs(dif).^2; f4=zeros(nrn,nan); if nrn<2*nrn_unit+1 error('图像行数不够');end if nan<2*nan_unit+1 error('图像列数不够');end for i=nrn_unit+1:nrn-nrn_unit for j=1:nan_protect a1=sum(sum(f3(i-nrn_unit:i+nrn_unit,nan-nan_unit+j:nan))); a2=sum(sum(f3(i-nrn_unit:i+nrn_unit,1:j+nan_unit))); c1=sum(sum(f3(i-nrn_protect:i+nrn_protect,nan-nan_protect+j:nan))); c2=sum(sum(f3(i-nrn_protect:i+nrn_protect,1:j+nan_protect))); f4(i,j)=a1+a2-c1-c2; end for j=nan_protect+1:nan_unit a1=sum(sum(f3(i-nrn_unit:i+nrn_unit,nan-nan_unit+j:nan))); a2=sum(sum(f3(i-nrn_unit:i+nrn_unit,1:j+nan_unit))); c1=sum(sum(f3(i-nrn_protect:i+nrn_protect,j-nan_protect:j+nan_protect))); f4(i,j)=a1+a2-c1; end for j=nan_unit+1:nan-nan_unit a1=sum(sum(f3(i-nrn_unit:i+nrn_unit,j-nan_unit:j+nan_unit))); c1=sum(sum(f3(i-nrn_protect:i+nrn_protect,j-nan_protect:j+nan_protect))); f4(i,j)=a1-c1; end for j=nan-nan_unit+1:nan-nan_protect a1=sum(sum(f3(i-nrn_unit:i+nrn_unit,j-nan_unit:nan))); a2=sum(sum(f3(i-nrn_unit:i+nrn_unit,1:j+nan_unit-nan))); c1=sum(sum(f3(i-nrn_protect:i+nrn_protect,j-nan_protect:j+nan_protect))); f4(i,j)=a1+a2-c1; end for j=nan-nan_protect+1:nan a1=sum(sum(f3(i-nrn_unit:i+nrn_unit,j-nan_unit:nan))); a2=sum(sum(f3(i-nrn_unit:i+nrn_unit,1:j+nan_unit-nan))); c1=sum(sum(f3(i-nrn_protect:i+nrn_protect,j-nan_protect:nan))); c2=sum(sum(f3(i-nrn_protect:i+nrn_protect,1:j+nan_protect-nan))); f4(i,j)=a1+a2-c1-c2; end end number=(nrn_unit*2+1)*(nan_unit*2+1)-(nrn_protect*2+1)*(nan_protect*2+1); f5=10^(threshold/10)*f4/number; f3(1:nrn_unit,:)=0; f3(nrn-nrn_unit+1:nrn,:)=0; [s1,s2]=find(f3>f5); coordinate_tar=[s1 s2];