- •Глава 1. Интерполирование
- •Постановка задачи интерполирования
- •1.2. Единственность интерполяционного многочлена
- •. Интерполяционный многочлен Лагранжа
- •1.4. Погрешность интерполяционной формулы Лагранжа
- •. Конечные разности и их свойства
- •1.6. Первая интерполяционная формула Ньютона
- •1.7. Вторая интерполяционная формула Ньютона
- •1.8. Интерполяционная формула Гаусса
- •1.10. Численное дифференцирование
- •Глава 2. Численное интегрирование
- •2.1. Общие замечания
- •Глава 3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1. Вводные замечания
- •3.2. Отделение корней
- •3.3. Метод половинного деления
- •3.5. Метод Ньютона (метод касательных)
- •3.6. Комбинированный метод
- •3.7. Метод итерации (метод последовательных приближений)
- •Глава 4. Решение систем линейных уравнений
- •4.1. Метод Гаусса (схема единственного деления)
- •4.2. Вычисление определителей
- •4.3. Вычисление обратной матрицы
- •4.4. Некоторые сведения из линейной алгебры
- •Рассмотрим систему линейных алгебраических уравнений
- •Вычитая второе равенство из первого, получаем
- •4.6. Приведение линейной системы к виду,
- •4.7. Метод Зейделя Пусть система линейных алгебраических уравнений приведена к виду
- •Глава 5. Обработка результатов наблюдений
- •5.1. Постановка задачи
- •5.2. Приближение функции, заданной таблично,
- •Глава 6. Численные методы решения
- •6.1. Методы решения задачи Коши. Вводные замечания
- •6.2. Решение дифференциальных уравнений
- •6.3. Метод Эйлера Пусть дано дифференциальное уравнение
- •6.4. Уточненный метод Эйлера Рассмотрим дифференциальное уравнение
- •6.6. Методы Рунге – Кутта Пусть дано дифференциальное уравнение
- •6.7. Метод Милна
- •6.8. Приближенное решение систем дифференциальных
- •6.9.Краевые задачи для обыкновенных дифференциальных
- •6.10. Решение краевых задач
- •6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
- •Глава 7. Численное решение дифференциальных уравнений в частных производных
- •7.1. Классификация дифференциальных уравнений
- •7.2. Метод сеток решения краевых задач
- •7.3. Погрешность аппроксимации дифференциальных уравнений
- •7.4. Аппроксимация граничных условий
- •7.6. Метод сеток для уравнений
3.5. Метод Ньютона (метод касательных)
Дано уравнение f(x)=0. Пусть найден отрезок , такой, что на концах его функция f(x) имеет разные знаки, т.е. . Пусть, кроме того, производные и на отрезке сохраняют знак. Пусть для определенности , , , при .
В точке проведем касательную к кривой y=f(x) (рис. 3.3).
У В0
0 а0 b2 b0 B1
x
A0
Рис. 3.3
Точку пересечения касательной с осью 0х принимаем за приближенное значение корня уравнения. Запишем уравнение касательной
.
Точку пересечения касательной с осью 0х обозначим через b1. Имеем
.
Принимая за конец нового отрезка можно повторить предыдущий шаг и найти :
и т.д.
. (3.4)
Докажем сходимость этого процесса. Рассмотрим последовательность
b0, b1, b2, …, bn, ….
Докажем, что все . Используем метод математической индукции. Прежде всего .
Пусть . Докажем ,что .
Положим
Применяя формулу Тейлора, получаем
.
где .
Так как и , получаем
.
Отсюда
, т.е. .
Доказано, что все и следовательно, . Из формулы (3.4) вытекает, что . Итак , последовательность b0, b1, b2, … монотонно убывающая, ограниченная.
Существует . Покажем, что с – это и есть корень уравнения . Действительно, переходя к пределу при в равенстве (3.4), получаем
.
Отсюда , т.е. .
Отметим, что если провести касательную к кривой y=f(x) в точке (рис. 3.3), получим точку , лежащую вне отрезка . Поэтому, применяя метод касательных,следует руководствоваться следующим правилом: касательная проводится на том же конце отрезка, где знаки функции и второй производной совпадают.
3.6. Комбинированный метод
Метод хорд и метод касательных дают приближение к корню с какого-либо одного конца отрезка; второй конец отрезка остается неподвижным. На практике удобнее использовать комбинированный метод, заключающийся в поочередном применении метода хорд и метода касательных.
Возьмем для определенности , , , при . Определим а1 по методу хорд и по методу касательных. Затем находим , применяя метод хорд на отрезке , и т.д. (рис. 3.4).
Комбинированный метод дает более быструю сходимость, чем метод хорд или метод касательных в отдельности. Кроме того, при применении комбинированного метода легко оценить погрешность результата, так как и находятся по разные стороны от корня. Отсюда, в частности, следует, что цифры, совпадающие у и , принадлежат точному корню .
у
0 а0 а1 а2
x
b2 b1 b0
Рис. 3.4
Четыре возможные комбинации знаков производных и определяют четыре типа расположения кривой y=f(x). условимся через b0 обозначать тот конец отрезка, на котором знак функции f(x) и ее второй производной совпадают (рис. 3.5).
Расчетные формулы комбинированного метода имеют вид:
где n=0, 1, 2,… .
Если корень уравнения требуется вычислить с точностью до , процесс вычисления корня можно прекращать в тот момент, когда . в ответа возьмем среднее арифметическое полученных значений аn и bn, т.е.
.
y
0 a0 b0 x 0 b0 a0 x
y
0 b0 a0 x 0 a0 b0 x
Рис. 3.5
Пример. Вычислить с точностью до действительный корень уравнения . вычисления в виде табл.3.2.