new_cfar用于sar检测

一种新的cfar检测算法，适用于sar图像，测试结果还不错。

d=0;fa=0;%d为发现信号的个数，fa为虚警个数
N=100;
faPos=ones（1200）;
for h=1:N
s1=2;
s2=4; 
SNR = 30;
noiseA1 = exprnd（s1179）;
noiseA2 = exprnd（s1180）;
clutter = exprnd（s2141）;
A=2/（s1^2）*10^（SNR/10）;   %目标幅值，由SNR确定
SigPos= 82;             %目标位置
sigA =zeros（1200）;
sigA（SigPos） = A;
plus = [noiseA1 clutter noiseA2];
indata = sigA + plus    ;%产生了目标回波
n=8;
Pf=1/N;
T=Pf^（（-1）/（2*n））-1;
thresholdsave=zeros（1200）;
for i=1:200
    if 8        cwin=[i-8:i-1 i+1:i+8];
    elseif i==1
        cwin=[i+1:i+16];
    elseif i>1&&i<=8
        cwin=[1:i-1 i+1:17];
    elseif i>=193&&i<200
        cwin=[184:i-1 i+1:200];
    elseif i==200
        cwin=[i-16:i-1];
    end
    Zsort=sort（indata（1cwin））;
    k=9;
    Z=16*Zsort（1k）;
    threshold=Z*T;
    thresholdsave（1i）=threshold;
         if indata（1i）>threshold&&i==SigPos
             d=d+1;
         elseif indata（1i）>threshold&&i~=SigPos
             fa=fa+1;
             faPos（1i）=0;
         end
end
end
num = 200*N;
Pd = d/N
Pfa= fa/num
plot（faPos（1:）‘r.‘）

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

     文件       1102  2011-06-22 16:13  new_CFAR\OS_cal_clutter.m

     文件       1232  2011-06-22 14:57  new_CFAR\OS_cal_clutter_draw1.m

     文件       1297  2011-06-21 19:47  new_CFAR\OS_cal_clutter_draw3.m

     文件        945  2011-06-22 15:10  new_CFAR\OS_cal_onlynoise.m

     文件       1105  2011-06-17 20:35  new_CFAR\OS_cal_onlynoise_draw3.m

     目录          0  2017-05-27 09:55  new_CFAR

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

                 5681                    6