Лаба 3 (Цепи Маркова и СМО) / L_3

clear
clc
close all

%% №1
P = [5 5 0 0; 0 5 5 0; 0 3 3 0; 0 0 3 3];

%% №2
MC = dtmc(P,"StateNames" , ["Healthy" "Unwell" "Sick" "Very sick"]);

%% №3
MCP = MC.P;
S=sum(MCP,2);

%% №4
figure('Name','Граф цепи маркова','NumberTitle','off')
graphplot(MC,'ColorEdges',true);

%% №5
P_cum = MCP;
for K = 1:4
for L=1:3
P_cum(K,L+1) = P_cum(K,L) + P_cum(K,L+1);
end
end
clear K L

%% №6/7/8
N=[200 1000 10000];
figure('Name','Моделирование поведения цепи Маркова','NumberTitle','off')
for i=1:3
subplot(3,1,i)
text=['Всего ',num2str(N(i)),' положений'];
title(text)
hold on
z = st(N(i),P_cum);
end

%% №10 (теоретический)

%% №12
lam=30;
mu=40;
n=6;
m=10;

p=lam/mu;
Psum=0;
for i=0:m
Psum=Psum+(p^i)/(factorial(i));
end
p0=Psum^(-1);
Potk=(p^(n+m))*p0/((n^m)*factorial(n));
Q=1-Potk;
A=lam*Q;

kzan=A/mu;
Loch=(p^(n+1))*(1-((p/n)^m)*(m+1-m*p/n))*p0/(n*factorial(n)*(1-p/n)^2);
Toch=Loch/lam;
Lsist=Loch+kzan;
Tsist=Lsist/lam;

% АЛИНА НАМ ОСТАЕТСЯ СДЕЛАТЬ 9 И 11-ЫЙ ПУНКТЫ
% но я совершенно не понимаю о чем речь, надо у Влада опять спрашивать,
% давай ты спросишь, мне кажется я уже успел и ЕА и его задолбать, но кто
% виноват, пусть методу напишут человеческим языком, я же не на иврите их
% прошу мне обьяснить задание