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

МФТИ 2012 Умнов АЕ АГ+ЛА Стр001-544

.pdf
Скачиваний:
39
Добавлен:
03.06.2015
Размер:
6.64 Mб
Скачать

Глава 6 . Системы линейных уравнений

199

что в силу теоремы 6.2.1 означает

det ( A B ) = det A × det B .

Теорема доказана.

§ 6.3. Разложение определителей

Выберем в

квадратной матрице n -го порядка

 

A

строки с но-

мерами i1 , i2

, ... ,ik и столбцы с номерами

j1 , j2 ,... , jk , где

1 k n .

 

 

 

 

Определение

Детерминант квадратной матрицы порядка k , обра-

6.3.1.зованной элементами, стоящими на пересечении

 

строк

i1 , i2 , ... ,ik

и столбцов j1 , j2 ,... , jk , назы-

 

вается

минором

k -го

порядка

и обозначается

 

j1, j2 ,..., jk

.

 

 

 

 

 

M i

,i

,...,i

 

 

 

 

1

2

k

 

 

 

 

 

 

 

 

 

Определение

Детерминант квадратной

матрицы

порядка n k ,

6.3.2.образованной элементами, остающимися после вы-

черкивания строк

i1 , i2 , ... ,ik

 

и столбцов

 

j1 , j2 ,... , jk , называется

минором,

дополнитель-

ным

к

 

j1 , j2 ,..., jk

,

и

обозначается

минору M i

,i

2

,...,i

k

 

 

 

 

 

1

 

 

 

 

 

 

 

 

j1, j2 ,..., jk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M i

,i ,...,i

k

.

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

Выберем в матрице A i -ю строку и j -й столбец, на пересече-

нии которых расположен элемент aij . Удалим из A выбранные

строку и столбец, рассмотрим квадратную матрицу A+ размера

(n 1) × (n 1) .

 

200

 

 

 

 

 

 

 

 

 

Аналитическая геометрия и линейная алгебра

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение

 

Детерминант матрицы

 

A+

 

называется дополни-

 

 

 

6.3.3.

 

 

 

 

 

 

 

 

 

тельным минором

 

ij

элемента αij .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

Сгруппируем в определении 6.1.2 –

детерминанта матрицы

 

A

 

 

 

 

 

 

 

 

 

все (n −1)!

слагаемых, содержащих элемент

αij , и вынесем

 

его

 

за

 

 

 

скобки. Получим выражение вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det

 

A

 

 

 

= αij Dij +K.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение

 

Число

Dij

 

называется алгебраическим дополнением

 

 

 

6.3.4.

 

 

 

 

 

 

 

 

 

элемента αij .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что по определению 6.1.2 имеют место равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

det

 

 

 

A

 

 

 

= αij Dij

= α kj Dkj j = [1, n], i = [1, n],

