Практика / Practic13

clear
close all
clc

lamda=1/6;
mu=1/2;
Q=mu/(mu+lamda);
A=lamda*Q;
p=lamda/(mu+lamda);
ro=lamda/mu;

matrix=ones(12,1);
for i=1:12
sv=randi(12);
switch sv
case 1
matrix(randi(12),1)=0;
case 2
matrix(randi(12),1)=0;
case 3
matrix(randi(12),1)=0;
end
end

matrix1=matrix;
for i=1:12
if matrix1(i,1)>0
matrix1(i,1)=NaN;
end
if matrix(i,1) == 0
matrix(i,1)=NaN;
end
end

t1=zeros(12,1);
e1 = exprnd(1,1,12);
for i = 2:12
while e1(i) <= e1(i-1)
e1(i)=e1(i-1)+exprnd(1);
end
t1(i,1)=t1(i-1,1)+e1(1,i);
end
t1=t1';

k=zeros(12,1);
for i=1:12
if matrix1(i,1)==0
k(i,1)=i;
end
end

k1=zeros(12:1);
for i=1:12
if k(i,1) > 0
k1(i,1)=t1(1,i);
end
end

Np = normalize(k1,'norm',1);
for i=1:length(Np)
if Np(i,1)==0
Np(i,1)=NaN;
end
end
Npl=zeros(1,length(Np));

matrix3=ones(1,12);
figure('Name','Одноканальная СМО с отказами','NumberTitle','off')
title('Одна из возможных реализаций')
stem(t1,matrix3,'k')
hold on
stem(t1,matrix1,'r','Marker','*')
if Npl ~= NaN
stem(Npl,Np);
end
hold on

matrix4=ones(12,1);
for i = 1:12
if matrix1(i,1) ~= 0
matrix4(i,1)=0;
else
matrix4(i,1)= NaN;
end
end
stem(t1,matrix4);
hold on
legend('На входе','Отказали','Этап на котором появилась','Приняли')
ylabel('Выходящие')
xlabel('Входящие')
xticks(t1);
xticklabels({'1','2','3','4','5','6','7','8','9','10','11','12'})

clear e1 i k k1 Np Npl t1