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

ЛЕКЦИЯ 01_Л.А

..pdf
Скачиваний:
95
Добавлен:
02.06.2015
Размер:
873.43 Кб
Скачать

ЛЕКЦИИ КАФЕДРЫ МАТЕМАТИКИ НИУ ВШЭ НН

21 ►

Покажем, что преобразования (1.10) можно осуществить при помощи умножения матрицы на некоторую квадратную матрицу K . Для определенности займемся преобразованиями 1 – 3 строк заданной матрицы.

Пусть A M m , n , так что определено произведение Em A, равное, очевидно, самой матрице A . В соответствии с равенством (1.8)’’ запишем

m

 

m

(Em A)i (Em )is as is as ai , i 1, m .

s 1

 

s 1

 

 

is

 

Такая интерпретация равенства Em A A в форме равенства строк с одинаковы-

ми номерами в матрицах Em A и A наталкивает на мысль искать квадратную матрицу

K M m , модифицируя подходящим образом единичную матрицу Em и показывает, что

нужное элементарное преобразование строк в A реализуется ее умножением на K слева.

1. Пусть в матрице K (m m) i я ( j я) строка есть

j я ( i я) строка из Em , а осталь-

ные строки – те же, что в Em : kis js , k js is , krs rs

при r i , r j для всех s 1, m .

Тогда

 

 

 

 

m

m

 

 

 

(K A)i

 

kis as js as

 

 

 

 

;

 

 

 

 

 

 

 

a j

 

 

 

 

 

 

 

s 1

s 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

m

 

 

 

 

 

k js as is as

ai

;

 

 

 

 

 

 

(K A) j

 

 

 

 

 

 

 

s 1

s 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

m

 

 

 

(K A)r

krs as rs as

ar

при

 

.

 

 

 

r i, r j

 

 

 

 

 

 

 

s 1

s 1

 

 

 

 

 

 

 

 

 

 

Следовательно, для того, чтобы поменять в матрице

A(m n)

местами две стро-

ки, достаточно умножить ее слева на матрицу K , которая получается из Em

при

помощи транспозиции строк с теми же номерами.

 

 

 

2. Пусть теперь в матрице K (m m) i я строка есть умноженная на число 0

i я

строка в Em , а остальные строки – те же, что в Em : kis is ,

krs rs

при r i , s 1, m .

 

 

Отсюда вытекает, что

 

 

 

 

 

 

 

m

m

 

 

 

(K A)i

 

is as is as

ai

;

 

 

 

 

 

 

 

s 1

s 1

 

 

 

 

 

Л Е К Ц И Я 1

Н.Н.БОБКОВ

ЛИНЕЙНАЯ АЛГЕБРА ◄ 22 ►

m

 

(K A)r rs as ar при

 

.

 

r i

 

 

 

s 1

Таким образом, чтобы умножить в матрице A(m n) строку на число, отличное

от нуля, достаточно умножить эту матрицу слева на матрицу K , которая получа-

ется из Em умножением в ней на это число строки с тем же номером25.

3. Пусть, наконец, в матрице K (m m)

i я строка есть сумма i

й и

j й строк единич-

ной матрицы Em , а остальные строки – те же, что в Em : kis is

js ,

i j ; krs rs

при

r i , s 1, m .

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

m

m

m

 

 

 

(K A)i

 

( is js ) as is as

js as

 

, i j ;

 

 

ai a j

 

 

 

 

 

 

s 1

s 1

s 1

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

(K A)r

rs as

ar

при

 

.

 

 

 

 

 

 

r i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s 1

 

 

 

 

 

 

 

 

 

 

 

Итак, чтобы прибавить в матрице

A(m n) к ее i й строке ее

j ю строку,

дос-

таточно умножить A слева на матрицу K , которая получается из Em прибавлени-

ем к ее i й строке ее j й строки.

Заметим, что описанное преобразование матрицы Em равносильно замене ее эле-

мента ij 0 , i j , единицей, поскольку

(Em )i (Em ) j (0 0

1

0 0) (0 0

1

0 0) , так что

 

i е место

 

j е место

 

(K )i (0 0 1 0 1 0) .

(i,i)

(i, j)

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

Л Е К Ц И Я 1

ЛЕКЦИИ КАФЕДРЫ МАТЕМАТИКИ НИУ ВШЭ НН

23 ►

Полученные общие результаты иллюстрируются ниже на простых примерах.

Примеры:

1 Умножение 1-й строки в A на число , 0 .

0

0

 

a11

a12

 

 

a11

 

0

1

0

 

 

a

a

22

 

 

 

a

21

 

 

 

 

 

 

21

 

 

 

 

 

 

0

0

1

 

 

a

a

 

 

 

a

 

 

 

 

 

 

31

 

32

 

 

 

 

31

 

 

 

 

 

 

 

