Решение системы линейных уравнений методом гаусса-жордана - онлайн калькулятор. Метод Гаусса-Жордана. Как найти обратную матрицу с помощью элементарных преобразований

Каждой системе линейных уравнений поставим в соответствие расширенную матрицу , полученную присоединением к матрице А столбца свободных членов:

Метод Жордана–Гаусса применяется для решения системы m линейных уравнений с n неизвестными вида:

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

Над строками расширенной матрицы осуществляем следующие элементарные преобразования:

1. перестановка двух строк ;

2. умножение строки на любое число, отличное от нуля ;

3. прибавление к одной строке другой строки, умноженной на некоторое число ;

4. отбрасывание нулевой строки (столбца) .

Пример 2.11. Решить методом Жордана–Гаусса системы линейных уравнений:

а ) Х 1 + Х 2 + 2Х 3 = -1

2Х 1 - Х 2 + 2Х 3 = -4

4Х 1 + Х 2 + 4Х 3 = -2

Решение: Составим расширенную матрицу:

Итерация 1

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

На этом первая итерация закончена.

Итерация 2

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

Итерация 3

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

откуда Х 1 = 1, Х 2 = 2, Х 3 = -2.

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

б ) Х 1 – Х 2 + Х 3 – Х 4 = 4

Х 1 + Х 2 + 2Х 3 +3Х 4 = 8

2Х 1 +4Х 2 + 5Х 3 +10Х 4 = 20

2Х 1 – 4Х 2 + Х 3 – 6Х 4 = 4

Решение: Расширенная матрица имеет вид:

Применяя элементарные преобразования, получим:

Исходная система эквивалентна следующей системе уравнений:

Х 1 – 3Х 2 – 5Х 4 = 0

2Х 2 + Х 3 + 4Х 4 = 4

Последние две строки матрицы A (2) являются линейно зависимыми.

Определение. Строки матрицы e 1 , e 2 ,…, e m называются линейно зависимыми , если существуют такие числа , не равные одновременно нулю, что линейная комбинация строк матрицы равна нулевой строке:

где 0 =(0, 0…0). Строки матрицы являются линейно независимыми , когда комбинация этих строк равна нулю тогда и только тогда, когда все коэффициенты равны нулю.

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

Теорема 2.3 (о ранге матрицы). Ранг матрицы равен максимальному числу её линейно независимых строк или столбцов, через которые линейно выражаются все остальные её строки (столбцы).

Ранг матрицы A (2) равен 2, т.к. в ней максимальное число линейно независимых строк равно 2 (это первые две строки матрицы).

Теорема 2.4 (Кронекера–Капели). Система линейных уравнений совместна и только тогда, когда ранг матрицы системы равен рангу расширенной матрицы этой системы.

1. Если ранг матрицы совместной системы равен числу переменных, т.е. r = n, то система имеет единственное решение.

2. Если ранг матрицы системы меньше числа переменных, т.е. r < n, то система неопределённая и имеет бесконечное множество решений.

В данном случае система имеет 4 переменных, а её ранг равен 2, следовательно, она имеет бесконечное множество решений.

Определение. Пусть r < n , r переменных x 1 , x 2 ,…, x r называются базисными , если определитель матрицы из коэффициентов при них (базисный минор ) отличен от нуля. Остальные n – r переменных называются свободными .

Определение. Решение системы, в котором все n – r свободных переменных равны нулю, называется базисным .

Совместная система m линейных уравнений с n переменными (m < n ) имеет бесконечное множество решений, среди которых базисных решений конечное число, не превосходящее , где .

В нашем случае , т.е. система имеет не более 6 базисных решений.

Общее решение имеет вид:

Х 1 = 3Х 2 +5Х 4

Х 3 = 4 – 2Х 2 – 4Х 4

Найдем базисные решения. Для этого полагаем Х 2 = 0, Х 4 = 0, тогда Х 1 =0, Х 3 = 4. Базисное решение имеет вид: (0, 0, 4, 0).

Получим другое базисное решение. Для этого в качестве свободных неизвестных примем Х 3 и Х 4 . Выразим неизвестные Х 1 и Х 2 через неизвестные Х 3 и Х 4:

Х 1 = 6 – 3/2Х 2 – Х 4

Х 2 = 2 – 1/2Х 3 – 2Х 4 .

Тогда базисное решение имеет вид: (6, 2, 0, 0).

Пример 2.12. Решить систему:

X 1 + 2X 2 – X 3 = 7

2X 1 – 3X 2 + X 3 = 3

4X 1 + X 2 – X 3 = 16

Решение.Преобразуем расширенную матрицу системы

