
От Пантелеева / Книги студентам / лИТЕРАТУРА / Тартышников
.pdfЕ. Е. Тыртышников |
37 |
|
|
6.4Группа обратимых матриц
Множество обратимых n ×n-матриц относительно операции умножения образует груп-
пу. Для доказательства есть все, кроме факта обратимости произведения обратимых матриц. Но это проверяется непосредственно: если A è B обратимы, то
(B−1A−1)(AB) = B−1(A−1A)B = B−1 I B = I,
(AB)(B−1A−1) = A(BB−1)A−1 = A−1 I A = I.
Отсюда
(AB)−1 = B−1 A−1.
Задача. Пусть H группа невырожденных верхних треугольных матриц порядка n, G группа всех обратимых матриц порядка n. Верно ли, что ABA−1 H для любых матриц B H è A G (в этом случае подгруппа H называется нормальным делителем группы G)?
6.5Обращение невырожденной матрицы
Квадратная матрица с отличным от нуля определителем называется невырожденной. Пусть Ae = [Aij] матрица порядка n, в которой элемент Aij есть алгебраическое
дополнение к элементу aij в матрице A. Матрица Ae> называется присоединенной äëÿ матрицы A.
Теорема. Если матрица A порядка n невырожденная, то она обратима и при этом
|
A−1 = |
1 |
Ae |
>, |
ãäå |
det A |
|||
|
A> присоединенная матрица для A. |
|
|
|
|
e |
|
|
|
Доказательство. Заметим, что
n |
|
0, |
k = i. |
( ) |
j=1 aijAkj = |
||||
X |
|
det A, |
k = i, |
|
|
|
6 |
|
Ïðè k = i равенство ( ) получается применением теоремы Лапласа при разложении определителя матрицы A ïî k-й строке. При k 6= i левая часть равенства ( ) представляет собой разложение по k-й строке определителя матрицы, полученной из A заменой k-й строки на i-ю. Такой определитель равен нулю как определитель матрицы с двумя одинаковыми строками. Например, если n = 3, i = 1, k = 2, òî
a11A21 + a12A22 + a13A23 = a11 |
− |
a32 |
a33 |
|
+ a12 |
a31 |
a33 |
|
+ a13 |
− |
a31 |
a32 |
|
|||||||||
|
|
|
|
|
|
a12 |
a13 |
|
|
|
a11 |
a13 |
|
|
|
|
a11 |
a12 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a11 |
a12 |
a13 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
= |
|
a |
|
a |
|
|
a |
|
= 0. |
|
|
|
|
|
|
|
|
|
|
|||
a11 |
a12 |
a13 |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
31 |
|
32 |
|
33 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Далее, |
n |
n |
|
|
aij(A>)jk, 1 ≤ i, k ≤ n. |
||
|
aijAkj = |
j=1 |
|
|
=1 |
e |
|
|
Xj |
X |

