4.5

Ещё о ранге

[20/100%]
Показать
LaTeX
Задача 4.5.1

Проверьте, что rk⁡(ATA)=rk⁡(A)=rk⁡(AAT)\operatorname {rk}\left(\mathbf{A}^T \mathbf{A}\right) = \operatorname {rk}\left(\mathbf{A}\right) = \operatorname {rk}\left(\mathbf{A}\mathbf{A}^T\right) для

A=(131−4−1−310262−8). \mathbf{A} = \begin{pmatrix} 1 & 3 & 1 & -4 \\ -1 & -3 & 1 & 0 \\ 2 & 6 & 2 & -8 \end{pmatrix}.
?
Задача 4.5.2

Определите dim⁡(N(A)∩R(B))\dim \left(N(\mathbf{A}) \cap R(\mathbf{B})\right) для

A=(−211−422000)иB=(131−4−1−310262−8). \mathbf{A} = \begin{pmatrix} -2 & 1 & 1 \\ -4 & 2 & 2 \\ 0 & 0 & 0 \end{pmatrix} \quad \text{и} \quad \mathbf{B} = \begin{pmatrix} 1 & 3 & 1 & -4 \\ -1 & -3 & 1 & 0 \\ 2 & 6 & 2 & -8 \end{pmatrix}.
?
Задача 4.5.3

Для матриц, заданных в упражнении 4.5.2, воспользуйтесь процедурой, описанной на стр.211, чтобы определить базис пространства N(A)∩R(B)N(\mathbf{A}) \cap R(\mathbf{B}).

Процедура: для Am×n\mathbf{A}_{m \times n} и Bn×p\mathbf{B}_{n \times p} базис пространства N(A)∩R(B)N(\mathbf{A}) \cap R(\mathbf{B}) строится следующим образом — находим базис {x→1,…,x→r}\left\{ \overrightarrow {x}_1, \ldots , \overrightarrow {x}_r\right\} пространства R(B)R(\mathbf{B}), полагаем Xn×r=(x→1∣x→2∣⋯∣x→r)\mathbf{X}_{n \times r} = (\overrightarrow {x}_1 \mid \overrightarrow {x}_2 \mid \cdots \mid \overrightarrow {x}_r), находим базис {v→1,…,v→s}\left\{ \overrightarrow {v}_1, \ldots , \overrightarrow {v}_s\right\} пространства N(AX)N(\mathbf{A}\mathbf{X}), и тогда B={Xv→1,Xv→2,…,Xv→s}\mathcal{B} = \left\{ \mathbf{X}\overrightarrow {v}_1, \mathbf{X}\overrightarrow {v}_2, \ldots , \mathbf{X}\overrightarrow {v}_s\right\} — базис пространства N(A)∩R(B)N(\mathbf{A}) \cap R(\mathbf{B}).

?
Задача 4.5.4

Если A1A2⋯Ak\mathbf{A}_1 \mathbf{A}_2 \cdots \mathbf{A}_k — произведение квадратных матриц, в котором некоторая Ai\mathbf{A}_i вырождена, объясните, почему всё произведение должно быть вырожденным.

?
Задача 4.5.5

Для A∈Rm×n\mathbf{A} \in \mathbb {R}^{m \times n} объясните, почему из ATA=0\mathbf{A}^T \mathbf{A} = \mathbf{0} следует A=0\mathbf{A} = \mathbf{0}.

?
Задача 4.5.6

Найдите rk⁡(A)\operatorname {rk}\left(\mathbf{A}\right) и все невырожденные подматрицы максимального порядка в

A=(2−114−218−41). \mathbf{A} = \begin{pmatrix} 2 & -1 & 1 \\ 4 & -2 & 1 \\ 8 & -4 & 1 \end{pmatrix}.
?
Задача 4.5.7

Возможно ли, чтобы rk⁡(AB)<rk⁡(A)\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) < \operatorname {rk}\left(\mathbf{A}\right) и rk⁡(AB)<rk⁡(B)\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) < \operatorname {rk}\left(\mathbf{B}\right) для одной и той же пары матриц?

