CA CFAR_radar 对单元平均CFAR门限的求解进行仿真

对单元平均CFAR门限的求解进行仿真，还分别求解了GOCA－CFAR以及SOCA－CFAR门限

clear all
sigma1=10^0.75;
sigma2=10^1.5;
sigma=[sigma1sigma2];
Lambda=sigma.^2;
N=200;
Pf=10^（-3）;
r_clutter=[exprnd（Lambda（1）1100）exprnd（Lambda（2）1N-100）];
r_target=[zeros（140）10^4zeros（146-41-1）10^4.5zeros（195-46-1）10^3.5...
    zeros（1150-95-1）10^7zeros（1200-150）];
r_signal=r_clutter+r_target;
hold on;
plot（1:N10.*log10（r_signal））;
xlabel（‘距离单元‘）;
ylabel（‘强度（dB）‘）;
grid;
T=10.*log10（-sigma.^2.*log（Pf））;
plot（1:100T（1）‘r-‘101:200T（2）‘r-‘）;
title（‘实际情况下CFAR检测门限‘）;
for i=19:182;
    S（1）=r_signal（1）;
    for u=2:16;
        S（u）=S（u-1）+r_signal（u）;
    end
    for v=17:200;
        S（v）=S（v-1）-r_signal（v-16）+r_signal（v）;% S（v）为v前面（算上第v个）16个单元的和。
    end
    T_c（i）=32*（Pf^（-1/32）-1）/32*（S（i-3）+S（i+18））;
end
plot（1:18210*log10（T_c）‘g‘）;
%------画左右两端的门限曲线-------
%for k=3:18;            % 定义SS（k）为位置k（包括k）之前的16个单元的和
%   SS（1）=0;
%  SS（2）=r_signal（1）;
%    SS（k）=SS（k-1）+r_signal（k-1）;
%    T_c1（k）=Pf^（-1/（16+k-3））*（SS（k-2）+S（k+18））
%end
%T_c1（1）=Pf^（-1/16）*S（19）;
%T_c1（2）=Pf^（-1/16）*S（20）;
%for h=183:198;         % 定义SSS（h）为位置h+3（包括h+3）之后16个单元的和
%    SSS（182）=S（200）;   % S（200）=r_signal（185）+...+r_signal（200）
%    SSS（h）=SSS（h-1）-r_signal（h+2）;
%    T_c2（h）=Pf^（-1/（16+200-h-2））*（SSS（h）+S（h-3））;
%end
%T_c2（199）=Pf^（-1/16）*S（196）;
%T_c2（200）=Pf^（-1/16）*S（197）;
%plot（1:1810*log10（T_c1）‘g‘1:20010*log10（T_c2）‘g‘）;
%-------------GOCA-CFAR----------------
for i=19:182;
    S（1）=r_signal（1）;
    for u=2:16;
        S（u）=S（u-1）+r_signal（u）;
    end
    for v=17:200;
        S（v）=S（v-1）-r_signal（v-16）+r_signal（v）;% S（v）为v前面（算上第v个）16个单元的和。
    end
    T_d1（i）=（Pf^（-1/16）-1）*（max（S（i-3）S（i+18）））;
end
plot（1:18210*log10（T_d1）‘y‘）;

%-------------SOCA-CFAR----------------
for i=19:182;
    S（1）=r_signal（1）;
    for u=2:16;
        S（u）=S（u-1）+r_signal（u）;
    end
    for v=17:200;
        S（v）=S（v-1）-r_signal（v-16）+r_signal（v）;% S（v）为v前面（算上第v个）16个单元的和。
    end
    T_d2（i）=（Pf^（-1/16）-1）*（min（S（i-3）S（i+18）））;
end
plot（1:18210*log10（T_d2）‘k‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2177  2014-05-20 00:20  radar2.m