Регулярні ланцюги Маркова та їх числові характеристики (реферат)

РЕФЕРАТ На тему:

Регулярні ланцюги Маркова та їх числові характеристики Як уже наголошувалося, регулярні ланцюги Маркова є частинним випадком ергодичних ланцюгів.

Регулярний ланцюг — це такий ланцюг, для якого за деякого значення k (кількість кроків) у матриці $${}^{\left(k\right)}$$ не буде нульових елементів. А це означає, що перехід системи з будь-якого стану до будь-якого іншого стану здійсниться за k кроків.

1. 1.Визначення стаціонарних (фінальних) імовірностей.

Нехай задано перехідну матрицю егулярного ланцюга Маркова. Тоді виконуються такі умови:

1) для будь-якого ймовірнісного вектора $$\stackrel{}{a}$$ вектор $$\stackrel{}{a}{}^{\left(k\right)}$$ прямує до вектора $$\stackrel{}{W}$$ за $$k->\mathrm{}-$$.

2) $$\stackrel{}{W}$$ — єдиний вектор, який має властивість.

$$\stackrel{}{W}=\stackrel{}{W}-$$ (39).

3) $$\stackrel{}{W}=\stackrel{}{W}=\stackrel{}{W}\text{.}$$.

Перша умова означає, що незалежно від імовірнісного вектора $$\stackrel{}{a}$$ (це може бути початковий вектор станів системи) для регулярних ланцюгів Маркова, які моделюються певною матрицею виконується рівність.

$$\underset{k->\mathrm{}}{\text{lim}}\stackrel{}{a}{}^{\left(k\right)}=\stackrel{}{W}\text{.}$$ (40).

Вектор $$\stackrel{}{W}$$ називають вектором стаціонарних (фінальних) імовірностей. Його компоненти задовольняють умову.

$$\underset{i=1}{\overset{N}{}}{w}_i=1\text{.}$$ (41).

Приклад 19. Дано однокрокову матрицю ймовірностей переходу, що є регулярним ланцюгом Маркова:

$$=\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,3}& \mathrm{0,2}\\ \mathrm{0,5}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,2}\\ \mathrm{0,}\text{15}& \mathrm{0,}\text{25}& \mathrm{0,3}& \mathrm{0,1}\\ \mathrm{0,}\text{05}& \mathrm{0,1}& \mathrm{0,}\text{15}& \mathrm{0,7}\end{array}\right)$$ ,.

а також вектор $$\stackrel{}{a}=\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,}\text{08}& \mathrm{0,}\text{32}& \mathrm{0,5}\end{array}\right)\text{.}$$.

Знайти стаціонарні (фінальні) ймовірності та визначити кількість кроків k, за яких регулярний процес вийде на стаціонарний режим.

Розв’язання. Скориставшись (40), дістанемо.

$$\stackrel{}{a}{}^{\left(k\right)}=\stackrel{}{W}->$$.

$$->\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,}\text{08}& \mathrm{0,}\text{32}& \mathrm{0,5}\end{array}\right){\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,3}& \mathrm{0,2}\\ \mathrm{0,5}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,2}\\ \mathrm{0,}\text{15}& \mathrm{0,}\text{25}& \mathrm{0,3}& \mathrm{0,3}\\ \mathrm{0,}\text{05}& \mathrm{0,1}& \mathrm{0,}\text{15}& \mathrm{0,7}\end{array}\right)}^k\text{.}$$.

За $$k=5$$

маємо:

.

$$\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,}\text{08}& \mathrm{0,}\text{32}& \mathrm{0,5}\end{array}\right){\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,3}& \mathrm{0,2}\\ \mathrm{0,5}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,2}\\ \mathrm{0,}\text{15}& \mathrm{0,}\text{25}& \mathrm{0,3}& \mathrm{0,3}\\ \mathrm{0,}\text{05}& \mathrm{0,1}& \mathrm{0,}\text{15}& \mathrm{0,7}\end{array}\right)}^5=\left(\begin{array}{cccc}\mathrm{0,}\text{208}& \mathrm{0,}\text{145}& \mathrm{0,}\text{204}& \mathrm{0,}\text{443}\end{array}\right)-$$.

за $$k=8$$

:

.

$$\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,}\text{08}& \mathrm{0,}\text{32}& \mathrm{0,5}\end{array}\right){\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,3}& \mathrm{0,2}\\ \mathrm{0,5}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,2}\\ \mathrm{0,}\text{15}& \mathrm{0,}\text{25}& \mathrm{0,3}& \mathrm{0,3}\\ \mathrm{0,}\text{05}& \mathrm{0,1}& \mathrm{0,}\text{15}& \mathrm{0,7}\end{array}\right)}^8=\left(\begin{array}{cccc}\mathrm{0,}\text{209}& \mathrm{0,}\text{145}& \mathrm{0,}\text{205}& \mathrm{0,}\text{441}\end{array}\right)-$$.

за $$k=\text{15}$$

:

.

