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

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

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

§ 5. Линейная зависимость столбцов матрицы

Дана матрица

3 2 2 14 —6 А = - 4 7 —1 12 4 5 2 8 24 3

14

Столбец Во = 12 можно линейно выразить через остальные

24

столбцы матрицы

2

 

 

2

 

—6

3

 

 

 

—4 , в 2 = 7

,

Я3 =

—1 и Bt —

4

5

2

 

 

8

 

3

Это значит, что существуют такие числа сц, а2,

а3

и а4, что В0 =

= eti • Bi + a2 • В2+ а 3 • В3+ а 4 • В4. Действительно,

 

 

 

В0 — 2• Bi + 3- В2 + 1*В3 + 0• В4,

 

 

так как

 

 

 

 

 

 

14

3

2

 

21

—6

12 = 2- —4 +3- 7

+ 1-

—1 + 0-

 

4

24

5

2

 

8 1

 

3

или

 

2-3

+ 3 - 2 +

1-2

+ 0- ( —6)

14

=

12

2 - ( - 4 ) + 3 .7 +

1.(—1)+0-4

24

 

2-5

+ 3 -2 -1 1-8-

+ 0 - 3

Оп р е д е л е н и е . Если один из столбцов матрицы можно ли­ нейно выразить через другие, то столбцы матрицы называются ли­ нейно зависимыми.

В рассмотренном примере столбцы матрицы линейно зависимы.

Оп р е д е л е н и е . Если никакой столбец линейно не выражается

через остальные, то такие столбцы называются линейно независи­ мыми.

В матрице

1 О О

0 4 0

2 5 3

столбцы линейно независимы. Действительно:

79

1)

нельзя найти такие числа cti и а2, что

 

 

1

0

0

ах-0 + a 2 -0

=

0

0

= a r 4

a 2‘ 0

= c^-4 4- a 2-0

a r 4

2

5

3

(х^ • 5 ~\~ct2■3

 

ctj ■5 -)- ct-2 • 3

так как при любых cti и а2 элемент первой строки суммы равен ну^

лю, а не единице;

 

 

и а2, что

2) нельзя найти такие числа ai

0

 

1

0

4

= ctr

0

+ a 2' 0

5

 

2

3

так как правая часть равна

ar l

О

сс^ • 2 -j- ct2 *3

и при любых Oi и о2 элемент второй строки суммы равен нулю, а не четырем;

3) третий столбец нельзя выразить через первые два, так как

1

0

 

ax. 1

«1- 0 + a2- 4

=

a =-4

2

5

 

ctj • 2 -j- g2 *5

Поэтому, чтобы в полученном столбце элементы первой и вто­ рой строк равнялись нулю, необходимо, чтобы ai = a2=0, но тогда элемент третьей строки также будет равен нулю, а не трем, как в

третьем столбце.

Итак, мы убедились, что ни один из столбцов рассматриваемой матрицы линейно не выражается через остальные. Это значит, что столбцы матрицы линейно независимы .

Можно доказать следующую теорему.

Т е о р е м а . Определитель D п-го порядка равен нулю тогда и только тогда, когда его столбцы линейно зависимы.

П р и м е р .

3

4

18

3

4

18

3

7

21

7

21

2

—2 —2 = 2-

1 —1

—1 =2

1

0

0

4 12 = 0 .

1

3

11

1

3

11

1

4

12

 

 

 

 

 

II

t

t

 

 

 

 

 

SO

Так как определитель равен нулю, то его столбцы линейно за висимые. Действительно,

18

3

+ 3-

4

—2 = 2-

2

—2

11

1

 

3

значит третий столбец линейно выражается через первые два.

§ 6. Условия существования решения системы линейных уравнений

В этом параграфе приводится условие, которому должны удов­ летворять коэффициенты системы линейных уравнений, чтобы систе­ ма была совместной.

Дана система т линейных уравнений с п неизвестными:

( au Х 1 + а Х2 х 2 + • • • + а Хп х п = Ь \ ,

ат\ Х 1 + ат2 х 2 + ■• • + От п Х п —■Ьт.

Обозначим через А матрицу

аи •• О-Хп

( 10)

атг - а т п

ичерез А расширенную матрицу

а1 1 - ■а\п bi

( И )

amV атп Ь,п

которая получается из матрицы А добавлением столбца свободных членов. Столбцы матрицы А обозначим через Кi, Кг,. . Кп, а стол­ бец свободных членов через В.

Рассмотрим матричное уравнение

Ki ‘xi -\~ Кг 'хг 4- • • • ~Ь Кц' хц В ,

( 12)

6—440

— til —

Э т о у р а в н е н и е м о ж н о з а п и с а т ь в д р у г о м в и д е :

 

Х 1 4“ °12 х 2 4" •

• 4- а 1п

b i

 

а и

х п

(1 3 )

 

 

 

сс

О- m l

* i + am2 *2 4- • * * ~Ь а т п %п

 

 

 

 

bn

 

Уравнение (12) и система

(9)

эквиваленты,

т. е. любое ре­

шение уравнения (12) является решением системы (9), и наоборот. Если система (9) несовместна, то и уравнение (12) не имеет решения.

Если X),

