Решение систем линейных уравнений методом Гаусса.
Одним из наиболее универсальных и эффективных методов решений линейных алгебраических систем является метод Гаусса, состоящий в последовательном исключении неизвестных.
Пусть дана система уравнений
,
, (3)
.
Процесс решения по методу Гаусса состоит из двух этапов. На первом этапе (прямой ход) система приводится к ступенчатому (в частности, треугольному) виду.
Приведенная ниже система имеет ступенчатый вид
,
,
,
где k ≤ n , aii ≠ 0 , . Коэффициенты aiiназываются главными элементами системы.
На втором этапе (обратный ход) идет последовательное определение неизвестных из этой ступенчатой системы.
Замечания: 1. Если ступенчатая система оказывается треугольной, т. е. k = n, то исходная система имеет единственное решение. Из последнего уравнения находим хn, из предпоследнего уравнения xn-1, далее поднимаясь по системе вверх, найдем все остальные неизвестные (xn-2, ... , x1).
2. На практике удобнее работать не с системой (3), а с расширенной ее матрицей, выполняя все элементарные преобразования над ее строками. Удобно, чтобы коэффициент а11 был равен 1 (уравнения переставить местами, либо разделить обе части уравнения на a11 ≠ 1 ).
Пример. Решить систему методом Гаусса:
,
,
.
Решение: Произведем элементарные преобразования над строчками расширенной матрицы системы:
~ ~ ~
Полученная матрица соответствует системе
x1 + x2 + x3 = 3 ,
x2 = 1 ,
x3 = 1 .
Осуществляя обратный ход, находим хз = 1, х2 = 1, х1 = 1.
Системы линейных однородных уравнений.
Пусть дана система линейных однородных уравнений
,
,
.
Очевидно, что однородная система всегда совместна (r(А) = r( )), она имеет нулевое (тривиальное) решение х1 = х2 = . . . = хn = 0.
При каких условиях однородная система имеет и ненулевые решения?
Теорема 4. Для того, чтобы система однородных уравнений имела ненулевые решения, необходимо и достаточно, чтобы ранг r ее основной матрицы был меньше числа п неизвестных, т. е. r < п.
Теорема 5. Для того, чтобы однородная система п линейных уравнений с п неизвестными имела ненулевые решения, необходимо и достаточно, чтобы ее определитель Δ был равен нулю, т. е. Δ = 0.
Если система имеет ненулевые решения, то Δ = 0. Ибо при Δ ≠ 0 система имеет только единственное, нулевое решение. Если же Δ = 0, то ранг r основной матрицы системы меньше числа неизвестных, т. е. r < n. И, значит, система имеет бесконечное множество (ненулевых) решений.
Пример. Решить систему
x1 – 2x2 + 4x3 = 0 ,
2x1 – 3x2 + 5x3 = 0 .
Решение: , r(A) = 2 ( ) , n = 3 .
Так как r < n , то система имеет бесчисленное множество решений. Найдем
их.
x1 – 2x2 = – 4x3 ,
2x1 – 3x2 = – 5x3 .
, . Стало быть,
, – общее решение.
Положив x3 = 0 , получаем одно частное решение: x1 = 0, х2 =0, x3 = 0. Положив х3 = 1, получаем второе частное решение: x1 = 2, х2 = 3, x3 = 1 и т. д.