Итак, уравнение, соответствующее третьей строке последней матрицы, противоречиво – оно привелось к неверному равенству 0 = –1, следовательно, данная система несовместна. Данный вывод можно также получить, если заметить, что ранг матрицы системы равен 2, тогда как ранг расширенной матрицы системы равен 3.

метод Гаусса–Жордана - один из наиболее известных и широко применяемых методов решения систем линейных уравнений. Матричный метод и метод Крамера обладают тем недостатком, что они не дают ответа в том случае, когда detA = 0, а определяют лишь единственное решение при detA неравном 0. Еще одним недостатком является то, что объем математических вычислений в рамках этих методов резко возрастает с ростом числа уравнений. Метод Гаусса практически свободен от этих недостатков.

Алгоритм метода Гаусса

1. На основании системы линейных уравнений составляем расширенную матрицу системы;
2. Приводим матрицу к "треугольному" виду;
3. Определяем ранги основной и расширенной матриц, и на основании этого делаем вывод о совместности системы и количестве допустимых решений;
4. В случае, если система имеет единственное решение производим обратную подстановку и находим его, если система имеет множество решений: выражаем базисные переменные через переменные которые могут принимать произвольные значения;

Комментарий к шагу 2 Метода Гаусса. Треугольной называют матрицу, в которой все элементы расположенные ниже главной диагонали равны нулю.

Для приведения исходной расширенной матрицы к треугольному виду используем следующие два свойства определителей:

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

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

На основании этих свойств определителей составим алгоритм преобразования матрицы к треугольному виду:

1. Рассматриваем строку i(начиная с первой). Если, элемент a i i равен нулю, меняем местами i-ю и i+1-ю строки матрицы. Знак определителя при этом изменится на противоположный. Если a 1 1 отличен от нуля - переходим к следующему шагу;
2. Для каждой строки j, ниже i-й находим значение коэффициента K j =a j i /a i i ;
3. Пересчитываем элементы всех строк j, расположенных ниже текущей строки i, с использованием соответствующих коэффициентов по формуле: a j k нов.=a j k -K j *a i k ; После чего, возвращаемся к первому шагу алгоритма и рассматриваем следующую строку, пока не доберемся до строки i=n-1, где n - размерность матрицы A
4. В полученной треугольной матрице расчитываем произведение всех элементов главной диагонали Пa i i , которое и будет являтся определителем;

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

4. Метод Жордана - Гаусса.

Схема с выбором главного элемента состоит в том, что требование неравенства нулю диагональных элементов akk, на которые происходит деление в процессе исключения, заменятся более жестким: из всех элементов К-го столба выбрать наибольший по модулю и переставить уравнения так, чтобы этот элемент оказался на месте элемента акк. Выбор главного элемента и связанная с ним перестановка строк необходимы в тех случаях, когда на каком-либо i-ом шаге акк=0 либо же акк очень мало по остальными элементами i- го столбца: при делении на такое «малое» акк будут получаться большие числа с большими абсолютными погрешностями, в результате чего решение может сильно исказиться.

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

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

1. В процессе исключений левая часть I –го уравнения системы обращается в нуль, а правая часть равна некоторому числу, отличному от нуля. т.е. 02+=bc0.

Это означает, что система не имеет решений, так как I – му уравнению не могут удовлетворять никакие значения неизвестных;

2. Левая и правая части I – го уравнения обращаются в нуль. Это означает, что I – ое уравнение является линейной комбинацией остальных, ему удовлетворяет любое найденное решение системы, поэтому оно может быть отброшено. В системе количество неизвестных больше количества уравнений и, следовательно, такая система имеет множество решений;

3. После того как все уравнения использованы для исключения неизвестных получено решение системы.

Таким образом, конечной целью преобразований Жордана-Гаусса является получение из заданной линейной системы

a11x1 + a12x2 + … + a1nxn = b1,n+1

a21x1 + a22x2 + … + a2nxn = b2,n+1

am1x1 + am2x2 + … + amnxn = bm.n+1

Здесь x1, x2, …, xn - неизвестные, которые надо определить. a11, a12, …, amn - коэффициенты системы - и b1, b2, … bm - свободные члены - предполагаются известными. Индексы коэффициентов (aij) системы обозначают номера уравнения (i) и неизвестного (j), при котором стоит этот коэффициент, соответственно.

Система (1) называется однородной, если все её свободные члены равны нулю (b1 = b2 = … = bm = 0), иначе - неоднородной.

Система (1) называется квадратной, если число m уравнений равно числу n неизвестных.

Решение системы (1) - совокупность n чисел c1, c2, …, cn, таких что подстановка каждого ci вместо xi в систему (1) обращает все ее уравнения в тождества.

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

Совместная система вида (1) может иметь одно или более решений.

Решения c1(1), c2(1), …, cn(1) и c1(2), c2(2), …, cn(2) совместной системы вида (1) называются различными, если нарушается хотя бы одно из равенств:

