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

ЛЕКЦИЯ 01_Л.А

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

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

11 ►

 

 

 

 

 

 

 

 

p

вынесен за знак суммирования по этому индексу, а сумма (bkj cjv ) есть элемент произве-

 

 

 

 

 

 

 

 

j 1

дения B на C с индексами « k, v », так что

 

 

n

p

 

n

 

 

 

 

 

A(BC) u v

, u 1, m; v 1,l .

 

 

 

 

 

au k (bkj cjv ) au k (BC)kv

k 1

j 1

 

k 1

 

 

 

 

 

 

 

 

В соответствии с определением (1.4) отсюда вытекает доказываемое матричное ра-

венство:

 

.

 

 

 

 

( AB)C A(BC)

 

 

 

 

 

 

 

 

 

 

 

p

n

Убедитесь непосредственным вычислением в том, что .

 

 

 

 

 

 

 

j 1

k 1

2 . A M m , n ; B,C M n , p A(B C) AB AC –умножение на матрицу слева дистрибутивно по отношению к сложению матриц (докажите самостоятельно).

Это же касается и умножения справа:

3 . A, B M m , n ; C M n , p ( A B)C AC BC (докажите самостоятельно).

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

мутативно. В самом деле, если произведение AB определено, то в общем случае произведение BA не только не равно AB , но даже может не существовать.

Пример:

A M 7 , 5 , B M 5 ,10

AB M 7 ,10 и BA , т.к. 10 7 15.

Определение

 

Матрица KAB

AB BA называется коммутатором матриц A и B . Если KAB O,

то матрицы A и B называются коммутирующими, или перестановочными.

Докажите, что KAB A, B M n , так что коммутировать могут только квадратные матрицы одинакового порядка. Приведите примеры таких матриц.

Задайте матрицу A M 2 и найдите все матрицы, перестановочные с ней.

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

Л Е К Ц И Я 1

Н.Н.БОБКОВ

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

ТРАНСПОНИРОВАНИЕ МАТРИЦЫ

Определение

Для всякой матрицы A M m , n транспонированная по отношению к ней матрица

C определяется условием

C AT ci j aji , i 1, n; j 1, m 16.

Из приведенного определения вытекает, что C AT M n , m , причем всякая строка

(всякий столбец) транспонированной матрицы состоит из элементов столбца (строки) матрицы A с тем же номером в их естественном порядке.

Примеры:

1). A 1

2

AT

1

 

 

 

 

 

 

.

 

 

 

 

 

(1 2 )

 

 

 

 

( 2 1)

2

 

 

 

 

 

 

 

a b

 

c

BT

 

a

d

 

 

 

2). B

 

 

b

e .

 

 

 

( 2 3)

d e

 

f

 

(3 2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

f

 

 

 

 

1

0

 

2

 

 

1

0

2

 

 

3). M 0

1

 

4 M T

0

1

4 M .

 

 

 

2 4

 

3

 

 

2

4 3

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

x

 

вектор-столбец размеров (n 1) . Тогда xT

x1 x2

xn вектор-

4). Пусть x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

строка размеров (1 n) . Ясно, что определено матричное умножение xT

на x , равное

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

n

 

 

xT x x1

 

 

xn x2

 

 

 

 

x2

 

x12 xn2 xk2 . Для n 2,3 эта сумма в соответствии

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

с теоремой Пифагора выражает квадрат длины вектора с декартовыми прямоугольными координатами x1 , x2 ; x1 , x2 , x3 соответственно.

16Транспонированная матрица иногда обозначается посредством A .

ЛЕ К Ц И Я 1

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

13 ►

Обобщая теорему Пифагора на случай произвольного числа координат, говорят о длине или так называемой евклидовой норме любого вектор-столбца x(n 1) , n :

 

 

 

n

1/2

 

 

(1.9)

x

 

xk2

 

(xT x)1/2

.

 

 

k 1

 

 

 

СВОЙСТВА ОПЕРАЦИИ ТРАНСПОНИРОВАНИЯ

1 . ( AT )T A (докажите).

2 . A B T AT BT .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть A,

B,C A B M m , n . Тогда CT M n , m

и

 

 

 

 

 

 

 

 

 

 

 

T

 

T

 

 

 

 

 

 

 

 

 

T

 

T

 

 

T

B

T

)i j .

 

( A B)

