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

RUDN-I

.pdf
Скачиваний:
11
Добавлен:
25.03.2016
Размер:
721.58 Кб
Скачать

1.3. Системы линейных уравнений. Основные понятия

21

Переменные xm+1, . . . , xn перенесем в правую часть и получим:

.

 

 

.

e

.

 

..

 

e .

 

 

 

x1

=

b1

a1(m+1)

xm+1

. . .

a1n

xn

 

x2

=

b2

a2(m+1)

· xm+1

. . .

a2n · xn

 

 

 

 

 

 

e

 

 

 

e

 

 

 

 

 

e ru

 

 

.

 

 

e

.

·

 

e

 

(2)

 

 

 

e.

 

 

.·

 

.

 

 

.

 

 

 

 

.

 

 

 

.

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

am(m+1) · xm+1

− . . . −

amn · xn.

 

xm

= bm

 

Переменные

 

m+1

 

e

 

n

 

 

 

 

 

 

 

.

 

 

x

 

 

,

. . . ,

x

 

,

 

 

 

 

 

 

|

 

 

 

 

 

 

0 matematem1 . . . 0 a2(r+1) . . . a2n

 

любые значения. Каждому набору свободных переменных отвечает од-

нозначно набор базисных переменных x1,

. . . ,

xm (по формулам (2)), а

значит, и решение системы (1).

 

 

 

 

 

 

 

 

Итак, в этом случае система имеет бесконечно много решений.

Замечание 4. В этом случае система (1) при любом столбце B правых частей разрешима и имеет бесконечно много решений. В частности, однородная система имеет не только нулевое решение.

Случай 4 (r < m, r < n).

 

 

 

 

 

 

 

 

 

 

1

0 . . .

0

||

a1(r+1)

. . . a1n

 

 

 

 

0

1 . . .

0

a2(r+1)

. . .

a2n

 

 

 

 

 

 

 

 

 

e

 

 

e

 

 

 

 

 

.. .. ..

..

 

e ..

..

.

e..

 

 

основная матрица

 

. .

. .

| .

 

.

 

A =

 

0

0 . . .

1

|

ar(r+1)

. . .

arn

 

преобразованной

 

 

 

 

 

 

 

 

 

 

 

— — —

 

|

 

 

 

— —

 

 

системы;

 

 

0

0 . . .

0

|

e 0

. . .

e0

 

 

 

 

 

 

 

e

 

 

 

 

 

.. .. .. ..

|

..

..

 

..

 

 

 

 

 

 

 

 

 

.

 

.

.

 

 

 

 

. . . .

 

 

 

 

 

 

 

0

0 . . .

0

 

0

. . .

0

 

 

 

 

 

 

 

1

0 . . .

0

a1(r+1)

. . .

a1n

 

b1

 

 

 

 

 

 

. .

 

 

. .

 

 

 

.

|.

b

 

 

 

 

 

 

 

.

 

.

 

 

e2

 

 

 

 

 

 

.

 

 

 

e

 

e

 

e

 

 

 

 

 

.. ..

 

.. .. ..

 

.. .. ..

 

 

 

 

 

 

 

 

 

 

e

 

 

e

 

 

 

 

 

(A

 

B) =

0

0

. . .

1

a

r(r+1)

. . . a

rn

 

b

r

 

 

 

 

 

 

 

 

 

|

 

 

 

|

 

 

— — — — —

 

 

 

 

 

 

 

 

 

 

— — — —

 

 

e

 

e

 

0 0 . . .

0

e 0

. . .

e0

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. .

 

 

. .

 

 

 

.

|.

 

 

 

 

 

 

 

. . ... . .

... . .

e

 

 

 

 

 

 

. .

 

 

. .

 

 

 

. .

 

 

 

 

 

 

 

0 0 . . .

0

 

0

. . .

0

 

bm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

e

 

 

 

— расширенная матрица преобразованной системы.

 

 

 

 

а) Еслиwwwb(r+1) = . . . = bm = 0, то последние (m − r) уравнений окажутся

e e