38 Лекция 6
Поэтому, в силу соотношений ( ),
AA > = |
det A ... |
det A |
|
= det A |
· |
I. |
e |
|
|
|
|
||
|
|
|
|
|
|
Используя теорему Лапласа для разложения определителя по k-му столбцу, находим
n
X
i=1
aijAik = |
0, |
k 6= j. |
|
|
det A, |
k = j, |
|
A >A = |
det A ... |
det A |
|
= det A |
· |
I. 2 |
e |
|
|
|
|
||
|
|
|
|
|
|
Задача. Докажите, что любую невырожденную матрицу можно сделать вырожденной, изменив лишь один из ее элементов.
6.6Правило Крамера
Теорема. Пусть A невырожденная матрица порядка n. Тогда система линейных алгебраических уравнений Ax = b имеет и притом единственное решение x c компонентами
xi = |
det Ai |
, 1 |
≤ i ≤ n, |
det A |
ãäå Ai матрица, получаемая из A заменой i-го столбца на b.
Доказательство. Согласно теореме об обращении невырожденной матрицы,
|
1 |
|
|
1 |
n |
|
|
x = A−1b = |
e |
>b |
Xj |
|
|||
|
|
|
|
||||
|
xi = det A |
≤ i ≤ n. |
|||||
det A A |
Ajibj, 1 |
||||||
|
|
|
|
|
|
=1 |
|
n
Остается заметить, что сумма P Ajibj есть разложение по i-му столбцу определителя
j=1
матрицы Ai. 2
6.7Определитель произведения матриц
Теорема. Определитель произведения квадратных матриц равен произведению их определителей.
Доказательство. Пусть A = [a1, . . . , an] матрица порядка n со столбцами a1, . . . , an è B = [bij]. Тогда любой столбец матрицы AB есть линейная комбинация столбцов матрицы A с коэффициентами из соответствующего столбца матрицы B:
" n |
n |
# |
X |
X |
|
AB = |
bi1 1 ai1 , . . . , |
bin n ain . |
i1=1 |
in=1 |
|
Е. Е. Тыртышников |
39 |
|
|
Используя линейность определителя по каждому столбцу, получаем
nn
XX
det(AB) = . . . bi11 . . . binn det[ai1 , . . . , ain ]
|
i1=1 in=1 |
= |
X bσ(1)1 . . . bσ(n)n det[aσ(1), . . . , aσ(n)] |
|
σ Sn |
= |
bσ(1)1 . . . bσ(n)n sgn(σ)! det A = det B · det A. 2 |
|
X |
|
σ Sn |
Задача. Пусть P обратимая матрица порядка n. Докажите, что для любых столбцов u, v Rn×1 выполняется равенство det(I + uv>) = det(I + (P u)(P −>v)>).
Задача. Докажите, что определитель трехдиагональной матрицы не изменится, если каждый наддиагональный элемент умножить, а каждый поддиагональный элемент поделить на одно и то же число.
6.8Обратимость и невырожденность
Теорема. Квадратная матрица обратима тогда и только тогда, когда она невырож-
денная.
Доказательство. Пусть для матрицы A существует обратная матрица A−1. Тогда AA−1 = I и, в силу теоремы об определителе произведения матриц,
det A · det A−1 = det I = 1 det A 6= 0.
Åñëè det A 6= 0, òî A обратима по теореме об обращении невырожденной матрицы. 2
Следствие. Столбцы матрицы A линейно независимы тогда и только тогда, когда det A 6= 0.
Задача. Пусть In è Im единичные матрицы порядка n è m. Докажите, что для любых матриц A размеров m×n è B размеров n×m из обратимости Im −AB вытекает обратимость In −BA. Докажите также, что обратимость каждой из этих матриц равносильна обратимости матрицы порядка m + n с блочным разбиением вида
Im |
A . |
B |
In |
Задача. Даны числа a è b такие, что 1 − ab 6= 0. Докажите, что матрица
|
1 |
a |
a2 |
... |
an |
|
|
b |
1 |
a |
... an−1 |
||
A = |
b2 |
b |
1 |
... an−2 |
|
|
|
... ... |
... |
... ... |
|||
|
|
|
|
|
|
|
|
bn |
bn−1 |
bn−2 |
... |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
обратима и обратная к ней матрица является трехдиагональной.
Задача. Пусть в условии предыдущей задачи a è b квадратные блоки порядка n, à 1 заменяется
единичной матрицей I того же порядка. Докажите, что если блок I −ab является обратимой матрицей,
то блочная матрица A с блоками порядка n обратима и при этом обратная к ней матрица является блочно-трехдиагональной.
40 |
Лекция 6 |
|
|

