Глава 3

Матричная алгебра

[78/100%]
Показать
LaTeX
§
Задача 3.2.1

Определите неизвестные величины в следующих выражениях.

?
(a)

3X=[0369]3\mathbf{X} = \begin{bmatrix} 0 & 3 \\ 6 & 9 \end{bmatrix}.

(b)

2[x+2y+330]=[36yz]T2\begin{bmatrix} x+2 & y+3 \\ 3 & 0 \end{bmatrix} = \begin{bmatrix} 3 & 6 \\ y & z \end{bmatrix}^{T}.

Задача 3.2.2

Определите, является ли каждая из следующих матриц симметричной, кососимметричной или ни той, ни другой.

?
(a)

[1−33−34−3330]\begin{bmatrix} 1 & -3 & 3 \\ -3 & 4 & -3 \\ 3 & 3 & 0 \end{bmatrix}.

(b)

[0−3−33013−10]\begin{bmatrix} 0 & -3 & -3 \\ 3 & 0 & 1 \\ 3 & -1 & 0 \end{bmatrix}.

(c)

[0−3−3−303−331]\begin{bmatrix} 0 & -3 & -3 \\ -3 & 0 & 3 \\ -3 & 3 & 1 \end{bmatrix}.

(d)

[120210]\begin{bmatrix} 1 & 2 & 0 \\ 2 & 1 & 0 \end{bmatrix}.

Задача 3.2.3

Постройте пример матрицы 3×33 \times 3 A\mathbf{A}, удовлетворяющей следующим условиям.

?
(a)

A\mathbf{A} одновременно симметрична и кососимметрична.

(b)

A\mathbf{A} одновременно эрмитова и симметрична.

(c)

A\mathbf{A} косоэрмитова.

Задача 3.2.4

Объясните, почему множество всех симметричных матриц размера n×nn \times n замкнуто относительно сложения матриц. То есть объясните, почему сумма двух симметричных матриц размера n×nn \times n снова является симметричной матрицей размера n×nn \times n. Замкнуто ли множество всех кососимметричных матриц размера n×nn \times n относительно сложения матриц?

?
Задача 3.2.5

Докажите, что каждое из следующих утверждений верно.

?
(a)

Если A=[aij]\mathbf{A} = [a_{ij}] кососимметрична, то ajj=0a_{jj} = 0 для каждого jj.

(b)

Если A=[aij]\mathbf{A} = [a_{ij}] косоэрмитова, то каждое ajja_{jj} — чисто мнимое число, т.е. кратное мнимой единице i\mathrm{i}.

(c)

Если A\mathbf{A} вещественна и симметрична, то B=iA\mathbf{B} = \mathrm{i}\mathbf{A} косоэрмитова.

Задача 3.2.6

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

?
(a)

Покажите, что A+AT\mathbf{A}+\mathbf{A}^{T} симметрична, а A−AT\mathbf{A}-\mathbf{A}^{T} кососимметрична.

(b)

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

Задача 3.2.7

Если A\mathbf{A} и B\mathbf{B} — две матрицы одинаковой формы, докажите, что каждое из следующих утверждений верно.

?
(a)

(A+B)∗=A∗+B∗(\mathbf{A}+\mathbf{B})^{*} = \mathbf{A}^{*}+\mathbf{B}^{*}.

(b)

(αA)∗=αˉA∗(\alpha \mathbf{A})^{*} = \bar{\alpha }\mathbf{A}^{*}.

Задача 3.2.8

Напомним постановку задачи о системе пружин из примера 3.2.1: для ряда пружин, соединяющих узлы, если узел ii смещён на xix_i единиц (положительное направление — влево, отрицательное — вправо), а соседние узлы соединены пружиной с коэффициентом жёсткости kk, то сила, действующая на внутренний узел со стороны двух соседних пружин, равна сумме произведений kk на относительное смещение с каждым из соседей; при этом, если полученные уравнения баланса сил организовать в линейную систему, коэффициентная матрица оказывается симметричной матрицей жёсткости K\mathbf{K}. Для двух пружин, соединяющих три узла, с общим коэффициентом жёсткости k1=k2=kk_1 = k_2 = k, эта конструкция даёт

K=k[1−10−12−10−11]. \mathbf{K} = k \begin{bmatrix} 1 & -1 & 0 \\ -1 & 2 & -1 \\ 0 & -1 & 1 \end{bmatrix}.

Используя соглашения примера 3.2.1, определите матрицу жёсткости для системы из nn одинаковых пружин с коэффициентом жёсткости kk, соединённых в ряд, аналогично изображённому на рисунке 3.2.1.

Рисунок 3.2.1Рисунок 3.2.1

?
§
Задача 3.3.1

Каждая из следующих функций отображает R2\mathbb {R}^{2} в R2\mathbb {R}^{2}. Определите, какие из них являются линейными функциями.

?
(a)

f(xy)=(x1+y)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} x \\ 1+y \end{pmatrix}.

(b)

f(xy)=(yx)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} y \\ x \end{pmatrix}.

(c)

f(xy)=(0xy)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ xy \end{pmatrix}.

(d)

f(xy)=(x2y2)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} x^{2} \\ y^{2} \end{pmatrix}.

(e)

f(xy)=(xsin⁡y)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} x \\ \sin y \end{pmatrix}.

(f)

f(xy)=(x+yx−y)f\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} x+y \\ x-y \end{pmatrix}.

Задача 3.3.2

Для x→=(x1x2⋮xn)\overrightarrow {x} = \begin{pmatrix} x_{1} \\ x_{2} \\ \vdots \\ x_{n} \end{pmatrix} и для констант ξi\xi_{i} убедитесь, что

f(x→)=ξ1x1+ξ2x2+⋯+ξnxn f(\overrightarrow {x}) = \xi _{1}x_{1}+\xi _{2}x_{2}+\cdots +\xi _{n}x_{n}

является линейной функцией.

?
Задача 3.3.3

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

?
Задача 3.3.4

Определите, какие из следующих трёх преобразований в R2\mathbb {R}^{2} являются линейными: (1) поворот точек против часовой стрелки на угол θ\theta; (2) отражение точек относительно оси xx; (3) проекция точек на прямую y=xy=x в перпендикулярном направлении.

Поворот, отражение и проекция в \mathbb {R}^{2}Поворот, отражение и проекция в \mathbb {R}^{2}

?
§
Задача 3.4.1

Определите матрицу, связанную с каждой из трёх линейных функций из R2\mathbb {R}^{2} в R2\mathbb {R}^{2}, описанных ниже. То есть определите aija_{ij}, такие что

f(p→)=f(x1x2)=(a11x1+a12x2a21x1+a22x2). f(\overrightarrow {p}) = f\begin{pmatrix} x_{1} \\ x_{2} \end{pmatrix} = \begin{pmatrix} a_{11}x_{1}+a_{12}x_{2} \\ a_{21}x_{1}+a_{22}x_{2} \end{pmatrix}.
?
Примечание.
?

Три линейные функции, о которых идёт речь, используются на протяжении всего этого раздела: поворот точек против часовой стрелки на угол θ\theta; отражение точек относительно оси xx; и проекция точек на прямую y=xy=x в перпендикулярном направлении.

Три линейных преобразования, используемые на протяжении раздела 3.4Три линейных преобразования, используемые на протяжении раздела 3.4

Задача 3.4.2

