
- •Раздел 6.
- •Раздел 6. Модели и алгоритмы решения задач численными методами с использованием математических пакетов Рекомендации по использованию учебного пособия
- •Тема 6.1. Элементы теории погрешностей
- •6.1.1. Точные и приближенные числа
- •6.1.2. Абсолютная и относительная погрешность
- •Тема 6.2. Методы решения нелинейных уравнений
- •6.2.1. Постановка задачи
- •Отделение корней (локализация корней);
- •Итерационное уточнение корней.
- •6.2.2. Отделение корней
- •6.2.2.1. Графическое отделение корней
- •6.2.2.2. Аналитическое отделение корней
- •6.2.3. Уточнение корней
- •6.2.3.1. Метод половинного деления
- •6.2.3.2. Метод итерации
- •6.2.3.3. Метод Ньютона (метод касательных)
- •6.2.3.4. Метод хорд
- •6.2.3.5. Сравнение методов решения нелинейных уравнений
- •6.2.4. Технология решения нелинейных уравнений средствами MathCad
- •Тема 6.3. Интерполяция функций
- •6.3.1. Постановка задачи
- •6.3.2. Интерполяционная формула Лагранжа
- •6.3.3. Интерполяционные формулы Ньютона
- •6.3.3.1. Конечные разности
- •6.3.3.2. Первая интерполяционная формула Ньютона
- •6.3.3.3. Вторая интерполяционная формула Ньютона
- •6.3.4. Сплайн – интерполяция
- •6.3.5. Сравнение интерполяционных многочленов по применению
- •6.3.6. Технология интерполяции функций в среде математических пакетов
- •Тема 6.4. Численное интегрирование
- •6.4.1. Постановка задачи
- •6.4.2. Метод прямоугольников
- •6.4.3. Формула трапеций
- •6.4.4. Формула Симпсона
- •6.4.5. Оценка погрешности численного интегрирования
- •6.4.6. Технология вычисления интегралов в среде математических пакетов
- •Тема 6.5. Методы решения обыкновенных дифференциальных уравнений
- •6.5.1. Постановка задачи
- •6.5.2. Метод Эйлера
- •6.5.3. Методы Рунге-Кутты
- •6.5.4. Решение оду n-го порядка
- •6.5.5. Сравнение методов решения оду
- •6.5.6. Технология решения обыкновенных дифференциальных уравнений средствами математических пакетов
- •6.6.2. Метод дихотомии
- •6.6.3. Метод золотого сечения
- •6.6.4. Сравнение методов
- •6.6.5. Технология решения задач одномерной оптимизации средствами математических пакетов
- •Тема 6.7. Аппроксимация функций
- •6.7.1. Постановка задачи аппроксимации
- •6.7.2. Метод наименьших квадратов
- •6.7.3. Технология решения задач аппроксимации функций средствами математических пакетов
- •Тема 6.8. Многомерная оптимизация
- •6.8.1. Постановка задачи и основные определения
- •6.8.2. Методы спуска
- •6.8.3. Метод градиентного спуска с дроблением шага
- •6.8.4. Метод наискорейшего спуска
- •6.8.5. Проблема оврагов. Метод покоординатного спуска
- •6.8.6. Технология решения задач многомерной оптимизации средствами математических пакетов
- •Список литературы
- •Тема 6.4. Численное интегрирование................................................71
- •Тема 6.5. Методы решения обыкновенных дифференциальных Уравнений............................................................................. 92
- •Тема 6.6. Одномерная оптимизация................................................ 115
- •Тема 6.7. Аппроксимация функций....................................................132
- •Тема 6.8. Методы многомерной оптимизации............................... 149
- •Список литературы.................................................................... 204
6.2.2. Отделение корней
Корень уравнения f(x)=0 считается отделенным (локализованным) на отрезке , если на этом отрезке данное уравнение не имеет других корней. Чтобы отделить корни уравнения, необходимо разбить область допустимых значений функции f(x) на достаточно узкие отрезки, в каждом из которых содержится только один корень. Существуют графический и аналитический способы отделения корней.
6.2.2.1. Графическое отделение корней
Графическое отделение корней основано на графическом способе решения уравнений – отыскании точек, в которых функция f(x) пересекает ось 0Х.
Пример 6.2.2-1. Отделить корни уравнения ln (x-1)2 – 0.5 = 0.
На
рис. 6.2.2-1 изображен график функции y
= ln
(x-1)2
– 0.5 , из которого следует, что
уравнение имеет два действительных
корня
[-1;0] и
[2;3].
Рис.6.2.2-1
В некоторых случаях удобно вначале преобразовать функцию f(x) к виду f(x)=g1(x) - g2(x), из которого, при условии f(x)=0, следует, что g1(x)=g2(x). При построении графиков y1=g1(x) и y2=g2(x) находят отрезки, содержащие точки пересечения этих графиков.
Пример 6.2.2-2. Отделить корни уравнения сos(x) – x + 1 = 0.
Приведем
исходное уравнение к виду сos(x)=
x
– 1. Построив графики функций
y1
= сos(x)
и y2
= х – 1 (рис. 6.2.2), выделим отрезок,
содержащий корень
[1;2].
Рис. 6.2.2-2
6.2.2.2. Аналитическое отделение корней
Аналитическое отделение корней основано на следующей теореме.
Если функция f(x) непрерывна и монотонна на отрезке [a;b] и принимает на концах отрезка значения разных знаков, то на отрезке [a;b] содержится один корень уравнения f(x)=0.
Действительно,
если условия теоремы выполнены, как это
имеет место на отрезке [a;b]
(рис. 6.2.2-3), то есть f(a)∙f(b)<0
и f'(x)>0
для xÎ
[a;b],
то график функции пересекает ось 0Х
только один раз и, следовательно, на
отрезке [a;b]
имеется один корень
уравнения f(x)
= 0. Аналогично можно доказать
единственность корня
на
отрезке [c;d],
на
[d;e]
и т.д
Рис. 6.2.2-3
Таким образом, для отделения корней нелинейного уравнения необходимо найти отрезки, в пределах которых функция монотонна и изменяет свой знак. Принимая во внимание, что непрерывная функция монотонна в интервалах между критическими точками, при аналитическом отделении корней уравнения можно рекомендовать следующий порядок действий:
установить область определения функции;
определить критические точки функции, решив уравнение f¢(x)=0;
составить таблицу знаков функции f(x) в критических точках и на границах области определения;
определить интервалы, на концах которых функция принимает разные знаки.
Пример 6.2.2-3. Отделить корни уравнения x - ln(x+2) = 0.
Область допустимых значений функции f(x) = x - ln(x+2) лежит в интервале (-2; ∞), найденных из условия x+2>0. Приравняв производную f¢(x)=1-1/(x+2) к нулю, найдем критическую точку хk= -1. Эти данные сведены в табл. 6.2.2-1 и табл. 6.2.2-2 знаков функции f(x).
Таблица 6.2.2-1 Таблица 6.2.2-.2
-
x
x→-2
-1
x→∞
x
-1.9
-1.1
-0.9
2.0
Sign(f(x))
+
-
+
Sign(f(x))
+
-
-
+
Уравнение x - ln(x+2) = 0 имеет два корня (-2;-1] и [-1; ∞) . Проверка знака функции внутри каждого из полученных полуинтервалов (табл.6.2.2) позволяет отделить корни уравнения на достаточно узких отрезках [-1.9;-1.1] и [-0.9;2.0].