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

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

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

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

219

По теореме 6.5.1 (о базисном миноре) последние m r урав- нений являются линейными комбинациями первых r уравне- ний, и, следовательно, их можно отбросить, поскольку они бу- дут тождественно удовлетворяться решениями первых r урав- нений.

В оставшихся уравнениях перенесем в правые части слагаемые, содержащие неизвестные ξr +1 , ξr +2 ,..., ξn .

Неизвестные ξ1 ,..., ξr называются основными (главными, за-

висимыми, базисными), а неизвестные ξr +1 ,..., ξn свободны-

ми (параметрическими, независимыми, небазисными).

α11ξ1 + α12 ξ2 + ... + α1r ξr = β1 − α1r +1ξr +1 ...− α1n ξn ,

α 21ξ1 + α22 ξ2

+ ...

+ α 2r ξr = β2 − α2r +1ξr +1 ...

− α 2n ξn ,

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

 

 

 

 

α r1ξ1 + α r 2 ξ2

+ ...

+ α rr ξr = βr − αr ,r +1ξr +1 ...

− α r ,n ξn .

Присвоим свободным неизвестным некоторые конкретные зна- чения ξr +1 = μ1 ,..., ξn = μ nr и рассчитаем по правилу Кра-

мера (теорема 6.4.1) соответствующие им значения основных неизвестных:

j -й столбец

 

 

 

 

α11

 

nr

 

 

α1r

 

 

 

 

 

 

...

β1 α1,r +k μ k

...

 

(6.7.1)

ξ j

=

1

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

det

...

...

...

 

...

...

,

 

M

 

 

 

 

 

α r1

 

nr

μ k

 

α rr

 

 

 

 

 

 

...

βr α r ,r +k

...

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

где j = [1, r] , а M

базисный минор.

 

 

 

 

220

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

 

2°. Заметим, что из соотношений (6.7.1), положив

 

 

 

μk = 0; k = [1, n r] ,

 

 

можно найти частное решение неоднородной системы (6.6.1).

Теперь рассмотрим однородную систему. По линейному свой- ству определителей (теорема 6.2.3) получаем выражения для значений неизвестных:

 

 

 

 

nr

 

 

 

 

 

 

 

 

 

 

ξ j = κ jk μ k ,

j = [1, r] ;

 

 

 

 

k =1

 

 

(6.7.2)

 

 

 

ξ r +i = μi ,

i = [1, n r] ,

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

κ jk

=

1

 

α11 ...

− α1,r +k ...

α1r

 

 

 

 

 

 

 

det

... ...

...

...

...

 

 

 

,

M

 

 

 

α r1 ...

− α r ,r +k ...

α rr

 

 

 

 

 

 

 

 

 

 

 

 

j = [1, r] , k = [1, n r].

j -й столбец

Наконец, в матричной форме соотношения (6.7.2) могут быть записаны в виде

ξ1

 

κ11

κ12

L κ1,nr

 

 

ξr +1

 

 

 

 

 

 

 

 

ξ2

=

κ21

κ 22

L κ 2,nr

 

 

ξr +2

 

 

(6.7.3)

L

 

L

L L L

 

 

L

 

 

 

ξr

 

κr1

κ r 2

L κr ,nr

 

 

ξn

 

 

 

 

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

 

 

 

221

 

 

ξ1

 

 

κ11

κ12

L κ1,nr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ2

 

 

κ21

κ 22

L κ2,nr

 

 

 

μ1

 

 

 

 

 

 

 

L

 

 

L

L

L

L

 

 

 

 

 

или

ξr

 

=

κr1

κr 2

L κr ,nr

 

 

 

μ 2

.

 

ξr +1

 

1

0

L

0

 

 

 

L

 

 

ξr +2

 

 

0

1

L 0

 

 

 

μ nr

 

 

 

L

 

 

L

L

L

L

 

 

 

 

 

 

 

ξn

 

 

0

0

L

1

 

 

 

 

 

3°. Полагая μ1 = 1, μ2 = μ3 = ... = μnr = 0 , получим решение

{ξ11 , ξ12 , ... , ξ1r , 1, 0, ... , 0} .

