Методы наивысшей алгебраической точности (Гаусса - Кристоффеля)

Материал из MachineLearning.

Перейти к: навигация, поиск

Содержание

Постановка задачи

Рассмотрим задачу поиска определённого интеграла вида

( 1 )

 I=\int_a^b{f(x)\rho(x)dx},\ \rho(x)>0,

где функция f(x) непрерывна на отрезке [a,b], а весовая функция \rho(x) непрерывна на интервале (a,b). Выразить интеграл через элементарные функции в общем случае не удаётся, поэтому обычно f(x) заменяют на некоторую аппроксимирующую функцию \varphi(x)\approx f(x). Она подбирается таким образом, чтобы интеграл от неё легко считался в элементарных функциях. Стандартный пример \varphi(x) - некоторый обобщённый интерполяционный многочлен. При этом f(x) заменяется линейным выражением, со значениями в узлах в качестве коэффициентов:

( 2 )

 f(x)=\sum_{i=0}^n {f(x_i)\varphi_i(x)} + r(x),

где функция r(x) - остаточный член аппроксимации. Подстановкой (2) в (1) получаем формулу численного интегрирования (квадратурную формулу):

( 3 )

 F=\sum_{i=0}^n {c_i f(x_i)} + R,

c_i = \int_a^b{\varphi(x)\rho(x)dx},\ \ R=\int_a^b{r(x)\rho(x)dx},

где x_i называются узлами, c_i - весами, а R - погрешностью или остаточным членом. Интеграл приближённо заменяется суммой, причём узлы и коэффициенты этой суммы не зависят от f(x). Таким образом, задача сводится к отысканию подходящих наборов узлов и весов, таких, чтобы обеспечить минимизацию погрешности R в приемлемое время.

Изложение метода

Анализ метода и ошибок

Числовой пример

Рекомендации программисту

Заключение

Список литературы

  • Н.Н.Калиткин. Численные методы М.: Наука, 1978.
Личные инструменты