Квадратура Гаусса
англ. Gaussian quadrature
Интеграл по n правильно выбранным точкам точен на многочленах степени 2n − 1: узлы — корни полиномов Лежандра, веса положительны; уже два узла «умеют» кубики.
Дегустация супа не требует съедания кастрюли: две правильно взятые ложки скажут о блюде больше, чем десяток случайных. Квадратура Гаусса выбирает ложки математически: и узлы, и веса подгоняются так, чтобы интеграл от многочлена степени считался точно всего по точкам. Для сравнения: правило трапеций с теми же узлами точно лишь на прямых, а Симпсону нужны точка ради точности на кубиках. Узлы — корни полиномов Лежандра, все лежат внутри отрезка и сгущаются к краям, веса положительны — сумма ведёт себя прилично.
Минимальный случай, : узлы , веса по единице, и . Проверка на : точный интеграл , приближение — точно. На — ноль и ноль. А вот : точно , Гаусс даёт — ошибка появляется ровно на четвёртой степени, как и обещано порогом . Переход к произвольному отрезку — линейной заменой переменной. Полная техника — в уроке численное интегрирование, поведение на некрасивых функциях — в источниках погрешностей, навык на задачах — в тренажёре численных методов.
Частые вопросы
Откуда берутся узлы Гаусса и почему они неровные?
Из подсчёта степеней свободы. У узлов и весов вместе параметров, и их хватает на точное интегрирование базиса от константы до . Система условий решается, и её ответ — корни полинома Лежандра. Равномерные узлы тратят параметры зря и страдают у концов отрезка, где функции обычно меняются резче; узлы Гаусса сгущаются именно там, выравнивая вклад. Цена — при смене числа узлов сетка строится заново.
Когда Гаусс проигрывает трапециям?
На некрасивых подынтегральных функциях. Высокий порядок обещан лишь для гладких, многочленоподобных функций; разрыв, пик или особенность у края обнуляют преимущество — нечему раскладываться в ряд. Адаптивные схемы на трапециях или Симпсоне сами сгущают сетку там, где локальная оценка ошибки велика, и честно сообщают о достигнутой точности. Плюс практический момент: сетки Гаусса разных порядков не вложены друг в друга, зато пары Гаусса–Кронрода дают сразу и результат, и оценку ошибки.
Решатель интегралов
Получить точное значение для сравнения
Численную квадратуру не с чем сверить без точного ответа: решатель возьмёт ваш интеграл аналитически.
Открыть решатель