Системы линейных уравнений. определители

Пусть дана система линейных уравнений

Коэффициенты при неизвестных составляют прямоугольную лицу

называемую матрицей системы. Первый индекс у коэффициента aij означает номер уравнения, второй - номер неизвестного, при котором стоит этот коэффициент. Коэффициенты Ь, , Ь гп называются свободными членами уравнений системы. Если свободные члены равны нулю, то система называется однородной , в противном случае - неоднородной. Матрицу

называют расширенной матрицей системы (2.1).

Решение системы (2.1) - это любой упорядоченный набор (ад, Х 2 , ? ??,х п) из п чисел, при подстановке которых в уравнения системы вместо соответствующих неизвестных каждое уравнение системы превращается в тождество. Система, не имеющая ни одного решения, называется несовместной , или противоречивой. Система, имеющая хотя бы одно решение, называется совместной.

Совместные системы подразделяют на определенные, обладающие единственным решением, и неопределенные, обладающие большим числом решений. Однородная система всегда совместна, так как имеет по крайней мере нулевое решение х - Х 2 - ... = х п = 0.

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

Над уравнениями системы обычно приходится проводить следующие элементарные преобразования:

  • 1) умножение обеих частей какого-либо уравнения на число, отличное от нуля;
  • 2) прибавление (вычитание) к одному уравнению другого, умноженного на некоторое число;
  • 3) перестановку уравнений;
  • 4) вычеркивание уравнений вида 0 х + 0 Х 2 + + 0 х п = 0, т.е. тождеств 0 = 0;
  • 5) перестановку неизвестных в системе уравнений.

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

Предположим, что коэффициент ац системы (2.1) отличен от нуля. Этого всегда можно добиться, переставляя в случае необходимости уравнения системы или неизвестные в ней и меняя нумерацию неизвестных. Умножим первое уравнение на а 2 /ац и вычтем из второго уравнения, затем на а^/ац и вычтем из третьего уравнения и т.д. Наконец, умножим первое уравнение на a m ja и вычтем из последнего уравнения. В результате неизвестное х будет исключено из всех уравнений, кроме первого, и система примет вид:

В системе (2.2) следует вычеркнуть уравнения вида 0 х + 0 Х 2 + ...+ +0 х п = 0, если такие появились. На этом первый шаг метода Гаусса заканчивается. Элемент Дц называют ведущим элементом этого шага.

Следующие шаги прямого хода метода Гаусса осуществляются аналогично. Так, на втором шаге при а 22 ^ 0 последовательно умножаем второе уравнение на а" 32 /а 22 , а! А2 /а! 22 , ..., a" m2 fa 22 и вычитаем его из 3-го, 4-го, ..., m-го уравнений. В результате неизвестное Х 2 исключается из всех уравнений, кроме 1-го и 2-го. На третьем шаге неизвестное хз исключается из всех уравнений, кроме первых трех, и т.д.

Возможно, что на некотором шаге прямого хода метода Гаусса встретится уравнение вида

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

Для упрощения записи в системе (2.4) штрихи над коэффициентами опущены. В ней не более т уравнений, т.е. г ^ т, так как некоторые уравнения, возможно, были приведены к виду 0 = 0 и вычеркнуты, и, очевидно также, что г ^ п.

При г = п система (2.4) имеет треугольный вид:

и в ней легко совершить обратный ход метода Гаусса. Для этого из последнего уравнения этой системы найдем значение неизвестного х п. Подставив его в предпоследнее уравнение, найдем значение.x n _i. Продолжая так далее, однозначно определим все неизвестные х, Х 2 , ..., х п. Следовательно, если система (2.1) при прямом ходе метода Гаусса сводится к системе треугольного вида, то такая система определенная, т.е. имеет единственное решение.

При г система (2.4) имеет вид трапеции. В ней неизвестные х, Х 2 , ..., х г принимают за главные, а неизвестные х+, х г+ 2 , ..., х п - за свободные. Свободные неизвестные могут принимать любые фиксированные значения. Полагая x r+ = 7 r +i, х г+ 2 = Ъ-+ 2 , , х п = 7 п, где 7r+i, 7г+2? , 7п - произвольные постоянные, и проведя в системе обратный ход метода Гаусса, получим формулы:

