Глава 7

Собственные значения и собственные векторы

[169/96%]
Показать
LaTeX
§
Задача 7.1.1

Найдите собственные значения и собственные векторы следующих матриц.

A=(−10−71411),B=(21684148−8−32−18),C=(3−250140−15), \mathbf{A} = \begin{pmatrix} -10 & -7 \\ 14 & 11 \end{pmatrix}, \quad \mathbf{B} = \begin{pmatrix} 2 & 16 & 8 \\ 4 & 14 & 8 \\ -8 & -32 & -18 \end{pmatrix}, \quad \mathbf{C} = \begin{pmatrix} 3 & -2 & 5 \\ 0 & 1 & 4 \\ 0 & -1 & 5 \end{pmatrix}, D=(063−151−124),E=(300030003). \mathbf{D} = \begin{pmatrix} 0 & 6 & 3 \\ -1 & 5 & 1 \\ -1 & 2 & 4 \end{pmatrix}, \quad \mathbf{E} = \begin{pmatrix} 3 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 3 \end{pmatrix}.

Для каких из них, если такие есть, не хватает собственных векторов в том смысле, что не существует полного линейно независимого набора?

?
Задача 7.1.2

Не выполняя вычисления собственных значений и собственных векторов, определите, какие из следующих векторов являются собственными для

A=(−9−6−2−4−8−6−3−12015853221712), \mathbf{A} = \begin{pmatrix} -9 & -6 & -2 & -4 \\ -8 & -6 & -3 & -1 \\ 20 & 15 & 8 & 5 \\ 32 & 21 & 7 & 12 \end{pmatrix},

и для тех, которые являются собственными векторами, укажите соответствующее собственное значение.

?
(a)

(−1,1,0,1)T(-1,1,0,1)^T.

(b)

(1,0,−1,0)T(1,0,-1,0)^T.

(c)

(−1,0,2,2)T(-1,0,2,2)^T.

(d)

(0,1,−3,0)T(0,1,-3,0)^T.

Задача 7.1.3

Объясните, почему собственные значения треугольных и диагональных матриц

T=(t11t12⋯t1n0t22⋯t2n⋮⋮⋱⋮00⋯tnn)иD=(λ10⋯00λ2⋯0⋮⋮⋱⋮00⋯λn) \mathbf{T} = \begin{pmatrix} t_{11} & t_{12} & \cdots & t_{1n} \\ 0 & t_{22} & \cdots & t_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & t_{nn} \end{pmatrix} \quad \text{и} \quad \mathbf{D} = \begin{pmatrix} \lambda _1 & 0 & \cdots & 0 \\ 0 & \lambda _2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda _n \end{pmatrix}

представляют собой попросту диагональные элементы — tiit_{ii} и λi\lambda_i.

?
Задача 7.1.4

Для T=(AB0C)\mathbf{T} = \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{0} & \mathbf{C} \end{pmatrix} докажите, что det⁡(T−λI)=det⁡(A−λI)det⁡(C−λI)\operatorname {det}\left(\mathbf{T}-\lambda \mathbf{I}\right) = \operatorname {det}\left(\mathbf{A}-\lambda \mathbf{I}\right)\operatorname {det}\left(\mathbf{C}-\lambda \mathbf{I}\right), и заключите, что σ(AB0C)=σ(A)∪σ(C)\sigma \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{0} & \mathbf{C} \end{pmatrix} = \sigma (\mathbf{A}) \cup \sigma (\mathbf{C}) для квадратных A\mathbf{A} и C\mathbf{C}.

?
Задача 7.1.5

Найдите собственные векторы D=diag⁡(λ1,λ2,…,λn)\mathbf{D} = \operatorname {diag}(\lambda_1, \lambda_2, \ldots , \lambda_n). В частности, каково собственное подпространство, соответствующее λi\lambda_i?

?
Задача 7.1.6

Докажите, что 0∈σ(A)0 \in \sigma (\mathbf{A}) тогда и только тогда, когда A\mathbf{A} — вырожденная матрица.

?
Задача 7.1.7

Объясните, почему очевидно, что

An×n=(n11⋯11n1⋯111n⋯1⋮⋮⋮⋱⋮111⋯n) \mathbf{A}_{n\times n} = \begin{pmatrix} n & 1 & 1 & \cdots & 1 \\ 1 & n & 1 & \cdots & 1 \\ 1 & 1 & n & \cdots & 1 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 1 & 1 & \cdots & n \end{pmatrix}

не имеет нулевого собственного значения, и, следовательно, почему A\mathbf{A} невырождена.

?
Задача 7.1.8

Объясните, почему собственные значения A∗A\mathbf{A}^{*}\mathbf{A} и AA∗\mathbf{A}\mathbf{A}^{*} вещественны и неотрицательны для любой A∈Cm×n\mathbf{A} \in \mathcal{C}^{m\times n}.

Подсказка: рассмотрите ∥Ax→∥22/∥x→∥22\left\| \mathbf{A}\overrightarrow {x}\right\|_2^2 / \left\| \overrightarrow {x}\right\|_2^2. Когда собственные значения A∗A\mathbf{A}^{*}\mathbf{A} и AA∗\mathbf{A}\mathbf{A}^{*} строго положительны?

?
Задача 7.1.9
?
(a)

Если A\mathbf{A} невырождена, и если (λ,x→)(\lambda , \overrightarrow {x}) — собственная пара для A\mathbf{A}, покажите, что (λ−1,x→)(\lambda^{-1}, \overrightarrow {x}) — собственная пара для A−1\mathbf{A}^{-1}.

(b)

Для всех α∉σ(A)\alpha \notin \sigma (\mathbf{A}) докажите, что x→\overrightarrow {x} является собственным вектором A\mathbf{A} тогда и только тогда, когда x→\overrightarrow {x} является собственным вектором (A−αI)−1(\mathbf{A}-\alpha \mathbf{I})^{-1}.

Задача 7.1.10
?
(a)

Покажите, что если (λ,x→)(\lambda , \overrightarrow {x}) — собственная пара для A\mathbf{A}, то (λk,x→)(\lambda^k, \overrightarrow {x}) — собственная пара для Ak\mathbf{A}^k при каждом натуральном kk.

(b)

Если p(x)=α0+α1x+α2x2+⋯+αkxkp(x) = \alpha_0 + \alpha_1 x + \alpha_2 x^2 + \cdots + \alpha_k x^k — произвольный многочлен, то определим p(A)p(\mathbf{A}) как матрицу

p(A)=α0I+α1A+α2A2+⋯+αkAk. p(\mathbf{A}) = \alpha _0 \mathbf{I} + \alpha _1 \mathbf{A} + \alpha _2 \mathbf{A}^2 + \cdots + \alpha _k \mathbf{A}^k.

Покажите, что если (λ,x→)(\lambda , \overrightarrow {x}) — собственная пара для A\mathbf{A}, то (p(λ),x→)(p(\lambda ), \overrightarrow {x}) — собственная пара для p(A)p(\mathbf{A}).

Задача 7.1.11

Объясните, почему теорема Гершгорина влечёт, что

A=(10−200120−410−100500) \mathbf{A} = \begin{pmatrix} 1 & 0 & -2 & 0 \\ 0 & 12 & 0 & -4 \\ 1 & 0 & -1 & 0 \\ 0 & 5 & 0 & 0 \end{pmatrix}

должна иметь по крайней мере два вещественных собственных значения. Подтвердите этот факт, вычислив собственные значения A\mathbf{A}.

Теорема Гершгорина: если nn кругов Гершгорина Gi={z∈C∣∣z−aii∣≤∑j≠i∣aij∣}G_i = \left\{ z \in \mathcal{C} \mid \left|z - a_{ii}\right| \le \sum_{j \ne i} \left|a_{ij}\right|\right\}, соответствующих строкам матрицы A\mathbf{A} размера n×nn \times n, таковы, что kk из них образуют связную область, не пересекающуюся с оставшимися n−kn-k кругами, то ровно kk собственных значений A\mathbf{A} лежат в связной области, образованной этими kk кругами, а оставшиеся n−kn-k собственных значений лежат в области, образованной объединением остальных n−kn-k кругов.

?
Задача 7.1.12

Если A\mathbf{A} нильпотентна (Ak=0\mathbf{A}^k = \mathbf{0} при некотором kk), объясните, почему trace⁡(A)=0\operatorname {trace}(\mathbf{A}) = 0.

Указание: Чему равно σ(A)\sigma (\mathbf{A})?

?
Задача 7.1.13

Если x→1,x→2,…,x→k\overrightarrow {x}_1, \overrightarrow {x}_2, \ldots , \overrightarrow {x}_k — собственные векторы матрицы A\mathbf{A}, отвечающие одному и тому же собственному значению λ\lambda, объясните, почему любая ненулевая линейная комбинация

v→=α1x→1+α2x→2+⋯+αnx→n \overrightarrow {v} = \alpha _1 \overrightarrow {x}_1 + \alpha _2 \overrightarrow {x}_2 + \cdots + \alpha _n \overrightarrow {x}_n

также является собственным вектором матрицы A\mathbf{A}, отвечающим собственному значению λ\lambda.

?
Задача 7.1.14

Объясните, почему собственный вектор квадратной матрицы A\mathbf{A} не может отвечать двум различным собственным значениям матрицы A\mathbf{A}.

?
Задача 7.1.15

Пусть σ(An×n)=σ(Bn×n)\sigma (\mathbf{A}_{n\times n}) = \sigma (\mathbf{B}_{n\times n}). Гарантирует ли это, что A\mathbf{A} и B\mathbf{B} имеют одинаковый характеристический многочлен?

?
Задача 7.1.16

Постройте примеры 2×22 \times 2, доказывающие следующие утверждения.

?
(a)

λ∈σ(A)\lambda \in \sigma (\mathbf{A}) и μ∈σ(B)\centernot  ⟹  λ+μ∈σ(A+B)\mu \in \sigma (\mathbf{B}) \centernot \implies \lambda +\mu \in \sigma (\mathbf{A}+\mathbf{B}).

(b)

λ∈σ(A)\lambda \in \sigma (\mathbf{A}) и μ∈σ(B)\centernot  ⟹  λμ∈σ(AB)\mu \in \sigma (\mathbf{B}) \centernot \implies \lambda \mu \in \sigma (\mathbf{A}\mathbf{B}).

Задача 7.1.17

Пусть {λ1,λ2,…,λn}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_n\right\} — собственные значения матрицы An×n\mathbf{A}_{n\times n}, и пусть (λk,c→)(\lambda_k, \overrightarrow {c}) — некоторая собственная пара.

?
(a)

Для λ∉σ(A)\lambda \notin \sigma (\mathbf{A}) объясните, почему (A−λI)−1c→=c→/(λk−λ)(\mathbf{A}-\lambda \mathbf{I})^{-1}\overrightarrow {c} = \overrightarrow {c}/(\lambda_k-\lambda ).

(b)

Для произвольного вектора d→n×1\overrightarrow {d}_{n\times 1} докажите, что собственные значения матрицы A+c→d→T\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^T совпадают с собственными значениями A\mathbf{A}, за исключением того, что λk\lambda_k заменяется на λk+d→Tc→\lambda_k + \overrightarrow {d}^T\overrightarrow {c}.

(c)

Как можно выбрать d→\overrightarrow {d}, чтобы гарантировать, что собственные значения матрицы A+c→d→T\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^T и матрицы A\mathbf{A} совпадают, за исключением того, что λk\lambda_k заменяется заданным числом μ\mu?

Задача 7.1.18

Пусть A\mathbf{A} — квадратная матрица.

?
(a)

Объясните, почему A\mathbf{A} и AT\mathbf{A}^T имеют одинаковые собственные значения.

(b)

Объясните, почему λ∈σ(A)  ⟺  λ‾∈σ(A∗)\lambda \in \sigma (\mathbf{A}) \iff \overline{\lambda } \in \sigma (\mathbf{A}^{*}).

Указание: Вспомните упражнение 6.1.8.

(c)

Следует ли из этих результатов, что λ∈σ(A)  ⟺  λ‾∈σ(A)\lambda \in \sigma (\mathbf{A}) \iff \overline{\lambda } \in \sigma (\mathbf{A}), когда A\mathbf{A} — квадратная матрица вещественных чисел?

(d)

Ненулевой вектор-строка y→∗\overrightarrow {y}^{*} называется левым собственным вектором матрицы A\mathbf{A}, если существует скаляр μ∈C\mu \in \mathcal{C}, такой что y→∗(A−μI)=0→\overrightarrow {y}^{*}(\mathbf{A}-\mu \mathbf{I}) = \overrightarrow {0}. Объясните, почему μ\mu должно быть собственным значением матрицы A\mathbf{A} в «правом» смысле этого термина, когда A\mathbf{A} — квадратная матрица вещественных чисел.

Задача 7.1.19

Рассмотрим матрицы Am×n\mathbf{A}_{m\times n} и Bn×m\mathbf{B}_{n\times m}.

?
(a)

Объясните, почему AB\mathbf{A}\mathbf{B} и BA\mathbf{B}\mathbf{A} имеют одинаковый характеристический многочлен, если m=nm=n.

Указание: Вспомните упражнение 6.2.16.

(b)

Объясните, почему характеристические многочлены матриц AB\mathbf{A}\mathbf{B} и BA\mathbf{B}\mathbf{A} не могут совпадать при m≠nm \ne n, а затем объясните, почему σ(AB)\sigma (\mathbf{A}\mathbf{B}) и σ(BA)\sigma (\mathbf{B}\mathbf{A}) совпадают, за возможным исключением нулевого собственного значения.

Задача 7.1.20

Если AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}, докажите, что A\mathbf{A} и B\mathbf{B} имеют общий собственный вектор.

Указание: Для λ∈σ(A)\lambda \in \sigma (\mathbf{A}) пусть столбцы матрицы X\mathbf{X} образуют базис пространства N(A−λI)N(\mathbf{A}-\lambda \mathbf{I}), так что (A−λI)BX=0(\mathbf{A}-\lambda \mathbf{I})\mathbf{B}\mathbf{X} = \mathbf{0}. Объясните, почему существует такая матрица P\mathbf{P}, что BX=XP\mathbf{B}\mathbf{X} = \mathbf{X}\mathbf{P}, а затем рассмотрите произвольную собственную пару матрицы P\mathbf{P}.

?
Задача 7.1.21

Для фиксированных матриц Pm×m\mathbf{P}_{m\times m} и Qn×n\mathbf{Q}_{n\times n} пусть T\mathbf{T} — линейный оператор на Cm×n\mathcal{C}^{m\times n}, определённый формулой T(A)=PAQ\mathbf{T}(\mathbf{A}) = \mathbf{P}\mathbf{A}\mathbf{Q}.

?
(a)

Покажите, что если x→\overrightarrow {x} — правый собственный вектор матрицы P\mathbf{P}, а y→∗\overrightarrow {y}^{*} — левый собственный вектор матрицы Q\mathbf{Q}, то x→y→∗\overrightarrow {x}\overrightarrow {y}^{*} является собственным вектором оператора T\mathbf{T}.

(b)

Объясните, почему trace⁡(T)=trace⁡(P)trace⁡(Q)\operatorname {trace}(\mathbf{T}) = \operatorname {trace}(\mathbf{P})\operatorname {trace}(\mathbf{Q}).

Задача 7.1.22

Пусть D=diag⁡(λ1,λ2,…,λn)\mathbf{D} = \operatorname {diag}(\lambda_1, \lambda_2, \ldots , \lambda_n) — вещественная диагональная матрица, такая что λ1<λ2<⋯<λn\lambda_1 < \lambda_2 < \cdots < \lambda_n, и пусть v→n×1\overrightarrow {v}_{n\times 1} — столбец из вещественных ненулевых чисел.

?
(a)

Докажите, что если α\alpha вещественно и отлично от нуля, то λi\lambda_i не является собственным значением матрицы D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T. Покажите, что собственные значения матрицы D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T на самом деле задаются решениями секулярного уравнения f(ξ)=0f(\xi )=0, определённого формулой

f(ξ)=1+α∑i=1nvi2λi−ξ. f(\xi ) = 1 + \alpha \sum _{i=1}^n \frac{v_i^2}{\lambda _i - \xi }.