K

 

 

A( m n ), m 3, n 2

 

 

 

 

a12

a22 . a32

Матрица K получена умножением 1-й строки в E3 на .

2 Перестановка 2-й и 3-й строк в матрице A .

a11

a12

 

 

a11

a12

 

 

1

0

0

 

a11

a12

 

a

a

 

 

a

a

 

 

 

0

0

1

 

 

a

21

a

 

21

22

 

r2 r3

31

32

 

 

 

 

 

 

 

 

 

22

 

a

a

 

a

a

 

 

 

0

1

0

 

 

a

a

 

31

32

 

 

21

22

 

 

 

 

 

 

 

 

31

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

Матрица K получена транспозицией 2-й и 3-й строк в E3 .

3 Прибавление к 1-й строке матрицы A ее 3-й строки.

a11

a12

 

 

a11 a31

a12 a32

 

 

1

0

1

 

a11

a

a

 

 

 

a

a

 

 

 

0

1

0

 

 

a

21

22

 

r r r

 

21

22

 

 

 

0

0

1

 

 

21

a

a

 

1 1 3

 

a

a

 

 

 

 

 

a

31

32

 

 

 

31

32

 

 

 

 

 

 

 

 

31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к1 й строке матрицы A

 

 

 

 

K

 

 

 

 

 

 

 

 

 

прибавилась ее 3 я строка

 

 

 

 

 

 

 

 

 

a12 a22 . a32

Матрица K получена заменой в E3 нуля с индексами (1,3) 26 на единицу.

Целесообразно проверить (сделайте это), что в этом примере умножение A слева

1

0

0

 

 

на K

0

1

0

приведет к прибавлению к ее 3 й строке

2 й строки, а умножение

 

0

1

1

 

 

1

0

0

 

слева на K

1

1

0

к прибавлению ее 1 й строки ко 2 й строке.

 

0

0

1

 

26Первый индекс в паре – номер строки, к которой прибавляется строка, номер которой стоит вторым в паре.

ЛЕ К Ц И Я 1

Н.Н.БОБКОВ

ЛИНЕЙНАЯ АЛГЕБРА ◄ 24 ►

Замечание

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

Одной из таких наиболее часто используемых композиций является прибавление к строке ai

матрицы A другой ее строки,

aj , j

i , умноженной на некоторое число . Очевидно, при

0 описанная операция не меняет строк

A , а при 0 она представляет собой компози-

цию следующих элементарных преобразований A : (1) умножение aj на :

aj

aj , (2)

прибавление к строке ai

строки

aj :

ai ai aj , (3) умножение

aj

на 1 / :

aj (1 / )aj .

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

Это замечание, конечно, в равной степени относится и к столбцам матрицы.

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

ОБРАТНАЯ МАТРИЦА И ЕЕ СВОЙСТВА. ОБРАТИМЫЕ И НЕОБРАТИМЫЕ МАТРИЦЫ

Определение

 

Пусть для квадратной матрицы A M n в множестве M n

существует матрица,

обозначаемая A 1 27 и обладающая свойствами

 

(1.11)

 

 

 

 

AA 1 A 1 A En

.

 

 

Тогда матрица A называется обратимой, а A 1 матрица,

обратная по отноше-

нию к A .

 

27 A 1 нерасчленимое обозначение, символическое возведение в степень « 1», сходное с обозначением обратной функции.

Л Е К Ц И Я 1

ЛЕКЦИИ КАФЕДРЫ МАТЕМАТИКИ НИУ ВШЭ НН

25 ►

В

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

A M n

в множестве M n нет такой матрицы B , что A B B A En , матрица A назы-

вается необратимой. В этом случае будем писать: A 1 .

ВЛекции 3 будет показано, что обратимость матрицы связана с отличием от нуля

ееопределителя (детерминанта) det A важной числовой характеристики квадратных

матриц. Матрицу A M n называют вырожденной при det A 0 и невырожденной в про-

тивном случае, то есть если det A 0 . Таким образом, обратимость и невырожденность матрицы, а также необратимость и вырожденность оказываются равносильными понятиями.

Свойство (1.11) является прямым аналогом связи действительного числа a и об-

ратного к нему числа a 1 1a . Известно, что среди действительных чисел лишь одно, а

именно a 0 , не имеет обратного. Очевидно, нулевая матрица O M n , служащая в M n

аналогом числа 0

в множестве ,

необратима, поскольку для любой матрицы

B M n

имеем O B B O O En , так что

O 1 . Однако, класс квадратных матриц из M n , не

имеющих обратной, более широк и не сводится к одной лишь нулевой матрице.

 

Например,

пусть A M n идемпотентная матрица:

A2 A и существует

A 1 . То-

гда A2 A 1 AA 1 ,

или A( AA 1 )

AA 1 , так что необходимо A E