Аналогично при μ1 = 0, μ2 = 1, μ3 = ... = μnr = 0 найдем решение {ξ12 , ξ22 , ... , ξ2r , 0 ,1, ... , 0}. И, продолжая этот про-

цесс, получим на последнем шаге при

μ1 = μ 2 = μ3 = ... = μ nr 1 = 0, μ nr = 1

решение {ξ1nr , ξn2r , ... , ξnr r , 0, 0, ... , 1} .

Полученные решения будем называть нормальными фунда-

ментальными решениями.

4°. Покажем теперь, что построенные n r частных решений однородной системы уравнений (6.6.1) являются линейно неза- висимыми. Действительно, записав эти решения как строки, получим матрицу вида

ξ11

ξ12

...

ξ1r

1

0 ...

0

 

ξ12

ξ22

...

ξr2

0

1 ...

0

. (6.7.4)

...

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

 

ξ1nr

ξ2nr

...

ξrnr

0

0 ...

1

 

222

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

 

Заметим, что ее ранг, с одной стороны, не меньше, чем

 

 

n r,

поскольку содержит ненулевой минор этого поряд-

ка, но, с другой стороны, не больше, чем число строк в этой матрице, равное n r, и потому ранг в точности равен

n r, что доказывает линейную независимость построен-

ных частных решений.

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

Определение Фундаментальной системой решений для системы

6.7.1.линейных уравнений (6.6.1) называется совокупность

любых n − rg A частных, линейно независимых решений однородной системы (6.6.1), где n число неизвестных в системе (6.6.1), а A ее основная

матрица. Матрица (6.7.4) называется фундаменталь-

ной.

Теорема Каждое частное решение однородной системы (6.6.1)

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

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

Пусть дано решение 1 , ξ2 ,..., ξn } однородной системы

(6.6.1). Рассмотрим матрицу размера (n r + 1) × n

ξ1

ξ2

...

ξr

ξr +1

ξr +2

...

ξn

 

 

 

 

 

ξ11

ξ12

...

ξ1r

1

0

...

0

 

 

 

ξ12

ξ22

...

ξr2

0

1

...

0

 

 

, (6.7.5)

...

... ... ...

...

... ... ...

 

 

 

ξ1nr

ξ2nr

...

ξrnr

0

0

...

1

 

 

 

 

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

223

 

ранг которой, с одной стороны, очевидно,

не меньше, чем

 

 

n r . С другой стороны, первые r столбцов этой матрицы

являются линейными комбинациями (заданными соотноше- ниями (6.7.3)) последних n r столбцов. Действительно, эти соотношения, связывающие значения свободных и основных переменных, одни и те же для всех строк матрицы (6.7.5), и по- тому в этой матрице каждый из первых r столбцов есть ли- нейная комбинация последних n r . Значит, ранг матрицы не превосходит n r , и, следовательно, равен в точности n r .

Наконец, по теореме 6.5.1 − о базисном миноре, который рас- полагается в последних r строках, первая строка матрицы (6.7.5) должна быть некоторой линейной комбинацией осталь- ных, и, следовательно, общее решение однородной системы (6.6.1) может быть записано в виде

 

ξ

1

 

 

 

ξ1

 

 

ξ 2

 

 

ξ nr

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

1

 

 

1

 

 

1

 

 

2

 

 

 

ξ1

 

 

ξ 2

 

 

ξ nr

 

 

 

 

 

 

2

 

 

2

 

 

2

 

 

...

 

 

 

...

 

 

...

 

 

...

 

 

ξ

r

 

= λ

1

ξ1

+ λ

2

ξ 2

+

... + λ nr

ξ nr

,

 

 

 

r

r

r

 

ξ r +1

 

 

 

1

 

 

0

 

 

0

 

 

ξ r +2

 

 

 

0

 

 

1

 

 

0

 

 

...

 

 

 

...

 

 

...

 

 

...

 

 

ξ n

 

 

 

0

 

 

0

 

 

1

 

где λi

i = [1, n r] произвольные константы.

 

 

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

 

 

 

 

 

 

 

 

