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

RUDN-I

.pdf
Скачиваний:
11
Добавлен:
25.03.2016
Размер:
721.58 Кб
Скачать

71

Тема 5. Алгебра матриц

Выше были определены операции сложения, вычитания матриц одного порядка, а также умножение матрицы на число. Здесь будет рассмотрена

еще одна важная операция — умножение матриц. Она играет ключевую

matematem

.

ru

 

роль в матричном методе исследования систем линейных уравнений и в матричном представлении действия линейных операторов (ему посвящена тема 16). Введено также понятие обратной матрицы, получен критерий ее существования и описаны способы вычисления. Изложен метод решения матричных уравнений и квадратных систем линейных уравнений с помощью обратной матрицы. В качестве следствия мы выводим формулы Крамера для решения квадратной системы линейных уравнений с невырожденной матрицей. В конце приведено дополнение (п. 5.8), которое содержит характеризацию определителя как полилинейной кососимметрической функции векторных аргументов.

5.1. Определение умножения матриц. Примеры

Пусть даны матрица-строка A = (a1 a2 . . . ap) и матрица-столбец той

 

b1

 

 

же длины B =

b2

 

...

.

 

 

bp

 

 

 

 

 

 

 

 

 

 

Определение 1. Произведением AB строки на столбец называется

число

.

p

 

 

 

 

 

 

 

AB = a1b1 + a2b2 + . . . + apbp = akbk

(1)

 

 

=1

 

 

 

k

 

Пусть теперь даны две матрицы A = aik размера (m × p) и B =

 

 

 

 

 

 

 

 

 

 

A1

 

 

 

 

 

bkj размера (p × n). Используем обозначение A

 

 

A2

 

 

 

 

 

=

 

 

...

 

 

, B

=

 

www

 

 

 

 

 

 

 

 

Am

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B1

B2 . . . Bn

,

где

A

i — строки матрицы

A, A

 

= (a

 

a

 

. . . a

 

),

[

 

]

 

 

i

 

 

i1

 

i2

 

ip

 

2 1
, B =
0 3
1 0
.
1 2
)
)
(n × n) · (n ×n) = (n × n),
| {z }

72

 

Тема 5. Алгебра матриц

 

 

 

 

b1j

.

Bj — столбцы матрицы B, Bj =

b2...j

 

bpj

 

 

 

 

 

 

 

Определение 2. Произведением C = AB матрицы A на матрицу B называется матрица C = cij размера (m × n), элементы которой

Умножение матриц подчинено правилу размеров, которое в определении

вычисляются по формуле

 

 

 

 

.

ru

 

j

 

 

 

 

(2)

cij = AiB , 1 ≤ i ≤ m,

1

 

 

 

≤ j ≤ n.

 

Итак, согласно определению (1)

 

 

 

 

 

 

 

 

b1j

 

 

 

 

 

 

 

b2j

(1)

p

 

 

 

cij = (ai1, ai2, . . . aip)

...

 

= k=1 aikbkj.

(3)

 

bpj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

matematem

 

 

 

(2) можно записать символически:

(m × p) · (p ×n) = (m × n)

(4)

|

 

{z

 

}

 

 

 

 

(длина строк 1-ой матрицы p должна равняться длине столбцов 2-ой

матрицы; число строк произведения равно числу строк 1-ой матрицы, а число столбцов произведения равно числу столбцов 2-ой матрицы).

Замечание 1. При умножении квадратных матриц одного порядка сно-

ва получим квадратную матрицу того же порядка.

www

A =

