ПК для ТЕБЯ

ПРОГРАММА РЕАЛИЗОВАНА В СРЕДЕ MATLAB.

Задача: Найти один корень заданного уравнения с относительной точностью 0.1%.. Рассчитать предельно-точное значение корня уравнения ?. В процессе решения уравнения методом простой итерации с погрешностью 10^-3 — 10^-5 вычислить отношение абсолютных погрешностей на соседних итерациях:

Теоретическая часть

Имеем уравнение f(x)=0, приведем его к удобному виду ?(x)=x, проделав, например, такие операции: f(x)=0 ? x+f(x)=x, где x+f(x)=?(x), получим, что ?(x)=x. Функция ?(x) называется итерационной функцией.

Для поиска корня следует выбрать начальное приближение x₀ такое, чтобы оно лежало как можно ближе к искомому корню ?. Последующие приближения вычисляются по итерационной формуле: x_n+1=?(x_n). Итерацией называется последовательность действий, в которой полученное значение на предыдущем шаге используется для поиска следующего значения на последующем шаге. Итерациями часто называют и сами полученные на каждых шагах значения x_n.

Теорема (о сходимости метода простой итерации): Пусть в некоторой  - окрестности корня ? функция ?(x) дифференцируема и удовлетворяет условию |?‘(x)|?q, где 0?q<1 — постоянная. Тогда независимо от выбора начального приближения из окрестности ?, итерационная последовательность не выходит из этой окрестности, метод сходится со скоростью геометрической прогрессии, и для метода справедлива оценка погрешности:

Критерий окончания итерационного процесса:

, где ? — заданная точность. Если величина 0<q?0,5, то очевидно, что при более простом критерии: |x_n - x_n-1|<?, метод будет сходиться быстрее.

Подготовка и решение задачи

Алгоритм:

Описание входной информации: xp — начальное приближение, e — относительная погрешность (по условию задачи e=0.001).

Описание выходной информации: xp — корень, который удовлетворяет условию, f(xp) — значение функции в найденном корне, n — количество проведенных приближений (итераций), ? — предельно-точное значение корня, Q — отношение абсолютных погрешностей, fip(ksi) — ?’(?), fip(xp) — ?’(x_n+1).

В программе ? определяется с помощью встроенной в MatLAB функции fzero(f,xp,eps). Где f- изначальная функция, xp-первое приближение, eps-точность (10^-6).

Ссылки для скачивания:
Метод простых итераций II способ
Вернуться назад