y |
0 |
F(a |
0 |
). |
a x a |
0 |
|
F(a0) |
. |
|
|||||||||
|
|
|
|
|
F (a0) |
||||
y F(a0) F (a0)(x a0).
|
Рис. 5 |
|
|
|
|||||
|
Начало |
|
|
|
|||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Ввод a, |
|
|
|
|||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
i = 0 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
аi 1 аi аi |
|
|
i = i +1 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Нет |
|
|
F(аi |
|
) |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
||||||
|
|
|
|
|
Да |
|
|
|
|
Вывод |
|
|
|
а, i |
|
|
|
||
|
|
|
|
||||||
Конец
Рис. 6
26
Рабочая формула метода касательных: ai 1 ai FF((aaii)),
ai 1 ai ai, i 0,1,2,3,...
Итерационный процесс продолжается до тех пор, пока F(x) не
станет меньше заданного числа : |
F(ai ) |
. При работе с этим ме- |
тодом возможна потеря корня, но при правильном применении метода он сходится быстро, 4-5 итераций дают погрешность 10-5, он используется также для уточнения значения корня [5]. Блок-схема алгоритма метода касательных представлена на рис. 6, где an – корень уравнения; n – число итераций; F(an) - значение функции в соответствующей точке.
2.4. Комбинированный метод хорд и касательных
Задача. Найти корень уравнения F(x) с заданной точностью . В этом случае используется одновременно методы касательных и
хорд. Приближение к корню происходит с двух сторон. Рассмотрим четыре случая, которые отвечают возможным комбинациям знаков
F (x) и F (x).
Из графиков, представленных на рис. 7, метод хорд применяется со стороны вогнутости, а метод касательных – со стороны выпуклости графика.
Совместное применение обоих методов дает сразу избыточное и недостаточное приближение. Применяя этот метод, мы предполагаем,
что |
|
|
|
|
F(x), F (x) |
и F (x) |
непрерывны на отрезке [a0, b0], причем F (x) |
||
и F |
|
|
|
|
(x) сохраняют свой знак. Известно, что сохранение знака у F (x) |
||||
говорит о монотонности |
|
означает, |
||
F(x), а сохранение знака у F (x) |
||||
что выпуклость кривой y F(x) при всех x [a0 , b0 ] обращена в одну сторону. Для удобства расчета обозначим через а0 тот конец отрезка [a0, b0], в котором знаки F (x) и F (x) совпадают.
Из возможных случаев рассмотрим случай первый. Пусть F(a)* F(b) 0 и F (x)* F (x) 0, т.е. знаки первой и второй производной совпадают. При решении уравнения каждая итерация заключается в следующем: из точки А проведем хорду, которая стягивает дугу АВ, и проведем касательную к дуге таким образом, чтобы точка пересечения касательной с осью ох оказалось внутри отрезка [a0, b0]. Хорда на графиках пересекает ось ох в точке b1,
27
Рис. 7
лежащей между точками b0 и искомым корнем , а касательная к дуге в точке А пересекает ось ох в точке а0, лежащей между точками а0 и искомым корнем уравнения ( рис. 8).
Полученное значение a1 и b1 дают новое приближение к корню. Приведем расчетные формулы для ai+1и bi+1, выведенные в п.2.1 и 2.2.
b |
b |
F(bi )(a0 bi ) |
, |
a |
|
a |
|
|
F(ai ) |
. |
|
|
|
|
|||||||
i 1 |
i |
F(a0 ) F(bi ) |
|
i 1 |
|
i |
|
F (ai ) |
||
Процесс нахождения ai+1 и bi+1 продолжается до тех пор, пока выполняется одно из следующих условий:
ai 1 bi 1 , где - заданная точность;
F(bi ) |
|
или |
|
F(ai ) |
; |
||||||
F |
|
a |
i 1 |
b |
i 1 |
|
. |
||||
|
|
|
|
|
|||||||
|
|
|
|
2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||
28
Рис. 8
Все округления при вычислениях следует производить в сторону от корня [4]. На рис. 9 представлена блок-схема комбинированного метода хорд и касательных, где n число итераций; аn , b n− значения приближения корня; F(аn) F(bn) − значения функции в данных точках.
2.5. Метод простой итерации (метод последовательных приближений)
Чтобы применить метод простой итерации для решения нелинейного уравнения F(x)=0, необходимо преобразовать его к следующему
виду: |
|
x (x). |
(2) |
Это преобразование (приведение уравнения к виду, удобному для итерации) можно выполнить различными способами; некоторые из них будут рассмотрены ниже. Функция называется итерационной функцией.
Выберем каким-либо образом приближенное значение корня х(0) и подставим его в правую часть уравнения (2). Получим значение
х(1) (х(0) ). Подставим теперь х(1) в правую часть уравнения (2),
имеем х(2) (х(1) ). Продолжая этот процесс неограниченно, получим последовательность приближений к корню, вычисляемых по формуле
x(n 1) |
(x(n)), |
n 0. |
(3) |
|
29 |
|
|
Если существует предел построенной последовательности
х |
lim x(n) , то, переходя к пределу в равенстве (3) |
и предполагая |
|||||
|
n |
|
|||||
функцию непрерывной, получим равенство |
|
||||||
|
|
х |
х |
( |
х |
) |
(4) |
|
Это значит, что |
– корень уравнения (2). |
|
||||
Начало
Ввод a, b,
i = 0
аi 1 аi аi |
|
|
|||
|
|
|
|||
bi 1 bi bi |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
с |
a b |
|
|
|
|
|
|
i = i +1 |
|||
|
|||||
2 |
|
|
|
||
|
|
|
|
|
|
Нет
F(c)
Да Вывод c, i
Конец
Рис. 9
Метод допускает простую геометрическую интерпретацию. Построим графики функций у = x и у = (x), (рис. 10,а и 10,б). Корнем уравнения у = (x) является абсцисса точки пересечения кривой
30