Для n=4n=4 и α>0\alpha >0 убедитесь, что график функции f(ξ)f(\xi ) имеет nn вертикальных асимптот при λi\lambda_i, причём ff возрастает от 00 до +∞+\infty на каждом интервале (λi,λi+1)(\lambda_i, \lambda_{i+1}), а ещё одна ветвь за пределами λn\lambda_n приближается к горизонтальной асимптоте на уровне 11, так что каждый интервал (λi,λi+1)(\lambda_i, \lambda_{i+1}) и луч (λn,∞)(\lambda_n, \infty ) содержат ровно один корень ξi\xi_i функции ff, и на этом основании заключите, что собственные значения матрицы D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T перемежаются (interlace) с собственными значениями D\mathbf{D}.

Рисунок 7.1.4Рисунок 7.1.4

(b)

Убедитесь, что (D−ξiI)−1v→(\mathbf{D}-\xi_i\mathbf{I})^{-1}\overrightarrow {v} является собственным вектором матрицы D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T, отвечающим собственному значению ξi\xi_i.

Задача 7.1.23

Пусть λ1,…,λn\lambda_1, \ldots , \lambda_n — корни многочлена p(λ)=λn+c1λn−1+c2λn−2+⋯+cnp(\lambda ) = \lambda^n + c_1\lambda^{n-1} + c_2\lambda^{n-2} + \cdots + c_n, и пусть τk=λ1k+λ2k+⋯+λnk\tau_k = \lambda_1^k + \lambda_2^k + \cdots + \lambda_n^k. Тождества Ньютона утверждают, что ck=−(τ1ck−1+τ2ck−2+⋯+τk−1c1+τk)/kc_k = -(\tau_1 c_{k-1} + \tau_2 c_{k-2} + \cdots + \tau_{k-1}c_1 + \tau_k)/k. Выведите эти тождества, выполнив следующие шаги:

?
(a)

Покажите, что p′(λ)=p(λ)∑i=1n(λ−λi)−1p'(\lambda ) = p(\lambda )\sum_{i=1}^n (\lambda -\lambda_i)^{-1} (логарифмическое дифференцирование).

(b)

Используя разложение в геометрический ряд для (λ−λi)−1(\lambda -\lambda_i)^{-1}, покажите, что при ∣λ∣>max⁡i∣λi∣\left|\lambda \right| > \max_i \left|\lambda_i\right|

∑i=1n1λ−λi=nλ+τ1λ2+τ2λ3+⋯ . \sum _{i=1}^n \frac{1}{\lambda -\lambda _i} = \frac{n}{\lambda } + \frac{\tau _1}{\lambda ^2} + \frac{\tau _2}{\lambda ^3} + \cdots .
(c)

Объедините эти два результата и приравняйте коэффициенты при одинаковых степенях λ\lambda.

Задача 7.1.24

Пусть характеристическое уравнение для A\mathbf{A} задано в виде λn+c1λn−1+c2λn−2+⋯+cn=0\lambda^n + c_1\lambda^{n-1} + c_2\lambda^{n-2} + \cdots + c_n = 0, и определим последовательность, полагая B0=I\mathbf{B}_0 = \mathbf{I} и

Bk=−trace⁡(ABk−1)kI+ABk−1для k=1,2,…,n. \mathbf{B}_k = -\frac{\operatorname {trace}(\mathbf{A}\mathbf{B}_{k-1})}{k}\mathbf{I} + \mathbf{A}\mathbf{B}_{k-1} \quad \text{для } k = 1, 2, \ldots , n.

Докажите, что при каждом kk

ck=−trace⁡(ABk−1)k. c_k = -\frac{\operatorname {trace}(\mathbf{A}\mathbf{B}_{k-1})}{k}.

Указание: используйте тождества Ньютона (упражнение 7.1.23) и вспомните упражнение 7.1.10(a).

?
§
Задача 7.2.1

Приведите A=(−8−61210)\mathbf{A} = \begin{pmatrix} -8 & -6 \\ 12 & 10 \end{pmatrix} к диагональному виду с помощью преобразования подобия, либо объясните, почему A\mathbf{A} нельзя диагонализовать.

?
Задача 7.2.2
?
(a)

Проверьте, что alg mult⁡A(λ)=geo mult⁡A(λ)\operatorname {alg\, mult}_{\mathbf{A}}(\lambda ) = \operatorname {geo\, mult}_{\mathbf{A}}(\lambda ) для каждого собственного значения матрицы

A=(−4−3−30−10665). \mathbf{A} = \begin{pmatrix} -4 & -3 & -3 \\ 0 & -1 & 0 \\ 6 & 6 & 5 \end{pmatrix}.
(b)

Найдите невырожденную матрицу P\mathbf{P} такую, что P−1AP\mathbf{P}^{-1}\mathbf{A}\mathbf{P} — диагональная матрица.

Задача 7.2.3

Покажите, что подобные матрицы не обязаны иметь одни и те же собственные векторы, приведя пример двух матриц, которые подобны, но имеют разные собственные подпространства.

?
Задача 7.2.4

λ=2\lambda =2 — собственное значение для A=(321020−2−30)\mathbf{A} = \begin{pmatrix} 3 & 2 & 1 \\ 0 & 2 & 0 \\ -2 & -3 & 0 \end{pmatrix}. Найдите alg mult⁡A(λ)\operatorname {alg\, mult}_{\mathbf{A}}(\lambda ), а также geo mult⁡A(λ)\operatorname {geo\, mult}_{\mathbf{A}}(\lambda ). Можно ли на основании этих результатов сделать какой-либо вывод о диагонализуемости A\mathbf{A}?

?
Задача 7.2.5

Если B=P−1AP\mathbf{B} = \mathbf{P}^{-1}\mathbf{A}\mathbf{P}, объясните, почему Bk=P−1AkP\mathbf{B}^k = \mathbf{P}^{-1}\mathbf{A}^k\mathbf{P}.

?
Задача 7.2.6

Вычислите lim⁡n→∞An\lim_{n\to \infty } \mathbf{A}^n для A=(7/51/5−11/2)\mathbf{A} = \begin{pmatrix} 7/5 & 1/5 \\ -1 & 1/2 \end{pmatrix}.

?
Задача 7.2.7

Пусть {x→1,x→2,…,x→t}\left\{ \overrightarrow {x}_1, \overrightarrow {x}_2, \ldots , \overrightarrow {x}_t\right\} — набор линейно независимых собственных векторов для An×n\mathbf{A}_{n\times n}, отвечающих соответственно собственным значениям {λ1,λ2,…,λt}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_t\right\}, и пусть X\mathbf{X} — произвольная матрица размера n×(n−t)n \times (n-t) такая, что Pn×n=(x→1∣⋯∣x→t∣X)\mathbf{P}_{n\times n} = (\overrightarrow {x}_1 \mid \cdots \mid \overrightarrow {x}_t \mid \mathbf{X}) невырождена. Докажите, что если P−1=(y→1∗⋮y→t∗Y∗)\mathbf{P}^{-1} = \begin{pmatrix} \overrightarrow {y}_1^{*} \\ \vdots \\ \overrightarrow {y}_t^{*} \\ \mathbf{Y}^{*} \end{pmatrix}, где y→i∗\overrightarrow {y}_i^{*} — строки, а Y∗\mathbf{Y}^{*} имеет размер (n−t)×n(n-t) \times n, то {y→1∗,y→2∗,…,y→t∗}\left\{ \overrightarrow {y}_1^{*}, \overrightarrow {y}_2^{*}, \ldots , \overrightarrow {y}_t^{*}\right\} — набор линейно независимых левых собственных векторов, отвечающих соответственно {λ1,λ2,…,λt}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_t\right\} (т.е. y→i∗A=λiy→i∗\overrightarrow {y}_i^{*}\mathbf{A} = \lambda_i \overrightarrow {y}_i^{*}).

?
Задача 7.2.8

Пусть A\mathbf{A} — диагонализуемая матрица, и пусть ρ(⋆)\rho (\star ) обозначает спектральный радиус. Докажите, что lim⁡k→∞Ak=0\lim_{k\to \infty } \mathbf{A}^k = \mathbf{0} тогда и только тогда, когда ρ(A)<1\rho (\mathbf{A}) < 1.

Примечание: этот результат справедлив и для недиагонализуемых матриц, как показано в другом месте текста.

?
Задача 7.2.9

Примените приём, используемый при доказательстве теоремы Шура о триангуляризации, чтобы построить ортогональную матрицу P\mathbf{P} такую, что PTAP\mathbf{P}^T\mathbf{A}\mathbf{P} верхнетреугольна для A=(13−916−11)\mathbf{A} = \begin{pmatrix} 13 & -9 \\ 16 & -11 \end{pmatrix}.

?
Задача 7.2.10

Проверьте теорему Кэли--Гамильтона для A=(1−4−48−11−8−885)\mathbf{A} = \begin{pmatrix} 1 & -4 & -4 \\ 8 & -11 & -8 \\ -8 & 8 & 5 \end{pmatrix}.

Указание: это матрица из примера 7.2.1 на стр. 507.

?
Задача 7.2.11

Поскольку сумма элементов каждой строки следующей симметричной матрицы A\mathbf{A} равна 44, ясно, что x→=(1,1,1,1)T\overrightarrow {x} = (1,1,1,1)^T является одновременно правым и левым собственным вектором, отвечающим λ=4∈σ(A)\lambda =4 \in \sigma (\mathbf{A}). Используя метод дефляции (перенормируйте x→\overrightarrow {x} до единичной длины y→\overrightarrow {y}, постройте ортогональную матрицу P=[y→∣X]\mathbf{P} = [\overrightarrow {y} \mid \mathbf{X}], у которой y→\overrightarrow {y} является первым столбцом, тогда оставшиеся собственные значения A\mathbf{A} в точности совпадают с собственными значениями меньшей матрицы XTAX\mathbf{X}^T\mathbf{A}\mathbf{X}), найдите оставшиеся собственные значения матрицы

A=(1021021121101102). \mathbf{A} = \begin{pmatrix} 1 & 0 & 2 & 1 \\ 0 & 2 & 1 & 1 \\ 2 & 1 & 1 & 0 \\ 1 & 1 & 0 & 2 \end{pmatrix}.
?
Задача 7.2.12

Объясните, почему AGi=GiA=λiGi\mathbf{A}\mathbf{G}_i = \mathbf{G}_i\mathbf{A} = \lambda_i \mathbf{G}_i для спектрального проектора Gi\mathbf{G}_i, отвечающего собственному значению λi\lambda_i диагонализируемой матрицы A\mathbf{A}.

?
Задача 7.2.13

Докажите, что A=c→n×1d→1×nT\mathbf{A} = \overrightarrow {c}_{n\times 1}\overrightarrow {d}_{1\times n}^T диагонализируема тогда и только тогда, когда d→Tc→≠0\overrightarrow {d}^T\overrightarrow {c} \ne 0.

?
Задача 7.2.14

Докажите, что A=(W00Z)\mathbf{A} = \begin{pmatrix} \mathbf{W} & \mathbf{0} \\ \mathbf{0} & \mathbf{Z} \end{pmatrix} диагонализируема тогда и только тогда, когда Ws×s\mathbf{W}_{s\times s} и Zt×t\mathbf{Z}_{t\times t} каждая диагонализируема.

?
Задача 7.2.15

Докажите, что если AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}, то A\mathbf{A} и B\mathbf{B} можно одновременно привести к треугольному виду унитарным преобразованием подобия --- т.е. U∗AU=T1\mathbf{U}^{*}\mathbf{A}\mathbf{U} = \mathbf{T}_1 и U∗BU=T2\mathbf{U}^{*}\mathbf{B}\mathbf{U} = \mathbf{T}_2 для некоторой унитарной матрицы U\mathbf{U}.

Указание: вспомните упражнение 7.1.20 вместе с ходом доказательства теоремы Шура о триангуляризации.

?
Задача 7.2.16

Для диагонализируемых матриц докажите, что AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A} тогда и только тогда, когда A\mathbf{A} и B\mathbf{B} можно одновременно диагонализировать --- т.е. P−1AP=D1\mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \mathbf{D}_1 и P−1BP=D2\mathbf{P}^{-1}\mathbf{B}\mathbf{P} = \mathbf{D}_2 для некоторой P\mathbf{P}.

Указание: если A\mathbf{A} и B\mathbf{B} коммутируют, то коммутируют и P−1AP=(λ1I00D)\mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \begin{pmatrix} \lambda_1\mathbf{I} & \mathbf{0} \\ \mathbf{0} & \mathbf{D} \end{pmatrix}, P−1BP=(WXYZ)\mathbf{P}^{-1}\mathbf{B}\mathbf{P} = \begin{pmatrix} \mathbf{W} & \mathbf{X} \\ \mathbf{Y} & \mathbf{Z} \end{pmatrix}.

?
Задача 7.2.17

Объясните, почему следующее <<доказательство>> теоремы Кэли--Гамильтона неверно: p(λ)=det⁡(A−λI)  ⟹  p(A)=det⁡(A−AI)=det⁡(0)=0p(\lambda ) = \operatorname {det}\left(\mathbf{A}-\lambda \mathbf{I}\right) \implies p(\mathbf{A}) = \operatorname {det}\left(\mathbf{A}-\mathbf{A}\mathbf{I}\right) = \operatorname {det}\left(\mathbf{0}\right) = 0.

?
Задача 7.2.18

Покажите, что собственные значения матрицы конечных разностей

A=(2−1−12−1⋱⋱⋱−12−1−12)n×n \mathbf{A} = \begin{pmatrix} 2 & -1 & & & \\ -1 & 2 & -1 & & \\ & \ddots & \ddots & \ddots & \\ & & -1 & 2 & -1 \\ & & & -1 & 2 \end{pmatrix}_{n\times n}

равны λj=4sin⁡2jπ2(n+1)\lambda_j = 4\sin^2\dfrac {j\pi }{2(n+1)}, 1≤j≤n1 \le j \le n.

Пример 7.2.5 устанавливает, что собственные значения этой матрицы конечных разностей равны λj=2−2cos⁡jπn+1\lambda_j = 2 - 2\cos \dfrac {j\pi }{n+1}, 1≤j≤n1 \le j \le n.

?
Задача 7.2.19

Пусть N=(01⋱⋱⋱10)n×n\mathbf{N} = \begin{pmatrix} 0 & 1 & & \\ & \ddots & \ddots & \\ & & \ddots & 1 \\ & & & 0 \end{pmatrix}_{n\times n}.

?
(a)

Покажите, что λ∈σ(N+NT)\lambda \in \sigma (\mathbf{N}+\mathbf{N}^T) тогда и только тогда, когда iλ∈σ(N−NT)\mathrm{i}\lambda \in \sigma (\mathbf{N}-\mathbf{N}^T).

Пример 7.2.5 показывает, что собственные значения N+NT\mathbf{N}+\mathbf{N}^T равны λj=2cos⁡(jπ/(n+1))\lambda_j = 2\cos (j\pi /(n+1)), 1≤j≤n1 \le j \le n.

(b)

Объясните, почему N+NT\mathbf{N}+\mathbf{N}^T невырождена тогда и только тогда, когда nn чётно.

(c)

Вычислите det⁡(N−NT)/det⁡(N+NT)\operatorname {det}\left(\mathbf{N}-\mathbf{N}^T\right)/\operatorname {det}\left(\mathbf{N}+\mathbf{N}^T\right), когда nn чётно.

Задача 7.2.20

Матрица Тёплица вида

C=(c0cn−1cn−2⋯c1c1c0cn−1⋯c2c2c1c0⋯c3⋮⋮⋮⋱⋮cn−1cn−2cn−3⋯c0)n×n \mathbf{C} = \begin{pmatrix} c_0 & c_{n-1} & c_{n-2} & \cdots & c_1 \\ c_1 & c_0 & c_{n-1} & \cdots & c_2 \\ c_2 & c_1 & c_0 & \cdots & c_3 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ c_{n-1} & c_{n-2} & c_{n-3} & \cdots & c_0 \end{pmatrix}_{n\times n}

называется циркулянтной матрицей. Если p(x)=c0+c1x+⋯+cn−1xn−1p(x) = c_0+c_1 x+\cdots +c_{n-1}x^{n-1} и {1,ξ,ξ2,…,ξn−1}\left\{ 1,\xi ,\xi^2,\ldots ,\xi^{n-1}\right\} — корни nn-й степени из единицы, то результаты упражнения 5.8.12 (с.379) гарантируют, что