(6.3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j =1

k =1

 

 

 

 

 

 

 

 

 

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

 

 

Справедливы равенства Dij = (−1)i + j

 

ij .

Теорема

M

6.3.1.

 

 

 

 

 

 

 

 

 

 

 

Доказательство.

 

 

 

 

 

1°. По определению детерминанта 6.1.2

 

 

 

 

det

 

A

 

 

 

= α11

(−1) Б(1,k2 ,k3 ,...,kn ) α 2k2

...α nkn

+ K ,

 

 

 

 

 

 

 

 

 

 

{1,k2 ,k3 ,...,kn }

 

 

 

 

то

есть

 

D11 =

(−1)Б(k2 ,...,kn ) α 2k2 α 3k3

...α nkn

, поскольку

 

 

 

 

 

 

 

{k2 ,...,kn }

 

 

 

 

очевидно, что Б(1, k2 , k3 ,..., kn ) = Б(k2 , k3 ,..., kn ) ,

 

Глава 6 . Системы линейных уравнений

 

 

 

201

 

но тогда выражение для

D11 совпадает с формулой определи-

 

 

 

теля матрицы порядка n 1 ,

получаемой из

 

A

 

вычеркива-

 

 

 

 

нием

 

 

первого столбца

и

первой строки.

Следовательно,

 

 

=

 

 

 

1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2°. Построим новую матрицу

 

A

 

 

 

, переместив элемент αij

 

мат-

 

 

 

 

 

 

рицы

 

 

A

 

в ее левый верхний угол, переставив i-ю строку на

 

 

 

 

 

первое

 

место, для чего потребуется i 1 перестановок, и пере-

 

ставим на первое место j-й столбец, что потребует

j 1 пере-

 

становок. Тогда определитель перестроенной матрицы

 

A

 

 

 

 

 

 

 

 

 

равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det A = (1)i 1+ j 1 det A = (1)i + j det A .

Согласно линейному свойству определителя (теорема 6.2.3) данное соотношение будет также выполняться и для каждого из его слагаемых, а значит, в силу формул (6.3.1) и для каждого

алгебраического дополнения Dij . Поэтому справедливо равен-

ство

D

= (1)i + j D.

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

11

 

 

 

 

 

 

 

 

 

 

3°. Наконец, очевидно, что значение дополнительного к

αij ми-

нора не зависит от положения αij

в матрице

 

 

 

A

 

 

 

,

и потому

 

 

 

 

 

 

ij

=

M 11 .

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

4°. Учитывая полученные соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij =

M 11 = D′ = (1)i+ j D ,

 

 

 

 

 

 

 

M

 

 

 

 

 

 

11

ij

 

приходим к равенству Dij = (1)i+ j M ij .

Теорема доказана.

202

Аналитическая геометрия и линейная алгебра

Следствие

Разложение определителя по i -у столбцу имеет вид

6.3.1.

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

det

 

 

A

 

= (-1)k +i aki

 

 

 

 

 

 

ik

 

 

 

 

 

M

 

 

или

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

+i M ki

 

 

 

 

 

 

ik .

 

det

 

 

A

= (-1)k

M

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

Для практических приложений особо полезной является

Теорема

Для любой квадратной матрицы

 

A

 

 

 

 

имеет место

 

 

 

6.3.2.

равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aij Dis = d js × D ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

где D = det

 

A

 

 

 

1 ,

j = s,

символ Кро-

 

 

 

 

 

 

 

 

и d js =

 

 

 

 

 

 

 

 

 

 

0 ,

j ¹ s

 

 

 

 

 

 

 

 

 

 

некера (см. § 2.2).

Доказательство.

По определению 6.3.4 алгебраического дополнения имеем det A = a1 j D1 j + a2 j D2 j +K + anj Dnj , то есть утвер-

ждение теоремы для случая j = s справедливо.

Пусть теперь j ¹ s . Тогда выражение

a1 j D1s + a2 j D2s + ... + anj Dns

можно рассматривать как разложение по s-му столбцу опреде-

лителя матрицы, у которой s-й столбец совпадает с j-м столб- цом. Но такой определитель равен нулю по следствию 6.2.2.

Теорема доказана.

 

Глава 6 .

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

 

 

 

 

 

 

 

 

 

 

 

 

203

 

Следствие

 

Если квадратная матрица

 

A

 

 

 

невырожденна, то

 

 

 

 

 

 

6.3.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

1

 

 

элементами ее обратной матрицы

 

 

являются

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(-1)i+ j

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

числа bij =

M

 

; i, j = [1, n] .

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство.

 

 

 

 

A

 

 

 

B

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем произведение матриц

 

 

 

и

 

 

элементы которых

 

 

 

 

 

 

 

ai j

 

и bi j ; i, j = [1, n] . Пусть

g pq

 

элемент произведения

 

 

A

 

и

 

B

 

, тогда, согласно определению 5.1.1

и теореме

 

 

 

 

 

6.3.2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

(-1) j +q

 

qj

1 n

 

 

 

 

 

 

 

 

M

 

 

 

g pq = a pj b jq =a pj

 

 

 

 

 

 

 

=

 

 

a pj Dqj =

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

j =1

j =1

 

 

 

 

 

 

 

 

 

 

 

 

D j =1

 

 

 

 

=

 

 

1

× D × d pq

= d pq ;

 

i, j = [1, n] .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогичное соотношение получается и для произведения

B A и по определению 1.1.4

A B = B A = E ,

но тогда, согласно определению 5.1.2 и лемме 5.1.1,

B = A 1 .

Следствие доказано.

Проверьте самостоятельно справедливость формулы (5.1.1).

Обозначим I = i1 + i2 + ... + ik и J = j1 + j2 + ... + jk , тогда оказывается справедливой обобщающая следствие 6.3.1

204

Аналитическая геометрия и линейная алгебра

Теорема

Для фиксированного набора столбцов j1,

j2 ,..., jk

6.3.3

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

 

 

 

 

 

 

 

(Лапласа).

I + J

 

 

 

 

j1 , j2 ,..., jk

 

det

 

=

(-1)

j1 , j2 ,..., jk

 

 

 

 

 

A

M i1 ,i2

,...,ik .

 

 

M i

,i ,...,i

 

 

 

 

 

 

 

1

2

k

 

 

 

{i1 ,i2 ,...,ik }

Отметим, что суммирование выполняется по всем возможным пе- рестановкам номеров строк i1 , i2 , ... ,ik .

Задача

6.3.1.

Решение.

Найти определитель матрицы n-го порядка:

 

x

a

a ...

a

 

 

a

x

a ...

a

 

Dn = det

a

a

x ...

a

.

 

... ... ... ... ...

 

 

a

a

a ...

x

 

 

 

 

 

 

 

1°. Заметим, что сумма элементов каждого столбца матрицы одинакова и равна x + a(n − 1) . Поэто-

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

следствия 6.2.4 и 6.2.3):

 

