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

9698

.pdf
Скачиваний:
0
Добавлен:
25.11.2023
Размер:
3.08 Mб
Скачать

1.3. Умножение матриц. Произведением матрицы A из m строк и k столбцов на матрицу B из k строк и n столбцов называется матрица

C = A × B

a11

a21

Kai1

Kam1

a12

a22

K

ai 2

K am2

K a2k

 

b11 K b1 j

K b1n

c

K a2k

 

11

 

 

 

 

 

L

K K

 

× b21

K b2 j

K b2n

= ci1

K aik

K K K K K

L

K K

 

K bkj

 

 

bk1

K bkn

 

 

 

 

 

 

cm1

K amk

 

 

 

 

L c1 j

L c1n

 

 

L L L L

 

L cij

L cin

 

,

 

L L L L

 

L cmj

 

 

 

L cmn

 

имеющая m строк и n столбцов, элемент которой cij , стоящий на пересечении строки i со столбцом j , равен сумме произведений элементов строки i матрицы A на соответствующие элементы столбца j матрицы B

cij = ai1b1 j + ai 2b2 j + K + aikbkj .

Обратим внимание на то, что число столбцов первого сомножителя

должно равняться числу строк второго сомножителя.

Пример.

1

2 5

6

1× 5 + 2 × 7

1× 6 + 2 ×8

19

22

,

 

 

×

 

=

×5

+ 4 × 7

 

=

 

 

3

4

7

8

3

3 × 6 + 4 ×8

43

50

 

5

6 1

2

5 ×1 + 6 ×3

5 × 2 + 6 × 4

23 34

 

 

×

 

=

 

+ 8 ×3

 

 

 

=

.

7

8 3

4

7 ×1

7 × 2 + 8 × 4

31 46

 

Отсюда видно, что A × B ¹ B × A . И, вообще, из существования произведения

A × B совсем не следует существование произведения B × A .

 

Единичная матрица порядка n

– это квадратная матрица

 

 

 

 

 

1

0 ...

0

0

 

 

 

 

 

 

 

0

1 ...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = ... ... ...

...

... ,

 

 

 

 

 

 

0

0 ...

1

0

 

 

 

 

 

 

 

0 ...

0

 

 

 

 

 

 

 

0

1

 

 

где на главной диагонали стоят единицы, а все остальные элементы равны нулю, играет роль «единицы» в алгебре матриц, ибо A × E = E × A , например:

10

a

a

 

1

0

a

a

 

11

12

 

×

 

= 11

12

.

a21

a22

0

1

a21

a22

1.4. Матрицы и линейные преобразования. Для «оправдания» определения операции умножения матриц рассмотрим связь между матрицами и линейными преобразованиями.

Пусть на плоскости задана декартова прямоугольная система координат xOy . Поставим в соответствие произвольной точке плоскости с

координатами (x, y) точку ( x, y) плоскости xOy′, координаты которой определяются следующими линейными соотношениями

x= a11 x + a12 yy¢ = a21x + a22 y ,

задающими так называемое линейное преобразование плоскости.

Очевидно, что такое преобразование однозначно определяется его матрицей

a

a

 

A = 11

12

.

 

a22

 

a21

 

Рассмотрим несколько примеров таких преобразований. Пусть преобразование задано следующим образом:

x= x

или, подробнее,

x= 1× x + 0 × y

 

 

.

y¢ = by

 

 

 

y¢ = 0

× x + b × y

 

1

0

 

 

 

Соответствующая ему матрица A =

β

. Это преобразование – « сжатие»

 

0

 

 

 

вдоль оси Oy (или к оси Ox ) с коэффициентом β > 1 (см. рис. 1.1).

y y

 

β

1

1

 

x

 

x

O

1

O

1

Рис. 1.1

11

 

x′ = αx

 

 

α

0

Преобразование

 

с матрицей A =

является «сжатием»

 

y′ = y

 

 

0

1

вдоль оси Ox (или к оси Oy ).

 

 

 

 

α

0

 

соответствует «деформации» плоскости

Матрица A =

β

 

0

 

 

 

 

 

 

одновременно вдоль оси

 

Ox и Oy .

 

 

 

Единичная матрица

 

1

0

 

задает тождественное преобразование

 

E =

 

 

 

 

 

0

1

 

 

 

