ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 183
Скачиваний: 1
Для применения этого способа проверки требуется определить ранг матриц [А у\ и А. Если ранги этих матриц равны между собой, урав нения совместны; в противном случае они несовместны.
Пример. Пусть даны уравнения Ах — у, в которых
2 |
6 |
4 |
2 ~ |
~ 1 |
4 |
15 |
14 |
7 |
и у — 6 |
2 |
9 |
10 |
5J |
- 5 _ |
Тогда, прибегнув к элементарным операциям со строками (см. параграф 7 главы VI), можно определить ранг матрицы А:
2 |
6 |
4 |
2~ |
А ^ 0 |
3 |
6 |
3 |
.0- |
0 |
0 |
0_ |
Следовательно, ранг матрицы А равен 2, |
или г (А) — 2. Расширенная |
||||
матрица будет выглядеть следующим образом: |
|||||
И У1 |
2 |
6 |
4 |
2 |
1 |
4 |
15 |
14 |
7 |
6 |
|
|
2 |
9 10 |
5 |
5 |
С помощью тех же элементарных операций эту матрицу можно приве сти к следующему виду:
2 |
6 |
И у] 9* 0 |
СО |
- 0 |
0 |
4 2 Г
6 3 4
0 0 0.
Таким образом, г [А у] = 2 = г (А), из чего можно заключить: урав нения Ах = у совместны.
В последующих параграфах этой главы будут изложены методы решения уравнений, предполагающие приведение матрицы Л к диаго нальному виду, поэтому мы приведем сейчас еще один способ проверки совместности уравнений, который может быть применен в сочетании с этими методами: допустим, что А представляет собой матрицу ранга г, в ней содержится р строк; через PAQ обозначим диагональную форму А, в таком случае уравнения Ах = у будут совместны тогда
итолько тогда, когда последние р—г элементов Ру равны нулю1.
1По предположению, PAQ представляет собой диагональную форму мат рицы А, тогда РА можно записать следующим образом: РА = J“y4rJ , где Аг — мат
рица, содержащая г линейно-независимых строк, а Р — произведение элемен тарных операторов соответствующего вида. Поэтому, если уравнения А $ = у совместны, то совместны уравнения РАх = Ру, а следовательно, совместной окажется и система j"^4глгj = Ру. Из этого непосредственно вытекает, что послед
ние р — г элементов Ру должны быть равны нулю. И наоборот, допустим, |
что |
последние р —г элементов Ру равны нулю. В таком случае уравнения |
х = |
— Ру совместны, а так как Я - 1 существует, то совместными окажутся и уравне ния Р ~ г j'Hj.j х = у, из этого следует, что Ах = у представляет собой систему
совместных уравнений.
177
С помощью такого подхода мы можем установить совместность тех или иных уравнений. Приведя матрицу А к диагональному виду PAQ, определяем г, т. е. ранг А; он равен числу ненулевых элементов в мат рице PAQ. Затем, зная р н г, мы можем исследовать последние р ~ г элементов Ру: если они равны нулю, то наши уравнения совместны, если, однако, хотя бы один из этих элементов отличен от нуля, уравне ния несовместны.
Пример. Выпишем диагональную форму матрицы PAQ уравнения
(2):
1 |
0 |
0 |
2 |
3 |
1 |
' |
— — |
1 |
0 |
1 |
1 |
1 |
|
2 |
|
|
|
|
|
|
—2 |
1 |
1 |
3 |
5 |
1 |
|
'l |
—— —2~ |
~2 |
0 |
0 |
|
|
2 |
|
] |
|
|
0 |
1 1 — |
0 |
0 |
||
— - |
|||||
|
|
0 |
2 |
0 |
|
о |
О |
0 |
1 |
1 |
Ранг матрицы А равен 2, поэтому для того, чтобы уравнения оказались совместными, последние р — г = 3—2 = 1 элементы Ру должны рав няться нулю. Поскольку же Ру имеет следующий вид:
1 |
0 0 |
14~ |
14 |
|
1 |
1 |
0 |
6 |
— 1 |
Ру-- 2 |
||||
—2 |
1 |
1 |
19 |
—3 |
можно утверждать, что уравнения Ах = у несовместны,
г) выводы
Мы показали, как решать совместные уравнения Ах = у с помощью формул (8) и (9). Требуется выбрать г линейно-независимых строк мат рицы А и представить их в виде [Лх Л 2], где Лх — квадратная невыро жденная матрица г-го порядка. В таком случае будет существовать и Л7 1, следовательно, могут быть использованы выражения (8) и (9). Такой путь решения предполагает, что мы можем выписать соответст вующую матрицу Лх. Для того чтобы отыскать Лх, необходимо, во-пер вых, найти г, т. е. ранг матрицы Л, во-вторых, выделить г линейно-не зависимых строк, а среди элементов этих строк — г линейно-незави симых столбцов, и, наконец, в-третьих, переставить (если это необхо димо) строки и столбцы матрицы Л и элементы векторов х и у таким образом, чтобы придать матрице удобную для расчленения форму.
Нетрудно проделать все это, если система состоит из небольшого числа уравнений; однако в тех случаях, когда число уравнений велико, любая из перечисленных операций может оказаться чрезвычайно тру доемкой. К счастью, существует метод, который позволяет избежать всех этих затруднений, он не требует ни вычисления ранга Л, ни
178
отыскания невырожденной подматрицы. Кроме того, этот метод может применяться и тогда, когда число уравнений оказывается больше или меньше числа неизвестных (матрица Л прямоугольна), в этом случае он будет лишь незначительно отличаться от способа, применяемого при равенстве числа уравнений числу неизвестных (матрица А квадратна). Такой метод предполагает использование матриц, называемых обоб щенными обратными матрицами.
3. ОБОБЩЕННЫЕ ОБРАТНЫЕ МАТРИЦЫ
а) ОПРЕДЕЛЕНИЯ
Пенроуз [6] показал, что для любой матрицы А существует некая единственная матрица К, удовлетворяющая следующим четырем условиям:
1) АКА = А; |
2) КАК = К\ |
3) (КА)' - /04; |
4) (АК)' = АК- |
Такую матрицу К Пенроуз назвал обобщенной обратной к А матрицей; он показал, что обобщенная обратная матрица всегда существует и однозначно определена при любой форме исходной матрицы А, кото рая может быть квадратной (вырожденной или невырожденной) или прямоугольной.
Пример. Пусть дана матрица А:
Г 1 О 2 П
1 2 О
Тогда матрица К, удовлетворяющая условиям 1—4, имеет следующий вид:
|
|
6 |
—2 |
—6 |
10' |
К |
66 |
0 |
—11 |
0 |
22 |
|
12 |
7 |
—12 —2 |
||
|
|
Пенроуз не только доказал существование и единственность матрицы, удовлетворяющей условиям 1—4, он показал также, как можно ею пользоваться при решении линейных уравнений1. Однако, если порядок матрицы А достаточно высок, построение обобщенной обратной матри цы К требует громоздких вычислений; между тем существуют различ ные модификации К, которые оказываются столь же полезными при
1 Этот вопрос рассматривался также в работах Мура [5] и Бьеммэра [1], об зор исследований в данной области содержится в работе Мэлика [4].
179
решении уравнений, но вычисляются более простым способом. Такие модификации обобщенной обратной матрицы могут носить различные названия. Рао [7] и Грэвилл [3], например, употребляют термин «мнимая обратная матрица», Уилкинсон [9] применяемую матрицу называет «эффективной обратной», а Голдмэн и Зилэн [2] именуют любую матри цу, удовлетворяющую условиям 1 , 2 и 3, «слабо обобщенной обратной матрицей». Термином «обобщенная обратная матрица» Рао [8] называет также любую матрицу Т, обеспечивающую решение совместных урав нений Ах = у в виде х = Ту. Далее, при доказательстве теоремы 1 будет показано, что, если такое решение существует, то имеет место АТА — А, т. е. тем самым предполагается условие 1. По этой причине, а также потому, что большая часть предлагавшихся в литературе обоб щенных обратных матриц Пенроуза удовлетворяет условию 1, мы при няли в этой книге определение «обобщенная обратная матрица» для обозначения любой матрицы G, которая удовлетворяет условию AGA = =А, т. е. для обозначения любой матрицы, удовлетворяющей первому из четырех условий Пенроуза*.
Из теоремы, которую мы упоминали, следует, что условия, пред полагаемые определением Рао [8], включены в наше определение обоб щенной обратной матрицы; то же самое можно сказать и о большинст ве других приводившихся определений, поскольку все они опираются на условие 1. Остальные условия (2, 3, 4), которые могут содержаться в том или ином определении, не противоречат нашему определению, а лишь дополнительно ограничивают его; самое сильное ограничение налагает одновременное включение всех четырех приведенных ус ловий. В этом случае мы можем воспользоваться тем определением обобщенной обратной матрицы, которое можно найти в работе Пен роуза, — «единственная обобщенная обратная матрица»1.
б) ПРОИЗВЕДЕНИЕ |
Я = GA |
|
|
|
|
|
|
||
По определению, G представляет собой любую матрицу, удовлетво |
|||||||||
ряющую |
условию |
AGA = А. Отсюда |
можно вывести |
два |
простых |
||||
свойства произведения Н = GA. Во-первых, ранг |
матрицы |
Я |
равен |
||||||
рангу |
А; |
ведь, |
поскольку Я = |
GA, |
г (Я) < |
г (А), |
а |
так |
как |
А = |
AGA |
= АН, |
г (А) <; г (Н). |
Следовательно, |
г (Н) = |
г (Л). |
|||
Во-вторых, Я2 = |
Я, потому что Я2 = GAGA = G (AGA) = |
GA = Я . |
Матрицы, обладающие таким свойством, часто называют идемпотентными: так, матрица Я идемпотентна, если Я2 = Я.
*В русском переводе книги Рао С. Р. «Линейные статистические методы и их применения» (М., «Наука», 1968) также принято сокращенное обозначение обоб щенной обратной матрицы — g-обратная матрица. — Прим, перев.
ХВ параграфе 9 главы V рассматривались левая и правая обратные матри цы; каждая из таких матриц удовлетворяет трем из четырех условий. Так, если
L — левая обратная матрица |
по |
отношению к Л и LA = /, тогда |
К = L |
и выполняются условия 1, 2 и 3. Аналогичным образом правая обратная |
матрица |
||
будет удовлетворять условиям |
1, 2 |
и 4. |
|
180