МатВектор

Command Palette

Search for a command to run...

🔢 Численные методы

Квадратура Гаусса

англ. Gaussian quadrature

Интеграл по n правильно выбранным точкам точен на многочленах степени 2n − 1: узлы — корни полиномов Лежандра, веса положительны; уже два узла «умеют» кубики.

Дегустация супа не требует съедания кастрюли: две правильно взятые ложки скажут о блюде больше, чем десяток случайных. Квадратура Гаусса выбирает ложки математически: и узлы, и веса подгоняются так, чтобы интеграл от многочлена степени считался точно всего по точкам. Для сравнения: правило трапеций с теми же узлами точно лишь на прямых, а Симпсону нужны точка ради точности на кубиках. Узлы — корни полиномов Лежандра, все лежат внутри отрезка и сгущаются к краям, веса положительны — сумма ведёт себя прилично.

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

узлов Гаусса — корни полинома Лежандра ; пара узлов с весами 1 уже точна на кубиках: и , впервые ошибается на

Частые вопросы

Откуда берутся узлы Гаусса и почему они неровные?

Из подсчёта степеней свободы. У узлов и весов вместе параметров, и их хватает на точное интегрирование базиса от константы до . Система условий решается, и её ответ — корни полинома Лежандра. Равномерные узлы тратят параметры зря и страдают у концов отрезка, где функции обычно меняются резче; узлы Гаусса сгущаются именно там, выравнивая вклад. Цена — при смене числа узлов сетка строится заново.

Когда Гаусс проигрывает трапециям?

На некрасивых подынтегральных функциях. Высокий порядок обещан лишь для гладких, многочленоподобных функций; разрыв, пик или особенность у края обнуляют преимущество — нечему раскладываться в ряд. Адаптивные схемы на трапециях или Симпсоне сами сгущают сетку там, где локальная оценка ошибки велика, и честно сообщают о достигнутой точности. Плюс практический момент: сетки Гаусса разных порядков не вложены друг в друга, зато пары Гаусса–Кронрода дают сразу и результат, и оценку ошибки.

пошагово

Решатель интегралов

Получить точное значение для сравнения

Численную квадратуру не с чем сверить без точного ответа: решатель возьмёт ваш интеграл аналитически.

Открыть решатель