1

1

1

...

1

 

 

a

x

a

...

a

 

Dn = ( x + a(n - 1)) × det

a

a

x

...

a

.

 

... ... ... ... ...

 

 

a

a

a

...

x

 

 

 

 

 

 

 

 

2°. Вычитая из каждой строки, начиная со второй, пер- вую строку, умноженную на a , получим

Глава 6 . Системы линейных уравнений

 

 

205

 

 

 

 

 

1

1

1

...

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

x - a

0

...

0

 

 

 

 

Dn

= (x + a(n -1)) × det

 

 

 

0

0

x - a

...

0

 

 

 

.

 

 

 

 

 

...

...

...

...

...

 

 

 

 

 

 

 

 

 

0

0

0

...

x - a

 

 

 

 

3°.

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

n 1 раз

следствие

 

6.3.1 для разложения определителя по первому

 

столбцу, приходим к выражению

 

 

 

 

 

 

n= ( x + a(n 1))( x a) n1 .

§6.4. Правило Крамера

Будем рассматривать систему

n линейных уравнений с

n неиз-

вестными:

 

 

α11ξ1 + α12 ξ 2

 

+ ... + α1n ξ n

= β1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ... + α 2n ξ n

= β 2 ,

 

 

 

 

 

 

α 21 ξ1 + α22 ξ2

 

 

 

 

 

 

 

...............................................