Используя умножение матриц, определите линейную функцию, получаемую в результате поворота с последующим отражением (используя поворот и отражение из упражнения 3.4.1).

?
Задача 3.4.3

Используя умножение матриц, определите линейную функцию, получаемую в результате сначала отражения, затем поворота и, наконец, проекции (используя отражение, поворот и проекцию из упражнения 3.4.1).

?
§
Задача 3.5.1

Для A=[1−230−544−38]\mathbf{A} = \begin{bmatrix} 1 & -2 & 3 \\ 0 & -5 & 4 \\ 4 & -3 & 8 \end{bmatrix}, B=[120437]\mathbf{B} = \begin{bmatrix} 1 & 2 \\ 0 & 4 \\ 3 & 7 \end{bmatrix} и C=[123]\mathbf{C} = \begin{bmatrix} 1 \\ 2 \\ 3 \end{bmatrix} вычислите следующие произведения, когда это возможно.

?
(a)

AB\mathbf{A}\mathbf{B}.

(b)

BA\mathbf{B}\mathbf{A}.

(c)

CB\mathbf{C}\mathbf{B}.

(d)

CTB\mathbf{C}^{T}\mathbf{B}.

(e)

A2\mathbf{A}^{2}.

(f)

B2\mathbf{B}^{2}.

(g)

CTC\mathbf{C}^{T}\mathbf{C}.

(h)

CCT\mathbf{C}\mathbf{C}^{T}.

(i)

BBT\mathbf{B}\mathbf{B}^{T}.

(j)

BTB\mathbf{B}^{T}\mathbf{B}.

(k)

CTAC\mathbf{C}^{T}\mathbf{A}\mathbf{C}.

Задача 3.5.2

Рассмотрим следующую систему уравнений:

2x1+x2+x3=3,4x1+2x3=10,2x1+2x2=−2. \begin{aligned} 2x_{1} + x_{2} + x_{3} & = 3, \\ 4x_{1} + 2x_{3} & = 10, \\ 2x_{1} + 2x_{2} & = -2. \end{aligned}
?
(a)

Запишите систему в виде матричного уравнения вида Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}.

(b)

Запишите решение системы в виде столбца s→\overrightarrow {s} и проверьте умножением матриц, что s→\overrightarrow {s} удовлетворяет уравнению Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}.

(c)

Запишите b→\overrightarrow {b} как линейную комбинацию столбцов матрицы A\mathbf{A}.

Задача 3.5.3

Пусть E=[100010301]\mathbf{E} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 3 & 0 & 1 \end{bmatrix}, и пусть A\mathbf{A} — произвольная матрица размера 3×33 \times 3.

?
(a)

Опишите строки матрицы EA\mathbf{E}\mathbf{A} через строки матрицы A\mathbf{A}.

(b)

Опишите столбцы матрицы AE\mathbf{A}\mathbf{E} через столбцы матрицы A\mathbf{A}.

Задача 3.5.4

Пусть e→j\overrightarrow {e}_{j} обозначает jj-й единичный столбец, содержащий 11 на jj-м месте и нули на всех остальных местах. Для общей матрицы An×n\mathbf{A}_{n \times n} опишите следующие произведения.

?
(a)

Ae→j\mathbf{A}\overrightarrow {e}_{j}.

(b)

e→iTA\overrightarrow {e}_{i}^{T}\mathbf{A}.

(c)

e→iTAe→j\overrightarrow {e}_{i}^{T}\mathbf{A}\overrightarrow {e}_{j}.

Задача 3.5.5

Пусть A\mathbf{A} и B\mathbf{B} — матрицы размера m×nm \times n. Если Ax→=Bx→\mathbf{A}\overrightarrow {x} = \mathbf{B}\overrightarrow {x} выполняется для всех столбцов x→\overrightarrow {x} размера n×1n \times 1, докажите, что A=B\mathbf{A} = \mathbf{B}.

?
Задача 3.5.6

Для A=[1/2α01/2]\mathbf{A} = \begin{bmatrix} 1/2 & \alpha \\ 0 & 1/2 \end{bmatrix} найдите lim⁡n→∞An\lim_{n \to \infty } \mathbf{A}^{n}.

?
Задача 3.5.7

Если Cm×1\mathbf{C}_{m \times 1} и R1×n\mathbf{R}_{1 \times n} — матрицы, состоящие соответственно из одного столбца и одной строки, то матричное произведение Pm×n=CR\mathbf{P}_{m \times n} = \mathbf{C}\mathbf{R} иногда называют внешним произведением C\mathbf{C} на R\mathbf{R}. Для согласованных матриц A\mathbf{A} и B\mathbf{B} объясните, как записать произведение AB\mathbf{A}\mathbf{B} в виде суммы внешних произведений, составленных из столбцов A\mathbf{A} и строк B\mathbf{B}.

?
Задача 3.5.8

Квадратная матрица U=[uij]\mathbf{U} = [u_{ij}] называется верхнетреугольной, если uij=0u_{ij} = 0 при i>ji > j — то есть все элементы ниже главной диагонали равны 00.

?
(a)

Если A\mathbf{A} и B\mathbf{B} — две верхнетреугольные матрицы размера n×nn \times n, объясните, почему произведение AB\mathbf{A}\mathbf{B} также обязано быть верхнетреугольным.

(b)

Если An×n\mathbf{A}_{n \times n} и Bn×n\mathbf{B}_{n \times n} верхнетреугольны, каковы диагональные элементы AB\mathbf{A}\mathbf{B}?

(c)

Матрица L\mathbf{L} называется нижнетреугольной, если ℓij=0\ell_{ij} = 0 при i<ji < j. Верно ли, что произведение двух нижнетреугольных матриц размера n×nn \times n снова является нижнетреугольным?

Задача 3.5.9

Если A=[aij(t)]\mathbf{A} = [a_{ij}(t)] — матрица, элементы которой являются функциями переменной tt, то производная матрицы A\mathbf{A} по tt определяется как матрица производных. То есть,

dAdt=[daijdt]. \frac{d\mathbf{A}}{dt} = \left[\frac{da_{ij}}{dt}\right].

Выведите правило дифференцирования произведения

d(AB)dt=dAdtB+AdBdt. \frac{d(\mathbf{A}\mathbf{B})}{dt} = \frac{d\mathbf{A}}{dt}\mathbf{B} + \mathbf{A}\frac{d\mathbf{B}}{dt}.
?
Задача 3.5.10

Вспомните построение из примера 3.5.2: для сети из nn узлов (например, городов, обслуживаемых авиакомпанией) матрица связности Cn×n=[cij]\mathbf{C}_{n \times n} = [c_{ij}] определяется как cij=1c_{ij} = 1, если существует прямой путь (рейс) из узла ii в узел jj, и cij=0c_{ij} = 0 в противном случае. Пусть Cn×n\mathbf{C}_{n \times n} — матрица связности, связанная с такой сетью из nn узлов, и пусть e→\overrightarrow {e} — столбец из всех единиц размера n×1n \times 1. В терминах сети опишите элементы каждого из следующих произведений.

?
(a)

Истолкуйте произведение Ce→\mathbf{C}\overrightarrow {e}.

(b)

Истолкуйте произведение e→TC\overrightarrow {e}^{T}\mathbf{C}.

Задача 3.5.11

