Адаптурная квадратура

Адаптивная квадратура автоматически подразделяет интервал интегрирования там, где подынтегральная функция сложна, используя оценки локальной ошибки для достижения требуемой точности с минимально возможным количеством вычислений функции.

Definition

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

Scope

Эта тема охватывает оценку локальной ошибки путем сравнения правил разных порядков или уровней уточнения, рекурсивное деление интервала пополам (адаптивная квадратура Симпсона и адаптивная квадратура Гаусса-Кронрода), глобальные бюджеты ошибок и критерии остановки, обработку сингулярностей и резких особенностей, а также разработку производственных автоматических интеграторов, таких как те, что представлены в библиотеке QUADPACK.

Core questions

• Как вычисляется локальная ошибка квадратурной оценки без знания точного интеграла?
• Как рекурсивное подразделение концентрирует усилия там, где подынтегральная функция изменяется больше всего?
• Какие критерии остановки надежно достигают требуемого допуска, избегая при этом напрасной работы?
• Как надежно обнаруживаются и обрабатываются интегрируемые сингулярности и разрывы?

Key theories

Оценка локальной ошибки и подразделение

Сравнение грубой и более точной (или более высокого порядка) оценки на подынтервале дает оценку локальной ошибки; если она превышает долю допуска, выделенную этому подынтервалу, подынтервал делится, и процедура повторяется, в противном случае его вклад принимается.

Глобально адаптивные стратегии

Вместо независимой обработки подынтервалов, глобально адаптивные интеграторы поддерживают очередь подынтервалов, упорядоченных по оценочной ошибке, и всегда уточняют наихудший, что эффективно обрабатывает локализованные сингулярности и лежит в основе процедур QUADPACK.

Mechanisms

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

Clinical relevance

Адаптивная квадратура — это то, на что полагаются универсальные процедуры интегрирования в научном программном обеспечении для предоставления результата с заданной пользователем точностью без необходимости анализа подынтегральной функции; она необходима для подынтегральных функций с пиками, поведением пограничного слоя или интегрируемыми сингулярностями, которые не поддаются фиксированному правилу, и лежит в основе автоматических интеграторов в широко используемых численных и статистических пакетах.

History

Автоматическое интегрирование с контролем ошибок развивалось в 1970-х и начале 1980-х годов, кульминацией чего стал пакет QUADPACK (1983), чьи адаптивные процедуры Гаусса-Кронрода с экстраполяцией стали де-факто стандартом и впоследствии были приняты, портированы или перереализованы во многих численных и статистических программных системах.

Key figures

• Robert Piessens
• Philip J. Davis
• Philip Rabinowitz

Seminal works

• davis1984
• piessens1983

Frequently asked questions

Как адаптивный интегратор узнает ошибку, если он не знает ответа?

Он оценивает локальную ошибку, сравнивая два приближения различной точности на одном и том же подынтервале — например, правило более высокого и более низкого порядка. Их разность аппроксимирует ошибку и направляет, где следует уточнять, даже если истинный интеграл неизвестен.

Когда адаптивная квадратура сталкивается с трудностями?

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