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

Линейная алгебра

.pdf
Скачиваний:
53
Добавлен:
05.06.2015
Размер:
1.13 Mб
Скачать

АЛГЕБРА МАТРИЦ (ПРОДОЛЖЕНИЕ)

Уже этих данных достаточно для того, чтобы сделать вывод: D2=1, следовательно,

E2 = D2 =1.

D1

Определяем D3

 

 

D =

 

λ 2

1

 

0

 

= (λ 2)3 ,

 

 

 

 

 

 

 

 

0

 

λ 2 1

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

λ 2

 

 

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

=

(λ 2)3

= (λ 2)3 .

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, канонической

формой

данной матрицы является следующая

λ -матрица:

 

 

 

 

 

1

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

0

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

0

0

(λ

3

 

 

 

 

 

 

 

2)

 

Матричным многочленом называется выражение вида

F (λ)= A0λS + A1λS 1 +... + AS ,

где λ – переменное; A0 , A1 ,..., AS – квадратные матрицыпорядка n с числовыми элементами. Если A0 0 , то S называют степенью матричного многочлена, n – порядком мат-

ричного многочлена.

Любую квадратичную λ -матрицу можно представить в виде матричного многочлена. Справедливо, очевидно, и обратное утверждение, т.е. любой матричный многочлен можно представить в виде некоторой квадратной λ -матрицы.

Справедливость данных утверждений со всей очевидностью вытекает из свойств операций над матрицами. Остановимся на следующих примерах:

Пример. Представить многочленную матрицу

 

λ2

1

 

 

λ +1

 

 

 

 

 

2

 

 

A(λ)=

λ

+1

λ

+ 2λ

 

 

 

+1

в виде матричного многочлена можно следующим образом

1

0

0

1

 

1 1

 

 

 

λ2

+

 

 

λ +

.

 

0

1

 

 

1

2

 

 

 

 

 

 

 

 

1 1

Пример. Матричный многочлен

2

0

1

1

3

5

G(λ)=

 

 

λ2

+

 

 

λ +

 

 

 

 

0

1

 

 

3

0

 

 

1

2

 

 

 

 

 

 

 

можно представить в виде следующей многочленной матрицы ( λ -матрицы)

 

2λ2 +λ +3

λ

5

 

 

 

 

2

 

.

G(λ)=

3λ +1

λ

2

 

 

 

 

Эта взаимозаменяемость матричных многочленов и многочленных матриц играет существенную роль в математическом аппарате методов факторного и компонентного анализа.

Матричные многочлены одинакового порядка можно складывать, вычитать и умножать аналогично обычным многочленам с числовыми коэффициентами. Следует, одна-

31

АЛГЕБРА МАТРИЦ (ПРОДОЛЖЕНИЕ)

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

Два матричных многочлена называются равными, если равны их коэффициенты, т.е. соответствующие матрицы при одинаковых степенях переменного λ .

Суммой (разностью) двух матричных многочленов F(λ) и G(λ) называется такой

матричный многочлен, у которого коэффициент при каждой степени переменного λ равен сумме (разности) коэффициентов при той же степени λ в многочленах F(λ) и G(λ).

Чтобы умножить матричный многочлен F(λ) на матричный многочлен G(λ), нужно каждый член матричного многочлена F(λ) умножить на каждый член матричного многочлена G(λ), сложить полученные произведения и привести подобные члены.

Степень матричного многочлена – произведения F(λ)G(λ) меньше или равна сум-

ме степеней сомножителей.

Операции над матричными многочленами можно осуществлять с помощью операций над соответствующими λ -матрицами.

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

Пример.

 

1 0

0

1

0

1

 

1 0

F(λ)=

 

λ +

 

 

G(λ)=

 

λ +

 

 

 

1 0

 

 

0

 

 

0

 

 

1 0

 

 

 

 

1

 

1

 

 

 

 

1

0 0 1

0 1

 

0 1

 

 

1 0 1

0

 

