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