Метод простых итераций
Материал из MachineLearning.
Строка 2: | Строка 2: | ||
Пусть есть функция <tex>y = f(x)</tex>.<br> | Пусть есть функция <tex>y = f(x)</tex>.<br> | ||
Требуется найти корень этой функции, то есть <tex>x</tex> при котором <tex>f(x)=0</tex><br> | Требуется найти корень этой функции, то есть <tex>x</tex> при котором <tex>f(x)=0</tex><br> | ||
- | Решение необходимо найти численно, то есть для реализации на ЭВМ. Для решения этой задачи предлагается использовать метод простых итераций | + | Решение необходимо найти численно, то есть для реализации на ЭВМ. Для решения этой задачи предлагается использовать метод простых итераций. |
== Метод простых итераций в общем виде == | == Метод простых итераций в общем виде == | ||
- | Заменеим исходное уравнение <tex>f(x)=0</tex> на эквивалентное <tex>g(x)=x</tex> | + | Заменеим исходное уравнение <tex>f(x)=0</tex> на эквивалентное <tex>g(x)=x</tex>,и будем строить итерации по правилу <tex>x_{n+1} = g(x_n)</tex>. Таким образом метод простой итерации - это одношаговый итерационный процесс. Для того, что бы начать данный процесс, необходимо знать начальное приближение <tex>x_1</tex>. Выясним условия сходимости метода. |
- | + | ||
- | + | ||
- | + | ||
- | Таким образом метод простой итерации - это одношаговый | + | |
===Сходимость метода простых итераций=== | ===Сходимость метода простых итераций=== | ||
Метод сходится, если при <tex>k \to \infty </tex> последовательность {<tex>x_n</tex>} имеет предел.<br> | Метод сходится, если при <tex>k \to \infty </tex> последовательность {<tex>x_n</tex>} имеет предел.<br> | ||
Строка 15: | Строка 11: | ||
===Метод релаксации=== | ===Метод релаксации=== | ||
+ | Так как для сходимости метода очень важен выбор функции <tex>g(x)</tex>, ее обычно берут вида <tex>g(x)=x+s(x)f(x)</tex>. Где <tex>s(x)</tex> не меняет знака на отрезке, на котором ищется корень функции.<br> | ||
Положим <tex>s(x) = c = const </tex> b и рассмотрим метод в этом случае.<br> | Положим <tex>s(x) = c = const </tex> b и рассмотрим метод в этом случае.<br> | ||
Тогда получим <tex>f(x_n) = \frac{x_{n+1}-x_{n}}{c}</tex>. <br> | Тогда получим <tex>f(x_n) = \frac{x_{n+1}-x_{n}}{c}</tex>. <br> |
Версия 08:33, 24 ноября 2008
Содержание |
Постановка задачи
Пусть есть функция .
Требуется найти корень этой функции, то есть при котором
Решение необходимо найти численно, то есть для реализации на ЭВМ. Для решения этой задачи предлагается использовать метод простых итераций.
Метод простых итераций в общем виде
Заменеим исходное уравнение на эквивалентное ,и будем строить итерации по правилу . Таким образом метод простой итерации - это одношаговый итерационный процесс. Для того, что бы начать данный процесс, необходимо знать начальное приближение . Выясним условия сходимости метода.
Сходимость метода простых итераций
Метод сходится, если при последовательность {} имеет предел.
Обозначим окресность точки радиуса .
Теорема. Если Липшиц непрерывна с константой на , то есть выполняется . При этом если также выполнено , то уравнение имеет решение на и метод простой итерации сходится к решению при любом выборе начального приближжения .
Метод релаксации
Так как для сходимости метода очень важен выбор функции , ее обычно берут вида . Где не меняет знака на отрезке, на котором ищется корень функции.
Положим b и рассмотрим метод в этом случае.
Тогда получим .
Числовые примеры
Рекомендации программисту
Заключение
Ссылки
Список литературы
- А.А.Самарский, А.В.Гулин. Численные методы. Москва «Наука», 1989.
- Н.Н.Калиткин. Численные методы. Москва «Наука», 1978.