F(λ)G(λ)=

 

 

 

λ2

+

 

 

λ

+

 

 

λ +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 0

 

 

 

 

 

1 0

0 1

0

1

0 1

 

 

1

0

 

+

0

1 1

0

0

1

 

1

1

1

0

 

 

 

 

 

 

=

 

λ2

+

 

 

λ +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0

1 1

0

0

1

 

1

1

 

0

 

 

 

 

С другой стороны F(λ)

и G(λ) можно записать в виде

 

 

 

 

 

λ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

λ

F(λ)=

 

 

 

 

 

 

 

и

 

 

 

 

 

 

G(λ)=

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

λ 1

λ2 +λ

 

0

1

1

1

1 0

 

F(λ)G(λ)=

 

 

 

 

 

=

 

 

λ2 +

 

 

λ +

 

 

 

 

 

 

 

 

λ2 λ

 

0

 

 

 

1

 

 

 

 

 

 

 

λ +1

 

1

 

1

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как умножение матриц не коммутативно, для матричных многочленов опреде-

ляются два деления с остатком – правое и левое.

 

 

 

 

 

 

 

 

 

Пусть даны два матричных многочлена порядка n

 

 

 

 

 

 

 

F(λ)= A λS + A λS 1 +... + A

 

 

 

 

 

G(λ)= B λt + B λt1

+... + B

 

0

1

 

 

S

 

 

 

 

 

 

 

0

1

 

 

t

где В0 – невырожденная матрица.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При делении F(λ)

на G(λ)

существует однозначно определенное правое частное

Q1 (λ) и правый остаток R1(λ)

F (λ)= Q1 (λ)G(λ)+ R1 (λ),

где степень R1 меньше степени G(λ), или R1 (λ)= 0 (деление без остатка), а также левое частное Q2 (λ) и левый остаток R2 (λ)

F (λ) = G(λ)Q2 (λ) + R2 (λ),

где степень R2 (λ) меньше степени G(λ), или R2 (λ)=0 (деление без остатка).

32

АЛГЕБРА МАТРИЦ (ПРОДОЛЖЕНИЕ)

Обобщённая теорема Безу. При делении матричного многочлена F(λ) на много-

член (λE A) правый остаток равен правому значению делимого F(λ) при

λ = A, т.е.

матрице

= A AS

+ A AS 1

 

= R ,

 

F

+... + A

(3.4.1)

(пр)

0

1

S

1

а левый остаток – левому значению делимого F(λ) при λ = A, т.е. матрице

 

F( лев) = AS A0 + AS 1 A0 +... + AS = R2

(3.4.2)

Доказательство. Доказательство справедливости обеих формул (3.4.1) и (3.4.2) осуществляется одинаково, непосредственной подстановкой. Докажем одну из них.

Итак, делимое – F (λ) , делитель – G = λE A , в качестве частного имеем многочлен

Q2 (λ) = A0λS 1 +(AA0 + A1)λS 2 +(A2 A0 + AA1 + A2 )λS 3 +... +

+(AS 1 A0 + AS 2 A1 +... + AS 1).

Определим произведение (λE A) Q2 (λ) :

A0λS +(AA0 + A1)λS 1 +(A2 A0 + AA1 + A2 )λS 2 +... +

+(AS 1 A0 + AS 2 A1 +... + AS 1)λ AA0λS 1 (A2 A0 + AA1)λS 2

(A3 A0 + A2 A1 + AA2 )λS 3 ... (AS A0 + AS 1 A1 +... + AAS 1) =

= A0λS + A1λS 1 + A2λS 2 +... + AS λ (AS A0 + AS 1 A1 +... + AAS 1),

т.е. ( AS A0 + AS 1 A1 +... + AAS 1 ) + F (λ) = (λE A)Q2 (λ), или F (λ) = (λE A)Q2 (λ) +( AS A0 + AS 1 A1 +... + AAS 1 ),

т.е. R2 = AS A0 + AS 1 A1 +... + AAS 1,