тождествами 0 = 0 и их можно отбросить. Далее, как в случае 3 (при

22

Тема 1.

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

m = r) переменные xr+1,

. . . , xn — свободные, могут принимать любые

значения; а базисные переменные x1, . . . , xr выражаются через них по формулам

 

 

x1 = b1

a1(r+1)

xr+1 . . .

a1n

xn

 

 

 

x2

= b2

a2(r+1)

· xr+1

. . .

a2n

· xn

 

решений

 

 

.

e.

.

·

− −

 

.·

 

 

.

.

 

e .

 

..

.

e .

 

бесконечно много.

.

.

 

.

 

 

 

.

 

 

 

e

 

e

 

 

 

e

 

 

 

ru

 

 

 

 

 

 

 

 

.

xr

= br

ar(r+1) · xr+1

− . . . −

arn · xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частности, однородная система имеет не только нулевые решения (для

В

 

 

e

 

e

 

 

 

e

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

нее B = B =

0

 

 

 

 

 

 

 

 

...

).

 

 

 

 

 

 

 

 

 

 

 

e

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вывод 2. При m = n имеетmatematemместо следующая альтернатива (альтер-

б) Если среди чисел br+1, . . . ,

bm хотя бы одно не равно 0, то система

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

e

e

Замечание 5. Из приведенных выше ответов следует, что однородная система (1) в случаях 1 и 2 имеет единственное решение. Это нулевое решение x1 = x2 = . . . = xn = 0. В случаях 3 и 4 решений бесконечно много.

В заключение сформулируем некоторые выводы, важные для дальнейшего.

Вывод 1. Если число неизвестных в однородной системе (1) больше числа уравнений. т.е. n > m, то реализуются случаи 3 или 4 и, значит, заведомо существуют. ненулевые решения однородной системы.

натива Фредгольма):

1)либо однородная система имеет только нулевое решение, а соответствующая неоднородная система разрешима (причем единственным образом) при любом столбце правых частей B;

2)либо существуют ненулевые решения однородной системы, а соответ-

ствующая неоднородная система совместна не при любом столбце правых частей B.www

Действительно, при m = n может возникнуть либо случай 1, и тогда реализуется первая возможность (см. замечание 1), либо случай 4, и тогда реализуется вторая возможность.

1.3. Системы линейных уравнений. Основные понятия

23

Замечание 6. Отметим, что при m ≠ n альтернатива Фредгольма не имеет места (см. замечания 2 и 3).

Пример. Найти все решения системы

 

 

 

ru

 

 

 

 

 

 

 

 

 

 

3x1 + 2x2 + x3 + 2x4 = 0

 

 

 

 

 

 

 

 

 

2x1

 

 

− x3 + x4 = 0

 

 

 

 

 

 

 

 

 

 

 

 

4x2

+ 5x3 + x4 = 0

 

 

 

Прямой ход метода Гаусса:

 

 

 

 

 

 

.

 

 

 

 

 

 

A =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

1

2

3

2

1

 

2

 

 

 

 

2 0

1 1

0

4/3

5/3

1/3

 

 

 

 

 

 

0 4

5 1

0

4

5

1

 

 

3

2

1

2

 

 

 

3

2

||

 

1

2

= A — ступенчатая матрица.

 

0

4

 

5

1

0

4

5

1

 

 

0 4 5 1

 

— — — — —

 

e

 

 

 

 

 

 

 

0 0

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

 

 

Обратный ход метода Гаусса:

 

 

 

 

 

 

 

 

 

 

 

 

A =

3 2

||

1 2

 

3 0

 

3/2 3/2

 

 

 

 

0 4

5 1

0 4 || −

5

 

1

 

 

 

e

 

 

— — — — —

 

 

— — —

 

 

 

 

 

0 0

 

0 0

0 0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

1

 

0

||

1/2 1/2

 

 

 

 

 

 

 

 

 

 

 

0

 

1

5/4

1/4

= A — простейшая матрица.

 

 

 

— —

 

 

 

 

b

 

 

 

 

 

 

 

 

0

 

0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

matematem

 

 

 

 

 

Преобразованная система имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

1/2x3 + 1/2x4 = 0

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

