
Лекция 14
.pdf
Лекция 14
6.2. Решение произвольной системы методом Гаусса
Рассмотрим произвольную систему линейных уравнений AX B и ее рас-
ширенную матрицу A |
|
|
B , т.е. матрицу, составленную из матрицы A и столбца |
||||||
|
|||||||||
свободных членов B : |
|
|
|
|
|
|
|
||
|
|
a 11 |
|
a 1n |
|
b 1 |
|
||
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
A |
|
B |
|
|
|
|
. |
||
|
|||||||||
|
|
|
|
a m n |
|
|
|
||
|
|
a m1 |
|
b m |
Если уравнение системы умножить на число, сложить с другим уравнением,
вычеркнуть одно из пропорциональных уравнений, поменять местами уравнения или неизвестные, то получим систему, равносильную исходной. Удобнее произ-
водить эти операции не с самими уравнениями, а с их коэффициентами, то есть с элементами расширенной матрицы.
Метод Гаусса (метод последовательного исключения неизвестных)
Прямой ход метода Гаусса состоит в приведении расширенной матрицы системы к матрице ступенчатого вида (с нулями ниже главной диагонали) с по-
мощью элементарных преобразований.
Обратный ход состоит в получении нулей выше главной диагонали. При этом на главной диагонали нужно получить ненулевые элементы (за исключением нулевых строк).
Будем использовать следующие элементарные преобразования:
умножения строки на число, отличное от нуля;
прибавления к строке другой строки;
вычеркивания одной из пропорциональных строк;
перестановки строк или столбцов местами;
вычеркивание нулевых строк.
Опишем подробнее прямой и обратный ход метода Гаусса.
1
Прямой ход. С помощью элементарных преобразований строк получим ну-
ли ниже главной диагонали сначала в первом столбце с помощью первой стро-
ки (это равносильно исключению x1 из всех уравнений, кроме первого); затем
получим нули ниже главной диагонали во втором столбце с помощью второй строки и т.д.
Обратный ход. Получим нули выше главной диагонали сначала в послед-
нем столбце с помощью последней строки, потом в предпоследнем столбце с
помощью предпоследней строки и т.д.. |
|
|
|
|
|
||||
|
Рассмотрим различные ситуации, которые могут возникнуть при примене- |
||||||||
нии метода Гаусса. |
|
|
|
|
|
|
|
||
|
Если в процессе элементарных преобразований мы получим строку из ну- |
||||||||
лей, |
а |
свободный |
член |
b 0 , |
то |
соответствующее |
уравнение |
||
0 x1 0 x2 |
0 xn b |
не имеет решение. Система несовместна. Заметим, что |
|||||||
за счет описанной строки ранг расширенной матрицы A |
|
B больше ранга мат- |
|||||||
|
рицы А.
Если в процессе элементарных преобразований мы получим строку расши-
ренной матрицы, всю состоящую из нулей, то эту строку вычеркнем (соответст-
вующее уравнение 0 x1 0 x2 0 xn 0 выполняется тождественно) и про-
должим метод элементарных преобразований.
Если получим ненулевую строку, но на главной диагонали будет нулевой элемент, то поменяем столбцы в матрице (а значит, неизвестные в уравнениях)
местами так, чтобы на главной диагонали оказался ненулевой элемент.
В итоге получим в расширенной матрице системы нули ниже и выше глав-
ной диагонали. При этом матрица системы приведется
|
|
|
1 |
0 |
|
0 |
d 1 |
|
|
|
|
|
|
|
|
|
|
|
|
либо к треугольному виду |
A |
|
B |
0 |
1 |
|
0 |
d 2 |
, |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
1 |
d n |
|
2
1 |
0 |
0 |
c 1r 1 |
c 1n |
d1 |
|
||||
|
0 |
1 |
0 |
c |
2 r 1 |
c |
2 n |
d |
2 |
|
либо к ступенчатому виду A | B |
|
|
|
|
|
|
. |
|||
|
|
|
|
|
|
|
|
|
||
|
0 |
0 |
1 |
cr r 1 |
cr n |
dr |
|
|||
|
|
В первом случае в результате элементарных преобразований получим рав-
носильную систему:
1 x1 0 x2 0 xn d1, |
|
|
x1 d1, |
|
|
|
|
|
|
|
||||||||||||||
|
|
1 x2 |
0 xn |
d2 , |
|
|
x2 d2 , |
|
|
|
|
|
|
|
||||||||||
0 x1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
............................................. |
|
|
|
|
|
|
|
|
........... |
|
|
|
|
|
|
|
|
|
|
|
|
|||
0 x 0 x |
2 |
1 x |
n |
d |
n |
, |
|
|
x |
n |
d |
n |
. |
|
|
|
|
|
|
|
||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Система имеет единственное решение, а rang A |
|
B rang A n . |
||||||||||||||||||||||
|
||||||||||||||||||||||||
Во втором случае rang A r n и в результате элементарных преобразова- |
||||||||||||||||||||||||
ний имеем равносильную систему: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1 x1 0 x2 0 xr c1 r 1xr 1 c1 n xn d1, |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 x |
1 x |
2 |
0 x |
r |
c |
2 r 1 |
x |
r 1 |
c |
2 n |
x |
n |
d |
2 |
, или |
|||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
................................................................................. |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 x2 |
1 xr |
cr r 1xr 1 cr n xn dr , |
|
|||||||||||||||||||
0 x1 |
|
|||||||||||||||||||||||
x1 d1 c1 r 1xr 1 c1 n xn , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2 c2 r 1xr 1 c2 n xn , |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
............................................... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dr cr r 1xr 1 |
cr n xn . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
xr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Получили общее решение системы уравнений; неизвестные x1, x 2 , , x r на- |
||||||||||||||||||||||||
зываются |
базисными неизвестными |
(число базисных неизвестных равно |
||||||||||||||||||||||
r rg A ), |
а неизвестные |
|
xr 1, , xn |
называются |
свободными неизвестными |
(их число равно n r ). Придавая свободным неизвестным числовые значения, по-
лучим из общего решения частные решения системы. Система совместна и имеет бесконечное множество решений. Итак, мы можем сделать вывод:
3
1). Если ранг расширенной матрицы больше ранга матрицы системы, то система несовместна.
2). Если ранг расширенной матрицы равен рангу матрицы системы, то сис-
тема совместна, при этом а) если ранг равен числу неизвестных, то система имеет единственное
решение;
б) если ранг меньше числа неизвестных, то система имеет бесконечно много решений.
Сформулируем без доказательства теорему.
Критерий Кронекера - Капелли. Система линейных уравнений совместна тогда и только тогда, когда ранг матрицы системы равен рангу расширенной мат-
рицы.
|
|
x1 2x2 x3 7, |
||||
|
|
|
3x3 |
10, |
||
Пример 6.2. Решить систему |
3x1 5x2 |
|||||
|
2x 7x |
|
x |
|
8, |
|
|
|
2 |
3 |
|||
|
1 |
|
|
|||
|
4x 11x |
2 |
x |
3 |
22. |
|
|
|
1 |
|
|
Решение. Составим расширенную матрицу системы A | B и обозначим ее
строки A1, A 2 , A3 , A 4 . С помощью элементарных преобразований строк получим нули ниже главной диагонали сначала в первом столбце с помощью первой стро-
ки (это равносильно исключению x1 из всех уравнений, кроме первого):
1 2 |
1 |
|
7 |
|
A1 |
|
|
|
|
|
A 1 |
|
|
1 2 |
1 |
|
7 |
|
||||
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
3A |
|
|
|
|||||||||
3 |
5 |
3 |
|
10 |
|
|
A2 |
|
|
A2 |
1 |
0 |
11 |
0 |
|
11 |
||||||
A | B 2 |
7 |
1 |
|
8 |
|
A |
|
|
~ |
A |
|
2 A |
|
|
0 |
3 |
3 |
|
6 |
. |
||
|
11 |
1 |
|
|
|
|
|
3 |
|
|
|
|
3 |
|
1 |
|
|
3 |
3 |
|
|
|
4 |
|
22 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
6 |
|||||
|
|
|
|
|
|
A |
4 |
|
A 4 4 A |
1 |
|
|
|
|
|
|
Разделим вторую строку на ( 11) . Третья и четвертая строки равны, сле-
довательно, одну из них можно вычеркнуть, например, четвертую. Третью строку можно сократить на 3. Затем с помощью второй строки получим нуль ниже глав-
ной диагонали во втором столбце:
4

