Методы приближенного решения уравнений. Метод последовательного поиска

Решаем приближенно уравнение F(x)=0 на отрезке [a,b] с точностью e. Пусть F(x) – непрерывная функция и на отрезке [a,b] выполняется достаточное условие существования корня . Если , то проверяем, какой из концов отрезка является корнем, в противном случае можно применить алгоритм последовательного поиска.

x=a; Fa=F(a);

do x=x+e; while F(x)*Fa>0;

Корнем объявляется x.

 

Методы приближенного решения уравнений. Метод вилки

Другое название метода – метод деления отрезка пополам. Решаем приближенно уравнение F(x)=0 на отрезке [a,b] с точностью e. Пусть F(x) – непрерывная функция и на отрезке [a,b] выполняется достаточное условие существования корня . Если , то проверяем, какой из концов отрезка является корнем, в противном случае можно применить алгоритм деления отрезка пополам.

Суть алгоритма заключается в следующем. Из условия задачи мы уже имеем решение задачи с точностью (в качестве решения можно взять любую точку отрезка [a,b], например точку a). Обозначим c=(a+b)/2 – середина отрезка [a,b], тогда возможны три варианта: либо F(c)=0 (тогда c – точное решение, в этом случае задача решена!), либо (тогда существует корень на отрезке [a,c], следовательно, мы имеем решение задачи с точностью ), либо, наконец, (тогда существует корень на отрезке [c,b], следовательно, мы получим решение задачи с точностью ).

Таким образом, алгоритм деления отрезка пополам позволяет за один шаг в два раза сократить область поиска решения уравнения. Применяя этот алгоритм раз, получаем решение уравнения с точностью e.

Методы приближенного решения уравнений. Метод простых итераций

Пусть – непрерывная и дифференцируемая (гладкая) функция. Требуется решить на с точностью уравнение . Заменим исходное уравнение на эквивалентное (то есть полученное из исходного применением инвариантных преобразований)

,

где . Обозначим левую часть нового соотношения , тогда уравнение будет выглядеть так: . Такая форма записи уравнения представляет собой итерационное соотношение. Возьмем наугад любое число из интервала – x0. С помощью итерационного соотношения можно получить последовательность чисел x1=Ф(x0); x2=Ф(x1); x3=Ф(x2) … и т.д. Последовательность может сходиться или расходиться. Если процесс сходится, то, очевидно, предельное значение последовательности является решением уравнения. Если полученное решение принадлежит отрезку , то это искомое решение.

Добиться сходимости можно, если подобрать такое , чтобы для некоторого числа выполнялось неравенство для всех x из .

 

 

 

Требуемая точность будет достигнута, как только

или

.