$$\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,}\text{08}& \mathrm{0,}\text{32}& \mathrm{0,5}\end{array}\right){\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,3}& \mathrm{0,2}\\ \mathrm{0,5}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,2}\\ \mathrm{0,}\text{15}& \mathrm{0,}\text{25}& \mathrm{0,3}& \mathrm{0,3}\\ \mathrm{0,}\text{05}& \mathrm{0,1}& \mathrm{0,}\text{15}& \mathrm{0,7}\end{array}\right)}^{\text{15}}=\left(\begin{array}{cccc}\mathrm{0,}\text{209}& \mathrm{0,}\text{145}& \mathrm{0,}\text{205}& \mathrm{0,}\text{441}\end{array}\right)\text{.}$$.

За $$k=8$$ і $$k=\text{15}$$ імовірності не змінилися. Отже, можна стверджувати, що при $$k=8$$ регулярний марковський процес досягає стаціонарних (фінальних) імовірностей, а саме: якщо далі збільшувати кількість кроків k, імовірності станів марковського процесу не змінюватимуться.

Отже, згідно зі щойно сказаним для регулярних ланцюгів Маркова маємо: $$\underset{k->\mathrm{}}{\text{lim}}{}^{\left(k\right)}=W,$$.

де.

$$W=\left(\begin{array}{ccccc}{w}_1& w_2& w_3& \dots & w_N\\ w_1& w_2& w_3& \dots & w_N\\ \dots & \dots & \dots & \dots & \dots \\ w_1& w_2& w_3& \dots & w_N\end{array}\right)\text{.}$$ (42).

Тут W є матрицею стаціонарних (фінальних) імовірностей, в якій усі рядки однакові.

Практично вже за $$k>=\text{10}$$ система стає стаціонарною.

Приклад 20. За даною однокроковою матрицею ймовірностей переходу для регулярного марковського процесу.

$$=\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)$$.

знайти матрицю W.

Розв’язання. Обчислимо $${}^k$$ для послідовних значень $$k=4-8-\text{10}\text{.}$$.

$${}^4={\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)}^4=\left(\begin{array}{cccc}\mathrm{0,}\text{071}& \mathrm{0,}\text{136}& \mathrm{0,}\text{378}& \mathrm{0,}\text{416}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{137}& \mathrm{0,}\text{377}& \mathrm{0,}\text{416}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{069}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{418}\end{array}\right)$$ ;

$${}^8={\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)}^8=\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)$$ ;

$${}^{\text{10}}={\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)}^{\text{10}}=\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)\text{.}$$.

Як бачимо, за $$k=8$$ марковський процес стає стаціонарним. Стаціонарним вектором буде $$\stackrel{}{W}=\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)\text{.}$$.

Застосовуючи результати, здобуті в прикладі 20, підтвердимо властивості.

$$\stackrel{}{W}=\stackrel{}{W}->$$.

$$->\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,4}\end{array}\right)\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,3}& \mathrm{0,5}\end{array}\right)=$$.

$$=\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{372}& \mathrm{0,4}\end{array}\right)-$$.

$$\stackrel{}{W}=\stackrel{}{W}=\stackrel{}{W}->$$.

$$->\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)=$$.

$$=\left(\begin{array}{cccc}0& \mathrm{0,4}& \mathrm{0,}\text{25}& \mathrm{0,}\text{35}\\ \mathrm{0,}\text{15}& \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,3}\\ \mathrm{0,}\text{02}& \mathrm{0,}\text{18}& \mathrm{0,}\text{42}& \mathrm{0,}\text{38}\\ \mathrm{0,1}& \mathrm{0,}\text{03}& \mathrm{0,}\text{37}& \mathrm{0,5}\end{array}\right)\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)=$$.

$$=\left(\begin{array}{cccc}\mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\\ \mathrm{0,}\text{07}& \mathrm{0,}\text{135}& \mathrm{0,}\text{378}& \mathrm{0,}\text{417}\end{array}\right)\text{.}$$.

Із (39) і (42) дістаємо систему рівнянь.

$$\{\begin{array}{c}\stackrel{}{W}=\stackrel{}{W}\\ w_i=\mathrm{1,}\end{array}$$.

з якої знаходимо значення елементів стаціонарного вектора $$\stackrel{}{W}=\left(\begin{array}{cccc}{w}_1& w_2& w_3\dots & w_N\end{array}\right)\text{.}$$.

Існування стаціонарних імовірностей для регулярних ланцюгів Маркова в попередніх міркуваннях було взято до уваги без доведення, але досить наочно проілюстровано конкретними прикладами.

Тепер на вищому рівні строгості висвітлимо питання про існування стаціонарного вектора $$\stackrel{}{W}$$ для регулярного процесу.

Із цією метою розглянемо стохастичну матрицю.

$$=\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,2}& \mathrm{0,5}& \mathrm{0,2}\\ \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}& \mathrm{0,2}\\ \mathrm{0,2}& \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}\\ \mathrm{0,1}& \mathrm{0,1}& \mathrm{0,5}& \mathrm{0,3}\end{array}\right),$$.

мінімальний елемент якої $$=\mathrm{0,1}$$, а також деякий вектор, наприклад.

$$\stackrel{}{x}=\left(\begin{array}{c}\mathrm{0,}\text{25}\\ \mathrm{0,}\text{05}\\ \mathrm{0,4}\\ \mathrm{0,3}\end{array}\right)$$ ,.

найменший і найбільший компонент якого відповідно такі: $$m_0=\mathrm{0,}\text{05}$$ і $$M_0=\mathrm{0,4}$$ .