Рассмотрим три резервуара, каждый из которых содержит VV галлонов рассола. Резервуары соединены в линию, и все краны открываются одновременно. Пресная вода поступает со скоростью rr галлонов/сек в верхнюю часть первого резервуара, rr галлонов/сек вытекает из нижней части и поступает в следующий резервуар, и так далее по цепочке — в каждый резервуар сверху поступает и снизу вытекает rr галлонов/сек.

Figure 3.5.2Figure 3.5.2

Пусть xi(t)x_{i}(t) обозначает количество фунтов соли в резервуаре ii в момент времени tt, и пусть

x→=(x1(t)x2(t)x3(t))иdx→dt=(dx1/dtdx2/dtdx3/dt). \overrightarrow {x} = \begin{pmatrix} x_{1}(t) \\ x_{2}(t) \\ x_{3}(t) \end{pmatrix} \qquad \text{и} \qquad \frac{d\overrightarrow {x}}{dt} = \begin{pmatrix} dx_{1}/dt \\ dx_{2}/dt \\ dx_{3}/dt \end{pmatrix}.

Предполагая, что в каждом резервуаре происходит полное непрерывное перемешивание, покажите, что

dx→dt=Ax→,гдеA=rV[−1001−1001−1]. \frac{d\overrightarrow {x}}{dt} = \mathbf{A}\overrightarrow {x}, \qquad \text{где} \qquad \mathbf{A} = \frac{r}{V}\begin{bmatrix} -1 & 0 & 0 \\ 1 & -1 & 0 \\ 0 & 1 & -1 \end{bmatrix}.
?
Примечание.
?

Используйте тот факт, что dxidt=скорость изменения=фунтовсек поступающих−фунтовсек убывающих\dfrac {dx_{i}}{dt} = \text{скорость изменения} = \dfrac {\text{фунтов}}{\text{сек}} \text{ поступающих} - \dfrac {\text{фунтов}}{\text{сек}} \text{ убывающих}.

§
Задача 3.6.1

Для блочных матриц

A=[100333100333122000]иB=[−1−10000−1−2−1−2−1−2], \mathbf{A} = \left[\begin{array}{c|cc|ccc} 1 & 0 & 0 & 3 & 3 & 3 \\ 1 & 0 & 0 & 3 & 3 & 3 \\ \hline 1 & 2 & 2 & 0 & 0 & 0 \end{array}\right] \qquad \text{и} \qquad \mathbf{B} = \left[\begin{array}{cc} -1 & -1 \\ \hline 0 & 0 \\ 0 & 0 \\ \hline -1 & -2 \\ -1 & -2 \\ -1 & -2 \end{array}\right],

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

?
Задача 3.6.2

Для всех матриц An×k\mathbf{A}_{n \times k} и Bk×n\mathbf{B}_{k \times n} покажите, что блочная матрица

L=(I−BAB2A−ABAAB−I) \mathbf{L} = \begin{pmatrix} \mathbf{I}-\mathbf{B}\mathbf{A} & \mathbf{B} \\ 2\mathbf{A}-\mathbf{A}\mathbf{B}\mathbf{A} & \mathbf{A}\mathbf{B}-\mathbf{I} \end{pmatrix}

обладает свойством L2=I\mathbf{L}^{2} = \mathbf{I}. Матрицы с таким свойством называются инволютивными, и они встречаются в криптографии.

?
Задача 3.6.3

Для матрицы

A=[1001/31/31/30101/31/31/30011/31/31/30001/31/31/30001/31/31/30001/31/31/3], \mathbf{A} = \begin{bmatrix} 1 & 0 & 0 & 1/3 & 1/3 & 1/3 \\ 0 & 1 & 0 & 1/3 & 1/3 & 1/3 \\ 0 & 0 & 1 & 1/3 & 1/3 & 1/3 \\ 0 & 0 & 0 & 1/3 & 1/3 & 1/3 \\ 0 & 0 & 0 & 1/3 & 1/3 & 1/3 \\ 0 & 0 & 0 & 1/3 & 1/3 & 1/3 \end{bmatrix},

найдите A300\mathbf{A}^{300}.

?
Задача 3.6.4

Для произвольной матрицы Am×n\mathbf{A}_{m \times n} покажите, что произведения A∗A\mathbf{A}^{*}\mathbf{A} и AA∗\mathbf{A}\mathbf{A}^{*} являются эрмитовыми матрицами.

?
Задача 3.6.5

Если A\mathbf{A} и B\mathbf{B} — симметричные матрицы, которые коммутируют, докажите, что произведение AB\mathbf{A}\mathbf{B} также симметрично. Если AB≠BA\mathbf{A}\mathbf{B} \neq \mathbf{B}\mathbf{A}, обязательно ли AB\mathbf{A}\mathbf{B} симметрично?

?
Задача 3.6.6

Докажите, что справедлив правый дистрибутивный закон: для согласованных матриц D\mathbf{D}, E\mathbf{E}, F\mathbf{F} выполняется (D+E)F=DF+EF(\mathbf{D}+\mathbf{E})\mathbf{F} = \mathbf{D}\mathbf{F}+\mathbf{E}\mathbf{F}.

?
Задача 3.6.7

Для каждой матрицы An×n\mathbf{A}_{n \times n} объясните, почему невозможно найти решение для Xn×n\mathbf{X}_{n \times n} в матричном уравнении

AX−XA=I. \mathbf{A}\mathbf{X}-\mathbf{X}\mathbf{A} = \mathbf{I}.
?
Задача 3.6.8

Пусть y→1×mT\overrightarrow {y}^{T}_{1 \times m} — строка неизвестных, а Am×n\mathbf{A}_{m \times n} и b→1×nT\overrightarrow {b}^{T}_{1 \times n} — известные матрицы.

?
(a)

Объясните, почему матричное уравнение y→TA=b→T\overrightarrow {y}^{T}\mathbf{A} = \overrightarrow {b}^{T} представляет собой систему из nn линейных уравнений с mm неизвестными.

(b)

Как решения для y→T\overrightarrow {y}^{T} в уравнении y→TA=b→T\overrightarrow {y}^{T}\mathbf{A} = \overrightarrow {b}^{T} связаны с решениями для x→\overrightarrow {x} в уравнении ATx→=b→\mathbf{A}^{T}\overrightarrow {x} = \overrightarrow {b}?

Задача 3.6.9

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

?
(a)

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

(b)

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

Задача 3.6.10

Вспомните приём приведения к блочно-треугольному виду из примера 3.6.7: если система Tn×nx→=b→\mathbf{T}_{n \times n}\overrightarrow {x} = \overrightarrow {b} имеет блочно-треугольную матрицу коэффициентов T=(AB0C)\mathbf{T} = \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{0} & \mathbf{C} \end{pmatrix} (где A\mathbf{A} имеет размер r×rr \times r, а C\mathbf{C} — размер (n−r)×(n−r)(n-r) \times (n-r)), а x→\overrightarrow {x}, b→\overrightarrow {b} разбиты на блоки соответствующим образом: x→=(x→1x→2)\overrightarrow {x} = \begin{pmatrix} \overrightarrow {x}_{1} \\ \overrightarrow {x}_{2} \end{pmatrix}, b→=(b→1b→2)\overrightarrow {b} = \begin{pmatrix} \overrightarrow {b}_{1} \\ \overrightarrow {b}_{2} \end{pmatrix}, то блочное умножение сводит Tx→=b→\mathbf{T}\overrightarrow {x} = \overrightarrow {b} к двум меньшим системам: Ax→1+Bx→2=b→1\mathbf{A}\overrightarrow {x}_{1}+\mathbf{B}\overrightarrow {x}_{2} = \overrightarrow {b}_{1} и Cx→2=b→2\mathbf{C}\overrightarrow {x}_{2} = \overrightarrow {b}_{2}: решите вторую относительно x→2\overrightarrow {x}_{2}, затем подставьте в первую и решите относительно x→1\overrightarrow {x}_{1}.