которые составляют общее решение системы (2.1). Из общего решения (2.С) при конкретных значениях 7 r +i, 7г+2, , 7п будут получаться частные решения системы (2.1). Так как каждое свободное неизвестное может принимать бесчисленное множество значений, система (2.1) при г т.е. в случае, когда она приводится к трапецеидальному виду, обладает бесчисленным множеством решений. Это справедливо для совместных систем, имеющих меньше уравнений, чем неизвестных, и, в частности, для однородных, имеющих меньше уравнений, чем неизвестных.

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

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

Решение. Оставляя в расширенной матрице системы

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

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

Вычеркивая здесь третью строку, придем к матрице

которая соответствует системе

Отсюда, совершая обратный ход метода Гаусса, найдем

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

Решение. Если в расширенной матрице системы

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

Строка (0 0 0 | - 5) соответствует уравнению 0 х + 0 х 2 + 0 хз = -5. Наличие такого уравнения указывает на несовместность рассматриваемой системы. ?

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

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


Последняя матрица этой цепочки соответствует системе

Полагая здесь хз = 73 (77 - произвольная постоянная) и проводя обратный ход метода Гаусса, получим общее решение:

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

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

целесообразно сначала из первого уравнения системы вычесть удвоенное третье, а остальные оставить без изменения. Тогда получим систему

в которой метод Гаусса проводится уже легко. Дополнительные преобразования совершаются также над матрицами.

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

Метод Гаусса – это просто! Почему? Известный немецкий математик Иоганн Карл Фридрих Гаусс еще при жизни получил признание величайшего математика всех времен, гения и даже прозвище «короля математики». А всё гениальное – просто! Кстати, портрет Гаусса красовался на купюре в 10 дойчмарок (до введения евро), и до сих пор Гаусс загадочно улыбается немцам с обычных почтовых марок.

Метод Гаусса прост тем, что для его освоения ДОСТАТОЧНО ЗНАНИЙ ПЯТИКЛАССНИКА. Про миноры и алгебраические дополнения можно на время забыть! Необходимо уметь складывать и умножать! Не случайно метод последовательного исключения неизвестных преподаватели часто рассматривают на школьных математических факультативах.

Парадокс, но у студентов метод Гаусса вызывает наибольшие сложности. Ничего удивительного – всё дело в методике, и мы постараемся в доступной форме рассказать об алгоритме метода.

Сначала немного систематизируем знания о системах линейных уравнений. Система линейных уравнений может:

1) Иметь единственное решение.

2) Иметь бесконечно много решений.

3) Не иметь решений (быть несовместной ).

Метод Гаусса – наиболее мощный и универсальный инструмент для нахождения решения любой системы линейных уравнений. Как мы помним, правило Крамера и матричный метод непригодны в тех случаях, когда система имеет бесконечно много решений или несовместна. А метод последовательного исключения неизвестных в любом случае приведет нас к ответу! На данном уроке мы вновь рассмотрим метод Гаусса для случая №1 (единственное решение системы), под ситуации пунктов №№ 2-3 отведена статья Несовместные системы и системы с общим решением. Заметим, что сам алгоритм метода во всех трёх случаях работает одинаково.

Вернемся к простейшей системе

И решим ее методом Гаусса.

На первом этапе запишем так называемую расширенную матрицу системы :

По какому принципу записаны коэффициенты, думаем, всем видно.

Примечание: Расширенная матрица системы получается из исходной с помощью «операции наращивания строк / столбцов». В данном случае матрицу нарастили за счёт столбца свободных членов исходной системы уравнений.

Примечание: Кроме перечисленных ранее 6-и алгебраических операций с матрицами и «операции наращивания» существует ещё «операция отбрасывания строк/столбцов». С помощью «операции отбрасывания строк/столбцов» составляют, например, подматрицы, определители которых являются минорами элементов матрицы.

Вертикальная черта внутри матрицы не несёт никакого математического смысла – это просто линия отчёркивания для удобства оформления.

Определение: Матрица системы – это матрица, составленная только из коэффициентов при неизвестных переменных системы линейных уравнений.

Определение: Расширенная матрица системы – это матрица системы, которую нарастили справа на столбец свободных членов.

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

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

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

Существуют следующие элементарные преобразования:

1) Строки матрицы можно переставлять местами. Например, в рассматриваемой матрице можно безболезненно переставить первую и вторую строки:

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

Рассмотрим, например матрицу . В данной матрице последние три строки пропорциональны, поэтому достаточно оставить только одну из них:

.

3) Если в матрице в ходе преобразований появилась нулевая строка, то ее также следует удалить . Рисовать не будем, понятно, нулевая строка – это строка, в которой одни нули .