Принагідно зазначимо, що в загальному випадку таких компонентів може бути кілька.

Знайдемо вектор

$$\stackrel{}{x}=\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,2}& \mathrm{0,5}& \mathrm{0,}\text{22}\\ \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}& \mathrm{0,2}\\ \mathrm{0,2}& \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}\\ \mathrm{0,1}& \mathrm{0,1}& \mathrm{0,5}& \mathrm{0,3}\end{array}\right)\left(\begin{array}{c}\mathrm{0,}\text{25}\\ \mathrm{0,}\text{05}\\ \mathrm{0,4}\\ \mathrm{0,3}\end{array}\right)=\left(\begin{array}{c}\mathrm{0,}\text{301}\\ \mathrm{0,3}\\ \mathrm{0,}\text{225}\\ \mathrm{0,}\text{32}\end{array}\right)$$ ,.

де $$m_1=\mathrm{0,}\text{225}$$ є найменшим компонентом цього вектора, а $$M_1=\mathrm{0,}\text{32}$$ — найбільшим.

Як бачимо, $$M_1<M_0,$$ $$m_1>m_0\text{.}$$ У загальному випадку $$M_1<=M_0,$$ $$m_1>=m_0\text{.}$$.

Потрібно довести, що виконується нерівність.

$$M_1-m_1<=\left(1-2\right)\left(M_0-m_0\right)\text{.}$$.

Побудуємо вектор $${\stackrel{}{x}}^{},$$ скориставшись для цього вектором $$\stackrel{}{x},$$ в якому замінимо всі компоненти, крім $$m_0=\mathrm{0,}\text{05},$$ на $$M_0=\mathrm{0,4}$$ :

$${\stackrel{}{x}}^{}=\left(\begin{array}{c}\mathrm{0,4}\\ \mathrm{0,}\text{05}\\ \mathrm{0,4}\\ \mathrm{0,4}\end{array}\right)\text{.}$$.

Отже, $$\stackrel{}{x}<={\stackrel{}{x}}^{}\text{.}$$ Далі знаходимо.

$${\stackrel{}{x}}^{}=\left(\begin{array}{cccc}\mathrm{0,1}& \mathrm{0,2}& \mathrm{0,5}& \mathrm{0,2}\\ \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}& \mathrm{0,2}\\ \mathrm{0,2}& \mathrm{0,3}& \mathrm{0,1}& \mathrm{0,4}\\ \mathrm{0,1}& \mathrm{0,1}& \mathrm{0,5}& \mathrm{0,3}\end{array}\right)\left(\begin{array}{c}\mathrm{0,4}\\ \mathrm{0,}\text{05}\\ \mathrm{0,4}\\ \mathrm{0,4}\end{array}\right)=\left(\begin{array}{c}\mathrm{0,}\text{338}\\ \mathrm{0,}\text{365}\\ \mathrm{0,}\text{295}\\ \mathrm{0,}\text{365}\end{array}\right)=\left(\begin{array}{c}\mathrm{0,2}\mathrm{0,}\text{05}+\mathrm{0,4}\left(1-\mathrm{0,2}\right)\\ \mathrm{0,1}\mathrm{0,}\text{05}+\mathrm{0,4}\left(1-\mathrm{0,1}\right)\\ \mathrm{0,3}\mathrm{0,}\text{05}+\mathrm{0,4}\left(1-\mathrm{0,3}\right)\\ \mathrm{0,1}\mathrm{0,}\text{05}+\mathrm{0,4}\left(1-\mathrm{0,1}\right)\end{array}\right)\text{.}$$ (а).

Як випливає з (а), компоненти вектора $${\stackrel{}{x}}^{}$$ подаються у загальному вигляді.

$${\text{am}}_0+\left(1-a\right)M_0=M_0-a\left(M_0-m_0\right),$$.

де $$a>=\text{.}$$ Зокрема, для першого елемента вектора $${\stackrel{}{x}}^{}$$ $$a=\mathrm{0,2}-$$ $$1-a=1-\mathrm{0,2}=\mathrm{0,8}$$ і т. д.

Оскільки $$\stackrel{}{x}<={\stackrel{}{x}}^{},$$ то $$\stackrel{}{x}<={\stackrel{}{x}}^{},$$ тобто кожний компонент вектора $$\stackrel{}{x}$$ не перевищує відповідного компонента вектора $${\stackrel{}{x}}^{}\text{.}$$ А звідси випливає.

$$M_1<=M_0-\left(M_0-m_0\right)\text{.}$$ (b).

Оскільки $$M_0-m_0>=\mathrm{0,}$$ $$>=\mathrm{0,}$$ то $$M_1<=M_0\text{.}$$.

Аналогічно для вектора $$\stackrel{}{x}$$ маємо.

$$-m_1<=-m_0-\left(-m_0+M_0\right),$$ ©.

звідки випливає, що $$m_1>=m_0\text{.}$$.

Додавши (b) і ©, дістанемо.

$$M_1-m_1<=M_0-m_0-2\left(M_0-m_0\right)->$$.

$$->M_1-m_1<=\left(1-2\right)\left(M_0-m_0\right)\text{.}$$.