Запишите следующую систему в виде Tn×nx→=b→\mathbf{T}_{n \times n}\overrightarrow {x} = \overrightarrow {b}, где T\mathbf{T} блочно-треугольна, а затем получите решение, решив две меньшие системы, как описано выше.

x1+x2+3x3+4x4=−1,2x3+3x4=3,x1+2x2+5x3+6x4=−2,x3+2x4=4. \begin{aligned} x_{1} + x_{2} + 3x_{3} + 4x_{4} & = -1, \\ 2x_{3} + 3x_{4} & = 3, \\ x_{1} + 2x_{2} + 5x_{3} + 6x_{4} & = -2, \\ x_{3} + 2x_{4} & = 4. \end{aligned}
?
Задача 3.6.11

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

?
(a)

\trace(ABC)=\trace(BCA)=\trace(CAB)\trace (\mathbf{A}\mathbf{B}\mathbf{C}) = \trace (\mathbf{B}\mathbf{C}\mathbf{A}) = \trace (\mathbf{C}\mathbf{A}\mathbf{B}).

(b)

\trace(ABC)\trace (\mathbf{A}\mathbf{B}\mathbf{C}) может отличаться от \trace(BAC)\trace (\mathbf{B}\mathbf{A}\mathbf{C}).

(c)

\trace(ATB)=\trace(ABT)\trace (\mathbf{A}^{T}\mathbf{B}) = \trace (\mathbf{A}\mathbf{B}^{T}).

Задача 3.6.12

Предположим, что Am×n\mathbf{A}_{m \times n} и x→n×1\overrightarrow {x}_{n \times 1} имеют вещественные элементы.

?
(a)

Докажите, что x→Tx→=0\overrightarrow {x}^{T}\overrightarrow {x} = 0 тогда и только тогда, когда x→=0→\overrightarrow {x} = \overrightarrow {0}.

(b)

Докажите, что \trace(ATA)=0\trace (\mathbf{A}^{T}\mathbf{A}) = 0 тогда и только тогда, когда A=0\mathbf{A} = \mathbf{0}.

§
Задача 3.7.1

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

?
(a)

[1213]\begin{bmatrix} 1 & 2 \\ 1 & 3 \end{bmatrix}.

(b)

[1224]\begin{bmatrix} 1 & 2 \\ 2 & 4 \end{bmatrix}.

(c)

[4−854−743−42]\begin{bmatrix} 4 & -8 & 5 \\ 4 & -7 & 4 \\ 3 & -4 & 2 \end{bmatrix}.

(d)

[123456789]\begin{bmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \end{bmatrix}.

(e)

[1111122212331234]\begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 2 & 2 & 2 \\ 1 & 2 & 3 & 3 \\ 1 & 2 & 3 & 4 \end{bmatrix}.

Задача 3.7.2

Найдите матрицу X\mathbf{X}, такую что X=AX+B\mathbf{X} = \mathbf{A}\mathbf{X}+\mathbf{B}, где

A=[0−1000−1000]иB=[122133]. \mathbf{A} = \begin{bmatrix} 0 & -1 & 0 \\ 0 & 0 & -1 \\ 0 & 0 & 0 \end{bmatrix} \qquad \text{и} \qquad \mathbf{B} = \begin{bmatrix} 1 & 2 \\ 2 & 1 \\ 3 & 3 \end{bmatrix}.
?
Задача 3.7.3

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

?
(a)

Если A\mathbf{A} содержит нулевую строку или нулевой столбец, то A\mathbf{A} вырождена.

(b)

Если A\mathbf{A} содержит две одинаковые строки или два одинаковых столбца, то A\mathbf{A} вырождена.

(c)

Если одна строка (или столбец) является кратной другой строке (или столбцу), то A\mathbf{A} обязательно вырождена.

Задача 3.7.4

Ответьте на каждый из следующих вопросов.

?
(a)

При каких условиях диагональная матрица невырождена? Опишите структуру обратной к диагональной матрице.

(b)

При каких условиях треугольная матрица невырождена? Опишите структуру обратной к треугольной матрице.

Задача 3.7.5

Если A\mathbf{A} невырождена и симметрична, докажите, что A−1\mathbf{A}^{-1} симметрична.

?
Задача 3.7.6

Если A\mathbf{A} — квадратная матрица, такая что I−A\mathbf{I}-\mathbf{A} невырождена, докажите, что

A(I−A)−1=(I−A)−1A. \mathbf{A}(\mathbf{I}-\mathbf{A})^{-1} = (\mathbf{I}-\mathbf{A})^{-1}\mathbf{A}.
?
Задача 3.7.7

Докажите, что если A\mathbf{A} имеет размер m×nm \times n, а B\mathbf{B} — размер n×mn \times m, причём AB=Im\mathbf{A}\mathbf{B} = \mathbf{I}_{m} и BA=In\mathbf{B}\mathbf{A} = \mathbf{I}_{n}, то m=nm = n.

?
Задача 3.7.8

Если A\mathbf{A}, B\mathbf{B} и A+B\mathbf{A}+\mathbf{B} каждая невырождена, докажите, что

A(A+B)−1B=B(A+B)−1A=(A−1+B−1)−1. \mathbf{A}(\mathbf{A}+\mathbf{B})^{-1}\mathbf{B} = \mathbf{B}(\mathbf{A}+\mathbf{B})^{-1}\mathbf{A} = (\mathbf{A}^{-1}+\mathbf{B}^{-1})^{-1}.
?
Задача 3.7.9

Пусть S\mathbf{S} — кососимметрическая матрица с вещественными элементами.

?
(a)

Докажите, что I−S\mathbf{I}-\mathbf{S} невырождена.

(b)

Если A=(I+S)(I−S)−1\mathbf{A} = (\mathbf{I}+\mathbf{S})(\mathbf{I}-\mathbf{S})^{-1}, покажите, что A−1=AT\mathbf{A}^{-1} = \mathbf{A}^{T}.

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

Для пункта (а): x→Tx→=0  ⟹  x→=0→\overrightarrow {x}^{T}\overrightarrow {x} = 0 \implies \overrightarrow {x} = \overrightarrow {0}.

Задача 3.7.10

Для матриц Ar×r\mathbf{A}_{r \times r}, Bs×s\mathbf{B}_{s \times s} и Cr×s\mathbf{C}_{r \times s}, таких что A\mathbf{A} и B\mathbf{B} невырождены, проверьте, что каждое из следующих утверждений верно.

?
(a)

(A00B)−1=(A−100B−1)\begin{pmatrix} \mathbf{A} & \mathbf{0} \\ \mathbf{0} & \mathbf{B} \end{pmatrix}^{-1} = \begin{pmatrix} \mathbf{A}^{-1} & \mathbf{0} \\ \mathbf{0} & \mathbf{B}^{-1} \end{pmatrix}.

(b)

(AC0B)−1=(A−1−A−1CB−10B−1)\begin{pmatrix} \mathbf{A} & \mathbf{C} \\ \mathbf{0} & \mathbf{B} \end{pmatrix}^{-1} = \begin{pmatrix} \mathbf{A}^{-1} & -\mathbf{A}^{-1}\mathbf{C}\mathbf{B}^{-1} \\ \mathbf{0} & \mathbf{B}^{-1} \end{pmatrix}.

Задача 3.7.11

Рассмотрим блочную матрицу (Ar×rCr×sRs×rBs×s)\begin{pmatrix} \mathbf{A}_{r \times r} & \mathbf{C}_{r \times s} \\ \mathbf{R}_{s \times r} & \mathbf{B}_{s \times s} \end{pmatrix}. Когда указанные обратные матрицы существуют, матрицы, определённые формулами

S=B−RA−1CиT=A−CB−1R \mathbf{S} = \mathbf{B}-\mathbf{R}\mathbf{A}^{-1}\mathbf{C} \qquad \text{и} \qquad \mathbf{T} = \mathbf{A}-\mathbf{C}\mathbf{B}^{-1}\mathbf{R}

называются дополнениями Шура матриц A\mathbf{A} и B\mathbf{B} соответственно.

?
(a)

Если A\mathbf{A} и S\mathbf{S} обе невырождены, проверьте, что

(ACRB)−1=(A−1+A−1CS−1RA−1−A−1CS−1−S−1RA−1S−1). \begin{pmatrix} \mathbf{A} & \mathbf{C} \\ \mathbf{R} & \mathbf{B} \end{pmatrix}^{-1} = \begin{pmatrix} \mathbf{A}^{-1}+\mathbf{A}^{-1}\mathbf{C}\mathbf{S}^{-1}\mathbf{R}\mathbf{A}^{-1} & -\mathbf{A}^{-1}\mathbf{C}\mathbf{S}^{-1} \\ -\mathbf{S}^{-1}\mathbf{R}\mathbf{A}^{-1} & \mathbf{S}^{-1} \end{pmatrix}.
(b)

Если B\mathbf{B} и T\mathbf{T} невырождены, проверьте, что

(ACRB)−1=(T−1−T−1CB−1−B−1RT−1B−1+B−1RT−1CB−1). \begin{pmatrix} \mathbf{A} & \mathbf{C} \\ \mathbf{R} & \mathbf{B} \end{pmatrix}^{-1} = \begin{pmatrix} \mathbf{T}^{-1} & -\mathbf{T}^{-1}\mathbf{C}\mathbf{B}^{-1} \\ -\mathbf{B}^{-1}\mathbf{R}\mathbf{T}^{-1} & \mathbf{B}^{-1}+\mathbf{B}^{-1}\mathbf{R}\mathbf{T}^{-1}\mathbf{C}\mathbf{B}^{-1} \end{pmatrix}.
Задача 3.7.12

Предположим, что A\mathbf{A}, B\mathbf{B}, C\mathbf{C} и D\mathbf{D} — матрицы размера n×nn \times n, такие что ABT\mathbf{A}\mathbf{B}^{T} и CDT\mathbf{C}\mathbf{D}^{T} каждая симметрична, а ADT−BCT=I\mathbf{A}\mathbf{D}^{T}-\mathbf{B}\mathbf{C}^{T} = \mathbf{I}. Докажите, что

ATD−CTB=I. \mathbf{A}^{T}\mathbf{D}-\mathbf{C}^{T}\mathbf{B} = \mathbf{I}.
?
§
Задача 3.8.1

Вспомните формулу Шермана–Моррисона: если An×n\mathbf{A}_{n \times n} невырождена, а c→\overrightarrow {c}, d→\overrightarrow {d} — столбцы размера n×1n \times 1, такие что 1+d→TA−1c→≠01+\overrightarrow {d}^{T}\mathbf{A}^{-1}\overrightarrow {c} \neq 0, то сумма A+c→d→T\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T} невырождена, и