4) Строку матрицы можно умножить (разделить) на любое число, отличное от нуля . Рассмотрим, например, матрицу . Здесь целесообразно первую строку разделить на –3, а вторую строку – умножить на 2: . Данное действие очень полезно, поскольку упрощает дальнейшие преобразования матрицы.

5) Это преобразование вызывает наибольшие затруднения, но на самом деле ничего сложного тоже нет. К строке матрицы можно прибавить другую строку, умноженную на число , отличное от нуля.

Рассмотрим нашу матрицу из практического примера: . Сначала распишем преобразование очень подробно.

Умножаем первую строку на (-2): , далее ко второй строке прибавляем первую строку, оставляя первую без изменений : . Теперь первую строку можно разделить «обратно» на (–2): .

Как видите, строка, которую ПРИБАВЛЯЛИ не изменилась . Всегда меняется строка, К КОТОРОЙ ПРИБАВЛЯЮТ .

На практике так подробно, конечно, не расписывают, а пишут короче:

Еще раз: ко второй строке прибавили первую строку, умноженную на (–2) . Умножают строку обычно устно или на черновике, при этом мысленный ход расчётов примерно такой:

«Переписываю матрицу и переписываю первую строку: »

«Сначала первый столбец. Внизу мне нужно получить ноль. Поэтому единицу вверху умножаю на –2: , и ко второй строке прибавляю первую: 2 + (–2) = 0.

Записываю результат во вторую строку: »

«Теперь второй столбец. Вверху –1 умножаю на –2: (-1∙(-2) = 2 ). Ко второй строке прибавляю первую: 1 + 2 = 3. Записываю результат во вторую строку:

»

«И третий столбец. Вверху –5 умножаю на –2: (-5∙(-2) = 10 ). Ко второй строке прибавляю первую: (–7 + 10 = 3 ). Записываю результат во вторую строку:

»

Пожалуйста, тщательно осмыслите этот пример и разберитесь в последовательном алгоритме вычислений, если вы это поняли, то метод Гаусса практически «в кармане». Но, конечно, над этим преобразованием мы еще поработаем.

Повторим: «Элементарные преобразования не изменяют решение системы»

ВНИМАНИЕ!: рассмотренные манипуляции нельзя использовать , если Вам предложено задание, где матрицы даны «сами по себе». Например, при «классических» действиях с матрицами что-то переставлять внутри матриц ни в коем случае нельзя!

Вернемся к нашей системе . Она уже почти решена.

Что просит Гаусс? Он говорит: «Запишите расширенную матрицу системы и с помощью элементарных преобразований приведите ее к ступенчатому виду ».

В данном случае для этого

(1) Ко второй строке прибавьте первую строку, умноженную на –2. Кстати, почему первую строку умножаем именно на –2? Для того чтобы внизу получить ноль, а значит, избавиться от одной переменной во второй строке.

(2) Разделите вторую строку на 3. Почему? Чтобы вторая строка давала сразу значение второй переменной.

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

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

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

Теперь систему нужно «раскрутить» в обратном направлении – снизу вверх, этот процесс называется обратным ходом метода Гаусса .

В нижнем уравнении у нас уже готовый результат: . Рассмотрим первое уравнение системы и подставим в него уже известное значение «игрек»:

Ответ:

Рассмотрим наиболее распространенную ситуацию, когда методом Гаусса требуется решить систему трёх линейных уравнений с тремя неизвестными.

Пример 1

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

Запишем расширенную матрицу системы:

Сейчас мы сразу нарисуем результат, к которому мы придём в ходе решения:

.

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

Сначала смотрим на левое верхнее число:

.

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

Теперь первая строка у нас останется неизменной до конца решения . Уже легче.

Единица в левом верхнем углу организована. Теперь нужно получить нули вот на этих местах:

Нули получаем как раз с помощью «трудного» преобразования. Сначала разбираемся со второй строкой (2, –1, 3, 13). Что нужно сделать, чтобы на первой позиции получить ноль? Нужно ко второй строке прибавить первую строку, умноженную на –2 . Мысленно или на черновике умножаем первую строку на –2: (–2, –4, 2, –18).

И последовательно проводим (опять же мысленно или на черновике) сложение, т. е. ко второй строке прибавляем первую строку, уже умноженную на –2 :

Результат записываем во вторую строку:

Аналогично разбираемся с третьей строкой (3, 2, –5, –1). Чтобы получить на первой позиции ноль, нужно к третьей строке прибавить первую строку, умноженную на –3 .