?
Задача 4.5.8

Верно ли, что rk⁡(AB)=rk⁡(BA)\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) = \operatorname {rk}\left(\mathbf{B}\mathbf{A}\right), когда оба произведения определены? Почему?

?
Задача 4.5.9

Объясните, почему rk⁡(AB)=rk⁡(A)−dim⁡(N(BT)∩R(AT))\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) = \operatorname {rk}\left(\mathbf{A}\right) - \dim \left(N(\mathbf{B}^T) \cap R(\mathbf{A}^T)\right).

?
Задача 4.5.10

Объясните, почему dim⁡N(Am×nBn×p)=dim⁡N(B)+dim⁡(R(B)∩N(A))\dim N(\mathbf{A}_{m \times n} \mathbf{B}_{n \times p}) = \dim N(\mathbf{B}) + \dim \left(R(\mathbf{B}) \cap N(\mathbf{A})\right).

?
Задача 4.5.11

Закон Сильвестра о нуль-мерности, сформулированный Джеймсом Дж. Сильвестром в 1884 году, утверждает, что для квадратных матриц A\mathbf{A} и B\mathbf{B}

max⁡{ν(A),ν(B)}≤ν(AB)≤ν(A)+ν(B), \max \left\{ \nu (\mathbf{A}), \nu (\mathbf{B})\right\} \leq \nu (\mathbf{A}\mathbf{B}) \leq \nu (\mathbf{A}) + \nu (\mathbf{B}),

где ν(⋆)=dim⁡N(⋆)\nu (\star ) = \dim N(\star ) обозначает нуль-мерность (размерность ядра).

?
(a)

Установите справедливость закона Сильвестра.

(b)

Покажите, что закон Сильвестра не выполняется для прямоугольных матриц, поскольку возможно ν(A)>ν(AB)\nu (\mathbf{A}) > \nu (\mathbf{A}\mathbf{B}). Возможно ли ν(B)>ν(AB)\nu (\mathbf{B}) > \nu (\mathbf{A}\mathbf{B})?

Задача 4.5.12

Для матриц Am×n\mathbf{A}_{m \times n} и Bn×p\mathbf{B}_{n \times p} докажите каждое из следующих утверждений:

?
(a)

rk⁡(AB)=rk⁡(A)\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) = \operatorname {rk}\left(\mathbf{A}\right) и R(AB)=R(A)R(\mathbf{A}\mathbf{B}) = R(\mathbf{A}), если rk⁡(B)=n\operatorname {rk}\left(\mathbf{B}\right) = n.

(b)

rk⁡(AB)=rk⁡(B)\operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) = \operatorname {rk}\left(\mathbf{B}\right) и N(AB)=N(B)N(\mathbf{A}\mathbf{B}) = N(\mathbf{B}), если rk⁡(A)=n\operatorname {rk}\left(\mathbf{A}\right) = n.

Задача 4.5.13

Выполните следующие вычисления, используя матрицы:

A=(122412.01)иb→=(121.01). \mathbf{A} = \begin{pmatrix} 1 & 2 \\ 2 & 4 \\ 1 & 2.01 \end{pmatrix} \quad \text{и} \quad \overrightarrow {b} = \begin{pmatrix} 1 \\ 2 \\ 1.01 \end{pmatrix}.
?
(a)

Найдите rk⁡(A)\operatorname {rk}\left(\mathbf{A}\right) и решите Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, используя точную арифметику.

(b)

Найдите rk⁡(ATA)\operatorname {rk}\left(\mathbf{A}^T \mathbf{A}\right) и решите ATAx→=ATb→\mathbf{A}^T \mathbf{A}\overrightarrow {x} = \mathbf{A}^T \overrightarrow {b} точно.

(c)

Найдите rk⁡(A)\operatorname {rk}\left(\mathbf{A}\right) и решите Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b} с 3-значной арифметикой.

