![]() |
Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | |
|
Постараемся повысить точность квадратурной формулы, используя неравномерное расположение узлов. Поставим задачу: построить квадратурную формулу
, (12.29)
которая при заданном n была бы точна для многочлена возможно большей степени. Такие формулы называют квадратурными формулами наивысшей алгебраической степени точности или формулами Гаусса.
Как обычно, обозначает весовую функцию, и
– непрерывная функция. Для удобства нумеруем узлы от 1 до n. Потребуем, чтобы формула (1) была точна для любого алгебраического многочлена степени m. Это эквивалентно требованию, чтобы формула была точной для функций:
Получаем условия:
(12.30)
Условия (12.30) представляют собой нелинейную систему из m+1 уравнений относительно неизвестных:
Число уравнений равно числу неизвестных, если
. Это число определяет наивысшую степень многочлена, для которого формула (12.29) является точной.
Выпишем отдельно уравнение для :
(12.31)
Уравнение (12.31) можно рассматривать как уравнение нормировки коэффициентов квадратурной формулы (12.29).
Рассмотрим примеры. Построим формулы Гаусса для разного числа узлов. Пусть имеем единичную весовую функцию (иначе говоря, весовая функция отсутствует) и стандартный, симметричный относительно начала координат, отрезок интегрирования
Пример 12. 11. n=1.
Получаем: Приходим к формуле прямоугольников:
Формула прямоугольников точна для любого многочлена первой степени.
Пример 12.12.n=2.
В данном случае m=3 и система (12.30) включает четыре уравнения:
Решаем систему. Во втором и в четвертом уравнениях перенесем второе слагаемое в правую часть и поделим четвертое уравнение на первое. Получим: , или
– узлы расположены симметрично. Учитывая симметрию, получаем из второго уравнения:
. Из первого уравнения следует, что
Наконец из третьего уравнения находим:
. Получаем квадратурную формулу
, (12.32)
точную для многочленов до третьей степени включительно. Обратим внимание, что из второго и четвертого уравнений (включающих нечетные степени x) следует только симметрия квадратурной формулы. Численные значения узлов и коэффициентов определяются из первого и третьего уравнений.
Интересен геометрический смысл полученной формулы. Через две заданные точки: – можно провести бесчисленное множество кривых третьего порядка
. Как следует из найденной квадратурной формулы, площади, заключенные между осью x и любой из этих кривых, и ограниченные вертикальными линиями
, одинаковы и равны
.
Формула (12.32) может быть получена также из геометрических соображений. Ненулевой вклад в интеграл вносят лишь четные функции:
и
. Интеграл в симметричных пределах от нечетных функций
равен нулю. Чтобы квадратурная формула правильно учитывала нулевой вклад в интеграл нечетных функций
, достаточно, чтобы формула была симметричной: узлы должны быть расположены симметрично
, а весовые коэффициенты должны быть равны
.
Опишем теперь площадь под кривой . Из соображений симметрии достаточно рассмотреть отрезок
. Разместим узел
таким образом, чтобы выполнялось равенство:
–
т.е. чтобы значение параболы в точке численно равнялось половине площади под кривой
(коэффициент
не играет роли). Отсюда
. Значение параболы в симметричном узле
даст вторую половину площади. Ясно, что те же самые отсчеты описывают также значение площади, ограниченной линией постоянной высоты
. Приходим, таким образом вновь к формуле (12.32).
![]() |
Пример 12.13. По формуле Гаусса с двумя узлами вычислим интеграл
и сравним найденное значение с точным. Точное значение интеграла легко вычисляется. После замены переменной имеем:
.
По формуле Гаусса . Поскольку формула Гаусса с двумя узлами точна для любого многочлена степени m =3, получаем точное значение.
Пример 12.14. Используя формулу Гаусса с двумя узлами, вычислим интеграл
.
Точное значение интеграла равно .
Чтобы применить формулу Гаусса, преобразуем отрезок интегрирования к стандартному отрезку . Полагаем:
.
Вновь получаем, что формула Гаусса дает точное значение интеграла:
.
Пример 12.15. Построим формулу Гаусса с тремя узлами.
Пусть по-прежнему Чтобы избежать необходимости решения системы из шести уравнений, используем геометрические соображения. Формула с тремя узлами точна для многочлена степени
. Поскольку интеграл от x в нечетной степени равен нулю, достаточно построить квадратурную формулу, точную для функций
. Эти функции показаны на рис. 12.8. Из соображений симметрии заключаем, что два узла должны быть расположены симметрично, и один узел должен находиться в центре:
Для описания функций
достаточно рассмотреть половину отрезка интегрирования
. Отсчет в точке
должен полностью описывать площадь под кривыми
. Получаем уравнения для нахождения
:
Поделив второе уравнение на первое, находим , или
. С учетом найденного значения из первого уравнения получаем
. Отсчет в точке
необходим для описания площади, ограниченной линией
. Из условия нормировки
находим .
Пример 12.16. Вычислим интеграл , используя формулу Гаусса с двумя и тремя узлами. Точное значение интеграла легко вычисляется:
. Чтобы привести отрезок интегрирования к стандартному отрезку [-1,1], заменим переменную интегрирования:
По формуле Гаусса с двумя узлами () получаем:
, и погрешность равна
.
Аналогично вычисляем по формуле Гаусса с тремя узлами:
, и погрешность равна 0.0014.
Сравним полученные результаты с результатами расчетов по формуле Симпсона с тремя узлами:
, погрешность равна 0.094.
Видим, что погрешность формулы Гаусса с тем же количеством узлов примерно в 70 раз меньше.
Дата публикования: 2015-04-07; Прочитано: 2369 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!