雷达系统仿真设计报告二 - 图文

雷达系统仿真设计报告

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 页