Лекция 7
7.1Разделение переменных и матрицы
При изучении функций от двух переменных особую роль играют функции с разделенными переменными f(x, y) = u(x)v(y) или суммы таких функций 1
f(x, y) = u1(x)v1(x) + . . . + ur(x)vr(y).
Пусть дана m × n-матрица A. Ее элемент aij
дискретных переменных i {1, . . . , m}, j {1, . . . , n}. В данном случае разделение переменных означает, что
aij = uivj, |
1 ≤ i ≤ m, 1 ≤ j ≤ n. |
|
Отсюда легко вывести, что A есть произведение столбца и строки: 2 |
||
|
u1 |
v1 |
A = uv>, |
u = u. .m. , v = |
.v.n. . |
7.2 Скелетное разложение |
|
|
Теперь предположим, что aij = ui1vj1 + . . . + uirvjr, |
1 ≤ i ≤ m, 1 ≤ j ≤ n. Â ýòîì |
|
случае A является суммой r матриц вида столбец на строку |
||
r |
uk = u.mk. . , vk = v.nk. . . |
|
A = k=1 ukvk>, |
||
X |
u1k |
v1k |
|
|
Это же равенство, записанное в виде произведения двух матриц
|
|
|
|
u11 |
. . . |
u1r |
v11 |
A = UV > = . . . . . . . . . . . . |
|||
um1 |
. . . |
umr |
v1r |
. . . vn1
. . . . . .
. . . unr
, U = [u1, . . . , ur], V = [v1, . . . , vr], ( )
называется скелетным разложением матрицы A. Оно означает, что каждый столбец матрицы A есть линейная комбинация столбцов матрицы U, а каждая строка матрицы A есть линейная комбинация строк матрицы V >. Отсюда сразу же вытекает
Теорема. Размерность линейной оболочки, натянутой на столбцы матрицы A, совпадает с размерностью линейной оболочки, натянутой на ее строки:
|
|
|
|
a1> |
|
dim L(a1 |
, . . . , an) = dim L(a1 |
, . . . , am), |
bm> |
||
A = [a1, . . . , an] = . . . . |
|||||
|
b |
b |
|
b |
|
|
|
|
|
a |
|
1Разделение переменных с большим успехом применяется для приближения функций общего вида. 2Матрица uv> иногда называется внешним произведением векторов u Cm è v Cn.
41
42 |
Лекция 7 |
|
|
Доказательство. |
Пусть столбцы матрицы U образуют базис линейной оболочки |
L(a1, . . . , an), à j-é |
столбец матрицы V > состоит из коэффициентов их линейной |
комбинации, дающей столбец aj. Тогда, очевидно, A = UV >. В силу предваряющего теорему замечания, размерность линейной оболочки строк матрицы A не выше числа V >, которое равно, по построению, размерности линейной оболочки столбцов матрицы A. Противоположное неравенство доказывается аналогично роль
столбцов и строк меняется транспонированием. 2
7.3Ранг матрицы
Размерность линейной оболочки столбцов (строк) матрицы иногда называется ее столб-
цовым (строчным) рангом. Поскольку столбцовый и строчный ранги совпадают, их об-
щее значение было бы естественно называть просто рангом матрицы. Из проведенного
нами доказательства этого факта вытекает также то, что это значение равно наименьшему числу матриц вида столбец на строку , дающих в сумме данную матрицу.
Однако, обычно дается другое определение: рангом матрицы называется наивыс-
ший порядок ее отличных от нуля миноров. Соответствующие минор и подматрица
называются базисным минором è базисной подматрицей. В силу уже установленной
эквивалентности обратимости и невырожденности, ранг матрицы равен наивысшему порядку обратимых подматриц в данной матрице. Обозначение: rankA.
Два очевидных свойства ранга матрицы
rankA ≤ min(m, n), rankA = rankA>.
Менее очевидно, что наивысший порядок отличных от нуля миноров матрицы совпадает с ее столбцовым и строчным рангом. Давайте это докажем.
7.4Окаймление обратимой подматрицы
Начнем с полезного вспомогательного предложения. Пусть матрица Q порядка k + 1 имеет блочный вид
Q = |
P v |
, P Rk×k, u, v Rk×1. |
u> c |
В этом случае Q называется окаймлением подматрицы P .
Если подматрица P обратима, а ее окаймление Q является необратимой матрицей, то последний столбец матрицы Q есть линейная комбинация первых k столбцов.
Доказательство. Используя правило умножения блочных матриц (см. Лекцию 1), легко проверить справедливость равенства
−u>P −1 |
1 |
u> |
c |
= |
0 |
γ |
, γ = c − u>P −1v. |
I |
0 |
P |
v |
|
P |
v |
|
Обозначим матрицу в правой части через M. Как произведение обратимой и необратимой матриц, M не может быть обратимой матрицей. Но она имеет блочно-треугольный вид, и если бы блоки P и γ были оба обратимы, то M имела бы обратную матрицу вида
M−1 = |
P −1 |
− |
P −1vγ |
− |
1 |
. |
0 |
γ−1 |
|
Е. Е. Тыртышников |
43 |
|
|
(Равенство MM−1 = I проверяется непосредственно.) Поскольку M не является обратимой матрицей, непременно
γ = c − u>P −1v = 0 c = u>P −1v.
Следовательно, |
u> |
(P −1v) = |
c . |
2 |
|
||||
|
P |
|
v |
|
7.5Теорема о базисном миноре
Теорема. Столбцы (строки), содержащие базисный минор, являются линейно независимыми, при этом любой столбец (любая строка) данной матрицы является их
линейной комбинацией.
Доказательство. Не ограничивая общности, предположим, что базисная подматрица P порядка k расположена в левом верхнем углу матрицы A размеров m × n. Таким
образом,
|
|
P |
|
vk+1 |
A = |
|
u> |
a |
k+1 k+1 |
k+1 |
|
|||
|
|
|
|
|
|
|
u. .m>. |
am. .k.+1 |
. . . |
vn |
. . . |
ak+1 n |
. . . . . . |
|
. . . |
amn |
, uk+1, . . . , um, vk+1, . . . , vn Rk×1.
По условию теоремы, любая подматрица порядка k + 1 âèäà
M = |
P |
vj |
, i, j > k, |
|
ui> |
aij |
|||
|
|
является необратимой. По лемме о необратимом окаймлении обратимой подматрицы, последний столбец в M есть линейная комбинация первых k столбцов. При этом ко-
эффициенты данной линейной комбинации не зависят от i (поскольку определяются вектором P −1 vj). Значит, j-й столбец матрицы A при j > k есть линейная комбинация первых k столбцов. Линейная независимость первых k столбцов доказывается следующим образом: пусть их линейная комбинация с коэффициентами α1, . . . , αk равна 0,
тогда |
P |
|
.α.1. |
|
= 0 |
|
|
.α.1. |
|
= 0. |
|
|
|
αk |
|
|
|
αk |
|
|
Утверждение теоремы относительно строк доказывается переходом к транспонированной матрице. 2
Матрица имеет ранг r тогда и только тогда, когда некоторый минор порядка r отличен от нуля, а все окаймляющие его миноры порядка r + 1 равны нулю.
Следствие 2.
Замечание. Теорема о базисном миноре не использует доказанной ранее теоремы о равенстве столбцового и строчного рангов. По существу, она дает еще одно доказательство этой теоремы.
Задача. Путь A n Ч n-матрица ранга k, а B любая невырожденная подматрица порядка k. Обозначим через R подматрицу размеров k Ч n, состоящую из строк матрицы A, содержащих подматрицу B, а через C подматрицу размеров n Ч k, состоящую из столбцов, содержащих B. Доказать,
÷òî
A = CB−1R.