(A+c→d→T)−1=A−1−A−1c→d→TA−11+d→TA−1c→. (\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T})^{-1} = \mathbf{A}^{-1} - \frac{\mathbf{A}^{-1}\overrightarrow {c}\overrightarrow {d}^{T}\mathbf{A}^{-1}}{1+\overrightarrow {d}^{T}\mathbf{A}^{-1}\overrightarrow {c}}.

Предположим, вам даны

A=[20−1−111−101]иA−1=[10101−1102]. \mathbf{A} = \begin{bmatrix} 2 & 0 & -1 \\ -1 & 1 & 1 \\ -1 & 0 & 1 \end{bmatrix} \qquad \text{и} \qquad \mathbf{A}^{-1} = \begin{bmatrix} 1 & 0 & 1 \\ 0 & 1 & -1 \\ 1 & 0 & 2 \end{bmatrix}.
?
(a)

Используя формулу Шермана--Моррисона, найдите обратную матрицу для B\mathbf{B}, полученной из A\mathbf{A} изменением элемента (3,2)(3,2) с 00 на 22.

(b)

Пусть C\mathbf{C} — матрица, совпадающая с A\mathbf{A} всюду, кроме того что c32=2c_{32}=2 и c33=2c_{33}=2. Используя формулу Шермана--Моррисона, найдите C−1\mathbf{C}^{-1}.

Задача 3.8.2

Предположим, что A\mathbf{A} и B\mathbf{B} — невырожденные матрицы, причём B\mathbf{B} получается из A\mathbf{A} заменой A∗j\mathbf{A}_{*j} на другой столбец b→\overrightarrow {b}. Используя формулу Шермана--Моррисона, выведите, что

B−1=A−1−(A−1b→−e→j)[A−1]j∗[A−1]j∗b→. \mathbf{B}^{-1} = \mathbf{A}^{-1} - \frac{(\mathbf{A}^{-1}\overrightarrow {b}-\overrightarrow {e}_{j})[\mathbf{A}^{-1}]_{j*}}{[\mathbf{A}^{-1}]_{j*}\overrightarrow {b}}.
?
Задача 3.8.3

Предположим, что матрица коэффициентов невырожденной системы Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b} обновляется, порождая другую невырожденную систему (A+c→d→T)z→=b→(\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T})\overrightarrow {z} = \overrightarrow {b}, где b→,c→,d→∈Rn×1\overrightarrow {b}, \overrightarrow {c}, \overrightarrow {d} \in \mathbb {R}^{n \times 1}, и пусть y→\overrightarrow {y} — решение Ay→=c→\mathbf{A}\overrightarrow {y} = \overrightarrow {c}. Покажите, что

z→=x→−y→d→Tx→1+d→Ty→. \overrightarrow {z} = \overrightarrow {x} - \frac{\overrightarrow {y}\overrightarrow {d}^{T}\overrightarrow {x}}{1+\overrightarrow {d}^{T}\overrightarrow {y}}.
?
Задача 3.8.4
?
(a)

Используя формулу Шермана--Моррисона, докажите, что если A\mathbf{A} невырождена, то A+αe→ie→jT\mathbf{A}+\alpha \overrightarrow {e}_{i}\overrightarrow {e}_{j}^{T} невырождена при достаточно малом α\alpha.

(b)

Используя пункт (а), докажите, что I+E\mathbf{I}+\mathbf{E} невырождена, когда все ϵij\epsilon_{ij} достаточно малы по модулю. Это альтернатива использованию рассуждения с рядом Неймана.

Задача 3.8.5

