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

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

Содержание

1 Постановка задачи
2 Численные методы решения уравнений
- 2.1 Метод простой итерации
  - 2.1.1 Применительно к СЛАУ
  - 2.1.2 Алгоритм
- 2.2 Метод Ньютона (метод касательных)
  - 2.2.1 Одномерный случай
  - 2.2.2 Многомерный случай
3 Литература
4 См. также

Постановка задачи

Рассмотрим методы численного решения уравнений и систем уравнений:

или

$\left\{ \begin{array}{lcr} f_1(x_1, x_2, \ldots, x_n) & =& 0 \\ \ldots & & \\ f_n(x_1, x_2, \ldots, x_n) & =& 0 \end{array}\right.$

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

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

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

Метод простой итерации

В основе метода заложено понятие сжимающего отображения. Определим терминологию:

Говорят, что функция осуществляет сжимающее отображение на , если

Тогда основная теорема будет выглядеть так:

Теорема Банаха (принцип сжимающих отображений).
Если — сжимающее отображение на , то:

— корень;
итерационная последовательность сходится к этому корню;
для очередного члена справедливо .

Поясним смысл параметра . Согласно теореме Лагранжа имеем:

Отсюда следует, что . Таким образом, для сходимости метода достаточно, чтобы

.........

и так далее, пока

Применительно к СЛАУ

Рассмотрим систему:

$\left\{ \begin{array}{ccc} a_{11} x_1 + \ldots + a_{1n} x_n & = & b_1 \\ \ldots & & \\ a_{n1} x_1 + \ldots + a_{nn} x_n & = & b_n \end{array}\right.$

Для неё итерационное вычисление будет выглядеть так:

$\left( \begin{array}{c} x_1\\ x_2\\ \vdots\\ x_n \end{array}\right)^{i+1} = \left( \begin{array}{cccc} a_{11}+1 & a_{12} & \ldots & a_{1n} \\ a_{21} & a_{22}+1 & \ldots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{n1} & a_{n2} & \ldots & a_{nn}+1 \end{array}\right) \left(\begin{array}{c} x_1\\ x_2\\ \vdots\\ x_n \end{array}\right)^{i}-\left(\begin{array}{c} b_1\\ b_2\\ \vdots\\ b_n \end{array}\right)$

Сходимость метода будет осуществлять

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

Решение уравнения cos(x)=x по методу простой итерации, очередная итерация: x_n+1=cos x_n, начальное приближение: x₁ = -1

Алгоритм

Условие преобразуется к виду , где — сжимающая
Задаётся начальное приближение и точность
Вычисляется очередная итерация
- Если , то и возврат к шагу 3.
- Иначе и остановка.

Решение уравнения f(x)=0 по методу Ньютона, начальное приближение: x₁=a.

Метод Ньютона (метод касательных)

Основная статья: Метод Ньютона

Одномерный случай

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

Воспользуемся тем, что , и получим окончательную формулу для :

С учётом этого сжимающая функция примет вид:

Тогда алгоритм нахождения численного решения уравнения сводится к итерационной процедуре вычисления:

Многомерный случай

Обобщим полученный результат на многомерный случай.

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

где .

Литература

Амосов А. А., Дубинский Ю. А., Копченова Н. П. Вычислительные методы для инженеров. — М.: Мир, 1998.
Бахвалов Н. С., Жидков Н. П., Кобельков Г. Г. Численные методы. — 8-е изд. — М.: Лаборатория Базовых Знаний, 2000.
Волков Е. А. Численные методы. — М.: Физматлит, 2003.
Коршунов Ю. М., Коршунов Ю. М. Математические основы кибернетики. — М.: Энергоатомиздат, 1972.
Калиткин Н. Н. Численные методы. — М.: Наука, 1978.

Ekb-oskab.ru

Прием лома металлов

Статьи