n

. Стало быть, единст-

 

 

 

 

 

 

 

En

En

 

 

 

 

 

 

 

 

 

 

 

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

Указанное обстоятельство имеет глубокое геометрическое истолкование. Действительно, выше уже отмечалось, что идемпотентные матрицы можно трактовать как некоторые проекторы. Но при проектировании, если оно не является тождественным преобразованием n мерного вектора в пространство n того же числа измерений28, безвозвратно теряется информация о некоторых координатах проектируемого вектора в том смысле, что по проекции невозможно однозначно восстановить его. В самом деле, одну и ту же проекцию на некоторую координатную ось или координатную плоскость имеет бесчис-

28 Общие определения понятий «вектор» и «число измерений» будут даны в Лекции 6, посвященной линейным (векторным) пространствам.

Л Е К Ц И Я 1

Н.Н.БОБКОВ ЛИНЕЙНАЯ АЛГЕБРА ◄ 26 ►

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

бы

как раз

умножением

координатного

столбца проекции слева на A 1 :

A 1

(A )

( A 1 A) E

 

 

. Именно из-за отсутствия такого обрат-

 

 

n

 

 

 

 

координатный

 

 

координатный столбец

 

 

столбец проекции

 

 

исходного вектора

 

ного отображения все идемпотентные матрицы в M n , кроме единичной, являются необ-

ратимыми.

Свойства обратной матрицы

1 .

Пусть матрица A M n

обратима и A 1 ее обратная матрица. Тогда (A 1 ) 1 A то

есть матрица A 1 тоже обратима, а ее обратная матрица совпадает с A .

 

Действительно, поскольку A 1 A En

и A A 1 En , то отсюда в соответствии с

определением (1.11) следует, что матрица A 1

имеет обратную матрицу, равную A . Тре-

буемое доказано.

 

 

Вывод: A и A 1 взаимно обратные матрицы.

2 . Пусть матрица A M n

обратима, A 1 ее обратная матрица и 0 . Тогда матрица

A тоже обратима, причем

( A) 1 1 A 1 .

 

В самом деле, учитывая, что AA 1 A 1 A En , находим

( A) ( 1 A 1 ) ( 1 )( AA 1 ) 1 En En ; ( 1 A 1 ) ( A) ( 1 )(A 1 A) 1 En En .

По определению это означает, что ( A) 1 1 A 1 .

3 . Пусть A M n и A 1 . Покажем, что тогда (AT ) 1 ( A 1 )T .

Вновь проверим выполнение равенств (1.11):

AT (A 1 )T

(A 1 A)T

ETn

En ;

( A 1 )T AT

(AA 1 )T

ETn

En .

Следовательно, ( AT ) 1 (A 1 )T , что и требовалось доказать.

Л Е К Ц И Я 1

ЛЕКЦИИ КАФЕДРЫ МАТЕМАТИКИ НИУ ВШЭ НН

27 ►

Итак, матрица, полученная транспонированием из обратимой матрицы, также обратима, причем ее обратная матрица получается транспонированием матрицы, обратной по отношению к взятой обратимой матрице. Матрицы AT и ( A 1 )T взаимно обратные.

4 . Пусть A и B обратимые матрицы, то есть A 1 , B 1 . Тогда и AB обратимая матри-

ца, причем обратной для нее служит матрица B 1 A 1 : ( AB) 1 B 1 A 1 .

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

( AB) (B 1 A 1 ) A (BB 1 ) A 1 A En A 1 ( A En ) A 1 AA 1 En ; (B 1 A 1 ) (AB) B 1 ( A 1 A) B B 1 En B (B 1 En ) B B 1 B En .

На основании (1.11) заключаем, что ( AB) 1 B 1 A 1 .

Как видно, обратимость матриц-множителей A, B влечет обратимость их произве-

дения AB и выполнение указанного равенства.

Обобщение свойства 4 . для любого конечного числа сомножителей выглядит так:

(1.11)’ ( A1 A2 Ak 1 Ak ) 1 Ak 1 Ak 11 A2 1 A1 1 (убедитесь в этом),

где Ai , i 1, k обратимые матрицы одного порядка.

Можно доказать (см. Лекцию 3), что обратимость произведения матриц влечет обратимость всех сомножителей. Равенства (1.11), (1.11)’ будут по-прежнему справедливы.

5 . Единственность обратной матрицы

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

Пусть для матрицы A M n существует две обратных матрицы, которые обозна-

чим как B1 и B2 . Составив произведение C B1 AB2 , усматриваем, что

C (B1 A)B2 En B2 B2 и C B1 ( AB2 ) B1En B1 .

Л Е К Ц И Я 1

Н.Н.БОБКОВ

ЛИНЕЙНАЯ АЛГЕБРА ◄ 28 ►

