Лабораторно-практическая работа № 7 (Вариант 17) - Расчет характеристик отрытых и замкнутых СМО, ручной счет и расчет с помощью MathLab
.docxФедеральное государственное бюджетное образовательное учреждение
высшего профессионального образования
«Петербургский государственный университет путей сообщения
Императора Александра I»
Кафедра «Высшая математика»
Лабораторно-практическая работа № 7
Расчет характеристик отрытых и замкнутых СМО
Ручной счет и расчет с помощью MathLab
Вариант 17
Выполнил: студент
факультета ТЭС
группы ПТМ-014
Алексеев А.Н.
Проверил: профессор
Боровских Ю.В.
Санкт-Петербург
2023
Открытая СМО типа MIMImIn
Дано:
Вычисление характеристик систем с помощью Matlab:
lambda=2; mu=4; m=3; n=3; P0=0.601
P=lambda/mu;
P_m=P/m;
N=m+n;
x=1:N;
x1=ones(1,N)*m;
xx=min(x,x1);
slave1=ones(1,N)*P;
slave2=slave1./xx;
slave=cumprod(slave2);
disp('Вероятность простоя системы');
P0=1/(1+sum(slave));
disp(P0);
for k=1:(m+n)
if k<=m
Pk(k)=P0*P^k/factorial(k);
else
Pk(k)=P0*P^k/(factorial(m)*m^(k-m));
end
end
disp(P0);
disp('Вероятности состояний 1,...,N');
Pk=slave*P0;
disp(Pk);
disp('Вероятность потери требования');
Potk=Pk(N);
disp(Potk);
disp('Вероятность ожидания перед обслуживанием');
Pw=sum(Pk(m:N-1));
disp(Pw);
disp('Вероятность обслуживания без ожидания');
Pww=P0+sum(Pk(1:(m-1)));
disp(Pww);
disp('Среднее число заявок в системе');
Qsr=(1:N)*Pk';
disp(Qsr);
xx1=x-x1;
disp('Среднее число заявок в очереди');
qsr=xx1((m+1):N)*(Pk((m+1):N))';
disp(qsr);
disp('Среднее число занятых приборов');
qs=xx*Pk';
disp(qs);
disp('Среднее время ожидания');
W=m*Pk(m)/2/mu*n*(n+1);
disp(W);
disp('Среднее время пребывания в системе');
T=W+1/mu;
disp(T);
t=0:0.01:10;
slv1=[];
slv2=[];
if n>1
for k=1:n-1
slv1=[slv1;(m*mu*P_m*t).^k./prod(1:k)];
slv2=[slv2;(m*mu*t).^k./prod(1:k)];
end
slv1=sum(slv1);
slv2=sum(slv2);
slv=slv1-slv2*P_m^n+1-P_m^n;
else
slv=1-P_m^n;
end
FW=1-Pk(m)/(1-P_m)*exp(-m*mu*t).*slv;
plot(t,FW);
Вероятность простоя системы
0.6061
Вероятности состояний 1,...,N
0.3030 0.0758 0.0126 0.0021 0.0004 0.0001
Вероятность потери требования
5.8456e-05
Вероятность ожидания перед обслуживанием
0.0151
Вероятность обслуживания без ожидания
0.6849
Среднее число заявок в системе
0.5030
Среднее число заявок в очереди
0.030
Среднее число занятых приборов
0.5000
Среднее время ожидания
0.7
Среднее время пребывания в системе
1.2
|
Вручную |
Matlab |
Вероятность простоя |
0.601 |
0.6050 |
Вероятность потери требования |
0.0000579 |
0.0000584 |
Вероятность ожидания перед обслуживанием |
0.0189 |
0.0151 |
Вероятность обслуживания без ожидания |
0.601 |
0.6849 |
Среднее число заявок в системе |
0.61 |
0.5030 |
Среднее число заявок в очереди |
0,0357 |
0.030 |
Среднее число занятых приборов |
0,494 |
0.5000 |
Среднее время ожидания |
|
0.7 |
Среднее время пребывания в системе |
|
1.2 |
Закрытая СМО типа MǀMǀmǀnǀS
Дано:
λ=7
μ=2
m=4
n=0
S=6
Вычисление характеристик систем вручную:
M=min{m+n, S}=4
Вычисление характеристик с помощью программы Matlab:
lambda=7; mu=2; m=4; n=0; S=6;
P=lambda/mu;
M=min(m+n,S);
x=1:M;
x1=ones(1,M-m)*m;
x2=S-x+1;
xx=[x(1:m),x1];
xx1=xx./x2;
slave1=ones(1,M)*P;
slave2=slave1./xx1;
slave=cumprod(slave2);
disp('Вероятность простоя');
P0=1/(1+sum(slave));
disp(P0);
disp('Вероятности состояний 1,...,M');
Pk=slave*P0;
disp(Pk);
disp('Вероятность потери требования');
PL=Pk(M);
disp(PL);
disp('Вероятность ожидания перед обслуживанием');
Pw=sum(Pk(m:M-1));
disp(Pw);
disp('Вероятность обслуживания без ожидания');
Pww=P0+sum(Pk(1:m-1));
disp(Pww);
disp('Среднее число заявок в системе');
Qsr=x*Pk';
disp(Qsr);
xq=x-xx;
disp('Среднее число заявок в очереди');
qsr=xq((m+1):M)*(Pk((m+1):M))';
disp(qsr);
disp('Среднее число занятых приборов');
qs=xx*Pk';
disp(qs);
disp('Среднее время ожидания');
W=(1:n)/m/mu.*(Pk((m+1):M));
disp(W(n));
disp('Среднее время пребывания в системе');
T=W(n)+1/mu;
disp(T);
t=0:0.01:10;
FW=1-Pw*exp(-(mu*t));
plot(t,FW);
Вероятность простоя
3.0173e-04
Вероятности состояний 1,...,M
0.0063 0.0554 0.2587 0.6792
Вероятность потери требования
0.6792
Вероятность ожидания перед обслуживанием
0
Вероятность обслуживания без ожидания
0.3208
Среднее число заявок в системе
3.6102
Среднее число заявок в очереди
0
Среднее число занятых приборов
3.6102
Среднее время ожидания
0.036
Среднее время пребывания в системе
0.536
|
Вручную |
Matlab |
Вероятность простоя |
0.0003 |
0.0003017 |
Вероятность потери требования |
0.6681 |
0.6792 |
Вероятность ожидания перед обслуживанием |
0 |
0 |
Вероятность обслуживания без ожидания |
0.3200 |
0.3208 |
Среднее число заявок в системе |
3.6102 |
3.6102 |
Среднее число заявок в очереди |
0.00001 |
0 |
Среднее число занятых приборов |
3.6 |
3.6102 |
Среднее время ожидания |
0.036 |
0.036 |
Среднее время пребывания в системе |
0,536 |
0,536 |