х2, ..., х п являются решением системы (9), т. е.

систе­

ма (9) совместна, то уравнение

(12)

обращается в тождество

 

 

Ki-Xi -f- К2 2

4~ • ■■4* Кп-хп В.

(14)

Следовательно, столбец В линейно выражается через столбцы

матрицы А

И, наоборот, если

известно, что В линейно выражается

через Ки К2, .... Кп, то это значит,

что существуют такие

числа

хи х2, .... л„, для которых имеет место (14).

Следовательно, исходя из эквивалентности системы (9) и урав­

нения (12), хи х2, ..., х п — решение системы (9).

Итак, система (9) совместна тогда и только тогда, когда стол­ бец свободных членов В линейно выражается через столбцы матри­

цы А.

матрицы с числом линейно независимых

Свяжем понятие ранга

столбцов.

 

 

 

Для этого рассмотрим систему

 

3х] -{- 2х2 —р 5х3 — Х \т-- 8,

|xi 4х2— Зх3 + Зх4 = — 2,

4х4 — 2 х2 -j- 2х3 4“ 2х4 = — 1.

3

2

5

1

Ранг матрицы А = 1 —4 —3

3 равен двум.

4 —2

2

2

Это значит, что имеется минор второго порядка не равный нулю, а все миноры третьего порядка равны нулю. (Проверьте это утверж­ дение путем вычисления окаймляющих миноров третьего порядка.)

На основании вышеприведенной теоремы это значит, что любые три столбца линейно зависимы, но имеются два столбца линейно Не­ зависимых.

82

У б е д и м с я в э т о м . Т а к к а к

 

 

 

6

 

 

7 К\

Кз

И

 

 

Kt — — ^

^ Кз,то

 

любые три столбца линейно зависимы,

но так как минор |j

_g|

= —14=И=0, то, следовательно, столбцы

Ki и Кз, составляющие

этот

минор, линейно независимы.

Выше рассмотренный пример приводит нас к важному выводу: ранг матрицы А равен максимальному числу линейно независимых столбцов этой матрицы.

Вопрос о совместности системы линейных уравнений полностью решается следующей теоремой.

Теорема Кронекера — Капелли.

Система линейных

(9) тогда и только тогда совместна,

когда ранг расширена

цы А равен рангу матрицы А.

 

Д о к а з а т е л ь с т в о . Пусть система совместна. Докажем, что

ранг матрицы А равен рангу расширенной матрицы.

 

Ранг матрицы А может быть либо равен рангу матрицы А, либо

на единицу больше последнего.

 

Если R — ранг матрицы А, то среди ее столбцов Кь ..., Кп мож­

но найти не более R линейно независимых столбцов. Обозначим их

К

K in ■ Эти столбцы будут линейно независимы и в матрице

h

R

А. Остальные столбцы матрицы А линейно выражаются через столб-

цы Kt ,..., Ki R.

Матрица А отличается от матрицы А только столбцом свобод­ ных членов. Так как система (9) по предположению совместна, то, как было установлено выше, столбец В линейно выражается через столбцы матрицы А. Это означает, что любой минор (Д+1)-го по­

рядка в матрице .4 равен нулю, т. е. ранги матриц А и А равны. Пусть, с другой стороны, ранг матрицы А равен рангу расширен­

ной матрицы А. Докажем, что система совместна. ЕслиД п,--, Ki^—

линейно независимые столбцы матрицы А, то эти же столбцы явля-

6*

83

ются линейно независимыми в матрице А. Любые /?+ 1 столбцов в матрице А линейно зависимы. Рассмотрим столбцы /<,■ , . . Ki^, В.

Они линейно зависимы. Из линейной независимости Ki

1

. . . . Ki сле-

 

R

дует, что столбец В линейно выражается через столбцы матрицы .4,

т. е. система совместна.

 

 

 

 

 

 

 

П р и м е р ы .

1) Решить систему

 

 

 

 

 

 

 

