Глава 8

Теория Перрона--Фробениуса для неотрицательных матриц

[33/97%]
Показать
LaTeX
§
Задача 8.2.1

Проверьте теорему Перрона, вычислив собственные значения и собственные векторы для

A=(723183129). \mathbf{A} = \begin{pmatrix} 7 & 2 & 3 \\ 1 & 8 & 3 \\ 1 & 2 & 9 \end{pmatrix}.

Найдите правый вектор Перрона pp, а также левый вектор Перрона qTq^T.

?
Задача 8.2.2

Убедитесь, что утверждения (8.2.2)--(8.2.5) действительно верны.

Утверждения (8.2.2)--(8.2.5) — это свойства корня Перрона и вектора Перрона, установленные в тексте для положительной матрицы An×n>0A_{n\times n} > 0 с r=ρ(A)r = \rho (A): (8.2.2) r>0r > 0; (8.2.3) r∈σ(A)r \in \sigma (A) (называемое корнем Перрона); (8.2.4) alg multA(r)=1\text{alg mult}_A(r) = 1; (8.2.5) существует собственный вектор x>0x > 0 такой, что Ax=rxAx = rx.

?
Задача 8.2.3

Приведите подробности, объясняющие, почему вектор Перрона определён однозначно.

?
Задача 8.2.4

Найдите корень Перрона и вектор Перрона для

A=(1−αβα1−β), A = \begin{pmatrix} 1 - \alpha & \beta \\ \alpha & 1 - \beta \end{pmatrix},

где α+β=1\alpha + \beta = 1 при α,β>0\alpha , \beta > 0.

?
Задача 8.2.5

Пусть An×n>0A_{n\times n} > 0 и ρ(A)=r\rho (A) = r.

?
(a)

Объясните, почему существует предел lim⁡k→∞(A/r)k\lim_{k\to \infty } (A/r)^k.

(b)

Объясните, почему lim⁡k→∞(A/r)k=G>0\lim_{k\to \infty } (A/r)^k = G > 0 является проектором на N(A−rI)N(A - rI) вдоль R(A−rI)R(A - rI).

(c)

Объясните, почему rk⁡(G)=1\operatorname {rk}\left(G\right) = 1.

Задача 8.2.6

Докажите, что если каждая сумма по строкам (или по столбцам) матрицы An×n>0A_{n\times n} > 0 равна ρ\rho, то ρ(A)=ρ\rho (A) = \rho.

?
Задача 8.2.7

Докажите, что если An×n>0A_{n\times n} > 0, то

min⁡i∑j=1naij≤ρ(A)≤max⁡i∑j=1naij. \min _i \sum _{j=1}^{n} a_{ij} \leq \rho (A) \leq \max _i \sum _{j=1}^{n} a_{ij}.
?
Примечание.
?

Вспомните факт (установленный ранее в книге), что если 0≤B≤C0 \leq B \leq C поэлементно, то ρ(B)≤ρ(C)\rho (B) \leq \rho (C).

Задача 8.2.8

Чтобы показать, в какой мере условие положительности не может быть ослаблено в теореме Перрона, постройте примеры квадратных матриц AA таких, что A≥0A \geq 0, но A≯0A \not> 0 (т.е. у AA есть хотя бы один нулевой элемент), с r=ρ(A)∈σ(A)r = \rho (A) \in \sigma (A), демонстрирующие справедливость следующих утверждений. Для разных утверждений можно использовать разные примеры.

?
(a)

rr может быть равно 00.

(b)

alg multA(r)\text{alg mult}_A(r) может быть больше 11.

(c)

index(r)\text{index}(r) может быть больше 11.

(d)

N(A−rI)N(A - rI) может не содержать положительный собственный вектор.

(e)

rr может быть не единственным собственным значением на спектральной окружности.

Задача 8.2.9

Установите минимаксный вариант формулы Коллатца--Виланда, утверждающий, что корень Перрона для A>0A > 0 задаётся как r=min⁡x∈Pg(x)r = \min_{x \in \mathcal{P}} g(x), где

g(x)=max⁡1≤i≤n[Ax]ixiиP={x∣x>0}. g(x) = \max _{1 \leq i \leq n} \frac{[Ax]_i}{x_i} \quad \text{и} \quad \mathcal{P} = \left\{ x \mid x > 0\right\} .
?
Задача 8.2.10