c1(1) = c1(2), c2(1) = c2(2), …, cn(1) = cn(2).

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

Решим следующую систему уравнений:

Запишем её в виде матрицы 3×4, где последний столбец является свободным членом:

Проведём следующие действия:

· К строке 2 добавим: -4 * Строку 1.

· К строке 3 добавим: -9 * Строку 1.

· К строке 3 добавим: -3 * Строку 2.

· Строку 2 делим на -2

· К строке 1 добавим: -1 * Строку 3.

· К строке 2 добавим: -3/2 * Строку 3.

· К строке 1 добавим: -1 * Строку 2.

В правом столбце получаем решение:

.

В методе Ньютона наблюдается ускорение сходимости процесса приближений. 5. Метод касательных (метод Ньютона) Метод касательных, связанный с именем И. Ньютона, является одним из наиболее эффективных численных методов решения уравнений. Идея метода очень проста. Возьмём производную точку x0 и запишем в ней уравнение касательной к графику функции f(x): y=f(x0)+ f ¢(x) (x-x0) (1.5) Графики...

Решения от численных методов расчёта. Для определения корней уравнения не требуется знания теорий групп Абеля, Галуа, Ли и пр. и применения специальной математической терминологии: колец, полей, идеалов, изоморфизмов и т.д. Для решения алгебраического уравнения n - ой степени нужно только умение решать квадратные уравнения и извлекать корни из комплексного числа. Корни могут быть определены с...

Математики тригонометрической подстановки и проверка эффективности разработанной методики преподавания. Этапы работы: 1. Разработка факультативного курса на тему: «Применение тригонометрической подстановки для решения алгебраических задач» с учащимися классов с углубленным изучением математики. 2. Проведение разработанного факультативного курса. 3. Проведение диагностирующей контрольной...

... «проявляется» лишь в процессе преобразований. Очевидность и «завуалированность» новой переменной мы рассмотрим на конкретных примерах во второй главе данной работы. 2. Возможности применения метода замены неизвестного при решении алгебраических уравнений В этой главе выявим возможности применения метода замены неизвестного при решении алгебраических уравнений в стандартных и нестандартных...

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

Размер матрицы: 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 X 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101

О методе

При решении системы линейных уравнений онлайн методом Гаусса выполняются следующие шаги.

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

Чтобы лучше всего понять принцип работы алгоритма Гаусса онлайн введите любой пример, выберите "очень подробное решение" и посмотрите его решение онлайн.

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

Тема 7

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

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

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

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

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

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

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

a 1 x 1 + a 2 x 2 + … a j x j + … a n x 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 1 , a i 2 , …, a in и свободным членом b i , т.е. имеет вид

a i1 x 1 + a i2 x 2 + … + a ij x j + … + a in x n = b i .

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

a 11 x 1 + a 12 x 2 + … + a 1j x j + … + a 1n x n = b 1

a 21 x 1 + a 22 x 2 + … + a 2j x j + … + a 2n x n = b 2

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

a i1 x 1 + a i2 x 2 + … + a ij x j + … + a in x n = b i (1)

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

a m1 x 1 + a m2 x 2 + … + a mj x j + … + a mn x n = b m

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

= b i , i = 1,…, m .

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

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

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

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

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

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

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

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

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

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

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

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

0 x 1 + 0 x 2 + … + 0 x j + … + 0 x n = b k .

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

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

= b 1 + b 2 .

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

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

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

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

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

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

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

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

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

Прямоугольная таблица чисел, составленная из коэффициентов 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 (0) вида (2) и выбор в ней ведущего элемента – любого ненулевого коэффициента при переменных из таблицы (2). Столбец и строка, на пересечении которых стоит ведущий элемент, называются ведущими . (Пусть выбран элемент a i 0 j 0 . Тогда i 0 – ая строка ведущая, j 0 - й столбец ведущий.) Переходим к основному этапу. Заметим, что часто ведущий элемент называют разрешающим .

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

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

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

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

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

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

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

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

Заключительный этап. Пусть проделано 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) представить в виде

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

x 1 +
=b (r) 1 ,

x 2 +
=b (r) 2 ,

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

x r +
=b (r) r ,

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

x i = b (r) i - a (r) ij x j , i = 1,…,r = m; x j – любое при j = (r+1),…,n. (5)

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

X b = {x=(x 1 , … ,x n) : x i = b (r) i - a (r) ij x 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 1 – 3 x 2 + 3 x 3 + 5 x 4 = -1,

3 x 1 + 4 x 2 - 2 x 3 + 6 x 4 = 2, (6)

5 x 1 – 4 x 2 + 6 x 3 + 10 x 4 = 2

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

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