
- •И. Н. Двойцова вычислительная математика
- •Введение
- •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.6. Оценка погрешностей простейших квадратурных формул
Теорема 4.1. Пусть функция дважды непрерывно дифференцируема на . Тогда для формул (4.5.3) и (4.5.4) справедливы следующие оценки погрешности:
(4.6.1)
(4.6.2)
Доказательство
Получим, например, формулу (4.6.1). По определению погрешности
С другой стороны, по формуле Тейлора
Тогда
Но
следовательно,
Аналогично выводится и формула (4.6.2).
Теорема 4.2. Пусть функция имеет на непрерывную четвертую производную. Тогда для формулы (4.5.6) справедлива следующая оценка:
(4.6.3)
Пример.
Вычислить интеграл
по формулам прямоугольников, трапеций
и парабол и оценить погрешность
вычислений.
-
числитель
знаменат.
0
0.50
1.100
0.453596
1.050
3.167423
0.143207
1/2
0.55
1.121
0.434782
1.125
3.202268
0.135773
1
0.60
1.144
0.413957
1.200
3.232039
0.128079
3/2
0.65
1.169
0.391072
1.275
3.256570
0.120088
2
0.70
1.196
0.366083
1.350
3.275723
0.111756
5/2
0.75
1.225
0.338946
1.425
3.289391
0.103042
3
0.80
1.256
0.309623
1.500
3.297495
0.093896
7/2
0.85
1.272
0.294370
1.575
3.299991
0.089203
4
0.90
1.324
0.244299
1.650
3.296865
0.074100
9/2
0.95
1.361
0.208261
1.725
3.288134
0.063337
5
1.00
1.400
0.169967
1.800
3.273848
0.051917
0.6 0.7 0.8 0.9 1.0
Оценим теперь погрешность вычислений только по теореме 4.1:
Вторая производная уже получается
настолько громоздкой, что оценка
погрешности становится затруднительной.
Тем не менее можно провести грубую
оценку сверху
Итак,
Так оно и есть, поскольку значения интеграла, полученные по формулам прямоугольников и трапеций, отличаются в четвертом знаке после запятой.