1 |
2 |
1 |
|
A | B ~ |
0 |
1 |
0 |
|
0 |
1 |
1 |
|
7 |
|
|
|
B |
|
|
B |
1 |
|||
1 |
|
|
2 |
||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B3 |
~
B1 |
|
1 |
2 |
1 |
|
7 |
|
|
|
||||||||
B |
|
|
0 |
1 |
0 |
|
1 |
. |
2 |
|
|
||||||
|
|
|
|
|
|
|
||
|
|
|
0 |
0 |
1 |
|
3 |
|
B3 B2 |
|
|||||||
|
|
1 |
2 |
1 |
|
7 |
|
C |
|
|||
|
|
|
|||||||||
|
|
|
C |
1 |
. |
||||||
Умножим третью строку на 1 , получим |
A | B ~ |
0 |
1 |
0 |
|
1 |
|
|
2 |
||
|
|
0 |
0 |
1 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
C3 |
|
|
|
|
|
|
|
|
|
|
|
Закончен прямой ход метода Гаусса, т. е. ниже главной диагонали все эле-
менты равны нулю, элементы главной диагонали равны единице. Перейдем к об-
ратному ходу метода Гаусса – получению нулей выше главной диагонали сначала в третьем столбце с помощью третьей строки, потом во втором столбце с по-
мощью второй строки:
C1
A | B ~ C2
C3
|
|
C1 C3 |
|
|
~ |
|
C2 |
|
|
||
|
|
|
C3 |
|
|
|
|
1 |
2 |
0 |
|
4 |
|
|
D |
|||
|
|
||||||||||
|
|
|
D |
1 |
|||||||
|
0 |
1 |
0 |
|
1 |
|
|
2 |
|||
|
|
0 |
0 |
1 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
D3 |
||
|
|
|
|
|
|
|
|
|
|
D1 2D2 |
|
|
~ |
|
D2 |
|
|
||
|
|
|
D3 |
|
|
|
|
1 |
0 |
0 |
|
2 |
|
|
|
|||||||
|
|
||||||
|
0 |
1 |
0 |
|
1 |
. |
|
|
|
||||||
|
0 |
0 |
1 |
|
3 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Обратный ход метода Гаусса закончен – ниже и выше главной диагонали все элементы равны нулю. В результате элементарных преобразований матрицы
A B , а значит и системы уравнений, получили равносильную систему:
1 x1 0 x2 0 x3 2, |
x1 2, |
||||||||
|
|
1 x2 |
0 x3 |
1, или |
x2 |
1, |
|||
0 x1 |
|||||||||
0 x |
0 x |
2 |
1 x |
3 |
3, |
x |
3 |
3. |
|
|
1 |
|
|
|
|
|
Итак, в результате элементарных преобразований на месте матрицы A мы получили единичную матрицу, а на месте столбца свободных членов – столбец-
решение.
|
|
x1 2x2 3x3 5x4 1, |
|||||||
|
|
|
|
|
|
22x4 1, |
|||
Пример 6.3. Решить систему |
x1 3x2 13x3 |
||||||||
|
3x 5x |
|
x |
|
2x |
|
5, |
||
|
|
2 |
3 |
4 |
|||||
|
1 |
|
|
|
|
||||
|
|
2x 3x |
2 |
4x |
3 |
7x |
4 |
4. |
|
|
|
1 |
|
|
|
|
Решение. Составим расширенную матрицу системы и обозначим ее строки
A1, A2 , A3, A4 . С помощью элементарных преобразований получим нули ниже
5
главной диагонали в первом столбце с помощью первой строки. С помощью вто-
рой строки получим нуль ниже главной диагонали во втором столбце:
|
|
1 |
2 3 |
5 |
|
1 |
|
A1 |
|
|
|
A1 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
A 2 |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
1 |
3 13 |
22 |
|
|
1 |
|
A 2 |
|
|
A1 |
|
|
|
|
|
|
||||||||
A| B |
|
|
|
|
|
|
|
|
|
|
|
|
~ |
A3 |
|
|
|
|
|
|
|
|
||||
|
|
|
3 5 1 |
|
2 |
|
|
|
5 |
|
|
A3 |
|
3A1 |
|
|
|
|
|
|
||||||
|
|
2 3 4 |
|
7 |
|
|
|
4 |
|
|
|
|
A 4 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
A 4 |
|
2 A1 |
|
|
|
|
|
|
|||||
1 2 |
|
3 5 |
|
1 |
|
B1 |
|
B1 |
|
1 |
2 |
3 5 |
|
1 |
||||||||||||
|
|
|
|
|||||||||||||||||||||||
|
0 |
1 |
|
10 17 |
|
2 |
|
|
B |
|
|
|
B |
|
|
|
0 |
1 |
10 17 |
|
2 |
|
||||
|
|
|
|
|
2 |
~ |
|
|
2 |
|
|
. |
||||||||||||||
|
0 |
1 |
|
10 17 |
|
2 |
|
|
B3 |
|
B3 B2 |
|
0 |
0 |
0 |
0 |
|
0 |
|
|||||||
0 |
1 |
|
10 17 |
|
2 |
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
0 |
|
0 |
||||||
|
|
|
B4 |
|
B4 B2 |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
Заметим, что rang A | B rang A 2 ; значит, система совместна. Вычеркнем
две нулевые строки и получим нуль выше главной диагонали во втором столбце:
1 |
2 |
3 5 |
5 |
C |
~ |
C 2C |
1 |
0 |
17 29 |
5 |
. |
A| B ~ 0 |
1 |
10 17 |
2 |
C12 |
1 C2 |
2 0 |
1 |
10 17 |
2 |
В результате элементарных преобразований матрицы A | B (или, что то же самое, системы уравнений) получили равносильную систему
1 x1 0 x2 |
17 x3 29 x4 5, |
. |
|
|
|
|
|
|
||
|
1 x2 |
|
|
|
|
|
|
|
|
|
0 x1 |
10 x3 17 x4 2. |
|
|
|
|
|
|
|||
Здесь |
x1 , x2 зависимые |
неизвестные, |
|
x3 , x4 свободные неизвестные. |
||||||
Общее решение имеет вид: x1 |
5 17x3 29x4 , |
|
|
|
||||||
|
|
x2 2 10x3 17x4. |
|
|
|
|||||
Придавая свободным неизвестным x3 , x4 |
числовые значения, из общего |
|||||||||
решения можно получить частные решения системы. |
|
|||||||||
|
|
|
|
|
2x x |
2 |
3x |
3 |
9, |
|
Пример 6.4. |
|
|
|
1 |
|
|
|
|||
Решить систему |
|
3x1 5x2 |
x3 |
4, |
||||||
|
|
|
|
|
4x 7x |
2 |
x |
3 |
5. |
|
|
|
|
|
|
1 |
|
|
|
Решение. Применим метод Гаусса:
6

