Преобразование Ганкеля - Hankel transform

В математика, то Преобразование Ганкеля выражает любую заданную функцию ж(р) как взвешенную сумму бесконечного числа Функции Бесселя первого рода Jν(кр). Все функции Бесселя в сумме имеют один и тот же порядок ν, но отличаются масштабным коэффициентом k вдоль р ось. Необходимый коэффициент Fν каждой функции Бесселя в сумме в зависимости от масштабного коэффициента k составляет преобразованную функцию. Преобразование Ханкеля - это интегральное преобразование и был впервые разработан математиком Герман Ганкель. Он также известен как преобразование Фурье – Бесселя. Так же, как преобразование Фурье для бесконечного интервала связана с Ряд Фурье на конечном интервале, поэтому преобразование Ганкеля на бесконечном интервале связано с Ряд Фурье – Бесселя на конечном интервале.

Определение

В Преобразование Ганкеля порядка функции ж(р) дан кем-то

где это Функция Бесселя первого вида заказа с . Обратное преобразование Ганкеля Fν(k) определяется как

что можно легко проверить, используя соотношение ортогональности, описанное ниже.

Область определения

Обращение преобразования Ганкеля функции ж(р) действует в каждой точке, в которой ж(р) непрерывна при условии, что функция определена в (0, ∞), кусочно-непрерывна и имеет ограниченную вариацию на каждом конечном подынтервале в (0, ∞), и

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

Альтернативное определение

Альтернативное определение гласит, что преобразование Ганкеля грамм(р) является[1]

Эти два определения связаны:

Если , тогда

Это означает, что, как и в предыдущем определении, преобразование Ханкеля, определенное таким образом, также является своим собственным обратным:

Очевидная область теперь имеет условие

но это можно продлить. Согласно приведенной выше ссылке, мы можем взять интеграл как предел, поскольку верхний предел стремится к бесконечности ( несобственный интеграл а не Интеграл Лебега ), и таким образом преобразование Ганкеля и его обратная работа для всех функций из L2 (0, ∞).

Преобразование уравнения Лапласа

Преобразование Ханкеля можно использовать для преобразования и решения Уравнение Лапласа выражается в цилиндрических координатах. При преобразовании Ганкеля оператор Бесселя превращается в умножение на .[2] В осесимметричном случае уравнение в частных производных преобразуется как

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

Ортогональность

Функции Бесселя образуют ортогональный базис относительно весового коэффициента р:[3]

Теорема Планшереля и теорема Парсеваля

Если ж(р) и грамм(р) таковы, что их преобразования Ганкеля Fν(k) и граммν(k) хорошо определены, то Теорема Планшереля состояния

Теорема Парсеваля, в котором говорится

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

Связь с многомерным преобразованием Фурье

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

Рассмотрим функцию из -мерный вектор р. Его -мерное преобразование Фурье определяется как

Чтобы переписать его в гиперсферических координатах, мы можем использовать разложение плоской волны на -мерные гиперсферические гармоники :[4]
где и - множества всех гиперсферических углов в -пространство и -Космос. Это дает следующее выражение для -мерное преобразование Фурье в гиперсферических координатах:
Если мы расширим и в гиперсферических гармониках:
преобразование Фурье в гиперсферических координатах упрощается до
Это означает, что функции с угловой зависимостью в виде гиперсферической гармоники сохраняют ее при многомерном преобразовании Фурье, в то время как радиальная часть подвергается преобразованию Ганкеля (с точностью до некоторых дополнительных факторов, таких как ).

Особые случаи

Преобразование Фурье в двух измерениях

Если двумерная функция ж(р) расширяется в многополюсный ряд,

то его двумерное преобразование Фурье имеет вид

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

Преобразование Фурье в трех измерениях

Если трехмерная функция ж(р) расширяется в многополюсный ряд над сферические гармоники,

то его трехмерное преобразование Фурье дается выражением

где
преобразование Ганкеля порядка .

Этот вид преобразования Ганкеля полуцелого порядка также известен как сферическое преобразование Бесселя.

Преобразование Фурье в d размеры (радиально-симметричный случай)

Если d-мерная функция ж(р) не зависит от угловых координат, то его d-мерное преобразование Фурье F(k) также не зависит от угловых координат и определяется выражением[5]

которое является преобразованием Ганкеля порядка до фактора .

2D-функции в ограниченном радиусе

Если двумерная функция ж(р) расширяется в многополюсный ряд а коэффициенты разложения жм достаточно гладкие вблизи начала координат и равны нулю вне радиуса р, радиальная часть ж(р)/рм может быть расширен до степенного ряда 1- (г / р) ^ 2:

такое, что двумерное преобразование Фурье ж(р) становится

где последнее равенство следует из §6.567.1.[6] Коэффициенты разложения жм, т доступны с дискретное преобразование Фурье техники:[7] если радиальное расстояние масштабируется с

коэффициенты ряда Фурье-Чебышева грамм появиться как

Использование повторного расширения

дает жм, т выражается в виде суммы граммм, дж.

Это одна из разновидностей техник быстрого преобразования Ганкеля.