{

 

x2 + 5/4x3 + 1/4x4 = 0,

 

 

 

 

 

 

 

 

 

 

x1

=

 

1/2x3

1/2x4 — ее решения.

 

 

(3)

 

 

 

 

{ x2 = 5/4x3 1/4x4

 

 

 

 

 

 

При этомwwwсвободные переменные x3, x4 могут принимать любые значения. Каждому набору свободных переменных x3, x4 отвечает по формуле

(3) определенный набор базисных переменных x1, x2, т.е. определенное решение системы.

24 Тема 1. Системы линейных уравнений. Метод Гаусса

Замечание 7. Изложенный метод решения системы (1) допускает элементарные преобразования 2-го рода (перестановку столбцов матрицы A при приведении матрицы A к ступенчатому виду). Для практического решения систем это не вполне удобно, т.к. связано с перенумерацией неизвестных. Обычно, необходимость в перестановке столбцов матрицы A не возникает. Исключительная ситуация, которая может потребовать перестановки столбцов, возникает, если на некотором шаге k в

ми в k-ой строке), 1 ≤ j1 < j2 < . . . jr ≤ n (причем номера могут идти не подряд). Обратный ход метода Гаусса дает тогда простейшую матрицу,

правом нижнем углу появится ненулевая матрица Ak, в которой один

 

 

 

 

 

 

 

 

 

 

 

ru

 

 

или несколько первых столбцов — нулевые. Если при этом избегать пе-

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

рестановки столбцов, то мы придем в итоге к ступенчатой матрице A

более общего вида, чем описан выше:

 

 

 

 

 

 

 

 

e

 

. . .

| a1j1 . . .

 

.a. 2.j2 ...... .. .. ..

 

.. .. ..

 

 

 

 

 

 

 

A ≈ A =

 

 

 

O

|

 

 

...

 

...

 

,

(4)

e

 

 

 

 

 

 

 

 

arjr

 

. . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

в которой элементы λk = akjk

̸= 0, k = 1,

2, . . . , r (их называют ведущи-

 

 

matematem

 

 

 

 

 

в которой akjk = 1, k = 1, 2, . . . , r, а элементы над ними равны нулю. В этом случае переменные xj1 , xj2 , . . . , xjr являются базисными, а

остальные свободными. Перенося свободные переменные в правые части уравнений, получим систему, эквивалентную (1), в которой из r уравнений базисные переменные xj1 , xj2 , . . . , xjr однозначно выражаются через свободные, а остальные уравнения (если r ≤ m − 1) имеют вид

e

0 · xj1 + 0 · xj2 + . . . + 0 · xjr = bk, k = r + 1, . . . , m.

3) существованиеwww"нулевого элемента" , т.е. матрицы O, такой что A+O = O + A = A для любой матрицы A,

Если все bk = 0, то получим тождества 0 = 0. Их можно отбросить

и мы получим бесконечно много решений системы, отвечающих любым

наборам свободныхe

.переменных. Если же среди чисел bk есть не равные

нулю, то система несовместна.

e

1.4. Теоретические упражнения к теме 1

1.1. Показать, что сложение матриц обладает следующими свойствами:

1) коммутативность: A + B = B + A,

2) ассоциативность: A + (B + C) = (A + B) + C,

1.4. Теоретические упражнения к теме 1

25

4) существование у каждой матрицы A ей противоположной −A, такой что A + (−A) = (−A) + A = O.

1.2. Показать, что умножение матрицы на число обладает следующими

свойствами:

 

ru

1)

1 · A = A,

.

2) α(βA) = (αβ)A.

1.3. Проверить свойства дистрибутивности:

1) α(A + B) = αA + αB,

2) (α + β)A = αA + βA.

 

1.4.Определим вычитание матриц: A − B := A + (−B). Показать, что 1) α(A − B) = αA − αB,

2) (α − β)A = αA − βA.

1.5.Показать, что транспонирование матрицы обладает следующими

свойствами:

1)(A + B)T = AT + BT ,

2)(αA)T = α · AT ,

3)(AT )T = A.