FnCFn−1=(p(1)0⋯00p(ξ)⋯0⋮⋮⋱⋮00⋯p(ξn−1)) \mathbf{F}_n\mathbf{C}\mathbf{F}_n^{-1} = \begin{pmatrix} p(1) & 0 & \cdots & 0 \\ 0 & p(\xi ) & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & p(\xi ^{n-1}) \end{pmatrix}

где Fn\mathbf{F}_n — матрица Фурье порядка nn. Проверьте эти факты для приведённого ниже циркулянта, вычислив его собственные значения и собственные векторы непосредственно:

C=(1010010110100101). \mathbf{C} = \begin{pmatrix} 1 & 0 & 1 & 0 \\ 0 & 1 & 0 & 1 \\ 1 & 0 & 1 & 0 \\ 0 & 1 & 0 & 1 \end{pmatrix}.
?
Задача 7.2.21

Предположим, что (λ,x→)(\lambda , \overrightarrow {x}) и (μ,y→∗)(\mu , \overrightarrow {y}^{*}) — правая и левая собственные пары матрицы A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} --- т.е. Ax→=λx→\mathbf{A}\overrightarrow {x}=\lambda \overrightarrow {x} и y→∗A=μy→∗\overrightarrow {y}^{*}\mathbf{A}=\mu \overrightarrow {y}^{*}. Объясните, почему y→∗x→=0\overrightarrow {y}^{*}\overrightarrow {x}=0 всякий раз, когда λ≠μ\lambda \ne \mu.

?
Задача 7.2.22

Рассмотрим A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n}.

?
(a)

Покажите, что если A\mathbf{A} диагонализуема, то существуют правый и левый собственные векторы x→\overrightarrow {x} и y→∗\overrightarrow {y}^{*}, отвечающие λ∈σ(A)\lambda \in \sigma (\mathbf{A}), такие что y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0, так что можно добиться y→∗x→=1\overrightarrow {y}^{*}\overrightarrow {x}=1.

(b)

Покажите, что не для всяких правого и левого собственных векторов x→\overrightarrow {x} и y→∗\overrightarrow {y}^{*}, отвечающих λ∈σ(A)\lambda \in \sigma (\mathbf{A}), выполняется y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0.

(c)

Покажите, что утверждение (а) может не выполняться, если A\mathbf{A} не диагонализуема.

Задача 7.2.23

Рассмотрим A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} с λ∈σ(A)\lambda \in \sigma (\mathbf{A}).

?
(a)

Докажите, что если λ\lambda простое, то y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0 для любой пары соответствующих правого и левого собственных векторов x→\overrightarrow {x} и y→∗\overrightarrow {y}^{*}, отвечающих λ\lambda, независимо от того, диагонализуема ли A\mathbf{A}.

Указание: Используйте разложение на ядровую и нильпотентную части (с. 397).

(b)

Покажите, что y→∗x→=0\overrightarrow {y}^{*}\overrightarrow {x}=0 возможно, когда λ\lambda не простое.

Задача 7.2.24

Для A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} с σ(A)={λ1,λ2,…,λk}\sigma (\mathbf{A}) = \left\{ \lambda_1,\lambda_2,\ldots ,\lambda_k\right\} покажите, что A\mathbf{A} диагонализуема тогда и только тогда, когда Rn=N(A−λ1I)⊕N(A−λ2I)⊕⋯⊕N(A−λkI)\mathbb {R}^{n} = N(\mathbf{A}-\lambda_1\mathbf{I}) \oplus N(\mathbf{A}-\lambda_2\mathbf{I}) \oplus \cdots \oplus N(\mathbf{A}-\lambda_k\mathbf{I}).

Указание: Вспомните упражнение 5.9.14.

?
Задача 7.2.25

Теорема Шура о триангуляризации гарантирует, что всякая квадратная матрица A\mathbf{A} унитарно подобна верхнетреугольной матрице — скажем, U∗AU=T\mathbf{U}^{*}\mathbf{A}\mathbf{U}=\mathbf{T}. Но даже если A\mathbf{A} вещественна, U\mathbf{U} и T\mathbf{T} могут оказаться комплексными, если A\mathbf{A} имеет комплексные собственные значения. Тем не менее, матрицы (и арифметику) можно сделать вещественными, если удовлетвориться блочно-треугольным результатом с блоками 2×22\times 2 или скалярными элементами на диагонали. Докажите, что для каждой A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} существует ортогональная матрица P∈Rn×n\mathbf{P} \in \mathbb {R}^{n\times n} и вещественные матрицы Bij\mathbf{B}_{ij} такие, что

PTAP=(B11B12⋯B1k0B22⋯B2k⋮⋮⋱⋮00⋯Bkk),где Bjj имеет размер 1×1 или 2×2. \mathbf{P}^T\mathbf{A}\mathbf{P} = \begin{pmatrix} \mathbf{B}_{11} & \mathbf{B}_{12} & \cdots & \mathbf{B}_{1k} \\ \mathbf{0} & \mathbf{B}_{22} & \cdots & \mathbf{B}_{2k} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{0} & \mathbf{0} & \cdots & \mathbf{B}_{kk} \end{pmatrix}, \quad \text{где } \mathbf{B}_{jj} \text{ имеет размер } 1\times 1 \text{ или } 2\times 2.

Если Bjj=[λj]\mathbf{B}_{jj} = [\lambda_j] имеет размер 1×11\times 1, то λj∈σ(A)\lambda_j \in \sigma (\mathbf{A}), а если Bjj\mathbf{B}_{jj} имеет размер 2×22\times 2, то σ(Bjj)={λj,λj‾}⊆σ(A)\sigma (\mathbf{B}_{jj}) = \left\{ \lambda_j, \overline{\lambda_j}\right\} \subseteq \sigma (\mathbf{A}).

?
Задача 7.2.26

Когда A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} диагонализуема преобразованием подобия S\mathbf{S}, матрица S\mathbf{S} может оказаться комплексной, если A\mathbf{A} имеет комплексные собственные значения. По аналогии с упражнением 7.2.25, можно остаться в области вещественных чисел, удовлетворившись блочно-диагональным результатом с элементами 1×11\times 1 или 2×22\times 2 на диагонали. Докажите, что если A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} диагонализуема с вещественными собственными значениями {ρ1,…,ρr}\left\{ \rho_1,\ldots ,\rho_r\right\} и комплексными собственными значениями {λ1,λ1‾,λ2,λ2‾,…,λt,λt‾}\left\{ \lambda_1,\overline{\lambda_1},\lambda_2,\overline{\lambda_2},\ldots ,\lambda_t,\overline{\lambda_t}\right\} с 2t+r=n2t+r=n, то существует невырожденная P∈Rn×n\mathbf{P} \in \mathbb {R}^{n\times n} и матрицы Bj∈R2×2\mathbf{B}_j \in \mathbb {R}^{2\times 2} такие, что

P−1AP=(D0⋯00B1⋯0⋮⋮⋱⋮00⋯Bt),где D=(ρ10⋯00ρ2⋯0⋮⋮⋱⋮00⋯ρr), \mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \begin{pmatrix} \mathbf{D} & \mathbf{0} & \cdots & \mathbf{0} \\ \mathbf{0} & \mathbf{B}_1 & \cdots & \mathbf{0} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{0} & \mathbf{0} & \cdots & \mathbf{B}_t \end{pmatrix}, \quad \text{где } \mathbf{D} = \begin{pmatrix} \rho _1 & 0 & \cdots & 0 \\ 0 & \rho _2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \rho _r \end{pmatrix},

а Bj\mathbf{B}_j имеет собственные значения λj\lambda_j и λj‾\overline{\lambda_j}.

?
Задача 7.2.27

Для A∈Cn×n\mathbf{A} \in \mathcal{C}^{n\times n} докажите, что x→∗Ax→=0\overrightarrow {x}^{*}\mathbf{A}\overrightarrow {x}=0 для всех x→∈Cn×1  ⟹  A=0\overrightarrow {x} \in \mathcal{C}^{n\times 1} \implies \mathbf{A}=\mathbf{0}. Покажите, что x→TAx→=0\overrightarrow {x}^T\mathbf{A}\overrightarrow {x}=0 для всех x→∈Rn×1\centernot  ⟹  A=0\overrightarrow {x} \in \mathbb {R}^{n\times 1} \centernot \implies \mathbf{A}=\mathbf{0}, даже если A\mathbf{A} вещественна.

?
§
Задача 7.3.1

Определите cos⁡A\cos \mathbf{A} для A=(−π/2π/2π/2−π/2)\mathbf{A} = \begin{pmatrix} -\pi /2 & \pi /2 \\ \pi /2 & -\pi /2 \end{pmatrix}.

?
Задача 7.3.2

Для матрицы A\mathbf{A} из примера 7.3.3 проверьте прямым вычислением, что eλ1tG1+eλ2tG2=P(eλ1t00eλ2t)P−1=eAt\mathrm{e}^{\lambda_1 t}\mathbf{G}_1 + \mathrm{e}^{\lambda_2 t}\mathbf{G}_2 = \mathbf{P}\begin{pmatrix} \mathrm{e}^{\lambda_1 t} & 0 \\ 0 & \mathrm{e}^{\lambda_2 t} \end{pmatrix}\mathbf{P}^{-1} = \mathrm{e}^{\mathbf{A}t}.

Матрица из примера 7.3.3: A=(0αβ−(α+β))\mathbf{A} = \begin{pmatrix} 0 & \alpha \\ \beta & -(\alpha +\beta ) \end{pmatrix}, где α,β>0\alpha , \beta > 0 — константы.

?
Задача 7.3.3

Объясните, почему sin⁡2A+cos⁡2A=I\sin^2 \mathbf{A} + \cos^2 \mathbf{A} = \mathbf{I} для диагонализируемой матрицы A\mathbf{A}.

?
Задача 7.3.4

Объясните, почему e0=I\mathrm{e}^{\mathbf{0}} = \mathbf{I} для любой квадратной нулевой матрицы.

?
Задача 7.3.5

Свойство спектрального отображения для диагонализируемых матриц гласит, что если f(A)f(\mathbf{A}) существует и {λ1,λ2,…,λn}\left\{ \lambda_1,\lambda_2,\ldots ,\lambda_n\right\} — собственные значения An×n\mathbf{A}_{n\times n} (с учётом кратностей), то {f(λ1),…,f(λn)}\left\{ f(\lambda_1),\ldots ,f(\lambda_n)\right\} — собственные значения f(A)f(\mathbf{A}).

?
(a)

Докажите это для диагонализируемых матриц.

(b)

Докажите это, когда бесконечный ряд f(z)=∑n=0∞cn(z−z0)nf(z) = \sum_{n=0}^\infty c_n(z-z_0)^n задаёт f(A)=∑n=0∞cn(A−z0I)nf(\mathbf{A}) = \sum_{n=0}^\infty c_n(\mathbf{A}-z_0\mathbf{I})^n.

Задача 7.3.6

Объясните, почему det⁡(eA)=etrace⁡(A)\operatorname {det}\left(\mathrm{e}^{\mathbf{A}}\right) = \mathrm{e}^{\operatorname {trace}(\mathbf{A})}.

?
Задача 7.3.7

Предположим, что для недиагонализируемых матриц Am×m\mathbf{A}_{m\times m} бесконечный ряд f(z)=∑n=0∞cn(z−z0)nf(z) = \sum_{n=0}^\infty c_n(z-z_0)^n используется для определения f(A)=∑n=0∞cn(A−z0I)nf(\mathbf{A}) = \sum_{n=0}^\infty c_n(\mathbf{A}-z_0\mathbf{I})^n. Не обращая внимания на вопросы сходимости, объясните, почему существует многочлен p(z)p(z) степени не выше m−1m-1 такой, что f(A)=p(A)f(\mathbf{A}) = p(\mathbf{A}).

?
Задача 7.3.8

Если f(A)f(\mathbf{A}) существует для диагонализируемой A\mathbf{A}, объясните, почему Af(A)=f(A)A\mathbf{A}f(\mathbf{A}) = f(\mathbf{A})\mathbf{A}. Что можно сказать, когда A\mathbf{A} не диагонализируема?

?
Задача 7.3.9

Объясните, почему eA+B=eAeB\mathrm{e}^{\mathbf{A}+\mathbf{B}} = \mathrm{e}^{\mathbf{A}}\mathrm{e}^{\mathbf{B}}, если AB=BA\mathbf{A}\mathbf{B}=\mathbf{B}\mathbf{A}. Приведите пример, показывающий, что eA+B\mathrm{e}^{\mathbf{A}+\mathbf{B}}, eAeB\mathrm{e}^{\mathbf{A}}\mathrm{e}^{\mathbf{B}} и eBeA\mathrm{e}^{\mathbf{B}}\mathrm{e}^{\mathbf{A}} могут все различаться, когда AB≠BA\mathbf{A}\mathbf{B} \ne \mathbf{B}\mathbf{A}.

Указание: Задачу 7.2.16 можно использовать для диагонализируемого случая. Для общего случая рассмотрите F(t)=e(A+B)t−eAteBt\mathbf{F}(t) = \mathrm{e}^{(\mathbf{A}+\mathbf{B})t} - \mathrm{e}^{\mathbf{A}t}\mathrm{e}^{\mathbf{B}t} и F′(t)\mathbf{F}'(t).

?
Задача 7.3.10

Покажите, что eA\mathrm{e}^{\mathbf{A}} является ортогональной матрицей всякий раз, когда A\mathbf{A} кососимметрична.

?
Задача 7.3.11

Некоторое электронное устройство состоит из набора переключающих схем, которые могут находиться либо во включённом (ON) состоянии, либо в выключенном (OFF) состоянии. Эти электронные переключатели могут менять состояние через регулярные промежутки времени, называемые тактами. Предположим, что в конце каждого такта 30%30\% переключателей, находящихся в состоянии OFF, переходят в состояние ON, а 90%90\% находящихся в состоянии ON возвращаются в состояние OFF.

?
(a)

Покажите, что устройство приближается к равновесию в том смысле, что доля переключателей в каждом состоянии со временем становится постоянной, и определите эти равновесные доли.

(b)

Независимо от начальных долей, примерно сколько тактов требуется устройству, чтобы стать практически стабильным?

Задача 7.3.12

Спектральный радиус матрицы A\mathbf{A} определяется как ρ(A)=max⁡λi∈σ(A)∣λi∣\rho (\mathbf{A}) = \max_{\lambda_i \in \sigma (\mathbf{A})} \left|\lambda_i\right|. Докажите, что если A\mathbf{A} диагонализируема, то

ρ(A)=lim⁡n→∞∥An∥1/nдля любой матричной нормы. \rho (\mathbf{A}) = \lim _{n\to \infty } \left\| \mathbf{A}^n\right\| ^{1/n} \quad \text{для любой матричной нормы.}

Замечание: Этот результат верен и для недиагонализируемых матриц, но доказательство на данном этапе более сложно.

?
Задача 7.3.13

Найдите доминирующую собственную пару для A=(723020−6−2−2)\mathbf{A} = \begin{pmatrix} 7 & 2 & 3 \\ 0 & 2 & 0 \\ -6 & -2 & -2 \end{pmatrix} методом степенных итераций.

?
Задача 7.3.14

Примените метод обратных степенных итераций, чтобы найти собственный вектор для каждого из собственных значений матрицы A\mathbf{A} из задачи 7.3.13.

?
Задача 7.3.15

Объясните, почему функция m(v→)m(\overrightarrow {v}), используемая при построении метода степенных итераций, не является непрерывной функцией, так что утверждения вроде m(x→n)→m(x→)m(\overrightarrow {x}_n) \to m(\overrightarrow {x}) при x→n→x→\overrightarrow {x}_n \to \overrightarrow {x} неверны. Тем не менее, если lim⁡n→∞x→n≠0→\lim_{n\to \infty } \overrightarrow {x}_n \ne \overrightarrow {0}, то lim⁡n→∞m(x→n)≠0\lim_{n\to \infty } m(\overrightarrow {x}_n) \ne 0.

?
Задача 7.3.16

Пусть H=(100−1−2−1021)\mathbf{H} = \begin{pmatrix} 1 & 0 & 0 \\ -1 & -2 & -1 \\ 0 & 2 & 1 \end{pmatrix}.

?
(a)

Примените «ванильную» QR-итерацию к H\mathbf{H}.

(b)

Примените QR-итерацию с одиночным сдвигом к H\mathbf{H}.

