Метод последовательного полного исключения неизвестных

Березнёва Т. Д.

Тема 7

СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ.

МЕТОД ГАУССА – ЖОРДАНА.»

( Учебная дисциплина “Введение в линейную алгебру и аналитическую геометрию”)

СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ.

МЕТОД ГАУССА – ЖОРДАНА.

Основные понятия

Уравнение с n переменными называется линейным, если все переменные (x₁,x₂, …x_n) входят в него в степени 1. Общий вид такого уравнения формально записывается следующим образом:

a₁x₁ + a₂x₂ + … a_jx_j + … a_nx_n = b, (*)

или

= b.

Величины a_j, j = 1,…,n, и b являются известными (заданными). Величины a_j называются коэффициентами при переменных (при неизвестных), а b - свободным членом.

Решением линейного уравнения (*) называется упорядоченный набор (, ,…, ) значений переменных, который при подстановке в уравнение (т.е. при заменеx_j на при всех j от1до n обращает его в тождество. Подчеркнем, что решение уравнения с n переменными всегда есть набор из n чисел и каждый такой набор из n чисел представляет собой одно решение. Очевидно, что если хотя бы один коэффициент при переменных не равен 0, то уравнение (*) имеет решение. В противном случае решение существует только при b = 0, и это все произвольные наборы из n чисел.

Рассмотрим одновременно m уравнений вида (*), т.е. систему m линейных алгебраических уравнений с n переменными. Пусть каждое i - е уравнение, i = 1,2,…,m, задается коэффициентами при переменных a_i₁, a_i₂, …, a_in и свободным членом b_i, т.е. имеет вид

a_i1x₁ + a_i2 x₂ + … + a_ij x_j + … + a_in x_n = b_i.

Тогда в общем виде система m линейных алгебраических уравнений с n переменными может быть записана в виде:

a₁₁x₁ + a₁₂ x₂ + … + a_1j x_j + … + a_1n x_n = b₁

a₂₁x₁ + a₂₂ x₂ + … + a_2j x_j + … + a_2n x_n = b₂

_{………………………………………………………………………………}

a_i1x₁ + a_i2 x₂ + … + a_ij x_j + … + a_in x_n = b_i (1)

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

a_m1x₁ + a_m2 x₂ + … + a_mj x_j + … + a_mn x_n = b_m

или, что то же самое,

= b_i, i = 1,…,m.

Если все свободные члены равны нулю, то система (1) называется однородной, т.е. имеет вид

= 0, i = 1,…,m, (1⁰)

в противном случае - неоднородной. Система (1⁰) является частным случает общей системы (1).

Решением системы уравнений (1) называется упорядоченный набор (, ,…, ) значений переменных, который при подстановке в уравнения системы (1) (т.е. при заменеx_j на , j = 1,…,n) все эти уравнения обращает в тождества, т.е. = b_i при всех i = 1,…,m.

Система уравнений (1) называется совместной, если у нее существует хотя бы одно решение. В противном случае система называется несовместной.

Совокупность всех решений системы уравнений (1) мы будем называть множеством ее решений и обозначать X^b(X⁰, если система однородная). Если система несовместна, то X^b = .

Основная задача теории систем линейных алгебраических уравнений состоит в том, чтобы выяснить, совместна ли система (1), и, если совместна, то описать множество всех её решений. Существуют методы анализа таких систем, которые позволяют описывать множество всех решений в случае совместных систем или убеждаться в несовместности в противном случае. Одним из таких универсальных методов является метод последовательного полного исключения неизвестных, или метод Гаусса - Жордана, который мы будем подробно изучать.

Прежде, чем переходить к описанию метода Гаусса - Жордана, приведем ряд полезных для дальнейшего определений и утверждений.

Две системы уравнений называются эквивалентными, если они имеют одно и то же множество решений. Другими словами, каждое решение одной системы является решением другой, и наоборот. Все несовместные системы считаются эквивалентными между собой.

Из определений эквивалентности и множества решений систем вида (1) сразу же вытекает справедливость следующих утверждений, которые мы сформулируем в виде теоремы.

Теорема 1. Если в системе (1) имеется уравнение с номером k, 1 k m, такое, что a_kj = 0 j, то

1) если соответствующий свободный член b_k 0, то система (1) несовместна;

2) если соответствующий свободный член b_k = 0, то k - е уравнение можно отбросить и получить систему из (m – 1) - го уравнения с n переменными, эквивалентную исходной.

Справедливость утверждений теоремы становится очевидной, если заметить, что k – е уравнение имеет вид