|
|
|
|
2 |
1 |
3 |
|
9 |
|
A 1 |
|
A |
2 A1 |
|
|
1 |
4 |
2 |
|
|
13 |
|
B 1 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
A |
|
B 3 |
5 |
1 |
|
4 A |
|
|
~ |
|
A |
|
|
3 |
5 1 |
|
|
4 B |
|
|
~ |
|||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
4 |
7 |
1 |
|
5 |
|
A 3 |
|
|
|
A 3 |
|
|
|
4 |
7 1 |
|
|
5 |
|
|
B 3 |
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
B 1 |
|
|
1 4 |
2 |
|
13 |
|
C1 |
~ |
|
|
C1 |
|
|
|
1 |
4 |
2 |
|
13 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
~ B |
|
3B |
|
|
0 7 |
7 |
|
35 C |
|
|
|
C |
|
7 |
|
0 |
1 1 |
|
5 . |
|
||||||||||||||
|
2 |
|
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
C3 |
|
C3 |
9 C2 |
|
7 |
|
|
|
|
|
|
|
|
||||||||||
B3 |
4B 1 |
0 9 |
9 |
|
57 |
|
|
|
|
0 |
0 0 |
|
12 |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
Последней строке соответствует уравнение 0 x1 0 x2 0 x3 12 , которое не имеет решения. Следовательно, и система не имеет решения.
7