cacfar单元平均恒虚警检测方法

cacfar单元平均恒虚警检测方法，简单易懂

clear all;clc;close all;
%%%%  仿真数据  %%%%
T=0.25;M=100;tc=T/M;fc=1/tc;%%% 一个扫频周期T 为250ms ； M为每个周期采样点100 ； fc采样频率
N=256;n=1:N;%% 一个相关积累时间64
%%%
fs=1;%%  信号频率
c_db=-35;c=10^（c_db/20）;%% 噪声基底
s_db=-30;s=10^（s_db/20）;%% 信号幅度
signal=s.*exp（-j*2*pi*fs*T*n）;
nosic=c.*random（‘Normal‘011N）;
x=nosic+signal;
%%%
Fx=abs（fftshift（fft（x）））;
freq=（-1/（2*T）:1/（T*N）:（1/（2*T）-1/（T*N）））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CA-CFAR  %%%%%%%%%%%%%%%%%%%%%%%%%
M=N;P_fa=10.^（-6）;
R=16;%参考单元数
n=R/2;
L_slipper=R+1;%滑窗长度
L_move=1;%滑窗间隔
L_num=floor（（M-L_slipper）/L_move）+1;%滑窗次数
Z=zeros（1L_num）;
for i=1:L_num
    for j=1:L_slipper
        Z（i）=Z（i）+Fx（（i-1）*L_move+j）;
    end
    Z（i）=Z（i）/L_slipper;
end
T=P_fa.^（-1/R）-1
S=Z.*T;
figure;plot（freq20.*log10（Fx/max（Fx））‘k‘）;
hold on;plot（freq（n+1:M-n）20.*log10（S/max（S）））;
xlabel（‘频率‘）;ylabel（‘ 归一化幅值/dB‘）;title（‘CA-CFAR‘）;set（gcf‘color‘‘white‘）;grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         991  2014-05-04 09:44  CACFAR_sn.m