|
чЯч. Учебнометодическое пособие по курсу Методы вычислений для студентов ivv курсов всех форм обучения
Многочлены Якоби (23) ортогональны на отрезке с весом и для любого многочлена степени меньшей
С помощью теоремы 3 получаем, что алгебраический порядок точности квадратурной формулы (22) равен . Квадратурная формула (22) содержит два параметра и , из нее могут быть получены специализированные квадратурные формулы, соответствующие распространенным видам степенных особенностей (см. [3]). В справочниках приведены квадратурные формулы Гаусса с другими весами.
1.7. Локально-интерполяционные (составные) квадратурные формулы
Для повышения точности квадратурных формул используют прием, идея которого восходит к римановым интегральным суммам. Отрезок интегрирования разбивают на некоторое число частичных отрезков, на каждом из которых применяют квадратурную формулу с небольшим числом узлов. В качестве параметра квадратурного процесса теперь используют число частичных отрезков.
Разобьем отрезок на частичных отрезков точками . Для вычисления интеграла на каждом частичном отрезке применим интерполяционную квадратурную формулу с узлами и
коэффициентами . Получим квадратурную формулу
Пусть и , тогда для погрешности квадратурной формулы (24) имеет место оценка
Квадратурная формула (24) называется локально-интерполяционной или составной.
Наиболее часто формула (24) используется в случае, когда отрезок разбит на частичные отрезки равной длины и на каждом частичном отрезке используется квадратурная формула Ньютона-Котеса с узлами. Из (24) получаем при , и локально-интерполяционную квадратурную формулу
где
Сумма абсолютных величин коэффициентов формулы (26)
не зависит от числа частичных отрезков .
Оценка погрешности квадратурной формулы (25) имеет вид
Из теоремы 2 и (27) следует, что квадратурный процесс, порожденный локально-интерполяционной квадратурной формулой (25), является сходящимся при (со скоростью на функциях из класса ).
Приведем простейшие составные квадратурные формулы, часто применяемые в практике. Правило трапеций
Если , то Правило Симпсона (парабол)
Если , то Замечание 5. Алгоритмы численного интегрирования, построенные на основе локально-интерполяционных квадратурных формул (25) имеют существенный недостаток – они насыщаемые. Насыщаемость проявляется в том, что асимптотическое представление погрешности формулы (25) имеет главный член. Отсюда следует неулучшаемость оценки погрешности, сколь бы ни была гладкой функция .
В зависимости от гладкости функции можно выписать любое заданное число членов асимптотического ряда, в который разлагается погрешность . Рассмотрим конкретный пример – правило трапеций. Если , то для погрешности квадратурной формулы имеет место представление
где и не зависит от . Из (28) и следует насыщаемость правила трапеций. Классом насыщения в данном случае является пространство .
Имеются простые способы преодоления дефекта локально интерполяционных квадратурных формул – их насыщаемости. Все они основаны на простом соображении, что у соответствующей линейной комбинации двух значений составной квадратурной формулы с различными, но кратными шагами, главный член погрешности исключается. Например, для правила трапеций в силу (28) , и мы получаем повышение порядка точности, если возьмем линейную комбинацию значений формулы для числа узлов и соответственно с коэффициентами 1 и – 4.
Пусть погрешность локально-интерполяционной квадратурной формулы (25) представима в виде
где и константа не зависит от . Тогда
Отсюда получаем
и, следовательно, с точностью до имеем
Если , то
Вычисление приближенной оценки погрешности по формуле (29) называется правилом Рунге.
Число
в (30) называется уточненным (экстраполированным) по Ричардсону приближенным значением интеграла (с погрешностью ).
Замечание 6. Используя этот прием, можно уничтожить и следующие члены асимптотического разложения погрешности квадратурной формулы. Однако целесообразнее применять квадратурные формулы, сразу приводящие к ненасыщаемым алгоритмам, например, составные формулы Гаусса. Отметим, что составные квадратурные формулы, основанные на формулах Гаусса с достаточно большим числом узлов, дают хорошие результаты как для очень гладких функций, так и для функций невыской гладкости.
Замечание 7. Каждая квадратурная формула рассчитывается на определенную гладкость подинтегральной функции. Например, для правила Симпсона погрешность , если . Если квадратурная формула имеет алгебраический порядок точности , то при ее применении можно рассчитывать получить «малую погрешность» только в том случае, когда имеет непрерывные производные до порядка, не меньшего . В противном случае погрешность вычисления интеграла может оказаться большой. Для увеличения порядка гладкости подинтегральную функцию представляют в виде двух слагаемых
которые выбирают так, чтобы: содержала все особенности или их главную часть и вычислялся точно; должна иметь непрерывные производные порядка, большего , для того, чтобы интеграл можно было вычислить с достаточной точностью с помощью выбранной квадратурной формулы. Приемы разложения (32) для конкретных классов подинтегральных функций изложены в [3].
1.8. Задание. Вычислить интеграл с точностью , используя правило Симпсона и составную квадратурную формулу Гаусса с пятью узлами. Оценить погрешность используемых квадратурных формул и определить число частичных отрезков разбиения, необходимое для достижения заданной точности вычисления интеграла.
Замечание 8. Обычно для вычисления интеграла с точностью используют итерационный процесс с последовательным удвоением числа частичных отрезков разбиения.
Если , то условием останова процесса является выполнение неравенства
при этом интеграл вычисляется по формуле (31). Варианты заданий
№ варианта
|
|
|
| № варианта
|
|
|
| 1
|
| 0
| 1
| 21
|
|
| 0,1
| 2
|
| 0
| 1
| 22
|
| 0,2
| 3
|
| 0
| 1
| 23
| 0,3
| 4
|
| 1
| 2
| 24
| 0,4
| 5
|
| 0
|
| 25
| 0,5
| 6
|
| 1
| 2
| 26
| 0,6
| 7
|
| 0
| 1
| 27
| 0,7
| 8
|
| 0
|
| 28
| 0,8
| 9
|
| 0
|
| 29
| 0,9
| 10
|
| 0
|
| 30
| 1,0
|
№ варианта
|
|
|
| № варианта
|
|
|
| 11
|
| 0
|
| 31
|
|
| 0,1
| 12
|
| 0
|
| 32
| 0,2
| 13
|
| 0
| 1
| 33
| 0,3
| 14
|
| 0
| 1
| 34
| 0,4
| 15
|
| 0
| 1
| 35
| 0,5
| 16
|
| 0
| 1
| 36
| 0,6
| 17
|
| 2
| 3
| 37
| 0,7
| 18
|
| 2
| 3
| 38
| 0,8
| 19
|
| 2
| 3
| 39
| 0,9
| 20
|
| 1
| 2
| 40
| 1,0
|
|
|
|