0x₁ + 0 x₂ + … + 0 x_j + … + 0 x_n = b_k.

Теорема 2. Если к одному уравнению системы (1) прибавить другое уравнение этой же системы, умноженное на любое число, то получится система уравнений, эквивалентная исходной системе.

Доказательство. Умножим, например, второе уравнение системы (1) на некоторое число и прибавим его к первому уравнению. В результате этого преобразования получим систему (1’), в которой все уравнения, начиная со второго, не изменились, а первое имеет следующий вид

= b₁ + b₂.

Очевидно, если какой-нибудь набор (, ,…, ) значений переменных обращает в тождества все уравнения системы (1), то он обращает в тождества и все уравнения системы (1’). Наоборот, решение (x’₁,x’₂,…,x’_j, …,x’_n) системы (1’) является также решением системы (1), так как система (1) получается из системы (1’) с помощью аналогичного преобразования, когда к первому уравнению системы (1’) прибавляется второе уравнение системы (1’), умноженное на число (- ).

Точно также доказывается и следующее утверждение.

Теорема 2’. Умножение произвольного уравнения системы (1) на любое число, отличное от нуля, переводит систему (1) в эквивалентную ей систему уравнений.

Теоремы 2 и 2’ дают два вида преобразований, которым подвергалась система (1), оставаясь эквивалентной:

а) умножение (или деление) произвольного уравнения системы (1) на любое число, отличное от нуля;

б) прибавление (или вычитание) к одному уравнению другого, умноженного на некоторое число.

Такие преобразования а) и б) называются элементарными преобразованиями системы уравнений (1).

Если к системе уравнений (1) несколько раз применить элементарные преобразования, то полученная в результате система, очевидно, также будет эквивалентна первоначальной.

Систему уравнений (1) можно записать в табличной форме:

x₁	x₂		x_j		x_n	b
a₁₁	a₁₂	…	a_1j	…	a_1n	b₁
a₂₁	a₂₂	…	a_2j	…	a_2n	b₂
…	…	…	…	…	…	…
a_i1	a_i2	…	a_ij	…	a_in	b_i
…	…	…	…	…	…	…
a_m1	a_m2	…	a_mj	…	a_mn	b_m

(2)

Прямоугольная таблица чисел, составленная из коэффициентов a_ij при неизвестных системы (1), называется матрицей системы (1) и обозначается A (в ней m строк и n столбцов), столбец свободных членов обозначается b. Прямоугольная таблица, составленная из коэффициентов a_ij при неизвестных и из столбца свободных членов b системы (1), называется расширенной матрицей системы (1) и обозначается (в ней m строк и (n+1) столбцов), т.е = (A, b). В i – ой строке матрицы содержатся все известные параметры, характеризующие i - ое уравнение системы (1), i = 1,…, m. В j – м столбце матрицы A содержатся все коэффициенты при неизвестном x_j, встречающиеся в системе (1).

Числа a_ij называются элементами матрицы А. Элемент a_ij находится в i - ой строкеи в j - м столбце матрицы А. Принято говорить, что элементa_ij находится на пересечении i - ой строки и j - го столбца матрицы А. Если все элементы строки (столбца) матрицы А (кроме одного) равны нулю, а ненулевой элемент равен единице, то такая строка (столбец) называется единичной (единичным).

Элементарным преобразованиям системы (1) соответствуют следующие элементарные преобразования таблицы (2):

а) умножение (или деление) всех элементов произвольной строки таблицы (2) на любое число, отличное от нуля,

б) прибавление (или вычитание) к одной строке (поэлементно) другой строки, умноженной на некоторое число.

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

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

Метод последовательного полного исключения неизвестных

(Метод Гаусса - Жордана)

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

Отметим принципиальное отличие предлагаемого метода решения систем линейных алгебраических уравнений от метода решения, скажем, стандартного квадратного уравнения. Оно решается с помощью хорошо известных формул, в которых неизвестные выражаются через коэффициенты уравнения. В случае общих систем линейных алгебраических уравнений мы таких формул не имеем и используем для отыскания решения метод итераций, или итеративный метод, или итерационный метод. Такие методы задают не формулы, а последовательность действий.

Метод Гаусса - Жордана представляет собой последовательную реализацию ряда однотипных больших шагов (или итераций). Это конкретный итерационный метод - один из многих методов итераций, предложенных для решения систем линейных алгебраических уравнений вида (1). Он состоит из начального этапа, основного этапа и заключительного этапа. Основной этап содержит повторяющиеся итерации – наборы однотипных действий.