1.6.Матрица A называется эквивалентной матрице B, если она получена из нее с помощью конечного числа элементарных преобразований. Показать, что это отношение эквивалентности обладает свойствами:

1) рефлексивности A ≈ A,

2) симметричности A ≈ B B ≈ A,

3) транзитивности A ≈ B, B ≈ C A ≈ C.

1.7.1) Показать, что matematem

 

 

 

2

 

 

 

.

1

x1

x12

 

1

0

0

 

 

x3

 

0

0

1

1 x3

 

1

x2

x22

 

0

1

0

при условии,wwwчто числа x1, x2, x3 разные. Рассмотреть также случаи совпадения двух из них; совпадения всех трех чисел.

2) Показать, что для любых различных чисел x1, x2, x3 и любых чисел y1, y2, y3 существует, причем единственный многочлен p(x) степени не выше 2, такой что p(xi) = yi, i = 1, 2, 3. Каков геометрический смысл этого утверждения?

1.8. Доказать следующие утверждения о связи решений однородной и неоднородной системы линейных уравнений.

1) Сумма решений однородной и неоднородной системы является решением неоднородной системы.

26 Тема 1. Системы линейных уравнений. Метод Гаусса

2) Разность двух решений неоднородной системы является решением однородной системы.

3) Пусть известно некоторое частное решение неоднородной системы.

 

 

 

 

ru

Тогда все ее решения можно получить, добавляя к нему решения одно-

родной системы.

 

 

.

 

 

 

matematem

 

www

.

 

 

 

 

 

 

 

 

 

27

Тема 2. Множества и отображения

Здесь рассматриваются некоторые общие понятия теории множеств

(сумма, произведение, разность, симметрическая разность, бинарные от-

ношения и, в частности, отношения эквивалентности). Обсуждаетсяru

од-

 

 

 

 

matematem

.

 

отображе-

 

 

 

 

 

 

 

 

 

 

 

 

(инъекции,

 

 

 

 

 

 

общим призна-

 

 

 

 

 

 

образующие

 

их

 

 

 

 

 

 

 

говорят, x

 

 

 

 

 

принадлежит

X.

 

 

1)

множество студентов данной группы A,

 

 

 

 

2)

множество студентов данного факультета B,

 

 

 

3)

множество студентов данного ВУЗа C

 

 

 

— это конечные множества;

 

 

 

 

4)

N = {1, 2, 3, . . .} — множество натуральных чисел,

 

 

5)

Z = {0, ±1, ±2, . . .} — множество целых чисел,

 

 

6)

 

 

.

 

 

 

 

 

R — множество действительных чисел

 

 

 

 

www

 

 

 

 

 

 

— это бесконечные множества.

 

 

 

 

 

 

 

 

 

 

если

a A

 

 

 

 

 

 

C, N

возможэлементы и

2) (a, b) = {x R : a < x < b} — интервал и т.д.

28 Тема 2. Множества и отображения

Пустое множество, множество без элементов . Например,

{x R : 2 < x < 1} = .

Равенство множеств: A = B, если A B и B A, т.е. A и B состоят

из одних и тех же элементов.

 

Определение 1. Сумма (объединение) C = A B определяется как

 

matematem

 

множество элементов, принадлежащих хотя бы одномуruиз множеств

A или B.

 

.

Определение 2. Пересечение (произведение) C = A ∩ B определяется как множество элементов, принадлежащих и A, и B.

Определение 3. Разность C = A \ B определяется как множество элементов из A, не принадлежащих B.

Определение 4. Симметрическая разность задается равенством

 

 

A B = (A \ B) (B \ A).

(1)

Иллюстрации:

 

 

 

A

B A

B A

B A

B

 

A B

A ∩ B

A \ B

A B

 

 

 

 

Рис. 1

Аналогично определяем сумму и пересечение любого (не обязательно конечного) семейства множеств:

(A B)www\ (B ∩ A). Итак, x A B x (A B) \ (B ∩ A), т.е.

C1

=

Aα

 

множество элементов, принадлежащих хотя бы

 

α

 

одному.

из Aα,

C2