Следствие

 

Общее решение неоднородной системы (6.6.1) может

6.7.1.быть дано формулой

224

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

ξ1

ξ2

...

ξr

ξr +1

ξr + 2

...

ξ n

 

ξ11

 

 

ξ12

 

 

 

 

 

 

 

ξ12

 

 

ξ 22

 

 

 

 

 

 

 

...

 

 

...

 

 

 

 

 

 

= λ1

ξ1r

+ λ

 

ξ r2

+ ...

+ λ n r

1

2

0

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

 

 

 

 

 

...

 

 

...

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

ξ10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ02

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

где

ξ0r

 

 

 

 

 

 

 

 

 

 

ξ0r +1

 

 

 

 

 

 

 

 

 

 

ξ0r +2

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

ξ0n

 

 

 

 

является некоторым частным родной системы (6.6.1), а числа произвольные константы.

ξ1n r

 

ξ10

 

ξ 2n r

 

ξ 02

 

...

 

...

 

ξ rn r

+

ξ 0r

,

0

ξ 0r +1

 

 

0

 

ξ 0r + 2

 

...

 

...

 

1

 

ξ 0n

 

решением неодно-

λi i = [1, n r]

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

 

 

x0

 

 

 

 

 

некоторое (найденное,

 

например, подбором) ча-

 

 

 

 

 

 

 

 

 

стное решение неоднородной системы (6.6.1), а

 

 

x

 

ее про-

 

 

 

извольное решение. Тогда по лемме 6.6.3 произвольное реше-

ние

 

однородной

системы (6.6.1)

 

y

 

 

 

 

представимо в виде

 

 

 

 

 

 

 

y

 

 

 

=

 

 

 

x

 

 

 

 

 

 

x0

 

. Откуда получаем

 

 

 

x

 

 

 

=

 

 

 

y

 

 

 

+

 

 

 

 

x0

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

225

Из теорем 6.7.1 и 6.7.2 непосредственно вытекает

 

Следствие Для того чтобы однородная система (6.6.1) с n < m

6.7.2.имела ненулевое решение, необходимо и достаточно,

чтобы ранг ее основной матрицы удовлетворял ус- ловию rg A < n .

В случае, когда основная матрица однородной сис- темы (6.6.1) квадратная, условие существования нетривиального решения равносильно равенству

det A = 0 .

Иное полезное для приложений условие совместности системы линейных уравнений дает

Теорема

Для того чтобы система (6.6.1) была совместной,

6.7.3необходимо и достаточно, чтобы каждое реше-

(Фредгольма).

 

 

 

y

 

=

 

η1

η2

 

 

...

 

ηm

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ние

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сопряженной

системы

 

α

 

 

η + α

 

 

η

 

 

+ ... + α

 

 

η

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

21

2

 

m1

m

 

 

 

 

 

 

 

 

 

 

 

 

11

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

η + α

 

 

η

 

 

+ ... + α

 

 

η

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

12

 

1

 

22

 

 

2

 

 

 

 

 

 

 

 

 

 

m 2

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

α

1n

η + α

2n

η

2

+ ... + α

mn

η

m

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(или в матричном виде

 

A

 

 

 

T

 

 

 

y

 

 

=

 

 

 

 

o

 

 

 

) удов-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

летворяло условию

 

βi ηi = 0 (или в матрич-

i=1

ном виде b T y = 0 ).

Доказательство необходимости.

Пусть система уравнений (6.6.1) совместна, то есть для каждо-

го ее решения x справедливо равенство b = A x .

 

226

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

Найдем

 

 

 

 

 

произведение

 

 

 

b

 

 

 

T

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

в предположении, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

T

 

 

 

y

 

 

 

=

 

 

 

 

 

o

 

 

 

 

 

 

 

. Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

T

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

= (

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

)T

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

=

 

 

x

 

 

 

T

 

 

 

 

 

A

 

 

 

T

 

 

 

y

 

 

 

=

 

 

 

x

 

 

 

T

 

 

 

o

 

 

 

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство достаточности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

b

 

 

 

T

 

 

 

 

y

 

 

 

= 0

 

 

 

 

 

 

 

 