Тоді для вектора $${}^k\stackrel{}{x}$$ матимемо максимальний компонент $$M_k$$ і мінімальний $$m_k\text{.}$$.

При $$k=\mathrm{1,}\mathrm{2,}\mathrm{3,}\text{.}\text{.}\text{.}$$ дістанемо $$M_1>=M_2>=M_3>=\text{.}\text{.}\text{.},$$ $$m_1<=m_2<=m_3<=\text{.}\text{.}\text{.},$$ $$M_k-m_k<=\left(1-2\right)\left(M_{k-1}-m_{k-1}\right)$$, оскільки.

$$\begin{array}{c}{M}_{k-1}-m_{k-1}<=\left(1-2\right)\left(M_{k-2}-m_{k-2}\right)\\ M_{k-2}-m_{k-2}<=\left(1-2\right)\left(M_{k-3}-m_{k-3}\right),\\ \dots \dots \dots \\ M_1-m_1<=\left(1-2\right)\left(M_0-m_0\right)\text{.}\end{array}$$ (d).

Тепер згідно з (d) дістанемо.

$$M_k-m_k<={\left(1-2\right)}^{k-1}\left(M_0-m_0\right)\text{.}$$ (e).

З огляду на те, що $$0<<\mathrm{0,5,}$$ за $$k->\mathrm{}$$ маємо такий граничний перехід: $${\left(1-2\right)}^{k-1}->\mathrm{0,}$$ звідки.

$$M_k-m_k->0->M_k->m_k\text{.}$$ (f).

Отже, при $$k->\mathrm{}$$ $${}^k\stackrel{}{x}$$ прямує до вектора, усі компоненти якого дорівнюють $${}_i\mathrm{:}$$.

$$\underset{k->\mathrm{}}{\text{lim}}{}^k=W,$$ (44).

де $$W=\left(\begin{array}{cccc}{w}_1& w_2& \dots & w_N\\ w_1& w_2& \dots & w_N\\ & & \dots & \\ w_1& w& \dots & w_N\end{array}\right)$$ — матриця стаціонарних імовірностей.

Для регулярних ланцюгів Маркова довготривале поводження системи (процесу) не залежить від початкового стану, з якого система почала функціювати.

Приклад 21. Дано матрицю однокрокового переходу, що описує регулярний процес переходу системи з одного стану до іншого $$\left(\mathrm{}=\left({}_1,{}_2,{}_3,{}_4\right)\right)$$ :

$$\begin{array}{c}\begin{array}{cccc}{}_1& {}_2& {}_3& {}_4\end{array}\\ =\begin{array}{c}{}_1\\ {}_2\\ {}_3\\ {}_4\end{array}\left(\begin{array}{cccc}\mathrm{0,4}& \mathrm{0,1}& \mathrm{0,2}& \mathrm{0,3}\\ \mathrm{0,2}& \mathrm{0,4}& \mathrm{0,3}& \mathrm{0,1}\\ \mathrm{0,5}& 0& \mathrm{0,2}& \mathrm{0,3}\\ \mathrm{0,1}& \mathrm{0,2}& \mathrm{0,1}& \mathrm{0,6}\end{array}\right)\end{array}$$ .

Знайти ймовірність того, що у стаціонарному режимі система перебуватиме у стані $${}_2,$$ якщо перший крок вона зробила зі стану $${}_1-$$ зі стану $${}_3\text{.}$$.

Розв’язання. Визначивши граничну матрицю.

$$W=\underset{k->\mathrm{}}{\text{lim}}{}^k=\left(\begin{array}{cccc}\mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\end{array}\right)$$ ,.

дістанемо.

$$\left(\begin{array}{cccc}1& 0& 0& 0\end{array}\right)\left(\begin{array}{cccc}\mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\end{array}\right)=\left(\begin{array}{cccc}\mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\end{array}\right)$$ ;

$$\left(\begin{array}{cccc}0& 0& 1& 0\end{array}\right)\left(\begin{array}{cccc}\mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\\ \mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\end{array}\right)=\left(\begin{array}{cccc}\mathrm{0,}\text{27}& \mathrm{0,}\text{171}& \mathrm{0,}\text{179}& \mathrm{0,}\text{38}\end{array}\right)\text{.}$$.

Отже, незалежно від того, з якого стану — $${}_1$$

чи.

$${}_3$$

— почала функціонувати система, ймовірність набути стан.

$${}_2$$

буде однакова й дорівнюватиме 0,171.

.

2. Знаходження фундаментальної матриці.

Щоб обчислити очікувану кількість кроків, що їх має здійснити система (процес) для регулярного марковського ланцюга, аби повернутися до певного фіксованого можливого стану, використовують фундаментальну матрицю, яку умовно позначають Z.

Нехай однокрокова матриця ймовірностей переходу системи з одного стану до іншого, а W — її гранична матриця.

Розглянемо матрицю $$Q=\left(-W\right)$$

і, урахувавши, що для граничної матриці W виконується рівність.

$$W^k=W$$.

$${\left(-W\right)}^k={}^k-W\text{.}$$.

Оскільки $$\underset{k->\mathrm{}}{\text{lim}}{}^k=W,$$

то.

$$\underset{k->\mathrm{}}{\text{lim}}{\left(-W\right)}^k=\underset{k->\mathrm{}}{\text{lim}}{Q}^k=0$$

