Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

.pdf
Скачиваний:
26
Добавлен:
23.02.2015
Размер:
285.47 Кб
Скачать

Определитель полураспавшейся матрицы (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: Умножение матриц