10 Умножение матриц
.pdfОпределитель полураспавшейся матрицы (3)
|
|
|
|
a32 |
: : : a3r |
|
f31 |
: : : f3q |
|
|
|
a32 |
: : : a3r |
f31 |
: : : f3q |
|
|||||||||||||
|
|
|
|
|
a22 |
: : : |
a2r |
|
f21 |
: : : |
f2q |
|
|
|
|
|
a12 |
: : : |
a1r |
f11 |
: : : |
f1q |
|
|
|||||
|
|
|
|
. |
. . |
. |
. |
. . |
. |
. . . |
. |
. . |
. |
. . . . |
. . . . |
|
|
|
. . . |
. . . . |
. . . . . |
. . . . |
. . . . |
. . . |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G |
= |
a11 |
|
|
|
|
|
: : : |
arr |
|
fr1 |
: : : |
frq |
|
|
a21 |
|
|
|
: : : |
arr fr1 |
: : : |
frq |
|
+ |
||||
|
|
ar2 |
|
|
|
|
|
|
|
|
ar2 |
|
|
|
|
||||||||||||||
j j |
|
|
|
|
0 |
|
: : : |
0 b11 |
: : : |
b1q |
|
|
|
0 |
: : : |
0 b11 |
: : : |
b1q |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
. . |
. |
. . |
. |
. |
. . . |
. |
. . |
. |
. . . . |
. . . . |
|
|
|
|
|
. . . . |
. . . . |
. . . . |
. . . . |
. . . . |
. . . |
|
|
|
|
|
|
|
|
0 |
|
: : : |
0 bq1 |
: : : |
bqq |
|
|
|
|
|
0 |
: : : |
0 bq1 |
: : : |
bqq |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a12 |
: : : |
|
a1r |
|
f11 : : : |
f1q |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a22 : : : |
|
a2r |
|
f21 : : : |
f2q |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
. |
. . . . |
. . . . |
. |
. . . |
. . . |
. |
. |
. . . . |
. . . . |
. . . . |
|
|
|
|
|
|
+ + ( |
|
|
)r |
+ |
1 |
|
|
|
|
|
: : : |
|
|
|
|
|
|
: : : |
|
|
|
|
|
|
||||
|
|
1 |
|
ar1 |
|
|
|
|
ar 1 r |
fr 1 1 |
|
= |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
ar 1 2 |
|
|
fr 1 q |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 |
: : : |
|
0 |
|
b11 |
: : : |
b1q |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
. |
. . . . |
. . . . |
. |
. . . |
. . . |
. |
. |
. . . . |
. . . . |
. . . . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
: : : |
|
0 |
|
bq1 |
: : : |
bqq |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=a11M11 jBj a21M21 jBj + + ( 1)r+1ar1Mr1 jBj =
=a11A11 jBj + a21A21 jBj + + ar1Ar1 jBj =
=(a11A11 + a21A21 + + ar1Ar1) jBj = jAj jBj:
Предложение доказано.
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (1)
Теперь мы можем доказать анонсированное выше утверждение об определителе произведения матриц.
Теорема 1
Если A = (aij ) и B = (bij ) квадратные матрицы одного и того же порядка, то jABj = jAj jBj.
Доказательство. Обозначим порядок матриц A и B через n, а матрицу ABчерез C. Пусть C = (cij ). Рассмотрим матрицу
|
0a21 |
a22 |
: : : a2n |
0 |
0 |
: : : |
0 1 |
|
||||||
|
|
a11 |
a12 |
: : : |
a1n |
0 |
0 |
: : : |
0 |
|
||||
|
Ba. . |
. . |
a. |
n.2. . |
:.:.:. .a. . |
. . . |
0. . . . |
0. . . |
:.:.:. . |
.0. . C |
: |
|||
|
= B n1 |
|
|
nn |
|
|
|
C |
||||||
D |
B |
|
1 0 |
: : : |
0 b11 |
|
: : : |
C |
|
|||||
|
B |
|
|
b12 |
b1nC |
|
||||||||
|
B |
|
|
|
1 |
: : : |
0 b21 |
|
: : : |
C |
|
|||
|
B |
0 |
|
|
b22 |
b2nC |
|
|||||||
|
B |
|
|
|
|
|
|
|
|
|
C |
|
||
|
B |
. . |
. . |
. |
. |
. . |
. . . . |
. . |
. . . |
. . . . |
. . . |
. . . . |
. . . C |
|
|
B |
0 0 |
: : : |
|
1 bn1 |
|
: : : |
C |
|
|||||
|
B |
|
|
bn2 |
bnnC |
|
||||||||
|
@ |
|
|
|
|
|
|
|
|
|
|
A |
|
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (2)
Матрица
|
|
0a12 |
a22 |
: : : |
an2 0 1 : : : |
0 1 |
|||||||
|
|
|
a11 |
a21 |
: : : |
an1 1 |
0 |
: : : |
0 |
||||
|
> |
Ba. |
1.n. . |
a. |
2.n. . |
:.:.:. . |
a. |
. . . |
.0. . . . |
0. . .:.:.:. . |
. .1. C |
||
|
= B |
|
|
|
|
nn |
|
|
C |
||||
D |
|
B |
0 |
0 |
: : : |
0 |
b11 b21 |
: : : |
C |
||||
|
|
B |
|
|
bn1C |
||||||||
|
|
B |
0 |
0 |
: : : |
0 |
b12 b22 |
: : : |
C |
||||
|
|
B |
|
|
bn2C |
||||||||
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
|
B |
. |
. . . |
. |
. . . |
. . . . |
. |
. . . |
. . . . . |
. . . |
. . . . |
. . . C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
|
@ |
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
0 |
0 : : : |
0 |
b1n b2n : : : |
bnn |
является полураспавшейся матрицей с диагональными блоками A и B.
Применяя предложение 9 из лекции 5 и предложение 1 из данной лекции, имеем jDj = jD>j = jA>j jB>j = jAj jBj, т. е.
jDj = jAj jBj: |
(2) |
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (3)
Прибавим к (n + 1)-му столбцу матрицы D ее первый столбец, умноженный на b11, второй, умноженный на b21, . . . , и, наконец, n-й, умноженный на bn1. Полученную матрицу обозначим через D0. Ясно, что во всех столбцах матрицы D0, кроме (n + 1)-го, стоят те же элементы, что и в матрице D. Если 1 6 i 6 n, то элемент, стоящий в i-й строке и
(n + 1)-м столбце матрицы D0, равен
ai1b11 + ai2b21 + + ainbn1 = ci1;
если же n + 1 6 i 6 2n, то этот элемент равен bi1 + bi1 = 0. Таким образом, матрица D0 имеет следующий вид:
|
|
0a21 |
a22 |
: : : a2n |
c21 |
0 |
: : : |
0 1 |
|
|||||||||
|
|
|
a11 |
a12 |
: : : |
a1n |
c11 |
0 |
: : : |
0 |
|
|||||||
|
0 |
Ba. |
n.1. . |
a. |
n.2. . |
:.:.:. .a. |
nn. . . |
c. |
. . . |
.0. . . |
:.:.:. . |
.0. .C |
: |
|||||
|
= B |
|
|
|
|
n1 |
|
C |
||||||||||
D |
|
B |
|
|
1 0 |
: : : |
0 0 |
b12 |
: : : |
C |
|
|||||||
|
|
B |
|
|
|
|
b1nC |
|
||||||||||
|
|
B |
|
|
|
1 |
: : : |
0 0 |
b22 |
: : : |
C |
|
||||||
|
|
B |
0 |
|
|
|
b2nC |
|
||||||||||
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
||
|
|
B |
. |
. |
. . |
. |
. |
. . |
. . . . |
. |
. |
. . |
. |
. . . |
. . . . |
. . . . |
. . .C |
|
|
|
B |
0 0 |
: : : |
|
|
1 0 |
|
: : : |
C |
|
|||||||
|
|
B |
|
|
bn2 |
bnnC |
|
|||||||||||
|
|
@ |
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (4)
Теперь прибавим к (n + 2)-му столбцу матрицы D0 ее первый столбец, умноженный на b12, второй, умноженный на b22, . . . , и, наконец, n-й, умноженный на bn2. Полученную матрицу обозначим через D00. Во всех столбцах матрицы D00, кроме (n + 2)-го, стоят те же элементы, что и в
матрице D0. Как и выше, проверяется, что элемент, стоящий в i-й строке и (n + 2)-м столбце матрицы D00, равен ci2, если 1 6 i 6 n, и 0, если
n + 1 6 i 6 2n. Таким образом,
|
|
0a21 |
a22 |
: : : a2n |
c21 |
c22 |
0 |
: : : |
0 1 |
|
|||||||||
|
|
|
a11 |
a12 |
: : : |
a1n |
c11 |
c12 |
0 |
: : : |
0 |
|
|||||||
|
00 |
Ba. . |
. |
.a. |
n.2. |
.:.:.:. .a. |
nn. . |
.c. |
n.1. |
.c. |
. . |
. .0. . |
.:.:.:. |
. .0. C |
: |
||||
|
= B n1 |
|
|
|
|
|
n2 |
|
C |
||||||||||
D |
|
B |
|
1 0 |
: : : |
0 0 |
0 |
b13 |
: : : |
C |
|
||||||||
|
|
B |
|
|
|
b1nC |
|
||||||||||||
|
|
B |
|
|
|
1 |
: : : |
0 0 |
0 |
b23 |
: : : |
C |
|
||||||
|
|
B |
0 |
|
|
|
b2nC |
|
|||||||||||
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
||
|
|
B |
. . |
. . |
. |
. |
. . |
. . . . |
. |
. |
. . |
. |
. . . |
. |
. . . |
. . . . |
. . . . |
. . C |
|
|
|
B |
0 0 |
: : : |
|
|
1 0 |
0 |
|
: : : |
C |
|
|||||||
|
|
B |
|
|
bn3 |
bnnC |
|
||||||||||||
|
|
@ |
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (5)
Изменим аналогичным образом (n + 3)-й столбец матрицы D00, затем (n + 4)-й столбец полученной матрицы и т. д. После того как будет изменен последний столбец, мы получим матрицу
|
|
0a21 |
a22 |
: : : a2n |
c21 |
c22 |
: : : c2n1 |
|
||||||||||||
|
|
|
a11 |
a12 |
: : : |
a1n |
c11 |
c12 |
: : : |
c1n |
|
|
||||||||
|
|
Ba. |
n.1. |
.a. |
n.2. |
.:.:.:. .a. |
. |
. .c. . . |
.c. |
n.2. |
.:.:.:. .c. |
. |
C |
: |
||||||
D |
= B |
|
: : : |
|
nn n1 |
|
: : : |
|
nnC |
|||||||||||
|
B |
|
|
1 |
0 |
0 |
0 |
0 |
0 |
C |
|
|||||||||
|
|
B |
|
|
: : : |
: : : |
C |
|
||||||||||||
|
|
B |
|
|
|
1 |
0 |
0 |
0 |
0 |
C |
|
||||||||
|
|
B |
0 |
|
|
|
C |
|
||||||||||||
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
||
|
|
B |
. |
. |
. . |
. |
. |
. . |
. . . . |
. |
. |
. . |
. . . . |
. |
. . . |
. . . . |
. |
. |
C |
|
|
|
B |
0 |
0 |
: : : |
|
|
1 |
0 |
0 |
: : : |
0 |
C |
|
||||||
|
|
B |
|
|
|
|
C |
|
||||||||||||
|
|
@ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
С учетом предложений 6 и 10 из лекции 5 получаем, что |
|
jDj = jD0j = jD00j = = jD j: |
|
Следовательно, |
|
jD j = jDj: |
(3) |
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (6)
Поменяем в матрице D местами сначала (n + 1)-й столбец с первым, затем (n + 2)-й столбец со вторым, . . . , наконец, последний столбец с n-м. В результате мы получим матрицу
|
0c21 |
c22 |
: : : c2n |
a21 |
a22 |
: : : |
a2n1 |
|
|||||||||
|
|
c11 |
c12 |
: : : |
c1n |
a11 |
a12 |
: : : |
a1n |
|
|
||||||
|
Bc. |
n.1. . |
c. |
n.2. |
.:.:.:. .c. . . .a. |
. |
. .a. |
. . |
.:.:.:. |
.a. |
. |
C |
: |
||||
|
|
|
|
||||||||||||||
D |
= B |
|
: : : |
nn |
n1 n2 |
: : : |
|
nnC |
|||||||||
B |
0 |
0 |
0 |
|
|
1 |
0 |
0 |
C |
|
|||||||
|
B |
: : : |
|
: : : |
C |
|
|||||||||||
|
B |
0 |
0 |
0 |
|
|
1 |
0 |
C |
|
|||||||
|
B |
|
0 |
|
|
C |
|
||||||||||
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
B |
. |
. . . |
. |
. . . |
. . . . |
. . . . |
. |
. |
. . |
. . |
. . |
. . . . |
. |
. |
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
@ |
0 |
0 |
: : : |
0 0 0 : : : 1 |
A |
|
||||||||||
|
|
|
|
Переходя от матрицы D к матрице D, мы сделали n перестановок столбцов. Применяя предложения 3 и 10 из лекции 5, имеем
j |
D |
j = ( 1)n jD j: |
(4) |
Б.М.Верников |
Лекция 10: Умножение матриц |
Определитель произведения матриц (7)
Матрица D является полураспавшейся матрицей с диагональными блоками C и E. Предложение 1 данной лекции и предложение 11 из
лекции 5 показывают, что jDj = jCj ( 1)n. Умножая обе части этого равенства на ( 1)n, имеем ( 1)n jDj = ( 1)2n jCj = jCj, т. е.
jCj = ( 1)n j |
D |
j: |
(5) |
Из равенств (2)–(5) вытекает, что
j Cj = ( 1)n jDj = ( 1)2n jD j = jD j = jDj = jAj jBj:
Теорема доказана.
Б.М.Верников |
Лекция 10: Умножение матриц |
Значение многочлена от квадратной матрицы (1)
Операции умножения и сложения матриц и умножения матрицы на число позволяют вычислять значения многочленов от квадратных матриц.
Очевидно, что если A квадратная матрица порядка n, то матрицы A2 = A A, A3 = A2 A, . . . , Ak = Ak 1 A, . . . существуют и являются
квадратными матрицами порядка n. В силу ассоциативности умножения матриц возведение матриц в степень обладает всеми теми же свойствами, что и возведение чисел в степень. А именно,
если A квадратная матрица, а m и n натуральные числа, то
Am An = Am+n, а (Am)n = Amn.
Определение
Если f (x) = anxn + an 1xn 1 + + a1x + a0 многочлен, а A
квадратная матрица, то значением многочлена f (x) от матрицы A
называется матрица
f (A) = anAn + an 1An 1 + + a1A + a0E;
где E единичная матрица того же порядка, что и A.
Б.М.Верников |
Лекция 10: Умножение матриц |
Значение многочлена от квадратной матрицы (2)
Можно рассматривать и матричные многочлены вида
F(x) = A0xr + A1xr 1 + + Ar 1x + Ar ;
где A0; A1; : : : ; Ar фиксированные квадратные матрицы порядка n. Значение такого многочлена от квадратной матрицы B порядка n вычисляется по формуле
F(B) = A0Br + A1Br 1 + + Ar 1B + Ar :
Б.М.Верников |
Лекция 10: Умножение матриц |