Вычислительная математика. Ч. 1. Асламова В.С - 7 стр.

UptoLike

13
1
1
1
1
2
2
2
2
y
1
y
2
Рис.1.
Теорема 1. Если непрерывная функция принимает значения разных знаков на
концах отрезка [a, b], т.е. f(a)*f(b) < 0, то внутри этого отрезка содержится
по меньшей мере один корень уравнения f(x)=0.
Корень заведомо будет единственным, если производная f '(x) существует и
сохраняет постоянный знак внутри интервала
(а, b).
Для отделения корней можно воспользоваться графическим путем. Для этого
представим уравнение
f(x)= 0 в виде x =
ϑ
(x), и строим графики функций у
1
= x
и y
2
=
ϑ
(x). Тогда абсциссы точек пересечения этих кривых и будут
действительными корнями уравнения
f(x) = 0.
Пример. Отделить корни уравнения
x
3
- 4x + 2 = 0.
Перепишем уравнение в виде
x(x
2
- 4) + 2 = 0; x = -2 / (x
2
- 4).
Построим графики функций
y
1
= x; y
2
= -2 / (x
2
- 4).
Абсциссы точек пересечения кривых
y
1
, y
2
дают приближенные значения
корней
x
2
0,5; x
2
1,6; x
3
-2,2 (рис.1).
Таким образом, функция имеет три корня, каждый из которых нужно
уточнить выбранным численным методом. Для уточнения корня выбираем
поочерёдно три интервала
x
1
[-2.5; -2.05], x
2
[0; 0.8], x
3
[1.2; 1.8].
14
2.2. Метод половинного деления
(метод бисекции, метод дихотомии)
Этот метод требует непрерывности функции
f(x), кроме того, на концах
отрезка [
a, b] функция должна иметь разные знаки.
В качестве начального приближения принимаем середину отрезка
2
ba
c
+
=
.
Для определения отрезка, содержащего корень, проверяем условие:
f(a)
f(с)<0 (2.1)
Можно сравнивать значения функции на концах отрезка [
с, b].
Если условие (2.1) выполняется, то корень лежит в отрезке [
a, c], поэтому
нужно передвинуть правую границу
b в точку c , т.е. b=c.
Если условие (2.1) не выполняется, то корень лежит в отрезке [
с, b] и следует
поменять левую границу
a=c.
Проверяем условие выхода из цикла: «длина отрезка меньше или равна
заданной точности?»:
|b-a|≤
ε
(2.2)
Если условие (2.2) не выполняется, то снова делим отрезок пополам,
определяем отрезок, содержащий корень и т.д. После каждой итерации отрезок,
содержащий корень уменьшается вдвое, т.е. после
n итераций он сократится в 2
n
раз. На рис.2 приведена графическая интерпретация метода половинного
деления, а на рис.3 – блок-схема метода половинного деления.
Метод половинного деления всегда сходится, то есть при его использовании
решение получается всегда, причём с заданной точностью. Требуемое обычно
большее число итераций по сравнению с некоторыми другими методами не
является препятствием к применению этого мет
ода, если каждое вычисление
значения функции
f(x) несложно.
a
b
f(x)
f(a)
f(b)
0
f(x)
Рис.2 Геометрическая интерпретация метода
половинного деления
c
f(c
0
)
b
c
c
f(c
1
)
х
    Теорема 1. Если непрерывная функция принимает значения разных знаков на                         2.2. Метод половинного деления
концах отрезка [a, b], т.е. f(a)*f(b) < 0, то внутри этого отрезка содержится
по меньшей мере один корень уравнения f(x)=0.                                                     (метод бисекции, метод дихотомии)
   Корень заведомо будет единственным, если производная f '(x) существует и
сохраняет постоянный знак внутри интервала (а, b).
   Для отделения корней можно воспользоваться графическим путем. Для этого         Этот метод требует непрерывности функции f(x), кроме того, на концах
                                                                                отрезка [a, b] функция должна иметь разные знаки.
представим уравнение f(x)= 0 в виде x = ϑ(x), и строим графики функций у1 = x
                                                                                                                                                           a+b
и y2 = ϑ(x). Тогда абсциссы точек пересечения этих кривых и будут                  В качестве начального приближения принимаем середину отрезка c =              .
действительными корнями уравнения f(x) = 0.                                                                                                                   2
     Пример. Отделить корни уравнения                                           Для определения отрезка, содержащего корень, проверяем условие:
                                   x3 - 4x + 2 = 0.                                                             f(a)⋅f(с)<0                                (2.1)
Перепишем уравнение в виде                                                         Можно сравнивать значения функции на концах отрезка [с, b].
                        x(x2 - 4) + 2 = 0; x = -2 / (x2 - 4).                      Если условие (2.1) выполняется, то корень лежит в отрезке [a, c], поэтому
                                                                                нужно передвинуть правую границу b в точку c , т.е. b=c.
Построим графики функций                                                           Если условие (2.1) не выполняется, то корень лежит в отрезке [с, b] и следует
                             y1 = x;    y2 = -2 / (x2 - 4).                     поменять левую границу a=c.
   Абсциссы точек пересечения кривых y1, y2 дают приближенные значения             Проверяем условие выхода из цикла: «длина отрезка меньше или равна
корней x2 ≈ 0,5; x2 ≈1,6; x3 ≈ -2,2 (рис.1).                                    заданной точности?»:
                                                                                                                  |b-a|≤ε                                    (2.2)
                                                      y1
                                                                                                       f(x)
                                                                                                                                     f(x)
                                                                                                       f(b)
                                                                                                      f(c1 )
                                        2                                                                                                       х
                                                                                                               a          c
                                        1                                                                0                     c c          b
                           2    1            1   2                                                    f(c0 )
                                        1
                                                                                                       f(a)
                                        2
                                                                                                   Рис.2 Геометрическая интерпретация метода
                               y2                                                                        половинного деления
                                                                                   Если условие (2.2) не выполняется, то снова делим отрезок пополам,
                                    Рис.1.                                      определяем отрезок, содержащий корень и т.д. После каждой итерации отрезок,
                                                                                содержащий корень уменьшается вдвое, т.е. после n итераций он сократится в 2n
   Таким образом, функция имеет три корня, каждый из которых нужно              раз. На рис.2 приведена графическая интерпретация метода половинного
уточнить выбранным численным методом. Для уточнения корня выбираем              деления, а на рис.3 – блок-схема метода половинного деления.
поочерёдно три интервала x1∈[-2.5; -2.05], x2∈[0; 0.8], x3∈[1.2; 1.8].             Метод половинного деления всегда сходится, то есть при его использовании
                                                                                решение получается всегда, причём с заданной точностью. Требуемое обычно
                                                                                большее число итераций по сравнению с некоторыми другими методами не
                                                                                является препятствием к применению этого метода, если каждое вычисление
                                                                                значения функции f(x) несложно.
                                        13                                                                             14