
книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов
.pdfНа основании свойства 8 имеем
п |
если |
i Ф / |
||
|
||||
*-•=1 |
если |
i — j |
||
|
|
|
||
Это означает, что |
|
|
|
|
CU |
". если |
i ф j, |
||
если |
i — j. |
|||
|
Поэтому A - B=E„, т. e.
В = A—1.
Таким образом, доказано, что любая неособенная квадратная матрица имеет обратную матрицу.
П р и м е р . Найти матрицу, обратную к матрице
3 4 1
А= 0 3 4 2 4 5
Р е ш е н и е . Вычислим опр. А и Aij:
|
|
|
3 4 1- |
= 23, |
|
|
опр. А — 0 3 4 |
||||
|
|
|
2 4 5 |
|
|
Ац |
|
|
1, |
А\ 2 |
-- (-- |
|
/li3 = |
( - U I+ 3- |
0 3 |
|
|
|
2 4 | — |
||||
А2 1 |
= (-D * + 1- |J |
J| = “ |
16’ |
А2 2 — (— |
|
|
/123 = (_1)з|-г. |
3 4 |
! = - |
||
|
2 4 |
||||
Ли = |
(-1 )» + ‘-|5 J|=13 . |
-4и = (-1 )*+ *-| 3 I |
^зз = (—1)3+ 3,1 q з | — 9.
=8 ,
=13,
— 12 и
69
Составим матрицы А и А '- |
|
|
|
|
||
—1 |
8 |
—6 |
1 — |
—1 —16 |
13 |
|
—16 |
13 |
—4 |
и А' = |
8 |
13 |
—12 |
13 |
—12 |
9 |
| |
—6 |
—4 |
9 |
поэтому искомая матрица имеет вид |
|
|
|
|
|
1 |
16 |
13 |
|
_ |
23 |
23 |
|
23 |
|
8 |
13 |
12 |
|
|
23 |
23 |
~ |
23 |
|
6 |
4 |
|
9 |
~ |
23 |
23 |
|
23 |
§3. Запись и решение системы линейных уравнений
вматричной форме
Рассмотрим систему |
« линейных уравнений с п неизвестными: |
||
( filial + |
«12*2 + |
■■• + |
а\пхп — bi, |
.......................................................... |
|
|
(5) |
( апХх\ + |
ап2 х 2 + |
• ■• + |
0-ППХП= Ьп. |
Обозначим через А матрицу |
коэффициентов системы (5) |
||
|
йц- • -а1п |
|
|
|
ап1 '" ’апп |
|
|
через X — матрицу-столбец неизвестных |
н |
Ь\
через В — матрицу-столбец свободных членов
Ьп
70
Вычислим произведение
Й11 • • |
Х 1 |
ЯцА'х + |
Ях2*2 • • • + |
a l n x,n |
|
d n v ' ’ а п п |
Xfi |
a n i x i + |
о-пгХг 4- • |
■+ |
а п п х п |
Это произведение является матрицей-столбцом, |
содержащей п |
||||
элементов. |
|
|
|
|
|
Рассмотрим матричное уравнение |
|
|
|
||
|
|
А - Х = В . |
|
|
(6) |
Это уравнение может быть записано так:
<2llx l |
+ |
Й12*2 "Ь" |
• + |
а \ пхп |
ь 1 |
ап1 х 1 |
+ |
<^П2Х 2 +• |
’ + |
&ппХп |
b,1 |
|
|
|
|
|
Написанное равенство двух матриц-столбцов имеет место в том случае, если каждый элемент одной из них равен соответствующему элементу другой, т. е. если справедлива система (5).
Это означает, что систему (5) из п линейных уравнений с п не известным можно заменить одним матричным уравнением (6). Ре шим это уравнение.
Если матрица А является неособенной, то можно найти матри цу А _1, обратную к матрице А.
Умножим левую и правую части уравнения (6) слева на Л-1, по лучаем
Л-1-(Л-Х) = A -i-B
или
(A -i'A )-X = Л-1-В,
но (Л-1 • А )—Еп, а Е п -Х ^Х , поэтому из последнего равенства полу чаем решение уравнения (6) в матричной форме:
X = Л-1-В.
71
Чтобы получить решение исходной системы, нужно вычислить произведение А ~ 1 • В и определить элементы вектора-столбца X. Они образуют искомое решение.
П р и м е р . Решить систему
2xi -j- Х2 |
— 4ха — 5, |
|
|
—xi |
+ хз — О, |
|
|
х 2 |
Зх3 — 2. |
|
|
|
2 1 —4 |
||
Р е ш е н и е . В этом примере |
Л —1 0 |
1 , х = |
|
5 |
0 |
1 |
3 |
|
|
|
и В — О
2
Вычислим А~'. Для этого находим опр. А и Л,-j:
|
опр. А = |
2 1 |
—4 |
|
|
|
||
|
-1 |
0 |
1 |
|
|
|
||
|
|
|
О |
1 |
3 |
|
|
|
— 1, |
А12 = 3, |
Ац — — 1, |
A2i |
— — 7, |
А2 2 = 6, |
|||
Лцз ~ |
— 2 , Лз1 = |
1, |
|
Л32 = |
2 и Л33 = |
1. |
||
Составляем Л и Л': |
|
|
|
|
|
|
|
|
|
—1 3 —1 |
и Л' = |
- 1 |
- 7 1 |
|
|||
|
—7 6 —2 |
3 |
6 2 |
|
||||
|
1 2 |
1 |
|
|
|
—1 |
- 2 1 |
|
Поэтому |
|
|
|
|
|
|
|
|
|
|
_1_ |
|
7_ |
_1_ |
|
||
|
|
|
5 |
|
5 |
5 |
|
|
|
|
|
|
|
6_ |
_ 2 |
|
|
|
|
|
|
|
5 |
5 |
|
|
|
|
J_ |
|
2___ 1^ |
|
|||
|
|
|
5 |
|
5 |
5 |
|
|
Х 1
х 2
Хз
72
Окончательно находим
или
Таким образом, решение системы линейных уравнений (5) сво дится к вычислению обратной матрицы А ~1 и произведения А~ 1 ■В.
Этим методом удобно пользоваться, когда |
приходится |
решать одну |
|||
и ту же систему с разными свободными членами. |
|
|
|||
П р и м е р . |
Рассмотрим системы |
|
|
|
|
2 * i + |
* 2 — 4 * з — 2 , |
2*1 -|- *2 — 4*3 —■8, |
|||
— *1 |
+ *3 = — 1, и |
— * i |
-|- *з = 4, |
||
|
*2 “Ь 3*з ~ 4 |
|
*2 |
3*з = |
— 7. |
|
|
2 |
|
8 |
через Вi и Вг\ |
членами. Обозначим матрицы-столбцы —1 |
И |
4 |
|||
|
|
4 |
|
—7 |
|
имеем:
|
|
|
1 |
|
7 |
1 |
|
|
|
|
_ |
5 |
~ |
5 |
5 |
|
2 |
X i |
|
|
3 |
|
6 |
2 |
|
|
х2 = |
A - ^ B i - - |
|
|
—1 |
||||
5 |
|
5 |
~~5 |
|
||||
*3 |
|
|
|
|
4 |
|||
|
|
|
1 |
|
2 |
1 |
|
|
|
|
~ |
5 |
~ |
5 |
5 |
|
|
|
|
|
|
|
1 |
|
|
9 |
|
|
|
|
|
— -4 |
|
5 |
|
|
|
|
|
|
5 |
|
|
|
|
|
( - D + T - 4 |
— |
|
8 |
|||
|
|
|
5 |
|||||
|
|
2 \ . . . . |
I |
|
|
4 |
||
|
|
тг |
8 |
|
5 ' |
4 |
|
5 |
поэтому |
9_ |
|
|
|
|
|
||
Xi |
~5 |
лг2 = — |
и *3 = — |
|
|
|||
|
|
|
|
|
|
|
|
|
Для второй системы находим |
1 |
|
7 |
1 |
|
|
||
|
|
|
|
|
|
|||
|
|
~ |
5 |
_ |
5 |
5 |
|
8 |
I *1 |
|
|
3 |
|
6 |
2 |
|
|
Х 2 | = А - 'В г = |
|
5 |
1Г |
5 |
|
4 |
||
I *3 |
|
|
—7 |
|||||
|
|
|
1 |
|
2 |
1 |
|
|
, 1\ |
/ |
7\ |
5 |
1 |
5 |
5 |
|
|
|
|
|
|
43 |
||||
Нт)-8+(-т)-4+т'(-7> |
~ |
5 |
||||||
3 |
6 |
2 |
|
|
|
|
|
34 |
т-8+т-4+т-(-7) |
|
= |
|
5 |
||||
|
|
|
|
|
|
|
|
|
-т)'8+(- т)-4+т-<-7) |
|
|
23 |
|||||
|
~ |
5 |
74
Поэтому |
43 |
34 |
23 |
|
|
|
|||
*i ~ |
5 |
х3 = — и х3 = |
* |
|
|
5 |
5 |
||
Определив один раз |
матрицу 7М, |
достаточно |
при изменении |
свободных членов вычислять произведение этой матрицы на матрицу-
столбец свободных членов. |
|
|
§ 4. Ранг |
матрицы |
|
Пусть дана матрица |
|
|
Ч ц - |
■ а ы |
|
а 2 1 ‘ |
• й ъ п |
|
Я m |
l ' |
' а т п |
типа тХп. Из элементов этой матрицы можно составлять определи
тели различных порядков. |
Возьмем произвольным образом k строк |
и k столбцов, 6^m in (т, |
п). |
Элементы, стоящие на пересечении этих строк и столбцов, обра зуют квадратную матрицу. Определитель этой матрицы называется минором k-ro порядка матрицы А.
Наибольший порядок миноров, который можно составить из эле
ментов матрицы А, |
равен |
£0 = min |
(tn,ri). |
|
|||
Для любого k< kn можно ппррлелить общее число миноров 6-го |
|||||||
порядка по формуле|С* • С „,/где |
С[— число сочетаний из i элемен |
||||||
тов по /. |
|
' |
|
|
|
|
|
Например, для матрицы |
|
|
|
|
|||
|
|
|
2 4 —5 6 |
|
|||
|
|
|
3 4 |
5 7 |
|
||
|
|
|
2 |
1 |
3 |
7 |
|
можно составить 18 миноров второго |
порядка |
||||||
(С\ ■С§= 6-3=18); например, |
|
|
|
|
|||
12 4] |
4, |
2 —5 |
= |
25, |
12 6 |
= — 4 и т. д.; |
|
| 3 4 J — |
|
3 |
5 |
|
|
|3 7 |
|
^— 75
четыре минора |
третьего |
порядка |
|
(Сл> • С |=4 -1 =4); например, |
||||||
2 4 - 5 |
= — 17, |
|
2 |
4 |
6 |
|
||||
3 4 |
5 |
|
3 |
4 7 |
|
|||||
2 |
1 |
3 |
|
|
|
2 |
|
1 |
7 |
|
Рассмотрим еще один пример: |
|
|
|
|
|
|||||
|
|
|
|
2 |
4 |
2 2 |
4 |
|
||
|
|
|
А = |
3 2 |
1 1 4 |
* |
||||
|
|
|
—1 4 2 4 2 |
|||||||
|
|
|
|
4 |
10 5 7 |
10 |
|
Все пять миноров четвертого порядка этой матрицы равны ну лю, но имеется минор третьего порядка
2 2 4 |
3 - 2 ) = |
0 0 |
—4 |
(_ 4 ).(_ 1 )ьИ . |
1 1 |
|
|
1 1 |
4 |
1 1 |
4 |
— — 8 , |
|||
2 4 |
2 |
|
2 4 |
2 |
|
2 4 |
|
|
|
|
|
не равный нулю. Порядок этого не равного нулю минора и есть ранг матрицы А.
Сорядок минора, не равного нулю.
Из этого определения следует, что ранг матрицы А равен R,максимальныйО п р е д е л е н и е . Рангом матрицы называется
если: |
имеется |
по крайней мере один минор матрицы А порядка R, |
||
1) |
||||
не равный нулю; |
матрицы А порядка |
Я+1 и выше равны нулю. |
||
2) |
все миноры |
|||
П р и м е р ы. Определить ранг данных матриц: |
||||
|
|
3 |
3 |
9 |
|
|
4 |
1 |
9 |
|
|
5 |
2 |
12 |
|
|
6 |
4 |
16 |
Р е ш е н и е . Вычислим все миноры третьего порядка:
3 |
3 |
9 |
3 |
3 |
9 |
= о, |
3 |
3 |
9 |
4 |
1 |
9 |
4 |
1 9 |
= 0, 4 1 9 |
5 |
2 12 |
= 0 и 5 |
2 |
12 |
|||||
5 |
2 |
12 |
6 |
4 |
16 |
|
6 |
4 |
16 |
6 |
4 |
16 |
76
Так как имеется минор второго порядка |
2 2 |
=24, |
не равный |
||||
4 16 |
|||||||
нулю, то ранг матрицы А равен двум. |
|
3 |
|
|
|
||
2 4 |
1 |
4 |
|
|
|
||
1 |
2 |
0,5 |
2 |
1,5 |
|
|
|
4 8 |
2 8 |
6 |
|
|
|
||
Р е ш е н и е . Все миноры |
третьего |
и второго |
порядков равны |
||||
нулю. У матрицы А имеется минор первого порядка, не |
равный ну |
лю, например |3|=3. Это означит, что ранг матрицы А равен еди
нице.
Можно доказать следующие свойства миноров матрицы.
1. Если все миноры k-ro порядка матрицы А равны нулю, то все миноры более высоких порядков также равны нулю.
2. Если минор D k-ro порядка не равен нулю, а все миноры (k+\)-ro порядка, полученные окаймлением этого минора, равны нулю, то и все миноры (ft-fl)-ro порядка равны нулю (значит ранг матрицы равен /г).
Рассмотрим второе свойство на примере. Пусть
2 |
4 |
1 |
3 4 |
1 |
А = 3 1 |
—4 1 2 |
3 . |
||
7 9 |
—2 7 10 5 |
Возьмем минор второго порядка
О = 11 f I = — 10^=0.
Вычислим все миноры третьего порядка, получаемые окаймлени ем минора D, таких миноров можно составить четыре:
2 4 |
1 |
2 4 3 |
2 4 4 |
2 4 1 |
3 1 —4 = 0, |
3 1 1 = 0, |
3 1 2 |
= 0 и 3 1 3 |
|
7 9 —2 |
7 9 7 |
7 9 10 |
7 9 5 |
Так как все эти миноры равны нулю, то на основании свойства 2 все остальные 16 миноров третьего порядка равны нулю. На пример,
3 4 |
1 — ( - 2 ) |
3 4 1 |
|||
1 |
2 |
3 |
1 |
2 |
3 |
7 |
10 5 ч- |
1 2 3 |
— //
4 1 3 4 1 3 4 1 3
1 —4 1 —И—0 = 1 —4 1 = 2 1 —4 1 |
||
9 —2 7 «_1 |
8 2 6 |
4 1 3 |
и т. д.
Таким образом, вместо вычисления 20 миноров третьего порядка достаточно вычислить 4, окаймляющие не равный нулю минор вто
рого порядка D.
Ранг матрицы А равен двум.
Для вычисления ранга матрицы целесообразно:
1)переходить от миноров меньших порядков, начиная с мино ров первого порядка, к минорам больших порядков;
2)если найден минор 6-го порядка D, не равный нулю, тогда
следует вычислять миноры (6+1)-го порядка, получающиеся окай млением минора D;
3) вычисление таких миноров нужно продолжать до тех пор, пока среди них не найдется минор, не равный нулю. В этом случае ранг матрицы не меньше, чем (6+1), и нужно повторить предыду щее правило. Если все эти миноры равны нулю, то ранг матрицы
равен 6. |
|
Определить ранг матрицы |
|
|
|||||||
Пр и м е р . |
|
|
|||||||||
|
|
|
А = |
1 2 |
5 |
1 2 |
|
|
|||
|
|
|
2 4 |
1 1 1 . |
|
|
|||||
|
|
|
|
|
3 6 6 2 3 |
|
|
||||
Р е ш е н и е . |
Минор |1|= 1 |
|
первого порядка не равен нулю. Рас |
||||||||
смотрим миноры второго порядка, окаймляющие этот минор: |
|||||||||||
|
|
|
| 2 4 | = ° ’ | 2 f | ^ |
8 ^ °' |
|
||||||
Мы нашли минор второго |
порядка, не равный нулю. Начинаем |
||||||||||
вычислять |
миноры |
третьего |
порядка, |
|
|
1151 |
|||||
окаймляющие минор L j : |
|||||||||||
|
1 |
5 |
2 |
1 |
5 |
1 |
|
|
1 |
5 |
2 |
|
2 1 4 = 0, |
2 1 1 = 0 и 2 1 1 |
|||||||||
|
3 6 6 |
3 6 2 |
|
|
3 6 3 |
||||||
Других |
миноров третьего |
порядка |
вычислять не нужно: из (7) |
||||||||
и (8) следует, что они равны нулю. |
Ранг матрицы равен двум. |
78