7.3

Функции диагонализируемых матриц

[17/88%]
Показать
LaTeX
Задача 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-итерацию с двойным сдвигом.