8.2

Положительные матрицы

[10/90%]
Показать
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}.
?