(

(

Действительно, пусть.

A

и

B — квадратные матрицы порядка n. Тогда

правило размеров (4) примет вид

т.е. C = AB — квадратная матрица порядка n.

Пример 1. Пусть

5.1. Определение умножения матриц. Примеры

73

Тогда

 

 

 

 

 

 

 

AB = (

2

1

1

0

(

3

2

0

3 )(

1

2 ) =

3

6 ),

BA = (

1

0

2

1

(

2

1

1

2 )(

0

3 ) =

2

7 ).

Вывод: вообще говоря, AB ̸= BA даже для квадратных матриц.

Если AB = BA, то матрицы A и B называют перестановочными, если

AB = −BA антиперестановочными.

 

.

ru

Пример 2. Пусть

 

 

 

 

 

 

 

 

 

a11

a12 . . . a1n

x1

 

a21

a22 . . . a2n

x2

 

A = ...

... ... ...

, X =

...

.

 

am1

am2 . . . amn

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда, правило размеров дает

 

 

 

 

 

 

(m × n) · (n ×1) = (m × 1),

 

 

т.е. AX — столбец

|

 

{z

 

}

 

 

 

 

 

 

 

 

a11

a12 . . . a1n

 

x1

 

 

 

 

a21

a22 . . . a2n

x2

 

 

 

... ... ... ...

...

=

 

 

 

am1 am2 . . . amn

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11x1

+ a12x2

+ . . . + a1nxn

 

 

a21x1

+ a22x2

+ . . . + a2nxn

.

matematem= ... ...

...

...

.

am1x1

+ am2x2

+ . . . + amnxn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечаниеwww2. Этот пример показывает, что систему m линейных уравнений с n неизвестными (1) п. 1.3 можно записать в форме одного матричного уравнения

AX = B.

Удобство такой записи отражает связь введенной операции умножения матриц с теорией линейных систем. Глубокая связь умножения матриц с теорией линейных отображений векторных пространств будет раскрыта в теме 15 (умножение матриц согласуется с умножением соответствующих линейных отображений).

74

 

 

 

 

 

 

 

 

Тема 5. Алгебра матриц

Пример 3. Пусть

 

 

 

 

 

 

 

 

 

 

 

A = (

0

1

 

 

2

 

0

 

 

Тогда

0 0 ) , B = (

0 0 ),

 

AB = (

 

 

 

 

 

 

 

 

 

 

 

0

1

2

0

 

 

0

0

 

cij = AiIj

0 0 )(

0 0 ) = (

0 0 )

 

= (ai1 . . . amatematemij . . . ain) 1 ← j-ый элемент

=

— пример делителей нуля:

A ̸= O, B

̸= O,

но

AB = O. ruОтметим, что в

этом примере BA ̸= O:

 

 

 

 

 

 

 

 

.

 

 

BA = (

2

 

0

0

1

 

 

0

2

 

 

0 0 )(

0 0 ) = (

 

0 0 ).

 

5.2. Свойства умножения матриц

1. Роль единичной матрицы.

Пусть I(m), I(n) — единичные матрицы порядков m и n. Тогда для матрицы A = aij размера (m × n) справедливы равенства

AI(n) = I(m)A = A.

(1)

Доказательство. Докажем, что AI(n) = A (второе — самостоятельно, см. упражнение 5.1). Обозначим AI(n) = C = cij . Правило размеров дает

 

 

 

(m × n) · (n ×n) = (m × n),

 

 

 

 

 

 

 

 

 

 

 

т.е.

C

имеет размер

(m

×

n)

 

Далее, по определению

 

 

 

 

.| {z }

 

 

 

.

 

 

 

 

 

 

 

0...

 

 

 

 

 

 

 

 

 

 

..

 

Пусть Awww= aij =

...

 

 

 

 

 

 

.

 

 

— матрица размера (m × p) со строками

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ai1 · 0 + . . . + aij · 1 + . . . + ain · 0 = aij.

Итак, cij = aij, 1 ≤ i ≤ m, 1 ≤ j ≤ n, т.е. I(n)A = C = A.

2. Строки произведения матриц.

 

 

 

 

 

A2

 

 

 

 

 

 

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Am

5.2. Свойства умножения матриц

 

 

 

 

 

 

 

75

 

 

 

 

 

 

 

B2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B1

 

 

 

 

 

 

 

 

Ai, i = 1, . . . , m, B

= bij

 

Bp

— матрица размера (p × n) со

 

 

 

= ...

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

 

ru

 

 

 

 

 

 

 

 

 

 

 

Cm

 

 

 

 

 

 

 

 

 

 

 

 

 

C2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

строками Bk, k = 1, . . . , p и пусть C = AB = ...

. Тогда справедливы

равенства

 

 

 

matematem

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

Ci

= AiB, i = 1, . . . , m,

 

 

(2)

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

(3)

 

 

 

 

Ci =

aikBk,

i = 1, . . . , m.

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

Доказательство. Докажем (2). Для Di

= AiB имеем по правилу разме-

ров:

 

 

 

 

(1 × p) · (p ×n) = (1 × n),

 

 

 

т.е. Di — строка (1 × n), Di |={z(d}i1, . . . , dij, . . . , din),

i = 1, . . . , m причем

по определению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

di1 = AiB1 = ci1, . . . , dij = AiBj = cij, . . . , din = AiBn = cin,

т.е. Ci = Di = AiB.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Докажем (3). Обозначим теперь Ei =

 

kp=1 aikBk — это линейная комби-

нация строк

B1, . . . , Bp,

т.е. тоже

строка (1

×

n). Далее,

 

 

 

.

 

 

 

 

 

 

 

 

Ei = ai1B1

+ ai2B2

+ . . . + aipBp =

 

 

 

 

 

 

 

 

www

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ai1(b11 b12 . . . b1n) + ai2

(b21 b22 . . . b2n) + . . . + aip(bp1 bp2 . . . bpn) =

= (ai1b11 + ai2b21 + . . . + aipbp1 ai1b12 + ai2b22 + . . . + aipbp2 . . .

. . . ai1b1n + ai2b2n + . . . + aipbpn) = (AiB1 AiB2

 

(2)

. . . AiBn) =

 

 

 

 

 

 

 

 

 

 

 

 

 

= (ci1 ci2 . . . cin) = Ci.

Итак, Ci = Ei =

 

kp=1 aikBk, i = 1, . . . , m т.е. верно (3).

 

Из формулы (2)

получим следующие следствия.

 

 

 

Следствие 1. Если i-ая строка матрицы A умножена на некоторое число, то i-ая строка матрицы C = AB умножается на то же число.

(B + C)A = BA + BC.
A(B + C) = AB + AC.

76

Тема 5. Алгебра матриц

Следствие 2. Если в матрице A поменять i-ую и j-ую строки, то в матрице C = AB также меняются местами i-ая и j-ая строки.

Следствие 3. Если к i-ой строке матрицы A прибавить j-ую с некоторым коэффициентом, то это же происходит со строками матрицы

C = AB.

Вывод. Элементарные преобразования строк матрицы A приводят к

тем же элементарным преобразованиям строк матрицы C = AB.

Иными словами, справедливо следующее утверждение..

ru

 

Лемма 1. Пусть матрица A получена из матрицы A элементарными

преобразованиями строк. Тогда матрица C = AB получается из C =

AB теми же элементарнымиe преобразованиями строк.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

e

 

 

 

 

 

 

3. Столбцы произведения матриц.

 

 

 

 

 

 

 

 

 

 

Пусть A =

 

aij

 

=

A1 A2 . . . Ap

 

— матрица размера (m

×

p) со

 

k

 

 

 

 

 

 

 

bij

= B

1

B

2

. . . B

n

 

 

столбцами

A

 

, k = 1, . . . , p, B =

 

 

 

 

— матри-

 

p

 

 

n

 

[

 

j ]

 

 

 

 

 

 

ца размера (

×

)

со столбцами

B , j = 1, . . . , n

и пусть

C = AB =

 

 

 

 

 

 

[

 

 

 

]

 

 

[ C1 C2 . . .

 

Cn ]. Тогдаj

справедливыj

равенства

 

 

 

 

(4)

 

 

 

 

 

 

 

 

C

= AB ,

j = 1, . . . , n,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cj =

bkjAk,

j = 1, . . . , n.

 

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

Упражнение 1. Вывести формулы (4), (5) (см. упражнение 5.7).

Из формулы (4) вытекает вывод: элементарные преобразования столб-

1) Пусть B и C — матрицы одного размера, причем умножение матриц AB, AC допускается правилом размеров. Тогда

цов матрицы B приводят к тем же элементарным преобразованиям

столбцов матрицы C = AB.

 

 

matematem

4. Дистрибутивность умножения относительно сложения.

www

.

 

(6)

2) Пусть B и C — матрицы одного размера, причем умножение BA, BC допускается правилом размеров. Тогда