(нульова матриця). Зі сказаного випливає, що матриця.

$$I-Q=I-\left(-W\right)$$

має обернену, адже.

$$\text{det}\left(I-\left(-W\right)\right)/=0\text{.}$$.

Таким чином, доведено, що коли однокроковою матрицею ймовірностей переходу для регулярного ланцюга Маркова, то матриця.

$$Z={\left(I-\left(-W\right)\right)}^{-1}$$.

є фундаментальною для цього ланцюга.

Приклад 22. Для даної матриці.

$$=\left(\begin{array}{cccc}\mathrm{0,}\text{12}& \mathrm{0,}\text{18}& \mathrm{0,}\text{45}& \mathrm{0,}\text{25}\\ \mathrm{0,}\text{04}& \mathrm{0,}\text{26}& \mathrm{0,}\text{55}& \mathrm{0,}\text{15}\\ \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,}\text{05}& \mathrm{0,4}\\ \mathrm{0,}\text{43}& \mathrm{0,}\text{17}& \mathrm{0,}\text{28}& \mathrm{0,}\text{12}\end{array}\right)$$.

регулярного ланцюга Маркова знайти фундаментальну матрицю Z.

Розв’язання. Насамперед знайдемо граничну матрицю.

$$W=\underset{k->\mathrm{}}{\text{lim}}{}^k=\left(\begin{array}{cccc}\mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\end{array}\right)$$

.

.

а далі матрицю.

$$\left(I-+W\right)=$$.

$$=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right)-\left(\begin{array}{cccc}\mathrm{0,}\text{12}& \mathrm{0,}\text{18}& \mathrm{0,}\text{45}& \mathrm{0,}\text{25}\\ \mathrm{0,}\text{04}& \mathrm{0,}\text{26}& \mathrm{0,}\text{55}& \mathrm{0,}\text{15}\\ \mathrm{0,2}& \mathrm{0,}\text{35}& \mathrm{0,}\text{05}& \mathrm{0,4}\\ \mathrm{0,}\text{43}& \mathrm{0,}\text{17}& \mathrm{0,}\text{28}& \mathrm{0,}\text{12}\end{array}\right)+\left(\begin{array}{cccc}\mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\\ \mathrm{0,}\text{199}& \mathrm{0,}\text{25}& \mathrm{0,}\text{31}& \mathrm{0,}\text{241}\end{array}\right)=$$.

$$=\left(\begin{array}{cccc}\mathrm{1,}\text{079}& \mathrm{0,}\text{07}& -\mathrm{0,}\text{14}& -\mathrm{0,}\text{009}\\ \mathrm{0,}\text{159}& \mathrm{0,}\text{99}& -\mathrm{0,}\text{24}& \mathrm{0,}\text{091}\\ -\mathrm{0,}\text{001}& -\mathrm{0,1}& \mathrm{1,}\text{26}& -\mathrm{0,}\text{159}\\ -\mathrm{0,}\text{231}& \mathrm{0,}\text{08}& \mathrm{0,}\text{03}& \mathrm{1,}\text{121}\end{array}\right)\text{.}$$.

Тоді.

$$Z={\left(I-\left(-W\right)\right)}^{-1}={\left(\begin{array}{cccc}\mathrm{1,}\text{079}& \mathrm{0,}\text{07}& -\mathrm{0,}\text{14}& -\mathrm{0,}\text{009}\\ \mathrm{0,}\text{159}& \mathrm{0,}\text{99}& -\mathrm{0,}\text{24}& \mathrm{0,}\text{091}\\ -\mathrm{0,}\text{001}& -\mathrm{0,1}& \mathrm{1,}\text{26}& -\mathrm{0,}\text{159}\\ -\mathrm{0,}\text{231}& \mathrm{0,}\text{08}& \mathrm{0,}\text{03}& \mathrm{1,}\text{121}\end{array}\right)}^{-1}=$$.

$$=\left(\begin{array}{cccc}\mathrm{0,}\text{941}& -\mathrm{0,}\text{0592}& \mathrm{0,}\text{0927}& \mathrm{0,}\text{0255}\\ -\mathrm{0,}\text{1668}& \mathrm{1,}\text{0452}& \mathrm{0,}\text{1820}& -\mathrm{0,}\text{0604}\\ \mathrm{0,}\text{0134}& \mathrm{0,}\text{0717}& \mathrm{0,}\text{8062}& \mathrm{0,}\text{1086}\\ \mathrm{0,}\text{2055}& -\mathrm{0,}\text{887}& -\mathrm{0,}\text{0155}& \mathrm{0,}\text{8987}\end{array}\right)\text{.}$$.

Як бачимо, елементи матриці Z набувають від'ємних значень.

Застосовуючи фундаментальну матрицю Z для регулярних ланцюгів Маркова, обчислюють числові характеристики, що описують поводження цих ланцюгів.

Наприклад, середнє значення кількості кроків $$f_j,$$

які необхідно здійснити системі, щоб уперше набути стану.

$${}_j$$

за умови, що перед цим вона перебувала у стані.

$${}_i$$

дорівнюватиме елементу.

$$m_{\text{ij}}$$

матриці.

.

$$M=\left(I-Z+{\text{EZ}}_{\text{dg}}\right)D$$

