- •Часть 2
- •Оглавление введение
- •Глава 1. Методы одномерной безусловной оптимизации § 1.1. Основные понятия
- •§ 1.2. Метод половинного деления
- •§ 1.3. Метод золотого сечения
- •§ 1.4. Метод Ньютона
- •Глава 2. Методы многомерной безусловной оптимизации § 2.1. Основные понятия
- •§ 2.2. Метод сопряженных направлений
- •§ 2.3. Метод наискорейшего спуска
- •§ 2.4. Метод Ньютона
- •Глава 3. Интерполирование и аппроксимация функций, заданных таблично § 3.1. Основные понятия
- •§ 3.2. Интерполяционный полином в форме Лагранжа
- •§ 3.3. Интерполяционный полином в форме Ньютона
- •§ 3.4. Сплайн-интерполяция таблично заданной функции
- •§ 3.5. Аппроксимация таблично заданных функций методом наименьших квадратов
- •Глава 4. Численное дифференцирование и интегрирование функций § 4.1. Основные понятия
- •§ 4.2. Численное дифференцирование
- •§ 4.3. Численное интегрирование
- •4.3.1. Квадратурные формулы прямоугольников, трапеций и Симпсона. Квадратурные формулы прямоугольников
- •Квадратурная формула трапеций
- •Квадратурная формула Симпсона
- •4.3.2. Автоматический выбор шага. Правило Рунге
- •Глава 5. Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений и их систем § 5.1. Основные понятия
- •§ 5.2. Метод Эйлера
- •§ 5.3. Явный 4-шаговый метод Адамса
- •§ 5.4. Метод Рунге Кутта четвертого порядка
- •Литература
4.3.1. Квадратурные формулы прямоугольников, трапеций и Симпсона. Квадратурные формулы прямоугольников
На каждом частичном отрезке разбиения [xk, xk + 1], k = 0, 1, 2, …, n – 1 воспользуемся квадратурной формулой вида
Ik f(xk)h, k = 0, 1, 2, …, n – 1. (4.3.9)
Тогда определенный интеграл от f(x) на всем отрезке [a, b] запишется в виде
(4.3.10)
Если на частичных отрезках воспользуемся квадратурной формулой
Ik f(xk + 1)h, k = 0, 1, 2, …, n – 1, (4.3.11)
получим следующую составную квадратурную формулу для отрезка [a, b]:
(4.3.12)
Наконец, в качестве квадратурной формулы на частичных отрезках разбиения можно использовать формулу
Ik f(xk + 1/2)h, k = 0, 1, 2, …, n – 1, (4.3.13)
где xk + ½ = (xk + xk + 1)/2 – середина k-го отрезка разбиения [xk, xk + 1], k = 0, 1, 2, …, n – 1.
В результате для отрезка [a, b] получим
(4.3.14)
Формулу (4.3.10) будем называть составной формулой левых прямоугольников, формулу (4.3.12) – составной формулой правых прямоугольников, формулу (4.3.14) – составной формулой центральных прямоугольников.
С геометрической точки зрения при использовании формул (4.3.9), (4.3.11) и (4.3.13) площадь под кривой y = f(x) заменяется на каждом частичном отрезке разбиения [xk, xk + 1] площадью прямоугольника шириной h и высотой, равной f(xk), f(xk + 1) или f(xk + 1/2) соответственно, чем и объясняется название этих формул.
Для квадратурных формул (4.3.9) и (4.3.11) – левых и правых прямоугольников – справедлива следующая оценка погрешности аппроксимации:
(4.3.15)
где
– максимум модуля первой производной
функции f(x)
на k-м частичном отрезке
разбиения [xk,
xk + 1],
k = 0, 1, 2, …, n
– 1.
Для составных формул (4.3.10) и (4.3.12) оценка погрешности аппроксимации:
(4.3.16)
где M1 – максимум модуля первой производной функции f(x) на отрезке [a, b].
Для формулы центральных прямоугольников (4.3.13) оценка погрешности аппроксимации
(4.3.17)
где
– максимум модуля второй производной
функции f(x)
на k-м частичном отрезке
разбиения [xk,
xk + 1],
k = 0, 1, 2, …, n
– 1.
Для составной формулы центральных прямоугольников
(4.3.18)
где M2 – максимум модуля второй производной функции f(x) на отрезке [a, b].
Считают, что формулы (4.3.10) и (4.3.12) имеют первый, формулы (4.3.9), (4.3.11) и (4.3.14) – второй, формула (4.3.13) – третий порядок точности по h.
Как следует из приведенных выше оценок погрешности аппроксимации, квадратурные формулы левых и правых прямоугольников точны для полиномов нулевого порядка, а квадратурная формула центральных прямоугольников – для полиномов до первого порядка включительно. Следует отметить, что вычисления с помощью квадратурной формулы центральных прямоугольников возможны, если значения f(xk + 1/2), k = 0, 1, 2, …, n – 1 либо могут быть вычислены, либо уже заданы в имеющейся таблице значений функции.
