Пример 1
.
Проверка:
Ответ: (3; -1).
Пример 2
Проверка:
Ответ: x=0,5; y=2; z=1,5.
Более удобным является так называемый метод Гаусса. Он применим и в более общем случае системы линейных уравнений, т. е. когда число уравнений не совпадает с числом неизвестных.
Итак, пусть дана система, содержащая m линейных уравнений с п неизвестными:
а11х1 + а12х2 + …+ а1nхn = b1;
а21х1 + а22х2 + …+ а2nхn = b2;
. ……………………………………
аm1х1 + аm2х2 + …+ аmnхn = bm
Метод Гаусса решения системы заключается в последовательном исключении переменных.
Схема единственного деления. Рассмотрим простейший вариант метода Гаусса, называемый схемой единственного деления.
Прямой ход состоит из n - 1 шагов исключения.
|
|
1-й шаг. Целью этого шага является исключение неизвестного x 1 из уравнений с номерами i = 2, 3, …, n. Предположим, что коэффициент a 11 ¹ 0. Будем называть его главным элементом 1- го шага.
Найдем величины
qi 1 = ai 1/ a 11 (i = 2, 3, …, n),
называемые множителями 1- го шага. Вычтем последовательно из второго, третьего, …, n- го уравнений системы первое уравнение, умноженное соответственно на q2 1, q 31, …, qn 1. Это позволит обратить в нуль коэффициенты при x 1 во всех уравнениях, кроме первого. В результате получим эквивалентную систему
a 11 x 1 + a 12 x 2 + a 13 x 3 + … + a 1 nxn = b 1 ,
a 22(1) x 2 + a 23(1) x 3 + … + a 2 n (1) xn = b 2(1) ,
a 32(1) x 2 + a 33(1) x 3 + … + a 3 n (1) xn = b 3(1),
...............
an 2(1) x 2 + an 3(1) x 3 + … + ann (1) xn = bn (1).
в которой aij (1) и bij (1) вычисляются по формулам
aij (1) = aij − qi 1 a 1 j , bi (1) = bi − qi 1 b 1.
2-й шаг. Целью этого шага является исключение неизвестного x 2 из уравнений с номерами i = 3, 4, …, n. Пусть a 22(1) ≠ 0, где a 22(1) – коэффициент, называемый главным (или ведущим) элементом 2- го шага. Вычислим множители 2-го шага
qi 2 = ai 2(1) / a 22(1) (i = 3, 4, …, n)
и вычтем последовательно из третьего, четвертого, …, n-о го уравнения системы второе уравнение, умноженное соответственно на q 32, q 42, …, qm 2. В результате получим систему
a 11 x 1 + a 12 x 2 + a 13 x 3 + … + a 1 nxn = b 1 ,
a 22(1) x 2 + a 23(1) x 3 + … + a 2 n (1) = b 2(1) ,
a 33(2) x 3 + … + a 3 n (2) xn = b 3(2) ,
...................
an 3(2) x 3 + … + ann (2) xn = bn (2) .
Здесь коэффициенты aij (2) и bij (2) вычисляются по формулам
aij (2) = aij (1) – qi 2 a 2 j (1), bi (2) = bi (1) – qi 2 b 2(1).
Аналогично проводятся остальные шаги. Опишем очередной k- й шаг.
k- й шаг. В предположении, что главный (ведущий) элемент k-о го шага akk ( k –1) отличен от нуля, вычислим множители k-го шага
|
|
qik = aik (k –1) / akk (k –1) (i = k + 1, …, n)
и вычтем последовательно из (k + 1)-го, …, n -го уравнений полученной на предыдущем шаге системы k -e уравнение, умноженное соответственно на qk +1, k, qk +2, k, …, qnk.
После (n - 1)-го шага исключения получим систему уравнений
a 11 x 1 + a 12 x 2 + a 13 x 3 + … + a 1 n xn = b 1 ,
a 22(1) x 2 + a 23(1) x 3 + … + a 2 n (1) xn = b 2(1) ,
a 33(2) x 3 + … + a 3 n (2) xn = b 3(2) ,
....................
ann (n –1) xn = bn (n –1) .
матрица A ( n -1) которой является верхней треугольной. На этом вычисления прямого хода заканчиваются.