Лекции по линейной алгебре / 5_метод Гаусса
.pdfЛекция 5. Метод Гаусса
План:
1.Метод Гаусса;
2.Теорема Кронекера-Капелли;
3.Примеры решения систем.
|
5.1. Решение систем методом Гаусса |
|
||||||
|
Рассмотрим систему |
m линейных уравнений с n |
неизвестными |
|||||
a11x1 |
a12 x2 |
... |
a1 j x j |
... |
a1n xn |
b1; |
|
|
a21x1 |
a22 x2 |
... |
a2 j x j |
... |
a2n xn |
b2 ; |
|
|
|
|
|
... |
|
|
|
(5.1) |
|
ai1x1 |
ai 2 x2 |
... |
aij x j |
... |
ain xn |
bi ; |
||
|
||||||||
|
|
|
... |
|
|
|
|
|
am1x1 |
am2 x2 |
... |
amj x j |
... |
amn xn |
bm . |
|
|
|
Метод |
|
Гаусса |
|
заключается в последовательном |
исключении |
переменных с помощью элементарных преобразований. С помощью таких преобразований система уравнений (5.1) приводится к равносильной системе ступенчатого или треугольного вида (прямой ход метода Гаусса). Из полученной системы, начиная с последних (по номеру) переменных, находятся все остальные переменные (обратный ход метода).
Пусть |
в системе (5.1) коэффициент при переменной x1 в |
первом |
|
уравнении |
a11 |
0 (путем перестановок уравнений системы всегда можно |
|
добиться, что a11 |
0 ). |
|
|
Шаг 1. Умножая первое уравнение на числа: a21 / a11, a31 / a11,..., |
am1 / a11 и |
прибавляя полученные уравнения соответственно ко второму, третьему,…,
m -му уравнению |
системы (5.1), исключается переменная x1 из всех |
||||||
последующих уравнений, начиная со второго. Т.е. получаем |
|
||||||
a11x1 |
a12 x2 |
... |
a1n xn |
b1, |
|
||
a(1) x ... |
a(1) x |
b(1) |
, |
|
|||
22 |
2 |
2n n |
2 |
|
|
||
..................................... |
(5.2) |
||||||
a(1) x ... |
a(1) x |
b(1) |
, |
||||
|
|||||||
i 2 |
2 |
in |
n |
i |
|
|
|
..................................... |
|
||||||
a(1) x ... |
a(1) x |
b(1) |
, |
|
|||
m2 |
2 |
mn n |
m |
|
|
где буквами с верхним индексом (1) обозначены новые коэффициенты, полученные после первого шага.
Шаг 2. Предположим, что a22(1) 0 (всегда можно путем перестановок уравнений системы получить данный факт). Умножая второе уравнение на
числа a(1) |
/ a(1) |
, |
a(1) |
/ a(1) |
,..., |
a(1) |
/ a(1) |
и прибавляя полученные уравнения |
32 |
22 |
|
42 |
22 |
|
m2 |
22 |
|
соответственно к третьему, |
четвертому,…., |
m -му уравнению системы (5.2), |
||||||||||||
исключим |
переменную |
x2 |
из |
всех |
последующих уравнений, начиная с |
|||||||||
третьего. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Продолжая |
процесс |
последовательного исключения переменных |
||||||||||||
x3 , x4 ,..., xr 1 , после (r |
1) -го шага получим систему: |
|||||||||||||
a11 x1 |
a12 x2 |
... |
a1r xr |
a1,r 1 xr |
1 |
... |
|
a1n xn |
b1 , |
|||||
a(1) x ... |
a(1) x |
a(1) x |
|
... |
a(1) x |
b(1) |
, |
|||||||
22 2 |
|
2r r 2,r 1 r 1 |
|
|
2n n |
2 |
|
|||||||
........................................................................ |
||||||||||||||
|
a(r |
1) x |
a(r |
1) x |
... |
a(r |
1) x |
|
b(r 1) |
, |
|
|||
|
rr |
r |
r ,r 1 r 1 |
|
|
rn |
n |
|
r |
|
(5.3) |
|||
..................................................................... |
||||||||||||||
|
||||||||||||||
|
|
|
|
0 |
b(r |
1) , |
|
|
|
|
|
|
|
|
|
|
|
|
|
r 1 |
|
|
|
|
|
|
|
|
........................................
0 bm(r 1) .
В |
последних (m |
r) |
уравнениях |
слева |
записаны |
выражения вида |
|||
0 x1 0 |
x2 |
... 0 xn . |
|
|
|
|
|
|
|
Если |
хотя бы |
одно |
из чисел |
b(r 1) |
,...,b(r 1) |
не |
равно нулю, то |
||
|
|
|
|
|
r 1 |
|
m |
|
|
соответствующее равенство противоречиво, и система (5.1) несовместна.
На практике преобразования Гаусса как правило проводят не с самими уравнениями, а с матрицей их коэффициентов (берут расширенную матрицу).
2x1 x2 6,
Пример1. Решить систему x1 2x2 x3 5, методом Гаусса.
3x1 4x2 2x3 13
Решение. Преобразования будем выполнять с расширенной матрицей системы
|
2 |
1 |
0 |
6 |
A | B |
1 |
2 |
1 |
5 . Приведем расширенную матрицу к |
|
3 |
4 |
2 |
13 |
|
|
|
|
|
ступенчатому виду (прямой ход метода Гаусса).
1) Для удобства преобразований поменяем первую и вторую строчки (это элементарное преобразование приводит к эквивалентной матрице):
1 |
2 |
1 |
5 |
2 |
1 |
0 |
6 ; |
3 |
4 |
2 |
13 |
|
|
|
|
2) Элемент, стоящий в первой строке и в первом столбце, будет ведущим. Элементы, находящиеся в первом столбце ниже ведущего необходимо обнулить. Для этого к элементам второй строки прибавляем удвоенные соответствующие элементы первой строки. Из элементов третьей
строки вычитаем соответствующие элементы первой строки, умноженные на 3. Получаем эквивалентную матрицу:
1 2 1 5
0 3 2 4 ;
0 10 1 2
3) Для упрощения преобразований умножим все элементы третьей строки на число 3:
1 2 1 5
0 3 2 4 ;
0 30 3 6
4) Для получения следующей "ступеньки" ведущим элементом выбираем элемент "-3", который находится во второй строке и во втором столбце. Поэтому необходимо обнулить элемент, находящийся под ним, для этого к элементам третьей строки прибавляем соответствующие элементы второй строки, умноженные на 10:
1 |
2 |
1 |
5 |
|
0 |
3 |
2 |
4 |
; |
0 |
0 |
17 |
34 |
|
|
|
|
|
|
5) Элементы третьей строки разделим на "-17":
1 2 1 5
0 3 2 4 .
0 0 1 2
Без комментариев процедура решения записывается в виде:
2 |
1 |
0 |
6 |
|
1 |
2 |
1 |
5 |
|
1 |
2 |
1 |
5 |
|
1 |
2 |
1 |
5 |
~ |
2 |
1 |
0 |
6 +2I ~ |
0 |
3 |
2 |
4 |
~ |
|
3 |
4 |
2 |
13 |
|
3 |
4 |
2 |
13 |
-3I |
0 |
10 |
1 |
2 |
*3 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
1 |
5 |
1 |
2 |
1 |
5 |
1 |
2 |
1 |
5 |
~ 0 |
3 |
2 |
4 |
~ 0 |
3 |
2 |
4 |
~ 0 |
3 |
2 |
4 . |
0 |
30 |
3 |
6 +10II |
0 |
0 |
17 |
34 :(-17) |
0 |
0 |
1 |
2 |
По полученной ступенчатой матрице записываем соответствующую систему:
x1 2x2 |
x3 |
5, |
3x2 |
2x3 |
4, |
x3 |
2. |
|
Выполняем обратный ход:
|
x3 |
2, |
|
|
|
|
|
x3 2, |
x1 |
3, |
|||
|
1 |
|
|
|
|
1 |
|
|
|||||
x |
|
2x 4 , |
|
x |
|
2 ( 2) 4 0, |
x |
0, |
|||||
|
|
|
|
|
|
|
|||||||
2 |
3 |
3 |
|
2 |
3 |
|
2 |
|
|||||
|
|
|
|
|
x3 |
2. |
|||||||
x1 |
2x2 |
x3 5, |
x1 |
2x2 |
x3 |
5 2 0 ( 2) 5 3, |
|||||||
|
|
||||||||||||
Ответ: Ответ: x1 |
3; x2 |
0; x3 |
2 . |
|
|
|
5.2. Теорема Кронекера-Капелли
Система линейных уравнений совместна тогда и только тогда, когда ранг матрицы системы равен рангу расширенной матрицы этой системы.
Пусть A основная матрица системы, A1 - расширенная матрица системы (5.1).
Результаты исследования решения системы (5.1) можно представить в виде схемы:
Система m линейных уравнений с
n неизвестными
r( A) r( A1 ) |
|
r( A) r( A1 ) r |
|
Система несовместна |
|
|
|
|
Система совместна |
||
|
|
|
|
r n r n
Система неопределенная |
Система определенная |
|
Единственное решение |
||
решений |
||
|
Вслучае неопределенной системы r переменных x1, x2 , ..., xr
называются основными (базисными) переменными, если определитель матрицы из коэффициентов при них отличен от нуля. Остальные n r переменные называются неосновными (свободными).
Решение системы (5.1), в котором все n r неосновных переменных равны нулю, называется базисным.
5.3. Примеры решения систем
|
|
|
|
|
|
|
2x1 |
3x2 |
x3 |
|
x4 |
5, |
|
|
|
|
|
|
|
|
|
|
а) Решить систему |
|
3x1 |
x2 |
2x3 |
|
x4 |
1, |
|
|
|
|
|
|
|
|
|
||||||
|
x1 |
2x2 |
3x3 |
4x4 |
6, |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
6x1 |
4x2 |
4x3 |
|
6x4 |
1. |
|
|
|
|
|
|
|
|
||
Решение. Решим систему методом Гаусса. |
|
|
|
|
|
|
|
|||||||||||||||
Прямой ход: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
3 |
1 |
1 |
5 |
|
1 |
2 |
3 |
4 |
6 |
|
-3I |
|
|
1 |
2 |
3 |
4 |
|
6 |
|
|
3 |
1 2 1 |
1 |
|
3 |
1 2 1 |
1 |
|
|
|
0 7 |
7 |
11 |
|
17 (-1) |
|
|||||||
~ |
|
-2I |
~ |
~ |
|
~ |
||||||||||||||||
1 |
2 |
3 |
4 |
6 |
2 |
3 |
1 |
1 |
5 |
|
0 |
1 |
7 |
7 |
|
7 |
(-1) |
|||||
|
|
-6I |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
6 |
4 4 6 |
1 |
|
6 |
4 4 6 |
1 |
|
|
|
0 |
8 14 18 |
|
35 (-1) |
|
||||||||
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
6 |
|
1 |
2 |
3 |
4 |
6 |
|
|
1 |
2 |
3 |
4 |
6 |
|
|
~ |
0 |
7 |
7 |
11 |
17 |
~ |
0 |
1 |
7 |
7 |
7 |
|
~ |
0 |
1 |
7 |
7 |
7 |
~ |
|
0 |
1 |
7 |
7 |
7 |
0 |
7 |
7 |
11 |
17 |
-7II |
0 |
0 |
42 |
38 |
32 |
|||||
|
|
|
|
|||||||||||||||||
|
0 |
8 |
14 |
18 |
35 |
|
0 |
8 |
14 |
18 |
35 |
-8II |
|
0 |
0 |
42 |
38 |
21 |
-III |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
6 |
~ |
0 |
1 |
7 |
7 |
7 . |
|
0 |
0 |
42 |
38 |
32 |
|
0 |
0 |
0 |
0 |
11 |
|
|
|
|
|
|
Ранг расширенной матрицы равен 4 (количество ненулевых строк в матрице ступенчатого вида). Ранг основной матрицы равен 3 (убираем последний столбец свободных членов и вычеркиваем последнюю нулевую строку).
Значит, по теореме Кронекера-Капелли система несовместна, т. е. не имеет решений.
Последняя строка полученной матрицы дает противоречивое уравнение 0 x1 0 x2 0 x3 0 x4 11 .
Ответ: система не имеет решения.
2x1 x2 3x3 2x4 4x5 |
1, |
б) Решить систему: 4x1 2x2 5x3 x4 7x5 2, 2x1 x2 x3 8x4 2x5 7.
Решение. Прямой ход метода Гаусса. Приводим расширенную матрицу системы к ступенчатому виду:
2 |
1 |
3 |
|
2 |
4 |
1 |
|
2 |
1 |
3 |
|
2 |
4 |
|
1 |
|
|
-2I ~ |
|
|
|
||||||||||||
4 |
2 |
5 |
1 |
|
7 |
2 |
~ 0 |
0 |
|
1 |
5 |
|
1 |
4 |
~ |
|
2 |
1 |
1 |
8 |
|
2 |
7 |
-I |
0 |
0 |
|
2 |
10 |
|
2 |
8 |
-2II |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
1 |
3 |
2 |
4 |
|
1 |
|
2x1 x2 |
x3 |
2x4 |
4x5 |
1, |
|
~ |
0 |
0 |
|
1 |
5 |
|
1 |
4 |
~ |
|
|||||
|
|
x3 |
5x4 x5 |
4. |
|
|
|||||||||
|
0 |
0 |
|
0 |
0 |
0 |
|
0 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||||||
Ранги основной и расширенной матриц равны 2. Значит, система |
|||||||||||||||
совместна. Так как r |
2 меньше числа переменных n 5 , |
то система имеет |
|||||||||||||
бесконечно много решений. |
|
|
|
|
|
|
|
||||||||
Количество |
базисных |
переменных |
равно |
рангу |
r 2 . В качестве |
базисных можно взять переменные: x1, x3 . Свободные переменные: x2 , x4 , x5 . Обратный ход метода Гаусса:
2x1 |
x2 |
x3 |
2x4 |
|
4x5 |
1, |
|
|
x3 |
5x4 |
x5 |
|
4, |
|
|
|
|
|
|
|
||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
x 5x x |
4, |
x |
x x 2x 4x 1 , |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
3 |
|
4 |
5 |
|
1 |
2 |
2 |
3 |
4 |
|
|
|
5 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
x3 |
|
5x4 |
x5 4, |
|
|
|
|
|
x3 |
|
|
5x4 |
x5 |
|
4, |
|
|
|||
x |
1 |
x 5x x 4 2x 4x 1 , |
x |
|
1 |
x |
|
|
3 |
x |
3 |
x |
3 |
. |
||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
1 |
2 |
2 |
|
4 |
5 |
4 |
5 |
|
|
1 |
2 |
|
2 |
2 |
|
4 |
2 |
5 |
2 |
|
Ответ: |
1 |
x |
3 |
x |
3 |
x |
3 |
; x ; 5x |
x |
4; x ; x |
, где |
x , x , x R . |
||||
|
|
|
|
|||||||||||||
|
2 |
2 |
2 |
4 |
2 |
5 |
2 |
2 |
4 |
5 |
4 |
5 |
|
2 |
4 |
5 |