44 |
Лекция 7 |
|
|
|
Задача. Докажите, что подматрица, расположенная на пересечении r линейно независимых строк |
èr линейно независимых столбцов матрицы ранга r, является невырожденной.
Задача. Известно, что A> = −A. Докажите, что ранг матрицы A число четное.
7.6Ранги и матричные операции
Утверждение 1. Ранг суммы матриц не превосходит суммы их рангов:
rank (A + B) ≤ rank A + rank B.
Доказательство. Очевидно, A è B должны иметь одинаковое число столбцов:
A = [a1, . . . , an], B = [b1, . . . , bn].
ßñíî, ÷òî L(a1 + b1, . . . , an + bn) L(a1, . . . , an, b1, . . . , bn). В меньшей линейной
оболочке выберем какую-либо систему векторов, образующую базис. Согласно лемме о дополнении до базиса, базис в большей линейной оболочке можно получить путем дополнения данной системы какими-то векторами из большей линейной оболочки. Поэтому
rank (A + B) = dim L(a1 + b1, . . . , an + bn) ≤ dim L(a1, . . . , an, b1, . . . , bn).
Пусть p = rank A, q = rank B, и предположим, не ограничивая общности, что базис в L(a1, . . . , an) образуют первые p векторов, а базис в L(b1, . . . , bn) первые q векторов.
Тогда L(a1, . . . , an, b1, . . . , bn) = L(a1, . . . , ap, |
b1, . . . , bq) |
|
dim L(a1, . . . , an, b1, . . . , bn) |
≤ p + q. |
2 |
Утверждение 2. Ранг произведения матриц не превосходит ранга каждого из сомножителей:
rank (AB) ≤ min (rank A, rank B).
Доказательство. Достаточно заметить, что каждый из столбцов матрицы AB является линейной комбинацией столбцов матрицы A. Поэтому линейная оболочка столбцов матрицы AB содержится в линейной оболочке столбцов матрицы A. Следовательно,
rank (AB) ≤ rank A. Далее,
rank (AB) = rank (AB)> = rank (B>A>) ≤ rank B> = rank B. 2
Утверждение 3. Ранг матрицы не изменяется при умножении ее слева или справа на обратимую матрицу.
Доказательство. Пусть B = P AQ, ãäå P è Q обратимые матрицы. 3 Â ñèëó
3Конечно, порядок P равен числу строк, а порядок Q числу столбцов матрицы A.
Е. Е. Тыртышников |
45 |
|
|
предыдущего утверждения, rank B ≤ rank A. В то же время, A = P −1BQ−1
rank A ≤ rank B. 2
Это утверждение полезно при вычислении ранга. Обычно это делается с помощью элементарных преобразований строк и столбцов, упрощающих вид матрицы путем исключения ее элементов; эти преобразования сводятся к умножению матрицы слева и справа на некоторые обратимые матрицы (специального вида) и поэтому сохраняют ранг.
Утверждение 4. При изменении k строк (столбцов) матрицы ее ранг не может измениться больше чем на k.
Доказательство. Достаточно рассмотреть случай k = 1. Пусть ранг равен r, а столб-
öû aj1 , . . . , ajr являются базисными. Если столбец aj совпадает с одним из этих столбцов, то при любом его изменении данная система столбцов будет иметь линейно неза-
висимую подсистему с числом столбцов r − 1. Поэтому ранг новой матрицы не меньше
÷åì r − 1 и, конечно, не больше r. Если столбец aj не сопадает ни с одним из выбран- ных базисных столбцов, то любой столбец новой матрицы линейно выражается через столбцы aj1 , . . . , ajr , aj, поэтому ее ранг не превышает r + 1 и не меньше чем r, òàê êàê
имеется система из r линейно независимых столбцов aj1 , . . . , ajr . 2
Следствие. При добавлении (изъятии) k столбцов ранг матрицы не может измениться более чем на k.
Задача. Пусть A è B матрицы ранга 1. Докажите, что если AB = BA 6= 0, то ранг матрицы
A + B не больше 1.
Задача. Заданы столбцы x, y Rn, причем x 6= 0. Докажите, что существует симметричная матрица A Rn×n такая, что rankA ≤ 2 è Ax = y.
Задача. Матрица A имеет r столбцов, а матрица B имеет r строк. Докажите, что
r≥ rank(A) + rank(B) − rank(AB).
7.7Однородная система линейных алгебраических уравнений
Система линейных алгебраических уравнений с нулевой правой частью
Ax = 0 ( )
называется однородной. Пусть в данной системе имеется m уравнений и n неизвестных. Тогда матрица коэффициентов A имеет размеры m × n. Рассмотрим A как систему столбцов A = [a1, . . . , an] и предположим, что ее ранг равен r. Не ограничивая общности, будем считать, что базисная подматрица в A расположена на первых r столбцах
будем называть их базисными. Отвечающие базисным столбцам компоненты решения x1, . . . , xr будем также называть базисными, а оставшиеся компоненты xr+1, . . . , xn свободными. Таким образом, вектор-решение имеет вид
x =
x1
. . .
xr
xr+1
.
. . .
xn

