Практика / Practic14

clc;
clear all;
close all;

tmax=1000;
lam=3*8;
mu=4;

p=lam/mu;
p0=(1+p+(p^2)/2+(p^3)/6)^(-1);
Potk=(p^3)*p0/6;
Q=1-Potk;
A=lam*Q;

i=0;
tz(1)=0;
while cumsum(tz)<tmax
i=i+1;
tz(i)=poissrnd(mu,1,1);
end
tz=cumsum(tz);
tz=unique(tz);

in1=0;
in2=0;
in3=0;
tr1=0;
tr2=0;
tr3=0;
i=1;
j=1;
k=1;
tot(1)=0;

for t=1:1:tmax
if t==tz(i)
if in1==0
in1=1;
tr1=tz(i)+poissrnd(lam,1,1);
i=i+1;
elseif in2==0
in2=1;
tr2=tz(i)+poissrnd(lam,1,1);
i=i+1;
elseif in3==0
in3=1;
tr3=tz(i)+poissrnd(lam,1,1);
i=i+1;
else
tot(j)=tz(i);
i=i+1;
j=j+1;
end
end
if t==tr1
tk(k)=tr1;
k=k+1;
in1=0;
end
if t==tr2
tk(k)=tr2;
k=k+1;
in2=0;
end
if t==tr3
tk(k)=tr3;
k=k+1;
in3=0;
end
end

ytz=zeros(length(tz),1);
ytot=zeros(length(tot),1);
ytk=ones(length(tk),1);
figure(1)
stem(tz,ytz)
hold on
stem(tot,ytot,'x')
hold on
stem(tk,ytk)
legend('Время поступления заявок','Время поступления необслуженных заявок','Время окончания обработки заявок')
xlabel('Время')
axis([0 tmax 0 1])

Potk1=length(tot)/length(tz);
Q1=1-Potk1;
A1=lam*Q1;
xlim([0 50]);