для

 

 

любого

решения системы линейных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнений

 

 

 

A

 

 

 

T

 

 

 

 

 

y

 

 

 

 

=

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

 

. Тогда общие решения систем ли-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нейных уравнений

 

A

 

 

 

T

 

 

 

 

 

 

y

 

 

=

 

 

 

 

 

 

 

o

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

T

 

 

 

 

 

 

 

y

 

 

 

 

 

=

 

 

 

 

 

 

o

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

T

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

совпадают, и для этих систем максимальное число линейно неза- висимых решений одинаково. Поэтому, согласно теоремам 6.7.1

и 6.7.2,

m − rg

 

 

 

A

 

 

 

T = m − rg

 

A

 

 

 

T

или rg

 

 

 

A

 

 

 

T = rg

 

A

 

 

 

 

T ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

но поскольку ранг матрицы не меняется при ее транспонирова-

нии, то имеет место равенство rg A = rg A| b , означающее в силу теоремы 6.6.1 совместность системы линейных уравнений

A x = b .

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

Альтернативное доказательство теоремы Фредгольма приведено в главе 10 (см. теоремы 10.6.4 и 10.6.5).

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

227

§ 6.8. Элементарные преобразования. Метод Гаусса

 

Практическое применение теорем 6.7.3 и 6.7.4 затрудняется тем, что заранее, как правило, неизвестно, совместна ли решаемая система. Определение же рангов основной и расширенной матриц независимо от поиска решений оказывается весьма нерациональной (с точки зре- ния расходования вычислительных ресурсов) процедурой.

Более эффективным вычислительным алгоритмом, позволяющим либо находить общее решение системы (6.6.1), либо устанавливать факт ее несовместности, является метод Гаусса.

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

Под наиболее простымвидом расширенной матрицы мы будем понимать верхнюю треугольную форму (т.е. случай, когда α ij = 0

при i > j ), для которой возможно рекуррентное нахождение неиз-

вестных путем лишь решения на каждом шаге процедуры линейного уравнения с одним неизвестным. Ниже приведен пример матрицы размера m × n (n > m) , имеющей верхнюю треугольную форму

a11

a12

a13 ...

a1,m2

a1,m1

a1,m

a1,m+1

...

a1n

 

0

a22

a23 ...

a2,m2

a2,m1

a2,m

a2,m+1

...

a2n

 

0

0

a33 ...

a3,m2

a3,m1

a3,m

a3,m +1

...

a3n

.

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

...

...

...

...

...

...

 

0

0

0 ...

0

am 1,m1

am1,m

am1,m+1

... am1,n

 

0

0

0 ...

0

0

amm

am,m+1

...

amn

 

Кэлементарным преобразованиям матрицы относятся:

-перестановка строк (перенумерация уравнений);

-перестановка столбцов основной матрицы (перенуме- рация неизвестных);

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

-удаление нулевой строки (исключение уравнений, тож- дественно удовлетворяющихся любыми значениями неизвестных);

-умножение строки на ненулевое число (нормирование уравнений);

-сложение строки с линейной комбинацией остальных строк с записью результата на место исходной строки (замена одного из уравнений системы следствием ее

уравнений, получаемым при помощи линейных опера- ций).

Решение неоднородной системы уравнений (равно как и ранг ее матрицы) не изменится также и при использовании любой комбина- ции элементарных операций.

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

-

перестановка столбцов с номерами i и j матрицы

 

 

A

 

 

 

размера

m × n осуществляется путем ее умно-

 

 

 

 

 

жения справа на матрицу

 

 

 

S

 

 

 

1 размера n × n , которая

 

 

 

 

 

 

 

 

 

 

 

в свою очередь получается из единичной матрицы n -

 

го порядка

 

 

 

E

 

 

 

путем перестановки в последней i -го

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и j -го столбцов;

-умножение i -й строки матрицы A на некоторое число l ¹ 0 осуществляется путем умножения A

слева на матрицу S 2 , которая получается из единич-

ной размера m × m матрицы E путем замены в по-

следней i -го диагонального элемента (равного едини-

це) на l ;