Для заданных матриц A\mathbf{A} и B\mathbf{B}, где A\mathbf{A} невырождена, объясните, почему A+ϵB\mathbf{A}+\epsilon \mathbf{B} также невырождена, когда вещественное число ϵ\epsilon ограничено достаточно малым интервалом вокруг начала координат. Иными словами, докажите, что малые возмущения невырожденных матриц также невырождены.

?
Задача 3.8.6

Выведите формулу Шермана--Моррисона--Вудбери: если An×n\mathbf{A}_{n \times n} невырождена, а C\mathbf{C}, D\mathbf{D} имеют размер n×kn \times k, причём (I+DTA−1C)−1(\mathbf{I}+\mathbf{D}^{T}\mathbf{A}^{-1}\mathbf{C})^{-1} существует, то

(A+CDT)−1=A−1−A−1C(I+DTA−1C)−1DTA−1. (\mathbf{A}+\mathbf{C}\mathbf{D}^{T})^{-1} = \mathbf{A}^{-1} - \mathbf{A}^{-1}\mathbf{C}(\mathbf{I}+\mathbf{D}^{T}\mathbf{A}^{-1}\mathbf{C})^{-1}\mathbf{D}^{T}\mathbf{A}^{-1}.
?
Примечание.
?

Вспомните упражнение 3.7.11 (формулы блочного обращения через дополнения Шура) и рассмотрите произведение (IC0I)(ACDT−I)(I0DTI)\begin{pmatrix} \mathbf{I} & \mathbf{C} \\ \mathbf{0} & \mathbf{I} \end{pmatrix}\begin{pmatrix} \mathbf{A} & \mathbf{C} \\ \mathbf{D}^{T} & -\mathbf{I} \end{pmatrix}\begin{pmatrix} \mathbf{I} & \mathbf{0} \\ \mathbf{D}^{T} & \mathbf{I} \end{pmatrix}.

Задача 3.8.7

Вспомните норму ∥A∥=max⁡i∑j∣aij∣\left\| \mathbf{A}\right\| = \max_{i}\sum_{j}\left|a_{ij}\right| (максимальную сумму модулей элементов строки) и связанное с ней число обусловленности κ(A)=∥A∥∥A−1∥\kappa (\mathbf{A}) = \left\| \mathbf{A}\right\| \left\| \mathbf{A}^{-1}\right\| невырожденной матрицы A\mathbf{A}, которое показывает, насколько плохо обусловлена A\mathbf{A} (чем больше κ\kappa, тем хуже обусловленность).

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

A=[1000−1000100−100−100−100300],B=[18−1−9−711111718],C=[122−4201−45−45−9481]. \mathbf{A} = \begin{bmatrix} 100 & 0 & -100 \\ 0 & 100 & -100 \\ -100 & -100 & 300 \end{bmatrix}, \qquad \mathbf{B} = \begin{bmatrix} 1 & 8 & -1 \\ -9 & -71 & 11 \\ 1 & 17 & 18 \end{bmatrix}, \qquad \mathbf{C} = \begin{bmatrix} 1 & 22 & -42 \\ 0 & 1 & -45 \\ -45 & -948 & 1 \end{bmatrix}.
?
Задача 3.8.8

Предположим, что элементы A(t)\mathbf{A}(t), x→(t)\overrightarrow {x}(t) и b→(t)\overrightarrow {b}(t) являются дифференцируемыми функциями вещественной переменной tt, такими что A(t)x→(t)=b→(t)\mathbf{A}(t)\overrightarrow {x}(t) = \overrightarrow {b}(t).

?
(a)

Предполагая, что A(t)−1\mathbf{A}(t)^{-1} существует, объясните, почему

dA(t)−1dt=−A(t)−1A′(t)A(t)−1. \frac{d\mathbf{A}(t)^{-1}}{dt} = -\mathbf{A}(t)^{-1}\mathbf{A}'(t)\mathbf{A}(t)^{-1}.
(b)

Выведите уравнение

x→′(t)=A(t)−1b→′(t)−A(t)−1A′(t)x→(t). \overrightarrow {x}'(t) = \mathbf{A}(t)^{-1}\overrightarrow {b}'(t) - \mathbf{A}(t)^{-1}\mathbf{A}'(t)\overrightarrow {x}(t).

Это показывает, что A−1\mathbf{A}^{-1} усиливает как изменение A\mathbf{A}, так и изменение b→\overrightarrow {b}, и тем самым подтверждает наблюдение, сделанное в задаче (3.8.8), о том, что чувствительность невырожденной системы к малым возмущениям напрямую связана с величиной элементов A−1\mathbf{A}^{-1}.

§
Задача 3.9.1

Некоторые сведения, необходимые для этого упражнения: для матрицы A\mathbf{A} размера m×nm \times n через EA\mathbf{E}_{\mathbf{A}} обозначается её (единственная) приведённая ступенчатая форма по строкам. Всякий раз, когда B\mathbf{B} можно получить из A\mathbf{A} последовательностью только элементарных операций над строками, мы пишем A∼rowB\mathbf{A} \overset {\text{row}}{\sim } \mathbf{B}, что выполняется тогда и только тогда, когда PA=B\mathbf{P}\mathbf{A} = \mathbf{B} для некоторой невырожденной P\mathbf{P}; аналогично A∼colB\mathbf{A} \overset {\text{col}}{\sim } \mathbf{B} (эквивалентность по столбцам) выполняется тогда и только тогда, когда AQ=B\mathbf{A}\mathbf{Q} = \mathbf{B} для некоторой невырожденной Q\mathbf{Q}, а A∼B\mathbf{A} \sim \mathbf{B} (эквивалентность, допускающая как операции над строками, так и над столбцами) выполняется тогда и только тогда, когда PAQ=B\mathbf{P}\mathbf{A}\mathbf{Q} = \mathbf{B} для невырожденных P\mathbf{P}, Q\mathbf{Q}.

Предположим, что A\mathbf{A} — матрица размера m×nm \times n.

?
(a)

Если [A∣Im][\mathbf{A}\mid \mathbf{I}_{m}] приводится по строкам к матрице [B∣P][\mathbf{B}\mid \mathbf{P}], объясните, почему P\mathbf{P} обязательно является невырожденной матрицей, такой что PA=B\mathbf{P}\mathbf{A} = \mathbf{B}.

(b)

Если [AIn]\left[\begin{smallmatrix} \mathbf{A} \\ \mathbf{I}_{n}\end{smallmatrix}\right] приводится по столбцам к [CQ]\left[\begin{smallmatrix} \mathbf{C} \\ \mathbf{Q}\end{smallmatrix}\right], объясните, почему Q\mathbf{Q} обязательно является невырожденной матрицей, такой что AQ=C\mathbf{A}\mathbf{Q} = \mathbf{C}.

(c)

Найдите невырожденную матрицу P\mathbf{P}, такую что PA=EA\mathbf{P}\mathbf{A} = \mathbf{E}_{\mathbf{A}}, где

A=[123424671236]. \mathbf{A} = \begin{bmatrix} 1 & 2 & 3 & 4 \\ 2 & 4 & 6 & 7 \\ 1 & 2 & 3 & 6 \end{bmatrix}.
(d)

Найдите невырожденные матрицы P\mathbf{P} и Q\mathbf{Q}, такие что PAQ\mathbf{P}\mathbf{A}\mathbf{Q} имеет ранговую нормальную форму (т.е. вид Nr=[Ir000]\mathbf{N}_{r} = \left[\begin{smallmatrix} \mathbf{I}_{r} & \mathbf{0} \\ \mathbf{0} & \mathbf{0}\end{smallmatrix}\right], где r=rk⁡(A)r = \operatorname {rk}\left(\mathbf{A}\right)).