2*1 -(“ *2

*3

*4

 

4 ,

 

 

 

—*1 + 2*2 +

*3

2*4 =

1,

( 15)

 

*1 -|- 3*2

*4

 

5 .

 

Р е ш е н и е .

 

2

 

 

 

 

 

 

1. Определяем ранг матрицы

2

1

 

—2 . Так как все миноры

—1

 

 

 

1 3

n

0 —1

 

третьего порядка

 

 

=

I

 

2 1I ,

не равен нулю,

равны нулю и минор и

 

_j g =5

то ранг матрицы А равен двум.

 

_

2 1 —1

144

 

 

 

—1 2

1

—2 1

• Он также

 

 

1

3

0 —1

5

 

равен двум. Поэтому система совместна.

 

первого

и второго урав

3.

Минор D состоит из коэффициентов

нений,

стоящих при первом и втором неизвестных, поэтому отбросим

уравнения, коэффициенты которых не входят в минор D, т. е. третье уравнение. Перенесем неизвестные *з и *4, коэффициенты которых не входят в минор D, в правую часть, получаем систему

 

 

 

2*1 ~\~

X 2 — 4 —f—*3 —

* 4,

(16)

 

 

 

-- *1

2*2 = 1

-- *3 “Ь 2*4 ,

 

 

 

 

Считая, что

|4 + *4

является столбцом свободных членов,

 

■*3 4- 2*4

найдем решение системы (16)

по правилу Крамера:

 

 

I 4 -р *з — *4 1 I

41

 

1

 

 

*! =

| 1 — *3 +

2*4 2 |

121

+

 

*з •

*4 =

 

D

 

5

5

 

 

 

 

 

 

 

 

 

 

3

*з —

4

 

 

 

 

 

 

£

7

х*

 

 

 

 

 

 

--

 

 

 

 

 

84

 

 

 

и

2 4

I 2 11

I 2 —1 [

2 4 -|- Х3 Xi

1 1 — Х'\ —j—2х,

I 1

I—1 —1|

| —1 2 I

х2 =

5

5 *3 +

I

D

6

1

3

 

= 5 “ Т Х з + Т Х4-

 

Таким образом, мы получили равенства:

 

7

3

 

4

 

*1 =

5 +

5

^ " Т

*

4'

*2

6

 

1

3

* 4.

с

 

г *3 +

г-

 

5

 

5

5

 

которые определяют общее решение заданной системы.

Задавая произвольные числовые значения неизвестным х3 и дг4, можно получить различные решения заданной системы. Например,

при

 

 

 

 

 

2

 

1

 

 

 

х3 = 40 и х4 =

20 полу чае м xj =

 

и х.> =

,

 

 

9 —

5 —

 

 

 

 

 

 

 

5

 

5

 

 

 

при х3

= xt = 0

 

7

и х2=

6

 

 

 

 

х1== —

— .

 

 

 

 

 

 

 

5

 

 

5

 

 

 

2)

Решить систему

 

 

 

 

 

 

 

 

 

 

2х3 -4* х2 — х3 — 2,

 

 

 

 

 

 

—Xl — 4X2 "i

Xj

 

1 ,

 

 

 

 

 

 

х4 — Зх2

— 3,

 

 

 

 

Р е ш е н и е . Ранг матрицы

2

 

1

■1

равен двум и минор

—1 —4

1

 

 

 

1

—3

0

 

2

1—12

D =

—4 1 = 3 f 0 .

Ранг расширенной

матрицы

—1 —4

11

 

—3 0

 

 

 

 

 

 

1 —3

0 3

 

 

 

 

 

 

 

 

также равен двум, и поэтому система совместна. Минор D состоит из элементов второго и третьего уравнений и второго и третьего столб­ цов, поэтому рассматриваем систему

85

 

 

f

--- 4X2 +

* 3

1

+

Х 1 >

 

 

 

{ — З х 2

 

— 3 — X i .

 

| 1 + * 1 11

n i l

 

1

I

l

l

 

13 0 I

 

1 о|

 

13 — * 1 0 1

 

 

1—

4

11 ~

3

+

 

3

 

Л 1 ~

 

I —

3

0 1

 

 

 

 

 

 

 

I — 4 1

 

Х \ 1

I — 4

1

1

 

1— 4

1 1

1— 3 3 — Х \ 1

1 — 3 3

!

I

— 3 — 1

—4 1

 

 

 

1

 

' г .

 

 

 

 

 

 

 

 

- 3

О

 

 

 

 

 

 

 

 

Способ составления общего решения системы линейных уравне­ ний, рассмотренный в предыдущем примере, можно сформулировать следующим образом.

1. Находим минор D максимального порядка не равный нулю в матрице А системы (9) и определяем ранг этой матрицы.

2.Определяем ранг расширенной матрицы.

3.Если эти ранги не равны между собой, то,система несовмест­

на — прекращаем вычисления.

4.Если ранги равны, то составляем систему, определитель ко­ торой совпадает с минором D.

Для этого из заданной системы отбрасываем все уравнения, ко­ эффициенты из которых не входят в минор D. В правую часть остав­

шихся уравнений переносим все неизвестные, коэффициенты которых

не входят

в минор

D. Если

R — ранг матрицы А,

то

получаем си­

стему, состоящую

из R уравнений и содержащую

R

неизвестных.

Остальные

n—R неизвестных

перенесены в правую

часть каждого

уравнения

и вместе с числами bt рассматриваются как свободные

члены.

-

5. Полученную систему решаем по формулам Крамера.

Таким

образом, мы получаем общее решение системы (9), поз­

воляющее выразить значение R неизвестных через другие п—R не­

известных.

Если R = n, то система имеет единственное решение.

§ 7. Вычисление определителей методом исключения (метод Гаусса)

Определители приходится вычислять не только при решении си­ стем линейных алгебраических уравнений. Многие технические задачи связаны с численным вычислением определителей высоких порядков.

86

Вычислять такие определители вручную очень трудно, а для порядка больше 10—15 практически невозможно. Метод исключения позволя­ ет использовать вычислительные машины для выполнения необходи­ мых вычислений. Этот метод основан на применении изложенных ра­ нее свойств определителей. Он заключается в преобразовании задан­ ного определителя к такому виду, чтобы на главной диагонали всюду стояли единицы, а все элементы, расположенные ниже главной диаго­ нали, были равны нулю.

Рассмотрим этот метод на примере вычисления следующего оп­

ределителя:

 

 

2

—3

7

D = 5

—1 8

3

4

6

Представим этот определитель в виде

 

 

 

 

 

 

 

1

с1 2

Cl3

 

 

 

 

D = b i bi Ь3 0 1

 

с23

 

 

 

 

 

 

0

0

 

1

 

 

тогда значение определителя

D равно произведению Ьф2 Ьз, так как

значение определителя с единицами по главной диагонали и нулями

ниже диагонали равно единице.

 

 

 

 

 

 

Нули в первом столбце во второй и третьей строках можно по­

лучить следующим

образом.

 

 

 

 

 

 

 

а)

элемент ац = 2 вынесем за знак определителя, для чего разде­

лим все элементы первой строки на два, тогда

 

 

 

 

D = 2-

1

—1,5

3,5

 

 

 

 

5 —1

 

8

 

 

 

 

 

 

3

4

 

6

 

 

 

б)