ij (C

 

)i j cji ( A B) ji aji bji ( A )i j (B

 

)i j ( A

 

 

 

 

Итак, (C

T

)i j

 

 

T

 

T

B

T

)i j , i 1, n ,

j 1, m , а это по определению оз-

 

 

( A B)

 

 

ij

( A

 

начает, что

( A B)T

AT BT

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 .

( AB)T BT AT .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

Пусть A M m , n , B M n , p AB M m , p ,( AB)i j aik bkj

, i 1, m ,

j 1, p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

Далее, BT M p , n ,

AT M n , m , так что определено произведение

BT AT M p , m ,

элементы которого выражаются следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

(BT AT )u v

(BT )uk ( AT )kv

bku avk

avk bku ( AB)vu

[(AB)T ]uv

 

для всех u 1, p ;

v 1, m .

 

 

k 1

 

 

 

 

k 1

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По определению равенства матриц это означает, что

( AB)T

BT AT

.

 

Л Е К Ц И Я 1

Н.Н.БОБКОВ

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

Замечание

Доказанное только что свойство операции транспонирования легко обобщается на любое число сомножителей. Так, например, ( ABC)T CT BT AT и т.п. (докажите).

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

СЛЕД КВАДРАТНОЙ МАТРИЦЫ

Определение

В квадратной матрице {ai j } A Mn совокупность элементов, индексы которых удовлетворяют уравнению i j , образует так называемую главную диагональ, идущую из левого верхнего угла в правый нижний. Совокупность элементов, для которых i j n 1, образует побочную диагональ, которая идет из правого верхнего угла в левый

нижний17:

главная диагональ

A

( n n )

a11

a21

a31

a n 2,1

a n 1,1

a n1

a12

a13

a1, n 2

a1, n 1

a1n

a22

a23

a 2, n 2

a2, n 1

a2n

a32

a33

a3, n 2

a3, n 1

a3n

 

 

 

 

 

 

 

 

 

a n 2,2

a n 2,3

an 2,n 2

an 2,n 1

a n 2,n

 

a n 1,2

a n 1,3

a n 1,n 2

a n 1,n 1

a n 1,n

 

 

an2

an3

a n,n 2

a n,n 1

an n

побочная диагональ

Определение

Следом квадратной матрицы называется сумма элементов ее главной диагонали:

n

tr A aii (обозначается также Sp A от «Spur» – след по-немецки.).

i 1

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

Л Е К Ц И Я 1

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

15 ►

Ниже приводятся основные свойства этой важной характеристики числовой квад-

ратной матрицы. Предполагается, что A, B M n .

Свойства следа матрицы

 

n

n

n

1 .

tr ( A) ( A)ii ( aii ) aii tr A при умножении матрицы на число ее след

 

i 1

i 1

i 1

также умножается на это число.

2 . tr ( A B) tr A tr B (докажите).

3 . tr AT tr A транспонирование квадратной матрицы, представляющее собой симметрич-

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

 

 

 

 

 

n

 

n

 

4 . tr ( AB) tr (BA) . В самом деле, ( AB)ij aik bkj ( AB)ii aik bki , а также (BA)ij

 

 

 

 

 

k 1

 

k 1

 

n

 

 

n

 

 

 

 

 

bik akj (BA)ii bik aki .Суммируя теперь

диагональные

элементы обеих матриц,

k 1

 

k 1

 

 

 

 

 

 

 

n

 

 

n n

n

n

n

получаем: tr ( AB) ( AB)ii aik bki , tr (BA) (BA)ii bik aki (i k)

 

 

i 1

 

 

i 1 k 1

i 1

i 1 k 1

n n

 

n

n

 

n n

 

 

 

bki aik

 

 

 

 

aik bki = tr ( AB) .

 

 

k 1 i 1

i 1

k 1

 

i 1 k 1

 

 

 

5 . tr ( AT A) tr (AAT ) 0 .

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

 

 

 

n

n

B AT . Далее, используя сделанные выше выкладки, получим tr (AT A) (AT )ik aki

 

 

 

i 1

k 1

n n

n

 

 

 

aki2

aki2

0 . Если элементы матрицы

A вещественные числа,

то равенство

i 1 k 1

i,k 1

 

 

 

здесь имеет место в том и только том случае, когда aki 0 ; k,i 1, n A O M n (до-

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

матрицы).

Л Е К Ц И Я 1

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

СПЕЦИАЛЬНЫЕ ВИДЫ КВАДРАТНЫХ МАТРИЦ

Определение

Квадратная матрица A называется:

Симметрической

 

 

если

AT

A

.

 

Кососимметрической (антисимметрической)

 

 

 

 

 

 

 

если

AT

A

.

Нижнетреугольной

если

 

при

 

.

ai j 0

i j

 

 

 

 

 

 

 

 

 

0

все элементы над главной

 

 

 

 

 

 

диагональю – нули

 

 

 

 

 

 

главная диагональ

 

 

Верхнетреугольной

если ai j 0

при i j .

0

все элементы под главной диагональю – нули

главная диагональ

 

 

Диагональной

если ai j 0

при i j .

0

0

все элементы вне главной диагонали – нули

главная диагональ

Запись B diag (1, 1,0,7) означает диагональную матрицу с перечисленными в скобках элементами, образующими ее главную диагональ:

1

0

0

0

 

diag (1, 1,0,7)

0

1

0

0

.

 

0

0

0

0

 

 

0

0

0

7

 

Л Е К Ц И Я 1

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

17 ►

Частные случаи диагональной матрицы – это нулевая матрица O M n и так назы-

 

 

 

1

0

 

0

 

 

 

 

 

ваемая единичная матрица:

En diag (1, ,1)

0

1

 

0

 

18.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n единиц

 

 

 

 

 

 

 

 

 

0

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Элементы единичной

матрицы в индексной

форме:

i j

1, i j

, где i 1, n ,

 

j

 

 

 

 

 

 

 

 

 

0, i

 

j 1, n . Объект i j называют символом Кронеккера.

 

 

 

 

 

 

 

Очевидно, матрица En

обладает свойством En

A

A En A и потому играет в

 

 

 

 

 

( n n )

 

 

 

 

 

 

множестве M n квадратных матриц порядка n ту же роль, которую играет в множестве

действительных чисел число 1. Роль числа 0

принадлежит в M n нулевой матрице, по-

скольку A O O A O.

 

 

 

 

Ортогональной

 

если

A AT AT A En

.

 

 

1, i j

.

Докажите, что если si , s j строки ортогональной матрицы, то si sTj i j

j

 

 

0, i

 

Говорят, что система строк, удовлетворяющая написанному условию, – ортонормирована: норма (1.7) каждой из строк равна 1, а произведение si sTj для не равных друг другу значений индексов i , j (аналог скалярного произведения векторов в геометрии) равно 0 . Напомним, что таким свойством обладают базисные векторы д.п.с.к.19, обозначае-

мые i , j и i , j , k соответственно в двумерном и трехмерном случаях.

Заметим, что прямоугольная матрица A M m , n , m n может обладать ровно од-

ним из свойств AAT Em или AT A En 20. Такие матрицы называются полуортогональ-

ными.

Дайте пример полуортогональной матрицы размера (2 3) .

18Обозначается также посредством In .

19Декартова прямоугольная система координат.

20Иначе говоря, обе матрицы AAT и AT A определены, но лишь одна из них – единичная. Можно доказать, что из этих двух квадратных матриц единичной всегда будет матрица меньшего порядка.

ЛЕ К Ц И Я 1

Н.Н.БОБКОВ

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

Как вы думаете, будет ли система столбцов ортогональной матрицы ортонормированной? Обоснуйте свою догадку.

Матрица B , удовлетворяющая равенству B2 A , где A, B M n , называется квадрат-

ным корнем из матрицы A и обозначается A1/2 . Матрица En1/2 , то есть квадратный корень из единичной матрицы, называется инволютивной матрицей.

Оказывается, что не для всех матриц A существует квадратный корень21

A1/2 , а ес-

ли существует, то он необязательно единственный.

 

 

 

 

 

 

1

0

 

2

1

1/2

Извлеките квадратный корень из матрицы A 0

2

. Существует ли

1

2

 

?

Матрицей перестановки22

 

если в любой

ее

строке, а

также в любом ее столбце все элементы, кроме одного, равного 1, нули.

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

Пример:

0 0

1

 

1

0

0

 

1 0

0

 

 

 

1

0

0

 

 

 

0

0

1

 

 

 

0

1

0

 

E3

,

 

 

 

 

 

 

 

0

1

0

 

r1 r2

 

0

1

0

 

r2 r3

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

или

0

0

1

 

1

0

0

 

1 0

0

 

 

 

 

 

1

0

0

 

 

 

0

0

1

 

 

 

0

1

0

 

E3

.

 

 

 

 

 

 

 

 

 

 

 

0

1

0

 

c1 c3

 

0

1

0

 

c2 c3

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Идемпотентной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если

A2 A

.

21Так было уже для действительных чисел.

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

ЛЕ К Ц И Я 1

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

19 ►

Пример:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

1

0

 

0

 

 

 

 

 

 

 

 

 

 

Рассмотрим вектор-столбец

b

 

и матрицу P

 

0

0

 

0

. Тогда

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

1 0

0

 

a

 

a

 

 

 

 

 

 

 

1

0

0

 

1

0

0

 

1

0

0

 

 

 

 

b

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

0 0

0

 

 

0

 

x

и

P2

 

0

0

0

 

0

 

0

0

0

0

0

P

, так что

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

0 0

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

0

0

 

 

 

0

0

0

 

 

 

0

0

0

 

 

 

 

 

c

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Px идемпотентная матрица. Глубокий смысл этих соотношений состоит в следующем: ес-

ли a,b, c координаты некоторого геометрического вектора u в выбранной системе коор-

динат Oxyz , то a,0,0 это координаты вектора ux , являющегося прямоугольной проекцией u на ось абсцисс. Таким образом, умножение матрицы Px на координатный столбец произвольного вектора u реализует преобразование этого вектора в его проекцию на указанную ось. Очевидно, повторное проектирование вектора ux даст вновь ux , поскольку этот

вектор уже лежит на оси Ox . В самом деле,

P

x

 

 

P (P ) (P

P ) P2

P

 

 

 

 

 

 

 

x

 

 

 

 

x

x

x

 

x

 

x

x

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обратно, если для любого вектора u

выполняется условие

 

 

, то матрица

 

Px

x x

Px идемпотентная. Действительно, Px x x

Px2

Px (Px2

Px ) O(3 1)

(0 0 0)T , откуда в силу произвольности вектора

u

вытекает,

 

что

Px2 Px

O(3 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(обоснуйте это утверждение), или

P2

P

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из идемпотентности матрицы Px следует, что любая ее натуральная степень совпа-

 

Px Pxn 1

 

 

дает с ней самой:

Pxn

Pxn 2 Px2 Pxn 2

 

Px

, n . Геометрически это озна-

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

 

 

1

0

0

 

 

По аналогичному принципу умножение матрицы P

 

0

1

0

 

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

xy

 

 

 

 

 

 

 

 

0

0

0

 

 

столбец произвольного вектора u реализует его прямоугольное проектирование на координатную плоскость Oxy :

Л Е К Ц И Я 1

Н.Н.БОБКОВ

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

 

 

1

0

0

 

a

 

a

 

 

 

P

 

0

1

0

 

b

 

 

b

 

 

xy

. И здесь не составляет большого труда убедиться в

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

 

c

 

0

 

 

 

справедливости равенства Px2y Pxy , означающего идемпотентность матрицы Pxy .

Преобразование вектора u в себя также можно считать проектированием в пространство Oxyz 3 . Очевидно, что оно реализуется единичной матрицей E3 , которая,

очевидно, идемпотентна: n E2n En .

Описанные преобразования векторов геометрического пространства называют опе-

раторами проектирования или проекторами23. В линейной алгебре проектором назы-

вают также произвольную идемпотентную симметрическую матрицу.

ЭЛЕМЕНТАРНЫЕ ПРЕОБРАЗОВАНИЯ МАТРИЦ И ИХ РЕАЛИЗАЦИЯ ПРИ ПОМОЩИ МАТРИЧНОГО УМНОЖЕНИЯ

Определение

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

(1.10)

1.Перестановка (транспозиция) двух строк матрицы.

2.Умножение строки матрицы на число, отличное от нуля24.

3.Сложение строки матрицы с некоторой другой ее строкой.

4.Те же операции со столбцами матрицы.

Действия 1 – 3 будем также называть элементарными преобразованиями строк матрицы, а действия 4 – элементарными преобразованиями ее столбцов. Как будет видно из дальнейшего изложения, элементарные преобразования играют важную роль во многих разделах линейной алгебры.

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

24Умножение в этом пункте и сложение в следующем, как всегда, поэлементные.

ЛЕ К Ц И Я 1

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