Мысленно или на черновике умножаем первую строку на –3: (–3, –6, 3, –27). И к третьей строке прибавляем первую строку, умноженную на –3 :

Результат записываем в третью строку:

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

Не нужно считать всё сразу и одновременно . Порядок вычислений и «вписывания» результатов последователен и обычно такой: сначала переписываем первую строку, и пыхтим себе потихонечку – ПОСЛЕДОВАТЕЛЬНО иВНИМАТЕЛЬНО :

.

А мысленный ход самих расчётов мы уже рассмотрели выше.

В данном примере это сделать легко, вторую строку делим на –5 (поскольку там все числа делятся на 5 без остатка). Заодно делим третью строку на –2, ведь чем меньше числа, тем проще решение:

На заключительном этапе элементарных преобразований нужно получить еще один ноль здесь:

Для этого к третьей строке прибавляем вторую строку, умноженную на –2 :

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

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

Теперь в действие вступает «обратный ход» метода Гаусса. Уравнения «раскручиваются» снизу вверх.

В третьем уравнении у нас уже готовый результат:

Смотрим на второе уравнение: . Значение «зет» уже известно, таким образом:

Пример 3

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

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

Поступим так:

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

Теперь слева вверху (–1), что нас вполне устроит. Кто хочет получить (+1), может выполнить дополнительное телодвижение: умножить первую строку на (–1), сменив у неё знак. Дальше алгоритм работает уже по накатанной колее:

.

(2) Ко второй строке прибавили первую строку, умноженную на 5. К третьей строке прибавили первую строку, умноженную на 3.

(3) Первую строку умножили на (–1). В принципе, это для красоты. У третьей строки также сменили знак и переставили её на второе место, таким образом, на второй «ступеньке у нас появилась нужная единица.

(4) К третьей строке прибавили вторую строку, умноженную на 2.

(5) Третью строку разделили на 3.

Заряжаем обратный ход, в оформлении примеров часто не переписывают саму систему, а уравнения «берут прямо из приведенной матрицы». Обратный ход, напоминаю, работает, снизу вверх. Да тут подарок получился:

Ответ: .

Пример 4

Решить систему линейных уравнений методом Гаусса

Это пример для самостоятельного решения, он несколько сложнее. Ничего страшного, если кто-нибудь запутается. Полное решение и образец оформления в конце урока. Ваш ход решения может отличаться от нашего хода решения.

Данный онлайн калькулятор находит решение системы линейных уравнений (СЛУ) методом Гаусса. Дается подробное решение. Для вычисления выбирайте количество переменных и количество уравнений. Затем введите данные в ячейки и нажимайте на кнопку "Вычислить."

x 1

+x 2

+x 3

x 1

+x 2

+x 3

x 1

+x 2

+x 3

=

=

=

Представление чисел:

Целые числа и (или) Обыкновенные дроби
Целые числа и (или) Десятичные дроби

Число знаков после десятичного разделителя

×

Предупреждение

Очистить все ячейки?

Закрыть Очистить