= α Aα

множество элементов, принадлежащих всем Aα.

Ясно, что C2 C1.

Пример 3. Показать, что A B = (A B) \ (B ∩ A).

(1)

Пусть x A B = (A \ B) (B \ A). Тогда x A \ B, либо x B \ A. Если x A \ B, то x A, x ̸B x A B, x ̸A ∩ B x

A B (A B) \ (B ∩ A).

2.2. Отношения эквивалентности. Фактор-множества

29

Обратно, пусть x (A B) \ (B ∩ A), т.е. x A B, x ̸B ∩ A. Либо x A, x ̸A ∩ B, т.е. x A \ B, либо x B, x ̸A ∩ B, т.е. x B \ A.

Итак, либо x A \ B, либо x B \ A, т.е. x (A \ B) (B \ A) = A B.

2)

A ∩ B = B ∩ A

}

 

B.

 

ru

 

Значит, (A B) \ (B ∩ A) A

 

 

 

Следовательно, A B = (A B) \ (B ∩ A).

 

 

Упражнение 1. Показать свойства сложения и умножения

 

1)

A B = B A

 

— коммутативность,

 

 

 

α

 

matematemα α α

 

 

3)

(A B) C = A (B C)

— ассоциативность,.

 

 

4) (A ∩ B) ∩ C = A ∩ (B ∩ C)

}

 

 

 

5)

(A B) ∩ C = (A ∩ C) (B ∩ C)

— дистрибутивность.

 

6) (A ∩ B) C = (A C) (B C)

}

 

 

(см. упражнение 2.1).

 

 

 

 

 

Принцип двойственности. Для любых множеств A1, A2, B справед-

ливы равенства

 

 

 

 

 

 

 

 

B \ (A1 A2) = (B \ A1) (B \ A2)

 

(2)

Докажем (2):

B \ (A1 ∩ A2) = (B \ A1) (B \ A2).

 

(3)

 

 

 

 

 

 

x B \ (A1 A2) x B, x ̸A1 A2 x B, x ̸A1, x ̸A2

 

x B \ A1, x B \ A2 x (B \ A1) (B \ A2).

 

Значит, верно (2).

 

 

 

 

 

Упражнение 2.

 

 

 

 

 

1)

Доказать равенство (3).

 

 

 

2)

Доказать обобщенный принцип двойственности:

 

 

 

.

 

 

 

 

 

 

B \ Aα = (B \ Aα), B \ Aα = (B \ Aα)

 

 

wwwC = Aα

 

{C =

Aα, Aα

Aβ = (α ̸= β)}.

(5)

(сумма и пересечение могут быть и бесконечного набора множеств Aα) (см. упражнение 2.2).

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

C = A B

{C = A B, A ∩ B = };

(4)

аналогично

 

 

 

 

α

α

30

Тема 2. Множества и отображения

2.2.Отношения эквивалентности. Фактор-множества

Определение 1. Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ru

X1, X2 — непустые множества. Произве-

дением (декартовым произведением) X1 × X2 называется множество

всех

упорядоченных пар (x1, x2), где x1

 

X1, x2

 

X2. Произведение

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X × X ≡ X

 

называют декартовым квадратом множества X.

 

 

 

matematem

 

 

Пример 1.

 

 

1

 

 

 

 

 

2} := R

2 .

 

1) R × R = {(x1, x2), xi

R

, i = 1,

 

— множество точек

плоскости.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

6

 

 

 

 

 

 

 

 

q(x1, x2)

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

R1

 

 

Рис. 1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

2) P = [a, b)×(c, d] = {(x1, x2) : a ≤ x < b, c < x ≤ d} — прямоугольник на плоскости.

 

 

R1 6

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

www

0 a

 

b

R1

Рис. 2

 

Br = {(x1, x2) : x12 + x22 ≤ r2} — цилиндр с кругом

3)

C = Br × [0, 1], где.

Br в основании (см. рис. 3).

 

 

 

 

 

 

 

 

 

 

 

 

6x3

 

 

 

 

 

 

 

 

q1

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

0

 

x2

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3

 

 

x1

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]