книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов
.pdf§ 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