(7)

Упражнение 2. Вывести формулы (6), (7) (см. упражнение 5.2).

5.2. Свойства умножения матриц

77

5. Ассоциативность умножения.

Пусть A, B, C — матрицы размеров (m × p), (p × k), (k × n) соответственно. Тогда

A(BC) = (AB)C. (8)

Замечание 1. Обе части (8) имеют смысл по правилам размеров: BC

Для доказательства (8) потребуется лемма о перестановке порядка суммирования в двойных суммах.

имеет размер (p×k)(k×n) = (p×n), A(BC) имеет размер (m×p)(p×n) =

(m × n). Таков же размер (AB)C.

.

ru

matematem

 

 

 

Лемма 2. В двойных суммах можно менять порядок суммирования:

для любых чисел αij,

i = 1, . . . , m,

j = 1, . . . , n

 

 

 

 

m

n

 

n

m

 

 

 

 

=1 (j=1 αij)

= j=1 (i=1 αij).

(9)

 

i

 

∑ ∑

 

 

Доказательство. Рассмотрим матрицу

 

 

 

 

 

 

α11

α12 . . .

α1n

 

 

 

 

 

α21

α22 . . . α2n

 

 

 

Λ =

...

 

... ...

...

.

 

 

 

 

αm1

αm2 . . . αmn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сумму всех ее элементов S(Λ) можем найти двумя способами.

 