Поскольку равенство матриц транзитивно, то отсюда вытекает, что B1 B2 , так что

обратная по отношению к A матрица – единственна.

Замечания

 

1. Пусть A, X M n

Тогда если AX En ( XA En ), то и XA En ( AX En ). Мы вернем-

ся к доказательству этого утверждения в Лекции 3, когда сможем связать вырожденность квадратной матрицы с линейной зависимостью системы ее столбцов (строк). Пока же заметим, что сформулированное свойство позволяет оставить в определении обратной матрицы лишь одно из равенств (1.11) и использовать только его при проверке правильности решения задач об обращении матриц.

2. Если A M m , n ,

а B M n , m , m n , то может оказаться, что ровно одно из произведе-

ний AB или BA равно единичной матрице соответствующего порядка, то есть

 

AB E

 

 

AB E

m .

 

 

m , либо

 

BA En

BA En

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

Приведите пример матриц указанного вида.

3. В ряде случаев задачу обращения заданной матрицы A M n удается решить, не привлекая систематических средств обращения и критериев обратимости. Например, для матрицы A En равенство En B En , B M n выполняется для матрицы B En . Отсюда в силу единственности обратной матрицы вытекает обратимость единичной матрицы и

справедливость равенства

En1 ( B En

.

 

 

 

Далее, решая систему скалярных уравнений

 

n

ij ; i, j 1, n , равносильную

 

aik akj

 

 

 

k 1

 

матричному уравнению A 2 En , найдем все квадратные корни из En , то есть построим множество инволютивных матриц заданного порядка. Очевидно, попутно будет решена и

задача обращения таких матриц, поскольку A 2 En A 1 A .

Примечательно, что общее решение сформулированной задачи нетривиально и не имеет простого вида A En , аналогичного решению x 1 соответствующего алгебраи-

ческого уравнения x2 1.

Л Е К Ц И Я 1

ЛЕКЦИИ КАФЕДРЫ МАТЕМАТИКИ НИУ ВШЭ НН

29 ►

Проверьте, что

A1

7

8

 

E1/2

2

и

A2

7

48

 

E1/2

2 , то есть

A1 , A2 инволютив-

 

6

7

 

 

1

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Алгоритм вычисления обратной матрицы

Обсуждаемый алгоритм, который иногда называют методом исключения, позволяет для заданной матрицы A M n найти обратную матрицу A 1 или установить, что матрица A необратима. Этот метод тесно связан с определенными выше элементарными преобразованиями матриц. Как вскоре будет показано, для решения поставленной задачи обращения квадратной матрицы достаточно использовать только элементарные преобразования ее строк29.

ТЕОРЕМА(критерий обратимости квадратной матрицы)

Для того, чтобы матрица A M n была обратимой, необходимо и достаточно, что-

бы ее можно было преобразовать в единичную матрицу En (редуцировать к En ) при по-

мощи конечного числа элементарных преобразований строк.

Доказательство:

1. Пусть заданную матрицу A M n можно преобразовать в матрицу En посредством элементарных преобразований строк в соответствии с изображенной ниже схемой, суть которой состоит в следующем. Позиции на главной диагонали заданной матрицы последовательно, начиная с левого верхнего угла, заполняются ненулевыми числами, которые преобразуются в единицы, с дальнейшим планомерным обнулением или «исключением» всех стоящих под ними элементов30.

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

29Конечно, возможность практической реализации метода ручным счетом во многом зависит от размера матрицы и образующих ее чисел. Метод часто используется при решении задач учебного характера при n d7 для матриц с целочисленными элементами.

30Аналогичная технология «исключения» применяется в методе Гаусса решения систем линейных уравнений, см. подробности в Лекции 2 и далее.

ЛЕ К Ц И Я 1

Н.Н.БОБКОВ

ЛИНЕЙНАЯ АЛГЕБРА ◄ 30 ►

Переход от каждой матрицы предыдущего шага к следующей, полученной из нее элементарным преобразованием строк, отмечается на схеме значком « »:

1 1 1

 

 

 

 

 

 

 

0

 

 

 

0 1

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

1

 

 

0

 

(1.12)

0

1

 

 

 

0

1

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

0

 

 

 

 

0

0

1

 

0

0

1

 

 

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

En .

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поясним этапы алгоритма приведения матрицы A к En более подробно:

создаем 1 на 1-м диагональном месте (индексы i j 1 ). Для этого можно переставить строки (если среди них есть начинающаяся с 1), прибавить к некоторой строке другую, умноженную на подходящее число или умножить первую строку на число, обратное ее первому элементу, если он отличен от нуля. При ручном счете рекомендуется избегать таких делений, осложняющих выкладки, до тех пор, пока это возможно;

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

создаем 1 на втором диагональном месте (индексы i j 2 );

Л Е К Ц И Я 1

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