(6.4.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ2 + ... + α nn ξ n

= βn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α n1 ξ1 + αn 2

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

α ji ξ i

= β j ;

j = [1, n] или же в матрич-

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ной форме

 

A

 

x

 

=

 

 

 

b

 

 

 

, где

квадратная матрица

 

A

 

имеет

 

 

 

 

 

 

 

 

 

компоненты α ji

, а столбцы

 

 

x

 

и

 

b

 

соответственно ξ i

и β j .

 

 

 

 

 

Определение

Упорядоченный набор чисел {ξ1 , ξ 2 ,..., ξ n } будем

6.4.1.

называть решением системы линейных уравнений,

если при подстановке этих чисел в каждое из уравне- ний системы мы получаем тождество.

206

Аналитическая геометрия и линейная алгебра

Имеет место

Теорема

Для того чтобы система линейных уравнений (6.4.1)

6.4.1имела единственное решение, необходимо и доста-

(правило

Крамера).

точно, чтобы D = det A ¹ 0 , и в этом случае реше-

ние данной системы будет иметь вид

xi = DDi ; i = 1,2,..., n ,

где Di определитель матрицы, получаемой из мат-

рицы A заменой ее i -го столбца на столбец сво-

бодных членов b :

 

a11

a12

...

b1

...

a1n

 

 

 

 

 

 

 

Di = det

a21

a22

...

b2

...

a2n

 

 

 

.

 

... ... ... ... ... ...

 

 

 

 

 

an1

an 2

...

bn

...

ann

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

i -й столбец

Доказательство.

1°. Получим вначале утверждение теоремы в предположении, что

x1

система (6.4.1) имеет решение x = x2 , то есть когда вы-

...

xn

полняются равенства

n

a ji xi = b j ; j = [1, n] .

i=1

 

Глава 6 . Системы линейных уравнений

 

207

 

Умножив последовательно для всех j = [1, n]

обе части этих

 

 

равенств на алгебраическое дополнение

D jk

и просуммиро-

 

вав результаты умножения по

j , получим

 

 

n

n

n

 

 

 

D jk (a ji xi ) = b j D jk

"k = [1, n].

 

j =1

i=1

j =1

 

 

Изменим порядок суммирования (то есть выполним пере- группировку слагаемых) в левой части этого равенства:

n

n

 

 

 

 

n

 

(a ji D jk )x i = b j D jk .

 

i=1 j =1

 

 

 

 

j =1

 

Но выражение в круглых скобках равно D×dik

(по теореме

6.3.2), поэтому, учитывая, что

 

n

 

 

 

 

 

n

 

b j D jk = Dk

и Ddik x i = x k Dk ,

j =1

 

 

 

 

 

i=1

 

получаем x k D = Dk

, k = [1, n] . Или окончательно

x k

=

Dk

 

 

"k = [1, n] .

(6.4.2)

D

 

 

 

 

 

 

Наконец, заметим, что из соотношений

 

 

x k D = Dk , k = [1, n]

 

в предположении D ¹ 0 следует, что решение,

определяемое

формулами (6.4.2), существует и единственно.

 

2°. Докажем теперь, что в условиях теоремы набор чисел

 

{ x i =

 

Dk

, i = [1, n] }

 

 

 

 

 

 

 

 

 

 

D

 

есть решение данной системы линейных уравнений. Убедим-

ся в этом, подставив значения xi в левые части исходной системы линейных уравнений (6.4.1):

208

 

Аналитическая геометрия и линейная алгебра

 

n

 

 

1

 

n

n

 

1

n

n

 

 

 

 

 

 

α ji

i

=

 

α ji (βk Dki ) =

βk (α ji Dki ) =

 

 

 

 

 

 

 

i=1

 

 

 

 

 

i=1

k =1

 

 

k =1

i=1

 

 

 

 

 

1

n

 

 

 

 

 

 

 

=

βk δkj

= β j , j = [1, n].

 

 

 

 

 

k =1

Для получения последнего равенства мы снова изменили по- рядок суммирования и воспользовались теоремой 6.3.2.

Теорема доказана.

§ 6.5. Ранг матрицы

Рассмотрим матрицу A размера m × n . Пусть

1 ≤ k ≤ min{m, n}.

Выберем в A k фиксированных столбцов и строк, на пересечении которых стоит матрица минора порядка k .

Пусть при данном k все миноры k -го порядка равны нулю, тогда будут равны нулю и все миноры порядка выше, чем k , поскольку каждый минор (k + 1) -го порядка представим в виде линейной ком-

бинации миноров порядка k (cм. следствие 6.3.1).

Определение Наивысший из порядков, отличных от нуля миноров

6.5.1.

матрицы A , называется рангом матрицы и обо-

значается rg A .

Определение Любой ненулевой минор матрицы, порядок которого

6.5.2.равен ее рангу, называется базисным минором.

Определение Столбцы (строки) матрицы, входящие в матрицу ба-

6.5.3.зисного минора, называются базисными.