Практика / Practic12

clc;
clear all;
close all;

N=10;
lam=8;
mu=5;

tot=0;
tlam=poissrnd(lam,N,1);
tlam=cumsum(tlam);
a=tlam(1);
tmu=poissrnd(mu,1,1);
tv(1)=a+tmu;
for i=2:N
if tlam(i)<=a+tmu
tot(i)=tlam(i);
else
tv(i)=a+tmu;
a=tlam(i);
tmu=poissrnd(mu,1,1);
end
end
tot(find(tot==0))=[];
tv(find(tv==0))=[];

ytlam=zeros(N,1);
ytot=zeros(length(tot),1);
ytv=ones(length(tv),1);
figure(1)
stem(tlam,ytlam)
hold on
stem(tot,ytot,'x')
hold on
stem(tv,ytv)
legend('Время поступления заявок','Время поступления необслуженных заявок','Время окончания обработки заявок')
xlabel('Время')