46 Лекция 7
Система ( ) равносильна равенству
x1a1 + . . . + xrar = −xr+1ar+1 − . . . − xnan. ( )
По теореме о базисной подматрице, столбцы ar+1, . . . , an принадлежат линейной обо- лочке столбцов a1, . . . , ar. Поэтому при любом выборе значений свободных неизвестных значения базисных неизвестных, удовлетворяющих равенству ( ), существуют и определяются однозначно. Таким образом, существуют n − r векторов вида
|
x.r. |
1. |
|
|
x.r. |
2. |
|
|
|
xr.n. |
−. |
r |
|
|
|||
|
|
x11 |
|
|
|
x12 |
|
|
|
|
x1 n−r |
|
|
||||
1 |
0 |
|
2 |
1 |
|
n−r |
|
0 |
|
|
|
||||||
v = |
. . . |
|
, v = |
. . . |
|
, . . . v |
= |
. . . |
|
|
, |
||||||
|
1 |
|
|
|
0 |
|
|
|
0 |
|
|
|
|||||
|
0 |
|
|
0 |
|
|
|
1 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
каждый из который является решением системы ( ):
Av1 = 0, . . . , Avn−r = 0.
Векторы v1, . . . , vn−r линейно независимы:
α1v1 + . . . + αn−rvn−r =
.. .
α1 α2
. . .
αn−r
= 0 α1 = . . . = αn−r = 0.
Кроме того, любая линейная комбинация этих векторов является решением системы ( ) и, более того, если x есть произвольное решение системы ( ), то
x = xr+1v1 + . . . + xnvn−r.
Таким образом, мы доказали следующее важное утверждение.
Теорема. Множество решений однородной системы линейных алгебраических уравнений ( ) совпадает с линейной оболочкой L(v1, . . . , vn−r) линейно независимых векторов
v1, . . . , vn−r.
Следствие. dim L(v1, . . . , vn−r) = n − r.
Линейно независимую систему решений w1, . . . , wk системы Ax = 0 называют фундаментальной системой, если ее линейная оболочка L(w1, . . . , wk) совпадает с множеством всех решений однородной системы Ax = 0.
Следствие. Число векторов в любой фундаментальной системе решений для Ax = 0 равно n − r, где r = rank A.
Для доказательства достаточно заметить, что линейная оболочка векторов фундаментальной системы решений имеет базис из построенных выше векторов v1, . . . , vn−r.
Задача. Даны матрицы A è B порядка n такие, что AB = 0 и при этом матрица A + B невырожденная. Доказать, что rankA + rankB = n.