Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов

.pdf
Скачиваний:
20
Добавлен:
29.10.2023
Размер:
7.59 Mб
Скачать

На основании свойства 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

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

Соседние файлы в папке книги из ГПНТБ