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




Лабораторная работа 4. РЕШЕНИЕ УРАВНЕНИЙ

Общие сведения

Огромное количество задач вычислительной математики связано с

решением нелинейных алгебраических уравнений, а также систем таких

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

возникает зачастую на промежуточных шагах, при реализации фрагментов

более сложных алгоритмов (к примеру, при расчетах дифференциальных

уравнений при помощи разностных схем и т. п.).

Численное решение нелинейного уравнения

Алгоритм приближенного решения уравнения f(x)=0 состоит из двух

этапов:

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

2. получения приближенного решения с заданной точностью с помощью функции root.

Если после многих итераций MathCAD не находит подходящего приближения, то появится сообщение

(отсутствует сходимость).

Пример 1. Решение уравнения с помощью функции ROOT

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

• уравнение не имеет корней;

• корни уравнения расположены далеко от начального приближения;

• выражение имеет комплексный корень, но начальное приближение было вещественным.

Чтобы установить причину ошибки, исследуйте график f(x). Он поможет выяснить наличие корней уравнения f(x) = 0 и, если они есть, то определить приблизительно их значения. Чем точнее выбрано начальное приближение корня, тем быстрее будет root сходиться.

Нахождение корней полинома/

Для нахождения корней выражения, имеющего вид:

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

Функция Polyroots(v) - возвращает корни полинома степени n. Коэффициенты полинома находятся в векторе v длины n + 1. Возвращает вектор длины n, состоящий из корней полинома.

Пример 2. Нахождение корней полинома

Решение систем уравнений.

Решение систем уравнений матричным методом.

Рассмотрим систему n линейных алгебраических уравнений относительно n неизвестных х1, х2, …, хn:

 

Система линейных уравнений может быть записана в матричном виде:

Ах = b,

где:

Если det A ≠ 0 то система или эквивалентное ей матричное уравнение имеет единственное решение.

Решение систем уравнений с помощью функции Lsolve Системы линейных уравнений удобно решать с помощью функции lsolve.

Функция lsolve(А, b) - возвращает вектор решения x такой, что Ах = b.

Пример 3. Решение системы уравнений

Запишем в матричном виде:

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

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

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

единичная матрица. Последний, (n + 1) столбец этой матрицы содержит

решение системы.

В MathCAD прямой и обратный ходы метода Гаусса выполняет функция

rref(A).

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

Формирование расширенной матрицы системы:

Приведение расширенной матрицы к ступенчатому виду

(прямой и обратный ходы метода Гаусса

 

Решение систем уравнений с помощью функций Find или Minner

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

выполнить следующее:

1. Задать начальное приближение для всех неизвестных, входящих в

систему уравнений. MathCAD решает систему с помощью итерационных

методов;

2. Напечатать ключевое слово Given. Оно указывает MathCAD, что далее

следует система уравнений;

 

3. Введите уравнения и неравенства в любом порядке. Используйте [Ctrl]=

для печати символа =. Между левыми и правыми частями неравенств может

стоять любой из символов <, >, ≥ и ≤;

4. Введите любое выражение, которое включает функцию Find,

например: х:= Find(х, у).

Ключевое слово Given, уравнения и неравенства, которые следуют за ним, и

какое - либо выражение, содержащее функцию Find, называют блоком

решения уравнений.

Пример 5. Решение системы уравнений с помощью функции Fin d

 

 



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2017-04-03 Нарушение авторских прав и Нарушение персональных данных


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