Инструкция ввода данных. Числа вводятся в виде целых чисел (примеры: 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 соответственно:

Делим каждую строку матрицы на соответствующий ведущий элемент (если ведущий элемент существует):

Подставив верхние выражения в нижние, получим решение.

Сегодня разбираемся с методом Гаусса для решения систем линейных алгебраических уравнений. О том, что это за системы, можно почитать в предыдущей статье, посвященной решению тех же СЛАУ методом Крамера. Метод Гаусса не требует каких-то специфических знаний, нужна лишь внимательность и последовательность. Несмотря на то что с точки зрения математики для его применения хватит и школьной подготовки, у студентов освоение этого метода часто вызывает сложности. В этой статье попробуем свести их на нет!

Метод Гаусса

Метод Гаусса – наиболее универсальный метод решения СЛАУ (за исключением ну уж очень больших систем). В отличие от рассмотренного ранее , он подходит не только для систем, имеющих единственное решение, но и для систем, у которых решений бесконечное множество. Здесь возможны три варианта.

  1. Система имеет единственное решение (определитель главной матрицы системы не равен нулю);
  2. Система имеет бесконечное множество решений;
  3. Решений нет, система несовместна.

Итак, у нас есть система (пусть у нее будет одно решение), и мы собираемся решать ее методом Гаусса. Как это работает?

Метод Гаусса состоит из двух этапов – прямого и обратного.

Прямой ход метода Гаусса

Сначала запишем расширенную матрицу системы. Для этого в главную матрицу добавляем столбец свободных членов.

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

Что можно делать:

  1. Можно переставлять строки матрицы местами;
  2. Если в матрице есть одинаковые (или пропорциональные) строки, можно удалить их все, кроме одной;
  3. Можно умножать или делить строку на любое число (кроме нуля);
  4. Нулевые строки удаляются;
  5. Можно прибавлять к строке строку, умноженную на число, отличное от нуля.

Обратный ход метода Гаусса

После того как мы преобразуем систему таким образом, одна неизвестная Xn становится известна, и можно в обратном порядке найти все оставшиеся неизвестные, подставляя уже известные иксы в уравнения системы, вплоть до первого.

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

Пример решения системы уравнений методом Гаусс

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

Сначала запишем расширенную матрицу:

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

Затем умножим 3-ую строку на (-1). Добавим 3-ую строку к 2-ой:

Умножим 1-ую строку на (6). Умножим 2-ую строку на (13). Добавим 2-ую строку к 1-ой:

Вуаля - система приведена к соответствующему виду. Осталось найти неизвестные:

Система в данном примере имеет единственное решение. Решение систем с бесконечным множеством решений мы рассмотрим в отдельной статье. Возможно, сначала Вы не будете знать, с чего начать преобразования матрицы, но после соответствующей практики набъете руку и будете щелкать СЛАУ методом Гаусса как орешки. А если Вы вдруг столкнетесь со СЛАУ, которая окажется слишком крепким орешком, обращайтесь к нашим авторам! вы можете, оставив заявку в Заочнике. Вместе мы решим любую задачу!

Раздел 3. Численные методы решения уравнений

Виды математических моделей (уравнений) в теории электрических цепей

1. - системы линейных алгебраических уравнений –

линейные цепи постоянного и синусоидального переменного (комплексный метод) тока.

2 . - системы нелинейных алгебраических или

трансцендентных уравнений – нелинейные цепи постоянного или синусоидального тока.

3. . системы нелинейных дифференциальных

уравнений первого порядка в обыкновенных производных – переходные процессы в нелинейных цепях.

Здесь F и ψ – вектор-функции, т.е. эквивалентно записи:

f 1 (X,b 1) = 0

f 2 (X,b 2) = 0

…………

f n (X,b n) = 0

а - записи:

ψ 1 (dX/dt,X,b 1 ,t) = 0

ψ 2 (dX/dt,X,b 2 ,t) = 0

…………………..

ψ n (dX/dt,X,b n ,t) = 0

Рассмотрим наиболее эффективные методы решения этих уравнений.

Численные методы решения систем линейных алгебраических уравнений (ЛАУ)

Метод Гаусса (исключения неизвестных)

Методы решения ЛАУ имеют важное значение, так как они применяются (итерационно) для решения более сложных уравнений.

Пусть система ЛАУ задана в виде:

,

где - квадратная матрица n – го порядка с ненулевыми диагональными элементами ; - вектор неизвестных; - вектор правых частей.

Алгоритм метода Гаусса состоит из прямого и обратного хода. Во время прямого хода осуществляется последовательное исключение неизвестных. Система приобретает вид:

Пересчет коэффициентов производится по формуле:

, где i, j = k+1, …n при исключение k -го неизвестного.

При этом столбец правых частей удобно рассматривать как n + 1 столбец матрицы коэффициентов , т.е. j = k+1, …n+1.

Обратный ход заключается в определении неизвестных, начиная с последнего уравнения где осталась одна неизвестная x n . Полученное значение x n подставляется в предыдущее уравнение и определяется x n -1 и т.д.

Для произвольного x k получается следующая формула:

где k = n, n -1,…1.

Трудоемкость метода Гаусса оценивается количеством выполняемых арифметических операций:

.

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



Основная идея метода разреженных матриц состоит в учете при вычислениях и хранении только ненулевых элементов матрицы . Степень разреженности матрицы характеризуется коэффициентом заполнения:

где n ннэ –число ненулевых элементов.

Существуют матрицы коэффициентов специального вида: ленточные, когда ненулевые элементы располагаются вдоль главной диагонали; и блочно-диагональные, когда вдоль главной диагонали располагаются ненулевые блоки. Еще встречаются блочно-диагональные с окаймлением.

Пример ленточной матрицы Пример блочно-диагональной матрицы


Пример блочно-диагональной матрицы с окаймлением

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

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