Пусть задана конкретная система линейных алгебраических уравнений (1). Это значит, что известны n, m, a_ij, b_i, i = 1,…,m; j = 1,…,n. Опишем предлагаемый метод решения этой системы.

Начальный этап включает в себя построение таблицы I⁽⁰⁾ вида (2) и выбор в ней ведущего элемента – любого ненулевого коэффициента при переменных из таблицы (2). Столбец и строка, на пересечении которых стоит ведущий элемент, называются ведущими. (Пусть выбран элемент a_i₀_j₀. Тогда i₀– ая строка ведущая, j₀- й столбец ведущий.) Переходим к основному этапу. Заметим, что часто ведущий элемент называют разрешающим.

Основной этап состоит из повторяющихся однотипных итераций с номерами k = 1, 2,…. Опишем подробно итерации метода Гаусса - Жордана.

К началу каждой итерации известна некоторая таблица I вида (2), в ней выбран ведущий (разрешающий) элемент и, соответственно, ведущий столбец и ведущая строка. Кроме того, имеется информация о том, какие строки и столбцы уже были ведущими. (Так, например, после начального этапа, т.е. на итерации 1 известны I⁽⁰⁾, ведущий (разрешающий) элемент a_i₀_j₀и i₀– ая строка ведущая, j₀- ой столбец ведущий.)

Итерация(с номером k) состоит из следующих действий.

1. Преобразование ведущего столбца (т.е. столбца, содержащего ведущий элемент) в единичный с 1 на месте ведущего элемента путем последовательного поэлементного вычитания ведущей строки (т.е. строки, содержащей ведущий элемент), умноженной на некоторые числа, из остальных строк таблицы. Сама ведущая строка преобразуется путем поэлементного деления ее на ведущий элемент.

2. Выписывается новая таблица I⁽^k⁾, (k - номер итерации), в которой все столбцы, которые были когда-либо ведущими, – единичные.

3. Проверяется, можно ли в таблице I⁽^k⁾ выбрать новый ведущий (разрешающий) элемент. По определению это любой ненулевой элемент, который стоит на пересечении строки и столбца, которые еще не были ведущими.

Если такой выбор возможен, то столбец и строка, на пересечении которых стоит ведущий (разрешающий) элемент, называются ведущими. Затем итерация повторяется с новой таблицей I⁽^k⁾, т.е. действия 1 – 3 повторяются с новой таблицей I⁽^k⁾. При этом строится новая таблица I⁽^k⁺¹⁾.

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

Заключительный этап. Пусть проделано r итераций, получена таблица I⁽^r⁾, состоящая из матрицы коэффициентов при переменных A⁽^r⁾ и столбца свободных членов b⁽^r⁾, и в ней нельзя выбрать новый ведущий элемент, т.е. метод остановился. Заметим, что метод обязательн о остановится за конечное число шагов, т.к. r не может быть больше min{m,n}.

Каковы варианты остановки метода? Что значит «нельзя выбрать новый ведущий элемент»? Это значит, что после r – ой итерации в матрице A⁽^r⁾ новой системы, эквивалентной системе (1), либо

а) все строки A⁽^r⁾ были ведущими, т.е. в каждой строке стоит одна и ровно одна единица, которая не стоит больше не в какой другой строке,

либо

б) остались строки в A⁽^r⁾, состоящие только из нулей.

Рассмотрим эти варианты.

а) В этом случае r = m, m n. Переставив строки и перенумеровав переменные (т.е. переставив столбцы), можно таблицу I⁽^r⁾ представить в виде

x₁	x₂		x_r	x_r+1		x_n	b
		…		a⁽^r⁾_1,r+1	…	a⁽^r⁾_1n	b⁽^r⁾₁
		…		a⁽^r⁾_2,r+1	…	a⁽^r⁾_2n	b⁽^r⁾₂
…	…	…	…	…	…	…	…
		…		a⁽^r⁾_i,r+1	…	a⁽^r⁾_in	b⁽^r⁾_i
…	…	…	…	…	…	…	…
		…		a⁽^r⁾_m,r+1	…	a⁽^r⁾_mn	b⁽^r⁾_m

(3)

Подчеркнем, что в таблице (3) каждая переменная с номером i, не превосходящим r, встречается только в одной строке. Таблица (3) соответствует системе линейных уравнений вида

x₁+ = b⁽^r⁾₁,

x₂+ = b⁽^r⁾₂,

………………………, (4)

x_r + = b⁽^r⁾_r,