(d)

Найдите ATA\mathbf{A}^T \mathbf{A}, ATb→\mathbf{A}^T \overrightarrow {b} и решение ATAx→=ATb→\mathbf{A}^T \mathbf{A}\overrightarrow {x} = \mathbf{A}^T \overrightarrow {b} с 3-значной арифметикой.

Задача 4.5.14

Докажите, что если элементы матрицы Fr×r\mathbf{F}_{r \times r} удовлетворяют ∑j=1r∣fij∣<1\sum_{j=1}^{r} \left|f_{ij}\right| < 1 для каждого ii (т.е. каждая сумма модулей элементов строки <1< 1), то I+F\mathbf{I} + \mathbf{F} невырождена.

?
Задача 4.5.15

Если A=(WXYZ)\mathbf{A} = \begin{pmatrix} \mathbf{W} & \mathbf{X} \\ \mathbf{Y} & \mathbf{Z} \end{pmatrix}, где rk⁡(A)=r=rk⁡(Wr×r)\operatorname {rk}\left(\mathbf{A}\right) = r = \operatorname {rk}\left(\mathbf{W}_{r \times r}\right), покажите, что существуют матрицы B\mathbf{B} и C\mathbf{C} такие, что

A=(WWCBWBWC)=(IB)W(IC). \mathbf{A} = \begin{pmatrix} \mathbf{W} & \mathbf{W}\mathbf{C} \\ \mathbf{B}\mathbf{W} & \mathbf{B}\mathbf{W}\mathbf{C} \end{pmatrix} = \begin{pmatrix} \mathbf{I} \\ \mathbf{B} \end{pmatrix} \mathbf{W} \begin{pmatrix} \mathbf{I} & \mathbf{C} \end{pmatrix}.
?
Задача 4.5.16

Для сходящейся последовательности {Ak}k=1∞\left\{ \mathbf{A}_k\right\}_{k=1}^{\infty } матриц пусть A=lim⁡k→∞Ak\mathbf{A} = \lim_{k \rightarrow \infty } \mathbf{A}_k.

?
(a)

Докажите, что если каждая Ak\mathbf{A}_k вырождена, то A\mathbf{A} вырождена.

(b)

Если каждая Ak\mathbf{A}_k невырождена, обязана ли A\mathbf{A} быть невырожденной? Почему?

Задача 4.5.17

Установите справедливость результата Фробениуса 1911 года, утверждающего, что если ABC\mathbf{A}\mathbf{B}\mathbf{C} существует, то

rk⁡(AB)+rk⁡(BC)≤rk⁡(B)+rk⁡(ABC). \operatorname {rk}\left(\mathbf{A}\mathbf{B}\right) + \operatorname {rk}\left(\mathbf{B}\mathbf{C}\right) \leq \operatorname {rk}\left(\mathbf{B}\right) + \operatorname {rk}\left(\mathbf{A}\mathbf{B}\mathbf{C}\right).
?
Задача 4.5.18

Если A\mathbf{A} имеет размер n×nn \times n, докажите, что следующие утверждения эквивалентны:

  1. N(A)=N(A2)N(\mathbf{A}) = N(\mathbf{A}^2).

  2. R(A)=R(A2)R(\mathbf{A}) = R(\mathbf{A}^2).

  3. R(A)∩N(A)={0→}R(\mathbf{A}) \cap N(\mathbf{A}) = \left\{ \overrightarrow {0}\right\}.

?
Задача 4.5.19

Пусть A\mathbf{A} и B\mathbf{B} — матрицы размера n×nn \times n, такие что A=A2\mathbf{A} = \mathbf{A}^2, B=B2\mathbf{B} = \mathbf{B}^2 и AB=BA=0\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A} = \mathbf{0}.

?
(a)

Докажите, что rk⁡(A+B)=rk⁡(A)+rk⁡(B)\operatorname {rk}\left(\mathbf{A} + \mathbf{B}\right) = \operatorname {rk}\left(\mathbf{A}\right) + \operatorname {rk}\left(\mathbf{B}\right).