Задача 3.9.2

Рассмотрим две матрицы

A=[220−13−1400−883]иB=[2−682514−13−9123]. \mathbf{A} = \begin{bmatrix} 2 & 2 & 0 & -1 \\ 3 & -1 & 4 & 0 \\ 0 & -8 & 8 & 3 \end{bmatrix} \qquad \text{и} \qquad \mathbf{B} = \begin{bmatrix} 2 & -6 & 8 & 2 \\ 5 & 1 & 4 & -1 \\ 3 & -9 & 12 & 3 \end{bmatrix}.
?
(a)

Эквивалентны ли A\mathbf{A} и B\mathbf{B}?

(b)

Эквивалентны ли A\mathbf{A} и B\mathbf{B} по строкам?

(c)

Эквивалентны ли A\mathbf{A} и B\mathbf{B} по столбцам?

Задача 3.9.3

Если A∼rowB\mathbf{A} \overset {\text{row}}{\sim } \mathbf{B}, объясните, почему базисные столбцы в A\mathbf{A} занимают в точности те же позиции, что и базисные столбцы в B\mathbf{B}.

?
Задача 3.9.4

Произведение элементарных матриц перестановки -- т.е. элементарных матриц типа I -- называется матрицей перестановки. Если P\mathbf{P} — матрица перестановки, объясните, почему P−1=PT\mathbf{P}^{-1} = \mathbf{P}^{T}.

?
Задача 3.9.5

Если An×n\mathbf{A}_{n \times n} — невырожденная матрица, какие (если таковые есть) из следующих утверждений верны?

?
(a)

A∼A−1\mathbf{A} \sim \mathbf{A}^{-1}.

(b)

A∼rowA−1\mathbf{A} \overset {\text{row}}{\sim } \mathbf{A}^{-1}.

(c)

A∼colA−1\mathbf{A} \overset {\text{col}}{\sim } \mathbf{A}^{-1}.

(d)

A∼I\mathbf{A} \sim \mathbf{I}.

(e)

A∼rowI\mathbf{A} \overset {\text{row}}{\sim } \mathbf{I}.

(f)

A∼colI\mathbf{A} \overset {\text{col}}{\sim } \mathbf{I}.

Задача 3.9.6

Какие (если таковые есть) из следующих утверждений верны?

?
(a)

A∼B  ⟹  AT∼BT\mathbf{A} \sim \mathbf{B} \implies \mathbf{A}^{T} \sim \mathbf{B}^{T}.

(b)

A∼rowB  ⟹  AT∼colBT\mathbf{A} \overset {\text{row}}{\sim } \mathbf{B} \implies \mathbf{A}^{T} \overset {\text{col}}{\sim } \mathbf{B}^{T}.

(c)

A∼rowB  ⟹  AT∼colBT\mathbf{A} \overset {\text{row}}{\sim } \mathbf{B} \implies \mathbf{A}^{T} \overset {\text{col}}{\sim } \mathbf{B}^{T}.

(d)

A∼rowB  ⟹  A∼B\mathbf{A} \overset {\text{row}}{\sim } \mathbf{B} \implies \mathbf{A} \sim \mathbf{B}.

(e)

A∼colB  ⟹  A∼B\mathbf{A} \overset {\text{col}}{\sim } \mathbf{B} \implies \mathbf{A} \sim \mathbf{B}.

(f)

A∼B  ⟹  A∼rowB\mathbf{A} \sim \mathbf{B} \implies \mathbf{A} \overset {\text{row}}{\sim } \mathbf{B}.

Задача 3.9.7

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

?
Задача 3.9.8

Если rk⁡(Am×n)=r\operatorname {rk}\left(\mathbf{A}_{m \times n}\right) = r, покажите, что существуют матрицы Bm×r\mathbf{B}_{m \times r} и Cr×n\mathbf{C}_{r \times n}, такие что A=BC\mathbf{A} = \mathbf{B}\mathbf{C}, где rk⁡(B)=rk⁡(C)=r\operatorname {rk}\left(\mathbf{B}\right) = \operatorname {rk}\left(\mathbf{C}\right) = r. Такое разложение называется разложением полного ранга.

?
Задача 3.9.9

Докажите, что rk⁡(Am×n)=1\operatorname {rk}\left(\mathbf{A}_{m \times n}\right) = 1 тогда и только тогда, когда существуют ненулевые столбцы u→m×1\overrightarrow {u}_{m \times 1} и v→n×1\overrightarrow {v}_{n \times 1}, такие что

A=u→v→T. \mathbf{A} = \overrightarrow {u}\overrightarrow {v}^{T}.
?
Задача 3.9.10

Докажите, что если rk⁡(An×n)=1\operatorname {rk}\left(\mathbf{A}_{n \times n}\right) = 1, то A2=τA\mathbf{A}^{2} = \tau \mathbf{A}, где τ=\trace(A)\tau = \trace (\mathbf{A}).

?
§
Задача 3.10.1

Пусть A=[1454182631630]\mathbf{A} = \begin{bmatrix} 1 & 4 & 5 \\ 4 & 18 & 26 \\ 3 & 16 & 30 \end{bmatrix}.

?
(a)

Найдите LU-множители матрицы A\mathbf{A}.

(b)

Используя LU-множители, решите Ax→1=b→1\mathbf{A}\overrightarrow {x}_{1} = \overrightarrow {b}_{1}, а также Ax→2=b→2\mathbf{A}\overrightarrow {x}_{2} = \overrightarrow {b}_{2}, где

b→1=(60−6)иb→2=(6612). \overrightarrow {b}_{1} = \begin{pmatrix} 6 \\ 0 \\ -6 \end{pmatrix} \qquad \text{и} \qquad \overrightarrow {b}_{2} = \begin{pmatrix} 6 \\ 6 \\ 12 \end{pmatrix}.
(c)

Используя LU-множители, найдите A−1\mathbf{A}^{-1}.

Задача 3.10.2

Пусть A\mathbf{A} и b→\overrightarrow {b} — матрицы

A=[1241736−12323−3202−26]иb→=(17334). \mathbf{A} = \begin{bmatrix} 1 & 2 & 4 & 17 \\ 3 & 6 & -12 & 3 \\ 2 & 3 & -3 & 2 \\ 0 & 2 & -2 & 6 \end{bmatrix} \qquad \text{и} \qquad \overrightarrow {b} = \begin{pmatrix} 17 \\ 3 \\ 3 \\ 4 \end{pmatrix}.
?
(a)

Объясните, почему матрица A\mathbf{A} не имеет LU-разложения.

(b)

Используя частичный выбор ведущего элемента, найдите матрицу перестановки P\mathbf{P}, а также LU-множители, такие что PA=LU\mathbf{P}\mathbf{A} = \mathbf{L}\mathbf{U}.

(c)

Используя информацию, содержащуюся в P\mathbf{P}, L\mathbf{L} и U\mathbf{U}, решите Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}.

Задача 3.10.3

Найдите все значения ξ\xi, при которых A=[ξ201ξ101ξ]\mathbf{A} = \begin{bmatrix} \xi & 2 & 0 \\ 1 & \xi & 1 \\ 0 & 1 & \xi \end{bmatrix} не имеет LU-разложения.