в которой каждая переменная с номером i, не превосходящим r, однозначно выражается через переменные x_r₊₁, …,x_n, коэффициенты матрицы a⁽^r⁾_ij, j = r+1,…,n, и свободный член b⁽^r⁾_i, представленные в таблице (3). На переменные x_r₊₁, …,x_n не накладываются никакие ограничения, т.е. они могут принимать любые значения. Отсюда произвольное решение системы, описываемой таблицей (3), или, что то же самое, произвольное решение системы (4), или, что то же самое, произвольное решение системы (1) имеет вид

x_i = b^(r)_i - a^(r)_ijx_j, i = 1,…,r = m; x_j – любое при j = (r+1),…,n. (5)

Тогда множество решений системы (1) можно записать как

X^b = {x=(x₁, …,x_n): x_i = b⁽^r⁾_i - a⁽^r⁾_ijx_j при i = 1,…, r = m; x_j – любое при j =(r+1),…,n.}.

б) В этом случае r < m, и существует хотя бы одна строка k, k > r, (предполагаем, что сделана перестановка строк и столбцов такая же, как в пункте а)) такая, что a⁽^r⁾_kj = 0 при всех j. Тогда, если соответствующий свободный член b⁽^r⁾_k не равен 0, то k - е уравнение не имеет решения, и, следовательно, вся система не имеет решения, т.е. система (1) несовместна.

Если же соответствующий b⁽^r⁾_k равен 0, то k - ое уравнение является лишним и его можно отбросить. Отбросив все такие уравнения, получим, что система (1) эквивалентна системе из r уравнений с n переменными, которая через r шагов записывается с помощью таблицы вида (3), в которой все строки были ведущими. Таким образом, мы пришли к рассмотренному выше случаю а) и можем выписать решение вида (5).

Метод Гаусса – Жордана описан полностью. За конечное число итераций система линейных алгебраических уравнений будет решена (если она совместна) или будет очевидно, что она несовместна (если она действительно несовместна).

Переменные, соответствующие ведущим (разрешающим) элементам, или стоящие в ведущих столбцах, принято называть базисными, а остальные переменные - свободными.

Обратим внимание на следующее.

1) Когдамы начинаем решать систему методом Гаусса - Жордана, мы можем не знать, совместна эта система или нет. Метод Гаусса - Жордана за конечное число итераций r даст ответ на этот вопрос. В случае совместной системы на основании последней таблицы выписывается общее решение исходной системы. В этом случае число базисных переменных обязательно равно номеру r последней итерации, т.е. числу выполненных итераций. Число r всегда не превосходит min{m,n},гдеm - число уравнений системы,а n - число переменных системы. Если r < n, то ( n – r) равно числу свободных переменных.

2) При записи общего решения не нужно перенумеровывать переменные, как это делалось для простоты понимания при описании Заключительного этапа. Это сделано для более ясного понимания.

3) При решении системы (1) методом Гаусса - Жордана базисными переменными будут только переменные, соответствующие столбцам, которые на каких-то итерациях выступали в роли ведущих, и наоборот, если на какой-то итерации столбец выступал в качестве ведущего, соответствующая ему переменная обязательно будет в числе базисных.

4) Если общее решение системы (1) содержит хотя бы одну свободную переменную, то эта система имеет бесконечно много частных решений, если же свободных переменных нет, то система имеет единственное решение, которое совпадает с общим решением.

5) Ведущие элементы могут быть выбраны на каждой итерации различным способом. Важно только то, что это ненулевые коэффициенты, стоящие на пересечении строки и столбца, которые до этого не были ведущими. Различный выбор ведущих элементов может дать различные записи множества решений. Однако, само множество решений при любой записи одно и то же.

Поясним работу метода на примерах.

Пример I. Решить следующую систему линейных алгебраических уравнений

2 x₁ – 3 x₂ + 3 x₃ + 5 x₄ = -1,

3 x₁ + 4 x₂ - 2 x₃ + 6 x₄ = 2, (6)

5 x₁ – 4 x₂ + 6 x₃ + 10 x₄ = 2

методом последовательного полного исключения неизвестных (методом Гаусса - Жордана).

Начальный этап. Сначала выпишем систему уравнений (6) в более удобной форме - в виде таблицы I⁽⁰⁾.

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
		-4

(0.1)

(0.2)

(0.3)

Первая строка таблицы I⁽⁰⁾ соответствует первому уравнению система (6), (будем ее, как и первое уравнение в (6), нумеровать (0.1)), вторая строка - второму уравнению системы (6) (будем ее нумеровать (0.2)), третья строка - третьему уравнению системы (6) (будем ее нумеровать (0.3)).