Задача 7.3.17

Покажите, что QR-итерация может не сходиться, на примере H=(001100010)\mathbf{H} = \begin{pmatrix} 0 & 0 & 1 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix}.

?
(a)

Сначала примените к H\mathbf{H} «обычную» QR-итерацию и посмотрите, что произойдёт.

(b)

Теперь попробуйте применить к H\mathbf{H} QR-итерацию с одинарным сдвигом.

(c)

Наконец, выполните для H\mathbf{H} QR-итерацию с двойным сдвигом.

§
Задача 7.4.1

Пусть An×n\mathbf{A}_{n\times n} диагонализуема, и пусть P=(x→1∣x→2∣⋯∣x→n)\mathbf{P} = \begin{pmatrix} \overrightarrow {x}_1 \mid \overrightarrow {x}_2 \mid \cdots \mid \overrightarrow {x}_n \end{pmatrix} — матрица, столбцами которой служит полный набор линейно независимых собственных векторов, соответствующих собственным значениям λi\lambda_i. Покажите, что решение уравнения u→′=Au→\overrightarrow {u}' = \mathbf{A}\overrightarrow {u}, u→(0)=c→\overrightarrow {u}(0) = \overrightarrow {c}, можно записать в виде

u→(t)=ξ1eλ1tx→1+ξ2eλ2tx→2+⋯+ξneλntx→n \overrightarrow {u}(t) = \xi _1 e^{\lambda _1 t}\overrightarrow {x}_1 + \xi _2 e^{\lambda _2 t}\overrightarrow {x}_2 + \cdots + \xi _n e^{\lambda _n t}\overrightarrow {x}_n

где коэффициенты ξi\xi_i удовлетворяют алгебраической системе Pξ→=c→\mathbf{P}\overrightarrow {\xi } = \overrightarrow {c}.

?
Задача 7.4.2

Используя только собственные значения, определите поведение решения уравнения u→′=Au→\overrightarrow {u}' = \mathbf{A}\overrightarrow {u}, u→(0)=c→\overrightarrow {u}(0) = \overrightarrow {c} при t→∞t \to \infty для каждой из следующих матриц.

?
(a)

A=(−1−20−3)\mathbf{A} = \begin{pmatrix} -1 & -2 \\ 0 & -3 \end{pmatrix}.

(b)

A=(1−203)\mathbf{A} = \begin{pmatrix} 1 & -2 \\ 0 & 3 \end{pmatrix}.

(c)

A=(1−21−1)\mathbf{A} = \begin{pmatrix} 1 & -2 \\ 1 & -1 \end{pmatrix}.

Задача 7.4.3

Рассмотрим два вида, сосуществующих в одной и той же среде, но конкурирующих за одни и те же ресурсы. Предположим, что численность каждого вида растёт пропорционально численности своего собственного вида, но убывает пропорционально численности конкурирующего вида — скажем, численность каждого вида растёт со скоростью, равной удвоенной его существующей численности, но убывает со скоростью, равной численности другой популяции. Предположим, что первоначально имеется 100 особей вида I и 200 особей вида II.

?
(a)

Определите численность каждого вида в любой последующий момент времени.

(b)

Определите, какой вид обречён на вымирание, и вычислите время до вымирания.

Задача 7.4.4

Рассмотрим два вида, выживающих в симбиотических отношениях в том смысле, что численность каждого вида убывает со скоростью, равной его собственной существующей численности, но растёт со скоростью, равной существующей численности другого вида.

?
(a)

Если первоначально имеется 200 особей вида I и 400 особей вида II, определите численность каждого вида в любой последующий момент времени.

(b)

Обсудите поведение каждого вида при t→∞t \to \infty.

§
Задача 7.5.1

Является ли A=(5+i−2i24+2i)\mathbf{A} = \begin{pmatrix} 5+i & -2i \\ 2 & 4+2i \end{pmatrix} нормальной матрицей?

?
Задача 7.5.2

Приведите примеры двух различных классов нормальных матриц, которые вещественны, но не симметричны.

?
Задача 7.5.3

Покажите, что A∈Rn×n\mathbf{A} \in \mathbb {R}^{n \times n} нормальна и имеет вещественные собственные значения тогда и только тогда, когда A\mathbf{A} симметрична.

?
Задача 7.5.4

Докажите, что собственные значения вещественной кососимметрической или косоэрмитовой матрицы обязательно являются чисто мнимыми числами (то есть кратными ii).

?
Задача 7.5.5

При попытке решить, что верно для матриц, а что нет, полезно мыслить в терминах следующих соответствий.

Эрмитовы матрицы⟷Вещественные числа (z‾=z).Косоэрмитовы матрицы⟷Мнимые числа (z‾=−z).Унитарные матрицы⟷Точки единичной окружности (z=eiθ). \begin{aligned} \text{Эрмитовы матрицы} & \longleftrightarrow \text{Вещественные числа } (\overline{z}=z). \\ \text{Косоэрмитовы матрицы} & \longleftrightarrow \text{Мнимые числа } (\overline{z}=-z). \\ \text{Унитарные матрицы} & \longleftrightarrow \text{Точки единичной окружности } (z=e^{i\theta }). \end{aligned}

Например, комплексная функция f(z)=(1−z)(1+z)−1f(z) = (1-z)(1+z)^{-1} отображает мнимую ось комплексной плоскости в точки единичной окружности, поскольку ∣f(z)∣2=1\left|f(z)\right|^2=1 всякий раз, когда z‾=−z\overline{z}=-z. Поэтому разумно предположить (как это сделал Кэли в 1846 году), что если A\mathbf{A} косоэрмитова (или вещественная кососимметрическая), то

f(A)=(I−A)(I+A)−1=(I+A)−1(I−A) f(\mathbf{A}) = (\mathbf{I}-\mathbf{A})(\mathbf{I}+\mathbf{A})^{-1} = (\mathbf{I}+\mathbf{A})^{-1}(\mathbf{I}-\mathbf{A})

унитарна (или ортогональна). Докажите, что это действительно так.

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

Выражение f(A)=(I−A)(I+A)−1=(I+A)−1(I−A)f(\mathbf{A}) = (\mathbf{I}-\mathbf{A})(\mathbf{I}+\mathbf{A})^{-1} = (\mathbf{I}+\mathbf{A})^{-1}(\mathbf{I}-\mathbf{A}) получило название преобразования Кэли.

Задача 7.5.6

Покажите на примере, что нормальная матрица может иметь полный независимый набор собственных векторов, не являющийся ортонормированным, а затем объясните, как любой полный независимый набор собственных векторов нормальной матрицы можно преобразовать в полный ортонормированный набор собственных векторов.

?
Задача 7.5.7

Постройте пример, показывающий, что обратное утверждение к «если A\mathbf{A} нормальна, то N(A−λiI)⊥N(A−λjI)N(\mathbf{A}-\lambda_i\mathbf{I}) \perp N(\mathbf{A}-\lambda_j\mathbf{I}) при λi≠λj\lambda_i \ne \lambda_j» неверно. Иными словами, покажите, что возможна ситуация, когда N(A−λiI)⊥N(A−λjI)N(\mathbf{A}-\lambda_i\mathbf{I}) \perp N(\mathbf{A}-\lambda_j\mathbf{I}) при всех i≠ji \ne j, но при этом A\mathbf{A} не является нормальной.

?
Задача 7.5.8

Объясните, почему треугольная матрица нормальна тогда и только тогда, когда она диагональна.

?
Задача 7.5.9

Используя результат упражнения 7.5.8, дайте альтернативное доказательство теоремы об унитарной диагонализации, которая утверждает, что A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} унитарно подобна диагональной матрице тогда и только тогда, когда A∗A=AA∗\mathbf{A}^{*}\mathbf{A} = \mathbf{A}\mathbf{A}^{*}.

?
Задача 7.5.10

Для нормальной матрицы A\mathbf{A} объясните, почему (λ,x→)(\lambda , \overrightarrow {x}) является собственной парой для A\mathbf{A} тогда и только тогда, когда (λ‾,x→)(\overline{\lambda }, \overrightarrow {x}) является собственной парой для A∗\mathbf{A}^{*}.

?
Задача 7.5.11

Чтобы проверить, поняли ли вы доказательство min-max-части теоремы Куранта--Фишера, постройте аналогичное доказательство для max-min-части.

Теорема Куранта--Фишера, min-max-часть: если A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} эрмитова с собственными значениями λ1≥λ2≥⋯≥λn\lambda_1 \ge \lambda_2 \ge \cdots \ge \lambda_n, то

λi=min⁡v→1,…,v→n−i∈Cnmax⁡x→⊥v→1,…,v→n−i∥x→∥2=1x→∗Ax→, \lambda _i = \min _{\overrightarrow {v}_1, \ldots , \overrightarrow {v}_{n-i} \in \mathbb {C}^{n}} \max _{\substack {\overrightarrow {x} \perp \overrightarrow {v}_1, \ldots , \overrightarrow {v}_{n-i} \\ \left\| \overrightarrow {x}\right\| _2 = 1}} \overrightarrow {x}^{*}\mathbf{A}\overrightarrow {x},

а max-min-часть утверждает аналогичную формулу

λi=max⁡dim⁡V=imin⁡y→∈V∥y→∥2=1y→∗Ay→. \lambda _i = \max _{\dim V = i} \min _{\substack {\overrightarrow {y} \in V \\ \left\| \overrightarrow {y}\right\| _2 = 1}} \overrightarrow {y}^{*}\mathbf{A}\overrightarrow {y}.
?
Задача 7.5.12

Теорема Куранта--Фишера имеет следующую альтернативную формулировку при 1<i<n1 < i < n:

λi=max⁡v→1,…,v→n−i∈Cnmin⁡x→⊥v→1,…,v→n−i∥x→∥2=1x→∗Ax→иλi=min⁡v→1,…,v→i−1∈Cnmax⁡x→⊥v→1,…,v→i−1∥x→∥2=1x→∗Ax→. \lambda _i = \max _{\overrightarrow {v}_1, \ldots , \overrightarrow {v}_{n-i} \in \mathbb {C}^{n}} \min _{\substack {\overrightarrow {x} \perp \overrightarrow {v}_1, \ldots , \overrightarrow {v}_{n-i} \\ \left\| \overrightarrow {x}\right\| _2 = 1}} \overrightarrow {x}^{*}\mathbf{A}\overrightarrow {x} \quad \text{и} \quad \lambda _i = \min _{\overrightarrow {v}_1, \ldots , \overrightarrow {v}_{i-1} \in \mathbb {C}^{n}} \max _{\substack {\overrightarrow {x} \perp \overrightarrow {v}_1, \ldots , \overrightarrow {v}_{i-1} \\ \left\| \overrightarrow {x}\right\| _2 = 1}} \overrightarrow {x}^{*}\mathbf{A}\overrightarrow {x}.

Чтобы проверить, действительно ли вы понимаете доказательство min-max-части теоремы Куранта--Фишера, адаптируйте его для доказательства приведённой выше альтернативной min-max-формулировки.

?
Задача 7.5.13
?
(a)

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

D=(eiθ10⋯00eiθ2⋯0⋮⋮⋱⋮00⋯eiθn). \mathbf{D} = \begin{pmatrix} e^{i\theta _1} & 0 & \cdots & 0 \\ 0 & e^{i\theta _2} & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & e^{i\theta _n} \end{pmatrix}.
(b)

Докажите, что всякая ортогональная матрица ортогонально подобна вещественной блочно-диагональной матрице вида

B=(±1⋱±1cos⁡θ1sin⁡θ1−sin⁡θ1cos⁡θ1⋱cos⁡θtsin⁡θt−sin⁡θtcos⁡θt). \mathbf{B} = \begin{pmatrix} \pm 1 & & & & & \\ & \ddots & & & & \\ & & \pm 1 & & & \\ & & & \cos \theta _1 & \sin \theta _1 & \\ & & & -\sin \theta _1 & \cos \theta _1 & \\ & & & & & \ddots \\ & & & & & & \cos \theta _t & \sin \theta _t \\ & & & & & & -\sin \theta _t & \cos \theta _t \end{pmatrix}.
§
Задача 7.6.1

Какие из следующих матриц являются положительно определёнными?

A=(1−1−1−151−115),B=(2068630808),C=(202062224). \mathbf{A} = \begin{pmatrix} 1 & -1 & -1 \\ -1 & 5 & 1 \\ -1 & 1 & 5 \end{pmatrix}, \quad \mathbf{B} = \begin{pmatrix} 20 & 6 & 8 \\ 6 & 3 & 0 \\ 8 & 0 & 8 \end{pmatrix}, \quad \mathbf{C} = \begin{pmatrix} 2 & 0 & 2 \\ 0 & 6 & 2 \\ 2 & 2 & 4 \end{pmatrix}.
?
Задача 7.6.2

Две массы m1m_1 и m2m_2 подвешены между тремя одинаковыми пружинами (с жёсткостью kk), как показано на рисунке 7.6.7. Каждая масса первоначально смещается из положения равновесия на некоторое горизонтальное расстояние и отпускается для свободных колебаний (предполагается, что вертикального смещения нет).

Рисунок 7.6.7Рисунок 7.6.7

?
(a)

Если xi(t)x_i(t) обозначает горизонтальное смещение массы mim_i от положения равновесия в момент времени tt, покажите, что Mx→′′=Kx→\mathbf{M}\overrightarrow {x}'' = \mathbf{K}\overrightarrow {x}, где

M=(m100m2),x→=(x1(t)x2(t)),иK=k(2−1−12). \mathbf{M} = \begin{pmatrix} m_1 & 0 \\ 0 & m_2 \end{pmatrix}, \quad \overrightarrow {x} = \begin{pmatrix} x_1(t) \\ x_2(t) \end{pmatrix}, \quad \text{и} \quad \mathbf{K} = k\begin{pmatrix} 2 & -1 \\ -1 & 2 \end{pmatrix}.

(Считайте силу, направленную влево, положительной.) Заметьте, что уравнение массы-жёсткости Mx→′′=Kx→\mathbf{M}\overrightarrow {x}'' = \mathbf{K}\overrightarrow {x} представляет собой матричную версию закона Гука F=kxF=kx, и K\mathbf{K} положительно определена.

(b)

Ищите решение в виде x→=eiθtv→\overrightarrow {x} = e^{i\theta t}\overrightarrow {v} для некоторого постоянного вектора v→\overrightarrow {v} и покажите, что это сводит задачу к решению алгебраического уравнения вида Kv→=λMv→\mathbf{K}\overrightarrow {v} = \lambda \mathbf{M}\overrightarrow {v} (при λ=−θ2\lambda = -\theta^2). Это называется обобщённой задачей о собственных значениях, поскольку при M=I\mathbf{M} = \mathbf{I} мы возвращаемся к обычной задаче о собственных значениях. Обобщённые собственные значения λ1\lambda_1 и λ2\lambda_2 являются корнями уравнения det⁡(K−λM)=0\operatorname {det}\left(\mathbf{K}-\lambda \mathbf{M}\right) = 0 --- найдите их при k=1k=1, m1=1m_1=1 и m2=2m_2=2 и опишите две моды колебаний.

(c)

Положите m1=m2=mm_1 = m_2 = m и примените технику, использованную в примере с колеблющимися бусинами, чтобы определить нормальные моды. Сравните результаты с результатами пункта (б).

Задача 7.6.3

Три массы m1m_1, m2m_2 и m3m_3 подвешены на трёх одинаковых пружинах (с жёсткостью kk), как показано ниже. Каждая масса первоначально смещается из положения равновесия на некоторое вертикальное расстояние, а затем отпускается для свободных колебаний.

Три массы, подвешенные последовательно на трёх одинаковых пружинахТри массы, подвешенные последовательно на трёх одинаковых пружинах

?
(a)

Если yi(t)y_i(t) обозначает смещение массы mim_i от положения равновесия в момент времени tt, покажите, что уравнение массы-жёсткости имеет вид My→′′=Ky→\mathbf{M}\overrightarrow {y}'' = \mathbf{K}\overrightarrow {y}, где

M=(m1000m2000m3),y→=(y1(t)y2(t)y3(t)),K=k(2−10−12−10−11) \mathbf{M} = \begin{pmatrix} m_1 & 0 & 0 \\ 0 & m_2 & 0 \\ 0 & 0 & m_3 \end{pmatrix}, \quad \overrightarrow {y} = \begin{pmatrix} y_1(t) \\ y_2(t) \\ y_3(t) \end{pmatrix}, \quad \mathbf{K} = k\begin{pmatrix} 2 & -1 & 0 \\ -1 & 2 & -1 \\ 0 & -1 & 1 \end{pmatrix}

