- •И. Н. Двойцова вычислительная математика
- •Введение
- •1. Элементарная теория погрешностей; вычислительные задачи, методы и алгоритмы
- •1.1. Источники и классификация погрешностей результата численного эксперимента
- •1.2. Погрешности чисел
- •1.3. Погрешности арифметических операций
- •1.4. Погрешности функций
- •1.5. Особенности машинной арифметики
- •1.6. Лабораторная работа № 1. Определение абсолютной и относительной погрешностей приближенных чисел. Оценка погрешностей результата
- •1.7. Корректность вычислительной задачи
- •1.8. Обусловленность вычислительной задачи
- •1.9. Вычислительные методы, их классификация
- •2. Приближение функций
- •2.1. Задача приближения функций
- •2.2. Интерполяция обобщенными многочленами
- •2.3. Полиномиальная интерполяция. Многочлен Лагранжа
- •2.4. Погрешность интерполяции
- •2.5. Конечные разности и их свойства
- •Доказательство
- •2.6. Разделенные разности и их свойства
- •2.7. Интерполяционный многочлен Ньютона
- •2.8. Вычислительная схема Эйткена
- •2.9. Лабораторная работа № 2. Интерполирование и экстраполирование данных. Интерполяционный многочлен Лагранжа
- •2.10. Интерполяционный многочлен Ньютона с конечными разностями
- •2.11. Лабораторная работа № 3. Интерполирование и экстраполирование данных. Интерполяционный многочлен Ньютона
- •3. Метод наименьших квадратов и специальные интерполяционные многочлены
- •3.1. Постановка задачи и вывод формул метода наименьших квадратов
- •3.2. Лабораторная работа № 4. Аппроксимация функции по методу наименьших квадратов
- •3.3. Глобальная полиномиальная интерполяция
- •3.4. Чувствительность интерполяционного многочлена к погрешностям входных данных
- •4. Численное дифференцирование и численное интегрирование
- •4.1. Простейшие формулы численного дифференцирования для первой производной
- •4.2. Формулы численного дифференцирования для второй производной
- •4.3. Формулы численного дифференцирования, основанные на интерполяции алгебраическими многочленами
- •4.4. Обусловленность формул численного дифференцирования
- •4.5. Простейшие квадратурные методы численного интегрирования
- •4.6. Оценка погрешностей простейших квадратурных формул
- •4.7. Квадратурные формулы интерполяционного типа
- •4.8. Квадратурные формулы Гаусса
- •4.9. Лабораторная работа № 8. Численное дифференцирование и численное интегрирование функций
- •Ирина Николаевна Двойцова Вычислительная математика
- •660049, Красноярск, пр. Мира , 82
4.5. Простейшие квадратурные методы численного интегрирования
В прикладных исследованиях, когда
возникает необходимость вычисления
и первообразной не существует, приходиться
интеграл считать численно. Наиболее
широко на практике используются
квадратурные формулы - приближенные
равенства вида
(4.5.1)
где
- некоторые точки из отрезка
- узлы квадратурной формулы (4.5.1),
- числовые коэффициенты, называемые
весами квадратурной формулы,
- целое число. Сумма
,
которая принимается за приближенное
значение интеграла, называется
квадратурной суммой.
Величина
называется погрешностью или остаточным
членом квадратурной формулы. Выведем
простейшие квадратурные формулы, исходя
из геометрической интерпретации
определенного интеграла:
...
Будем интерпретировать интеграл
как площадь криволинейной трапеции,
ограниченной графиком функции
,
осью абсцисс и прямыми
Разобьем отрезок
на элементарные отрезки
точками
Интеграл
представится таким образом:
(4.5.2)
Введем обозначения:
Формула
центральных прямоугольников. Заменим
приближенно площадь элементарной
криволинейной трапеции площадью
прямоугольника, основанием которого
является отрезок
,
а высота равна значению
.
Тогда сразу получается элементарная
квадратурная формула прямоугольников
Суммируя
по всему отрезку
,
получим
(4.5.3)
Совершенно аналогично можно получить
формулы
и
,
которые называются квадратурными
формулами левых и правых прямоугольников.
Их точность
,
тогда как точность формулы (4.5.3)
.
Ф
ормула
трапеций. Соединив точки
,
получим формулу трапеций. Заменим
площадь элементарной криволинейной
трапеции площадью построенной фигуры.
Тогда
,
а итоговая формула примет вид
(4.5.4)
Ф
ормула
парабол (Симпсона).
Если площадь элементарной криволинейной
трапеции заменить площадью фигуры,
расположенной под параболой, проходящей
через точки
,
то получим приближенное равенство
- интерполяционный многочлен второй
степени с узлами
Для этих точек справедлива формула
. (4.5.5)
Действительно, легко проверить, что
,
,
.
Кроме того, формула (4.5.5) относительно
представляет уравнение второй степени.
Подставим теперь точки
,
,
в уравнение (4.5.5). Получим
Таким
образом, точки
,
,
параболы
удовлетворяют уравнению (4.5.5).
Проинтегрируем полученную формулу по
отрезку
.
Тогда
.
Применяя полученную формулу на каждом элементарном отрезке, получим квадратурную формулу парабол (Симпсона)
(4.5.6)