Заметьте, что N={x∣x≥0 при x≠0}\mathcal{N} = \left\{ x \mid x \geq 0 \text{ при } x \neq 0\right\} используется в максиминном варианте формулы Коллатца--Виланда, а P={x∣x>0}\mathcal{P} = \left\{ x \mid x > 0\right\} — в минимаксном варианте из предыдущего упражнения. Приведите пример матрицы A>0A > 0, показывающий, что r≠min⁡x∈Ng(x)r \neq \min_{x \in \mathcal{N}} g(x), если g(x)g(x) определена как

g(x)=max⁡1≤i≤nxi≠0[Ax]ixi. g(x) = \max _{\substack {1 \leq i \leq n \\ x_i \neq 0}} \frac{[Ax]_i}{x_i}.
?
§
Задача 8.3.1

Пусть A=(010303020)A = \begin{pmatrix} 0 & 1 & 0 \\ 3 & 0 & 3 \\ 0 & 2 & 0 \end{pmatrix}.

?
(a)

Покажите, что AA неприводима.

(b)

Найдите корень Перрона и вектор Перрона для AA.

(c)

Найдите число собственных значений на спектральной окружности матрицы AA.

Задача 8.3.2

Предположим, что индекс импримитивности 5×55 \times 5 неотрицательной неприводимой матрицы AA равен h=3h = 3. Объясните, почему AA обязательно вырождена, причём alg multA(0)=2\text{alg mult}_A(0) = 2.

?
Задача 8.3.3

Предположим, что AA — неотрицательная матрица, обладающая положительным спектральным радиусом и соответствующим положительным собственным вектором. Обязана ли из этого следовать неприводимость AA?

?
Задача 8.3.4

Не вычисляя собственные значения или характеристический многочлен, объясните, почему σ(Pn)={1,ω,ω2,…,ωn−1}\sigma (P_n) = \{ 1, \omega , \omega^2, \ldots , \omega^{n-1}\}, где ω=e2πi/n\omega = e^{2\pi i/n} для

