книги / Матричные и асимптотические методы в теории линейных систем
..pdfопределитель которой отличен от нуля, во-вторых, опре делитель обратной матрицы тоже отличен от нуля и, в-тре тьих, определители матриц А и Л-1 — взаимно обратные числа:
и- I—TV
Итак, обратную может иметь только невырожденная мат рица.
Всякая невырожденная матрица имеет обратную. В са мом деле, как это следует из равенств (5.1), обратной для матрицы А является, например, матрица
A -l = -fa -A e. |
(6.3) |
Каждая матрица может иметь только одну обратную. Действительно, пусть матрица С удовлетворяет, как и мат
рица Л-1» условиям (6.1). Тогда
С = СЕ = С (АА~') = (СА) А~' = ЕА~1 = А~\
что и доказывает единственность обратной матрицы. Матрица, обратная произведению двух невырожденных
матриц Л и В, определяется равенством
(.АВГ' = Вг'А~'. |
(6.4) |
В самом деле, умножая обе части равенства
(ЛВ)"1АВ = Е
справа на произведение обратных В-1Л-1, сразу получаем (6.4).
§7. Транспонирование матрицы
ипереход к сопряженной матрице
Рассмотрим |
прямоугольную |
матрицу Л = |
(atj) с |
раз |
|
мерами т X п. |
|
(а]/) с размерами п X т называется тран |
|||
Матрица Л' = |
|||||
спонированной |
по |
отношению |
к матрице Л, |
если |
аы = |
aik‘ |
|
' |
|
|
|
Матрица Л* = (а*/) с размерами' п X т называется со пряженной (или эрмитово сопряженной) по отношению к
матрице Л, если аы = aik, где а(к — число, комплексно
сопряженное элементу aik. В частности, для скалярной ве личины а, которую можно рассматривать как матрицу с размерами 1 X 1, я* — ~а.
Операции транспонирования и перехода к сопряженной матрице обладают следующими легко доказываемыми свой ствами:
1) |
(А')' = А, |
(А*)* = А, |
|
2) |
(А + В )' = А' + В', |
(А + В)* = |
А* + В*, |
3) |
(аЛ)' = аА', |
(осА)* = |
осА*, |
4) |
(АВУ = В’А \ |
(АВ)* = В*А*. |
|
Кроме того, для квадратных матриц |
|||
5) |
(Л” 1)' = |
( Л 'Г \ |
(Л-1)* = (Л*)-1, |
6) |
det А' = |
det A, |
det А* = det А. |
Если матрица А совпадает со своей транспонированной (А = А'), то матрица А называется симметрической. Ес ли же А совпадает со своей сопряженной (А = Л*), то она называется эрмитовой. Симметрическими и эрмитовыми матрицами могут быть только квадратные матрицы.
§ 8. Блочные матрицы
Прямоугольную матрицу
ап |
#12 |
а\п |
#21 |
#22 |
#2п |
Р-т\ |
U-m2 |
атп |
горизонтальными и вертикальными линиями можно рас сечь на прямоугольные клетки (блоки):
( s < m , / < л ) .
Каждый из блоков (субматриц) Ац представляет собой некоторую прямоугольную матрицу (и, в частности, число)
с размерами mt X nt; например,
где
А = |
« п |
: fll2 |
«13 |
«14 |
«2J • «23 |
«23 |
«24 |
||
|
^4со СЗ |
'■«82 |
«33 «34i |
|
^11 «11» |
^ 1 2 |
— («12 |
«13 а и ) > ^ 2 1 |
д |
__ |
/« 2 2 |
«23 |
«24 \ |
|
22 “ |
U . |
«зз |
« s i / |
В частности, матрица может быть рассечена только го ризонтальными или только вертикальными линиями. При этом блочная матрица будет иметь соответственно вид
А = |
или А = {Аг As |
А,). |
Сокращенно блочную матрицу обозначают так:
А — (Лар)5>*.
Рассмотрим две матрицы А и В одинаковых размеров и с одинаковым разбиением на блоки, т. е.
|
Л. |
А 12 |
|
/В ц |
В13 |
в и |
А = |
Л21 |
Л2а |
|
в = | Ви |
|
B^t |
|
|
|
|
|
||
|
,Л51 |
A S2 |
А и / |
\ в ,, |
Bs2 |
К |
и матрицы Лсф и В а р |
имеют одинаковые |
размеры |
т а X лр |
|||
(а = 1, 2, ..., |
s; р = 1, 2, ..., /). Тогда в соответствии с пра |
|||||
вилом сложения матриц |
|
|
|
|||
|
|
( ^ и + |
^ 12 + |
В13 |
+ |
&\ty |
|
|
^ 21 ~Ь ^ 2 1 ^ 2 2 "Н ^ 2 2 |
|
|
A s \ -}- В 51 А$2 -f- B S2
Таким образом, операция сложения над блочными мат рицами одинаковых размеров с одинаковым разбиением на
блоки производится формально так, как если бы вместо бло ков стояли числовые элементы.
Для того чтобы правило умножения матриц можно бы ло перенести на блочные матрицы, необходимо, чтобы все горизонтальные размеры блоков в первом сомножителе сов падали с соответствующими вертикальными размерами бло ков во втором сомножителе. Иными словами, если
Л ц ^12
^21 ^22
Д>1 AS2
и, кроме того, число столбцов блока Ааь равно числу строк блока Вар (а = 1, 2, ..., s; р = 1, 2, ...» и\ б = 1, 2, ..., {), то возможно перемножение матриц А и В формально так, как если бы вместо блоков стояли числовые элементы:
где |
|
АВ = |
С = |
(СаР), |
t |
|
|
|
|
|
|
|
|
|
= |
Лаб5бр |
= |
1, 2, |
. . . , sj Р = 1, 2, . . . , ц). |
6=1
Квадратная матрица, у которой все элементы, располо женные под (над) главной диагональю, равны нулю, назы вается верхней (нижней) треугольной матрицей.
Аналогичные понятия вводятся и для блочных матриц. Блочная матрица
' Аа |
А12 |
А\р |
Л* |
А22 |
Л2р |
А = [ |
Ар2 |
Арр |
\КApi |
называется верхней (нижней) квазитреугольной матрицей,
если все диагональные блоки и сама матрица А — квадрат ные матрицы, а все недиагональные блоки, расположенные под (над) диагональными блоками,— нулевые матрицы.
Блочная матрица А = {Ац) называется квазидиагональной, если все диагональные блоки и сама матрица А — квад ратные матрицы, а недиагональные блоки — нулевые мат рицы. Квазидиагональная матрица является частного вида квазитреугольной матрицей.
Определитель квазитреугольной матрицы А = (Ац)ррсвя зан с определителями диагональных блоков соотношением
р |
(8.1) |
d e M = П etAjjd . |
|
Докажем это. |
матрицу |
Рассмотрим сначала квазитреугольную |
По определению |
|
|
|
d e tA = |
2 |
( - l ) ' * 1’*......atkla2ll, |
|
|
k^^kj (Izjfcj) |
|
|
Так как А12 = 0, то из всех произведений |
яnk„ |
могут быть не равны нулю только те, в которых индексы klt k2, ..., kr принадлежат множеству 1, 2, ..., г. Вследствие
этого остальные индексы |
kr+if kr+2» •••» kn могут |
при |
нимать значения только из |
множества г -f 1 , г + 2 , |
...» п. |
В этих условиях число транспозиций элементов, необходи мых для приведения перестановки 1 , 2 , ..., п к расположе
нию klt k2, |
kn, равно сумме числа транспозиций элемен |
||||
тов, необходимых для |
приведения перестановки 1 , 2 , |
г |
|||
к расположению klt k2, ..., k, |
и числа транспозиций, необ |
||||
ходимых для приведения перестановки r - f |
1 , г -f 2 , |
п |
|||
к расположению 6г+ь |
Ь+2» |
kn: |
|
|
|
t (ftj, k2>• • • |
, Йд) |
it^2*• • • |
i ^r) “h ^2 (^ + 1» |
• • • » |
|
Учитывая это, находим |
|
|
|||
det Л — |
V |
|
| |
kr)“W* |
-I */|j ^ |
|
2 |
|
|
||
ft,. .... |
*/+1' |
—**„=г+1 |
|
|
|
|
*0*1 |
<МЛ |
|
|
|
|
К fltife, |
в/Л^Дг-Н ftr.j_| |
|
|
|
|
2 |
|
(__1 )*' (fc.......Vai*, |
flrftr Й |
|
|
ft,..... ft/==^ |
|
|
|
|
|
ft£<#ft/ |
|
|
|
|
к |
S |
( - 1 )'* ............‘”’a,+i *,+■ |
Onftrt. |
||
|
kr+l' ■■* |
|
|
|
|
|
ft^fcft/ </*Л |
|
|
||
Отсюда |
|
det A — det Au det Л22. |
|
(8.2) |
|
|
|
|
Рассматривая в общем случае матрицу
согласно (8.2) будем иметь det А = det Аи det Л22. Матри-
да А22 снова квазитреугольная. Проделав над ней ту же операцию, получим
det А = det Ап det Л22 det Л83.
После р — 1 таких шагов придем к соотношению (8.1). Таким же путем может быть доказано равенство (8 . 1 ) при менительно в верхней квазитреугольной матрице.
§ 9. Линейные преобразования и матрицы |
|
|||
Пусть т величин у1%у21 .... ут выражаются |
линейно и |
|||
однородно через п других величин хх, х2, ...» х„: |
|
|||
Ух = |
а 11х 1 + aUx i + |
+ |
пХп, |
|
У2 = |
^21^1 4" 022Х24* |
4* ^2пХ>пл |
(9.1) |
|
Ут ~ |
a m lX i + й т 2Х 2 -f- |
4* |
ОтпХп . |
|
Преобразование величин xlt х2, ...» хп в величины у1% Уъ* •••» Ут посредством равенств (9.1) называется линейным преобразованием.
Система равенств (9.1) эквивалентна одному матричному равенству
/У г \ |
/ ап |
|
\ |
( xi \ |
Уг |
а21 |
fl22 |
02я |
X* |
|
|
|
|
9 |
т |
0>т\ |
Ощ2 |
Отп |
|
в чем легко убедиться, выполнив умножение матриц в пра вой части этого равенства и приравняв друг другу соответ ствующие элементы матриц, расположенных слева и справа от знака равенства. Обозначая
Ух |
|
а1Х |
а11 |
а\п |
Уг |
А = |
й21 |
а92 |
о-гп |
|
|
|
9 |
|
1т |
'Т! |
flml |
От2 |
Omnj |
можно вместо (9 . 1) записать коротко: |
|
|
|
|
|
у = Ах. |
|
|
(9.2) |
Таким образом, линейное преобразование (9.1) одно значно определяет матрицу А и, обратно, всякая матрица А определяет некоторое линейное преобразование одной си
стемы чисел в другую.
Допустим, что xlt х2, .... хп в свою очередь выража
ются через величины zlt z2, ...» zp |
посредством равенств |
||||
x i = |
b Xxz x + |
^12г 2 + |
4 - |
b \pzpt |
|
Х2 = |
k 21Zx + |
^2222 + |
+ |
^2pZp, |
(9.3) |
|
|
|
|
|
|
х„ = bn\Zx |
ьп2гг + |
4 *bnpZp. , |
|
Можно |
уъ |
у2, |
...» Ут непосредственно выразить |
через |
zlt zz, ..., |
zp. |
Для |
этого нужно с помощью равенств |
(9.3) |
исключить Ху xZ}...» *п из равенств (9.1). В результате полу
чим |
|
ух = спгг + с{гЧ + |
|
+ |
|
|
|||
|
|
|
P z p > |
|
|||||
|
|
y z = с 21г г + с ггг ъ 4* |
|
+ c2pZp, |
(9.4) |
||||
|
|
|
|
|
|
|
|
||
где |
|
д т = |
C/nlzi 4 *cm2za 4" |
|
4” CtnpZp,. |
|
|||
|
|
|
|
|
|
|
|
||
СИ= |
2 |
aikbki |
(i = |
1 , 2 , |
, |
m; |
/ = 1 , 2 , |
rt* |
|
|
A = |
I |
|
|
|
|
|
(9.5) |
|
В самом деле, учитывая , что (см. (9.1) и (9.3)) |
|||||||||
|
|||||||||
|
|
|
fl |
|
|
P |
|
|
|
|
|
У1 = 2! |
^A |
|
21 |
, |
|
||
|
|
|
A = 1 |
|
|
|
/ = 1 |
|
|
последовательно |
получаем |
|
|
|
|
||||
n |
p |
n |
p |
|
|
|
|
||
= 2 |
= I |
2 |
== 2 |
2 aikbkjZj= |
|
|
|||
A |
/ = 1 |
A=1 j=J |
|
|
|
|
|||
|
|
|
|
P |
|
/ n |
\ |
P |
|
|
|
|
|
= 2 |
( |
2 |
Q/A^A/ )ZI |
— 2 со г/* |
|
|
|
|
|
/=i \ A=I |
/ |
/=| |
откуда и следуют соотношения (9.4), (9.5).
Эту операцию можно выполнить также, используя матПолагая
|
|
( ьп |
Ь1г |
К |
г = |
в = |
( 6а1 |
^22 |
Ьгр |
|
|
|||
|
|
V -„1 |
ЬП2 |
Ьпр |
вместо (9.3) будем |
иметь |
|
|
|
|
х = Bz. |
|
(9.6) |
|
Подставляя (9.6) в (9.2), |
получим |
|
||
|
у = ABz = Cz, |
(9.7) |
где С — матрица с размерами т X р, элементы которой, в соответствии с правилом умножения матриц, определяются формулой (9 5).
Допустим, что квадратная матрица А порядка /г, опре деляющая линейное преобразование
ffi = Q i\x | 4" а & х 2 4" 4” ainXn (i = 1 , 2 , ri), (9.8)
— невырожденная матрица. Тогда система алгебраических
уравнений (9.8) может |
быть разрешена относительно xlt |
х2, .... хп при любых ylt |
у2, ...» у,;, согласно правилу Кра |
мера |
|
Х1 ~
Отсюда
|
а п |
# i /—1 |
9 i |
#1Я -1 |
а \ п |
|
1 |
#21 |
# 2 / —1 |
Уъ |
#2 Я -1 |
#2п |
|
M l |
|
|
|
|
|
|
|
#п! |
#п /—I |
Уп |
#п /+1 |
#пд |
|
X j = |
- щ - 2 |
А ц у ь |
( / = 1 , 2 , |
. . . , п ), |
(9.9) |
где Alf — алгебраическое дополнение элемента ац матри цы А.
В матричной записи (9.9) принимает вид
* = \А А 9 или, если учитывать (6.3),
х = А~ху.
Этот же результат немедленно следует и из матричного равенства
Ах = у
после умножения обеих частей этого равенства слева на мат рицу А ~\
Г л а в а II
ВЕКТОРЫ, ВЕКТОРНЫЕ ПРОСТРАНСТВА, ЛИНЕЙНЫЕ ОПЕРАТОРЫ И МАТРИЦЫ
§ 1. Векторы и векторное пространство
Пусть дано некоторое множество R элементов j , y , z , ...
и числовое поле ZK.
Множество R элементов х , y , z , ... называется линейным пространством, если введены операции сложения элементов и умножения элемента из R на число из d i, т. е.
а) каждым двум элементам х у у £ R поставлен в соот ветствие элемент х + у 6 R, называемый суммой элементов х н у ,
б) каждому элементу х £ R и каждому числу X £ di поставлен в соответствие элементов X х € /?, называемый
произведением числа X на элемент х ,
иэти операции удовлетворяют постулатам:
1)+ у = у + х (коммутативность);
2) {х + у») + z = х + (у + г) (ассоциативность);
3) существует нулевой элемент 0 в R такой, что произве дение числа 0 на любой х £ R равно элементу 0:
0л: = 0;
4)1 х = х\
5)аф л;) = (сф) л:;
6)(а - f Р) х = ах + рлг;
7)а (х -f- у) — ах + &У-
Вдальнейшем элементы х,у),2, ... мы будет называть
векторами, а пространство/? — линейным векторным про• странством или просто векторным пространством.
Пр и м е р ы .
1.В геометрии, физике, механике рассматриваются на правленные отрезки (векторы), для которых операции сложения и умножения на число (действительное) вводятся