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

Поиск:

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

⇐ Предыдущая 1 2 345 6 7 8 9 10 Следующая ⇒

Исходное уравнение (3.1) преобразуем к эквивалентному уравнению:

x =

(x).

(3.8)

Пусть известно начальное приближение (полученное, например, на этапе отделения корней): x = x ₀. Подставим его в правую часть (3.8) и получим новое приближение: x₁ = (x₀). Повторяя эту процедуру, будем иметь в общем виде на некотором k-м шаге:

x_k = (x_k-1) .

В качестве условия окончания вычислительного процесса можно взять выполнение неравенства: ½x_k - x_k-1½ < .

Значение x_k, удовлетворяющее ему, и есть корень уравнения F(x, a₁, a₂, ..., a_k) = 0.

Геометрическая интерпретация этого метода приведена на рис.3.8, 3.9. Здесь x^* - истинное, искомое значение корня; x₀ - начальное приближение к корню; x₁, x₂, x₃ - очередные итерации.

Рис.3.8.

Рис.3.9.

При использовании этого метода возникает вопрос о его сходимости. Дело в том, что при некоторых условиях расстояние между истинным корнем и приближениями к нему может возрастать с каждой новой итерацией, как это показано на рис.3.10, 3.11.

Рис.3.10.

Рис.3.11.

Условием сходимости метода простых итераций является выполнение в окрестности искомого корня неравенства:

(x)½ < 1

(3.9)

Это условие является достаточным, т.е. если оно выполняется, то процесс обязательно сходится; если же условие (3.9) не выполняется или выполняется не во всех точках

x₀, x₁, x₂, ..., x_k, ... ,

то заранее сказать что-либо конкретное о сходимости нельзя.

Итак, для решения уравнения F(x) = 0методом простых итераций надо преобразовать его к уравнению вида x = (x) так, чтобы выполнялось условие ½ (x)½ < 1. Сходимость к истинному корню будет тем быстрее, чем ближе к единице значение (x).

Существует более или менее универсальный способ преобразования уравнения F(x, a₁, a₂, ..., a_k) = 0к виду x = (x):

F(x) = 0

C ^. F(x) = 0

C ^. F(x) + x = x

(3.15)

Здесь C - некоторый параметр, выбираемый из условия сходимости процесса.

При использовании преобразования (3.15) условием окончания вычислительного процесса является выполнение неравенства

В программе необходимо указывать функцию F(x) и вводить вычисленный заранее параметрС и значение допустимой погрешности . Программа должна осуществлять не более 100 итераций. Если за 100 итераций не достигнута требуемая точность, то программа выводит сообщение об отсутствии сходимости и прекращает работу.

⇐ Предыдущая 1 2 345 6 7 8 9 10 Следующая ⇒

Дата добавления: 2015-02-12; просмотров: 1223 | Нарушение авторских прав | Изречения для студентов