Методы наивысшей алгебраической точности (Гаусса - Кристоффеля)
Материал из MachineLearning.
м |
(категория, шаблон) |
||
Строка 21: | Строка 21: | ||
была точна для любого многочлена степени <tex>2n-1</tex>. Покажем, как находятся узлы и веса этих формул. | была точна для любого многочлена степени <tex>2n-1</tex>. Покажем, как находятся узлы и веса этих формул. | ||
- | Будем считать, что вес положителен <tex>\rho(x)>0</tex>, непрерывен на <tex>(a,b)</tex>, может обращаться в нуль или бесконечность на концах отрезка так, чтобы существовал <tex>\int_a^b\rho(x)dx</tex>. Известно{{ | + | Будем считать, что вес положителен <tex>\rho(x)>0</tex>, непрерывен на <tex>(a,b)</tex>, может обращаться в нуль или бесконечность на концах отрезка так, чтобы существовал <tex>\int_a^b\rho(x)dx</tex>. Известно{{источник?}}, что при выполнении этих условий существует полная система алгебраических многочленов <tex>P_m(x)</tex>, ортогональных на <tex>[a,b]</tex> с заданным весом: |
{{eqno|4}} | {{eqno|4}} | ||
<p align="center"> <tex>\int_a^b P_k(x)P_m(x)\rho(x)dx=\delta_k_m ||P_k(x)||^2_L_2</tex> </p> | <p align="center"> <tex>\int_a^b P_k(x)P_m(x)\rho(x)dx=\delta_k_m ||P_k(x)||^2_L_2</tex> </p> | ||
Строка 65: | Строка 65: | ||
{{stub}} | {{stub}} | ||
[[Категория:Численное интегрирование]] | [[Категория:Численное интегрирование]] | ||
+ | [[Категория:Учебные задачи]] |
Версия 20:14, 19 октября 2008
Содержание |
Постановка задачи
Рассмотрим задачу поиска определённого интеграла вида
где функция непрерывна на отрезке
, а весовая функция
непрерывна на интервале
.
Выразить интеграл через элементарные функции в общем случае не удаётся, поэтому обычно
заменяют на некоторую аппроксимирующую функцию
. Она подбирается таким образом, чтобы интеграл от неё легко считался в элементарных функциях. Стандартный пример
- некоторый обобщённый интерполяционный многочлен. При этом
заменяется линейным выражением, со значениями в узлах в качестве коэффициентов:
где функция - остаточный член аппроксимации. Подстановкой (2) в (1) получаем формулу численного интегрирования (квадратурную формулу):
,
где называются узлами,
- весами, а
- погрешностью или остаточным членом. Интеграл приближённо заменяется суммой, причём узлы и коэффициенты этой суммы не зависят от
.
Таким образом, задача сводится к отысканию подходящих наборов узлов и весов, таких, чтобы обеспечить минимизацию погрешности
в приемлемое время.
Изложение метода
Параметрами формулы (3) являются узлы и веса. В известных формулах численного интегрирования, таких как формулы трапеций, Симпсона, принято фиксировать положение узлов и по ним находить веса. Таким образом в них не полностью используются возможности общей формулы. Логично предположить, что выбор оптимального положения узлов приведёт к улучшению работы метода.
Итак, формула (3) с n узлами содержит параметров, столько же коэффициентов у многочлена степени
. Значит, параметры можно подобрать так, чтобы квадратурная формула (3)
была точна для любого многочлена степени . Покажем, как находятся узлы и веса этих формул.
Будем считать, что вес положителен , непрерывен на
, может обращаться в нуль или бесконечность на концах отрезка так, чтобы существовал
. Известно[источник?], что при выполнении этих условий существует полная система алгебраических многочленов
, ортогональных на
с заданным весом:
Все нули этих многочленов действительны и расположены на интервале .
Составим по узлам интегрирования многочлен n-й степени . Функция
при
является многочленом степени не выше 2n-1. Следовательно, для неё формула Гаусса-Кристоффеля точна. Тогда получим:
так как . Значит многочлен
ортогонален всем многочленам
степени
.
Если разложить в ряд по рассматриваемым ортогональным многочленам и подставить этот ряд в условие ортогональности (5), то получим:
,
,
т.е. все коэффициенты разложения при
. Это значит, что
с точностью до численного множителя совпадает с
. Значит, узлами формулы Гаусса-Кристоффеля являются нули многочленов соответствующей степени
, ортогональных на
с весом
.
Веса интегрирования нетрудно определить, если узлы уже найдены. Функция
есть многочлен степени n-1, т.е. для неё формула Гаусса-Кристоффеля точна. Подставляя её в формулу (3) и учитывая, что эта формула равна нулю во всех узлах, кроме m-го, получим веса формулы Гаусса-Кристоффеля:
из которого видно, что все веса положительны. Подставляя в формулу Гаусса , получим соотношение:
,
из которого следует равномерная ограниченность весов.
Формулы Гаусса-Кристоффеля называют также формулами наивысшей алгебраической точности, поскольку для произвольного многочлена степени выше формула (3) с
узлами уже не может быть точной.
Различные случаи
Рассмотрим некоторые частные случаи:
- Собственно формула Гаусса соответствует
. Линейным преобразованием аргумента можно перейти к отрезку
. На нём ортогональны с единичным весом многочлены Лежандра. Если обозначить их узлы и соответствующие веса через
, то обратным линейным преобразованием можно получить узлы и веса для произвольного отрезка
,
.
получаем формулу средних. Погрешность формулы Гаусса (приводится без вывода) пропорциональна той производной, которая соответствует низшей неучтённой формуле аргумента; верхняя граница погрешности равна
.
.
- Формула Эрмита позволяет интегрировать на отрезке
с весом
. При этих условиях ортогональны многочлены Чебышева первого рода
. Соответствующие узлы и веса интегрирования равны:
Отметим, что веса во всех узлах одинаковы. На произвольный отрезок эти узлы и веса преобразуются так же, как в формуле Гаусса. Погрешность формулы Эрмита не превышает
.
Анализ метода и ошибок
Числовой пример
Рекомендации программисту
Заключение
Список литературы
- Н.Н.Калиткин. Численные методы М.: Наука, 1978.