Метод Гаусса.

Практическое значение правила Крамера для решения системы n линейных уравнений с п неизвестными невелико, так как при его применении приходится вычислять п +1 определителей n-го порядка: , x1, x2, …,xn. Более удобным является так называемый метод Гаусса. Он применим и в более общем случае системы линейных уравнений, т. е. когда число уравнений не совпадает с числом неизвестных.

Итак, пусть дана система, содержащая m линейных уравнений с п неизвестными:

а11х1 + а12х2 + …+ а1nхn = b1;

а21х1 + а22х2 + …+ а2nхn = b2; (6)

. ……………………………………

аm1х1 + аm2х2 + …+ аmnхn = bm

Требуется найти все решения системы уравнений (6). Будем производить над системой элементарные преобразования: исключение из системы уравнения вида

0х1 + 0х2 + …+ 0хn = 0 (7)

и прибавление к обеим частям одного из уравнений системы соответствующих частей другого уравнения, умноженных на любое число .

Очевидно, что если мы проделаем над уравнениями системы (6) любое из приведенных выше преобразований, то получим систему, равносильную исходной. При необходимости систему (6) будем подвергать еще одному виду преобразований - перенумерации переменных и уравнений. Идея этого преобразования заключается в следующем. Если, например, возникает необходимость, чтобы в каком-то уравнении системы (например, в k -м) неизвестная x1 стояла на первом месте, то в результате перенумерации соответствующее уравнение запишется в виде

akix1 +... + ak2x2 + … + ak1xi+... + aknxn = bk,

т. е. вместо прежней неизвестной хi мы будем писать х1, а вместо x1 - хi Метод Гаусса решения системы (6) заключается в последовательном исключении переменных.

Если среди уравнений системы есть хотя бы одно уравнение вида

0xl + 0x2+... + 0xn= b, (8)

причем b 0, то совершенно очевидно, что ни одна система значений х1, х2..., хп не удовлетворяет этому уравнению, а следовательно, и системе в целом, поэтому система несовместна.

Пусть теперь система (6) не содержит уравнений вида (7) или (8). Это значит, что в каждом уравнении системы хотя бы один из коэффициентов отличен от нуля. Пусть a110 (в противном случае, применив элементарные преобразования, мы сможем добиться, чтобы первый коэффициент первого уравнения был отличен от нуля). Оставив первое уравнение без изменения, исключим из всех уравнений системы (6), начиная со второго, неизвестную х1. Для этого из второго уравнения вычтем первое, умноженное на a21/a11, затем из третьего уравнения вычтем также первое, но уже умноженное на a31/a11, и так до последнего уравнения. В результате этих преобразований мы получим равносильную систему

а11х1 + а12х2 + … + а1nхn = b1;

а?22х2 + …+ а?2nхn = b?2;

………………………… (9)

а?m2х2 + …+ а?mnхn = b?m

Заметим, что в системе (9) число уравнений может быть и меньше m, так как среди них могут оказаться уравнения вида (7), которые, как мы условились ранее, можно отбросить.

Пусть а22 0. Применим те же самые рассуждения и исключим из последних п - 2 уравнений системы (9) неизвестную х2 путем вычитания из третьего уравнения второго, умноженного на a?32/a?22, из четвертого уравнения -- второго, умноженного на a?34/a?22 и т. д. В результате получим систему

а11х1 + а12х2 + а13х3 + …+ а1nхn = b1;

а?22х2 + а?23х3 + …+ а?2nхn = b?2;

а?22х2 + а?23х3 + …+ а?2nхn = b?2;

а??33х3 + …+ а?3хn = b?3;

……………………………

а?m3х3 +…+а?mnхn = b?m

.

Продолжая этот процесс, систему (6) приведем к равносильной системе вида

c11х1 + c12х2 + c13х3 + …+ c1kхk + …+ c1nхn = d1;

c22х2 + c23х3 + …+ c2kхk + …+ c2nхn = d2;

c33х3 + …+ c3kхk + …+ c3nхn = d3; (10)

………………………………………

ckkхk + …+cknхn = dk.

в которой коэффициенты c11, c22,..., ckk отличны от нуля.

Может оказаться, что в процессе преобразования на каком-то шаге в полученной системе окажется уравнение вида (8). В этом случае система (7) не имеет решений. Предположим теперь, что среди уравнений полученной системы нет уравнения вида (8). Тогда для решения системы (6) необходимо решить систему (9), что не составляет особого труда. Рассмотрим два возможных случая.

1. k = n (это частный случай, когда число уравнений совпадает с числом неизвестных). Тогда последнее уравнение системы (10) имеет вид сппхп = dn, откуда хп = dn /cnn. Подставив это значение в предпоследнее уравнение системы (7), имеющее вид

cn-1n-1xn-1 + cn-1nxn= dn-1, найдем значение неизвестной xn-1 и т. д.; наконец, из первого уравнения найдем неизвестную x1 Таким образом, в случае k = п система уравнений (6) имеет единственное решение.

2. k < n. Тогда из последнего уравнения системы (10), найдем неизвестную xk, выраженную через неизвестные хk+1, хk+2,... xn:

xk = (dkk - ck k+1xk+1 - … - cknxn)

Подставив это значение неизвестной в предпоследнее уравнение системы (10), найдем выражение для неизвестной хk-1,и т. д.; наконец, подставив значения неизвестных хk, хk-1,... x2 в первое уравнение системы (10), получим выражение для неизвестной x1. В результате указанная система уравнений (6) приводится к виду

x1 = d?1 + c?1 k+1xk+1 + …+ c?1nxn;

x2 = d?2 + c?2 k+1xk+1 + …+ c?2nxn; (11)

………………………………………

xk = d?k + c?k k+1 xk+1 + …+ cknxn.

Неизвестные хk+1, хk+2, …, хп называются свободными. Им можно придать различные значения и затем из системы (6) найти значения неизвестных х1, х2, …, хk. Таким образом, в случае k < п совместная система уравнений (6) имеет бесчисленное множество решений.

Заметим, что если в процессе приведения системы (6) к системе (11) была произведена перенумерация неизвестных, то в системе (11) необходимо вернуться к их первоначальной нумерации.

На практике процесс решения системы уравнений облегчается тем, что указанным выше преобразованиям подвергают не саму систему, а матрицу, составленную из коэффициентов уравнений системы (6) и их свободных членов. При этом каждому элементарному преобразованию, проведенному над системой (6), соответствует преобразование над матрицей (12): вычеркивание строки, все элементы которой состоят из нулей, прибавление к элементам некоторой строки соответствующих элементов другой строки, умноженных на некоторое число, и перестановка двух столбцов матрицы (12).

Пример 1. Решить методом Гаусса систему уравнений

x1 - 2x2 + x3 + x4 = -1;

3x1 + 2x2 - 3x3 - 4x4 = 2;

2x1 - x2 + 2x3 - 3x4 = 9;

x1 + 3x2 - 3x3 - x4 = -1