雷达目标恒虚警处理

主要是雷达目标检测……是单元平均恒虚警检测

clear all;
close all;
clc;
%%%%%%%%%%%产生瑞利噪声%%%%%%%%%%%%%%%%%%%%%%%
sigma=2;
t=1e-3;
fs=1e6;
ts=1/fs;
t1=0.05e-3:1/fs:0.2e-3-1/fs;
n=length（t1）;
rand（‘state‘0）;
u=rand（1n）;
rayleigh_noise=sqrt（2*log2（1./u））*sigma;

%%%%%%%%%%%%产生目标回波%%%%%%%%%%%%%%%%%%%
N=t/ts;
s_pc_1=[zeros（1100）1zeros（1N-101）];
noise=rand（1N）;
rayleigh_clutter=[zeros（150）rayleigh_noisezeros（1N-200）];
s_pc=s_pc_1+0.1*rayleigh_clutter+0.1*noise;
figureplot（（0:ts:t-ts）s_pc）
xlabel（‘t（单位:s）‘）title（‘叠加了瑞利分布杂波、热噪声的目标回波‘）;

%%%%%%%%%%%%%慢门限恒虚警处理%%%%%%%%%%%%%%%%%%%
 cfar_result=zeros（1N）;
 cfar_result（11）=s_pc（11）;
 for i=2:N
     cfar_result（i）=s_pc（1i）/mean（s_pc（11:i））;
 end
 figureplot（（0:ts:t-ts）cfar_result）
xlabel（‘t（单位:s）‘）title（‘采用慢门限处理结果‘）;

%%%%%%%%%%%%%%%快门限恒虚警处理%%%%%%%%%%%%%
cfar_k_result=zeros（1N）;

cfar_k_result（11）=s_pc（11）/（sqrt（2）/pi*mean（s_pc（12:17）））;
 
%第l点恒虚警处理时噪声均值由其后面的l6点的噪声决定
for i=2:16
%第2点到第16点的恒虚警处理的噪声均但由其前面和后面l6点的噪声共同决定
    noise_mean=sqrt（2）/pi*（mean（s_pc（11:i-1））+mean（s_pc（1i+1:i+16）））/2;
    cfar_k_result（1i）=s_pc（1i）/noise_mean;
end
for i=17:N-17
%正常的数据点的恒虚警处理的噪声均值由其前面和后面各26点的噪声中的决定
    noise_mean=sqrt（2）/pi*max（mean（s_pc（1i-16:i-1））mean（s_pc（1i+1:i+16）））;
    cfar_k_result（1i）=s_pc（1i）/noise_mean;
end
for i=N-16:N-1
%例数第l6点到倒数第2点恒虚警处理的噪声均值由其前面16和后面的噪声共同决定
    noise_mean=sqrt（2）/pi*（mean（s_pc（1i-16:i-1））+mean（s_pc（1i+1:N）））/2;
    cfar_k_result（1i）=s_pc（1i）/noise_mean;
end
%最后-点的桓虚警处理的噪声均值由其前面l6点的噪声决定
cfar_k_result（1N）=s_pc（1N）/（sqrt（2）/pi*mean（s_pc（1N-16:N-1）））;
figureplot（0:ts:t-tscfar_k_result）xlabel（‘t（单位:s）‘）title（‘采用快门限恒虚警处理处理结果‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1879  2007-10-22 20:52  ex9_14.m

----------- ---------  ---------- -----  ----

                 2097                    2