Pn=(010⋯0001⋯0⋮⋮⋱⋱⋮00⋯01100⋯0). P_n = \begin{pmatrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & 0 & 1 \\ 1 & 0 & 0 & \cdots & 0 \end{pmatrix}.
?
Задача 8.3.5

Определите, является ли

A=(0120000070200000920400010) A = \begin{pmatrix} 0 & 1 & 2 & 0 & 0 \\ 0 & 0 & 0 & 7 & 0 \\ 2 & 0 & 0 & 0 & 0 \\ 0 & 9 & 2 & 0 & 4 \\ 0 & 0 & 0 & 1 & 0 \end{pmatrix}

приводимой или неприводимой.

?
Задача 8.3.6

Определите, является ли матрица AA из упражнения 8.3.5 примитивной или импримитивной.

?
Задача 8.3.7

Матрица Sn×n≥0S_{n\times n} \geq 0, у которой суммы по строкам не превосходят 11, причём хотя бы одна сумма по строке строго меньше 11, называется субстохастической матрицей.

?
(a)

Объясните, почему ρ(S)≤1\rho (S) \leq 1 для любой субстохастической матрицы.

(b)

Докажите, что ρ(S)<1\rho (S) < 1 для любой неприводимой субстохастической матрицы.

Задача 8.3.8

Неотрицательная матрица, у которой каждая сумма по строке равна 11, называется стохастической матрицей (некоторые говорят "построчно стохастической"). Докажите, что если An×nA_{n\times n} неотрицательна и неприводима, причём r=ρ(A)r = \rho (A), то AA подобна rPrP для некоторой неприводимой стохастической матрицы PP.

?
Примечание.
?

Рассмотрите D=(p10⋯00p2⋯0⋮⋮⋱⋮00⋯pn)D = \begin{pmatrix} p_1 & 0 & \cdots & 0 \\ 0 & p_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & p_n \end{pmatrix}, где pkp_k — компоненты вектора Перрона для AA.

Задача 8.3.9

Виландт построил матрицу

Wn=(010⋯0001⋯0⋮⋮⋱⋱⋮00⋯01110⋯0) W_n = \begin{pmatrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & 0 & 1 \\ 1 & 1 & 0 & \cdots & 0 \end{pmatrix}

чтобы показать, что Wnn2−2n+2>0W_n^{n^2-2n+2} > 0, но [Wnn2−2n+1]11=0[W_n^{n^2-2n+1}]_{11} = 0. Проверьте, что это верно при n=4n = 4.

?
Задача 8.3.10

В модели популяции Лесли объясните, что происходит с вектором f(t)f(t) при t→∞t \to \infty в зависимости от того, выполняется ли r<1r < 1, r=1r = 1 или r>1r > 1.

?
Примечание.
?

Модель популяции Лесли разбивает популяцию самок на возрастные группы G1,…,GnG_1, \ldots , G_n, где bkb_k — коэффициент рождаемости, а sks_k — коэффициент выживаемости для группы GkG_k, а f(t)=(f1(t),…,fn(t))Tf(t) = (f_1(t), \ldots , f_n(t))^T — вектор размеров групп в момент времени tt, удовлетворяющий f(t+1)=Lf(t)f(t+1) = Lf(t) для матрицы Лесли LL; r=ρ(L)r = \rho (L) — корень Перрона матрицы LL, а p,qp, q — соответствующие векторы Перрона для LL и LTL^T.

Задача 8.3.11

Используя характеристическое уравнение, покажите, что матрица Лесли (из модели популяции Лесли) примитивна даже если b1=0b_1 = 0 (при условии, что все остальные bkb_k и sks_k положительны).

?
Задача 8.3.12

Матрица A∈Rn×nA \in \mathbb {R}^{n \times n} называется существенно положительной, если AA неприводима и aij≥0a_{ij} \geq 0 для всех i≠ji \neq j. Докажите, что каждое из следующих утверждений эквивалентно тому, что AA существенно положительна.

?
(a)

Существует такое α∈R\alpha \in \mathbb {R}, что A+αIA + \alpha I примитивна.

(b)

etA>0e^{tA} > 0 для всех t>0t > 0.

Задача 8.3.13

Пусть AA — существенно положительная матрица, как определено в задаче 8.3.12. Докажите, что каждое из следующих утверждений верно.

?
(a)

AA имеет собственную пару (ξ,x)(\xi , x), где ξ\xi вещественно и x>0x > 0.

(b)

Если λ\lambda — любое собственное значение AA, отличное от ξ\xi, то Re(λ)<ξ\mathrm{Re}\left(\lambda \right) < \xi.

(c)

ξ\xi увеличивается при увеличении любого элемента AA.

Задача 8.3.14

Пусть A≥0A \geq 0 — неприводимая матрица, и пусть aij(k)a_{ij}^{(k)} обозначает элементы AkA^k. Докажите, что AA примитивна тогда и только тогда, когда

ρ(A)=lim⁡k→∞[aij(k)]1/k. \rho (A) = \lim _{k\to \infty } \left[a_{ij}^{(k)}\right]^{1/k}.
?
§
Задача 8.4.1

Найдите стационарное распределение для

P=(1/4003/43/81/43/801/31/61/61/3001/21/2). P = \begin{pmatrix} 1/4 & 0 & 0 & 3/4 \\ 3/8 & 1/4 & 3/8 & 0 \\ 1/3 & 1/6 & 1/6 & 1/3 \\ 0 & 0 & 1/2 & 1/2 \end{pmatrix}.

Представляет ли это стационарное распределение предельное распределение в обычном смысле или только в смысле Чезаро?

?
Задача 8.4.2

Дважды стохастической матрицей называется неотрицательная матрица Pn×nP_{n\times n}, у которой суммы всех строк, а также суммы всех столбцов равны 11. Для неприводимой nn-состояний цепи Маркова, у которой матрица переходов дважды стохастична, какова доля времени в долгосрочной перспективе, проводимая в каждом состоянии? Какой вид имеют lim⁡k→∞(I+P+⋯+Pk−1)/k\lim_{k\to \infty } (I + P + \cdots + P^{k-1})/k и lim⁡k→∞Pk\lim_{k\to \infty } P^k (если он существует)?

?
Примечание.
?

Цель этого упражнения — показать, что дважды стохастические матрицы не представляют большого интереса с точки зрения теории цепей Маркова. Однако существует интересный теоретический результат (принадлежащий Дж. Биркгофу, 1946 г.), утверждающий, что множество n×nn \times n дважды стохастических матриц образует выпуклый многогранник в Rn×n\mathbb {R}^{n \times n} с матрицами перестановок в качестве вершин.

Задача 8.4.3

Объясните, почему rk⁡(I−P)=n−1\operatorname {rk}\left(I-P\right) = n-1 для любой неприводимой стохастической матрицы Pn×nP_{n\times n}. Приведите пример, показывающий, что это не обязательно верно для приводимых стохастических матриц.

?
Задача 8.4.4

Докажите, что левый перроновский вектор для неприводимой стохастической матрицы Pn×nP_{n\times n} (n>1n > 1) задаётся формулой

πT=1∑i=1nPi(P1,P2,…,Pn), \pi ^T = \frac{1}{\sum _{i=1}^{n} P_i}\left(P_1, P_2, \ldots , P_n\right),

где PiP_i — ii-й главный минор порядка n−1n-1 в I−PI-P.

?
Примечание.
?

Чему равно [adj⁡(A)]A[\operatorname {adj}(A)]A, если AA вырождена?

Задача 8.4.5

Пусть Pn×nP_{n\times n} — неприводимая стохастическая матрица, и пусть Qk×kQ_{k\times k} — главная подматрица матрицы I−PI - P, где 1≤k<n1 \leq k < n. Докажите, что ρ(Q)<1\rho (Q) < 1.

?
Задача 8.4.6

Пусть Pn×nP_{n\times n} — неприводимая стохастическая матрица, и пусть Qk×kQ_{k\times k} — главная подматрица матрицы I−PI - P, где 1≤k<n1 \leq k < n. Объясните, почему QQ является M-матрицей.

?
Примечание.
?

M-матрицей называется вещественная невырожденная матрица Qk×kQ_{k\times k}, такая что [Q]ij≤0[Q]_{ij} \leq 0 для всех i≠ji \neq j и Q−1≥0Q^{-1} \geq 0; эквивалентно, QQ является M-матрицей тогда и только тогда, когда существуют матрица B≥0B \geq 0 и вещественное число s>ρ(B)s > \rho (B), такие что Q=sI−BQ = sI - B.

Задача 8.4.7

Пусть Pn×nP_{n\times n} (n>1n > 1) — неприводимая стохастическая матрица. Объясните, почему все главные миноры порядка 1≤k<n1 \leq k < n в I−PI - P положительны.

?
Задача 8.4.8

Используйте те же предположения, что и для отказоустойчивой системы, описанной в разобранном примере с двумя независимыми контролями (пример 8.4.5), но возьмите три контроля, AA, BB и CC, вместо двух. Определите среднее время до отказа, начиная с трёх исправных контролей, с двух исправных и одного непроверенного контроля, а также с трёх непроверенных контролей.

?
Примечание.
?

Пример, на который даётся ссылка: система имеет два независимых контроля, которые могут предотвратить её разрушение, активируемые в дискретные моменты времени; система «под контролем», если хотя бы один контроль исправен, и разрушена, если все контроли отказывают одновременно. Контроль, который был исправен при одной активации, с вероятностью 90%90\% надёжен при следующей; контроль, который отказал и был заменён, надёжен при следующей активации лишь с вероятностью 60%60\%.

Задача 8.4.9

Мышь помещают в одну камеру ящика, изображённого на рисунке ниже, а кошку — в другую камеру. Каждую минуту двери камер открывают ровно настолько, чтобы позволить переход из одной камеры в соседнюю. В половине случаев, когда двери открываются, кошка не покидает занимаемую ею камеру. То же верно и для мыши. Когда кошка или мышь перемещается, дверь, через которую она проходит, выбирается случайным образом.

Рисунок 8.4.1: план ящика с мышью (камера # 1 соединена и с # 2, и с # 3, а # 2 соединена с # 3).Рисунок 8.4.1: план ящика с мышью (камера # 1 соединена и с # 2, и с # 3, а # 2 соединена с # 3).

?
(a)

Объясните, почему кошка и мышь в конце концов обязательно окажутся в одной и той же камере, и определите ожидаемое число шагов до этого момента.

(b)

Определите вероятность того, что кошка поймает мышь в камере jj для каждого j=1,2,3j = 1, 2, 3.