(47).

.

що міститься на перетині і-го рядка та j-го стовпця.

Тут Е — одинична матриця (усі її елементи дорівнюють одиниці):

$$$$ $$E=\left(\begin{array}{ccccc}1& 1& 1& \dots & 1\\ 1& 1& 1& \dots & 1\\ \dots & \dots & \dots & \dots & \dots \\ \dots & \dots & \dots & \dots & \dots \\ 1& 1& 1& \dots & 1\end{array}\right)-$$.

$$Z_{\text{dg}}$$ — діагональна матриця, діагональні елементи якої збігаються з діагональними елементами матриці Z;

D — діагональна матриця, елементи якої - $$d_{\text{ii}}=\frac{1}{w_i},$$ де $$w_i$$ ;

і-й компонент стаціонарної ймовірності стаціонарного вектора $$\stackrel{}{W}\text{.}$$.

Дисперсія числа $$f_j$$ обчислюється так:

$$D\left(f_j\right)=M\left(2Z_{\text{dg}}D-I\right)+2\left(\text{ZM}-E{\left(\text{ZM}\right)}_{\text{dg}}\right)-M(f_i{)}_{\text{sq}}\text{.}$$.

Приклад 23. Статистична обробка спостережень метеослужби, здійснюваних улітку для певного міста України, дала такі результати:

• кщо певний день був теплим і безхмарним, то ймовірність, що така сама погода лишиться й наступного дня, дорівнює 0,6- імовірність того, що вона зміниться на вітряну погоду, дорівнює 0,25, а на дощову — 0,15;

• оли ж погода була вітряна, то ймовірність того, що вона такою самою і залишиться наступного дня, дорівнює 0,32, а ймовірність того, що вона зміниться на тиху сонячну погоду, дорівнює 0,46, і на дощову — 0,22;

• якщо погода була дощова, то ймовірність того, що вона не зміниться й наступного дня, дорівнює 0,26, а ймовірність того, що зміниться на вітряну або сонячну тиху погоду, становить відповідно 0,29 і 0,45.

Виконати такі завдання:

1) побудувати матрицю однокрокового переходу p>

2) визначити матриці $$M(f_j),$$ $$D(f_j)\text{.}$$.

Розв’язання. 1) Запишемо матрицю

$$\begin{array}{c}\begin{array}{ccc}{}_1& {}_2& {}_3\end{array}\\ =\begin{array}{c}{}_1\\ {}_2\\ {}_3\end{array}\left(\begin{array}{ccc}\mathrm{0,6}& \mathrm{0,}\text{25}& \mathrm{0,}\text{15}\\ \mathrm{0,}\text{46}& \mathrm{0,}\text{32}& \mathrm{0,}\text{22}\\ \mathrm{0,}\text{29}& \mathrm{0,}\text{45}& \mathrm{0,}\text{26}\end{array}\right),\end{array}$$.

де $${}_1$$, $${}_2$$ і $${}_3$$ — відповідно сонячна, вітряна та дощова погода.

2) Щоб визначити $$M\left(f_j\right)$$ та $$D\left(f_j\right)$$, потрібна фундаментальна матриця Z, яка подається через граничну матрицю.

$$W=\underset{k->\mathrm{}}{\text{lim}}{}^k=\left(\begin{array}{ccc}\mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\\ \mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\\ \mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\end{array}\right)\text{.}$$.

Тепер маємо:

$$\left(I-+W\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right)-\left(\begin{array}{ccc}\mathrm{0,6}& \mathrm{0,}\text{25}& \mathrm{0,}\text{15}\\ \mathrm{0,}\text{46}& \mathrm{0,}\text{32}& \mathrm{0,}\text{22}\\ \mathrm{0,}\text{29}& \mathrm{0,}\text{45}& \mathrm{0,}\text{26}\end{array}\right)+\left(\begin{array}{ccc}\mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\\ \mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\\ \mathrm{0,}\text{497}& \mathrm{0,}\text{31}& \mathrm{0,}\text{193}\end{array}\right)=$$.

$$=\left(\begin{array}{ccc}\mathrm{0,}\text{897}& \mathrm{0,}\text{06}& \mathrm{0,}\text{043}\\ \mathrm{0,}\text{037}& \mathrm{0,}\text{99}& -\mathrm{0,}\text{027}\\ \mathrm{0,}\text{207}& -\mathrm{0,}\text{14}& \mathrm{0,}\text{933}\end{array}\right),$$.

звідки.

$$Z={\left(I-+W\right)}^{-1}={\left(\begin{array}{ccc}\mathrm{0,}\text{897}& \mathrm{0,}\text{06}& \mathrm{0,}\text{043}\\ \mathrm{0,}\text{037}& \mathrm{0,}\text{99}& -\mathrm{0,}\text{027}\\ \mathrm{0,}\text{207}& -\mathrm{0,}\text{14}& \mathrm{0,}\text{933}\end{array}\right)}^{-1}=\left(\begin{array}{ccc}\mathrm{1,}\text{131}& -\mathrm{0,}\text{0762}& -\mathrm{0,}\text{0543}\\ -\mathrm{0,}\text{0493}& \mathrm{1,}\text{018}& \mathrm{0,}\text{0317}\\ -\mathrm{0,}\text{2582}& \mathrm{0,}\text{1696}& \mathrm{1,}\text{0886}\end{array}\right)\text{.}$$.

