Метод Гаусса (численное интегрирование)

Численное интегрирование функции $x^{5}+6x^{2}+1$ методом гаусса-3

Метод Гаусса — метод численного интегрирования, позволяющий повысить алгебраический порядок точности методов на основе интерполяционных формул путём специального выбора узлов интегрирования без увеличения числа используемых значений подынтегральной функции. Метод Гаусса позволяет достичь максимальной для данного числа узлов интегрирования алгебраической точности.

Например, для двух узлов можно получить метод 3-го порядка точности

$I\approx {\frac {b-a}{2}}\left(f\left({\frac {a+b}{2}}-{\frac {b-a}{2{\sqrt {3}}}}\right)+f\left({\frac {a+b}{2}}+{\frac {b-a}{2{\sqrt {3}}}}\right)\right)\,$,

тогда как для равноотстоящих узлов метода выше 2-го порядка получить невозможно. В общем случае, используя $n$ точек, можно получить метод с порядком точности $2n-1$. Значения узлов метода Гаусса по $n$ точкам являются корнями полинома Лежандра степени $n$. Значения весов вычисляются по формуле $a_{i}={\frac {2}{(1-x_{i}^{2})\,[P_{n}'(x_{i})]^{2}}}$, где $P_{n}'$ — первая производная полинома Лежандра.

Для $n=3$ узлы и веса имеют следующие значения: $x_{1,3}=\pm {\sqrt {0.6}},x_{2}=0$, веса : $a_{1,3}={\frac {5}{9}},a_{2}={\frac {8}{9}}$.

(Полином определен на отрезке $[-1,1]$).

Наиболее известен метод Гаусса по пяти точкам.

Метод Гаусса — Кронрода

Недостаток метода Гаусса состоит в том, что он не имеет лёгкого (с вычислительной точки зрения) пути оценки погрешности полученного значения интеграла. Использование правила Рунге при дроблении отрезка интегрирования требует вычисления подынтегральной функции примерно в таком же числе точек, не давая при этом практически никакого выигрыша в точности, в отличие от простых методов, где точность увеличивается в несколько раз при каждом новом разбиении. Кронродом был предложен следующий метод оценки значения интеграла

$I\approx \sum _{i=1}^{n}a_{i}\,f(x_{i})+\sum _{i=1}^{n+1}b_{i}\,f(y_{i})$,

где $x_{i}$ — узлы метода Гаусса по $n$ точкам, а $3n+2$ параметров $a_{i}$, $b_{i}$, $y_{i}$ подобраны таким образом, чтобы порядок точности метода был равен $3n+1$. Тогда для оценки погрешности можно использовать эмпирическую формулу:

$\Delta =\left(200|I-I_{G}|\right)^{1.5}$,

где $I_{G}$ — приближённое значение интеграла, полученное методом Гаусса по $n$ точкам. Библиотеки gsl и SLATEC для вычисления определённых интегралов содержат подпрограммы, использующие метод Гаусса — Кронрода по 15, 21, 31, 41, 51 и 61 точкам.

См. также

• Метод прямоугольников
• Метод трапеций
• Список квадратурных формул

Литература

1. Болтачев Г.Ш. Численные методы в теплофизике. Курс лекций Лекция 3: Численное интегрирование