1 способ:

 

 

 

 

 

 

 

 

a) сложим все элементы i-ой строки, i = 1, . . . , m,

 

b) сложим полученные суммы по всем строкам.

 

 

Итак,

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

Si) =

αij

— сумма элементов i-ой строки, i = 1, . . . , m,

 

 

=1

 

 

 

 

 

 

 

 

j

 

m

 

m

n

αij).

 

 

S(Λ) =

=1 (Si)) = i=1

(j=1

(10)

 

 

 

i

 

∑ ∑

 

2 способ:www

a) сложим все элементы j-го столбца, j = 1, . . . , n, b) сложим полученные суммы по всем столбцам.

78

 

 

 

 

 

 

 

 

 

 

Тема 5.

Алгебра матриц

Итак,

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ru

 

 

 

m

 

 

 

 

 

 

 

 

 

Sj) =

 

 

 

 

 

 

 

 

 

 

 

 

 

αij — сумма элементов j-го столбца, j = 1, . . . , n,

 

 

=1

 

j

(

 

)

∑ ∑

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

m

 

 

 

 

 

 

S(Λ) =

=1

Sj)

= j=1 (i=1 αij).

 

(11)

 

 

 

 

 

 

 

 

 

 

 

.

 

 

Это и означает, что E matematem= G требовалось.

 

 

Сопоставляя (10) и (11), получаем (9).

 

 

 

 

Доказательство свойства 5. Введем обозначения

 

 

 

 

 

 

D = AB,

E = (AB)C = DC,

 

(12)

 

 

 

 

F = BC,

G = A(BC) = AF.

 

(13)

Нужно показать, что E = G. Из (12) по правилу умножения матриц

eij = ∑k

dikckj (т.к. E = DC),

dik = ∑l

ailblk (т.к. D = AB).

Подставляя, получим

 

 

 

 

 

 

 

 

 

eij = ∑k

(l

ailblj)ckj = k

l

 

 

(9)

 

 

 

ailblkckj =

 

ail (∑k

blkckj).

 

 

 

 

.

 

 

= l

k

ailblkckj = ∑l

 

 

 

 

 

 

 

 

 

 

l ailflj = gij т.к. G = AF .

Но k blkckj = flj т.к. F = BC. Итак, eij =

www

 

, что и

 

 

 

 

 

 

 

 

 

 

 

 

6. Транспонирование матриц.

Пусть A, B — матрицы размера (m × p) и (p × n), соответственно. Тогда

 

 

 

(AB)T = BT AT .

 

 

 

 

 

 

(14)

Замечание 2.

 

 

AB

 

 

(m

×

n)

 

(AB)

T

 

По правилу размеров T

 

имеет размер

×

 

 

T, тогда

 

имеет размер (n × T ).TДалее,

B

имеет размер (

n

p

),

A

 

имеет

 

 

 

m

 

 

 

 

 

 

 

размер (p × m), так что B A

имеет смысл, и ее размер (m × n).

 

Упражнение 3. Доказать формулу (14) (см. упражнение 5.2).

5.3. Свойства умножения квадратных матриц

79

5.3.Свойства умножения квадратных матриц

Обозначим Mn — множество всех квадратных матриц порядка n. Для A, B Mn по правилу размеров C = AB Mn:

