Метод Гаусса для решения систем линейных уравнений
Система линейных алгебраических уравнений
Понятие системы линейных алгебраических уравнений
Система уравнений – это условие, состоящее в одновременном выполнении нескольких уравнений относительно нескольких переменных. Системой линейных алгебраических уравнений (далее – СЛАУ), содержащей m уравнений и n неизвестных, называется система вида:
где числа aij называются коэффициентами системы, числа bi – свободными членами, aij и bi (i=1,…, m; b=1,…, n) представляют собой некоторые известные числа, а x1,…, xn – неизвестные. В обозначении коэффициентов aij первый индекс i обозначает номер уравнения, а второй j – номер неизвестного, при котором стоит этот коэффициент. Подлежат нахождению числа xn. Такую систему удобно записывать в компактной матричной форме: AX=B. Здесь А – матрица коэффициентов системы, называемая основной матрицей;
– вектор-столбец из неизвестных xj.
– вектор-столбец из свободных членов bi.
Произведение матриц А*Х определено, так как в матрице А столбцов столько же, сколько строк в матрице Х (n штук).
Расширенной матрицей системы называется матрица A системы, дополненная столбцом свободных членов
Решение системы линейных алгебраических уравнений
Решением системы уравнений называется упорядоченный набор чисел (значений переменных), при подстановке которых вместо переменных каждое из уравнений системы обращается в верное равенство.
Решением системы называется n значений неизвестных х1=c1, x2=c2,…, xn=cn, при подстановке которых все уравнения системы обращаются в верные равенства. Всякое решение системы можно записать в виде матрицы-столбца
Система уравнений называется совместной, если она имеет хотя бы одно решение, и несовместной, если она не имеет ни одного решения.
Совместная система называется определенной, если она имеет единственное решение, и неопределенной, если она имеет более одного решения. В последнем случае каждое ее решение называется частным решением системы. Совокупность всех частных решений называется общим решением.
Решить систему – это значит выяснить, совместна она или несовместна. Если система совместна, найти ее общее решение.
Две системы называются эквивалентными (равносильными), если они имеют одно и то же общее решение. Другими словами, системы эквивалентны, если каждое решение одной из них является решением другой, и наоборот.
Преобразование, применение которого превращает систему в новую систему, эквивалентную исходной, называется эквивалентным или равносильным преобразованием. Примерами эквивалентных преобразований могут служить следующие преобразования: перестановка местами двух уравнений системы, перестановка местами двух неизвестных вместе с коэффициентами у всех уравнений, умножение обеих частей какого-либо уравнения системы на отличное от нуля число.
Система линейных уравнений называется однородной, если все свободные члены равны нулю:
Однородная система всегда совместна, так как x1=x2=x3=…=xn=0 является решением системы. Это решение называется нулевым или тривиальным.
Метод исключения Гаусса
2.1 Сущность метода исключения Гаусса
Классическим методом решения систем линейных алгебраических уравнений является метод последовательного исключения неизвестных – метод Гаусса (его еще называют методом гауссовых исключений). Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе ступенчатого (или треугольного) вида, из которого последовательно, начиная с последних (по номеру) переменных, находятся все остальные переменные.
Процесс решения по методу Гаусса состоит из двух этапов: прямой и обратный ходы.
Прямой ход.
На первом этапе осуществляется так называемый прямой ход, когда путём элементарных преобразований над строками систему приводят к ступенчатой или треугольной форме, либо устанавливают, что система несовместна. А именно, среди элементов первого столбца матрицы выбирают ненулевой, перемещают его на крайнее верхнее положение перестановкой строк и вычитают получившуюся после перестановки первую строку из остальных строк, домножив её на величину, равную отношению первого элемента каждой из этих строк к первому элементу первой строки, обнуляя тем самым столбец под ним.
После того, как указанные преобразования были совершены, первую строку и первый столбец мысленно вычёркивают и продолжают пока не останется матрица нулевого размера. Если на какой-то из итераций среди элементов первого столбца не нашёлся ненулевой, то переходят к следующему столбцу и проделывают аналогичную операцию.
На первом этапе (прямой ход) система приводится к ступенчатому (в частности, треугольному) виду.
Приведенная ниже система имеет ступенчатый вид:
,
где
Коэффициенты aii называются главными (ведущими) элементами системы.
1-й шаг.
Будем считать, что элемент (если a11=0, переставим строки матрицы так, чтобы a11 не был равен 0. Это всегда возможно, т. к. в противном случае матрица содержит нулевой столбец, ее определитель равен нулю и система несовместна).
Преобразуем систему, исключив неизвестное х1 во всех уравнениях, кроме первого (используя элементарные преобразования системы). Для этого умножим обе части первого уравнения на и сложим почленно со вторым уравнением системы (или из второго уравнения почленно вычтем первое, умноженное на ). Затем умножим обе части первого уравнения на и сложим с третьим уравнением системы (или из третьего почленно вычтем первое, помноженное на ). Таким образом, последовательно умножаем первую строку на число и прибавляем к i-й строке, для i=2, 3, …, n.
Продолжая этот процесс, получим эквивалентную систему:
Здесь – новые значения коэффициентов при неизвестных и свободные члены в последних m-1 уравнениях системы, которые определяются формулами:
Таким образом, на первом шаге уничтожаются все коэффициенты, лежащие под первым ведущим элементом a110, на втором шаге уничтожаются элементы, лежащие под вторым ведущим элементом а22(1) (если a22(1)0) и т.д. Продолжая этот процесс и дальше, мы, наконец, на (m-1) шаге приведем исходную систему к треугольной системе.
Если в процессе приведения системы к ступенчатому виду появятся нулевые уравнения, т.е. равенства вида 0=0, их отбрасывают. Если же появится уравнение вида то это свидетельствует о несовместности системы.
На этом прямой ход метода Гаусса заканчивается.
Обратный ход.
На втором этапе осуществляется так называемый обратный ход, суть которого заключается в том, чтобы выразить все получившиеся базисные переменные через небазисные и построить фундаментальную систему решений, либо, если все переменные являются базисными, то выразить в численном виде единственное решение системы линейных уравнений.
Эта процедура начинается с последнего уравнения, из которого выражают соответствующую базисную переменную (она в нем всего одна) и подставляют в предыдущие уравнения, и так далее, поднимаясь по «ступенькам» наверх.
Каждой строчке соответствует ровно одна базисная переменная, поэтому на каждом шаге, кроме последнего (самого верхнего), ситуация в точности повторяет случай последней строки.
Примечание: на практике удобнее работать не с системой, а с расширенной ее матрицей, выполняя все элементарные преобразования над ее строками. Удобно, чтобы коэффициент a11 был равен 1 (уравнения переставить местами, либо разделить обе части уравнения на a11).
Примеры решения СЛАУ методом Гаусса
В данном разделе на трех различных примерах покажем, как методом Гаусса можно решить СЛАУ.
Пример 1. Решить СЛАУ 3-го порядка.
Обнулим коэффициенты при во второй и третьей строчках. Для этого домножим их на 2/3 и 1 соответственно и сложим с первой строкой:
Теперь обнулим коэффициент при в третьей строке, домножив вторую строку на 6 и вычитая из неё третью:
В результате мы привели исходную систему к треугольному виду, тем самым закончив первый этап алгоритма.
На втором этапе разрешим полученные уравнения в обратном порядке. Имеем:
из третьего;
из второго, подставив полученное ;
из первого, подставив полученные и .
В случае, если число уравнений в совместной системе получилось меньше числа неизвестных, то тогда ответ будет записываться в виде фундаментальной системы решений.
Пример 2. Решить неопределенную СЛАУ 4-го порядка:
В результате элементарных преобразований над расширенной матрицей системы
исходная система свелась к ступенчатой, где количество уравнений меньше, чем количество неизвестных:
Поэтому общее решение системы: x2=5×4–13×3–3; x1=5×4–8×3–1. Если положить, например, что x3=0, x4=0, то найдем одно из частных решений этой системы x1=-1, x2=-3, x3=0, x4=0.
Пример 3. Решить СЛАУ 4-ого порядка.
Условие:
baza-referat.ru
Метод Гаусса для решения систем линейных уравнений
Система линейных алгебраических уравнений
Понятие системы линейных алгебраических уравнений
Система уравнений – это условие, состоящее в одновременном выполнении нескольких уравнений относительно нескольких переменных. Системой линейных алгебраических уравнений (далее – СЛАУ), содержащей m уравнений и n неизвестных, называется система вида:
где числа aij называются коэффициентами системы, числа bi – свободными членами, aij и bi (i=1,…, m; b=1,…, n) представляют собой некоторые известные числа, а x1,…, xn – неизвестные. В обозначении коэффициентов aij первый индекс i обозначает номер уравнения, а второй j – номер неизвестного, при котором стоит этот коэффициент. Подлежат нахождению числа xn. Такую систему удобно записывать в компактной матричной форме: AX=B. Здесь А – матрица коэффициентов системы, называемая основной матрицей;
– вектор-столбец из неизвестных xj.
– вектор-столбец из свободных членов bi.
Произведение матриц А*Х определено, так как в матрице А столбцов столько же, сколько строк в матрице Х (n штук).
Расширенной матрицей системы называется матрица A системы, дополненная столбцом свободных членов
Решение системы линейных алгебраических уравнений
Решением системы уравнений называется упорядоченный набор чисел (значений переменных), при подстановке которых вместо переменных каждое из уравнений системы обращается в верное равенство.
Решением системы называется n значений неизвестных х1=c1, x2=c2,…, xn=cn, при подстановке которых все уравнения системы обращаются в верные равенства. Всякое решение системы можно записать в виде матрицы-столбца
Система уравнений называется совместной, если она имеет хотя бы одно решение, и несовместной, если она не имеет ни одного решения.
Совместная система называется определенной, если она имеет единственное решение, и неопределенной, если она имеет более одного решения. В последнем случае каждое ее решение называется частным решением системы. Совокупность всех частных решений называется общим решением.
Решить систему – это значит выяснить, совместна она или несовместна. Если система совместна, найти ее общее решение.
Две системы называются эквивалентными (равносильными), если они имеют одно и то же общее решение. Другими словами, системы эквивалентны, если каждое решение одной из них является решением другой, и наоборот.
Преобразование, применение которого превращает систему в новую систему, эквивалентную исходной, называется эквивалентным или равносильным преобразованием. Примерами эквивалентных преобразований могут служить следующие преобразования: перестановка местами двух уравнений системы, перестановка местами двух неизвестных вместе с коэффициентами у всех уравнений, умножение обеих частей какого-либо уравнения системы на отличное от нуля число.
Система линейных уравнений называется однородной, если все свободные члены равны нулю:
Однородная система всегда совместна, так как x1=x2=x3=…=xn=0 является решением системы. Это решение называется нулевым или тривиальным.
Метод исключения Гаусса
2.1 Сущность метода исключения Гаусса
Классическим методом решения систем линейных алгебраических уравнений является метод последовательного исключения неизвестных – метод Гаусса (его еще называют методом гауссовых исключений). Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе ступенчатого (или треугольного) вида, из которого последовательно, начиная с последних (по номеру) переменных, находятся все остальные переменные.
Процесс решения по методу Гаусса состоит из двух этапов: прямой и обратный ходы.
Прямой ход.
На первом этапе осуществляется так называемый прямой ход, когда путём элементарных преобразований над строками систему приводят к ступенчатой или треугольной форме, либо устанавливают, что система несовместна. А именно, среди элементов первого столбца матрицы выбирают ненулевой, перемещают его на крайнее верхнее положение перестановкой строк и вычитают получившуюся после перестановки первую строку из остальных строк, домножив её на величину, равную отношению первого элемента каждой из этих строк к первому элементу первой строки, обнуляя тем самым столбец под ним.
После того, как указанные преобразования были совершены, первую строку и первый столбец мысленно вычёркивают и продолжают пока не останется матрица нулевого размера. Если на какой-то из итераций среди элементов первого столбца не нашёлся ненулевой, то переходят к следующему столбцу и проделывают аналогичную операцию.
На первом этапе (прямой ход) система приводится к ступенчатому (в частности, треугольному) виду.
Приведенная ниже система имеет ступенчатый вид:
,
где
Коэффициенты aii называются главными (ведущими) элементами системы.
1-й шаг.
Будем считать, что элемент (если a11=0, переставим строки матрицы так, чтобы a11 не был равен 0. Это всегда возможно, т. к. в противном случае матрица содержит нулевой столбец, ее определитель равен нулю и система несовместна).
coolreferat.com
Метод Гаусса для решения систем линейных алгебраических уравнений (СЛАУ)
Теоретическая часть
Метод Гаусса представляет собой обобщение способа подстановки и состоит в последовательном исключении неизвестных до тех пор, пока не останется одно уравнение с одним неизвестным.
При этом матрица СЛАУ приводится треугольному виду, где ниже главной диагонали располагаются только нули.
Приведение матрицы к треугольному виду называется прямым ходом метода Гаусса. Обратный ход начинается с решения последнего уравнения и заканчивается определением первого неизвестного.
Имеем Ax=b, где A=[aij] — матрица размерности n?n, det A>0, b=(a1, n+1, …, an, n+1)T.
В предположении, что a11>0, первое уравнение системы (1):
делим на коэффициент a11, в результате получаем уравнение
Затем из каждого из остальных уравнений вычитается первое уравнение, умноженное на соответствующий коэффициент ai1. В результате эти уравнения преобразуются к виду
Первое неизвестное оказалось исключенным из всех уравнений, кроме первого. Далее предполагаем, что и исключаем неизвестное x2 из всех уравнений, начиная со второго, и т.д. В результате последовательного исключения неизвестных система уравнений преобразуется в систему уравнений с треугольной матрицей (2):
Совокупность проведенных действий называется прямым ходом метода Гаусса.
Из n—го уравнения системы (2) определяем
Реализация прямого хода требует арифметических операций, а обратного — арифметических операций.
Алгоритм
Описание входной информации: Str (Stlb) — количество строк (столбцов) в расширенной матрице, Ab [i, j] — матрица A (i — строки, j — столбцы).
Описание выходной информации: Ans — матрица-ответ.
Ссылки для скачивания:
Метод Гаусса на MatLab7
Метод Гаусса на Pascal ABC
Вернуться назад
pcfu.ru