Для визначення $$M\left(f_j\right)$$ знаходимо:

$${\text{EZ}}_{\text{dg}}=\left(\begin{array}{ccc}1& 1& 1\\ 1& 1& 1\\ 1& 1& 1\end{array}\right)\left(\begin{array}{ccc}\mathrm{1,}\text{131}& 0& 0\\ 0& \mathrm{1,}\text{018}& 0\\ 0& 0& \mathrm{1,}\text{0886}\end{array}\right)=\left(\begin{array}{ccc}\mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\\ \mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\\ \mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\end{array}\right)-$$.

$$I-Z+{\text{EZ}}_{\text{dg}}=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right)-\left(\begin{array}{ccc}\mathrm{1,}\text{131}& \mathrm{0,}\text{0762}& -\mathrm{0,}\text{0543}\\ -\mathrm{0,}\text{0493}& \mathrm{1,}\text{018}& \mathrm{0,}\text{0317}\\ -\mathrm{0,}\text{2582}& \mathrm{0,}\text{1696}& \mathrm{1,}\text{0886}\end{array}\right)+$$.

$$+\left(\begin{array}{ccc}\mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\\ \mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\\ \mathrm{1,}\text{131}& \mathrm{1,}\text{018}& \mathrm{1,}\text{0886}\end{array}\right)=\left(\begin{array}{ccc}1& \mathrm{1,}\text{094}& \mathrm{1,}\text{14}\\ \mathrm{1,}\text{18}& 1& \mathrm{1,}\text{054}\\ \mathrm{1,}\text{389}& \mathrm{0,}\text{848}& 1\end{array}\right)-$$.

$$D=\left(\begin{array}{ccc}\mathrm{2,}\text{012}& 0& 0\\ 0& \mathrm{3,}\text{226}& 0\\ 0& 0& \mathrm{5,}\text{181}\end{array}\right)\text{.}$$.

Остаточно маємо:

$$M\left(f_i\right)=\left(\begin{array}{ccc}1& \mathrm{1,}\text{094}& \mathrm{1,}\text{14}\\ \mathrm{1,}\text{18}& 1& \mathrm{1,}\text{054}\\ \mathrm{1,}\text{389}& \mathrm{0,}\text{848}& 1\end{array}\right)\left(\begin{array}{ccc}\mathrm{2,}\text{012}& 0& 0\\ 0& \mathrm{3,}\text{226}& 0\\ 0& 0& \mathrm{5,}\text{181}\end{array}\right)=$$.

$$=\left(\begin{array}{ccc}\mathrm{2,}\text{012}& \mathrm{3,}\text{606}& \mathrm{5,}\text{905}\\ \mathrm{2,}\text{374}& \mathrm{3,}\text{296}& \mathrm{5,}\text{46}\\ \mathrm{2,}\text{795}& \mathrm{2,}\text{795}& \mathrm{5,}\text{18}\end{array}\right)\text{.}$$.

Отже, якщо сьогодні, наприклад, погода тиха й сонячна, то середня кількість днів до найближчого вітряного дня становитиме 3,6, а до похмурої дощової погоди — 6 (5,905).

Коли ж дощова погода, то середня кількість днів до тихої сонячної погоди наближено дорівнює 3 (2,795).

Щоб визначити $$D\left(f_i\right)$$, знайдемо:

$$\text{ZM}=\left(\begin{array}{ccc}\mathrm{1,}\text{131}& -\mathrm{0,}\text{0762}& -\mathrm{0,}\text{0543}\\ -\mathrm{0,}\text{0493}& \mathrm{1,}\text{018}& \mathrm{0,}\text{0317}\\ -\mathrm{0,}\text{2582}& \mathrm{0,}\text{1696}& \mathrm{1,}\text{0886}\end{array}\right)\left(\begin{array}{ccc}\mathrm{2,}\text{012}& \mathrm{3,}\text{606}& \mathrm{5,}\text{905}\\ \mathrm{2,}\text{374}& \mathrm{3,}\text{296}& \mathrm{5,}\text{46}\\ \mathrm{2,}\text{795}& \mathrm{2,}\text{795}& \mathrm{5,}\text{18}\end{array}\right)=$$.

$$=\left(\begin{array}{ccc}\mathrm{1,}\text{943}& \mathrm{3,}\text{675}& \mathrm{5,}\text{981}\\ \mathrm{2,}\text{406}& \mathrm{3,}\text{266}& \mathrm{5,}\text{431}\\ \mathrm{2,}\text{926}& \mathrm{2,}\text{671}& \mathrm{5,}\text{04}\end{array}\right)-$$.

$${\left(\text{ZM}\right)}_{\text{dg}}=\left(\begin{array}{ccc}\mathrm{1,}\text{943}& 0& 0\\ 0& \mathrm{3,}\text{266}& 0\\ 0& 0& \mathrm{5,}\text{04}\end{array}\right)-$$.