что и требовалось доказать.

Следствие. F(λ) делится справа (слева) на многочлен (λE A) тогда и только то-

гда, когда F

 

= R (F

 

 

= R )

равно 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(пр)

1

( лев)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Показать, что матричный многочлен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(λ)

 

 

λ2 2λ

2

 

1

0

+

2

 

0

 

 

0

2

 

 

 

 

=

 

 

 

 

 

=

λ2

 

 

 

 

λ

+

 

 

 

 

 

 

 

 

 

 

 

λ

 

2λ

+

 

 

 

 

 

1

 

2

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

2

0

1

 

 

 

 

 

 

 

 

делится на матричный многочлен (λE A) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

 

, слева без остатка.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. В самом деле, справедливо равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (λ)= (λE A)Q(λ)+ R1 , где R1 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ 2

1

 

 

 

λ

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

λE A =

 

 

 

; Q(λ)=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

λ +

 

 

 

 

 

0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Подсчитаем значение левого остатка по теореме Безу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

21 0

 

 

0

2

=

 

 

 

 

 

 

 

 

 

 

 

 

R = A2

 

+

A

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

1 2

1

0

2

1 2

0

 

0

2

=

0

0

 

 

 

 

=

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

1

 

1

2

 

 

 

 

 

2

 

 

 

0

 

 

 

 

 

1

1 0

0

 

1

 

 

0

 

 

 

0

 

33

АЛГЕБРА МАТРИЦ (ПРОДОЛЖЕНИЕ)

3.5.Задания для самостоятельной работы по главе 3

3.1.Найти обратную матрицу

 

1

2

3

4

 

 

2

3

1

2

 

 

 

 

1

1

1

1

.

 

 

 

1

0

2

6

 

 

 

3.2. Найти обратную матрицу порядка n

 

 

 

1

1

 

1 ...

 

 

1

 

 

 

 

 

0

1

 

1 ...

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

1 ...

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

... ...

 

 

 

 

 

 

 

...

 

 

...

 

 

 

 

 

0

0

 

0 ...

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3.3. Найти обратную матрицу порядка n

n 1

 

 

 

 

n

 

 

1

2

3

 

4

...

 

 

 

 

 

 

0

1

2

 

3

...

n 2

 

 

n 1

 

 

 

 

 

 

 

0

0

1

 

2

...

n 3

 

 

n

 

2

.

 

 

 

 

 

 

...

...

 

 

 

 

...

 

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

 

 

 

 

 

 

0

0

0

0

...

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

0

0

0

0

...

0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

3.4. Найти обратную матрицу порядка n

 

 

 

 

 

 

 

 

 

 

 

1

1

 

1 ...

 

 

1

 

 

 

 

 

1

0

 

1 ...

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

0 ...

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

... ...

 

 

 

 

 

 

 

...

 

 

...

 

 

 

 

 

1

1

 

1 ...

 

 

0

 

 

 

 

 

 

 

 

 

 

 

3.5. Найти обратную матрицу

cosα

sinα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

cosα

 

 

 

 

 

 

 

 

 

 

sinα

 

 

 

 

 

 

 

 

3.6. Найти обратную матрицу порядка (n+1)

 

 

 

n

 

 

 

 

 

1

a

a

2

a

3

...

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

a

a2

...

 

an1

 

 

 

 

0

0

1

a

...

a

n2

.

 

 

 

 

 

 

 

 

 

... ... ...

...

...

 

 

...

 

 

 

 

 

0

0

0

0

...

 

 

1

 

 

 

 

 

 

 

 

 

 

34

АЛГЕБРА МАТРИЦ (ПРОДОЛЖЕНИЕ)

3.7. Найти обратную матрицу порядка n

 

2

1

0

0

...

0

 

 

1

2

1

0

...

0

 

 

 

 

0

1

2

1

...

0

.

 

 

 

 

 

...

...

 

... ... ...

...

 

0