Связь с преобразованиями Фурье и Абеля

Преобразование Ханкеля является одним из членов Цикл FHA интегральных операторов. В двух измерениях, если мы определим А как Преобразование Абеля оператор F как преобразование Фурье оператор и ЧАС как оператор преобразования Ганкеля нулевого порядка, то частный случай теорема о проекции для кругово-симметричных функций утверждает, что

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

Числовая оценка

Простой и эффективный подход к численной оценке преобразования Ханкеля основан на наблюдении, что оно может быть представлено в виде свертка логарифмической заменой переменных[8]

В этих новых переменных преобразование Ханкеля имеет вид
где
Теперь интеграл можно вычислить численно с помощью сложность с помощью быстрое преобразование Фурье. Алгоритм может быть дополнительно упрощен, если использовать известное аналитическое выражение для преобразования Фурье :[9]
Оптимальный выбор параметров зависит от свойств , в частности его асимптотика при и .

Этот алгоритм известен как «квази-быстрое преобразование Ганкеля» или просто «быстрое преобразование Ганкеля».

Поскольку он основан на быстрое преобразование Фурье в логарифмических переменных, должен быть определен на логарифмической сетке. Для функций, определенных на однородной сетке, существует ряд других алгоритмов, в том числе простой квадратура, методы, основанные на теорема о проекции, и методы, использующие асимптотическое разложение функций Бесселя.[10]

Некоторые пары преобразований Ханкеля

[11]

Можно выразить в терминах эллиптические интегралы.[12]

Kп(z) это модифицированная функция Бесселя второго рода.K(z) это полный эллиптический интеграл первого рода.

Выражение

совпадает с выражением для Оператор Лапласа в полярные координаты (k, θ) применительно к сферически симметричной функции F0(k).

Преобразование Ганкеля Многочлены Цернике по сути являются функциями Бесселя (Noll 1976):

даже для пм ≥ 0.

Смотрите также

Рекомендации

  1. ^ Луи де Бранж (1968). Гильбертовы пространства целых функций. Лондон: Прентис-Холл. п.189. ISBN  978-0133889000.
  2. ^ Пуларикас, Александр Д. (1996). Справочник преобразований и приложений. Бока-Ратон, Флорида: CRC Press. ISBN  0-8493-8342-0. OCLC  32237017.
  3. ^ Понсе де Леон, Дж. (2015). «Возвращаясь к ортогональности функций Бесселя первого рода на бесконечном интервале». Европейский журнал физики. 36 (1): 015016. Bibcode:2015EJPh ... 36a5016P. Дои:10.1088/0143-0807/36/1/015016.
  4. ^ Эйвери, Джеймс Эмиль, автор. Гиперсферические гармоники и их физические приложения. ISBN  978-981-322-930-3. OCLC  1013827621.CS1 maint: несколько имен: список авторов (связь)
  5. ^ Фарис, Уильям Г. (2008-12-06). «Радиальные функции и преобразование Фурье: Notes for Math 583A, Fall 2008» (PDF). Университет Аризоны, факультет математики. Получено 2015-04-25.
  6. ^ Градштейн, И. С .; Рыжик И. М. (2015). Цвиллинджер, Даниэль (ред.). Таблица интегралов, серий и продуктов (Восьмое изд.). Академическая пресса. п. 687. ISBN  978-0-12-384933-5.
  7. ^ Секада, Хосе Д. (1999). «Численная оценка преобразования Ханкеля». Комп. Phys. Comm. 116 (2–3): 278–294. Bibcode:1999CoPhC.116..278S. Дои:10.1016 / S0010-4655 (98) 00108-8.
  8. ^ Зигман, А. Э. (1977-07-01). «Квази быстрое преобразование Ганкеля». Письма об оптике. 1 (1): 13. Bibcode:1977 ОптЛ .... 1 ... 13S. Дои:10.1364 / ol.1.000013. ISSN  0146-9592. PMID  19680315.
  9. ^ Талман, Джеймс Д. (октябрь 1978 г.). «Числовые преобразования Фурье и Бесселя в логарифмических переменных». Журнал вычислительной физики. 29 (1): 35–48. Bibcode:1978JCoPh..29 ... 35T. Дои:10.1016/0021-9991(78)90107-9. ISSN  0021-9991.
  10. ^ Кри, M.J .; Кости, П.Дж. (июль 1993 г.). «Алгоритмы для численной оценки преобразования Ханкеля». Компьютеры и математика с приложениями. 26 (1): 1–12. Дои:10.1016/0898-1221(93)90081-6. ISSN  0898-1221.
  11. ^ Папулис, Афанасиос (1981). Системы и преобразования с приложениями в оптике. Флорида, США: Издательство Кригер. С. 140–175. ISBN  978-0898743586.
  12. ^ Kausel, E .; Ирфан Байг, М. М. (2012). «Преобразование Лапласа произведений функций Бесселя: посещение более ранних формул» (PDF). Квартал прикладной математики. 70: 77–97. Дои:10.1090 / s0033-569x-2011-01239-2. HDL:1721.1/78923.