probint3.m

此程序给出了不同情况的雷达距离测量仿真平台

% Probability of Intercept Simulation
% Using Three Windows
% ------------------------------------

clear;clc;clf;

% Input Window Parameters

s=[];
T1=6.13;      	% Emitter 1 Period in s	
T2=1.05;    	% Emitter 2 Period in s
T3=.043;       % Receiver Scan Period in s
t1=.0333;	   % Intercept 1 Pulse in s
t2=.00278;  	% Intercept 2 Pulse in s
t3=.0021;      % Receiver Open Time in s
df1=100*t1/T1;	% Duty Factor 1 in % 
df2=100*t2/T2;	% Duty Factor 2 in %
df3=100*t3/T3;	% Duty Factor 3 in %
n=200;         % Number of Trials

% Simulate Probability of Intercept

for i=1:n;
p=[];
c=0;
r1=rand;r3=rand;
while length(p)<1,
   t=linspace(c*T2,c*T2+t2,100);
   x=square(2*pi*(t/T1+r1),df1);
   xx=(x+1)/2;
   y=square(2*pi*(t/T3+r3),df3);
   yy=(y+1)/2;
   z=xx.*yy;
   p=find(z);
   tot=c*T2;
   if tot>5000,break,end;
   c=c+1;
end;

% Find Probability of Intercept

s(i)=c*T2;
phase(i)=r1;
end;
[hits,time]=hist(s,600);
cumulative=cumsum(hits);
yz=cumulative/n;

% Plot Probability of Intercept

x1=time;
stairs(x1,yz);
%axis([0,500,0,1]);
grid;
title('Probability of Main Beam Intercept');
xlabel('Intercept Duration-seconds');
ylabel('Cumulative Probability of Intercept');
hold on;

% Calculate Probability of Intercept

POI95=.95;
To=(T1*T2*T3)/(t1*t2+t2*t3+t1*t3);	% Mean Intercept Time
Tj=-To*log(1-POI95);                % Time POI=.95
tx=0:1:Tj;
POI=1-exp(-tx/To);
plot(tx,POI);
hold off;