Ортогональная матрица
Квадратная матрица А называется ортогональной, если
ААТ=АТА = Е, (1.59)
т.е. А - ортогональна, если ее транспонированная матрица является обратной исходной. Таким образом, каждая ортогональная матрица обратима.
Так как (/1Г) = Д, то из (1.59) вытекает, что обращение ортогональной матрицы есть ортогональная матрица.
Примеры ортогональных матриц:
/ cos a sin сп / 0,8 — 0,6 /1 0
V-sina cos а/ 0,6 -0,8/’ 17
Отметим, что единичная матрица любого порядка является ортогональной.
Справедливы следующие утверждения
1. Произведение ортогональных матриц есть ортогональная матрица. Действительно, пусть матрицы А Ортогональны, тогда
Ат = А'1, Вт = В~' и, значит,
(АВУ = ВТАТ = В-1 А-1 = {АВ}-1 (1.60)
Формулы (1.60) доказывают утверждение.
2. Произведение двух ортогональных матриц есть ортогональная матрица.
Докажем, пусть А иВ - ортогональные матрицы одного порядка. Так как (АВ)Т =ВТАТ, то (АВ)г(АВ) = ВгАгАВ = Вт(а'А}В = ВгЕВ = ВтВ = Е, т.е. (АВ) т (АВ) = Е. Следовательно, АВ - ортогональная матрица.
Равенство Д7 - А1 выражает необходимое и достаточное условие ортогональности матрицы А.
В самом деле, если А - ортогональная матрица, то ДГД = Е. С другой стороны, А-1 А = Е. Из этих двух равенств следует, что Ат = А 1.
Обратно, пусть Ат - А1, тогда отсюда и из равенства А1 А = Е следует Ar А = Е, т.е. А - ортогональная матрица.
3. Матрица, полученная транспонированием ортогональной матрицы -ортогональна.
Действительно, если А - ортогональная матрица, то (иг) Ат = ААТ = АА-1 = Е, или (аг^Ат =Е. Но это означает, что Ат -ортогональная матрица.
4. Матрица, обратная ортогональной матрице, - ортогональна.
В самом деле, если А - ортогональная матрица, то Ат = А’, AJA-E. Далее, (д”7) А~1=АА~1=Е, или (д”7) А~]=Е, что означает, А~] -ортогональная матрица.
Определитель матрицы
Определитель - числовая характеристика квадратной матрицы. Для того, что бы обусловить естественную природу такой числовой характеристики, чуть- чуть забежим вперед, и решим в общем виде систему двух уравнений с двумя переменными, известным нам еще со школы, методом исключения неизвестных.
И так, пусть задана система:
(аг1х + а12у = Ьг
{а21х + а22у = Ь2
Исключим переменную у. Для этого умножим первое уравнение на а22, второе на (—а12):
Г а.ц(і22х 4- сі^2а.22у — Ьга22
(—а21а12х — сі22а12у = ~^2аі2
Теперь сложим уравнения:
(а11а22_а21а12)х = ^1а22 — ^2а12
И выразим переменную х:
^1«22 ^2 «12
«11«22—«21«12
Далее, в точности так же, исключим переменную у. Первое уравнение умножаем на (~а12), второе на а1г. И дальше опять складываем уравнения.
С—а11а21х — а12а21у — —Ь±а21 I а1±а21х + аца22у = Ь2аг1
у(а11а22~а12а21) — ^2а11 - ^1«21 ^2 «11 ~ ^1а21 а11а22~а12а21
Проанализируем, что получилось. Знаменатели дробей (для х и у) одинаковы. Если построить матрицу из коэффициентов системы:
Л = /«11 «12
К«21 «22/
то знаменатель есть произведение элементов матрицы, стоящих на главной диагонали минус произведение элементов, стоящих на побочной диагонали. Вот эту структуру и назвали определителем (detriminant).
Обозначение:
detA^11 а2
1«21 «221
Таким образом, имеем
|«11 «121
1«21 «22 I “ 61116122 “ «21«12-
Далее, проанализируем числители для х и для у. И тут такая же структура, т.е. определители. Для переменной х, обозначим его как Дь для у- как Д2.
Тогда л ^1 «12 , ,
^1_ h л ~ ^1а22 “ ^2«12 «2 «22
Л «11 ^1 К. I
^2— л к — «11^2 “ «21^1 «21 «2
И так, мы вышли на понятие определителя и вдобавок вывели формулы Крамера (об этом подробнее позже).
Определителем (детерминантом) квадратной матрицы второго порядка (или - определителем второго порядка) называется число, равное разности произведений элементов, стоящих на главной диагонали матрицы и произведений элементов побочной диагонали:
Рисунок 16
Для определителя матрицы А употребляются следующие обозначения:|Л|, A, detA, det(aij').
Аналогично, рассмотрев уже систему трех уравнений с тремя переменными, приходим к определителю третьего порядка (попробуйте самостоятельно):
Рисунок 17
Правило вычисления определителя по этой схеме называют правилом Сюрроса.
Для того чтобы определиться со знаками слагаемых при вычислении определителя произвольного, п — го порядка нам необходимо познакомиться с элементарными сведениями теории перестановок.
Перестановкой из п натуральных чисел 1, 2, 3, ..., п называется любое их расположение, взятое в определенном порядке.
Произвольную перестановку из п чисел будем записывать в виде (а1Ла2, где каждое - одно из чисел 1,2, ..., п и cq Ф о, при і ф j.
Две перестановки (ах, а2>..., ап) и (fii>P2,называются равными, если «і і при/= 1, и. В противном случае перестановки называют неравными (различными).
Число различных перестановок из чисел 1, 2, 3, ..., п, равно произведению 1 ? 2 • 3 • ... • (n—1) -п, которое обозначается п! («эн факториал»).
Говорят, что два числа образуют инверсию в перестановке, если бо'лыпее число стоит левее ме'ныпего.
Обозначение: k(alt а2,..., ап) - число инверсий в перестановке (ai,a2,...,an).
Например, рассмотрим перестановку (3,4,1,2,7,10). Перечислим пары чисел, образующих инверсию: (3,1), (3,2), (4,1), (4,2).
То есть, к (3,4,1,2,7,10) = 4.
Способ вычисления числа инверсий основан на утверждении. к(а1,а2, ...,ап) = к± + к2+... +/сп-і,гДе fc/(i = l, л-1) - число чисел, стоящих перед числом «j в перестановке, полученной из данной зачеркиванием чисел, меныпихар
Например, найдем используя это утверждение, количество инверсий в перестановке
(5,1,3,2,4,6)
Перед числом 1 стоит одно число, следовательно, кг = 1.
Зачеркивая единицу, получаем
(5,4-,3,2,4,6)
Дальше перед числом 2 стоит два числа, значит к2 = 2. И зачеркиваем двойку:
(5,4,3,2,4,6)
Перед числом 3, одно число, то есть к3 = 1, и зачеркиваем 3.
(5,4,3,2,4,6)
Перед числом 4 стоит одно число, значит к4 = 1, и зачеркиваем четверку.
(5,4,3,2,4,6)
Дальше перед пятеркой нет чисел к5 = 0:
(5,4,3,2,4,6)
И к6 = 0.
Таким образом,
fc(5,l,3,2,4,6) = к1 + к2 + к3 + к4 + к5 + к6 = 142 + 1 + 1 + 0 + 0 = 5.
Очевидно, если числа в перестановке записаны в порядке возрастания, то число инверсий равно нулю.
Если число инверсий в перестановке четное, то она называется четной, в противном случае - нечетной.
Считают, две перестановки имеют разный характер четности, если одна из них - четная, а другая - нечетная.
Очевидны следующие утверждения:
- 1. Пусть в некоторой перестановке сделана транспозиция. Она равна произведению нечетного числа транспозиций соседних элементов.
- 2. При транспозиции соседних элементов число инверсий в перестановке меняется на одну единицу.
- 3. Если в перестановке сделать транспозицию соседних элементов, то четность перестановки изменится на противоположную.
- 4. Любая транспозиция изменяет четность перестановки на противоположную.
- 5. Число четных перестановок п элементов равно числу нечетных перестановок.
- 6. Любая перестановка может быть получена из любой другой посредством нескольких транспозиций.
Определитель п-го порядка
Определителем квадратной матрицы называется алгебраическая сумма всевозможных произведений элементов этой матрицы, взятых по одному из каждой строки и по одному из каждого столбца. Сомножители в каждом слагаемом записываются в порядке следования строк, номера же столбцов образуют перестановки; слагаемые, соответствующие четным перестановкам, берутся со знаком «плюс», соответствующие нечетным - со знаком «минус».
|
Дл |
я определителя третьего |
порядка, в соответствии с этим |
|
определением: «11 «12 ... «1и |
= ^-^rWaljl - а2І2 - ...? а„1п (/) |
|
|
с |
ІеіЛ = а21 а22 а2п О-пІ ^п2 ^пп «11 «12 «13 |
|
|
а21 а22 а23 = У азі аз2 азз (aitZ2 ап |
(-1)*^.....an)alaia2a2 ...апап, ) |
|
где (alt а2, ...,ап) пробегает все перестановки чисел 1, 2,..., п.
Далее, множитель (—l)fc(ai-a2. ".«п) равен +1, если (alt а2,ап) четная перестановка, и равен - 1, если нечетная.
Понятие определителя имеет смысл для матриц с элементами из любого ассоциативного коммутативного кольца и, в частности, из любого поля.
Вычислить определитель высокого порядка, пользуясь определением, задача не из легких. Поэтому на практике вначале преобразуют определитель на основании свойств, затем пользуются теоремой Лапласа или одним из методов, предлагаемых вашему вниманию далее.
Теорема Лапласа
Введем несколько понятий.
(Цц а12 ... а1п
^21 а.?.2 J ? квадратная матрица порядка п.
ап1 ап2 ••• апп/
Минором порядка к для этой матрицы называется определитель матрицы, составленной из элементов, находящихся на пересечении некоторых выбранных к строк и к столбцов:
аагрг аагрк
а«к01 ?" аЯк0к
Минором, дополнительным к данному минору порядка к, называется минор порядка п — к, матрица которого получается из исходной посредством вычеркивания строк и столбцов, содержащих данный минор порядка к.
Алгебраическим дополнением к данному минору называется дополнительный минор С множителем (—l)ai+ "+ak+/?i+ - +/?fc.
Теорема 1.25 (Лапласа)
Пусть в матрице определителя выбраны к строк. Определитель равен сумме произведений всех миноров порядка к, составленных из этих строк, на их алгебраические дополнения.
Например, вычислим определитель
теорему Лапласа
четвёртого порядка, используя
Д=
- 1 -1
- 2
- 2
- 2
- 2 О
- 1
- 3
- 1 -1
- 1
О
- 1
- 3 -1
Выберем любое, например, две строки. Пусть это будут первая и вторая строка.
Тогда
Ж -1І-ІЛ И ДНА И 71-
И ДГІ2 71-1? ?1-1? ?1=
= (2 + 2)-(1 — 3) — (И-3)-(0 — 3) + (1 — 0)-(0 + 1) — (2 — 6)-(-2 - 6) + (2 - 0) ? (2 + 2) - (3 - 0) • (2 - 0) = 4 ? (-2) - 4 ? (-3) +11- (-4) ? (-8) + 2 ? 4 — 3 ? 2 = —8 + 12 + 1 — —32 + 8 — 6 = 7.
Частный случай теоремы Лапласа
Рассмотрим ступенчатую матрицу А:
|
/ «и «21 |
«12 «22 |
|
|
|
||
|
А = |
«ml |
«m2 |
- «тт |
0 |
0 |
|
|
ат+1,1 |
«т+1,2 |
"• «т+1,т |
«m+lrn+І |
«т+1,л |
||
|
ап1 |
«л2 |
- «лт |
«и,т+1 |
«тт , |
||
|
Если к определителю |
ступенчатой матрицы |
применить |
||||
|
Лапласа, то |
||||||
|
«11 |
«12 |
- «14 |
«т+1,т+1 |
«т+1,п |
||
|
det А = |
«21 |
«22 |
«24 . |
|||
|
«ml |
«m2 |
«тт |
«п,т+1 |
«п,п |
||
теорему
Например, вычислим определитель 7-го порядка.
|
1 |
2 |
0 0 |
0 |
0 |
0 | ||||||||||
|
3 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
3 |
0 |
0 |
||||
|
7 |
8 |
0 |
1 |
3 |
0 |
0 |
11 |
9 і |
2 |
1 |
0 |
0 |
0 |
||
|
1 |
3 |
2 |
1 |
0 |
0 |
0 |
4 ' |
3 |
0 |
2 |
0 |
0 |
|||
|
5 |
6 |
3 |
0 |
2 |
с |
1 0 |
3 и |
5 |
8 |
6 |
1 |
0 |
|||
|
3 |
7 |
5 |
8 |
6 |
1 |
0 |
0 |
2 |
-1 |
3 |
5 |
||||
|
1 |
-1 |
0 |
2 |
-1 |
3 |
! 5 |
|||||||||
|
0 |
1 |
3 |
|||||||||||||
|
=1з |
4І |
|
1 0 |
|
•Із |
^1 |
|= |
-2- |
(- |
13) ? |
?5 = |
130. |
|||
На основании теоремы Лапласа можно утверждать, что поскольку любой определитель и-го порядка сводится к вычислению определителей, например, второго или третьего порядков, то все свойства для определителей второго и третьего порядков справедливы и для определителя любого порядка.
Введем еще понятия.
Минором элементна^ заданного определителя называется определитель, полученный из данного вычеркиванием той строки и того столбца, которым принадлежит данный элемент. Минор элемента aik обозначается Mik.
Например, минором элемента а21 определителя
«и
«21
«31
«12 «13
(1-61)
«22 «23
«32 «33
является определитель второго порядка а,-, Мгї = 12 13
«32 «33
Алгебраическим дополнением злементааік заданного определителя называется его минор, взятый со знаком (—l)i+k и обозначается т.е.
Так, алгебраическим дополнением элемента а21 определителя (1.61) является определитель М2/, взятый со знаком минус, т.е.
«12 «13
Л7=(-1)2+1М27=-
«32 «33
Далее, мы рассмотрим свойства определителей. Доказательство некоторых из них предлагаем читателю произвести самостоятельно, используя для этого определители третьего порядка. Для чего вам необходимо будет только лишь в соответствии с условием, записать определитель, раскрыть его, используя правило Сюрроса, и сравнить с исходным (неизмененным определителем.)