Первый столбец таблицы образуют коэффициенты системы (6) при неизвестной x₁, второй столбец - коэффициенты системы (6) при неизвестной x₂, аналогично поясняются третий и четвертый столбцы таблицы I⁽⁰⁾. Пятый столбец (b) таблицы составляют правые части уравнений (свободные члены) системы (6).

В качестве ведущего (разрешающего) элемента может быть выбран любой ненулевой элемент из первых четырех столбцов. Выберем в качестве ведущего элемент 2, стоящий на пересечении 1-ой строки и первого столбца. Тогда ведущая строка – первая, ведущий столбец - первый. Перейдем к основному этапу.

Основной этап. Известны таблица I⁽⁰⁾, ведущий элемент 2, стоящий на пересечении первой строки и первого столбца, ведущая строка – первая, ведущий столбец - первый.

Итерация 1. Подробно опишем преобразования таблицы I⁽⁰⁾в новую таблицу I⁽¹⁾, в которой первый столбец – единичный с 1 на месте ведущего элемента. Записывать новую таблицу будем под таблицей I⁽⁰⁾.

Сначала опишем преобразование, в результате которого элемент, расположенный непосредственно под ведущим элементом и равный 3, перейдет в нулевой элемент. Для этого нужно из строки (0.2) вычесть (поэлементно) ведущую строку (0.1), умноженную на 3/2. Тогда строка (0.2) заменится на строку

17/2

- 13/2

- 3/2

7/2

Поскольку в данном случае все элементы второй строки можно умножить на число 2 и ничего с точки зрения множества решений не изменится, то более удобно совершить еще одно элементарное преобразование (умножить все элементы строки на 2) и в таблицу I⁽¹⁾ следует поместить в качестве второй такую строку

-13

-3

(1.2)

Очевидно, такое преобразование таблицы I⁽⁰⁾соответствует элементарным преобразованиям системы (6), в результате которых система (6) превращается в эквивалентную систему, описываемую уравнениями (0.1), (1.2) и (0.3). Итак, в результатепреобразования второй строки таблицы I⁽⁰⁾мы получили новую вторую строку (вторую строку новой таблицы I⁽¹⁾), в которой в ведущем столбце (в нашем случае в первом столбце) вместо числа 3 стоит число 0.

Преобразуем третью строку таблицы I⁽⁰⁾ так, чтобы вместо 5 стоял 0. Для этого нужно из строки (0.3) вычесть (поэлементно) ведущую строку (0.1), умноженную на 5/2. Получим

-3

-5

(1.3)

Из трех строк таблицы I⁽⁰⁾непреобразованной осталась только ведущая (первая) строка. Все элементы ведущей строки (0.1) делим на ведущий элемент (т.е. в нашем конкретном случае - на число 2) и полученную новую строку записываем на ее место (т.е. на место (1.1)) в таблице I⁽¹⁾.

-3/2

3/2

5/2

-1/2

(1.1)

Таблица I⁽¹⁾будет иметь вид (три последние строки, выделенные курсивом)

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
		-4
I⁽¹⁾	1	-3/2	3/2	5/2	-1/2
	0	17	-13	-3	7
	0	7	-3	-5	9

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

(1.3)

Закончив заполнение таблицы I⁽¹⁾, мы проверяем, можем ли мы в ней выбрать новый ведущий элемент. По определению это любой ненулевой элемент, стоящий на пересечении строки и столбца, которые еще не были ведущими. В нашей ситуации это любой ненулевой элемент, который не стоит в первой строке и в первом столбце.

Выберем в качестве ведущего в таблице I⁽¹⁾ элемент (-3), стоящий на пересечении 3-его столбца и 3 – ей строки. Тогда ведущая строка – третья, ведущий столбец - третий.

Прежде чем переходить к следующей (второй) итерации, сделаем ряд важных замечаний, которые можно повторять дословно или с необходимыми корректировками после остальных итераций:

1. Итерация имеет вполне определенный (свой) ведущий элемент и, следовательно, определенные (свои) ведущие столбец и строку. Рабочей таблицей первой итерации была таблица I⁽⁰⁾, итоговой таблицей первой итерации - таблица I⁽¹⁾.

