Рассмотрим СЛАУ
(3.17)
Запишем расширенную прямоугольную матрицу коэффициентов системы
(3.17):
. (3.18)
Среди элементов матрицы
aij (i,j = 1, ...n) выберем наибольший по модулю, называемый главным элементом. Пусть им будет, например, элемент apq. Строка, содержащая главный элемент, называется главной строкой.
Далее вычисляем множители
mi = aiq / apq для всех i неравных p.Затем преобразуем матрицу (3.18) следующим образом: из каждой i-ой неглавной строки вычитаем почленно главную строку, умноженную на mi. В результате получим матрицу, у которой все элементы q-го столбца за исключением apq, равны 0. Отбрасывая этот столбец и главную строку, получим новую матрицу M1 с числом строк и столбцов на 1 меньше.
Над матрицей М
1 повторяем те же операции, после чего получим матрицу M2 и т.д. Таким образом продолжаем до тех пор, пока не получим матрицу, содержащую одну строку из двух элементов, которую тоже считаем главной.
Затем объединим все главные строки, начиная с последней. После некоторой перестановки они образуют треугольную матрицу, эквивалентную исходной. На этом заканчивается этап вычислений, называемый прямым ходом. Решив систему с полученной треугольной матрицей коэффициентов, найдём последовательно значения неизвестных
xi(i = 1, 2, ..., n). На этом заканчивается обратный ход.
Смысл выбора главного элемента состоит в том, чтобы сделать возможно меньшими числа
mi и тем самым уменьшить погрешность вычислений.
Пример. Решить систему уравнений.
Решение. Начинаем действовать так же, как и в предыдущем примере. Пусть .Тогда прямым ходом метода Гаусса исключаем из второго и третьего уравнений переменную :
Коэффициент при во втором уравнении был обозначен нами в предыдущем пункте через :
Если теперь , то можно исключить переменную из последнего полученного уравнения:
Коэффициент при в этом уравнении был обозначен в предыдущем пункте через . Если выражение
отлично от нуля, то последнее уравнение системы однозначно разрешимо относительно :
Подставляем найденное значение в уравнение
получаем:
Наконец, подставляя найденные значения для и в первое уравнение системы, находим значение :
Ответ. При условиях
система имеет единственное решение.
Итак, мы получили общие формулы решения системы уравнений. Правда, мы разобрали пока только случай единственности решения. Условия существования единственного решения получились в виде набора неравенств на коэффициенты системы. Оказывается, что не все эти неравенства являются существенными для единственности, можно доказать, что необходимым и достаточным в каждом примере является последнее полученное неравенство: для системы