умножим первую строку на числа —5 и —3 и прибавим соот­

ветственно ко второй и третьей строкам; тогда

 

 

 

= 2-

,-- г_ 1 —1,5 3,5

=2-

1

—1,5

3,5

 

L 5 —1

8

0

6,5 —9,5

 

 

(—5)

 

 

 

 

0

8,5

—4,5

 

(-3)1-----►3 4 6

 

 

 

 

 

Таким образом получен определитель, у которого в первом столб­

це ниже диагонали стоят нули;

 

 

 

 

 

 

в)

«сделаем нули» во втором столбце.

 

 

 

 

 

 

-

87

-

 

 

 

 

Для этого вынесем за знак определителя элемент 022—6,5, по­

лучим

1

- 1 ,5

3,5

 

 

;

Dt

2-6,5- О

1

—1,46

 

О

8,5

—4,5

 

г) чтобы «сделать нуль» во втором столбце ниже диагонали,

вычтем из третьей строки вторую,

умноженную на «—8,5», получаем

Д , =

2-6,5-

-1 ,5

3,5

;

1

—1,46

 

 

О

7,92

 

д) вынесем за

знак определителя элемент я33=7,92; получим

1

- 1 ,5

3,5

О

1

—1,46 -2-6,5-7,92 = 2-6,5-7,92 = 102,96.

О

0

1

Промежуточные результаты, которые получились в составленных выше определителях, целесообразно заносить в специальную таблицу

(табл.

1). Эта таблица состоит из 12 столбцов и 9 строк (для опреде­

лителя

л-го порядка — 2л + 6 столбцов

и л2 строк).

Эта таблица разделена на три зоны.

Каждая зона состоит из трех

строк.

Первые четыре столбца называются дополнительными. Каж­

дый нечетный столбец, начиная с пятого, называется основным, а все остальные — вспомогательными. Таким образом, составленная табли­ ца содержит четыре дополнительных, четыре основных и четыре вспо­ могательных столбца.

Каждый основной столбец и стоящий справа от него вспомога­ тельный называются связанными.

Элементы рассматриваемого определителя записываются в пер­ вых трех основных столбцах в строках в I зоне.

Для определителя п-то порядка таблица состоит из л2 строк, разделенных на л зон, которые соответствуют строкам опреде­ лителя.

Таблица содержит л+1 основных, л+1 вспомогательных и четы­ ре дополнительных столбца. Элементы определителя записываются в первых л основных столбцах I зоны.

В последнем основном столбце (предпоследнем столбце в табли­ це) записываются контрольные числа, равные сумме элементов исход­ ного определителя, расположенных в соответствующих строках.

В строках И зоны основных столбцов записываются элементы первого промежуточного определителя и контрольные суммы его эле-

-86

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