2. В результате выполнения итерации ведущий столбец превращается в единичный (на месте ведущего элемента появляется единица (она выделена полужирным шрифтом), а на местах остальных элементов ведущего столбца - нули, тем самым в результате итерации мы исключили неизвестную из всех уравнений системы кроме одного. Поскольку на каждой итерации полностью исключается своя неизвестная (переменная), становится понятным название метода Гаусса-Жордана как метода последовательно полного исключения неизвестных.

3. Таблица I⁽¹⁾ соответствует системе уравнений, эквивалентной исходной системе (6), которой соответствует таблица I⁽⁰⁾, т.к. преобразования, которым мы подвергали строки таблицы I⁽⁰⁾, представляют собой элементарные преобразования исходной системы (6).

4. Наиболее трудоемкую часть итерации составляют преобразования табличных строк.

5. Первую строку таблицы I⁽¹⁾ на число 2 умножать нельзя, т.к. на месте элемента, который был ведущим, должна стоять единица.

Итерация 2. Опишем преобразования таблицы I⁽¹⁾в новую таблицу I⁽²⁾, в которой третий столбец (ведущий) будет единичным. Записывать новую таблицу будем под таблицей I⁽¹⁾.

Далее, как и на предыдущей (первой) итерации, будем подвергать все строки элементарным преобразованиям.

Сначала переведем в нулевой элемент 3/2 ведущего столбца таблицы I⁽¹⁾. Для этого ведущую строку (1.3) умножим на 1/2 и прибавим к строке (1.1). Получим первую строку (2.1) новой таблицы I⁽²⁾.

(2.1)

Преобразуем теперь вторую строку таблицы I⁽¹⁾ так, чтобы вместо (-13) стоял 0. Для этого нужно из строки (1.2) вычесть (поэлементно) ведущую строку (1.3), умноженную на 13/3. Получим

-40/3

56/3

-32

В данном случае естественно все элементы этой строки умножить на 3/8 и в таблицу I⁽²⁾ поместить в качестве второй следующую строку

-5

-12

(2.2)

Из трех строк таблицы I⁽¹⁾непреобразованной осталась только ведущая (третья) строка. Все элементы ведущей строки (1.3) делим на ведущий элемент (т.е. в нашем конкретном случае - на число -3) и полученную новую строку записываем на ее место (т.е. на место (2.3)) в таблице I⁽²⁾.

-7/3

5/3

-3

(2.3)

Таблица I⁽²⁾будет иметь вид (три последние строки, выделенные курсивом)

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
		-4
I⁽¹⁾		-3/2	3/2	5/2	-1/2
			-13	-3
			-3	-5
I⁽²⁾		2	0	0	4
	0	-5	0		-12
	0	-7/3		5/3	-3

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

(1.3)

(2.1)

(2.2)

(2.3)

Закончив заполнение таблицы I⁽²⁾, проверяем, можно ли мы выбрать новый ведущий элемент – ненулевой элемент, стоящий на пересечении строки и столбца, которые еще не были ведущими. В нашей ситуации это любой ненулевой элемент, который не стоит в первой и третьей строке и в первом и третьем столбце. Выберем в качестве ведущего в таблице I⁽²⁾ элемент (7), стоящий на пересечении 4 - ого столбца и 2 – ей строки. Тогда ведущая строка – вторая, ведущий столбец- четвертый.

Итерация 3. Преобразуем четвертый столбец (ведущий) в единичный и получим новую таблицу I⁽³⁾. Записываем новую таблицу под таблицей I⁽²⁾.

В первой строке уже стоит 0 в четвертом столбце. Поэтому строку (2.1), не меняя, переносим в таблицу I⁽³⁾ как (3.1). С помощью элементарного преобразования переведем теперь в нулевой элемент 5/3 ведущего столбца таблицы I⁽²⁾. Для этого ведущую строку (2.2) умножим на 5/21 и прибавим к строке (2.3). Получим третью строку (3.3) новой таблицы I⁽³⁾:

-8/7

-1/7

(3.3)

Все элементы ведущей строки (2.2) делим на ведущий элемент (т.е. в нашем конкретном случае - на число 7) и полученную новую строку

-5/7

-12/7

(3.2)

записываем на ее место (т.е. на место (3.2)) в таблице I⁽³⁾.

Таблица I⁽³⁾будет иметь вид (три последние строки, выделенные курсивом)

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
		-4
I⁽¹⁾		-3/2	3/2	5/2	-1/2
			-13	-3
			-3	-5
I⁽²⁾
		-5			-12
		-7/3		5/3	-3
I⁽³⁾		2	0	0	4
	0	-5/7	0		-12/7
	0	-8/7		0	-1/7

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

(3.3)

(2.1)

(2.2)

(2.3)

(3.1)

(3.2)

(3.3)

Теперь проверяем, можем ли мы выбрать в I⁽³⁾ новый ведущий элемент – ненулевой элемент, стоящий на пересечении строки и столбца, которые еще не были ведущими. В нашей ситуации в таблице I⁽³⁾уже нельзя выбрать новый ведущий элемент, т.к. в ней всего три строки, и все они уже были ведущими на предыдущих итерациях. Следовательно, таблица I⁽³⁾(и третья итерация) является последней (конечной). Переходим к Заключительному этапу.

Заключительный этап. В таблице I⁽³⁾ имеется три единичных столбца, каждый из которых на одной из итераций был ведущим. Соответствующие переменные (в нашем случае это переменные x₁, x_3,x₄) являются базисными, а остальные переменные (в нашем случае это переменная x₂) - свободные.

Запишем систему линейных алгебраических уравнений, соответствующую последней таблице I⁽³⁾.

x₁ + 2 x₂ = 4

(-5/7) x₂ + x₄ =-12/7 (7)

(-8/7) x₂ + x₃ = - 1/7.

Система (7) эквивалентна системе (6), т.е. имеет то же самое множество решений. Выразим базисные переменные через свободную. С учетом того, что свободная переменная может принимать любое значение, получим

x₁ = 4 - 2 x₂, x₃ = -1/7 + (8/7) x₂, x₄ = -12/7 + (5/7) x₂, x₂ - любое. (8)

Любой набор x = (x₁, x₂, x₃, x₄), удовлетворяющий (8), представляет собой решение системы (7)и, следовательно, эквивалентной ей исходной системы (6). Множество решений X^b = { x = (x₁, x₂, x₃, x₄): x₁ = 4 - 2x₂, x₃ = -1/7 + (8/7)x₂, x₄ = -12/7 + (5/7 x₂, x₂ – любое}.

Придавая свободной переменной всевозможные значения, будем из набора (8) получать различные частные решения системы (6). Например, полагая x₂ = 0, получим частное решение (4, 0, - 1/7, - 12/7). Общее решение (8) системы (6) содержит (описывает единой формулой) вс е ее частные решения.

Рассмотренный пример показал один возможный вариант окончания метода, при котором система уравнений оказалась совместной, и было выписано общее решение. Приведем примеры, демонстрирующие другие возможные варианты остановки метода. При этом покажем, как метод Гаусса – Жордана применяется для одновременного решения нескольких систем, отличающихся только свободными членами (правыми частями).

Пример 2. Решить две следующие системы, которые отличаются одна от другой только правыми частями

2 x₁ – 3 x₂ + 3 x₃ + 5 x₄ = -1,

3 x₁ + 4 x₂ - 2 x₃ + 6 x₄ = 2, (9)

4 x₁ + 6 x₂ - 6 x₃ - 10 x₄ = 3

2 x₁ – 3 x₂ + 3 x₃ + 5 x₄ = -1,

3 x₁ + 4 x₂ - 2 x₃ + 6 x₄ = 2, (10)

4 x₁ + 6 x₂ - 6 x₃ - 10 x₄ = 2.

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

Начальный этап. Перепишем две системы в виде таблицы I⁽⁰⁾, в которой будут фигурировать два столбца правых частей b⁽⁹⁾ и b⁽¹⁰⁾ (а не один, как в предыдущем примере).

	x₁	x₂	x₃	x₄	b⁽⁹⁾	b⁽¹⁰⁾
I⁽⁰⁾		-3			-1	-1
			-2
	-4		-6	-10

(0.1)

(0.2)

(0.3)

Выберем, как и в примере 1, в качестве ведущего элемент 2, стоящий на пересечении 1-ой строки и 1 - го столбца. Тогда ведущая строка – первая, ведущий столбец - первый. Перейдем к основному этапу.

Итерация 1. Построим новую таблицу I⁽¹⁾, в которой первый столбец – единичный с 1 на месте ведущего элемента. Записывать новую таблицу будем под таблицей I⁽⁰⁾. Очевидно, что первая и вторая строка в ней совпадают с соответствующими строками в таблице I⁽¹⁾ примера 1. Для того, чтобы получить 0 на месте (-4) в третьей строке, умножим строку (0.1) на 2 и прибавим (почленно) к строке (0.3). Таблица I⁽¹⁾будет иметь вид (три последние строки, выделенные курсивом)

	x₁	x₂	x₃	x₄	B⁽⁹⁾	b⁽¹⁰⁾
I⁽⁰⁾		-3			-1	-1
			-2
	-4		-6	-10
I⁽¹⁾		-3/2	3/2	5/2	-1/2	-1/2
	0	17	-13	-3	7	7
	0	0	0	0	1	0

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

(1.3)

Таблица I⁽¹⁾ с первым столбцом свободных членов соответствует системе уравнений, эквивалентной системе (9), а с вторым столбцом свободных членов - системе уравнений, эквивалентной системе (10). Для системы (9) последнее уравнение имеет вид

0 x₁ + 0 x₂ + 0 x₃ + 0 x₄ = 1,

что показывает несовместность системы (9).

Для системы (10) последнее уравнение имеет вид

0 x₁ + 0 x₂ + 0 x₃ + 0 x₄ = 0.

Отсюда следует, что его можно отбросить, и что система (10) эквивалентна системе из 2 – х уравнений с 4 – мя неизвестными, представляемой строками (1.1), (1.2) и вторым столбцом свободных членов в таблице I⁽¹⁾. Мы можем продолжать решать систему (10), используя следующую таблицу I⁽¹⁾

I⁽¹⁾		-3/2	3/2	5/2	-1/2
	0	17	-13	-3	7

(1.1)

(1.2)

или

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
	-4		-6	-10
I⁽¹⁾		-3/2	3/2	5/2	-1/2
	0	17	-13	-3	7

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

В таблице I⁽¹⁾ можно выбрать в качестве ведущего любой ненулевой элемент, стоящий во второй строке и в любом, кроме первого, столбце. Выберем в качестве ведущего в таблице I⁽¹⁾ элемент (-3), стоящий на пересечении 4 - ого столбца и 2 – ой строки. Тогда ведущая строка – вторая, ведущий столбец- четвертый.

Итерация 2. Преобразуем таблицу I⁽¹⁾в новую таблицу I⁽²⁾, в которой четвертый столбец (ведущий) будет единичным. Для этого преобразуем первую строку таблицы I⁽¹⁾ так, чтобы вместо 5/2 стоял 0. Для этого нужно из строки (1.1) вычесть (поэлементно) ведущую строку (1.2), умноженную на (-5/6). Получим

38/3

-28/3

16/3

(2.1)

Все элементы ведущей строки (1.2) делим на ведущий элемент (т.е. в нашем конкретном случае - на число (-3)) и полученную новую строку записываем на ее место (т.е. на место (2.2)) в таблице I⁽²⁾.

Таблица I⁽²⁾будет иметь вид (две последние строки, выделенные курсивом)

	x₁	x₂	x₃	x₄	b
I⁽⁰⁾		-3			-1
			-2
	-4		-6	-10
I⁽¹⁾		-3/2	3/2	5/2	-1/2
			-13	-3
I⁽²⁾	1	38/3	-28/3	0	16/3
	0	-17/3	13/3	1	-7/3

(0.1)

(0.2)

(0.3)

(1.1)

(1.2)

(2.1)

(2.2)

Обе строки уже были ведущими, поэтому мы не можем в таблице I⁽²⁾ выбрать ведущий элемент. Переходим к Заключительному этапу.

Заключительный этап. Все строки были ведущими, поэтому система, определяемая таблицей I⁽²⁾, и, следовательно, система (10) разрешима. При этом переменные x₁ и x₄– базисные, а x₂и x₃– свободные. Запишем систему линейных алгебраических уравнений, соответствующую последней таблице I⁽²⁾.

x₁ + (38/3) x₂ - (28/3) x₃= 16/3

x₄ - (17/3) x₂+ (13/3) x₃ = - 7/3 (11)

Из (11) можно однозначным образом выразить базисные переменные через свободные и правую часть. На свободные переменные при этом нет никаких ограничений, т.е. они могут принимать любые значения. Таким образом, множеством решений системы (11) и, следовательно, системы (10) является следующее множество X^b,

X^b = { x = (x₁, x₂, x₃, x₄): x₁ = 16/3 - (38/3) x₂ + (28/3) x₃, x₄= - 7/3 + (17/3) x₂

- (13/3) x₃, x₂и x₃– любые.}

Примеры 1 и 2 проиллюстрировали работу метода Гаусса – Жордана и все возможные способы остановки метода. Повторим основные выводы, которые можно сделать из изложенного.

1) Метод остановится через конечное число итераций r. При этом выяснится, является система совместной или несовместной.

2) Если система совместна, то число базисных переменных равно числу итераций r, а число свободных переменных – (n – r). Все базисные переменные однозначно выражаются через свободные переменные и соответствующие свободные члены.

3) Если система совместна и общее решение содержит хотя бы одну свободную переменную, то число решений бесконечно. Если же свободных переменных нет (n = r), то решение единственно.