Каково современное состояние в сильно колебательных интегральных вычислениях?


23

Каково современное состояние приближения сильно колебательных интегралов как в одном измерении, так и в более высоких измерениях с произвольной точностью?


Это плохо ... пока нет общего метода ... Просто много попыток, но ожидаем, что они терпят неудачу время от времени ... Некоторые статьи утверждают, что у них есть джекпот, но когда это звучит слишком хорошо, чтобы быть правдой ... это так.

@ Гиги: Добро пожаловать в SciComp! Ваш комментарий немного расплывчатый; Не могли бы вы пояснить, почему вы считаете, что современное состояние в приближении сильно колебательных интегралов является плохим?
Джефф Оксберри

Что ж, это действительно правда, что в вычислении сильно колебательных интегралов еще нет «волшебной пули», но мы обходимся тем, что имеем, и всегда благодарны, если они работают.
JM

Ответы:


19

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

Существует ряд эффективных методов для квадратуры колебательных интегралов. Есть методы, подходящие для конечных колебательных интегралов, и есть методы для бесконечных колебательных интегралов.

Для бесконечных колебательных интегралов два наиболее эффективных метода - это метод Лонгмана и модифицированная двойная экспоненциальная квадратура, обусловленная Урой и Мори. (Но см. Также эти две статьи Арье Изерл.)

Метод Лонгмана основан на преобразовании колебательного интеграла в чередующийся ряд путем разделения интервала интегрирования, а затем суммирования чередующегося ряда методом преобразования последовательности. Например, при интегрировании колебательного интеграла вида

0f(t)sintdt

каждый превращает это в переменную сумму

k=0kπ(k+1)πf(t)sintdt

Члены этой чередующейся суммы вычисляются с помощью некоторого квадратурного метода, такого как схема Ромберга или квадратура Гаусса. В оригинальном методе Лонгмана использовалось преобразование Эйлера , но современные реализации заменяют Эйлера более мощными методами ускорения сходимости, такими как преобразование Шенкса или преобразование Левина .

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

Для конечных осциллирующих интегралов Писсенс (один из авторов QUADPACK) и Брандерс в двух работах подробно описывают модификацию квадратуры Кленшоу-Кертиса (то есть построение полиномиального чебышевского разложения неосциллирующей части подынтегрального выражения). Метод Левина , с другой стороны, использует метод коллокации для квадратуры. (Мне сказали, что теперь есть более практичная версия старого режима ожидания, метод Филона, но у меня нет опыта работы с ним.)


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


11

Помимо «многомерного против одномерного» и «конечного диапазона против бесконечного диапазона», важной категоризацией для методов является «один конкретный тип осциллятора (обычно тип Фурье: , и т. Д.). или типа Бесселя: и т. д.) против более общего осциллятора ( или даже более общего осциллятора ) ".exp ( i t ) J 0 ( t ) exp ( i g ( t ) ) w ( t )sin(t)exp(it)J0(t)exp(ig(t))w(t)

Во-первых, колебательные методы интеграции ориентированы на конкретные генераторы. Как сказал Дж.М. , наиболее известными являются метод Филона и метод Кленшоу-Кертиса (эти два тесно связаны) для интегралов с конечным диапазоном, а также методы, основанные на экстраполяции рядов, и метод двойной экспоненты Оуры и Мори для интегралов с бесконечным диапазоном.

Совсем недавно были найдены некоторые общие методы. Два примера:

  1. Метод Левина, основанный на коллокации, для любого ( Левин 1982 ) или более поздней версии для любого осциллятора определяемого линейным ОДУ ( Левин 1996 как связанный с JM ). Mathematica использует метод Левина для интегралов, не охватываемых более специализированными правилами.w ( t )exp(ig(t))w(t)

  2. Метод Huybrechs и Vandewalle, основанный на аналитическом продолжении по сложному пути, где подынтегральное выражение не является осциллирующим ( Huybrechs and Vandewalle 2006 ).

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

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

Создание автоматических подпрограмм для более общих методов может быть затруднено в большинстве языков программирования (C, Python, Fortran и т. Д.), В которых вы обычно ожидаете программировать свой интеграл как функцию / подпрограмму и передавать его в подпрограмму интегратора, потому что чем больше общие методы должны знать структуру подынтегрального выражения (какие части выглядят осциллирующими, какой тип осциллятора и т. д.) и не могут рассматривать его как «черный ящик».


Бумага Хайбрехса / Вандевалле - это то, чего я еще не видел, так что +1 за это. Похоже, что это похоже на исследование, проведенное Temme и другими для оценки специальных функций, за исключением того, что асимптотические разложения не участвуют в Huybrechs / Vandewalle. Кроме того, я думаю, что несколько решателей применили аналогичный подход к первой задаче из трехзначной задачи Трефетена.
JM

2

Вы также можете проверить работу Marnix Van Daele и соавторов. Смотрите, например, это и это .

Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.