Page 23 - 4861

S1=0;
for j=1:N+1
if j==1
F=1;
else
F=F*(j-1);
end
S1=S1+alfa^(j-1)/F;
end
S2=0;
for m=1:n
pr=1;
for k=1:m
pr=pr*(N+k*hama);
end
S2=S2+alfa^m/pr;
end
PS=S1+(alfa^N/F)*S2;
P0=1/PS;
%Обчислення Pj для значень j від 1 до N
fk=1;
for j=1:N
fk=fk*j;
P(j)=alfa^j*P0/fk;
end
%Обчислення Pj для значень j від N+1 до N+n
for j=N+1:N+n
fs=1;
for k=1:j-N
fs=fs*(N+k*hama);
end
P(j)=(alfa^j*P0)/(F*fs);
end
%Перевірка точності обчислень ергодичних значень Pj
[M,m]=size(X);
for i=1:m
if i==1
D(1)=(P0-X(M,1))*100/P0;
Z(1)=X(M,1);H(1)=P0;
else
D(i)=(P(i-1)-X(M,i))*100/P(i-1);
Z(i)=X(M,i);H(i)=P(i-1);
end
end
%Таблиця точності обчислень
Y=[Z' H' D'];

Лістинг 8.2 – Файл - функція формування матриці із правих частин рівняння (8.2)

function F=rksol(t,x,A)
%Формування матриці правих частин ДР
F=A*x;

Результатом розв'язку системи диференціальних рівнянь (8.2) чисельним методом є графіки
ймовірностей P , j , 0 10 (рис. 8.1) заходження системи в одному із станів j .
j
Аналіз графіків дає можливість зробити висновок про зміну ймовірностей P в часі. Видно,
j

22