0

0

0

...

2

 

 

 

3.8. Как изменится обратная матрица A1, если в данной матрице A: а) переставить i-ую и j-ую строки?

б) i-ую строку умножить на число с, не равное нулю?

в) к i-ой строке прибавить j-ую, умноженную на число с, или совершить аналогичное преобразование столбцов?

3.9. Найти матрицу A

1

, обратную для матрицы

E

k

U

, где

E

 

и E

– еди-

 

 

 

 

 

 

 

 

A =

O

E

 

 

k

 

l

 

 

 

 

 

 

l

 

 

 

 

 

ничные матрицы соответственно порядков k и l, U – произвольная матрица порядка k ×l , а все остальные элементы равны нулю.

3.10. Показать, что операция транспонирования матрицы обладает свойствами:

а) (A + B)= A′+ B;

б) (AB)= BA; в) (cA)= cA;

г) (A1 )= (A)1 ,

где с – число, а А и В – матрицы.

3.11.Доказать, что если А и В – симметрические квадратные матрицы одинакового порядка, то матрица C = A B A B ... A B A является симметрической.

3.12.Показать, что для любой матрицы Вматрица A = B Bявляется симметрической.

3.13. Квадратная матрица A = (aij ) порядка n называется ортогональной, если

A A′ = E , где Е – единичная матрица. Показать, что для ортогональности квадратной матрицы А необходимо и достаточно любое из следующих условий:

а) столбцы А образуют ортонормированную систему, т.е.

n

akiakj =δi j , k =1

где δi j – символ Кронекера, обозначающий 1 при i=j и 0 при i j ; б) строки А образуют ортонормированную систему, т.е.

n

aik a jk =δij .

k=1

3.14.Доказать, что ранг суммы двух матриц не больше суммы их рангов.

3.15.Доказать, что если ранг матрицы А равен r, то минор d, стоящий на пересечении любых r линейно независимых строк и r линейно независимых столбцов этой матрицы, отличен от нуля.

35

РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Глава 4. Решение системы линейных уравнений

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

Системой m линейных уравнений с n неизвестными называется система m алгебраических уравнений первой степени вида

 

 

 

a11x1 + a12 x2 +... + a1n xn =b1,

 

 

 

 

a21x1 + a22 x2 +... + a2n xn =b2

,

 

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

(4.1.1)

 

 

где x j , j =

 

am1x1 + am2 x2 +... + amn xn =bm ,

 

– неизвестные, подлежащие определению;

1,n

aij

– числа, называемые коэффициентами при неизвестных;

bi

– числа, называемые свободными членами.

 

Решением системы уравнений (4.1.1) называется совокупность n чисел α1,α2 ,...,αn

таких, что если в каждое уравнение системы вместо неизвестных подставить эти числа (α1 вместо x1 , α2 вместо x2 ,...,αn вместо xn ), то все уравнения обратятся в тождества.

Если система линейных уравнений (4.1.1) имеет хотя бы одно решение, то она называется совместной. В противном случае система называется несовместной.

Совместная система, имеющая единственное решение, называется определенной, а система, имеющая более одного решения – неопределенной.

Две системы линейных уравнений называются эквивалентными, если любое решение каждой из них является одновременно решением и другой системы.

Две произвольные несовместные системы считаются эквивалентными.

Системе линейных уравнений (4.1.1) поставим в соответствие матрицу A = (aij )m,n и

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

~

A

a11 = a21

a...

m1

a

...

a

b

 

 

 

12

 

 

1n

1

 

 

a22

...

a2n

b2

 

,

...

... ...

...

 

 

 

 

a

m2

...

a

mn

b

 

 

 

 

 

m

 

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

4.2. Методы решения системы n линейных уравнений с n неизвестными

Рассмотрим систему n линейных уравнений с n неизвестными

a11x1 + a12 x2 +... + a1n xn = b1,

 

a21x1 + a22 x2 +... + a2n xn = b2

,

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

