%本程序研究在單個脈沖發現概率一定的條件下，被動導引頭對目標輻射源發現概率隨時間的變化。
%采用蒙特卡羅過程仿真方法；
%發現目標輻射源：連續發現5個目標輻射源脈沖；
%跟蹤丟失：連續丟失8個目標輻射源脈沖。
clc;clear;
a=0.2;
b=0.9;
ab=15;
Tpd=[];
for Pd=a:(b-a)/ab:b;%設置單個脈沖發現概率；
Tr=1;tt=[];
 for i=1:1000
 count=0;T=0;
    while(1)
        rn1=rand(1,1);
       if rn1>Pd
         count=count+1;
       else
         count=0;
       end
         T=T+Tr;
         if T>500*Tr
             break;
         end
       if count==5
          break;
       end
    end
    tt=[tt,T];
 end
 Tpd=[Tpd,ceil(mean(tt))];
end
Pd=a:(b-a)/ab:b
plot(Pd,Tpd*3e-3*50)
grid on
xlabel('單個脈沖發現概率')
% ylabel('重復周期數目')
% ylabel('發現輻射源的時間(s)')
ylabel('無人機所飛行的距離(m)')