(k33=1k_{33}=1 --- это не ошибка!).

(b)

Покажите, что K\mathbf{K} положительно определена.

(c)

Найдите нормальные моды при m1=m2=m3=mm_1=m_2=m_3=m.

Задача 7.6.4

Диагонализуя квадратичную форму 13x2+10xy+13y213x^2+10xy+13y^2, покажите, что повёрнутый график 13x2+10xy+13y2=7213x^2+10xy+13y^2=72 представляет собой эллипс в каноническом виде.

Рисунок 7.2.1 (стр. 505) показывает этот эллипс: в системе координат xyxy он расположен наклонно, с уравнением 13x2+10xy+13y2=7213x^2+10xy+13y^2=72; поворот системы координат xyxy против часовой стрелки на 45∘45^{\circ } в систему координат uvuv устраняет перекрёстный член, и уравнение эллипса становится u2/9+v2/4=1u^2/9+v^2/4=1.

?
Задача 7.6.5

Предположим, что A\mathbf{A} --- вещественная симметричная матрица. Объясните, почему знаки главных элементов (пивотов) в LDU\mathbf{L}\mathbf{D}\mathbf{U}-разложении матрицы A\mathbf{A} показывают инерцию матрицы A\mathbf{A}.

?
Задача 7.6.6

Рассмотрим квадратичную форму

f(x→)=19(−2x12+7x22+4x32+4x1x2+16x1x3+20x2x3). f(\overrightarrow {x}) = \frac{1}{9}\bigl(-2x_1^2+7x_2^2+4x_3^2+4x_1x_2+16x_1x_3+20x_2x_3\bigr).
?
(a)

Найдите симметричную матрицу A\mathbf{A}, такую что f(x→)=x→TAx→f(\overrightarrow {x}) = \overrightarrow {x}^T\mathbf{A}\overrightarrow {x}.

(b)

Диагонализуйте квадратичную форму с помощью разложения LDLT\mathbf{L}\mathbf{D}\mathbf{L}^T и определите сигнатуру (инерцию) матрицы A\mathbf{A}.

(c)

Является ли эта форма положительно определённой?

(d)

Проверьте правильность найденной выше инерции, вычислив собственные значения матрицы A\mathbf{A}.

(e)

Проверьте закон инерции Сильвестра, придумав преобразование конгруэнтности C\mathbf{C} и вычислив инерцию матрицы CTAC\mathbf{C}^T\mathbf{A}\mathbf{C}.

Задача 7.6.7

Объясните, почему каждая невырожденная матрица A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} может быть единственным образом разложена в виде A=RU\mathbf{A} = \mathbf{R}\mathbf{U}, где R\mathbf{R} — эрмитова положительно определённая матрица, а U\mathbf{U} — унитарная матрица. Это матричный аналог полярной формы комплексного числа z=reiθz = re^{i\theta }, r>0r>0, поскольку эрмитовы положительно определённые матрицы размера 1×11\times 1 — это положительные вещественные числа, а унитарные матрицы размера 1×11\times 1 — это точки на единичной окружности.

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

Указание: сначала объясните, почему R=(AA∗)1/2\mathbf{R} = (\mathbf{A}\mathbf{A}^{*})^{1/2}.

Задача 7.6.8

Объясните, почему попытки получить более точные приближения к решению задачи Дирихле, дискретизированной в виде Lu→=g→\mathbf{L}\overrightarrow {u}=\overrightarrow {g} (дискретная лапласова система из примера 7.6.2), путём использования более мелких сеток с большим числом узлов приводят ко всё более плохо обусловленной линейной системе Lu→=g→\mathbf{L}\overrightarrow {u}=\overrightarrow {g}.

Дискретизация в примере 7.6.2 даёт собственные значения дискретного лапласиана L\mathbf{L} в виде λij=4−2cos⁡(iπn+1)−2cos⁡(jπn+1)\lambda_{ij} = 4 - 2\cos \bigl(\frac{i\pi }{n+1}\bigr) - 2\cos \bigl(\frac{j\pi }{n+1}\bigr) для i,j=1,2,…,ni,j = 1, 2, \ldots , n, где nn — число внутренних узлов сетки на одной стороне.

?
Задача 7.6.9

Для заданной функции ff уравнение ∇2u=f\nabla^2 u = f называется уравнением Пуассона. Рассмотрим уравнение Пуассона на квадрате в двумерном пространстве с граничными условиями Дирихле. То есть

∂2u∂x2+∂2u∂y2=f(x,y)приu(x,y)=g(x,y) на границе. \frac{\partial ^2 u}{\partial x^2} + \frac{\partial ^2 u}{\partial y^2} = f(x,y) \quad \text{при} \quad u(x,y) = g(x,y) \text{ на границе.}

Дискретизируем задачу, наложив на квадрат регулярную сетку с n2n^2 внутренними узлами, отстоящими друг от друга на расстояние hh, как в примере 7.6.2. Пусть fij=f(xi,yj)f_{ij} = f(x_i,y_j), и определим вектор f→\overrightarrow {f} как f→=(f11,f12,…,f1n∣f21,f22,…,f2n∣⋯∣fn1,fn2,…,fnn)T\overrightarrow {f} = (f_{11}, f_{12}, \ldots , f_{1n} \mid f_{21}, f_{22}, \ldots , f_{2n} \mid \cdots \mid f_{n1}, f_{n2}, \ldots , f_{nn})^T. Покажите, что дискретизация уравнения Пуассона приводит к системе линейных уравнений вида Lu→=g→−h2f→\mathbf{L}\overrightarrow {u} = \overrightarrow {g}-h^2\overrightarrow {f}, где L\mathbf{L} — дискретный лапласиан, а u→\overrightarrow {u} и g→\overrightarrow {g} такие же, как описано в примере 7.6.2.

?
Задача 7.6.10

Как определено в упражнении 5.8.15 и обсуждено в упражнении 7.8.11, произведение Кронекера (иногда называемое тензорным произведением или прямым произведением) матриц Am×n\mathbf{A}_{m \times n} и Bp×q\mathbf{B}_{p \times q} — это матрица размера mp×nqmp \times nq

A⊗B=(a11Ba12B⋯a1nBa21Ba22B⋯a2nB⋮⋮⋱⋮am1Bam2B⋯amnB). \mathbf{A} \otimes \mathbf{B} = \begin{pmatrix} a_{11}\mathbf{B} & a_{12}\mathbf{B} & \cdots & a_{1n}\mathbf{B} \\ a_{21}\mathbf{B} & a_{22}\mathbf{B} & \cdots & a_{2n}\mathbf{B} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1}\mathbf{B} & a_{m2}\mathbf{B} & \cdots & a_{mn}\mathbf{B} \end{pmatrix}.

Проверьте, что если In\mathbf{I}_n — единичная матрица размера n×nn\times n, и если

An=(2−1−12−1⋱⋱⋱−12−1−12)n×n \mathbf{A}_n = \begin{pmatrix} 2 & -1 & & \\ -1 & 2 & -1 & \\ & \ddots & \ddots & \ddots \\ & & -1 & 2 & -1 \\ & & & -1 & 2 \end{pmatrix}_{n \times n}

— матрица конечных разностей nn-го порядка, то дискретный лапласиан задаётся формулой

Ln2×n2=(In⊗An)+(An⊗In). \mathbf{L}_{n^2 \times n^2} = (\mathbf{I}_n \otimes \mathbf{A}_n) + (\mathbf{A}_n \otimes \mathbf{I}_n).

Таким образом, мы получаем изящную матричную связь между конечно-разностными приближениями одномерного и двумерного лапласианов. Эта формула ведёт к простому альтернативному выводу формулы (7.6.8) --- см. упражнение 7.8.12. Как можно догадаться, дискретный трёхмерный лапласиан имеет вид

Ln3×n3=(In⊗In⊗An)+(In⊗An⊗In)+(An⊗In⊗In). \mathbf{L}_{n^3 \times n^3} = (\mathbf{I}_n \otimes \mathbf{I}_n \otimes \mathbf{A}_n) + (\mathbf{I}_n \otimes \mathbf{A}_n \otimes \mathbf{I}_n) + (\mathbf{A}_n \otimes \mathbf{I}_n \otimes \mathbf{I}_n).
?
§
Задача 7.7.1

Может ли индекс нильпотентной матрицы размера n×nn\times n превышать nn?

?
Задача 7.7.2

Определите все возможные жордановы формы N\mathbf{N} для нильпотентной матрицы размера 4×44\times 4.

?
Задача 7.7.3

Объясните, почему число блоков размера i×ii\times i или больше в жордановой форме для нильпотентной матрицы задаётся выражением rk⁡(Li−1)−rk⁡(Li)\operatorname {rk}\left(\mathbf{L}^{i-1}\right) - \operatorname {rk}\left(\mathbf{L}^i\right).

?
Задача 7.7.4

Для нильпотентной матрицы L\mathbf{L} индекса kk пусть Mi=R(Li)∩N(L)M_i = R(\mathbf{L}^i) \cap N(\mathbf{L}). Докажите, что Mi⊆Mi−1M_i \subseteq M_{i-1} для каждого i=0,1,…,ki=0,1,\ldots ,k.

?
Задача 7.7.5

Докажите, что R(Lk−1)∩N(L)=R(Lk−1)R(\mathbf{L}^{k-1}) \cap N(\mathbf{L}) = R(\mathbf{L}^{k-1}) для всех нильпотентных матриц L\mathbf{L} индекса k>1k>1. Иными словами, докажите, что Mk−1=R(Lk−1)M_{k-1} = R(\mathbf{L}^{k-1}).

?
Задача 7.7.6

Пусть L\mathbf{L} — нильпотентная матрица индекса k>1k>1. Докажите, что если столбцы B\mathbf{B} образуют базис для R(Li)R(\mathbf{L}^i) при i≤k−1i \le k-1, и если {v→1,v→2,…,v→s}\left\{ \overrightarrow {v}_1, \overrightarrow {v}_2, \ldots , \overrightarrow {v}_s\right\} — базис для N(LB)N(\mathbf{L}\mathbf{B}), то {Bv→1,Bv→2,…,Bv→s}\left\{ \mathbf{B}\overrightarrow {v}_1, \mathbf{B}\overrightarrow {v}_2, \ldots , \mathbf{B}\overrightarrow {v}_s\right\} — базис для MiM_i.

?
Задача 7.7.7

Найдите P\mathbf{P} и N\mathbf{N} такие, что P−1LP=N\mathbf{P}^{-1}\mathbf{L}\mathbf{P} = \mathbf{N} имеет жорданову форму, где

L=(3321−2−1−1−11−101−5−4−3−2). \mathbf{L} = \begin{pmatrix} 3 & 3 & 2 & 1 \\ -2 & -1 & -1 & -1 \\ 1 & -1 & 0 & 1 \\ -5 & -4 & -3 & -2 \end{pmatrix}.
?
Задача 7.7.8

Определите жорданову форму для следующей нильпотентной матрицы 8×88\times 8.

L=(41301574613−54−39−19−9−6−8−2−496212101−6−5−3−21−100−32−24−13−6−2−5−1−2−10−7−20−303−2−4−3−2−10−1−101712632321). \mathbf{L} = \begin{pmatrix} 41 & 30 & 15 & 7 & 4 & 6 & 1 & 3 \\ -54 & -39 & -19 & -9 & -6 & -8 & -2 & -4 \\ 9 & 6 & 2 & 1 & 2 & 1 & 0 & 1 \\ -6 & -5 & -3 & -2 & 1 & -1 & 0 & 0 \\ -32 & -24 & -13 & -6 & -2 & -5 & -1 & -2 \\ -10 & -7 & -2 & 0 & -3 & 0 & 3 & -2 \\ -4 & -3 & -2 & -1 & 0 & -1 & -1 & 0 \\ 17 & 12 & 6 & 3 & 2 & 3 & 2 & 1 \end{pmatrix}.
?
Задача 7.7.9

Докажите, что если N\mathbf{N} — жорданова форма для нильпотентной матрицы L\mathbf{L}, как описано ниже, то для любого набора ненулевых скаляров {ζ1,ζ2,…,ζt}\left\{ \zeta_1, \zeta_2, \ldots , \zeta_t\right\} матрица L\mathbf{L} подобна матрице N~\tilde{\mathbf{N}} вида

N~=(ζ1N10⋯00ζ2N2⋯0⋮⋮⋱⋮00⋯ζtNt). \tilde{\mathbf{N}} = \begin{pmatrix} \zeta _1\mathbf{N}_1 & 0 & \cdots & 0 \\ 0 & \zeta _2\mathbf{N}_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \zeta _t\mathbf{N}_t \end{pmatrix}.

Иными словами, единицы на наддиагонали блоков Ni\mathbf{N}_i являются искусственными, поскольку на наддиагональ любого Ni\mathbf{N}_i можно принудительно поместить любое ненулевое значение. Важным в «жордановой структуре» матрицы L\mathbf{L} является число и размеры нильпотентных жордановых блоков (или цепочек), а не значения, стоящие на наддиагоналях этих блоков.

Жорданова форма для нильпотентной матрицы: всякая нильпотентная матрица Ln×n\mathbf{L}_{n\times n} индекса kk подобна блочно-диагональной матрице P−1LP=N=diag⁡(N1,N2,…,Nt)\mathbf{P}^{-1}\mathbf{L}\mathbf{P} = \mathbf{N} = \operatorname {diag}(\mathbf{N}_1, \mathbf{N}_2, \ldots , \mathbf{N}_t), в которой каждый блок Nj\mathbf{N}_j — нильпотентная матрица с единицами на наддиагонали и нулями в остальных местах.

?
§
Задача 7.8.1

Найдите жорданову форму следующей матрицы, различные собственные значения которой σ(A)={0,−1,1}\sigma (\mathbf{A}) = \left\{ 0,-1,1\right\}. Не пугайтесь размера A\mathbf{A}.

A=(−4−5−31−201−2473−130−120−1000000−112−420−31−8−14−51−601−4474−33−1−342−2−25−304−167302003). \mathbf{A} = \begin{pmatrix} -4 & -5 & -3 & 1 & -2 & 0 & 1 & -2 \\ 4 & 7 & 3 & -1 & 3 & 0 & -1 & 2 \\ 0 & -1 & 0 & 0 & 0 & 0 & 0 & 0 \\ -1 & 1 & 2 & -4 & 2 & 0 & -3 & 1 \\ -8 & -14 & -5 & 1 & -6 & 0 & 1 & -4 \\ 4 & 7 & 4 & -3 & 3 & -1 & -3 & 4 \\ 2 & -2 & -2 & 5 & -3 & 0 & 4 & -1 \\ 6 & 7 & 3 & 0 & 2 & 0 & 0 & 3 \end{pmatrix}.
?
Задача 7.8.2

Для матрицы A=(301−41−2−40−1)\mathbf{A} = \begin{pmatrix} 3 & 0 & 1 \\ -4 & 1 & -2 \\ -4 & 0 & -1 \end{pmatrix}, используя стандартную процедуру построения жорданова базиса, постройте невырожденную матрицу P\mathbf{P} такую, что P−1AP=J\mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \mathbf{J} имеет жорданову форму.

?
Задача 7.8.3

Объясните, почему index⁡(λ)≤alg mult⁡(λ)\operatorname {index}(\lambda ) \le \operatorname {alg\ mult}(\lambda ) для каждого λ∈σ(An×n)\lambda \in \sigma (\mathbf{A}_{n\times n}).

?
Задача 7.8.4

Объясните, почему index⁡(λ)=1\operatorname {index}(\lambda ) = 1 тогда и только тогда, когда λ\lambda — полупростое собственное значение.

?
Задача 7.8.5

Докажите, что каждая квадратная матрица подобна своей транспонированной.

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

Указание: рассмотрите «матрицу обращения» R=(11\iddots1)\mathbf{R} = \begin{pmatrix} & & & 1 \\ & & 1 & \\ & \iddots & & \\ 1 & & & \end{pmatrix}, полученную обращением порядка строк (или столбцов) единичной матрицы I\mathbf{I}.

Задача 7.8.6

