
- •Методические указания
- •Часть 1
- •1. Численное решение нелинейных уравнений.
- •1.1. Метод деления отрезка пополам.
- •1.2. Метод Ньютона (метод касательных).
- •1.3. Метод простой итерации.
- •2. Методы решения систем линейных алгебраических уравнений.
- •2.1. Метод Гаусса.
- •2.2. Метод обратной матрицы.
- •2.3. Метод прогонки.
- •2.4. Метод простой итерации (метод Якоби).
- •2.5. Метод Зейделя.
- •3. Численные методы решения систем нелинейных уравнений.
- •3.1. Метод простой итерации (метод Якоби).
- •3.2. Метод Зейделя.
- •3.3. Метод Ньютона.
- •Литература
1.2. Метод Ньютона (метод касательных).
Суть метода состоит
в том, что на
-й
итерации в точке
строится касательная к кривой
и ищется точка пересечения касательной
с осью абсцисс (рис. 1.6).
Если задан интервал изоляции корня
,
то за начальное приближение
принимается тот конец отрезка, на котором
. (1.1)
Уравнение
касательной, проведенной к кривой
в точке
с координатами
и
,
имеет вид:
(1.2)
|
Рис. 1.6. Метод касательных. |
За следующее
приближение корня
примем абсциссу точки пересечения
касательной с ocьюOX.
Из (1.2) при
,
получим
(1.3)
При этом необходимо,
чтобы
.
Аналогично могут
быть найдены и следующие приближения
как точки пересечения с осью абсцисс
касательных, проведенных в точках
,
и т.д. Формула для
-го
приближения имеет вид:
(1.4)
Для завершения
итерационного процесса можно использовать
условия
или
.
Объем вычислений
в методе Ньютона больше, чем в других
методах, поскольку приходится находить
значение не только функции
,
но и ее производной. Однако скорость
сходимости здесь значительно выше.
Пример 1.2.Решить уравнениена отрезке
методом Ньютона c точностью
.
Решение. Определим
первые и вторые производные заданной
функции
:
;
.
Проверим выполнение условия сходимости
на концах заданного интервала:
- не выполняется,
- выполняется. За начальное приближение
корня можно принять
.
Находим первое приближение:
.
Аналогично находится второе приближение:
.
Третье приближение:
.
Так как
,
итерационный процесс заканчивается.
Таким образом, приближенным решением
данного уравнения является
.
На рис. 1.7 приведена программа решения данного уравнения методом Ньютона. В качестве исходных данных вводятся начальное приближение и точность вычисления.
CLS DEF FNF(X)=X^3+X-1 DEF FNP(X)=3*X+1 INPUT X, E 1 X=X-FNF(X)/FNP(X) PRINT X, FNF(X) IF ABS(FNF(X)/FNP(X))>E THEN 1 END
|
Рис. 1.7. Программа нахождения корней методом Ньютона на языке QUICK BASIC. |
Пример 1.3.Решить уравнениена отрезке
методом Ньютона c точностью
с помощью программыExcel.
|
A |
B |
C |
D |
1 |
x |
F(x) |
F'(x) |
погрешность |
2 |
1,00000 |
|
|
|
3 |
0,75000 |
1,00000 |
4,00000 |
0,25000 |
4 |
0,68605 |
0,17188 |
2,68750 |
0,06395 |
5 |
0,68234 |
0,00894 |
2,41198 |
0,00371 |
6 |
0,68233 |
0,00003 |
2,39676 |
0,00001 |
Рис. 1.8. Решение уравнения методом Ньютона с помощью программы Excel. |
Порядок решения (рис. 1.8).
Ввести в ячейки A1:D1 заголовки столбцов.
В ячейку A2 – значение начального приближения
В ячейку B3 – формулу функции =A2^3+A2-1
В ячейку C3 – формулу производной функции =3*A2^2+1
В ячейку A3 – формулу первого приближения =A2-B3/C3
В ячейку D3 – погрешность =ABS(A3-A2)
Выделить ячейки A3:D3 и скопировать формулы в соседние ячейки расположенных ниже строк A4:D4, A5:D5, и т.д. при помощи маркера заполнения. Каждая новая строка содержит результаты очередного приближения.
В столбце A найти значение корня, соответствующее заданной точности.
Приближенное
решение данного уравнения
содержится в ячейкеA6
(погрешность
в ячейкеD6).