(4.2.1)

 

an1x1 + an2 x2 +... + ann xn = bn .

Определитель |A| матрицы А называется определителем системы (4.2.1).

Теорема Крамера. Если определитель |A| системы (4.2.1) отличен от нуля, то система совместна и имеет единственное решение.

36

РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Доказательство. Пусть система (4.2.1) совместна и α1,α2 ,...,αn – одно из ее решений. Тогда получим n тождеств:

a11α1 + a12α2 +... + a1nαn = b1,

 

 

a21α1 + a22α2 +... + a2nαn = b2 ,

(4.2.2)

 

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

 

 

 

an1α1 + an2α2 +... + annαn = bn.

 

 

Умножим обе части первого из равенств (4.2.2) на алгебраическое дополнение А

,

 

1 j

 

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

и т.д. и обе час-

2 j

 

 

ти n-ого равенства – на А . Складывая левые и правые части полученных выражений,

nj

 

придем к следующему равенству:

 

(a11 А1 j + a21 А2 j +... + an1 Аnj )α1 +(a12 А1 j + a22 А2 j +... + an2 Аnj )α2 +

 

+... +(a1 j А1 j + a2 j А2 j +... + anj Аnj )α j +... +(a1n А1 j + a2n А2 j +... +

(4.2.3)

+ ann Аnj )αn = b1 А1 j +b2 А2 j +... +bn Аnj

Коэффициент при α j равен определителю |A| системы (4.2.1), коэффициент при α1,α2 ,...,α j1,α j+1,...,αn равен нулю, а правая часть равенства (4.2.3) является определителем,

полученным из определителя|A| путем замены j-гостолбца столбцом свободных членов. Обозначим данный определитель через

Aj = b1 A1 j +b2 A2 j +... +bn Anj

Тогда равенство (4.2.3) примет вид: | А|α j =| Aj |, откуда

 

| Aj

|

 

 

 

 

 

 

 

 

α j =

, j =1, n

 

 

 

 

(4.2.4)

| A |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из формулы (4.2.4) следует, что если система (4.2.1) совместна, то она обладает

единственным решением.

 

 

 

 

 

 

 

 

 

 

Формулы (4.2.4) называются формулами Крамера.

| Aj |

 

 

 

 

Непосредственной подстановкой

значений α j =

, j =

 

, во

все уравнения

1, n

| A |

 

 

 

 

 

 

 

 

 

 

системы убедимся в том, что они образуют ее решение:

n

| Aj

|

 

1

 

n

 

 

 

1

 

 

 

n

n

 

 

1

 

n

n

aij

 

 

 

=

 

 

 

aij | Aj | =

 

 

 

 

 

aij bk Akj

=

 

 

 

 

bk aij Akj .

| A

|

 

 

 

 

 

 

 

 

 

 

 

 

 

j=1

 

 

| A | j=1

 

 

 

| A | j=1

k=1

 

 

| A | k=1

j=1

При i = k ,

n

 

 

 

 

 

 

 

 

 

 

при i k ,

n

 

 

 

 

= 0 .

 

aij Akj =| A |;

 

 

 

aij Akj

 

 

j =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

Таким образом, получим

| Aj |

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

aij

 

 

=

 

 

 

 

bi | A |= bi , i =

1, n

 

 

 

 

 

 

 

 

 

 

| A |

| A

|

 

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

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

Пример. Решить систему линейных уравнений методом Крамера:

x1 + 2x2 + 3x3 = 7, x1 3x2 + 2x3 =5, x1 + x2 + x3 =3.

37

РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Решение. Вычислим определитель | A |,| A1 |,| A2 |,| A3 |:

| А|=

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 3 2

=9 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| А

 

 

7

2

3

 

 

 

 

 

 

 

1

0

1

 

 

 

 

 

1

1

 

= −(5 14) =9 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|=

 

5

3

2

 

=

 

 

 

14

0

5

 

= −

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

5

 

 

 

 

 

 

 

 

 

 