Докажите теорему Кэли--Гамильтона с помощью жордановой формы; то есть докажите, что каждая матрица A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} удовлетворяет своему собственному характеристическому уравнению.

?
Задача 7.8.7

Докажите, что если λ\lambda — собственное значение матрицы A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} такое, что index⁡(λ)=k\operatorname {index}(\lambda ) = k и alg mult⁡A(λ)=m\operatorname {alg\ mult}_{\mathbf{A}}(\lambda ) = m, то dim⁡N((A−λI)k)=m\dim N\bigl((\mathbf{A}-\lambda \mathbf{I})^k\bigr) = m. Верно ли также, что dim⁡N((A−λI)m)=m\dim N\bigl((\mathbf{A}-\lambda \mathbf{I})^m\bigr) = m?

?
Задача 7.8.8

Пусть λj\lambda_j — собственное значение матрицы A\mathbf{A} с index⁡(λj)=kj\operatorname {index}(\lambda_j) = k_j. Докажите, что если Mi(λj)=R((A−λjI)i)∩N(A−λjI)M_i(\lambda_j) = R\bigl((\mathbf{A}-\lambda_j\mathbf{I})^i\bigr) \cap N(\mathbf{A}-\lambda_j\mathbf{I}), то

{0→}=Mkj(λj)⊆Mkj−1(λj)⊆⋯⊆M0(λj)=N(A−λjI). \left\{ \overrightarrow {0}\right\} = M_{k_j}(\lambda _j) \subseteq M_{k_j-1}(\lambda _j) \subseteq \cdots \subseteq M_0(\lambda _j) = N(\mathbf{A}-\lambda _j\mathbf{I}).
?
Задача 7.8.9

Объясните, почему (A−λjI)ix→=b→(λj)(\mathbf{A}-\lambda_j\mathbf{I})^i\overrightarrow {x} = \overrightarrow {b}(\lambda_j) обязательно должна быть совместной системой всякий раз, когда λj∈σ(A)\lambda_j \in \sigma (\mathbf{A}) и b→(λj)∈Si(λj)\overrightarrow {b}(\lambda_j) \in S_i(\lambda_j), где Si(λj)S_i(\lambda_j) — расширяющее множество, используемое при построении базиса для N(A−λjI)N(\mathbf{A}-\lambda_j\mathbf{I}) из базиса для Mi(λj)=R((A−λjI)i)∩N(A−λjI)M_i(\lambda_j) = R\bigl((\mathbf{A}-\lambda_j\mathbf{I})^i\bigr) \cap N(\mathbf{A}-\lambda_j\mathbf{I}), так что каждый b→(λj)∈Si(λj)\overrightarrow {b}(\lambda_j) \in S_i(\lambda_j) лежит в Mi(λj)M_i(\lambda_j).

?
Задача 7.8.10

Распространяется ли результат упражнения 7.7.5 на ненильпотентные матрицы? То есть, если λ∈σ(A)\lambda \in \sigma (\mathbf{A}) с index⁡(λ)=k>1\operatorname {index}(\lambda ) = k > 1, верно ли, что Mk−1=R((A−λI)k−1)M_{k-1} = R\bigl((\mathbf{A}-\lambda \mathbf{I})^{k-1}\bigr)?

?
Задача 7.8.11

Как было определено в упражнении 5.8.15 и упомянуто в упражнении 7.6.10, произведение Кронекера (иногда называемое тензорным произведением или прямым произведением) матриц Am×n\mathbf{A}_{m\times n} и Bp×q\mathbf{B}_{p\times q} — это mp×nqmp\times nq матрица

A⊗B=(a11Ba12B⋯a1nBa21Ba22B⋯a2nB⋮⋮⋱⋮am1Bam2B⋯amnB). \mathbf{A} \otimes \mathbf{B} = \begin{pmatrix} a_{11}\mathbf{B} & a_{12}\mathbf{B} & \cdots & a_{1n}\mathbf{B} \\ a_{21}\mathbf{B} & a_{22}\mathbf{B} & \cdots & a_{2n}\mathbf{B} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1}\mathbf{B} & a_{m2}\mathbf{B} & \cdots & a_{mn}\mathbf{B} \end{pmatrix}.
?
(a)

Предполагая согласованность размеров, установите следующие свойства.

  • A⊗(B⊗C)=(A⊗B)⊗C\mathbf{A} \otimes (\mathbf{B} \otimes \mathbf{C}) = (\mathbf{A} \otimes \mathbf{B}) \otimes \mathbf{C}.

  • A⊗(B+C)=(A⊗B)+(A⊗C)\mathbf{A} \otimes (\mathbf{B}+\mathbf{C}) = (\mathbf{A}\otimes \mathbf{B}) + (\mathbf{A}\otimes \mathbf{C}).

  • (A+B)⊗C=(A⊗C)+(B⊗C)(\mathbf{A}+\mathbf{B}) \otimes \mathbf{C} = (\mathbf{A}\otimes \mathbf{C}) + (\mathbf{B}\otimes \mathbf{C}).

  • (A1⊗B1)(A2⊗B2)⋯(Ak⊗Bk)=(A1⋯Ak)⊗(B1⋯Bk)(\mathbf{A}_1\otimes \mathbf{B}_1)(\mathbf{A}_2\otimes \mathbf{B}_2)\cdots (\mathbf{A}_k\otimes \mathbf{B}_k) = (\mathbf{A}_1\cdots \mathbf{A}_k)\otimes (\mathbf{B}_1\cdots \mathbf{B}_k).

  • (A⊗B)∗=A∗⊗B∗(\mathbf{A}\otimes \mathbf{B})^{*} = \mathbf{A}^{*}\otimes \mathbf{B}^{*}.

  • rk⁡(A⊗B)=rk⁡(A)rk⁡(B)\operatorname {rk}\left(\mathbf{A}\otimes \mathbf{B}\right) = \operatorname {rk}\left(\mathbf{A}\right)\operatorname {rk}\left(\mathbf{B}\right).

  • Для оставшихся свойств предположите, что A\mathbf{A} — матрица размера m×mm\times m, а B\mathbf{B} — размера n×nn\times n.

  • trace⁡(A⊗B)=trace⁡(A)trace⁡(B)\operatorname {trace}(\mathbf{A}\otimes \mathbf{B}) = \operatorname {trace}(\mathbf{A})\operatorname {trace}(\mathbf{B}).

  • (A⊗In)(Im⊗B)=A⊗B=(Im⊗B)(A⊗In)(\mathbf{A}\otimes \mathbf{I}_n)(\mathbf{I}_m\otimes \mathbf{B}) = \mathbf{A}\otimes \mathbf{B} = (\mathbf{I}_m\otimes \mathbf{B})(\mathbf{A}\otimes \mathbf{I}_n).

  • det⁡(A⊗B)=det⁡(A)mdet⁡(B)n\operatorname {det}\left(\mathbf{A}\otimes \mathbf{B}\right) = \operatorname {det}\left(\mathbf{A}\right)^m\operatorname {det}\left(\mathbf{B}\right)^n.

  • (A⊗B)−1=A−1⊗B−1(\mathbf{A}\otimes \mathbf{B})^{-1} = \mathbf{A}^{-1}\otimes \mathbf{B}^{-1}.

(b)

Пусть собственные значения Am×m\mathbf{A}_{m\times m} обозначены через λi\lambda_i, а собственные значения Bn×n\mathbf{B}_{n\times n} — через μj\mu_j. Докажите следующее.

  • Собственные значения A⊗B\mathbf{A}\otimes \mathbf{B} — это mnmn чисел {λiμj}i=1,j=1m,n\left\{ \lambda_i\mu_j\right\}_{i=1,j=1}^{m,n}.

  • Собственные значения (A⊗In)+(Im⊗B)(\mathbf{A}\otimes \mathbf{I}_n) + (\mathbf{I}_m\otimes \mathbf{B}) — это {λi+μj}i=1,j=1m,n\left\{ \lambda_i+\mu_j\right\}_{i=1,j=1}^{m,n}.

Задача 7.8.12

Используя пункт (б) упражнения 7.8.11 вместе с результатом упражнения 7.6.10, постройте альтернативный вывод формулы (7.6.8). А именно, покажите, что n2n^2 собственных значений дискретного лапласиана Ln2×n2\mathbf{L}_{n^2 \times n^2}, описанного в примере 7.6.2, задаются формулой

λij=4(sin⁡2(iπ2(n+1))+sin⁡2(jπ2(n+1))),i,j=1,2,…,n. \lambda _{ij} = 4\left(\sin ^2\left(\frac{i\pi }{2(n+1)}\right) + \sin ^2\left(\frac{j\pi }{2(n+1)}\right)\right), \quad i,j=1,2,\ldots ,n.
?
Примечание.
?

Указание: вспомните упражнение 7.2.18.

Задача 7.8.13

Определите собственные значения трёхмерного дискретного лапласиана, используя формулу из упражнения 7.6.10, которая утверждает, что

Ln3×n3=(In⊗In⊗An)+(In⊗An⊗In)+(An⊗In⊗In). \mathbf{L}_{n^3\times n^3} = (\mathbf{I}_n\otimes \mathbf{I}_n\otimes \mathbf{A}_n) + (\mathbf{I}_n\otimes \mathbf{A}_n\otimes \mathbf{I}_n) + (\mathbf{A}_n\otimes \mathbf{I}_n\otimes \mathbf{I}_n).
?
§
Задача 7.9.1

Озеро ii в замкнутой системе из трёх озёр равного объёма VV первоначально содержит cic_i фунтов загрязняющего вещества. Если вода в системе циркулирует со скоростями (галлонов/сек), указанными на рисунке 7.9.2, найдите количество загрязняющего вещества в каждом озере в момент времени t>0t>0 (предполагая непрерывное перемешивание), а затем определите загрязнённость каждого озера в долгосрочной перспективе.

Рисунок 7.9.2Рисунок 7.9.2

?
Задача 7.9.2

Предположим, что A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} имеет собственные значения λi\lambda_i с index⁡(λi)=ki\operatorname {index}(\lambda_i)=k_i. Объясните, почему ii-й спектральный проектор задаётся формулой

Gi=fi(A),гдеfi(z)={1при z=λi,0иначе. \mathbf{G}_i = f_i(\mathbf{A}), \quad \text{где} \quad f_i(z) = \begin{cases} 1 & \text{при } z=\lambda _i, \\ 0 & \text{иначе.} \end{cases}

Спектральное разложение f(A)f(\mathbf{A}): для A\mathbf{A} с σ(A)={λ1,…,λs}\sigma (\mathbf{A}) = \left\{ \lambda_1, \ldots , \lambda_s\right\}, ki=index⁡(λi)k_i = \operatorname {index}(\lambda_i), и ff такой, что f(λi),f′(λi),…,f(ki−1)(λi)f(\lambda_i), f'(\lambda_i), \ldots , f^{(k_i-1)}(\lambda_i) существуют для каждого λi\lambda_i,

f(A)=∑i=1s∑j=0ki−1f(j)(λi)j!(A−λiI)jGi. f(\mathbf{A}) = \sum _{i=1}^s \sum _{j=0}^{k_i-1} \frac{f^{(j)}(\lambda _i)}{j!}(\mathbf{A}-\lambda _i\mathbf{I})^j\mathbf{G}_i.
?
Задача 7.9.3

Объясните, почему каждый спектральный проектор Gi\mathbf{G}_i можно выразить в виде многочлена от A\mathbf{A}.

?
Задача 7.9.4

Если σ(An×n)={λ1,λ2,…,λs}\sigma (\mathbf{A}_{n\times n}) = \left\{ \lambda_1, \lambda_2, \ldots , \lambda_s\right\} с ki=index⁡(λi)k_i = \operatorname {index}(\lambda_i), объясните, почему

Ak=∑i=1s∑j=0ki−1(kj)λik−j(A−λiI)jGi. \mathbf{A}^k = \sum _{i=1}^s \sum _{j=0}^{k_i-1} \binom {k}{j}\lambda _i^{k-j}(\mathbf{A}-\lambda _i\mathbf{I})^j\mathbf{G}_i.

Спектральное разложение f(A)f(\mathbf{A}): f(A)=∑i=1s∑j=0ki−1f(j)(λi)j!(A−λiI)jGif(\mathbf{A}) = \sum_{i=1}^s \sum_{j=0}^{k_i-1} \dfrac {f^{(j)}(\lambda_i)}{j!}(\mathbf{A}-\lambda_i\mathbf{I})^j\mathbf{G}_i.

?
Задача 7.9.5

При соглашении (kj)=0\binom {k}{j} = 0 для j>kj>k, объясните, почему для жорданова блока размера k×kk\times k с собственным значением λ\lambda

(λ1⋱⋱⋱1λ) ⁣ ⁣k=(λk(k1)λk−1(k2)λk−2⋯(km−1)λk−m+1λk(k1)λk−1⋱⋮⋱⋱(k2)λk−2λk(k1)λk−1λk)m×m. \begin{pmatrix} \lambda & 1 & & \\ & \ddots & \ddots & \\ & & \ddots & 1 \\ & & & \lambda \end{pmatrix}^{\! \! k} = \begin{pmatrix} \lambda ^k & \binom {k}{1}\lambda ^{k-1} & \binom {k}{2}\lambda ^{k-2} & \cdots & \binom {k}{m-1}\lambda ^{k-m+1} \\ & \lambda ^k & \binom {k}{1}\lambda ^{k-1} & \ddots & \vdots \\ & & \ddots & \ddots & \binom {k}{2}\lambda ^{k-2} \\ & & & \lambda ^k & \binom {k}{1}\lambda ^{k-1} \\ & & & & \lambda ^k \end{pmatrix}_{m\times m}.

Функции жордановых блоков: для жорданова блока J⋆\mathbf{J}_{\star } размера k×kk\times k с собственным значением λ\lambda, и для ff такой, что f(λ),f′(λ),…,f(k−1)(λ)f(\lambda ), f'(\lambda ), \ldots , f^{(k-1)}(\lambda ) существуют,

f(J⋆)=(f(λ)f′(λ)f”(λ)2!⋯f(k−1)(λ)(k−1)!f(λ)f′(λ)⋱⋮⋱⋱f”(λ)2!f(λ)f′(λ)f(λ)). f(\mathbf{J}_{\star }) = \begin{pmatrix} f(\lambda ) & f'(\lambda ) & \dfrac {f”(\lambda )}{2!} & \cdots & \dfrac {f^{(k-1)}(\lambda )}{(k-1)!} \\ & f(\lambda ) & f'(\lambda ) & \ddots & \vdots \\ & & \ddots & \ddots & \dfrac {f”(\lambda )}{2!} \\ & & & f(\lambda ) & f'(\lambda ) \\ & & & & f(\lambda ) \end{pmatrix}.
?
Задача 7.9.6

Найдите eAe^{\mathbf{A}} для A=(628−22−2002)\mathbf{A} = \begin{pmatrix} 6 & 2 & 8 \\ -2 & 2 & -2 \\ 0 & 0 & 2 \end{pmatrix}.

?
Задача 7.9.7

Для f(z)=z−14f(z) = \sqrt[4]{z-1} найдите f(A)f(\mathbf{A}), если A=(−3−8−95119−1−21)\mathbf{A} = \begin{pmatrix} -3 & -8 & -9 \\ 5 & 11 & 9 \\ -1 & -2 & 1 \end{pmatrix}.

?
Задача 7.9.8
?
(a)

Объясните, почему каждая невырожденная A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} имеет квадратный корень.

(b)

Дайте необходимые и достаточные условия существования A\sqrt{\mathbf{A}}, когда A\mathbf{A} вырождена.

Задача 7.9.9

Докажите, что если (λ,x→)(\lambda , \overrightarrow {x}) — собственная пара для A\mathbf{A}, то (f(λ),x→)(f(\lambda ), \overrightarrow {x}) — собственная пара для f(A)f(\mathbf{A}) всякий раз, когда f(A)f(\mathbf{A}) существует. Следует ли отсюда также, что alg mult⁡A(λ)=alg mult⁡f(A)(f(λ))\operatorname {alg\ mult}_{\mathbf{A}}(\lambda ) = \operatorname {alg\ mult}_{f(\mathbf{A})}(f(\lambda ))?

?
Задача 7.9.10

Пусть ff определена на A\mathbf{A}, и пусть λ∈σ(A)\lambda \in \sigma (\mathbf{A}). Приведите пример или объяснение того, почему следующие утверждения не обязательно верны.

