d=0;fa=0;%d為發(fā)現(xiàn)信號的個數(shù)，fa為虛警個數(shù)
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????;%產(chǎn)生了目標回波
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