3

1

1

 

 

 

 

 

 

 

3

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| А

 

 

 

 

1 7

3

 

 

 

 

 

 

 

1

 

 

 

 

7

3

 

 

 

 

2

1

 

= 4 4 = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|=

 

 

 

1 5 2

 

=

 

0 2 1

=

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 2

 

 

 

 

 

 

 

 

 

 

 

 

1 3 1

 

 

 

 

 

0 4 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| А

 

 

1

2

7

 

 

 

 

1

2

7

 

 

 

5

 

 

2

 

= 20 2 =18 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|=

 

1 3 5

 

=

 

0 5 2

 

=

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 4

 

 

 

 

 

1

1 3

 

 

 

0 1 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда x1 =1, x2 = 0, x3 = 2.

Решение системы линейных уравнений с определителем |A|, отличным от нуля, можно найти с помощью обратной матрицы. Для этого запишем систему (4.2.1) в виде матричного уравнения

где A = (a

 

 

 

 

 

 

 

 

 

 

АХ=В

 

 

(4.2.5)

ij

)

n,n

; X = (x

, x

2

,..., x

n

)T , B = (b ,b ,...,b )T .

 

 

 

 

 

1

 

 

1

 

2

 

n

 

 

 

Решение матричного уравнения (4.2.5) имеет вид

 

 

 

 

 

 

 

 

 

 

X = A1B

 

 

(4.2.6)

Пример. Решить систему линейных уравнений с помощью обратной матрицы

 

 

 

 

 

 

 

 

 

x1 + 2x2 + 3x3 = 2

 

 

 

 

 

 

 

 

 

 

x1 3x2 + 2x3 = 0

 

 

 

 

 

 

 

 

 

 

x1 + x2 + x3 =1

 

 

Решение. Вычислим для матрицы

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

 

 

 

 

 

 

 

 

 

 

 

А= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

ее обратную матрицу

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

5

1

 

13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А1 =

 

1

2

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

4

1

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определим неизвестную матрицу-столбец Х:

 

 

5

1

13

 

2

 

 

 

3

1/ 3

 

Х = А1В =

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

,

 

1

2 1

 

 

0

 

=

 

3

= 1/ 3

9

9

 

 

4

1

5

 

 

1

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1/ 3

 

откуда x1 =1/ 3, x2 =1/ 3, x3 =1/ 3.

38

РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Формулы Крамера (4.2.4) могут быть получены из выражения (4.2.6). Действитель-

но, запишем матричное равенство

x1

xM2 =xn

X = A1B в развернутом виде:

 

 

 

A

A

...

A

b

 

 

 

 

11

21

 

n1

 

1

 

1

 

 

A12

A22

...

An2

b2

.

 

 

 

 

 

... ...

 

 

 

| A |

 

 

M

... ...

 

 

 

 

 

A

A

...

A

b

 

 

 

 

1n

2n

 

nn

 

n

 

Из полученного выражения непосредственно следуют формулы Крамера:

 

1

n

xi =

Akibk , i =

1, n

.

 

 

| A | k=1

4.3. Теорема Кронекера-Карелли

Теорема. Система линейных уравнений (4.1.1) совместна тогда и только тогда, ко-

гда = ~ . r( A) r( A)

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

Необходимость. Пусть система (4.1.1) совместна и пусть числа α1,α2 ,...,αn – одно из ее решений. Подставляя эти числа вместо неизвестных в систему (4.1.1), получим m тож-

деств, которые показывают, что последний столбец матрицы ~ является линейной комби-

A

нацией всех остальных столбцов, взятых соответственно с коэффициентами α1,α2 ,...,αn .

Всякий другой столбец матрицы ~ входит и в матрицу А. Поэтому максимальное число

A

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

~

~

A

совпадает. Следовательно, r( A) = r( A) .

Достаточность. Пусть дано, что = ~ = . Отсюда следует, что максимальное r(A) r( A) r

