radar.m

一些常被用于教学或者参考的概率论的实例的源代码

%radar.m
%
%
max=input('Number of total observations? ');
p=input('Minimum decimal fraction of total observations to declare a target? ');
ps=input('0/1 signal threshold probability? ');
ops=1-ps;
pn=input('0/1 noise threshold probability? ');
opn=1-pn;
clear pdv;
clear pfav;
for n=1:max
   m=floor(n*p)+1;
   PD=0;
   PFA=0;
   for k=m:n
      factor=binomial(n,k);
      PD=PD+factor*(ps^k)*ops^(n-k);
      PFA=PFA+factor*(pn^k)*opn^(n-k);
   end
   pdv(n)=PD;
   pfav(n)=PFA;
end
x=1:max;
plot(x,pdv,'-',x,pfav,':')
legend('Probability of Detection','Probability of False Alarm')
xlabel('number of observations')
ylabel('probability')
title('Probabilities of Detection and False Alarm for Radar Binary Integration')
figure(1)