$$E{\left(\text{ZM}\right)}_{\text{dg}}=\left(\begin{array}{ccc}1& 1& 1\\ 1& 1& 1\\ 1& 1& 1\end{array}\right)\left(\begin{array}{ccc}\mathrm{1,}\text{943}& 0& 0\\ 0& \mathrm{3,}\text{266}& 0\\ 0& 0& \mathrm{5,}\text{04}\end{array}\right)=\left(\begin{array}{ccc}\mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\\ \mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\\ \mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\end{array}\right)-$$.

$$\text{ZM}-E{\left(\text{ZM}\right)}_{\text{dg}}=\left(\begin{array}{ccc}\mathrm{1,}\text{943}& \mathrm{3,}\text{675}& \mathrm{5,}\text{981}\\ \mathrm{2,}\text{406}& \mathrm{3,}\text{266}& \mathrm{5,}\text{431}\\ \mathrm{2,}\text{926}& \mathrm{2,}\text{671}& \mathrm{5,}\text{04}\end{array}\right)-\left(\begin{array}{ccc}\mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\\ \mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\\ \mathrm{1,}\text{943}& \mathrm{3,}\text{266}& \mathrm{5,}\text{04}\end{array}\right)=$$.

$$=\left(\begin{array}{ccc}0& \mathrm{0,}\text{409}& \mathrm{0,}\text{941}\\ \mathrm{0,}\text{463}& 0& \mathrm{0,}\text{391}\\ \mathrm{0,}\text{983}& -\mathrm{0,}\text{596}& 0\end{array}\right)\text{.}$$.

$$2\left(\text{ZM}-E(\text{ZM}{)}_{\text{dg}}\right)=\left(\begin{array}{ccc}0& \mathrm{0,}\text{818}& \mathrm{1,}\text{882}\\ \mathrm{0,}\text{926}& 0& \mathrm{0,}\text{782}\\ \mathrm{1,}\text{966}& -\mathrm{1,}\text{192}& 0\end{array}\right)$$ ;

$$\begin{array}{c}M\left(2Z_{\text{dg}}D-I\right)=\left(\begin{array}{ccc}\mathrm{2,}\text{012}& \mathrm{3,}\text{606}& \mathrm{5,}\text{905}\\ \mathrm{2,}\text{374}& \mathrm{3,}\text{296}& \mathrm{5,}\text{46}\\ \mathrm{2,}\text{795}& \mathrm{2,}\text{795}& \mathrm{5,}\text{18}\end{array}\right)(\left(\begin{array}{ccc}\mathrm{2,}\text{262}& 0& 0\\ 0& \mathrm{2,}\text{036}& 0\\ 0& 0& \mathrm{2,}\text{1772}\end{array}\right)x\\ \end{array}$$ $$\begin{array}{c}\left(\begin{array}{ccc}\mathrm{2,}\text{012}& 0& 0\\ 0& \mathrm{3,}\text{226}& 0\\ 0& 0& \mathrm{5,}\text{18}\end{array}\right)-\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right))=\left(\begin{array}{ccc}\mathrm{2,}\text{012}& \mathrm{3,}\text{606}& \mathrm{5,}\text{905}\\ \mathrm{2,}\text{374}& \mathrm{3,}\text{296}& \mathrm{5,}\text{46}\\ \mathrm{2,}\text{795}& \mathrm{2,}\text{795}& \mathrm{5,}\text{18}\end{array}\right)x\\ \end{array}$$.

$$\left(\begin{array}{ccc}\mathrm{3,}\text{551}& 0& 0\\ 0& \mathrm{5,}\text{568}& 0\\ 0& 0& \text{10},\text{278}\end{array}\right)=\left(\begin{array}{ccc}\mathrm{7,}\text{145}& \text{20},\text{078}& \text{60},\text{692}\\ \mathrm{8,}\text{43}& \text{12},\text{352}& \text{56},\text{118}\\ \mathrm{9,}\text{925}& \text{15},\text{563}& \text{53},\text{24}\end{array}\right)\text{.}$$.

Отже, остаточно дістаємо:

$$D\left(f_j\right)=M\left(2Z_{\text{dg}}D-I\right)+2\left(\text{ZM}-E{\left(\text{ZM}\right)}_{\text{dg}}\right)-M(f_i{)}_{\text{sq}}=$$.

$$\begin{array}{c}=\left(\begin{array}{ccc}\mathrm{7,}\text{145}& \text{20},\text{078}& \text{60},\text{692}\\ \mathrm{8,}\text{43}& \text{12},\text{352}& \text{56},\text{118}\\ \mathrm{9,}\text{925}& \text{15},\text{563}& \text{53},\text{24}\end{array}\right)-\left(\begin{array}{ccc}\mathrm{4,}\text{048}& \text{13}& \text{34},\text{869}\\ \mathrm{5,}\text{636}& \text{10},\text{864}& \text{29},\text{812}\\ \mathrm{7,}\text{812}& \mathrm{7,}\text{812}& \text{26},\text{832}\end{array}\right)=\\ =\left(\begin{array}{ccc}\mathrm{3,}\text{097}& \mathrm{7,}\text{896}& \text{27},\text{705}\\ \mathrm{3,}\text{72}& \mathrm{1,}\text{488}& \text{27},\text{089}\\ \mathrm{4,}\text{079}& \mathrm{6,}\text{599}& \text{26},\text{408}\end{array}\right)\text{.}\end{array}$$.