- •Численные методы
- •Содержание
- •Введение
- •Тема №1. Решение уравнений
- •1. Основные сведения
- •1.1 Отделение корней уравнения
- •Решение
- •1.2 Нахождение корней уравнения
- •Решение
- •Решение
- •Метод касательных (метод Ньютона)
- •2. Лабораторная работа №1 «Решение уравнений».
- •3. Контрольные вопросы
- •Тема №2. Решение систем линейных уравнений
- •1. Основные сведения
- •1.1 Прямые методы
- •Решение
- •Решение
- •1.2 Итерационные методы Метод простой итерации
- •Решение
- •Метод Зейделя
- •2. Лабораторная работа №2 «Решение систем линейных уравнений»
- •3. Контрольные вопросы
- •Тема №3. Интерполяция и аппроксимация функций
- •Основные сведения
- •Интерполяция
- •1.1.1. Глобальная полиномиальная интерполяция
- •Интерполяционный полином Лагранжа
- •Решение
- •Интерполяционный многочлен Ньютона
- •Решение
- •Решение
- •Локальная (кусочно-полиномиальная) интерполяция
- •Решение
- •Решение
- •Интерполяция сплайнами
- •1.2. Аппроксимация. Метод наименьших квадратов
- •Решение
- •2. Лабораторная работа № 3 «Интерполяция и аппроксимация функций»
- •3. Контрольные вопросы
- •3.1. Интерполяция
- •Аппроксимация
- •Тема №4. Дифференцирование функций
- •Основные сведения
- •Вычисление первой производной
- •Вычисление второй производной
- •Дифференцирование с использованием интерполяционных формул Дифференцирование на основе интерполяционной формулы Ньютона
- •Решение
- •Дифференцирование на основе многочлена Лагранжа
- •Лабораторная работа №4 «Дифференцирование функций».
- •Контрольные вопросы.
- •Тема №5 Интегрирование функций
- •Основные сведения
- •1.1. Общие принципы численного интегрирования.
- •Простейшие квадратурные формулы прямоугольников, трапеций и Симпсона.
- •Формулы прямоугольников
- •Формула трапеций
- •Формула Симпсона (формула парабол)
- •Решение
- •1.3. Квадратурные формулы Ньютона – Котеса.
- •1.4. Квадратурные формулы Гаусса.
- •Решение
- •2. Лабораторная работа №5 «Интегрирование функций».
- •3. Контрольные вопросы.
- •Тема №6 Решение обыкновенных дифференциальных уравнений
- •1. Основные сведения.
- •1.1. Понятие о численном решении задачи Коши для дифференциального уравнения первого порядка.
- •1.2. Метод Эйлера. Модифицированный метод Эйлера.
- •Решение
- •1.3. Методы Рунге - Кутта.
- •Решение
- •Лабораторная работа №6 «Решение обыкновенных дифференциальных уравнений».
- •3. Контрольные вопросы. Метод Эйлера
- •Решение
- •Решение
Решение
Корни данного уравнения локализованы на отрезках [-3;-2], [-1;0], [0;1].
На отрезке
преобразуем уравнение к виду
.
Максимальное
значение производной функции
на
отрезке
равно
:
,
а значит,
,
и итерационный процесс будет сходящимся.
На отрезках [0;1] и [-1;0] производная функции неограниченна, и указанное преобразование уравнения (1.1) не подходит.
На отрезках [0;1] и [-1;0] приемлемы следующие преобразования исходного уравнения к виду (1.4):
на
[-1;0] и
на [0;1] .
Одним из способов приведения уравнения (1.1) к виду , допускающему решение его методом простой итерации, является его преобразование в форме
(1.7)
где m — отличная от нуля константа.
В этом случае
.
Дифференцируя,
получим
.
Далее подбираем m так, чтобы выполнялось условие
на
всем промежутке изменения переменной
.
Из последнего неравенства получаем:
.
Если
для всех x
и
,
то
.
Если
и
,
то
.
Значение q, входящее в (1.5) , определяется из условия
.
Так,
если
принять
,
то для
; если
,
то для
.
Пример №4. Привести уравнение
к итерационному виду (1.4).
Решение
Корень данного уравнения принадлежит отрезку [1;1,5] (см. пример 1).
Можно преобразовать уравнение к виду
.
Однако
расчеты показывают, что значения
производной функции
на отрезке [1;1,5] по абсолютной величине
больше двух:
,
и, следовательно, итерационный процесс не будет сходящимся.
Преобразуем данное уравнение к виду
.
В этом случае
и
.
Подберем m так, чтобы значение производной по абсолютной величине было меньше единицы.
Так
как производная
отрицательна и ее наименьшее значение
на правом конце отрезка
,
то можно принять
и рассматривать уравнение
,
пригодное для решения методом простой итерации.
Метод касательных (метод Ньютона)
Рис. №6 иллюстрирует сущность метода Ньютона.
Пусть - первое приближение корня
уравнения (1.1):
.
Проведем в точке
((
,
)
касательную к кривой
.
За второе приближение корня
примем точку пересечения этой касательной
с осью абсцисс.
Аналогично получим следующее приближение корня, используя касательную к кривой, проведенную в точке ( ,
).В результате получим последовательной точек
,
,…,
,…
для которых справедливо
следующее утверждение.
Пусть
-
простой корень уравнения
,
в некоторой окрестности которого функция
дважды непрерывно дифференцируема.
Тогда
найдется такая малая
-окрестность
корня
,
что при произвольном выборе начального
приближения
из этой окрестности итерационная
последовательность метода Ньютона
, ,…, ,…
не выходит за пределы этой -окрестности корня, и справедлива оценка
.
(1.8)
Таким образом, при выборе начального приближения из достаточно малой окрестности корня метод Ньютона сходится квадратично.
Это означает, что на каждой итерации число верных цифр приближения примерно удваивается.
Так
как уравнение касательной к кривой
в точке (
,
)
имеет вид
,
то,
полагая в нем
,
получим расчетную
формулу метода Ньютона
(если
):
.
(1.9)
Используя (1.9), можно получить оценку погрешности метода Ньютона. На практике обычно используется следующая простая оценка:
.
(
1.10)