?
(a)

f(A)f(\mathbf{A}) подобна A\mathbf{A}.

(b)

geo mult⁡A(λ)=geo mult⁡f(A)(f(λ))\operatorname {geo\ mult}_{\mathbf{A}}(\lambda ) = \operatorname {geo\ mult}_{f(\mathbf{A})}(f(\lambda )).

(c)

index⁡A(λ)=index⁡f(A)(f(λ))\operatorname {index}_{\mathbf{A}}(\lambda ) = \operatorname {index}_{f(\mathbf{A})}(f(\lambda )).

Задача 7.9.11

Объясните, почему Af(A)=f(A)A\mathbf{A}f(\mathbf{A}) = f(\mathbf{A})\mathbf{A} всякий раз, когда f(A)f(\mathbf{A}) существует.

?
Задача 7.9.12

Объясните, почему функция ff определена на A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} тогда и только тогда, когда ff определена на AT\mathbf{A}^T, а затем докажите, что f(AT)=f(A)Tf(\mathbf{A}^T) = f(\mathbf{A})^T. Почему нельзя использовать ( ⋅ )∗(\, \cdot \, )^{*} вместо ( ⋅ )T(\, \cdot \, )^T?

?
Задача 7.9.13

Используя приём записи h(z)=p(f1(z),f2(z))h(z) = p(f_1(z), f_2(z)) для полиномиального тождества pp, тождественно обращающегося в нуль, так что h(A)=p(f1(A),f2(A))=0h(\mathbf{A}) = p(f_1(\mathbf{A}), f_2(\mathbf{A})) = \mathbf{0}, установите следующие тождества.

?
(a)

eAe−A=Ie^{\mathbf{A}}e^{-\mathbf{A}} = \mathbf{I} для всех A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n}.

(b)

eαA=(eA)αe^{\alpha \mathbf{A}} = \bigl(e^{\mathbf{A}}\bigr)^{\alpha } для всех α∈C\alpha \in \mathbb {C} и A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n}.

(c)

eiA=cos⁡A+isin⁡Ae^{i\mathbf{A}} = \cos \mathbf{A} + i\sin \mathbf{A} для всех A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n}.

Задача 7.9.14
?
(a)

Покажите, что если AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}, то eA+B=eAeBe^{\mathbf{A}+\mathbf{B}} = e^{\mathbf{A}}e^{\mathbf{B}}.

(b)

Приведите пример, показывающий, что в общем случае eA+B≠eAeBe^{\mathbf{A}+\mathbf{B}} \ne e^{\mathbf{A}}e^{\mathbf{B}}.

Задача 7.9.15

Найдите интерполяционный многочлен Эрмита p(z)p(z), т.е. многочлен наименьшей степени, совпадающий со значениями и производными ff в каждом собственном значении вплоть до порядка на единицу меньше его индекса, такой что p(A)=eAp(\mathbf{A}) = e^{\mathbf{A}} для A=(321−3−2−1−3−2−1)\mathbf{A} = \begin{pmatrix} 3 & 2 & 1 \\ -3 & -2 & -1 \\ -3 & -2 & -1 \end{pmatrix}.

?
Задача 7.9.16

Теорема Кэли--Гамильтона утверждает, что каждая A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} удовлетворяет своему собственному характеристическому уравнению, и это гарантирует, что An+j\mathbf{A}^{n+j} (j=0,1,2,…j=0,1,2,\ldots) может быть выражена как многочлен от A\mathbf{A} степени не выше n−1n-1. Поскольку f(A)f(\mathbf{A}) всегда является многочленом от A\mathbf{A}, теорема Кэли--Гамильтона гарантирует, что f(A)f(\mathbf{A}) может быть выражена как многочлен от A\mathbf{A} степени не выше n−1n-1. Такой многочлен можно определить всякий раз, когда f(j)(λi)f^{(j)}(\lambda_i), j=0,1,…,ai−1j=0,1,\ldots ,a_i-1 существует для каждого λi∈σ(A)\lambda_i \in \sigma (\mathbf{A}), где ai=alg mult⁡(λi)a_i = \operatorname {alg\ mult}(\lambda_i). Стратегия та же, что использовалась для нахождения интерполяционного многочлена Эрмита, согласованного с индексом собственного значения, за исключением того, что вместо kik_i используется aia_i. Если мы можем найти многочлен p(z)=α0+α1z+⋯+αn−1zn−1p(z) = \alpha_0+\alpha_1 z+\cdots +\alpha_{n-1}z^{n-1} такой, что для каждого λi∈σ(A)\lambda_i \in \sigma (\mathbf{A})

p(λi)=f(λi),p′(λi)=f′(λi),…,p(ai−1)(λi)=f(ai−1)(λi), p(\lambda _i) = f(\lambda _i), \quad p'(\lambda _i) = f'(\lambda _i), \quad \ldots , \quad p^{(a_i-1)}(\lambda _i) = f^{(a_i-1)}(\lambda _i),

то p(A)=f(A)p(\mathbf{A}) = f(\mathbf{A}). Почему? Эти уравнения представляют собой линейную систему n×nn\times n с αi\alpha_i в качестве неизвестных, и по той же причине, что и раньше, решение всегда возможно.

?
(a)

Какие преимущества и недостатки имеет этот подход по сравнению с подходом на основе индекса из упражнения 7.9.15?

(b)

Используйте этот метод, чтобы найти многочлен p(z)p(z) такой, что p(A)=eAp(\mathbf{A}) = e^{\mathbf{A}} для A=(321−3−2−1−3−2−1)\mathbf{A} = \begin{pmatrix} 3 & 2 & 1 \\ -3 & -2 & -1 \\ -3 & -2 & -1 \end{pmatrix}. Сравните с упражнением 7.9.15.

Задача 7.9.17

Покажите, что если ff — функция, определённая на

A=(αβγ0αβ00α)=αI+βN+γN2,гдеN=(010001000), \mathbf{A} = \begin{pmatrix} \alpha & \beta & \gamma \\ 0 & \alpha & \beta \\ 0 & 0 & \alpha \end{pmatrix} = \alpha \mathbf{I}+\beta \mathbf{N}+\gamma \mathbf{N}^2, \quad \text{где} \quad \mathbf{N} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix},

то f(A)=f(α)I+βf′(α)N+(γf′(α)+β2f”(α)2!)N2f(\mathbf{A}) = f(\alpha )\mathbf{I} + \beta f'(\alpha )\mathbf{N} + \left(\gamma f'(\alpha ) + \dfrac {\beta^2 f”(\alpha )}{2!}\right)\mathbf{N}^2.

?
Задача 7.9.18

Если h(z)=f(g(z))h(z) = f(g(z)), где ff и gg — функции, такие что g(A)g(\mathbf{A}) и f(g(A))f(g(\mathbf{A})) каждая существует, то h(A)=f(g(A))h(\mathbf{A}) = f(g(\mathbf{A})). Однако доказывать это простой заменой ``zz на A\mathbf{A}'' недопустимо. Один из способов доказать, что h(A)=f(g(A))h(\mathbf{A}) = f(g(\mathbf{A})), — показать, что h(J⋆)=f(g(J⋆))h(\mathbf{J}_{\star }) = f(g(\mathbf{J}_{\star })) для типичного жорданова блока, а затем воспользоваться тем, что f(A)=Pf(J)P−1f(\mathbf{A}) = \mathbf{P}f(\mathbf{J})\mathbf{P}^{-1} строится поблочно из жордановых блоков J\mathbf{J}. Проделайте это для жорданова блока 3×33\times 3 --- обобщение на блоки k×kk\times k аналогично. То есть пусть h(z)=f(g(z))h(z) = f(g(z)), и с помощью упражнения 7.9.17 докажите, что если g(J⋆)g(\mathbf{J}_{\star }) и f(g(J⋆))f(g(\mathbf{J}_{\star })) каждая существует, то

h(J⋆)=f(g(J⋆))для J⋆=(λ100λ100λ). h(\mathbf{J}_{\star }) = f\bigl(g(\mathbf{J}_{\star })\bigr) \quad \text{для } \mathbf{J}_{\star } = \begin{pmatrix} \lambda & 1 & 0 \\ 0 & \lambda & 1 \\ 0 & 0 & \lambda \end{pmatrix}.
?
Задача 7.9.19

Докажите, что если Γi\Gamma_i — простой замкнутый контур, охватывающий λi∈σ(A)\lambda_i \in \sigma (\mathbf{A}), но не охватывающий остальные собственные значения A\mathbf{A}, то ii-й спектральный проектор равен

Gi=12πi∮Γi(ξI−A)−1 dξ=12πi∮ΓiR(ξ) dξ. \mathbf{G}_i = \frac{1}{2\pi i}\oint _{\Gamma _i} (\xi \mathbf{I}-\mathbf{A})^{-1}\, d\xi = \frac{1}{2\pi i}\oint _{\Gamma _i} \mathbf{R}(\xi )\, d\xi .
?
Задача 7.9.20

Для f(z)=z−1f(z) = z^{-1} проверьте, что f(A)=A−1f(\mathbf{A}) = \mathbf{A}^{-1} для любой невырожденной A\mathbf{A}.

?
Задача 7.9.21

Если Γ\Gamma — простой замкнутый контур, охватывающий все собственные значения невырожденной матрицы A\mathbf{A}, чему равно значение 12πi∮Γξ−1(ξI−A)−1 dξ\dfrac {1}{2\pi i}\oint_{\Gamma } \xi^{-1}(\xi \mathbf{I}-\mathbf{A})^{-1}\, d\xi?

?
Задача 7.9.22

Функция обращения f(z)=z−1f(z) = z^{-1} не определена на вырожденных матрицах, но обобщённая функция обращения