число линейно независимых столбцов матриц А и ~ совпадает и равно r. Для определен-

A

ности предположим, что первые r столбцов матриц А и ~ линейно независимы, а осталь-

A

ные (n-r) столбцов является их линейными комбинациями. Выражая последний столбец матрицы А как линейную комбинацию первых r столбцов, получим:

ai1 α1 + ai2 α2 +... + air αr = bi , i =1, m

или

ai1 α1 + ai2 α2 +... + air αr + ai,r +1 0 +... + ai,n 0 = bi , i =1, m

откуда следует, что числа α1,α2 ,...,αr ,0,...,0 являются решением системы (4.1.1), т.е. сис-

тема (4.1.1) совместна. Теорема доказана.

На основании теоремы Кронекера-Капелли имеем:

1.

~

Если r( A) r( A) , то система (4.1.1) несовместна;

2.

~

Если r( A) = r( A) = r , то система (4.1.1) совместна.

Пусть для определенности базисный минор порядка r расположен в верхнем левом углу матрицы А. Тогда первые r строк матрицы А линейно независимы, а остальные ее строки являются линейной комбинацией первых r строк. Но это означает, что первые r уравнений системы (4.1.1) линейно независимы, а остальные (m-r) ее уравнений являются их линейными комбинациями. Поэтому достаточно решить систему r уравнений; решения такой системы будут, очевидно, удовлетворять и остальным (m-r) уравнениям.

При этом возможны два случая:

39

РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

1. r = n . Тогда систему, состоящую из первых r уравнений системы (4.1.1)

a11x1 + a12 x2 +... + a1n xn = b1, a21x1 + a22 x2 +... + a2n xn = b2 ,

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

ar1x1 + ar 2 x2 +... + arn xn = br .

можно решить, например, по правилу Крамера. В этом случае система имеет единственное решение, т.е. система совместна и определена;

2. r < n . Рассмотрим первые r уравнений системы (4.1.1). Оставив в левых частях первые r неизвестных, перенесем остальные в правые части. Получим систему:

a11x1 +... + a1r xr =b1 a1,r+1xr+1 ... a1n xn , a21x1 +... + a2r xr =b2 a2,r+1xr+1 ... a2n xn ,

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

ar1x1 +... + arr xr =br ar,r+1xr+1 ... arn xn .

Очевидно, что полученная система и, следовательно, система (4.1.1) являются совместными и неопределенными.

Таким образом, если = ~ , то система (4.1.1) совместна (определенная или r( A) r( A)

неопределенная), если < ~ , то система (4.1.1) несовместна. r( A) r( A)

Если в системе n линейных уравнений с n неизвестными определитель системы ра-

вен нулю, то < . Тогда если = ~ , то система является совместной и неопре- r( A) n r( A) r( A)

деленной. Если < ~ , то система несовместна. r( A) r( A)

Теорема Кронекера-Капелли устанавливает необходимое и достаточное условие совместности системы (4.1.1), но не дает способа нахождения решения этой системы. Рассмотрим методЖордана-Гаусса– методрешениясистемыm линейныхуравненийсn неизвестными.

4.4. Метод Жордана-Гаусса

Метод Жордана-Гаусса основан на элементарных преобразованиях (п.3.2) строк расширенной матрицы

a11

~ = a21

A ...

am1

a

...

a

b

 

 

12

 

 

1n

1

 

a22

...

a2n

b2

 

...

... ...

...

 

 

a

m2

...

a

mn

b

 

 

 

 

m

системы (4.1.1).

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

Пусть дана система m линейных уравнений с n неизвестными. Применяя элементарные преобразования, построим эквивалентную систему специального вида. Для этого выберем в качестве первого уравнений одно из тех уравнений системы, где коэффициент при х1 отличен от нуля. Не нарушая общности, предположим, что а11 0 . Тогда первым

уравнением системы будет уравнение

a11 x1 + a12 x2 +... + a1n xn = b1 .

40