(b)

Докажите, что rk⁡(A)+rk⁡(I−A)=n\operatorname {rk}\left(\mathbf{A}\right) + \operatorname {rk}\left(\mathbf{I} - \mathbf{A}\right) = n.

Задача 4.5.20

Для A∈Rm×n\mathbf{A} \in \mathbb {R}^{m \times n} такой, что rk⁡(A)=r\operatorname {rk}\left(\mathbf{A}\right) = r, пусть A=BC\mathbf{A} = \mathbf{B}\mathbf{C} — разложение A\mathbf{A} полного ранга, в котором Bm×r\mathbf{B}_{m \times r} — матрица базисных столбцов из A\mathbf{A}, а Cr×n\mathbf{C}_{r \times n} — матрица ненулевых строк из EA\mathbf{E}_{\mathbf{A}} (см. упражнение 3.9.8). Матрица, определённая как

A†=CT(BTACT)−1BT \mathbf{A}^{\dagger } = \mathbf{C}^T \left(\mathbf{B}^T \mathbf{A} \mathbf{C}^T\right)^{-1} \mathbf{B}^T

называется псевдообратной матрицей Мура–Пенроуза для A\mathbf{A}. Некоторые авторы называют A†\mathbf{A}^{\dagger } псевдообратной или обобщённой обратной матрицей для A\mathbf{A}.

?
(a)

Объясните, почему матрица BTACT\mathbf{B}^T \mathbf{A} \mathbf{C}^T невырождена.

(b)

Убедитесь, что x→=A†b→\overrightarrow {x} = \mathbf{A}^{\dagger } \overrightarrow {b} является решением нормальных уравнений ATAx→=ATb→\mathbf{A}^T \mathbf{A} \overrightarrow {x} = \mathbf{A}^T \overrightarrow {b} (а также Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, когда эта система совместна).

(c)

Покажите, что общее решение для ATAx→=ATb→\mathbf{A}^T \mathbf{A} \overrightarrow {x} = \mathbf{A}^T \overrightarrow {b} (а также Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, когда эта система совместна) можно записать как

x→=A†b→+(I−A†A)h→, \overrightarrow {x} = \mathbf{A}^{\dagger } \overrightarrow {b} + \left(\mathbf{I} - \mathbf{A}^{\dagger } \mathbf{A}\right) \overrightarrow {h},

где h→\overrightarrow {h} — «свободный» вектор-параметр из Rn×1\mathbb {R}^{n \times 1}.

(d)

Если r=nr = n, объясните, почему A†=(ATA)−1AT\mathbf{A}^{\dagger } = \left(\mathbf{A}^T \mathbf{A}\right)^{-1} \mathbf{A}^T.

(e)

Если A\mathbf{A} квадратная и невырожденная, объясните, почему A†=A−1\mathbf{A}^{\dagger } = \mathbf{A}^{-1}.

(f)

Убедитесь, что A†=CT(BTACT)−1BT\mathbf{A}^{\dagger } = \mathbf{C}^T \left(\mathbf{B}^T \mathbf{A} \mathbf{C}^T\right)^{-1} \mathbf{B}^T удовлетворяет уравнениям Пенроуза:

AA†A=A,(AA†)T=AA†,A†AA†=A†,(A†A)T=A†A. \mathbf{A}\mathbf{A}^{\dagger }\mathbf{A} = \mathbf{A}, \qquad \left(\mathbf{A}\mathbf{A}^{\dagger }\right)^T = \mathbf{A}\mathbf{A}^{\dagger }, \qquad \mathbf{A}^{\dagger }\mathbf{A}\mathbf{A}^{\dagger } = \mathbf{A}^{\dagger }, \qquad \left(\mathbf{A}^{\dagger }\mathbf{A}\right)^T = \mathbf{A}^{\dagger }\mathbf{A}.