À×´ïϵͳ·ÂÕæÉè¼Æ±¨¸æ
xlabel('t£¨µ¥Î»£ºs£©'),title('¹Ì¶¨ÔÓ²¨¶þ´Î¶ÔÏûºó½á¹û'); clear s_MTI_result_1;
% %%%%%%%%%MTD´¦ÀíºÍÇóÄ£%%%%%%%%%%%%%%%%%%% s_MTD_result_1=zeros(N_echo_frame-2,L); s_MTD_result=zeros(1,L); for i=1:L
s_MTD_result_1(:,i)=fft(s_MTI_result(:,i),16); end for i=1:L
s_MTD_result(i)=max(abs(s_MTD_result_1(:,i))); end
subplot(234),plot(0:ts:T_frame-ts,s_MTD_result),
xlabel('t£¨µ¥Î»£ºs£©'),title('MTD´¦ÀíºóÇóÄ£½á¹û(ÐźÅ×î´óͨµÀ)'); clear s_MTD_result_1;
%?AR´¦Àí%%%%%%%%%%%%%%%%
cfar_result(1,1)=s_MTD_result(1,1)/(sqrt(2)/pi*mean(s_MTD_result(1,2:17)));%
µÚ1µãºãÐ龯´¦ÀíʱÔëÉù¾ùÖµÓÉÆäºóÃæµÄ16µãµÄÔëÉù¾ö¶¨
for i=2:16 %%%%µÚ2µãµ½µÚ16µãµÄºãÐ龯µÄÔëÉù¾ùÓÉÆäÇ°ÃæºÍºóÃæµÄ16µãµÄÔë
Éù¹²Í¬¾ö¶¨
noise_mean=sqrt(2)/pi*(mean(s_MTD_result(1,1:i-1))+mean(s_MTD_result(1,i+1:i+16)))/2;
cfar_result(1,i)=s_MTD_result(1,i)/noise_mean; end
µÚ 25 Ò³ ¹² 26 Ò³
À×´ïϵͳ·ÂÕæÉè¼Æ±¨¸æ
for i=17:L-17
%Õý³£µÄÊý¾ÝµãºãÐ龯´¦ÀíµÄÔëÉù¾ùÖµÓÉÆäÇ°ÃæºÍºóÃæ¸÷16µãÔëÉù¹²Í¬¾ö¶¨
noise_mean=sqrt(2)/pi*max(mean(s_MTD_result(1,i-16:i-1)),mean(s_MTD_result(1,i+1:i+16)));
cfar_result(1,i)=s_MTD_result(1,i)/noise_mean; end
for i=L-16:L-1
%µ¹ÊýµÚ16µãµ½µ¹ÊýµÚ2µãºãÐ龯´¦ÀíµÄÔëÉù¾ùÖµÓÉÆäÇ°Ãæ16µãºÍºóÃæµÄÔëÉù¹²Í¬
¾ö¶¨
noise_mean=sqrt(2)/pi*(mean(s_MTD_result(1,i-16:i-1))+mean(s_MTD_result(1,(i+1):L)))/2;
cfar_result(1,i)=s_MTD_result(1,i)/noise_mean; end
%%%×îºóÒ»µãµÄºãÐ龯´¦ÀíµÄÔëÉù¾ùÖµÓÉÆäÇ°ÃæµÄ16µãµÄÔëÉù¾ö¶¨
cfar_result(1,L)=s_MTD_result(1,L)/(sqrt(2)/pi*mean(s_MTD_result(1,L-16:L-1)));
s_result=zeros(1,L); for i=1:L
if s_MTD_result(1,i)>=cfar_result(1,i) s_result(1,i)=s_MTD_result(1,i); else s_result(1,i)=0; end end
subplot(235),plot((0:ts:T_frame-ts),cfar_result),xlabel('t(µ¥Î»):s'),title('ºãÐ龯µÄÃÅÏÞÖµ');
subplot(236),plot((0:ts:T_frame-ts),s_result),xlabel('t(µ¥Î»):s'),title('ºãÐ龯´¦Àí½á¹û');
µÚ 26 Ò³ ¹² 26 Ò³