3.5

Умножение матриц

[11/100%]
Показать
LaTeX
Задача 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{ убывающих}.