?
Задача 3.10.4

Если A\mathbf{A} — невырожденная матрица, обладающая LU-разложением, докажите, что ведущий элемент, возникающий после (k+1)(k+1) шагов стандартного метода Гаусса с использованием только операций III типа, задаётся формулой

pk+1=ak+1,k+1−c→TAk−1b→, p_{k+1} = a_{k+1,k+1} - \overrightarrow {c}^{T}\mathbf{A}_{k}^{-1}\overrightarrow {b},

где Ak\mathbf{A}_{k} и Ak+1=(Akb→c→Tak+1,k+1)\mathbf{A}_{k+1} = \begin{pmatrix} \mathbf{A}_{k} & \overrightarrow {b} \\ \overrightarrow {c}^{T} & a_{k+1,k+1} \end{pmatrix} — главные ведущие подматрицы порядков kk и k+1k+1 соответственно. Используя это, выведите, что все ведущие элементы должны быть отличны от нуля, если для A\mathbf{A} существует LU-разложение.

?
Задача 3.10.5

Если A\mathbf{A} — матрица, содержащая только целочисленные элементы, и все её ведущие элементы равны 11, объясните, почему A−1\mathbf{A}^{-1} также должна быть целочисленной матрицей.

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

Этот факт можно использовать для построения случайных целочисленных матриц, обладающих целочисленными обратными, случайным образом порождая целочисленные матрицы L\mathbf{L} и U\mathbf{U} с единичными диагоналями и затем составляя произведение A=LU\mathbf{A} = \mathbf{L}\mathbf{U}.

Задача 3.10.6

Рассмотрим трёхдиагональную матрицу T=[β1γ100α1β2γ200α2β3γ300α3β4]\mathbf{T} = \begin{bmatrix} \beta_{1} & \gamma_{1} & 0 & 0 \\ \alpha_{1} & \beta_{2} & \gamma_{2} & 0 \\ 0 & \alpha_{2} & \beta_{3} & \gamma_{3} \\ 0 & 0 & \alpha_{3} & \beta_{4} \end{bmatrix}.

?
(a)

Считая, что T\mathbf{T} обладает LU-разложением, убедитесь, что оно задаётся формулами

L=[1000α1/π11000α2/π21000α3/π31],U=[π1γ1000π2γ2000π3γ3000π4], \mathbf{L} = \begin{bmatrix} 1 & 0 & 0 & 0 \\ \alpha _{1}/\pi _{1} & 1 & 0 & 0 \\ 0 & \alpha _{2}/\pi _{2} & 1 & 0 \\ 0 & 0 & \alpha _{3}/\pi _{3} & 1 \end{bmatrix}, \qquad \mathbf{U} = \begin{bmatrix} \pi _{1} & \gamma _{1} & 0 & 0 \\ 0 & \pi _{2} & \gamma _{2} & 0 \\ 0 & 0 & \pi _{3} & \gamma _{3} \\ 0 & 0 & 0 & \pi _{4} \end{bmatrix},

где числа πi\pi_{i} порождаются рекуррентной формулой

π1=β1иπi+1=βi+1−αiγiπi. \pi _{1} = \beta _{1} \qquad \text{и} \qquad \pi _{i+1} = \beta _{i+1} - \frac{\alpha _{i}\gamma _{i}}{\pi _{i}}.

(Это верно для трёхдиагональных матриц произвольного размера, что делает вычисление LU-множителей таких матриц очень простым.)

(b)

Применяя приведённую выше рекуррентную формулу, найдите LU-разложение матрицы

T=[2−100−12−100−12−100−11]. \mathbf{T} = \begin{bmatrix} 2 & -1 & 0 & 0 \\ -1 & 2 & -1 & 0 \\ 0 & -1 & 2 & -1 \\ 0 & 0 & -1 & 1 \end{bmatrix}.
Задача 3.10.7

Матрица An×n\mathbf{A}_{n \times n} называется ленточной матрицей (band matrix), если aij=0a_{ij} = 0 при ∣i−j∣>w\left|i-j\right| > w для некоторого положительного целого числа ww, называемого шириной ленты (bandwidth). Иными словами, ненулевые элементы A\mathbf{A} ограничены полосой из ww диагональных линий выше и ниже главной диагонали. Например, трёхдиагональные матрицы имеют ширину ленты один, а диагональные матрицы имеют ширину ленты нуль. Если A\mathbf{A} — невырожденная матрица с шириной ленты ww и если A\mathbf{A} имеет LU-разложение A=LU\mathbf{A} = \mathbf{L}\mathbf{U}, то L\mathbf{L} наследует нижнюю ленточную структуру A\mathbf{A}, а U\mathbf{U} наследует верхнюю ленточную структуру в том смысле, что L\mathbf{L} имеет «нижнюю ширину ленты» ww, а U\mathbf{U} имеет «верхнюю ширину ленты» ww. Проиллюстрируйте, почему это так, используя типичную матрицу 5×55 \times 5 с шириной ленты w=2w = 2.

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

Постройте пример невырожденной симметричной матрицы, не обладающей LU- (или LDU-) разложением.

(b)

Постройте пример невырожденной симметричной матрицы, обладающей LU-разложением, но не являющейся положительно определённой.

Задача 3.10.9

Напомним про LDU-разложение: вынесение диагональных ведущих элементов из верхнетреугольного множителя U\mathbf{U} LU-разложения A=LU\mathbf{A} = \mathbf{L}\mathbf{U} даёт A=LDU′\mathbf{A} = \mathbf{L}\mathbf{D}\mathbf{U}', где D=\diag(u11,u22,…,unn)\mathbf{D} = \diag (u_{11}, u_{22}, \ldots , u_{nn}) — диагональная матрица ведущих элементов, а U′\mathbf{U}' — верхнетреугольная матрица с единицами на диагонали (переобозначая далее U′\mathbf{U}' как U\mathbf{U}, так что и L\mathbf{L}, и U\mathbf{U} в LDU-разложении имеют единичную диагональ).

?
(a)

Найдите LDU-множители для A=[1454182631630]\mathbf{A} = \begin{bmatrix} 1 & 4 & 5 \\ 4 & 18 & 26 \\ 3 & 16 & 30 \end{bmatrix} (это та же матрица, что использовалась в упражнении 3.10.1).

(b)

Докажите, что если матрица обладает LDU-разложением, то LDU-множители определены однозначно.

(c)

Если A\mathbf{A} симметрична и обладает LDU-разложением, объясните, почему оно обязательно имеет вид A=LDLT\mathbf{A} = \mathbf{L}\mathbf{D}\mathbf{L}^{T}.

Задача 3.10.10

Напомним: симметричная матрица A\mathbf{A}, обладающая LU-разложением, в котором все ведущие элементы положительны, называется положительно определённой, и такая матрица однозначно раскладывается как A=RTR\mathbf{A} = \mathbf{R}^{T}\mathbf{R}, где R\mathbf{R} — верхнетреугольная матрица с положительными диагональными элементами; это разложение Холецкого матрицы A\mathbf{A}, а R\mathbf{R} — множитель Холецкого матрицы A\mathbf{A}.

Объясните, почему A=[123281231227]\mathbf{A} = \begin{bmatrix} 1 & 2 & 3 \\ 2 & 8 & 12 \\ 3 & 12 & 27 \end{bmatrix} положительно определена, а затем найдите множитель Холецкого R\mathbf{R}.

?