{ x′ = x

плоскости y′ = y .

В качестве упражнения напишите формулы преобразования, которое произвольную точку «переводит» в точку, симметричную ей относительно начала координат, и получите соответствующую матрицу.

Аналогичным образом можно интерпретировать линейные преобразования для трех переменных

x′ = a11x + a12 y + a13 z

′ = + +

y a21x a22 y a23 z .z′ = a31x + a32 y + a33 z

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

x′ = b11x + b12 y

y′ = b21 x + b22 y

с матрицей B =

b

b

 

 

переводит точку (x, y) в точку ( x, y) , а линейное

11

12

 

 

 

b21

b22

 

 

 

 

преобразование

 

 

 

 

x′′ = a x′ + a y

 

 

 

 

 

 

 

 

 

 

 

11

12

 

 

 

 

 

 

 

 

 

 

 

 

y′′ = a21x′ + a22 y

с матрицей A =

a

a

 

 

переводит точку ( x, y) в точку ( x′′, y′′) . Эти два

11

12

 

 

a21

a22

 

 

 

преобразования, проведенные друг за другом, переводят точку (x, y) в точку

( x′′, y′′) . Очевидно, что результат этих преобразований тоже линейное преобразование. Какова его матрица?

12

Выразим координаты точки ( x¢¢, y¢¢) через координаты исходной точки

(x, y)

x¢¢ = a

(b x + b y) + a

 

(b x + b y ) = (a b + a b

) x + (a b + a b

) y

 

11

11

12

12

 

21

22

11

11

12

21

 

11

12

12

22

y¢¢ = a

(b x + b y ) + a (b x + b y) = (a b + a b ) x + (a b + a b ) y

 

21

11

12

 

22

21

22

21

11

22

21

21

12

22

 

22

 

 

 

Таким образом, матрица результирующего преобразования равна произведению матриц A × B исходных преобразований, взятых в соответствующем порядке, т.е.

a

a

 

b

b

 

a b

+ a b

a b

+ a b

 

11

12

 

× 11

12

 

= 11 11

12 21

11 12

12 22

.

a21

a22

b21

b22

a21b11

+ a22b21

a21b12

+ a22b22

 

Введенные выше матричные операции обладают следующими

свойствами:

( A× B) ×C = A×( B ×C)

a( A× B) = (aA) × B = A×(a × B) ( A + B) ×C = A ×C + B ×C .

Определим операцию транспонирования матрицы как замену в ней строк на соответствующие столбцы и наоборот

a

a K a

 

 

a

a K a

 

11

 

12

1n

 

 

11

21

m1

 

a21

a22 K a2n

 

T

a12

a22 K am 2

A = K K

K

,

A

= K K K

.

a

a

m2

K a

 

 

a

a

K a

 

m1

 

mn

 

 

1n

2n

mn

 

 

 

 

 

 

 

 

 

 

 

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

 

 

2

 

2

 

2

 

 

x1

 

T

x

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

x1 + x2 + x3 = ( x1 x2 x3 ) x2

= X X ,

 

 

 

 

 

 

 

 

 

 

 

 

X = x2

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax

2

+ 2bxy + cy

2

=

( x

a b

x

 

= X

T

 

a b

 

x

 

 

y )

 

 

 

AX , A =

 

 

X =

.

 

 

 

 

 

 

b c

y

 

 

 

 

b c

 

y

13

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

Лекция 2. Правило Крамера и определители матриц

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

Бертран Рассел (1872-1970гг.)

2.1. Системы двух уравнений с двумя неизвестными. Уже при решении уравнения первой степени

a × x = b

возможны три случая:

если a ¹ 0 , то уравнение имеет единственное решение x = b a ;

если a = 0 и b ¹ 0 , то уравнение не имеет решения;

если a = 0 и b = 0 , то уравнение имеет бесчисленное множество решений, т.к. равенство 0 × x = 0 выполняется при любых значениях x .

Оказывается, что такие же случаи имеют место при решении системы

уравнений с большим числом неизвестных. В стандартной записи система двух линейных уравнений с двумя неизвестными имеет вид

 

 

a1 x + b1 y = c1

,

(2.1)

 

 

a2 x + b2 y = c2

 

 

 

 

где коэффициенты a1 , a2

и b1 , b2 при неизвестных x и y , а также правые

части уравнений c1 и c2

заданные действительные числа.

x = α , y = β,

Решением системы (2.1)

называется такая пара чисел

которая оба уравнения системы (2.1) обращает в тождества. В этом случае говорят, что пара чисел (α,β) «удовлетворяет» системе уравнений (2.1), а система уравнений называется совместной. Если система уравнений не имеет решений, то она называется несовместной. Эта терминология относится к системам произвольного числа уравнений.

Существует несколько способов нахождения решения системы, если оно существует. Проделаем следующие операции:

 

a x + b y = c

 

× b

 

-

1

1

1

 

2 .

a2 x + b2 y = c2

 

× b1

 

 

14

 

 

 

В результате получим

x(a1b2 a2b1 ) = c1b2 c2b1 .

Если число a1b2 - a2b1 ¹ 0 , то

x =

 

c1b2

c2b1

 

 

(2.2)

a1b2

a2b1

 

 

Аналогичным образом, исключая из системы неизвестное x

(проверьте!),

получим

 

a2c1

 

 

y =

a1c2

.

(2.3)

 

 

 

 

a1b2

a2b1

 

Формулы (2.2) и (2.3) достаточно громоздки, и от их вида «в глазах рябит». Вот здесь и приходят на помощь те «удачные обозначения», о которых говорил Б. Рассел. Выпишем коэффициенты при неизвестных в виде следующей матрицы:

a

b

 

A = 1

1

.

a2

b2

 

Число

( A) = a1b2 a2b1

называется определителем матрицы A и

обозначается следующим образом:

 

 

 

 

 

 

a1

b1

 

 

= ( A) = a b a b .

 

 

 

 

 

a2

b2

 

 

1

2

2

1

 

 

 

 

 

 

 

 

Для удобства в обозначение определителя включена матрица, которой он соответствует.

Вычислим ещё два определителя матриц, полученных из матрицы A путем замены ее первого и соответственно второго столбца на столбец правой части системы, т.е.

x

=

c1

b1

= c b c b

,

y

=

a1

c1

= a c

c a

2

.

 

c2

b2

1 2 2 1

 

 

a2

c2

1 2

1

 

 

 

 

 

 

 

 

 

 

 

Тогда формулы (2.2) и (2.3) в случае, когда D ¹ 0 , можно представить в виде

15

 

 

c1

b1

 

 

D

 

 

 

 

 

a1

c1

 

 

 

Dy

 

x =

 

c

 

b

 

=

 

 

y =

 

 

a

 

c

 

 

=

 

 

 

2

2

 

 

x

,

 

 

 

2

2

 

 

 

.

 

 

 

 

 

 

 

 

 

a1

b1

 

 

 

 

 

a1

b1

 

 

D

 

 

 

 

D

 

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

Эти формулы связаны с именем швейцарского математика Г. Крамера (1704-1752гг.). Сама идея определителя появилась ещё у Г. Лейбница

(1646-1716гг.), а название дал в 1812г. О. Коши (1789-1857гг.). Это название довольно естественно: определители , D x и y полностью

определяются системой уравнений и, наоборот, по заданным определителям (в принятых обозначениях) можно «восстановить» систему уравнений.

Исследуем полученные формулы. Будем предполагать, что не все коэффициенты при неизвестных одновременно равны нулю. В противном случае система (2.1) примет вид

0 × x + 0

× y = c

(2.4)

 

1 .

0 × x + 0

× y = c2

 

Тогда, если c1 = c2 = 0 , то решением является любая пара чисел;

если же

хотя бы одно из чисел c1 или c2 отлично от нуля, то система (2.4) не имеет решений и называется несовместной.

Учитывая принятые обозначения для определителей, имеем два

уравнения

 

 

 

D × x = Dx

,

 

D × y = Dy

 

 

при условии, что не все входящие в эти уравнения определители равны нулю, или одно из этих уравнений, если все эти определители равны нулю.

Отсюда получаем следующие утверждения:

· если определитель ¹ 0, то система имеет единственное решение,

определяемое по формулам Крамера

x =

x

,

y =

Dy

;

D

D

 

 

 

 

· если определитель

= 0 и хотя бы один из определителей D x или

yотличен от нуля, то система несовместна (пусть, для определённости,

=0 и D x ¹ 0 , тогда получим противоречивое равенство 0 × x = Dx ¹ 0 );

16

· если все три определителя системы равны нулю, то система имеет бесчисленное множество решений: x – любое, а y выражается через x из

любого уравнения исходной системы или наоборот: y

любое, а x

выражается через y из любого уравнения системы.

 

Отметим, что определитель обращается в нуль, когда

a1b2 a2b1 = 0 ,

т.е. элементы какой-то строки (столбца) равны нулю или элементы строк (столбцов) пропорциональны

a1 = b1 . a2 b2

С учетом этого переформулируем приведенные выше утверждения в терминах коэффициентов системы:

· если коэффициенты при неизвестных не пропорциональны, т.е.

a1 ¹ b1 или a1b2 - a2b1 ¹ 0 , то система имеет единственное решение; a2 b2

· если коэффициенты при неизвестных пропорциональны, но не пропорциональны свободным членам,

a1

=

b1

¹

c1

(D = 0, Dx ¹ 0) , или Dy ¹ 0 или оба не равны нулю,

a2

b2

c2

 

 

 

то система несовместна; · если коэффициенты при неизвестных пропорциональны и

пропорциональны свободным членам,

a1

=

b1

=

c1

( = 0,

x = 0,

y = 0 ) ,

a

 

b

 

2

 

 

c

2

 

 

 

 

2

 

 

 

 

 

то система имеет бесчисленное множество решений.

Эти утверждения становятся очевидными, если вспомнить, что каждое уравнение первой степени задаёт прямую в плоскости xOy , и решить систему (2.1), значит найти общие точки двух прямых. В первом случае прямые пересекаются в одной точке, координаты которой и будут решением системы, во втором – прямые параллельны (решения нет) и, наконец, в третьем случае прямые совпадают, и координаты любой точки на этой общей прямой будут решением системы.

2.2. Системы трех уравнений с тремя неизвестными. Выше мы показали, как вычисляется определитель матрицы второго порядка

D( A) =

a1

b1

= a1b2 a2b1 .

 

a2

b2

 

Для вычисления определителя матрицы третьего порядка

17

a11

a12

A = a

a

22

21

 

 

a32

a31

a13 a23 a33

нам потребуются понятия минора и алгебраического дополнения.

Пусть мы выбрали некоторый элемент aij (этот элемент стоит на пересечении строки i со столбцом j ; пара (i j) его «адрес» в матрице).

Вычеркнем теперь строку i и столбец j , где стоит этот элемент, и получим определитель второго порядка. Его называют дополнительным минором элемента aij (или просто минором) и обозначают M ij .

Алгебраическим дополнением элемента aij назовем величину

Aij = (−1)i+ j Mij .

Таким образом, алгебраическое дополнение – это минор с соответствующим знаком (+) или (–). Знаки алгебраических дополнений элементов определителя легко определяются просто «глядя» на этот определитель. При этом полезно иметь в виду следующую схему

+− +

− + − ,

+− +

где знаком (+) отмечены места тех элементов, алгебраические дополнения которых равны минорам, взятым с их собственными знаками.

Теперь мы можем дать правило вычисления определителя

матрицы третьего порядка

a11 a12 a13

a21 a22 a23 = ( A) = ai1 Ai1 + ai 2 Ai 2 + ai3 Ai3 , i = 1, 2,3 ,

a31 a32 a33

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

Оказывается, что разложение по столбцам дает тот же результат:

a11 a12 a13

a21 a22 a23 = ( A) = a1 j A1 j + a2 j A2 j + a3 j A3 j , j = 1, 2,3 . a31 a32 a33

18

Где гарантия, что в результате будет получаться одно и то же число? Оказывается, что справедливо следующее утверждение.

Для каждой квадратной матрицы порядка n сумма произведений

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

Вот это число и называется определителем, а способов его вычисления существует много и тот, который даётся введённым выше определением, просто один из них. Этот способ называется разложением по элементам строки или столбца.

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

( A) = a11a22a33 + a12a23a31 + a21a32a13 a31a22a13 a21a12 a33 a32a23a11 .

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

Решение системы трех уравнений с тремя неизвестными

a x + b y + c z = d

1 1 1 1

a2 x + b2 y + c2 z = d2a3 x + b3 y + c3 z = d3

выражается через определители третьего порядка по аналогичным формулам Крамера:

 

 

d1

b1

c1

 

 

 

 

 

 

a1

d1

c1

 

 

 

 

 

a1

b1

d1

 

 

 

 

 

 

 

d2

b2

c1

 

 

 

 

 

 

a2

d2

c2

 

 

 

 

 

a2

b2

d2

 

 

 

 

 

x =

 

d3

b3

c3

 

=

x

 

; y =

 

a3

d3

c3

 

=

y

; z =

 

 

a3

b3

d3

 

 

=

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

b1

c1

 

 

 

 

 

a1

b1

c1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

b1

c1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

c1

 

 

 

 

 

 

 

a2

b2

c1

 

 

 

 

 

 

a

b c

 

 

 

 

 

 

a3

b3

c3

 

 

 

 

 

 

 

a3

b3

c3

 

 

 

 

 

 

 

2

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a b c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

3

 

 

 

 

 

 

в предположении, что определитель матрицы системы

 

 

не равен нулю.

Отметим, что

x , y ,

z

определители матриц, получаемых из матрицы

системы путем замены столбца коэффициентов при соответствующей

неизвестной на столбец правых частей. Если

= 0 и существует хотя бы

один из определителей

x , y , z , отличный от

нуля,

то

система

несовместна. Рассмотрение оставшегося случая

=

x =

y =

z = 0 мы

отложим до введения понятия ранга матрицы.

19

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