Page 23 - 4861
P. 23

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
   18   19   20   21   22   23   24   25   26   27   28