g(z)={z−1если z≠0,0если z=0, g(z) = \begin{cases} z^{-1} & \text{если } z \ne 0, \\ 0 & \text{если } z = 0, \end{cases}

определена на всех квадратных матрицах. Из упражнения 7.9.20 ясно, что если A\mathbf{A} невырожденна, то g(A)=A−1g(\mathbf{A}) = \mathbf{A}^{-1}, так что g(A)g(\mathbf{A}) — естественный способ распространить понятие обращения на вырожденные матрицы. Объясните, почему g(A)=ADg(\mathbf{A}) = \mathbf{A}^D является дрейзиновской обратной матрицей A\mathbf{A}, а не обязательно псевдообратной Мура--Пенроуза A†\mathbf{A}^{\dagger }.

?
Задача 7.9.23

Пусть A\mathbf{A} — вырожденная матрица, и пусть Γ\Gamma — простой замкнутый контур, содержащий все собственные значения A\mathbf{A}, кроме λ1=0\lambda_1=0, которое не лежит ни внутри, ни на Γ\Gamma. Докажите, что

12πi∮Γξ−1(ξI−A)−1 dξ=AD \frac{1}{2\pi i}\oint _{\Gamma } \xi ^{-1}(\xi \mathbf{I}-\mathbf{A})^{-1}\, d\xi = \mathbf{A}^D

есть дрейзиновская обратная матрица A\mathbf{A}.

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

Подсказка: теорема Коши--Гурса утверждает, что если функция ff аналитична во всех точках внутри и на простом замкнутом контуре Γ\Gamma, то ∮Γf(z) dz=0\oint_{\Gamma } f(z)\, dz = 0.

§
Задача 7.10.1

Какие из следующих матриц сходящиеся, а какие суммируемые?

A=(−1/23/2−3/210−1/21−11/2),B=(010001100),C=(−1−2−3/2121113/2). \mathbf{A} = \begin{pmatrix} -1/2 & 3/2 & -3/2 \\ 1 & 0 & -1/2 \\ 1 & -1 & 1/2 \end{pmatrix}, \quad \mathbf{B} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & 0 & 0 \end{pmatrix}, \quad \mathbf{C} = \begin{pmatrix} -1 & -2 & -3/2 \\ 1 & 2 & 1 \\ 1 & 1 & 3/2 \end{pmatrix}.
?
Задача 7.10.2

Для матриц из упражнения 7.10.1 вычислите предел каждой сходящейся матрицы и вычислите предел Чезаро для каждой суммируемой матрицы.

Если A\mathbf{A} сходится или суммируема к G\mathbf{G}, и если I−A=BC\mathbf{I}-\mathbf{A} = \mathbf{B}\mathbf{C} --- факторизация полного ранга, то G=I−B(CB)−1C\mathbf{G} = \mathbf{I} - \mathbf{B}(\mathbf{C}\mathbf{B})^{-1}\mathbf{C}.

?
Задача 7.10.3

Проверьте, что

x→(k)=Akx→(0)иx→(k)=Akx→(0)+∑j=0k−1Ak−j−1b→(j)для k=1,2,3,… \overrightarrow {x}(k) = \mathbf{A}^k\overrightarrow {x}(0) \quad \text{и} \quad \overrightarrow {x}(k) = \mathbf{A}^k\overrightarrow {x}(0) + \sum _{j=0}^{k-1}\mathbf{A}^{k-j-1}\overrightarrow {b}(j) \quad \text{для } k=1,2,3,\ldots

действительно являются соответствующими решениями разностных уравнений x→(k+1)=Ax→(k)\overrightarrow {x}(k+1) = \mathbf{A}\overrightarrow {x}(k) (однородная система) и x→(k+1)=Ax→(k)+b→(k)\overrightarrow {x}(k+1) = \mathbf{A}\overrightarrow {x}(k)+\overrightarrow {b}(k) (неоднородная система).

?
Задача 7.10.4

Найдите предельный вектор для следующей игры в напёрстки, вычислив сначала предел Чезаро G\mathbf{G} с помощью факторизации полного ранга.

Игра в напёрстки: горошина изначально находится под одним из четырёх напёрстков и после каждого хода перемещается под один из (не более трёх) соседних напёрстков с равной вероятностью, что описывается 4×44\times 4 столбцово-стохастической матрицей перехода A\mathbf{A}, для которой p→(k)=Akp→(0)\overrightarrow {p}(k) = \mathbf{A}^k\overrightarrow {p}(0) задаёт распределение вероятностей по положениям напёрстка после kk ходов. Если I−A=BC\mathbf{I}-\mathbf{A} = \mathbf{B}\mathbf{C} — факторизация полного ранга, то предел Чезаро равен G=I−B(CB)−1C\mathbf{G} = \mathbf{I}-\mathbf{B}(\mathbf{C}\mathbf{B})^{-1}\mathbf{C}, а предельный вектор ожидаемых долей равен p→=Gp→(0)\overrightarrow {p} = \mathbf{G}\overrightarrow {p}(0).

?
Задача 7.10.5

Проверьте, что

x→(k)=Akx→(0)иx→(k)=Akx→(0)+∑j=0k−1Ak−j−1b→(j)при k=1,2,3,… \overrightarrow {x}(k) = \mathbf{A}^k\overrightarrow {x}(0) \quad \text{и} \quad \overrightarrow {x}(k) = \mathbf{A}^k\overrightarrow {x}(0) + \sum _{j=0}^{k-1}\mathbf{A}^{k-j-1}\overrightarrow {b}(j) \quad \text{при } k=1,2,3,\ldots

действительно являются соответствующими решениями разностных уравнений x→(k+1)=Ax→(k)\overrightarrow {x}(k+1) = \mathbf{A}\overrightarrow {x}(k) (однородная система) и x→(k+1)=Ax→(k)+b→(k)\overrightarrow {x}(k+1) = \mathbf{A}\overrightarrow {x}(k)+\overrightarrow {b}(k) (неоднородная система).

?
Задача 7.10.6

Докажите, что если существует матричная норма, для которой ∥A∥<1\left\| \mathbf{A}\right\| < 1, то lim⁡k→∞Ak=0\lim_{k\to \infty }\mathbf{A}^k = \mathbf{0}.

?
Задача 7.10.7

Исследуя итерационную матрицу, сравните сходимость метода Якоби и метода Гаусса--Зейделя для каждой из следующих матриц коэффициентов при произвольной правой части. Объясните, почему это показывает, что ни один из методов нельзя считать универсально предпочтительным по сравнению с другим.

A1=(12−2111221),A2=(2−11222−1−12). \mathbf{A}_1 = \begin{pmatrix} 1 & 2 & -2 \\ 1 & 1 & 1 \\ 2 & 2 & 1 \end{pmatrix}, \quad \mathbf{A}_2 = \begin{pmatrix} 2 & -1 & 1 \\ 2 & 2 & 2 \\ -1 & -1 & 2 \end{pmatrix}.
?
Задача 7.10.8

Пусть A=(2−10−12−10−12)\mathbf{A} = \begin{pmatrix} 2 & -1 & 0 \\ -1 & 2 & -1 \\ 0 & -1 & 2 \end{pmatrix} (матрица конечных разностей из примера 1.4.1).

?
(a)

Проверьте, что A\mathbf{A} удовлетворяет специальным условиям, гарантирующим справедливость формулы оптимальной релаксации SOR: det⁡(αD−L−U)=det⁡(αD−βL−β−1U)\operatorname {det}\left(\alpha \mathbf{D}-\mathbf{L}-\mathbf{U}\right) = \operatorname {det}\left(\alpha \mathbf{D}-\beta \mathbf{L}-\beta^{-1}\mathbf{U}\right) для всех вещественных α\alpha и β≠0\beta \ne 0, а итерационная матрица Якоби HJ\mathbf{H}_J имеет вещественные собственные значения, причём ρ(HJ)<1\rho (\mathbf{H}_J)<1, где D\mathbf{D}, L\mathbf{L}, U\mathbf{U} — диагональная, строго нижнетреугольная и строго верхнетреугольная части матрицы A\mathbf{A}.

(b)

Определите оптимальный параметр релаксации SOR, если известно, что ωopt=21+1−ρ2(HJ)\omega_{\text{opt}} = \dfrac {2}{1+\sqrt{1-\rho^2(\mathbf{H}_J)}} и ρ(Hωopt)=ωopt−1\rho (\mathbf{H}_{\omega_{\text{opt}}}) = \omega_{\text{opt}}-1.

(c)

Найдите асимптотические скорости сходимости для методов Якоби, Гаусса--Зейделя и оптимального SOR.

(d)

Используя x→(0)=(1,1,1)T\overrightarrow {x}(0) = (1,1,1)^T и b→=(2,4,6)T\overrightarrow {b} = (2,4,6)^T, выполните несколько шагов методов Якоби, Гаусса--Зейделя и оптимального SOR для решения системы Ax→=b→\mathbf{A}\overrightarrow {x}=\overrightarrow {b}, пока не станет виден характер сходимости.

Задача 7.10.9

Докажите, что если ρ(Hω)<1\rho (\mathbf{H}_{\omega }) < 1, где Hω\mathbf{H}_{\omega } — итерационная матрица метода SOR, то 0<ω<20 < \omega < 2.

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

Указание: используйте det⁡(Hω)\operatorname {det}\left(\mathbf{H}_{\omega }\right), чтобы показать, что ∣λk∣≥∣1−ω∣\left|\lambda_k\right| \ge \left|1-\omega \right| для некоторого λk∈σ(Hω)\lambda_k \in \sigma (\mathbf{H}_{\omega }).

Задача 7.10.10

Покажите, что спектральный радиус итерационной матрицы Якоби для дискретного лапласиана Ln2×n2\mathbf{L}_{n^2\times n^2}, описанного в примере 7.6.2, равен ρ(HJ)=cos⁡(π/(n+1))\rho (\mathbf{H}_J) = \cos \bigl(\pi /(n+1)\bigr).

?
Задача 7.10.11

Рассмотрим скалярную последовательность {α1,α2,α3,…}\left\{ \alpha_1, \alpha_2, \alpha_3, \ldots \right\} и связанную с ней последовательность средних Чезаро {μ1,μ2,μ3,…}\left\{ \mu_1, \mu_2, \mu_3, \ldots \right\}, где μn=(α1+α2+⋯+αn)/n\mu_n = (\alpha_1+\alpha_2+\cdots +\alpha_n)/n. Докажите, что если {αn}\left\{ \alpha_n\right\} сходится к α\alpha, то {μn}\left\{ \mu_n\right\} также сходится к α\alpha.

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

Как и для скаляров, векторная последовательность {v→n}\left\{ \overrightarrow {v}_n\right\} в конечномерном пространстве сходится к v→\overrightarrow {v} тогда и только тогда, когда для каждого ε>0\varepsilon >0 найдётся натуральное число N=N(ε)N=N(\varepsilon ), такое что ∥v→n−v→∥<ε\left\| \overrightarrow {v}_n-\overrightarrow {v}\right\| <\varepsilon для всех n≥Nn\ge N, причём неважно, какая норма используется. Поэтому ваше доказательство должно быть справедливо и для векторов (и матриц).

Задача 7.10.12

Для матриц с неположительными внедиагональными элементами (ZZ-матриц) докажите, что следующие утверждения эквивалентны.

?
(a)

A\mathbf{A} является MM-матрицей.

(b)

Все ведущие главные миноры A\mathbf{A} положительны.

(c)

A\mathbf{A} имеет LU\mathbf{L}\mathbf{U}-разложение, причём и L\mathbf{L}, и U\mathbf{U} являются MM-матрицами.

(d)

Существует вектор x→>0→\overrightarrow {x}>\overrightarrow {0} такой, что Ax→>0→\mathbf{A}\overrightarrow {x}>\overrightarrow {0}.

(e)

Каждый aii>0a_{ii}>0 и AD\mathbf{A}\mathbf{D} диагонально доминантна для некоторой диагональной матрицы D\mathbf{D} с положительными диагональными элементами.

(f)

Ax→≥0→\mathbf{A}\overrightarrow {x} \ge \overrightarrow {0} влечёт x→≥0→\overrightarrow {x} \ge \overrightarrow {0}.

Задача 7.10.13

Предположим, что λ∈σ(A)\lambda \in \sigma (\mathbf{A}), и пусть M1=A−λI\mathbf{M}_1 = \mathbf{A}-\lambda \mathbf{I}. Следующая процедура даёт значение index⁡(λ)\operatorname {index}(\lambda ): разложим M1=B1C1\mathbf{M}_1 = \mathbf{B}_1\mathbf{C}_1 как факторизацию полного ранга, положим M2=C1B1\mathbf{M}_2 = \mathbf{C}_1\mathbf{B}_1, разложим M2=B2C2\mathbf{M}_2 = \mathbf{B}_2\mathbf{C}_2 как факторизацию полного ранга, положим M3=C2B2\mathbf{M}_3 = \mathbf{C}_2\mathbf{B}_2, и так далее. В общем случае Mi=Ci−1Bi−1\mathbf{M}_i = \mathbf{C}_{i-1}\mathbf{B}_{i-1}, где Mi−1=Bi−1Ci−1\mathbf{M}_{i-1} = \mathbf{B}_{i-1}\mathbf{C}_{i-1} --- факторизация полного ранга.

?
(a)

Объясните, почему эта процедура должна в конечном счёте дать матрицу Mk\mathbf{M}_k, которая либо невырождена, либо нулевая.

(b)

Докажите, что если kk --- наименьшее натуральное число такое, что Mk−1\mathbf{M}_k^{-1} существует или Mk=0\mathbf{M}_k = \mathbf{0}, то

index⁡(λ)={k−1если Mk невырождена,kесли Mk=0. \operatorname {index}(\lambda ) = \begin{cases} k-1 & \text{если } \mathbf{M}_k \text{ невырождена}, \\ k & \text{если } \mathbf{M}_k = \mathbf{0}. \end{cases}
Задача 7.10.14

Используя процедуру из задачи 7.10.13, найдите индекс каждого собственного значения A=(−3−8−95119−1−21)\mathbf{A} = \begin{pmatrix} -3 & -8 & -9 \\ 5 & 11 & 9 \\ -1 & -2 & 1 \end{pmatrix}.

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

Подсказка: σ(A)={4,1}\sigma (\mathbf{A}) = \left\{ 4,1\right\}.

Задача 7.10.15

Пусть A\mathbf{A} --- матрица, заданная в задаче 7.10.14.

?
(a)

Найдите жорданову форму для A\mathbf{A}.

(b)

Для произвольной функции ff, определённой на A\mathbf{A}, найдите интерполяционный полином Эрмита и опишите f(A)f(\mathbf{A}).

Задача 7.10.16

Пусть дана матрица Bn×n\mathbf{B}_{n\times n} ранга rr такая, что index⁡(B)≤1\operatorname {index}(\mathbf{B}) \le 1 (т.е. index⁡(λ=0)≤1\operatorname {index}(\lambda =0)\le 1), жорданова форма для B\mathbf{B} имеет вид (000Cr×r)=P−1BP\begin{pmatrix} 0 & 0 \\ 0 & \mathbf{C}_{r\times r} \end{pmatrix} = \mathbf{P}^{-1}\mathbf{B}\mathbf{P}, так что B=P(000C)P−1\mathbf{B} = \mathbf{P}\begin{pmatrix} 0 & 0 \\ 0 & \mathbf{C} \end{pmatrix}\mathbf{P}^{-1}, где C\mathbf{C} невырождена. Это означает, что B\mathbf{B} принадлежит алгебраической группе G\mathcal{G} относительно умножения матриц, и обратный к B\mathbf{B} элемент в G\mathcal{G} есть B#=P(000C−1)P−1\mathbf{B}^{\# } = \mathbf{P}\begin{pmatrix} 0 & 0 \\ 0 & \mathbf{C}^{-1} \end{pmatrix}\mathbf{P}^{-1}; B#\mathbf{B}^{\# } называется групповым обратным для B\mathbf{B}. Докажите, что если lim⁡k→∞Ak\lim_{k\to \infty }\mathbf{A}^k существует и B=I−A\mathbf{B} = \mathbf{I}-\mathbf{A}, то

lim⁡k→∞Ak=I−BB#. \lim _{k\to \infty }\mathbf{A}^k = \mathbf{I}-\mathbf{B}\mathbf{B}^{\# }.

Иными словами, предельную матрицу можно охарактеризовать как разность двух единичных элементов --- I\mathbf{I} является единицей в мультипликативной группе невырожденных матриц, а BB#\mathbf{B}\mathbf{B}^{\# } является единичным элементом в мультипликативной группе, содержащей B\mathbf{B}.

Пределы степеней: для A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n}, lim⁡k→∞Ak\lim_{k\to \infty }\mathbf{A}^k существует тогда и только тогда, когда ρ(A)<1\rho (\mathbf{A})<1, либо ρ(A)=1\rho (\mathbf{A})=1, причём λ=1\lambda =1 --- единственное собственное значение на единичной окружности, и λ=1\lambda =1 полупростое.

?
Задача 7.10.17

Если Mn×n\mathbf{M}_{n\times n} — групповая матрица (т.е. если index⁡(M)≤1\operatorname {index}(\mathbf{M}) \le 1), то групповую обратную матрицу M\mathbf{M} можно охарактеризовать как единственное решение M#\mathbf{M}^{\# } уравнений MM#M=M\mathbf{M}\mathbf{M}^{\# }\mathbf{M} = \mathbf{M}, M#MM#=M#\mathbf{M}^{\# }\mathbf{M}\mathbf{M}^{\# } = \mathbf{M}^{\# } и MM#=M#M\mathbf{M}\mathbf{M}^{\# } = \mathbf{M}^{\# }\mathbf{M}. Более того, некоторые авторы используют именно эти уравнения для определения M#\mathbf{M}^{\# }. Используя эту характеризацию, покажите, что если M=BC\mathbf{M} = \mathbf{B}\mathbf{C} — произвольное разложение M\mathbf{M} на множители полного ранга, то M#=B(CB)−2C\mathbf{M}^{\# } = \mathbf{B}(\mathbf{C}\mathbf{B})^{-2}\mathbf{C}. В частности, если M=U1DV1∗\mathbf{M} = \mathbf{U}_1\mathbf{D}\mathbf{V}_1^{*} — разложение полного ранга, полученное из сингулярного разложения, то

M#=U1D−1/2(V1∗U1)−2D−1/2V1∗=U1D−1(V1∗U1)−2V1∗=U1(V1∗U1)−2D−1V1∗. \mathbf{M}^{\# } = \mathbf{U}_1\mathbf{D}^{-1/2}(\mathbf{V}_1^{*}\mathbf{U}_1)^{-2}\mathbf{D}^{-1/2}\mathbf{V}_1^{*} = \mathbf{U}_1\mathbf{D}^{-1}(\mathbf{V}_1^{*}\mathbf{U}_1)^{-2}\mathbf{V}_1^{*} = \mathbf{U}_1(\mathbf{V}_1^{*}\mathbf{U}_1)^{-2}\mathbf{D}^{-1}\mathbf{V}_1^{*}.
?
§
Задача 7.11.1

Найдите минимальный многочлен для A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix}.

?
Задача 7.11.2

Найдите минимальный многочлен вектора b→=(−1,1,1)T\overrightarrow {b} = (-1,1,1)^T относительно матрицы A\mathbf{A} из упражнения 7.11.1.

?
Задача 7.11.3

Используя метод Крылова, найдите характеристический многочлен матрицы A\mathbf{A} из упражнения 7.11.1.

?
Задача 7.11.4

Какова жорданова форма матрицы, минимальный многочлен которой равен m(x)=(x−λ)(x−μ)2m(x) = (x-\lambda )(x-\mu )^2, а характеристический многочлен — c(x)=(x−λ)2(x−μ)4c(x) = (x-\lambda )^2(x-\mu )^4?

?
Задача 7.11.5

Используя метод вычисления минимального многочлена матрицы через ортогонализованную крыловскую последовательность (последовательно ортогонализуя по Граму--Шмидту I,A,A2,…\mathbf{I}, \mathbf{A}, \mathbf{A}^2, \ldots в скалярном произведении Фробениуса ⟨X∣Y⟩=trace⁡(XTY)\langle \mathbf{X} \mid \mathbf{Y} \rangle = \operatorname {trace}(\mathbf{X}^T\mathbf{Y}) до тех пор, пока не окажется, что степень Ak\mathbf{A}^k является линейной комбинацией предыдущих ортогонализованных членов U0,…,Uk−1\mathbf{U}_0, \ldots , \mathbf{U}_{k-1}, причём коэффициенты находятся решением треугольной системы Rc→=R−1\mathbf{R}\overrightarrow {c} = \mathbf{R}^{-1} относительно αi\alpha_i), найдите минимальный многочлен для

A=(−7−48−8−4−14−4−16−817−16−6−36−5). \mathbf{A} = \begin{pmatrix} -7 & -4 & 8 & -8 \\ -4 & -1 & 4 & -4 \\ -16 & -8 & 17 & -16 \\ -6 & -3 & 6 & -5 \end{pmatrix}.
?
Задача 7.11.6

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

?
Задача 7.11.7

Покажите, что две матрицы могут иметь одинаковые минимальный и характеристический многочлены, не будучи подобными, рассмотрев A=(N00N)\mathbf{A} = \begin{pmatrix} \mathbf{N} & \mathbf{0} \\ \mathbf{0} & \mathbf{N} \end{pmatrix} и B=(N000)\mathbf{B} = \begin{pmatrix} \mathbf{N} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}, где N=(0100)\mathbf{N} = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}.

?
Задача 7.11.8

Докажите, что если A\mathbf{A} и B\mathbf{B} — циклические матрицы с одинаковым характеристическим многочленом, то A\mathbf{A} подобна B\mathbf{B}.

?
Задача 7.11.9

Используя алгоритм Ланцоша, найдите ортогональную матрицу P\mathbf{P} такую, что PTAP=T\mathbf{P}^T\mathbf{A}\mathbf{P} = \mathbf{T} трёхдиагональна, где A=(211121112)\mathbf{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}.

?
Задача 7.11.10

Начиная с x→0=0→\overrightarrow {x}_0 = \overrightarrow {0}, примените алгоритм сопряжённых градиентов для решения Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, где A=(211121112)\mathbf{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix} и b→=(400)\overrightarrow {b} = \begin{pmatrix} 4 \\ 0 \\ 0 \end{pmatrix}.

?
Задача 7.11.11

Используя алгоритм Арнольди, найдите ортогональную матрицу Q\mathbf{Q} такую, что QTAQ=H\mathbf{Q}^T\mathbf{A}\mathbf{Q} = \mathbf{H} имеет верхнюю форму Хессенберга, где A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix}.

?
Задача 7.11.12

Используя GMRES, решите Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b} для A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix} и b→=(121)\overrightarrow {b} = \begin{pmatrix} 1 \\ 2 \\ 1 \end{pmatrix}.

?