(n × n) · (n ×n) = (n × n).

| {z }

Поэтому, кроме общих свойств умножения матриц, умножение квадрат-

 

 

 

matematem

ru

 

 

ных матриц имеет ряд специфических свойств.

 

 

 

1. Роль "единицы" при умножении квадратных матриц.

из Mn играет

единичная матрица порядка n I = I(n) Mn:

 

 

 

 

 

 

 

 

AI = IA = A, A Mn.

 

 

 

 

 

2. Для A Mn, m = 0, 1, 2, . . . имеют смысл степени

 

 

(1)

A0 = I,

A1 = A,

A2 = AA, . . . , Am = Am−1A,

m = 2, 3, . . .

Очевидно, что

 

 

Am = A · A · . . . · A, m ≥ 2.

 

 

 

 

 

Степени

 

|

 

 

{z

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

m раз

 

 

 

 

 

 

 

Замечание 1.

 

 

 

матрицы коммутируют между собой:

 

 

AmAp = ApAm = Am+p,

m, p N0 = {0, 1, 2, . . .}

 

(2)

3. Имеют смысл многочлены от матриц A Mn.

 

 

 

 

Пусть

 

 

p(λ) = p0 + p1λ + . . . + pkλk

 

 

 

 

 

 

 

 

 

 

 

— многочлен степени k,

p0, p1, . . . , pk R. Полагаем

 

 

 

(3)

www

p(A) = p0I + p1A + . . . + pkAk Mn.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эту матрицу называют многочленом от матрицы A.

 

 

 

 

Пример.

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(λ) = λ2 + λ − 2, A = (

1

1

 

 

 

 

 

 

0 2 ) ,

 

 

 

 

p(A) = A2 + A − 2I = (

1

1

 

 

1

1

+ (

1

1

2 (

1

0

0 2 )(

0 2 )

0 2 )

0 1 ) =

( 0

4 )

+ (

0 2 ) + (

0 2 ) = (

0 0 ) .

 

 

1

1

 

 

 

1

 

1

 

 

2

0

 

0

0

 

 

Этот результат не случаен (см. упражнение 5.13)!

4. Определитель произведения квадратных матриц.

= det A = sign σ · det A.

80

Тема 5. Алгебра матриц

Теорема 1. При умножении квадратных матриц их определители перемножаются, т.е. если A, B Mn, то

det AB = det A · det B.

ru

(4)

 

Для доказательства теоремы 1 потребуются некоторые вспомогательные

утверждения.

.

 

Обозначим C = AB Mn. Согласно формуле (5) для столбцов матрицы C справедливы равенстваmatematem

n

Cj = bkjAk, j = 1, . . . , n,

k=1

т.е. Cj есть линейная комбинация столбцов A1, . . . , An матрицы A, причем коэффициенты ее — элементы j-го столбца матрицы B. Эти формулы при разных значениях j будут использованы совместно, поэтому индексы суммирования в них следует обозначать разными буквами, например

 

 

 

 

kj

 

 

 

 

 

 

n

 

 

 

 

 

Cj =

bkjjAkj , j = 1, . . . , n

 

(5)

 

 

 

 

=1

 

 

(причем k1, . . . , kn независимо друг от друга меняются от 1 до n).

 

Лемма 1. Пусть A =

A1 A2 . . . An Mn, σ

= (k1, k2, . . . , kn)

 

Sn

— перестановка

чисел (1, 2, . . . , n), A(σ) = Ak1

Ak2 . . . Akn

 

 

[

]

 

Mn — матрица с переставленными столбцами.[Тогда

 

]

 

 

 

det A(σ) = sign σ · det A.

 

(6)

Доказательство леммы.

из A транспозицией (заменой местами) двух столбцов. По свойству 2

1) Если σ есть транспозиция, то sign σ = 1. При этом A(σ) получается

www

случае

определителей в этом.

det A(σ)

2) В общем случае перестановка σ есть произведение нескольких транспозиций, причем их число четно, если σ — четная, и нечетно, если σ — нечетная (см. теорему 2 п. 3.2). При этом A(σ) получается из A соответствующими транспозициями столбцов. При каждой из них знак определителя меняется. В результате, для четной перестановки σ он сменится четное число раз, т.е. останется прежним:

det A(σ) = det A = sign σ · det A, если σ — четная.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]