雷达系统仿真设计报告
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 页