Решение систем линейных уравнений методом Гаусса.
Пусть нам требуется найти решение системы из n
линейных уравнений с n
неизвестными переменными
определитель основной матрицы которой отличен от нуля.
Суть метода Гаусса состоит в последовательном исключении неизвестных переменных: сначала исключается x 1 из всех уравнений системы, начиная со второго, далее исключается x 2 из всех уравнений, начиная с третьего, и так далее, пока в последнем уравнении останется только неизвестная переменная x n . Такой процесс преобразования уравнений системы для последовательного исключения неизвестных переменных называется прямым ходом метода Гаусса . После завершения прямого хода метода Гаусса из последнего уравнения находитсяx n , с помощью этого значения из предпоследнего уравнения вычисляется x n-1 , и так далее, из первого уравнения находится x 1 . Процесс вычисления неизвестных переменных при движении от последнего уравнения системы к первому называется обратным ходом метода Гаусса .
Кратко опишем алгоритм исключения неизвестных переменных.
Будем считать, что , так как мы всегда можем этого добиться перестановкой местами уравнений системы. Исключим неизвестную переменную x 1
из всех уравнений системы, начиная со второго. Для этого ко второму уравнению системы прибавим первое, умноженное на , к третьему уравнению прибавим первое, умноженное на , и так далее, к n-ому
уравнению прибавим первое, умноженное на . Система уравнений после таких преобразований примет вид
где , а .
К такому же результату мы бы пришли, если бы выразили x 1 через другие неизвестные переменные в первом уравнении системы и полученное выражение подставили во все остальные уравнения. Таким образом, переменная x 1 исключена из всех уравнений, начиная со второго.
Далее действуем аналогично, но лишь с частью полученной системы, которая отмечена на рисунке
Для этого к третьему уравнению системы прибавим второе, умноженное на , к четвертому уравнению прибавим второе, умноженное на , и так далее, к n-ому
уравнению прибавим второе, умноженное на . Система уравнений после таких преобразований примет вид
где , а . Таким образом, переменная x 2
исключена из всех уравнений, начиная с третьего.
Далее приступаем к исключению неизвестной x 3
, при этом действуем аналогично с отмеченной на рисунке частью системы
Так продолжаем прямой ход метода Гаусса пока система не примет вид
С этого момента начинаем обратный ход метода Гаусса: вычисляем x n из последнего уравнения как , с помощью полученного значения x n находим x n-1 из предпоследнего уравнения, и так далее, находим x 1 из первого уравнения.
Пример.
Решите систему линейных уравнений методом Гаусса.
Сегодня разбираемся с методом Гаусса для решения систем линейных алгебраических уравнений. О том, что это за системы, можно почитать в предыдущей статье, посвященной решению тех же СЛАУ методом Крамера. Метод Гаусса не требует каких-то специфических знаний, нужна лишь внимательность и последовательность. Несмотря на то что с точки зрения математики для его применения хватит и школьной подготовки, у студентов освоение этого метода часто вызывает сложности. В этой статье попробуем свести их на нет!
Метод Гаусса
Метод Гаусса – наиболее универсальный метод решения СЛАУ (за исключением ну уж очень больших систем). В отличие от рассмотренного ранее , он подходит не только для систем, имеющих единственное решение, но и для систем, у которых решений бесконечное множество. Здесь возможны три варианта.
- Система имеет единственное решение (определитель главной матрицы системы не равен нулю);
- Система имеет бесконечное множество решений;
- Решений нет, система несовместна.
Итак, у нас есть система (пусть у нее будет одно решение), и мы собираемся решать ее методом Гаусса. Как это работает?
Метод Гаусса состоит из двух этапов – прямого и обратного.
Прямой ход метода Гаусса
Сначала запишем расширенную матрицу системы. Для этого в главную матрицу добавляем столбец свободных членов.
Вся суть метода Гаусса заключается в том, чтобы путем элементарных преобразований привести данную матрицу к ступенчатому (или как еще говорят треугольному) виду. В таком виде под (или над) главной диагональю матрицы должны быть одни нули.
Что можно делать:
- Можно переставлять строки матрицы местами;
- Если в матрице есть одинаковые (или пропорциональные) строки, можно удалить их все, кроме одной;
- Можно умножать или делить строку на любое число (кроме нуля);
- Нулевые строки удаляются;
- Можно прибавлять к строке строку, умноженную на число, отличное от нуля.
Обратный ход метода Гаусса
После того как мы преобразуем систему таким образом, одна неизвестная Xn становится известна, и можно в обратном порядке найти все оставшиеся неизвестные, подставляя уже известные иксы в уравнения системы, вплоть до первого.
Когда интернет всегда под рукой, можно решить систему уравнений методом Гаусса онлайн . Достаточно лишь вбить в онлайн-калькулятор коэффициенты. Но согласитесь, гораздо приятнее осознавать, что пример решен не компьютерной программой, а Вашим собственным мозгом.
Пример решения системы уравнений методом Гаусс
А теперь - пример, чтобы все стало наглядно и понятно. Пусть дана система линейных уравнений, и нужно решить ее методом Гаусса:
Сначала запишем расширенную матрицу:
Теперь займемся преобразованиями. Помним, что нам нужно добиться треугольного вида матрицы. Умножим 1-ую строку на (3). Умножим 2-ую строку на (-1). Добавим 2-ую строку к 1-ой и получим:
Затем умножим 3-ую строку на (-1). Добавим 3-ую строку к 2-ой:
Умножим 1-ую строку на (6). Умножим 2-ую строку на (13). Добавим 2-ую строку к 1-ой:
Вуаля - система приведена к соответствующему виду. Осталось найти неизвестные:
Система в данном примере имеет единственное решение. Решение систем с бесконечным множеством решений мы рассмотрим в отдельной статье. Возможно, сначала Вы не будете знать, с чего начать преобразования матрицы, но после соответствующей практики набъете руку и будете щелкать СЛАУ методом Гаусса как орешки. А если Вы вдруг столкнетесь со СЛАУ, которая окажется слишком крепким орешком, обращайтесь к нашим авторам! вы можете, оставив заявку в Заочнике. Вместе мы решим любую задачу!
Одним из универсальных и эффективных методов решения линейных алгебраических систем является метод Гаусса , состоящий в последовательном исключении неизвестных.
Напомним, две системы называются эквивалентными (равносильными), если множества их решений совпадают. Другими словами, системы эквивалентны, если каждое решение одной из них является решением другой и наоборот. Эквивалентные системы получаются приэлементарных преобразованиях уравнений системы:
умножение обеих частей уравнения на число отличное от нуля;
прибавление к некоторому уравнению соответствующих частей другого уравнения, умноженных на число отличное от нуля;
перестановка двух уравнений.
Пусть дана система уравнений
Процесс решения этой системы по методу Гаусса состоит из двух этапов. На первом этапе (прямой ход) система с помощью элементарных преобразований приводится к ступенчатому , илитреугольному виду, а на втором этапе (обратный ход) идет последовательное, начиная с последнего по номеру переменного, определение неизвестных из полученной ступенчатой системы.
Предположим, что коэффициент данной
системы
,
в противном случае в системе первую
строку можно поменять местами с любой
другой строкой так, чтобы коэффициент
прибыл отличен от нуля.
Преобразуем систему, исключив неизвестное во всех уравнениях, кроме первого. Для этого умножим обе части первого уравнения наи сложим почленно со вторым уравнением системы. Затем умножим обе части первого уравнения наи сложим с третьим уравнением системы. Продолжая этот процесс, получим эквивалентную систему
Здесь
– новые значения коэффициентов и
свободных членов, которые получаются
после первого шага.
Аналогичным образом, считая главным
элементом
,
исключим неизвестноеиз всех уравнений системы, кроме первого
и второго. Продолжим этот процесс, пока
это возможно, в результате получим
ступенчатую систему
,
где
,
,…,–
главные элементы системы
.
Если в процессе приведения системы к
ступенчатому виду появятся уравнения
,
т. е. равенства вида
,
их отбрасывают, так как им удовлетворяют
любые наборы чисел
.
Если же при
появится уравнение вида,
которое не имеет решений, то это
свидетельствует о несовместности
системы.
При обратном ходе из последнего уравнения
преобразованной ступенчатой системы
выражается первое неизвестное
через все остальные неизвестные
,
которые называютсвободными
.
Затем выражение переменнойиз последнего уравнения системы
подставляется в предпоследнее уравнение
и из него выражается переменная
.
Аналогичным образом последовательно
определяются переменные
.
Переменные
,
выраженные через свободные переменные,
называютсябазисными
(зависимыми).
В результате получается общее решение
системы линейных уравнений.
Чтобы найти частное решение
системы, свободным неизвестным
в общем решении придаются произвольные
значения и вычисляются значения
переменных
.
Технически удобнее подвергать элементарным преобразованиям не сами уравнения системы, а расширенную матрицу системы
.
Метод Гаусса - универсальный метод,
который позволяет решать не только
квадратные, но и прямоугольные системы,
в которых число неизвестных
не равно числу уравнений
.
Достоинство этого метода состоит
также в том, что в процессе решения мы
одновременно исследуем систему на
совместность, так как, приведя расширенную
матрицу
к ступенчатому виду, легко определить
ранги матрицыи расширенной матрицы
и применитьтеорему Кронекера -
Капелли
.
Пример 2.1 Методом Гаусса решить систему
Решение
. Число уравнений
и число неизвестных
.
Составим расширенную матрицу системы, приписав справа от матрицы коэффициентов столбец свободных членов.
Приведём матрицу к треугольному виду; для этого будем получать «0» ниже элементов, стоящих на главной диагонали с помощью элементарных преобразований.
Чтобы получить «0» во второй позиции первого столбца, умножим первую строку на (-1) и прибавим ко второй строке.
Это преобразование запишем числом (-1) против первой строки и обозначим стрелкой, идущей от первой строки ко второй строке.
Для получения «0» в третьей позиции первого столбца, умножим первую строку на (-3) и прибавим к третьей строке; покажем это действие с помощью стрелки, идущей от первой строки к третьей.
.
В полученной матрице, записанной второй в цепочке матриц, получим «0» во втором столбце в третьей позиции. Для этого умножили вторую строку на (-4) и прибавили к третьей. В полученной матрице вторую строку умножим на (-1), а третью - разделим на (-8). Все элементы этой матрицы, лежащие ниже диагональных элементов - нули.
Так как , система является совместной и определенной.
Соответствующая последней матрице система уравнений имеет треугольный вид:
Из последнего (третьего) уравнения
.
Подставим во второе уравнение и получим
.
Подставим
и
в первое уравнение, найдём
.
Данный онлайн калькулятор находит решение системы линейных уравнений (СЛУ) методом Гаусса. Дается подробное решение. Для вычисления выбирайте количество переменных и количество уравнений. Затем введите данные в ячейки и нажимайте на кнопку "Вычислить."
|
Представление чисел:
Целые числа и (или) Обыкновенные дробиЦелые числа и (или) Десятичные дроби
Число знаков после десятичного разделителя
×
Предупреждение
Очистить все ячейки?
Закрыть Очистить
Инструкция ввода данных. Числа вводятся в виде целых чисел (примеры: 487, 5, -7623 и т.д.), десятичных чисел (напр. 67., 102.54 и т.д.) или дробей. Дробь нужно набирать в виде a/b, где a и b (b>0) целые или десятичные числа. Примеры 45/5, 6.6/76.4, -7/6.7 и т.д.
Метод Гаусса
Метод Гаусса − это метод перехода от исходной системы линейных уравнений (при помощи эквивалентных преобразований) к системе, которая решается проще, чем исходная система.
Эквивалентными преобразованиями системы линейных уравнений являются:
- перемена местами двух уравнений в системе,
- умножение какого-либо уравнения в системе на ненулевое действительное число,
- прибавление к одному уравнению другого уравнения, умноженного на произвольное число.
Рассмотрим систему линейных уравнений:
(1) |
Запишем систему (1) в матричном виде:
Ax=b | (2) |
(3) |
A -называется матрица коэффициентов системы, b − правая часть ограничений, x − вектор переменных, которую нужно найти. Пусть rang(A )=p .
Эквивалентные преобразования не меняют ранг матрицы коэффициентов и ранг расширеннной матрицы системы. Не меняется также множество решений системы при эквивалентных преобразованиях. Суть метода Гаусса заключается в приведении матрцы коэффициентов A к диагональному или ступенчатому.
Построим расшренную матрицу системы:
На следующем этапе обнуляем все элементы столбца 2, ниже элемента . Если данный элемент нулевой, то эту строку меняем местами со строкой, лежащий ниже данной строки и имеющий ненулевой элемент во втором столбце. Далее обнуляем все элементы столбца 2 ниже ведущего элемента a 22 . Для этого сложим строки 3, ... m со строкой 2, умноженной на −a 32 /a 22 , ..., −a m2 /a 22 , соответственно. Продолжая процедуру, получим матрицу диагонального или ступенчатого вида. Пусть полученная расширенная матрица имеет вид:
(7) |
Так как rangA=rang (A|b ), то множество решений (7) есть (n−p )− многообразие. Следовательно n−p неизвестных можно выбрать произвольно. Остальные неизвестные из системы (7) вычисляются так. Из последнего уравнения выражаем x p через остальные переменные и вставляем в предыдущие выражения. Далее из предпоследнего уравнения выражаем x p−1 через остальные переменные и вставляем в предыдущие выражения и т.д. Рассмотрим метод Гаусса на конкретных примерах.
Примеры решения системы линейных уравнений методом Гаусса
Пример 1. Найти общее решение системы линейных уравнений методом Гаусса:
Обозначим через a ij элементы i -ой строки и j -ого столбца.
Исключим элементы 1-го столбца матрицы ниже элемента a 1 1 . Для этого сложим строки 2,3 со строкой 1, умноженной на -2/3,-1/2 соответственно:
Делим каждую строку матрицы на соответствующий ведущий элемент (если ведущий элемент существует):
Подставив верхние выражения в нижние, получим решение.
Пусть задана система линейных алгебраических уравнений, которую необходимо решить (найти такие значения неизвестных хi, что обращают каждое уравнение системы в равенство).
Мы знаем, что система линейных алгебраических уравнений может:
1) Не иметь решений (бытьнесовместной
).
2) Иметь бесконечно много решений.
3) Иметь единственное решение.
Как мы помним,правило Крамера и матричный методнепригодны в тех случаях, когда система имеет бесконечно много решений или несовместна. Метод Гаусса – наиболее мощный и универсальный инструмент для нахождения решения любой системы линейных уравнений , который в каждом случае приведет нас к ответу! Сам алгоритм метода во всех трёх случаях работает одинаково. Если в методах Крамера и матричном необходимы знания определителей, то для применения метода Гаусса необходимо знание только арифметических действий, что делает его доступным даже для школьников начальных классов.
Преобразования расширенной матрицы (это матрица системы - матрица, составленная только из коэффициентов при неизвестных, плюс столбец свободных членов) системы линейных алгебраических уравнений в методе Гаусса:
1) с троки матрицыможно переставлять местами.
2) если в матрице появились (или есть) пропорциональные (как частный случай – одинаковые) строки, то следуетудалить из матрицы все эти строки кроме одной.
3) если в матрице в ходе преобразований появилась нулевая строка, то ее также следует удалить .
4) строку матрицы можноумножить (разделить) на любое число,отличное от нуля.
5) к строке матрицы можноприбавить другую строку, умноженную на число , отличное от нуля.
В методе Гаусса элементарные преобразования не меняют решение системы уравнений.
Метод Гаусса состоит из двух этапов:
- «Прямой ход» - с помощью элементарных преобразований привести расширенную матрицу системы линейных алгебраических уравнений к «треугольному» ступенчатому виду: элементы расширенной матрицы, расположенные ниже главной диагонали, равны нулю (ход «сверху-вниз»). Например, к такому виду:
Для этого выполним следующие действия:
1) Пусть мы рассматриваем первое уравнение системы линейных алгебраических уравнений и коэффициент при х 1 равен К. Второе, третье и т.д. уравнения преобразуем следующим образом: каждое уравнение (коэффициенты при неизвестных, включая свободные члены) делим на коэффициент при неизвестном х 1 , стоящий в каждом уравнении, и умножаем на К. После этого из второго уравнения (коэффициенты при неизвестных и свободные члены) вычитаем первое. Получаем при х 1 во втором уравнении коэффициент 0. Из третьего преобразованного уравнения вычитаем первое уравнение, так до тех пор, пока все уравнения, кроме первого, при неизвестном х 1 не будут иметь коэффициент 0.
2) Переходим к следующему уравнению. Пусть это будет второе уравнение и коэффициент при х 2 равен М. Со всеми «нижестоящими» уравнениями поступаем так, как описано выше. Таким образом, «под» неизвестной х 2 во всех уравнениях будут нули.
3) Переходим к следующему уравнению и так до тех пора, пока не останется одна последняя неизвестная и преобразованный свободный член.
- «Обратный ход» метода Гаусса – получение решения системы линейных алгебраических уравнений (ход «снизу-вверх»). Из последнего «нижнего» уравнения получаем одно первое решение – неизвестную х n . Для этого решаем элементарное уравнение А*х n = В. В примере, приведенном выше, х 3 = 4. Подставляем найденное значение в «верхнее» следующее уравнение и решаем его относительно следующей неизвестной. Например, х 2 – 4 = 1, т.е. х 2 = 5. И так до тех пор, пока не найдем все неизвестные.
Пример.
Решим систему линейных уравнений методом Гаусса, как советуют некоторые авторы:
Запишем расширенную матрицу системы и с помощью элементарных преобразований приведем ее к ступенчатому виду:
Смотрим на левую верхнюю «ступеньку». Там у нас должна быть единица. Проблема состоит в том, что в первом столбце единиц нет вообще, поэтому перестановкой строк ничего не решить. В таких случаях единицу нужно организовать с помощью элементарного преобразования. Обычно это можно сделать несколькими способами. Поступим так:
1 шаг
. К первой строке прибавляем вторую строку, умноженную на –1. То есть, мысленно умножили вторую строку на –1 и выполнили сложение первой и второй строки, при этом вторая строка у нас не изменилась.
Теперь слева вверху «минус один», что нас вполне устроит. Кто хочет получить +1, может выполнить дополнительное действие: умножить первую строку на –1 (сменить у неё знак).
2 шаг . Ко второй строке прибавили первую строку, умноженную на 5. К третьей строке прибавили первую строку, умноженную на 3.
3 шаг . Первую строку умножили на –1, в принципе, это для красоты. У третьей строки также сменили знак и переставили её на второе место, таким образом, на второй «ступеньке у нас появилась нужная единица.
4 шаг . К третьей строке прибавили вторую строку, умноженную на 2.
5 шаг . Третью строку разделили на 3.
Признаком, который свидетельствует об ошибке в вычислениях (реже – об опечатке), является «плохая» нижняя строка. То есть, если бы у нас внизу получилось что-нибудь вроде (0 0 11 |23) , и, соответственно, 11x 3 = 23, x 3 = 23/11, то с большой долей вероятности можно утверждать, что допущена ошибка в ходе элементарных преобразований.
Выполняем обратный ход, в оформлении примеров часто не переписывают саму систему, а уравнения «берут прямо из приведенной матрицы». Обратный ход, напоминаю, работает «снизу вверх». В данном примере получился подарок:
x 3 = 1
x 2 = 3
x 1 + x 2 – x 3 = 1, следовательно x 1 + 3 – 1 = 1, x 1 = –1
Ответ :x 1 = –1, x 2 = 3, x 3 = 1.
Решим эту же систему по предложенному алгоритму. Получаем
4 2 –1 1
5 3 –2 2
3 2 –3 0
Разделим второе уравнение на 5, а третье – на 3. Получим:
4 2 –1 1
1 0.6 –0.4 0.4
1 0.66 –1 0
Умножим второе и третье уравнения на 4, получим:
4 2 –1 1
4 2,4 –1.6 1.6
4 2.64 –4 0
Вычтем из второго и третьего уравнений первое уравнение, имеем:
4 2 –1 1
0 0.4 –0.6 0.6
0 0.64 –3 –1
Разделим третье уравнение на 0,64:
4 2 –1 1
0 0.4 –0.6 0.6
0 1 –4.6875 –1.5625
Умножим третье уравнение на 0,4
4 2 –1 1
0 0.4 –0.6 0.6
0 0.4 –1.875 –0.625
Вычтем из третьего уравнения второе, получим «ступенчатую» расширенную матрицу:
4 2 –1 1
0 0.4 –0.6 0.6
0 0 –1.275 –1.225
Таким образом, так как в процессе вычислений накапливалась погрешность, получаем х 3 = 0,96 или приблизительно 1.
х 2 = 3 и х 1 = –1.
Решая таким образом, Вы никогда не запутаетесь в вычислениях и не смотря на погрешности вычислений, получите результат.
Такой способ решения системы линейных алгебраических уравнений легко программируем и не учитывает специфические особенности коэффициентов при неизвестных, ведь на практике (в экономических и технических расчетах) приходиться иметь дело именно с нецелыми коэффициентами.
Желаю успехов! До встречи на занятиях! Репетитор Дмитрий Айстраханов .
сайт, при полном или частичном копировании материала ссылка на первоисточник обязательна.