Функция Бесселя

редактировать

Семейства решений связанных дифференциальных уравнений Функции Бесселя - это радиальная часть форм колебаний кругового барабана.

Функции Бесселя, сначала выступ математиком Даниэлем Бернулли, объединенные Фридрихом Бесселем, являются каноническими решениями y (x) дифференциального уравнения Бесселя

Икс 2 d 2 ydx 2 + xdydx + (x 2 - α 2) y = 0 {\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + x {\ frac {dy} {dx}} + \ left (x ^ {2} - \ alpha ^ {2} \ right) y = 0}{\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + x {\ frac {dy} {dx}} + \ left (x ^ {2} - \ alpha ^ {2} \ right) y = 0}

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

Наиболее важные случаи, когда α является целым числом или полуцелым числом. Функции Бесселя для целого числа α также известны как цилиндрические функции или цилиндрические гармоники, они появляются в решении уравнения Лапласа в цилиндрические координаты. Сферические функции Бесселя с полуцелым числом α получаются, когда уравнение Гельмгольца решается в сферических координатах.

Содержание
  • 1 Функции Бесселя
  • 2 Определения
    • 2.1 Функции Бесселя первого рода: J α
      • 2.1.1 Интегралы Бесселя
      • 2.1.2 Связь с гипергеометрическими рядами
      • 2.1.3 Связь с полиномами Лагерра
    • 2.2 Функции Бесселя второго рода: Y α
    • 2.3 Функции Ганкеля: H α, H α
    • 2.4 Модифицированные функции Бесселя: I α, K α
    • 2.5 Сферические функции Бесселя: j n, y n
      • 2.5.1 Производящая функция
      • 2.5.2 Дифференциальные отношения
    • 2.6 Сферические функции Ганкеля: h n, h n
    • 2.7 Функции Риккати - Бесселя: S n, C n, ξ n, ζ n
  • 3 Asy Мптотические формы
  • 4 Полнообластные приближения с элементарными функциями
  • 5 Свойства
    • 5.1 Рекуррентные соотношения
  • 6 Теорема умножения
  • 7 Нули функции Бесселя
    • 7.1 Гипотеза Бурже
    • 7.2 Численные подходы
  • 8 См. Также
  • 9 Примечания
  • 10 Ссылки
  • 11 Внешние ссылки
Приложения функций Бесселя

Уравнение Бесселя возникает при нахождении разделимых решений уравнения уравнения Лапласа и уравнение Гельмгольца в цилиндрических или сферических координатах. Поэтому функции Бесселя особенно важны для многих задач распространение волн и статических потенциалов. При решении задач в цилиндрических системах получаются функции Бесселя целого порядка (α = n); в сферических задачах получаются полуцелые порядки (α = n + 1/2). Например:

Также функции Бесселя появляются в других задачах, таких как обработка сигналов (например, см. FM -синтез, окно Кайзера или фильтр Бесселя ).

Определения

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

ТипПервый видВторой вид
Функции Бесселя
Модифицированные функции Бесселя
Функции ХанкеляH. α= J α + iY αH. α= J α - iY α
Сферические функции Бесселяjnyn
Сферические функции Ганкеляh. n= j n + iy nh. n= j n - iy n

Функции Бесселя второго рода и сферические функции Бесселя второго рода иногда обозначаются N n и n n соответственно, а не Y n и y n.

Функции Бесселя первого рода: J α

Функции Бесселя первого рода, обозначенные как J α (x), являются решениями дифференциального уравнения Бесселя, которые являются конечными в начале координат (x = 0) для целого или положительного α и расходятся, когда x приближается к нулю для отрицательного нецелого числа α. Функцию можно определить по ее разложению в ряду вокруг x = 0, которое можно найти, применив метод Фробениуса к уравнению Бесселя:

J α (x) = ∑ m = 0 ∞ (- 1) мм! Γ (м + α + 1) (Икс 2) 2 м + α, {\ Displaystyle J _ {\ альфа} (х) = \ сумма _ {м = 0} ^ {\ infty} {\ frac {(-1) ^ {м}} {м! \ Gamma (m + \ alpha +1)}} {\ left ({\ frac {x} {2}} \ right)} ^ {2m + \ alpha},}{ \ displaystyle J _ {\ alpha} (x) = \ sum _ {m = 0} ^ {\ infty} {\ frac {(-1) ^ {m}} {m! \ Gamma (m + \ alpha +1)} } {\ left ({\ frac {x} {2}} \ right)} ^ {2m + \ alpha},}

где Γ (z) - это гамма -функция , сдвинутое обобщение функций факториала на нецелочисленные значения. Функция Бесселя первого рода является целой функцией, если α является целым числом, в результате этого многозначная функция с сингулярностью в нуле. Графики функций Бесселя выглядят примерно как осциллирующие синусоидальные или косинусные функции, которые убивают пропорционально x - 1 2 {\ displaystyle x ^ {- {\ frac {1} {2}}}}{\ displaystyle x ^ {- {\ frac {1} {2}}}} (см. также их асимптотические формы ниже), хотя их корни в общем случае не периодичны, за исключением асимптотики для больших x. (Ряд указывает, что −J 1 (x) является производной от J 0 (x), так же как −sin x является производной от cos x; в более общем смысле, производная из J n (x) может быть выражено через J n ± 1 (x) тождествами ниже.)

График функции Бесселя первого рода, J α (x), для целых порядков α = 0, 1, 2

Для нецелых α функций J α (x) и J −α (x) линейно независимы и, следовательно, являются решения двух дифференциальных уравнений. С другой стороны, для целого порядка n справедливо следующее соотношение (гамма-функция имеет простые полюсы при каждом из неположительных целых чисел):

J - n (x) = (- 1) n J n (Икс). {\ displaystyle J _ {- n} (x) = (- 1) ^ {n} J_ {n} (x).}{\ displaystyle J _ {- n} (x) = (- 1) ^ {n} J_ {n} (x).}

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

Интегралы Бесселя

Другое определение функции Бесселя для целых значений, возможно с использованием интегрального представления:

J n (x) = 1 π ∫ 0 π cos ⁡ (n τ - x sin ⁡ τ) d τ. {\ displaystyle J_ {n} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} \ cos (n \ tau -x \ sin \ tau) \, d \ tau.}{\ displaystyle J_ {n} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ { \ пи} \ соз (п \ тау -x \ грех \ тау) \, д \ тау.}

Другое интегральное представление:

J n (x) = 1 2 π ∫ - π π ei (x sin ⁡ τ - n τ) d τ. {\ displaystyle J_ {n} (x) = {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} e ^ {i (x \ sin \ tau -n \ tau)} \, d \ tau.}{\ displaystyle J_ {n} ( x) = {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} e ^ {i (x \ sin \ tau -n \ tau)} \, d \ tau. }

Это был подход, который использовал Бессель, и из этого определения он вывел несколько функций функций. Определение может быть расширено на нецелые порядки с помощью одного из интегралов Шлефли, для Re (x)>0:

J α (x) = 1 π ∫ 0 π cos ⁡ (α τ - x sin ⁡ τ) d τ - sin ⁡ α π π ∫ 0 ∞ e - x sinh ⁡ t - α tdt. {\ displaystyle J _ {\ alpha} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} \ cos (\ alpha \ tau -x \ sin \ tau) \, d \ tau - {\ frac {\ sin \ alpha \ pi} {\ pi}} \ int _ {0} ^ {\ infty} e ^ {- x \ sinh t- \ alpha t} \, dt.}{\ displaystyle J _ {\ alpha} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} \ cos (\ alpha \ tau -x \ sin \ tau) \, d \ tau - {\ frac {\ sin \ alpha \ pi} {\ pi}} \ int _ {0} ^ {\ infty} e ^ {- x \ sinh t - \ alpha t} \, dt.}

Связь с гипергеометрическим рядом

Функции Бесселя можно выразить через обобщенный гипергеометрический ряд as

J α (x) = (x 2) α Γ (α + 1) 0 F 1 (α + 1; - х 2 4). {\ displaystyle J _ {\ alpha} (x) = {\ frac {\ left ({\ frac {x} {2}} \ right) ^ {\ alpha}} {\ Gamma (\ alpha +1)}} \; _ {0} F_ {1} \ left (\ alpha +1; - {\ frac {x ^ {2}} {4}} \ right).}{\ displaystyle J _ {\ alpha} (x) = {\ frac {\ left ({\ frac {x} {2}} \ right) ^ {\ alpha}} {\ Гамма (\ alpha +1)}} \; _ {0} F_ {1} \ left (\ alpha +1; - {\ frac {x ^ {2}} {4}} \ right).}

Это связано с разработкой функций Бесселя в терминах функции Бесселя - Клиффорда.

Связь с полиномами Лагерра

В терминах полиномов Лагерра Lkи произвольно выбранные программы t функция Бесселя может быть выражена как

J α (x) (x 2) α знак равно e - t Γ (α + 1) ∑ k = 0 ∞ L k (α) (x 2 4 t) (k + α k) tkk!. {\ displaystyle {\ frac {J _ {\ alpha} (x)} {\ left ({\ frac {x} {2}} \ right) ^ {\ alpha}}} = {\ frac {e ^ {- t}} {\ Gamma (\ alpha +1)}} \ sum _ {k = 0} ^ {\ infty} {\ frac {L_ {k} ^ {(\ alpha)} \ left ({\ frac {x ^ {2}} {4t}} \ right)} {\ binom {k + \ alpha} {k}}} {\ frac {t ^ {k}} {k!}}.}{\ displaystyle {\ frac {J _ {\ alpha} (x)} {\ left ({\ frac {x} {2}} \ right) ^ {\ alpha} }} = {\ frac {e ^ {- t}} {\ Gamma (\ alpha +1)}} \ sum _ {k = 0} ^ {\ infty} {\ frac {L_ {k} ^ {(\ альфа)} \ left ({\ frac {x ^ {2}} {4t}} \ right)} {\ binom {k + \ alpha} {k}}} {\ frac {t ^ {k}} {k! }}.}

Функции Бесселя второго kind: Y α

Функции Бесселя второго рода, обозначаемые Y α (x), иногда вместо этого обозначаемые N α (x), являются решениями дифференциала Бесселя уравнение, которое имеет особенность в начало координат (x = 0) и является многозначным. Иногда их называют функции Вебера, так как они были введены H. М. Вебер (1873), а также функции Неймана после Карла Неймана.

График функции Бесселя второго рода, Y α (x), для целочисленных порядков α = 0, 1, 2

Для нецелого числа α Y α (x) связано с J α (x) по

Y α (x) = J α (x) cos ⁡ (α π) - J - α (x) sin ⁡ (α π). {\ Displaystyle Y _ {\ alpha} (x) = {\ frac {J _ {\ alpha} (x) \ cos (\ alpha \ pi) -J _ {- \ alpha} (x)} {\ sin ( \ alpha \ pi)}}.}{\ displaystyle Y _ {\ alpha} (x) = {\ frac {J _ {\ alpha } (х) \ соз (\ альфа \ пи) -J _ {- \ альфа} (х)} {\ грех (\ альфа \ пи)}}.}

В случае целого порядка n функция достигается путем взятия предела, когда нецелое число α к n:

Y n (x) = lim α → n Y α (Икс). {\ displaystyle Y_ {n} (x) = \ lim _ {\ alpha \ to n} Y _ {\ alpha} (x).}Y_ {n} (x) = \ lim _ {\ alpha \ to n} Y _ {\ alpha} (x).

Если n - неотрицательное целое число, мы имеем ряд

Y n (г) знак равно - (г 2) - N π ∑ К знак равно 0 N - 1 (N - К - 1)! к! (z 2 4) k + 2 π J n (z) ln ⁡ z 2 - (z 2) n π ∑ k = 0 ∞ (ψ (k + 1) + ψ (n + k + 1)) (- z 2 4) кк! (п + к)! {\ displaystyle Y_ {n} (z) = - {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {- n}} {\ pi}} \ sum _ {k = 0 } ^ {n-1} {\ frac {(nk-1)!} {k!}} \ left ({\ frac {z ^ {2}} {4}} \ right) ^ {k} + {\ frac {2} {\ pi}} J_ {n} (z) \ ln {\ frac {z} {2}} - {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {n}} {\ pi}} \ sum _ {k = 0} ^ {\ infty} (\ psi (k + 1) + \ psi (n + k + 1)) {\ frac {\ left (- { \ frac {z ^ {2}} {4}} \ right) ^ {k}} {к! (n + k)!}}}{\ displaystyle Y_ {n} (z) = - {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {- n}} {\ pi}} \ sum _ {k = 0} ^ { n-1} {\ frac {(nk-1)!} {k!}} \ left ( {\ frac {z ^ {2}} {4}} \ right) ^ {k} + {\ frac {2} {\ pi}} J_ {n} (z) \ ln {\ frac {z} {2 }} - {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {n}} {\ pi}} \ sum _ {k = 0} ^ {\ infty} (\ psi ( k + 1) + \ psi (n + k + 1)) {\ frac {\ left (- {\ frac {z ^ {2}} {4}} \ right) ^ {k}} {k! (n + к)!}}}

где ψ (z) {\ displaystyle \ psi (z)}\ psi (z) - дигамма-функция, логарифмическая производная от гамма-.

Существует также соответствующая функция интегральной формулы (для Re (x)>0):

Y n (x) = 1 π ∫ 0 π sin ⁡ (x sin ⁡ θ - n θ) d θ - 1 π ∫ 0 ∞ (ent + (- 1) ne - nt) е - х з тдт. {\ displaystyle Y_ {n} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} \ sin (x \ sin \ theta -n \ theta) \, d \ theta - {\ frac {1} {\ pi}} \ int _ {0} ^ {\ infty} \ left (e ^ {nt} + (- 1) ^ {n} e ^ {- nt} \ right) e ^ {- x \ sh t} \, dt.}{\ displaystyle Y_ {n} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} \ sin (x \ sin \ theta -n \ theta) \, d \ theta - {\ frac {1} {\ pi}} \ int _ {0} ^ {\ infty} \ left (e ^ {nt} + (- 1) ^ {n} e ^ {- nt} \ right) e ^ {- x \ sinh t} \, dt.}

Yα(x) необходимо как второе линейно независимое решение уравнения Бесселя, когда α является целым числом. Но Y α (x) имеет большее значение, чем это. Его можно рассматривать как «естественного» партнера J α (x). См. Также подраздел о функциях Ганкеля ниже.

Когда α является целым числом, кроме того, как было аналогично для функций первого рода, справедливо следующее соотношение:

Y - n (x) = (- 1) n Y n (Икс). {\ displaystyle Y _ {- n} (x) = (- 1) ^ {n} Y_ {n} (x).}{\ displaystyle Y _ {- n} (x) = (- 1) ^ {n} Y_ {n} (x).}

И J α (x), и Y α (x) - голоморфные функции от x на комплексной плоскости, разрезанной вдоль отрицательной действительной оси. Когда α является целым числом, функции Бесселя J являются целыми функциями от x. Если x удерживается фиксированным на ненулевом значении, то Бесселя выполняет целые функции от α.

Функции Бесселя второго рода, когда одним из основных примеров решения типа в теореме Фукса.

Функции Ганкеля: H. α, H. α

Еще одна важная формулировка двух линейно независимых решения уравнения Бесселя являются функции Ганкеля первого и второго рода, H. α(x) и H. α(x), определяемые как

H α (1) (Икс) знак равно J α (Икс) + я Y α (Икс), ЧАС α (2) (Икс) = J α (Икс) - я Y α (Икс), {\ Displaystyle {\ begin {выровнено} H _ {\ alpha} ^ { (1)} (x) = J _ {\ alpha} (x) + iY _ {\ alpha} (x), \\ H _ {\ alpha} ^ {(2)} (x) = J _ {\ alpha} (x) -iY _ {\ alpha} (x), \ end {align}}}{\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (x) = J _ {\ alpha} (x) + iY _ {\ alpha} (x), \\ H _ {\ alpha} ^ {(2)} (x) = J _ {\ alpha} (x) -iY _ {\ альфа} (х), \ end {align}}}

, где i - мнимая единица. Эти линейные комбинации также известны как функции Бесселя третьего рода ; они являются двумя линейно независимыми решениями дифференциального уравнения Бесселя. Они названы в честь Германа Ганкеля.

. Эти формы линейной комбинации представлены многочисленным на первый взгляд простых свойств, такими как асимптотические формулы или интегральные представления. Здесь «простой» появление фактора формы e. Таким образом, можно думать, что функция Бесселя второго рода естественным образом возникает как мнимая часть функций Ганкеля.

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

Используя предыдущие соотношения, их можно выразить как

H α (1) (x) = J - α (x) - e - α π i J α (x) i sin ⁡ α π, H α (2) (x) = J - α (x) - e α π i J α (x) - i sin ⁡ α π. {\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (x) = {\ frac {J _ {- \ alpha} (x) -e ^ {- \ alpha \ pi i } J_ {\ alpha} (x)} {i \ sin \ alpha \ pi}}, \\ H _ {\ alpha} ^ {(2)} (x) = {\ frac {J _ {- \ alpha } (x) -e ^ {\ alpha \ pi i} J _ {\ alpha} (x)} {- i \ sin \ alpha \ pi}}. \ end {align}}}{\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (x) = {\ frac {J _ {- \ alpha} (x) -e ^ {- \ alpha \ pi i} J _ {\ alpha} (x)} {i \ sin \ alpha \ pi}}, \\ H _ {\ alpha} ^ {(2)} (x) = {\ frac {J _ {- \ alpha } (x) -e ^ {\ alpha \ pi i} J _ {\ alpha} (x)} {- i \ sin \ alpha \ pi}}. \ end {align}}}

Если α - целое число, предел имеет быть рассчитанным. Следующие соотношения действительны, независимо от того, являются ли α целым или нет:

H - α (1) (x) = e α π i H α (1) (x), H - α (2) (x) = е - α π i H α (2) (x). {\ displaystyle {\ begin {align} H _ {- \ alpha} ^ {(1)} (x) = e ^ {\ alpha \ pi i} H _ {\ alpha} ^ {(1)} (x), \\ H _ {- \ alpha} ^ {(2)} (x) = e ^ {- \ alpha \ pi i} H _ {\ alpha} ^ {(2)} (x). \ End {align}}}{\ displaystyle {\ begin {align} H _ {- \ alpha} ^ {(1)} (x) = e ^ {\ alpha \ pi i} H _ {\ alpha} ^ {(1)} (x), \\ H _ {- \ alpha} ^ {(2)} (x) = e ^ {- \ альфа \ пи я} Н _ {\ альфа} ^ {(2)} (х). \ конец {выровнено}}}

В частности, если α = m + 1/2 с ma неотрицательным целым числом, из приведенных выше указанных прямо следует, что

J - (m + 1 2) (x) = (- 1) m + 1 Y m + 1 2 (х), Y - (м + 1 2) (х) = (- 1) м Дж м + 1 2 (х). {\ displaystyle {\ begin {align} J _ {- (m + {\ frac {1} {2}})} (x) = (- 1) ^ {m + 1} Y_ {m + {\ frac {1} {2}}} (x), \\ Y _ {- (m + {\ frac {1} {2}})} (x) = (- 1) ^ {m} J_ {m + {\ frac {1}) {2}}} (х). \ End {align}}}{\ displaystyle {\ begin {align} J _ {- (m + {\ frac {1} {2}})} (x) = (- 1) ^ {m + 1} Y_ {m + {\ frac {1} {2} }}} (x), \\ Y _ {- (m + {\ frac {1} {2}})} (x) = (- 1) ^ {m} J_ {m + {\ frac {1} {2 }}} (x). \ end {align}}}

Они полезны при разработке сферических функций Бесселя (см. Ниже).

Функции Ганкеля допускают следующие интегральные представления для Re (x)>0:

H α (1) (x) = 1 π i ∫ - ∞ + ∞ + i π ex sinh ⁡ t - α tdt, ЧАС α (2) (Икс) знак равно - 1 π я ∫ - ∞ + ∞ - я π ex sinh ⁡ T - α tdt, {\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (x) = {\ frac {1} {\ pi i}} \ int _ {- \ infty} ^ {+ \ infty + i \ pi} e ^ {x \ sinh t- \ alpha t} \, dt, \\ H _ {\ alpha} ^ {(2)} (x) = - {\ frac {1} {\ pi i}} \ int _ {- \ infty} ^ {+ \ infty -i \ pi} e ^ {x \ sinh t- \ alpha t} \, dt, \ end {align}}}{\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (x) = {\ frac {1} {\ pi i}} \ int _ {- \ infty} ^ {+ \ infty + i \ pi} e ^ {x \ sinh t- \ alpha t} \, dt, \\ H _ {\ alpha} ^ {(2)} (x) = - {\ frac {1} {\ pi i}} \ int _ {- \ infty} ^ {+ \ infty -i \ pi} e ^ {x \ sinh t- \ alpha t} \, dt, \ end {align}}}

где пределы интегрирования укажите интегрирование вдоль контура , который можно выбрать следующим образом: от −∞ до 0 по отрицательной действительной оси, от 0 до ± iπ по мнимой оси и от ± iπ до + ∞ ± iπ по контуру, параллельному действительной оси.

Модифицированные функции Бесселя: I α, K α

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

I α (x) знак равно я - α J α (ix) знак равно ∑ м знак равно 0 ∞ 1 м! Γ (m + α + 1) (x 2) 2 m + α, K α (x) = π 2 I - α (x) - I α (x) sin ⁡ α π, {\ displaystyle {\ begin {выровнено) } I _ {\ alpha} (x) = i ^ {- \ alpha} J _ {\ alpha} (ix) = \ sum _ {m = 0} ^ {\ infty} {\ frac {1} {m ! \, \ Gamma (m + \ alpha +1)}} \ left ({\ frac {x} {2}} \ right) ^ {2m + \ alpha}, \\ K _ {\ alpha} (x) = {\ frac {\ pi} {2}} {\ frac {I _ {- \ alpha} (x) -I _ {\ alpha} (x)} {\ sin \ alpha \ pi}}, \ end { align}}}{\ displaystyle {\ begin {align} I _ {\ alpha} (x) = i ^ {- \ alpha} J _ {\ alpha} (ix) = \ sum _ {m = 0} ^ {\ infty} {\ frac {1} {m! \, \ Gamma (m + \ alpha +1)}} \ left ({\ frac {x} {2}} \ right) ^ {2m + \ alpha}, \\ K _ {\ alpha} (x) = {\ frac {\ pi} {2}} {\ frac {I _ {- \ alpha} (x) -I _ {\ alpha} (x)} {\ sin \ альфа \ пи}}, \ конец {выровнено}}}

когда α не целое; когда α - целое число, используется предел. Они выбраны с действительными значениями для действительных и положительных аргументов x. Таким образом, разложение в ряд для I α (x) аналогично разложению для J α (x), но без переменного (-1) множителя.

K α {\ displaystyle K _ {\ alpha}}{\ displaystyle K _ {\ alpha}} можно выразить через функции Ганкеля:

K α = {π 2 i α + 1 H α (1) (ix) - π < arg ⁡ x ≤ π 2 π 2 ( − i) α + 1 H α ( 2) ( − i x) − π 2 < arg ⁡ x ≤ π {\displaystyle K_{\alpha }={\begin{cases}{\frac {\pi }{2}}i^{\alpha +1}H_{\alpha }^{(1)}(ix)-\pi <\arg x\leq {\frac {\pi }{2}}\\{\frac {\pi }{2}}(-i)^{\alpha +1}H_{\alpha }^{(2)}(-ix)-{\frac {\pi }{2}}<\arg x\leq \pi \end{cases}}}{\ Displaystyle K _ {\ альфа} = {\ begin {case} {\ frac {\ pi} {2}} i ^ {\ alpha +1} H _ {\ alpha} ^ {(1)} (ix) - \ pi <\ arg x \ leq {\ frac {\ pi} {2}} \\ {\ frac {\ pi} {2}} (- i) ^ {\ alpha +1} H _ {\ alpha} ^ {(2)} (- ix) - {\ frac {\ pi} {2}} <\ arg x \ leq \ pi \ end {cases}}}

Мы можем выразить первую и вторую функции Бесселя через модифицированные функции Бесселя (они верны, если −π < arg z ≤ π/2):

J α (iz) = e α i π 2 I α (z), Y α (iz) знак равно е (α + 1) я π 2 я α (z) - 2 π е - α я π 2 К α (z). {\ Displaystyle {\ begin {выровнено} J _ {\ альфа} (iz) = e ^ {\ frac {\ alpha i \ pi} {2}} I _ {\ alpha} (z), \\ Y _ {\ alpha} (iz) = e ^ {\ frac {(\ alpha +1) i \ pi} {2}} I _ {\ alpha} (z) - {\ frac {2} {\ pi}} e ^ {- {\ frac {\ alpha i \ pi} {2}}} K_ {\ alpha} (z). \ end {align}}}{\ displaystyle {\ begin {align} J _ {\ alpha} (iz) = e ^ {\ frac {\ alpha i \ pi} {2}} I _ {\ alpha} (z), \\ Y _ {\ alpha} (iz) = e ^ {\ frac {(\ alpha +1) i \ pi} { 2}} I _ {\ alpha} (z) - {\ frac {2} {\ pi}} e ^ {- {\ fra с {\ альфа я \ пи} {2}}} К _ {\ альфа} (г). \ конец {выровнено}}}

Iα(x) и K α (x) - два линейно независимых решения модифицированного уравнения Бесселя :

x 2 d 2 ydx 2 + xdydx - (x 2 + α 2) y = 0. {\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + x {\ frac {dy} {dx}} - \ left (x ^ {2} + \ alpha ^ {2} \ right) y = 0.}{\ displaystyle x ^ { 2} {\ frac {d ^ {2} y} {dx ^ {2}}} + x {\ frac {dy} {dx}} - \ left (x ^ {2} + \ alpha ^ {2} \ вправо) y = 0.}

В отличие от обычных функций Бесселя, которые колеблются как функции дей ствительного аргумента, I α и K α - это экспоненциально растущие и убывающие функции соответственно. Как и обычная функция Бесселя J α, функция I α стремится к нулю при x = 0 при α>0 и конечна при x = 0 при α = 0. Аналогично, K α расходится при x = 0 с особенностью логарифмического типа для K 0 и ½Γ (| α |) (2 / x) в случае.

Модифицированные функции Бесселя первого рода, I α (x), для α = 0, 1, 2, 3 Модифицированные функции Бесселя второго рода, K α (x), для α = 0, 1, 2, 3

. Две интегральные формулы для модифицированных функций Бесселя (для Re (x)>0):

I α (x) = 1 π ∫ 0 π ex cos ⁡ θ cos Α θ d θ - sin ⁡ α π π ∫ 0 ∞ e - x ch ⁡ t - α tdt, K α (x) = ∫ 0 ∞ e - x ch ⁡ t ch ⁡ α tdt. {\ displaystyle {\ begin {align} I _ {\ alpha} (x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} e ^ {x \ cos \ theta } \ cos \ alpha \ theta \, d \ theta - {\ frac {\ sin \ alpha \ pi} {\ pi}} \ int _ {0} ^ {\ infty} e ^ {- x \ ch t- \ альфа t} \, dt, \\ K _ {\ alpha} (x) = \ int _ {0} ^ {\ infty} e ^ {- x \ ch t} \ cosh \ alpha t \, dt. \ end {выровнены}}}{\ displaystyle {\ begin {align} I _ {\ alpha} ( x) = {\ frac {1} {\ pi}} \ int _ {0} ^ {\ pi} e ^ {x \ cos \ theta} \ cos \ alpha \ theta \, d \ theta - {\ frac {\ sin \ alpha \ pi} {\ pi}} \ int _ {0} ^ {\ infty} e ^ {- x \ ch t- \ alpha t} \, dt, \\ K _ {\ alpha} (x) = \ int _ {0} ^ {\ infty} e ^ {- x \ ch t} \ cosh \ alpha t \, dt. \ end {align}}}

Функции Бесселя можно описать как преобразование Фурье степеней квадратичных функций. Например:

2 K 0 (ω) = ∫ - ∞ ∞ e i ω t t 2 + 1 d t. {\ displaystyle 2 \, K_ {0} (\ omega) = \ int _ {- \ infty} ^ {\ infty} {\ frac {e ^ {я \ omega t}} {\ sqrt {t ^ {2} +1}}} \, dt.}{\ displaystyle 2 \, K_ {0} (\ omega) = \ int _ {- \ infty} ^ {\ infty} {\ frac {е ^ {я \ omega t}} {\ sqrt {t ^ {2} +1}}} \, dt.}

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

Модифицированные функции Бесселя K 1/3 и K 2/3 могут быть представлены в виде быстро сходящихся интегралов

K 1 3 (ξ) = 3 ∫ 0 ∞ exp ⁡ (- ξ (1 + 4 x 2 3) 1 + x 2 3) dx, K 2 3 (ξ) = 1 3 ∫ 0 ∞ 3 + 2 x 2 1 + x 2 3 exp ⁡ (- ξ (1 + 4 x 2 3) 1 + x 2 3) dx. {\ displaystyle {\ begin {align} K _ {\ frac {1} {3}} (\ xi) = {\ sqrt {3}} \ int _{0} ^ {\ infty} \ exp \ left (- \ xi \ left (1 + {\ frac {4x ^ {2}} {3}} \ right) {\ sqrt {1 + {\ frac {x ^) {2}} {3}}}} \ right) \, dx, \\ K _ {\ frac {2} {3}} (\ xi) = {\ frac {1} {\ sqrt {3}} } \ int _ {0} ^ {\ infty} {\ frac {3 + 2x ^ {2}} {\ sqrt {1 + {\ frac {x ^ {2}} {3}}}}} \ exp \ left (- \ xi \ left (1 + {\ frac {4x ^ {2}} {3}} \ right) {\ sqrt {1 + {\ frac {x ^ {2}} {3}}}} \ вправо) \, dx. \ End {align}}}{\ displaystyle {\ begin {align} K _ {\ frac {1} {3}} (\ xi) = { \ sqrt {3}} \ int _ {0} ^ {\ infty} \ exp \ left (- \ xi \ left (1 + {\ frac {4x ^ {2}} {3}} \ right) {\ sqrt {1 + {\ frac {x ^ {2}} {3}}}} \ right) \, dx, \\ K _ {\ frac {2} {3}} (\ xi) = {\ frac {1 } {\ sqrt {3}}} \ int _ {0} ^ {\ infty} {\ frac {3 + 2x ^ {2}} {\ sqrt {1 + {\ frac {x ^ {2}} {3 }}}}} \ exp \ left (- \ xi \ left (1 + {\ frac {4x ^ {2}} {3}} \ right) {\ sqrt {1 + {\ frac {x ^ {2} } {3}}}} \ right) \, dx. \ End {align}}}

модифицированная функция Бесселя второго также вызывалась под именами (сейчас редко):

Сферические функции Бесселя: j n, y n Сферические функции Бесселя первого вида, j n (x), для n = 0, 1, 2 Сферические функции Бесселя второго рода, y n (x), для n = 0, 1, 2

При решении уравнения Гельмгольца в сферической координате путем разделения альтернативное уравнение имеет вид

x 2 d 2 ydx 2 + 2 xdydx + (x 2 - n (n + 1)) y = 0. {\ displaystyle x ^ {2} {\ гидроразрыв {d ^ {2} y} {dx ^ {2}}} + 2x {\ frac {dy} {dx}} + \ left (x ^ {2} -n (n + 1) \ right) y = 0.}{\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + 2x {\ frac {dy} {dx}} + \ слева (x ^ {2} -n (n + 1) \ right) y = 0.}

Два линейно независимых решения этого уравнения сферические функции Бесселя jnи y n и с вязаны с обычными функциями Бесселя J n и Y nby

jn (x) = π 2 x J n + 1 2 (x), yn (x) = π 2 x Y n + 1 2 (x) = (- 1) n + 1 π 2 х Дж - п - 1 2 (х). {\ displaystyle {\ begin {align} j_ {n} (x) = {\ sqrt {\ frac {\ pi} {2x}}} J_ {n + {\ frac {1} {2}}} (x), \\ y_ {n} (x) = {\ sqrt {\ frac {\ pi} {2x}}} Y_ {n + {\ frac {1} {2}}} (x) = (- 1) ^ {n + 1} {\ sqrt {\ frac {\ pi} {2x}}} J _ {- n - {\ frac {1} {2}}} (x). \ end {align}}}{\ displaystyle {\ begin {align} j_ {n} (x) = {\ sqrt {\ frac {\ pi} {2x}}} J_ {n + {\ frac {1} {2}}} (x), \\ y_ {n} (x) = {\ sqrt {\ frac {\ pi} {2x}}} Y_ {n + {\ frac {1} {2}}} (x) = (- 1) ^ {n + 1} {\ sqrt {\ frac {\ pi} {2x}}} J _ {- n - {\ frac {1} {2}}} (x). \ End {align}}}

ynтакже обозначается n n или η n; некоторые авторы эти функции являются сферическими функциями Неймана .

Сферические функции Бесселя также могут быть записаны как (формулы Рэлея )

jn (x) = (- x) n (1 xddx) n sin ⁡ хх, уn (х) = - (- х) п (1 хддх) п соз ⁡ хх. {\ displaystyle {\ begin {align} j_ {n} (x) = (- x) ^ {n} \ left ({ \ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {n} {\ frac {\ sin x} {x}}, \\ y_ {n} (x) = - (- x) ^ {n} \ left ({\ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {n} {\ frac {\ cos x} {x} }. \ end {align}}}{\ displaystyle {\ begin {a ligned} j_ {n} (x) = (- x) ^ {n} \ left ({\ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {n} { \ frac {\ sin x} {x}}, \\ y_ {n} (x) = - (- x) ^ {n} \ left ({\ frac {1} {x}} {\ frac {d } {dx}} \ right) ^ {n} {\ frac {\ cos x} {x}}. \ end {align}}}

Первая сферическая функция Бесселя j 0 (x) также известна как (ненормализованная) функция sinc. Первые сферические функции Бесселя являются:

j 0 (x) = sin ⁡ xx. J 1 (x) = sin ⁡ xx 2 - cos ⁡ xx, j 2 (x) = (3 x 2 - 1) sin ⁡ xx - 3 соз ⁡ хх 2, j 3 (х) = (15 х 3 - 6 х) грех хх - (15 х 2 - 1) соз х хх {\ displaystyle {\ begin {align} j_ {0} (x) = {\ frac { \ sin x} {x}}. \\ j_ {1} (x) = {\ frac {\ sin x} {x ^ {2}}} - {\ frac {\ cos x} { x}}, \\ j_ {2} (x) = \ left ({\ frac {3} {x ^ {2}}} - 1 \ right) {\ frac {\ sin x} {x}} - {\ frac {3 \ cos x} {x ^ {2}}}, \\ j_ {3} (x) = \ left ({\ frac {15} {x ^ {3}}} - {\ frac {6} {x}} \ right) {\ frac {\ sin x} {x}} - \ left ({\ frac {15} {x ^ {2}}} - 1 \ right) {\ frac {\ cos x} {x}} \ end {align}}}{\ displaystyle {\ begin {align} j_ {0} (x) = {\ frac {\ sin x} {x}}. \\ j_ {1} (x) = {\ frac {\ sin x} {x ^ {2}}} - {\ frac {\ cos x} {x}}, \\ j_ {2} (x) = \ left ({\ frac {3} {x ^ {2}}} - 1 \ right) {\ frac {\ sin x} {x}} - {\ frac {3 \ cos x} {x ^ {2}} }, \\ j_ {3} (x) = \ left ({\ frac {15} {x ^ {3}}} - {\ frac {6} {x}} \ right) {\ frac {\ sin x} {x}} - \ left ({\ frac {15} {x ^ {2}}} - 1 \ right) {\ frac {\ cos x} {x}} \ end {align}}}

и

y 0 (x) = - j - 1 (x) = - cos ⁡ xx, y 1 (x) = j - 2 ( x) = - cos ⁡ xx 2 - sin ⁡ xx, y 2 (x) = - j - 3 (x) = (- 3 x 2 + 1) cos ⁡ xx - 3 sin ⁡ xx 2, y 3 (x) знак равно j - 4 (x) = (- 15 x 3 + 6 x) cos ⁡ xx - (15 x 2 - 1) sin ⁡ xx. {\ displaystyle {\ begin {align} y_ {0} (x) = - j _ {- 1} (x) = - {\ frac {\ cos x} {x}}, \\ y_ {1} ( x) = j _ {- 2} (x) = - {\ frac {\ cos x} {x ^ {2}}} - {\ frac {\ sin x} {x}}, \\ y_ {2 } (x) = - j _ {- 3} (x) = \ left (- {\ frac {3} {x ^ {2}}} + 1 \ right) {\ frac {\ cos x} {x }} - {\ frac {3 \ sin x} {x ^ {2}}}, \\ y_ {3} (x) = j _ {- 4} (x) = \ left (- {\ frac { 15} {x ^ {3}}} + {\ frac {6} {x}} \ right) {\ frac {\ cos x} {x}} - \ left ({\ frac {15} {x ^ { 2}}} -1 \ right) {\ frac {\ sin x} {x}}. \ End {align}}}{\ displaystyle {\ begin {align} y_ {0} (x) = - j _ {- 1} (x) = - {\ frac {\ cos x} {x}}, \\ y_ {1} (x) = j _ {- 2} (x) = - {\ frac {\ cos x} {x ^ {2}}} - {\ frac {\ sin x} {x}}, \\ y_ {2} (x) = - j _ {- 3} (x) = \ left (- {\ frac {3} {x ^ {2}}} + 1 \ right) {\ frac {\ cos x} {x}} - {\ frac {3 \ sin x} {x ^ {2} }}, \\ y_ {3} (x) = j_ {-4} (x) = \ left (- {\ frac {15} {x ^ {3}}} + {\ frac {6} {x}} \ right) {\ frac {\ cos x} {x }} - \ left ({\ frac {15} {x ^ {2}}} - 1 \ right) {\ frac {\ sin x} {x}}. \ end {align}}}

Производящая функция

Сферические функции Бесселя имеют производящие функции

1 z cos ⁡ (z 2 - 2 zt) знак равно ∑ N знак равно 0 ∞ tnn! J N - 1 (Z), 1 Z грех ⁡ (Z 2 - 2 Z T) знак равно ∑ N знак равно 0 ∞ T N N! у н - 1 (г). {\ displaystyle {\ begin {align} {\ frac {1} {z}} \ cos \ left ({\ sqrt {z ^ {2} -2zt}} \ right) = \ sum _ {n = 0} ^ {\ infty} {\ frac {t ^ {n}} {n!}} j_ {n-1} (z), \\ {\ frac {1} {z}} \ sin \ left ({\ sqrt {z ^ {2} -2zt}} \ right) = \ sum _ {n = 0} ^ {\ infty} {\ frac {t ^ {n}} {n!}} y_ {n-1} ( z). \ end {align}}}{\ displaystyle {\ begin {align} {\ frac {1} {z}} \ cos \ left ({\ sqrt {z ^ {2} -2zt}} \ right) = \ sum _ {n = 0} ^ {\ infty} {\ frac {t ^ { n}} {n!}} j_ {n-1} (z), \\ {\ frac {1} {z}} \ sin \ left ({\ sqrt {z ^ {2} -2zt}} \ right) = \ sum _ {n = 0} ^ {\ infty} {\ frac {t ^ {n}} {n!}} y_ {n-1} (z). \ end {align}}}

Дифференциальные отношения

Далее f n представляет собой любое из j n, y n, h. n, h. nдля n = 0, ± 1, ± 2,...

(1 zddz) m (zn + 1 fn (z)) = zn - m + 1 fn - m (z), (1 zddz) m (z - nfn (z)) = (- 1) mz - n - mfn + m (z). {\ displaystyle {\ begin {align} \ left ({\ frac {1} {z}} {\ frac {d} {dz}} \ right) ^ {m} \ left (z ^ {n + 1} f_ {n} (z) \ right) = z ^ {nm + 1} f_ {nm} (z), \\\ left ({\ frac {1} {z}} {\ frac {d} {dz} } \ right) ^ {m} \ left (z ^ {- n} f_ {n} (z) \ right) = (- 1) ^ {m} z ^ {- nm} f_ {n + m} ( z). \ end {align}}}{\ displaystyle {\ begin {align} \ left ({\ frac {1} {z}} {\ frac {d} {dz}} \ right) ^ {m} \ left (z ^ {n + 1} f_ {n} (z) \ right) = z ^ {n-m + 1} f_ {nm} (z), \\\ left ({\ frac {1} {z}} {\ frac { d} {dz}} \ right) ^ {m} \ left (z ^ {- n} f_ {n} (z) \ right) = (- 1) ^ {m} z ^ {- nm} f_ { п + м} (г). \ конец {выровнено}}}

Сферические функции Ганкеля: h. n, h. n

Существуют также сферические функции-аналоги Ханкеля:

hn (1) (x) = jn (х) + iyn (x), hn (2) (x) = jn (x) - iyn (x). {\ displaystyle {\ begin {align} h_ {n} ^ {(1)} (x) = j_ {n} (x) + iy_ {n} (x), \\ h_ {n} ^ {(2)} (х) = j_ {n} (x) -iy_ {n} (x). \ end {align}}}{\ displaystyle {\ begin {выровнено} h_ {n} ^ {(1)} (x) = j_ {n} (x) + iy_ {n} (x), \\ h_ {n} ^ {(2)} (x) = j_ {n} (x) -iy_ {n} (x). \ end {выравнивается}}}

На самом деле есть простые выражения в замкнутой форме для функций Бесселя полуцелочисленный порядок в терминах стандартных тригонометрических функций и, следовательно, для сферических функций Бесселя. В частности, для неотрицательных целых чисел n:

h n (1) (x) = (- i) n + 1 e i x x ∑ m = 0 n i m m! (2 х) м (п + м)! (п - м)!, {\ displaystyle h_ {n} ^ {(1)} (x) = (- i) ^ {n + 1} {\ frac {e ^ {ix}} {x}} \ sum _ {m = 0} ^ {n} {\ frac {i ^ {m}} {м! \, (2x) ^ {m}}} {\ frac {(n + m)!} {(Nm)!}},}{\ displaystyle h_ {n} ^ {(1)} (x) = (- i) ^ {n + 1} {\ frac {e ^ {ix}} {x}} \ sum _ {m = 0} ^ {n} {\ frac {i ^ {m}} {m! \, (2x) ^ {m}}} {\ frac {(n + m)!} {(нм)!} },}

и h. n- это комплексно-сопряженное значение этого (для действительного x). Отсюда следует, например, что j 0 (x) = sin x / x и y 0 (x) = −cos x / x и так далее. мультипольном разложении электромагнитного поля.

Функции Риккати - Бесселя - Бесселя: S n мультипольном разложении электромагнитного поля

., C n, ξ n, ζ n

Riccati - функции Бесселя незначительно отличаются от сферических функций Бесселя:

S n (x) = xjn (x) = π x 2 J n + 1 2 (x) C n (x) = - xyn (x) = - π x 2 Y n + 1 2 (x) ξ n (x) = xhn (1) (x) = π x 2 H n + 1 2 (1) (x) = S n (x) - i C n (x) ζ n (x) = xhn (2) (x) = π x 2 H n + 1 2 ( 2) (Икс) знак равно SN (Икс) + я С N (Икс) {\ Displaystyle {\ begin {Выровнено} S_ {п} (х) = xj_ {п} (х) = {\ sqrt {\ frac {\ pi x} {2}}} J_ {n + {\ frac {1} {2}}} (x) \\ C_ {n} (x) = - xy_ {n} (x) = - { \ sqrt {\ frac {\ pi x} {2}}} Y_ {n + {\ frac {1} {2}}} (x) \\\ xi _ {n} (x) = xh_ {n} ^ {(1)} (x) = {\ sqrt {\ frac {\ pi x} {2}}} H_ {n + {\ frac {1} {2}}} ^ {(1)} (x) = S_ {n} (x) -iC_ {n} (x) \\\ zeta _ {n} (x) = xh_ {n} ^ { (2)} (x) = {\ sqrt {\ frac {\ pi x} {2}}} H_ {n + {\ frac {1} {2}}} ^ {(2)} (x) = S_ {n} (x) + iC_ {n} (x) \ end {align}}}{\ displaystyle {\ begin {align} S_ {n} (x) = xj_ {n} (x) = {\ sqrt {\ frac {\ pi x} {2}}} J_ {n + {\ frac { 1} {2}}} (x) \\ C_ {n} (x) = - xy_ {n} (x) = - {\ sqrt {\ frac {\ pi x} {2}}} Y_ {n + {\ frac {1} {2}}} (x) \\\ xi _ {n} (x) = xh_ {n} ^ {(1)} (x) = {\ sqrt {\ frac {\ pi x} {2}}} H_ {n + {\ frac {1} {2}}} ^ {(1)} (x) = S_ {n} (x) -iC_ {n} (x) \\\ zeta _ {n} (x) = xh_ {n} ^ {(2)} (x) = {\ sqrt {\ frac {\ pi x} {2}}} H_ {n + {\ frac {1} {2 }}} ^ {(2)} (x) = S_ {n} (x) + iC_ {n} (x) \ end {align}}}

Они удовлетворяют дифференциальному уравнению

x 2 d 2 ydx 2 + (x 2 - n (n + 1)) y Знак равно 0. {\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + \ left (x ^ {2} -n (n + 1) \ right) y = 0.}{\ displaystyle x ^ {2} {\ frac {d ^ {2} y} {dx ^ {2}}} + \ left (x ^ {2} -n (n + 1) \ справа) y = 0.}

Например, такого рода дифференциальное уравнение появляется в квантовой механике при решении радиальной составляющей уравнения уравнения Шредингера с гипотетическим цилиндрическим бесконечным потенциальным барьером. Это дифференциальное уравнение и решения Риккати - Бесселя также задает задаче рассеяния электромагнитных волн сферой, известной как рассеяние Ми после первого решения Ми (1908). См. Например, Du (2004) для последних разработок и ссылок.

После Дебая (1909) иногда используется обозначение ψ n, χ n вместо S n, C n.

Асимптотические формы

Функции Бесселя имеют следующие асимптотические формы. Для малых аргументов 0 < z ≪ √α + 1, one obtains, when α is not a negative integer:

J α (z) ∼ 1 Γ (α + 1) (z 2) α. {\ Displaystyle J _ {\ alpha} (z) \ sim {\ frac {1} {\ Gamma (\ alpha +1)}} \ left ({\ frac {z} {2}} \ right) ^ {\ alpha}.}{\ displaystyle J _ {\ alpha} (z) \ sim {\ frac {1} {\ Gamma (\ alpha +1)}} \ left ({\ frac {z} {2}} \ right) ^ {\ alpha}.}

Когда α - отрицательное целое число, мы имеем

J α (z) ∼ (- 1) α (- α)! (2 z) α. {\ displaystyle J _ {\ alpha} (z) \ sim {\ frac {(-1) ^ {\ alpha}} {(- \ alpha)!}} \ left ({\ frac {2} {z}} \) ^ {\ alpha}.}{\ displaystyle J _ {\ alpha} (z) \ sim {\ frac {(-1) ^ {\ alpha}} {(- \ alpha)!}} \ left ({\ frac {2} {z}} \ right) ^ {\ alpha}.}

Для функций Бесселя второго имеет три случая справа:

Y α (z) ∼ {2 π (ln ⁡ (z 2) + γ), если α = 0 - Γ (α) π (2 z) α + 1 Γ (α + 1) (z 2) α cot ⁡ (α π), если α не является целым неположительным числом (один член доминирует, если α не является мнимым), - (- 1) α Γ (- α) π (z 2) α, если α - целое отрицательное число, {\ displaystyle Y _ {\ alpha} (z) \ sim {\ begin {cases} {\ dfrac {2} {\ pi}} \ left (\ ln \ left ({\ dfrac {z} {2}} \ right) + \ gamma \ right) {\ text {if}} \ alpha = 0 \\ - {\ dfrac {\ Gamma (\ alpha)} {\ pi}} \ left ({\ dfrac {2} {z}} \ right) ^ {\ alpha} + {\ dfrac {1} {\ Gamma (\ alpha +1) }} \ left ({\ dfrac {z} {2}} \ right) ^ {\ alpha} \ cot (\ alpha \ pi) {\ text {if}} \ alpha {\ text {не является неположительным целым числом (один член доминирует, если}} \ alpha {\ text {не является воображаемым)}}, \\ - {\ dfrac {(-1) ^ {\ alpha} \ Gamma (- \ alpha)} {\ pi}} \ left ({\ dfrac {z} {2}} \ right) ^ {\ alpha} {\ text {if}} \ alpha {\ text { - отрицательное целое число,}} \ end {cases}}}{\ displaystyle Y _ {\ alpha } (z) \ sim {\ begin {cases} {\ dfrac {2} {\ pi}} \ left (\ ln \ left ({\ dfrac {z} {2}} \ right) + \ gamma \ right) {\ text {if}} \ alpha = 0 \\ - {\ dfrac {\ Gamma (\ alpha)} {\ pi}} \ left ({\ dfrac {2} {z}} \ right) ^ {\ alpha} + {\ dfrac {1} {\ Gamma (\ alpha +1)}} \ left ({\ dfrac {z} {2}} \ right) ^ {\ alpha} \ cot (\ alpha \ pi) {\ text {if}} \ alpha {\ text {не является целым неположительным числом (один член доминирует, если только}} \ alpha {\ text {не являетсявоображаемым)}}, \\ - {\ dfrac {(-1) ^ {\ а lpha} \ Gamma (- \ alpha)} {\ pi}} \ left ({\ dfrac {z} {2}} \ right) ^ {\ alpha} {\ text {if}} \ alpha {\ text { является отрицательным целым числом,}} \ end {cases}}}

где γ - константа Эйлера - Маскерони (0,5772...).

Для больших вещественных аргументов z ≫ | α - 1/4 | нельзя записать истинную асимптотику для функций Бесселя и второго рода (если α не является полуцелым числом ), потому что они имеют нули до бесконечности, что должно было бы точно соответствовать любому асимптотическому разложению. Однако для данного значения arg z можно написать уравнение, содержащее член порядка | z |:

J α (z) = 2 π z (cos ⁡ (z - α π 2 - π 4) + e | Im ⁡ (z) | O (| z | - 1)) для | arg ⁡ z | < π, Y α ( z) = 2 π z ( sin ⁡ ( z − α π 2 − π 4) + e | Im ⁡ ( z) | O ( | z | − 1)) for | arg ⁡ z | < π. {\displaystyle {\begin{aligned}J_{\alpha }(z)={\sqrt {\frac {2}{\pi z}}}\left(\cos \left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)+e^{\left|\operatorname {Im} (z)\right|}\mathrm {O} \left(|z|^{-1}\right)\right){\text{for }}\left|\arg z\right|<\pi,\\Y_{\alpha }(z)={\sqrt {\frac {2}{\pi z}}}\left(\sin \left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)+e^{\left|\operatorname {Im} (z)\right|}\mathrm {O} \left(|z|^{-1}\right)\right){\text{for }}\left|\arg z\right|<\pi.\end{aligned}}}{\ displaystyle {\ begin {align} J_ { \ alpha} (z) = {\ sqrt {\ frac {2} {\ pi z}}} \ left (\ cos \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right) + e ^ {\ left | \ operatorname {Im} (z) \ right |} \ mathrm {O} \ left (| z | ^ {- 1} \ right) \ right) {\ text {for}} \ left | \ arg z \ right | <\ pi, \\ Y _ {\ alpha} (z) = {\ sqrt {\ frac {2} {\ pi z} }} \ left (\ sin \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right) + e ^ {\ left | \ operatorname { Im} (z) \ right |} \ mathrm {O} \ left (| z | ^ {- 1} \ right) \ right) {\ text {for}} \ left | \ arg z \ right | <\ пи. \ конец {выровнено}}}

(При α = 1/2 члены в этих формулах полностью выпадают; см. Сферические функции Бесселя выше.) Даже если эти уравнения верны, для комплексного z могут быть доступны лучшие приближения. Например, J 0 (z), когда z находится рядом с отрицательной действительной линией, лучше аппроксимируется

J 0 (z) ≈ - 2 π z cos ⁡ (z + π 4) {\ displaystyle J_ {0} (z) \ приблизительно {\ sqrt {\ frac {-2} {\ pi z}}} \ cos \ left (z + {\ frac {\ pi} {4}} \ right)}{\ displaystyle J_ {0} (z) \ приблизительно {\ sqrt {\ frac {-2 } {\ pi z}}} \ cos \ left (z + {\ frac {\ pi} {4}} \ right)}

чем на

J 0 (z) ≈ 2 π z cos ⁡ (z - π 4). {\ displaystyle J_ {0} (z) \ приблизительно {\ sqrt {\ frac {2} {\ pi z}}} \ cos \ left (z - {\ frac {\ pi} {4}} \ right). }{\ displaystyle J_ {0} (z) \ приблизительно {\ sqrt {\ frac {2} {\ pi z}}} \ cos \ left (z - {\ frac {\ pi} {4}} \ right).}

Асимптотики для функций Ганкеля:

H α (1) (z) ∼ 2 π zei (z - α π 2 - π 4) для - π < arg ⁡ z < 2 π, H α ( 2) ( z) ∼ 2 π z e − i ( z − α π 2 − π 4) for − 2 π < arg ⁡ z < π. {\displaystyle {\begin{aligned}H_{\alpha }^{(1)}(z)\sim {\sqrt {\frac {2}{\pi z}}}e^{i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}-\pi <\arg z<2\pi,\\H_{\alpha }^{(2)}(z)\sim {\sqrt {\frac {2}{\pi z}}}e^{-i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}-2\pi <\arg z<\pi.\end{aligned}}}{\ displaystyle {\ begin {align} H _ {\ alpha} ^ {(1)} (z) \ sim {\ sqrt { \ frac {2} {\ pi z}}} e ^ {i \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right)} {\ text {for}} - \ pi <\ arg z <2 \ pi, \\ H _ {\ alpha} ^ {(2)} (z) \ sim {\ sqrt {\ frac {2} {\ pi z}}} e ^ {- i \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right)} {\ text {для }} - 2 \ pi <\ arg z <\ pi. \ End {align}}}

Их можно распространить на другие значения arg z с использованием уравнения, связывающих H. α(ze) и H. α(ze) с H. α(z) и H. α(z).

Интересно, что хотя функция Бесселя первого - это среднее значение двух функций Ганкеля, J α (z) не является асимптотическим средним средним этих двух асимптотических форм, когда z отрицательно (потому что одна или другая функция) другая не будет там будет правильно, в зависимости от используемого arg z). Но асимптотики для функций Ганкеля позволяют нам записывать асимптотики для функций Бесселя первого и второго рода для комплексных (не действующих) z, пока | z | стремится к бесконечности при постоянном фазовом угле arg z (с использованием квадратного корня с положительной действительной частью):

J α (z) ∼ 1 2 π zei (z - α π 2 - π 4) для - π < arg ⁡ z < 0, J α ( z) ∼ 1 2 π z e − i ( z − α π 2 − π 4) for 0 < arg ⁡ z < π, Y α ( z) ∼ − i 1 2 π z e i ( z − α π 2 − π 4) for − π < arg ⁡ z < 0, Y α ( z) ∼ − i 1 2 π z e − i ( z − α π 2 − π 4) for 0 < arg ⁡ z < π. {\displaystyle {\begin{aligned}J_{\alpha }(z)\sim {\frac {1}{\sqrt {2\pi z}}}e^{i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}-\pi <\arg z<0,\\J_{\alpha }(z)\sim {\frac {1}{\sqrt {2\pi z}}}e^{-i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}0<\arg z<\pi,\\Y_{\alpha }(z)\sim -i{\frac {1}{\sqrt {2\pi z}}}e^{i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}-\pi <\arg z<0,\\Y_{\alpha }(z)\sim -i{\frac {1}{\sqrt {2\pi z}}}e^{-i\left(z-{\frac {\alpha \pi }{2}}-{\frac {\pi }{4}}\right)}{\text{for }}0<\arg z<\pi.\end{aligned}}}{\ displaystyle {\ begin {выровнено } J _ {\ alpha} (z) \ sim {\ frac {1} {\ sqrt {2 \ pi z}}} e ^ {i \ left (z - {\ frac {\ alpha \ pi} {2}) } - {\ frac {\ pi} {4}} \ right)} {\ text {for}} - \ pi <\ arg z <0, \\ J _ {\ alpha} (z) \ sim {\ frac {1} {\ sqrt {2 \ pi z}}} e ^ {- i \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right)} {\ text {for}} 0 <\ arg z <\ pi, \\ Y _ {\ alpha} (z) \ sim -i {\ frac {1} {\ sqrt {2 \ pi z }}} e ^ {i \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right)} {\ text {for}} - \ pi <\ arg z <0, \\ Y _ {\ alpha} (z) \ sim -i {\ frac {1} {\ sqrt {2 \ pi z}}} e ^ {- i \ left (z - {\ frac {\ alpha \ pi} {2}} - {\ frac {\ pi} {4}} \ right)} {\ text {for}} 0 <\ arg z <\ pi. \ end { выровнено}}}

Для модифицированных функций Бесселя Ханкель разработал также асимптотические разложения (с большим аргументом) :

I α (z) ∼ ez 2 π z (1 - 4 α 2 - 1 8 z + (4 α 2 - 1) (4 α 2 - 9) 2! (8 z) 2 - (4 α 2 - 1) (4 α 2 - 9) (4 α 2 - 25) 3! (8 z) 3 + ⋯) для | arg ⁡ z | < π 2, K α ( z) ∼ π 2 z e − z ( 1 + 4 α 2 − 1 8 z + ( 4 α 2 − 1) ( 4 α 2 − 9) 2 ! ( 8 z) 2 + ( 4 α 2 − 1) ( 4 α 2 − 9) ( 4 α 2 − 25) 3 ! ( 8 z) 3 + ⋯) for | arg ⁡ z | < 3 π 2. {\displaystyle {\begin{aligned}I_{\alpha }(z)\sim {\frac {e^{z}}{\sqrt {2\pi z}}}\left(1-{\frac {4\alpha ^{2}-1}{8z}}+{\frac {\left(4\alpha ^{2}-1\right)\left(4\alpha ^{2}-9\right)}{2!(8z)^{2}}}-{\frac {\left(4\alpha ^{2}-1\right)\left(4\alpha ^{2}-9\right)\left(4\alpha ^{2}-25\right)}{3!(8z)^{3}}}+\cdots \right){\text{for }}\left|\arg z\right|<{\frac {\pi }{2}},\\K_{\alpha }(z)\sim {\sqrt {\frac {\pi }{2z}}}e^{-z}\left(1+{\frac {4\alpha ^{2}-1}{8z}}+{\frac {\left(4\alpha ^{2}-1\right)\left(4\alpha ^{2}-9\right)}{2!(8z)^{2}}}+{\frac {\left(4\alpha ^{2}-1\right)\left(4\alpha ^{2}-9\right)\left(4\alpha ^{2}-25\right)}{3!(8z)^{3}}}+\cdots \right){\text{for }}\left|\arg z\right|<{\frac {3\pi }{2}}.\end{aligned}}}{\ displaystyle {\ begin {align} I _ {\ alpha} (z) \ sim {\ frac {e ^ {z}} { \ sqrt {2 \ pi z}}} \ left (1 - {\ frac {4 \ alpha ^ {2} -1} {8z}} + {\ frac {\ left (4 \ alpha ^ {2} -1 \ right) \ left (4 \ alpha ^ {2} -9 \ right)} {2! (8z) ^ {2}}} - {\ frac {\ left (4 \ alpha ^ {2} -1 \ right) \ left (4 \ alpha ^ {2} -9 \ right) \ left (4 \ alpha ^ {2} -25 \ right)} {3! (8z) ^ {3}}} + \ cdots \ right) {\ text {for}} \ left | \ arg z \ right | <{\ frac {\ pi} {2}}, \\ K _ {\ alpha} (z) \ sim {\ sqrt {\ frac { \ pi} {2z}}} e ^ {- z} \ left (1 + {\ frac {4 \ alpha ^ {2} -1} {8z}} + {\ frac {\ left (4 \ alpha ^ { 2} -1 \ right) \ left (4 \ alpha ^ {2} -9 \ right)} {2! (8z) ^ {2}}} + {\ frac {\ left (4 \ alpha ^ {2} -1 \ вправо) \ влево (4 \ альф a ^ {2} -9 \ right) \ left (4 \ alpha ^ {2} -25 \ right)} {3! (8z) ^ {3}}} + \ cdots \ right) {\ text {для }} \ left | \ arg z \ right | <{\ frac {3 \ pi} {2}}. \ end {align}}}

Когда α = 1/2, все члены, кроме первого, обращаются в нуль, и мы имеем

I 1 2 (z) = 2 π z sinh ⁡ (z) ∼ e z 2 π z для | arg ⁡ z | < π 2, K 1 2 ( z) = π 2 z e − z. {\displaystyle {\begin{aligned}I_{\frac {1}{2}}(z)={\sqrt {\frac {2}{\pi z}}}\sinh(z)\sim {\frac {e^{z}}{\sqrt {2\pi z}}}{\text{for }}\left|\arg z\right|<{\tfrac {\pi }{2}},\\K_{\frac {1}{2}}(z)={\sqrt {\frac {\pi }{2z}}}e^{-z}.\end{aligned}}}{\ displaystyle {\ begin {align} I _ {\ frac {1} {2}} (z) = {\ sqrt {\ frac {2} {\ pi z}}} \ sinh (z) \ sim {\ frac {e ^ {z}} {\ sqrt {2 \ pi z}}} {\ текст {for}} \ left | \ arg z \ right | <{\ tfrac {\ pi} {2}}, \\ K _ {\ frac {1} {2}} (z) = {\ sqrt {\ frac {\ pi} {2z}}} e ^ {- z}. \ end {align}}}

Для малых аргументов 0 < |z| ≪ √α + 1, we have

I α (z) ∼ 1 Γ (α + 1) (z 2) α, K α (z) ∼ {- ln ⁡ (z 2) - γ, если α = 0 Γ (α) 2 (2 z) α, если α>0 {\ displaystyle {\ begin {align} I _ {\ alpha} (z) \ sim {\ frac {1} {\ Gamma (\ alpha +1)}} \ left ({\ frac {z} {2}} \ right) ^ {\ alpha}, \\ K _ {\ alpha} (z) \ sim {\ begin {cases} - \ ln \ left ({\ dfrac {z} {2}} \ right) - \ gamma {\ text {if}} \ alpha = 0 \\ {\ frac {\ Gamma (\ alpha)} {2}} \ left ({\ dfrac {2 } {z}} \ right) ^ {\ alpha} {\ text {if}} \ alpha>0 \ end {cases}} \ end {align}}}{\displaystyle {\begin{aligned}I_{\alpha }(z)\sim {\frac {1}{\Gamma (\alpha +1)}}\left({\frac {z}{2}}\right)^{\alpha },\\K_{\alpha }(z)\sim {\begin{cases}-\ln \left({\dfrac {z}{2}}\right)-\gamma {\text{if }}\alpha =0\\{\frac {\Gamma (\alpha)}{2}}\left({\dfrac {2}{z}}\right)^{\alpha }{\text{if }}\alpha>0 \ end {cases} } \ end {align}}}
Приближение полной области с элементарными функциями

Очень хорошее приближение (ошибка ниже 0,3% {\ displaystyle 0.3 \%}{\ displaystyle 0.3 \%} максимального значения 1) функции Бесселя J 0 {\ displaystyle J_ {0}}J_ {0} для произвольного значения аргумента x может быть получено с помощью элементарных функций путем объединения тригонометрическое приближение, работающее для меньших значений x с выражением, содержащим функцию ослабленного косинуса, действительную для больших аргументов с использованием функции плавного перехода 1 1 + (x / 7) 20 {\ displaystyle {\ frac {1} {1} + (x / 7) ^ {20}}}}{\ displaystyle {\ frac {1} {1+ (x / 7) ^ {20} }}} т.е.

J 0 (x) ≈ [1 6 + 1 3 cos ⁡ x 2 + 1 3 cos ⁡ 3 x 2 + 1 6 cos ⁡ x] 1 1 + (x / 7) 20 + 2 π | х | cos ⁡ [x - π 4 sign ⁡ (x)] [1 - 1 1 + (x / 7) 20]. {\ displaystyle J_ {0} (x) \ приблизительно \ left [{\ frac {1} {6}} + {\ frac {1} {3}} \ cos {\ frac {x} {2}} + { \ frac {1} {3}} \ cos {\ frac {{\ sqrt {3}} x} {2}} + {\ frac {1} {6}} \ cos x \ right] {\ frac {1 } {1+ (x / 7) ^ {20}}} + {\ sqrt {\ frac {2} {\ pi | x |}}} \ cos \ left [x - {\ frac {\ pi} {4 }} \ operatorname {sgn} (x) \ right] \ left [1 - {\ frac {1} {1+ (x / 7) ^ {20}}} \ right].}{\ displaystyle J_ {0} (x) \ приблизительно \ left [{\ frac {1} {6}} + {\ frac {1} {3}} \ cos {\ frac {x} {2}} + {\ frac {1} {3}} \ cos {\ frac {{\ sqrt {3}} x} {2}} + {\ frac {1} {6}} \ cos x \ right] {\ frac {1} {1+ (x / 7) ^ {20}}} + {\ sqrt {\ frac {2} {\ pi | x |}}} \ cos \ left [x - {\ frac {\ pi} {4}} \ operatorname {sgn} (x) \ right] \ left [1 - {\ frac { 1} {1+ (x / 7) ^ {20}}} \ right].}
Свойства

Для целочисленного порядка α = n, J n часто определяется через ряд Лорана для производящей функции:

e (x 2) (t - 1 t) Знак равно ∑ N = - ∞ ∞ J N (x) tn {\ displaystyle e ^ {\ left ({\ frac {x} {2}} \ right) \ left (t - {\ frac {1} {t}}) \ right)} = \ sum _ {n = - \ infty} ^ {\ infty} J_ {n} (x) t ^ {n}}{\ displaystyle e ^ {\ left ({\ frac {x} {2}} \ right) \ left (t - {\ frac {1} {t}} \ right)} = \ sum _ {n = - \ infty} ^ {\ inft y} J_ {n} (x) t ^ {n}}

подход, используемый P. А. Хансен в 1843 году. (Это можно обобщить на нецелочисленный порядок с помощью контурного интегрирования или других методов.) Еще одним важным соотношением для целочисленных порядков является разложение Якоби – Ангера :

eiz соз ⁡ ϕ = ∑ N = - ∞ ∞ в J N (z) ein ϕ {\ displaystyle e ^ {iz \ cos \ phi} = \ sum _ {n = - \ infty} ^ {\ infty} i ^ {n} J_ {n} (z) e ^ {in \ phi}}{\ displaystyle e ^ {iz \ cos \ phi} = \ sum _ {n = - \ infty} ^ {\ infty} я ^ {n} J_ {n} (z) e ^ {in \ phi}}

и

e ± iz sin ⁡ ϕ = J 0 (z) + 2 ∑ n = 1 ∞ J 2 n (z) cos ⁡ (2 n ϕ) ± 2 i ∑ n = 0 ∞ J 2 n + 1 (z) sin ⁡ ((2 n + 1) ϕ) {\ displayst yle e ^ {\ pm iz \ sin \ phi} = J_ {0} (z) +2 \ sum _ {n = 1} ^ {\ infty} J_ {2n} (z) \ cos (2n \ phi) \ pm 2i \ sum _ {n = 0} ^ {\ infty} J_ {2n + 1} (z) \ sin ((2n + 1) \ phi)}{\ displaystyle e ^ {\ pm iz \ sin \ phi} = J_ {0} (z) +2 \ sum _ {n = 1} ^ {\ infty} J_ {2n} (z) \ cos (2n \ phi) \ pm 2i \ sum _ {n = 0} ^ {\ infty} J_ {2n + 1} (z) \ sin ((2n + 1) \ phi)}

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

В более общем смысле, ряд

f (z) = a 0 ν J ν (z) + 2 ⋅ ∑ k = 1 ∞ ak ν J ν + k (z) {\ displaystyle f (z) = a_ {0} ^ {\ nu} J _ {\ nu} (z) +2 \ cdot \ sum _ {k = 1} ^ {\ infty} a_ {k} ^ {\ nu} J _ {\ nu + k} (z)}{\ displaystyle f (z) = a_ {0} ^ {\ nu} J _ {\ nu} (z) +2 \ cdot \ sum _ {k = 1} ^ {\ infty} a_ {k} ^ {\ nu} J _ {\ nu + k} (z)}

называется разложением Неймана функции f. Коэффициенты при ν = 0 имеют явный вид

a k 0 = 1 2 π i ∫ | z | знак равно cf (z) О К (z) dz {\ displaystyle a_ {k} ^ {0} = {\ frac {1} {2 \ pi i}} \ int _ {| z | = c} f (z) O_ {k} (z) \, dz}{\ displaystyle a_ {k} ^ {0} = {\ frac {1} {2 \ pi i}} \ int _ {| z | = c} f (z) O_ {k} (z) \, dz}

где O k - многочлен Неймана.

Выбранные функции допускают специальное представление

f (z) = ∑ k = 0 ∞ ak ν J ν + 2 K (z) {\ displaystyle f (z) = \ sum _ {k = 0} ^ {\ infty} a_ {k} ^ {\ nu} J _ {\ nu + 2k} (z)}{\ displaystyle f (z) = \ sum _ {k = 0} ^ {\ infty} a_ {k} ^ {\ nu} J_ {\ nu + 2k} (z)}

с

ak ν = 2 (ν + 2 k) ∫ 0 ∞ f (z) J ν + 2 k (z) zdz {\ displaystyle a_ {k} ^ {\ nu} = 2 ( \ nu + 2k) \ int _ {0} ^ {\ infty} f (z) {\ frac {J _ {\ nu + 2k} (z)} {z}} \, dz}{\ displaystyle a_ {k} ^ {\ nu} = 2 (\ nu + 2k) \ int _ {0} ^ {\ infty} f (z) {\ frac {J _ {\ nu + 2k} (z)} {z }} \, dz}

из-за ортогональности отношение

∫ 0 ∞ J α (z) J β (z) dzz = 2 π sin ⁡ (π 2 (α - β)) α 2 - β 2 {\ displaystyle \ int _ {0} ^ {\ infty } J _ {\ alpha} (z) J _ {\ beta} (z) {\ frac {dz} {z}} = {\ frac {2} {\ pi}} {\ frac {\ sin \ left ({\ frac {\ pi} {2}} (\ alpha - \ beta) \ right)} {\ alpha ^ {2} - \ beta ^ {2}}}}{\ displaystyle \ int _ {0} ^ {\ infty} J _ {\ alpha} (z) J _ {\ beta} (z) {\ frac {dz} { z}} = {\ frac {2} {\ pi}} {\ frac {\ sin \ left ({\ frac {\ pi} {2}} (\ alpha - \ beta) \ right)} {\ alpha ^ {2} - \ beta ^ {2}}}}

В общем, если f имеет точку ветвления около начала координат такой природы, что

f (z) = ∑ k = 0 ak J ν + k (z) {\ displaystyle f (z) = \ sum _ {k = 0} a_ {k} J_ { \ nu + k} (z)}{\ displaystyle f (z) = \ sum _ {k = 0} a_ {k} J _ {\ nu + k} (z)}

, тогда

L {∑ k = 0 ak J ν + к} (s) знак равно 1 1 + s 2 ∑ К знак равно 0 ак (s + 1 + s 2) ν + к {\ displaystyle {\ mathcal {L}} \ left \ {\ sum _ {k = 0} a_ {k} J _ {\ nu + k} \ right \} (s) = {\ frac {1} {\ sqrt {1 + s ^ {2}}}} \ sum _ {k = 0} {\ frac { a_ {k}} {\ left (s + {\ sqrt {1 + s ^ {2}}} \ right) ^ {\ nu + k}}}}{\ displaystyle {\ mathcal {L}} \ left \ {\ sum _ {k = 0} a_ {k} J _ {\ nu + k} \ right \} (s) = {\ frac {1} {\ sqrt {1+ s ^ {2}}}} \ sum _ {k = 0} {\ frac {a_ {k}} {\ left (s + {\ sqrt {1 + s ^ {2}}} \ right) ^ {\ nu + k}}}}

или

∑ k = 0 ak ξ ν + К знак равно 1 + ξ 2 2 ξ L {е} (1 - ξ 2 2 ξ) {\ displaystyle \ sum _ {k = 0} a_ {k} \ xi ^ {\ nu + k} = {\ frac { 1+ \ xi ^ {2}} {2 \ xi}} {\ mathcal {L}} \ {f \} \ left ({\ frac {1- \ xi ^ {2}} {2 \ xi}} \ right)}{\ displaystyle \ sum _ {k = 0} a_ {k} \ xi ^ {\ nu + k} = {\ frac {1+ \ xi ^ {2}} {2 \ xi}} {\ mathcal {L}} \ {f \} \ left ({\ frac {1- \ xi ^ {2}} {2 \ xi}} \ right)}

где L {f} {\ displaystyle {\ mathcal {L}} \ {f \}}{\ displaystyle {\ mathcal {L}} \ {f \}} - это преобразование Лапласа для f.

Другой способ определения функций Бесселя - это метод Пуассона формула представления и формула Мелера-Сонина:

J ν (z) = (z 2) ν Γ (ν + 1 2) π ∫ - 1 1 eizs (1 - s 2) ν - 1 2 ds = 2 (z 2) ν ⋅ π ⋅ Γ (1 2 - ν) ∫ 1 ∞ sin ⁡ zu (u 2 - 1) ν + 1 2 du {\ displaystyle {\ begin {align} J _ {\ nu} (z) = {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {\ nu}} {\ Gamma \ left (\ nu + {\ frac {1} {2}} \ right) {\ sqrt {\ pi}}}} \ int _ {- 1} ^ {1} e ^ {izs} \ left (1-s ^ {2} \ справа) ^ {\ nu - {\ frac {1} {2}}} \, ds \\ [5px] = {\ frac {2} {{\ left ({\ frac {z} {2}} \ right)} ^ {\ nu} \ cdot {\ sqrt {\ pi}} \ cdot \ Gamma \ left ({\ frac {1} {2}} - \ nu \ right)}} \ int _ {1} ^ {\ infty} {\ frac {\ sin zu} {\ left (u ^ {2} -1 \ right) ^ {\ nu + {\ frac {1} {2}}}}} \, du \ end { align}}}{\ displaystyle {\ begin {align} J _ {\ nu} (z) = {\ frac {\ left ({\ frac {z } {2}} \ right) ^ {\ nu}} {\ Gamma \ left (\ nu + {\ frac {1} {2}} \ right) {\ sqrt {\ pi}}}} \ int _ { -1} ^ {1} e ^ {izs} \ left (1-s ^ {2} \ right) ^ {\ nu - {\ frac {1} {2}}} \, ds \\ [5px] = {\ frac {2} {{\ left ({\ frac {z} {2}} \ right)} ^ {\ nu} \ cdot {\ sqrt {\ pi}} \ cdot \ Gamma \ left ({\ frac {1} {2}} - \ nu \ right)}} \ int _ {1} ^ {\ infty} {\ frac {\ sin zu} {\ left (u ^ {2} -1 \ right) ^ {\ nu + {\ frac {1} {2}}}}} \, du \ end {align}}}

где ν>−1/2 и z ∈ C . Эта формула особенно полезна при работе с преобразованием Фурье.

уравнение Бесселя становится эрмитовым (самосопряженным), если оно делится на x, решения должны удовлетворять ортогональности для соответствующих граничных условий. В частности, отсюда следует, что:

∫ 0 1 x J α (xu α, m) J α (xu α, n) dx = δ m, n 2 [J α + 1 (u α, m)] 2 знака равно δ м, N 2 [J α ′ (u α, m)] 2 {\ displaystyle \ int _ {0} ^ {1} xJ _ {\ alpha} \ left (xu _ {\ alpha, m} \ right) J_ {\ alpha} \ left (xu _ {\ alpha, n} \ right) \, dx = {\ frac {\ delta _ {m, n}} {2}} \ left [J _ {\ альфа +1} \ left (u _ {\ alpha, m} \ right) \ right] ^ {2} = {\ frac {\ delta _ {m, n}} {2}} \ left [J _ {\ alpha} '\ left (u_ {\ alpha, m} \ right) \ right] ^ {2}}{\displaystyle \int _{0}^{1}xJ_{\alpha }\left(xu_{\alpha,m}\right)J_{\alpha }\left(xu_{\alpha,n}\right)\,dx={\frac {\delta _{m,n}}{2}}\left[J_{\alpha +1}\left(u_{\alpha,m}\right)\right]^{2}={\frac {\delta _{m,n}}{2}}\left[J_{\alpha }'\left(u_{\alpha,m}\right)\right]^{2}}

где α>−1, δ m, n - дельта Кронекера, а u α, m - m-й ноль из J α (x). Это соотношение ортогональности можно использовать для извлечения коэффициентов в рядх Фурье - Бесселя, где функция раскрывается в базисе функций J α (xu α, m) для фиксированного α и переменного м.

Аналогичное соотношение для сферических функций Бесселя следует немедленно:

∫ 0 1 x 2 j α (xu α, m) j α (xu α, n) dx = δ m, n 2 [j α + 1 (U α, м)] 2 {\ displaystyle \ int _ {0} ^ {1} x ^ {2} j _ {\ alpha} \ left (xu _ {\ alpha, m} \ right) j _ { \ alpha} \ left (xu _ {\ alpha, n} \ right) \, dx = {\ frac {\ delta _ {m, n}} {2}} \ left [j _ {\ alpha +1} \ left (u_ {\ alpha, m} \ right) \ right] ^ {2}}{\ displaystyle \ int _ {0} ^ {1} x ^ {2} j _ {\ alpha} \ left (xu _ {\ alpha, m} \ right) j _ {\ alpha} \ left (xu _ {\ alpha, n} \ right) \, dx = {\ frac {\ delta _ {m, n}} {2}} \ left [j_ { \ alpha +1} \ left (u _ {\ alpha, m} \ right) \ right] ^ {2}}

Если определить функцию товарного вагона x, которая зависит от малого параметра ε, как:

f ε (Икс) знак равно ε rect ⁡ (Икс - 1 ε) {\ Displaystyle f _ {\ varepsilon} (x) = \ varepsilon \ OperatorName {rect} \ left ({\ frac {x-1} {\ varepsilon}} \ right)}{\ displaystyle f _ {\ varepsilon} (x) = \ varepsilon \ operatorname {rect} \ left ({\ frac {x-1} {\ varepsilon}} \ справа)}

(где rect - это функция прямоугольника ), затем его преобразование Ханкеля (любого заданного порядка α>−1/2), g ε (k) приближается к J α (k), когда ε приближается к нулю для любого заданного k. И наоборот, преобразование Ханкеля (того же порядка) для g ε (k) есть f ε (x):

∫ 0 ∞ k J α (kx) g ε (к) dk знак равно е ε (Икс) {\ Displaystyle \ int _ {0} ^ {\ infty} kJ _ {\ alpha} (kx) g _ {\ varepsilon} (k) \, dk = f _ {\ varepsilon} ( x)}{\ displaystyle \ int _ {0} ^ {\ infty} kJ _ {\ alpha} (kx) g _ {\ varepsilon} (k) \, dk = f _ {\ varepsil on} (x)}

который равен нулю везде, кроме около 1. Когда ε приближается к нулю, правая часть приближается к δ (x - 1), где δ - дельта-функция Дирака. Это допускает предел (в распределительном смысле):

∫ 0 ∞ k J α (kx) J α (k) dk = δ (x - 1) {\ displaystyle \ int _ {0 } ^ {\ infty} kJ _ {\ alpha} (kx) J _ {\ alpha} (k) \, dk = \ delta (x-1)}{\ displaystyle \ int _ {0} ^ {\ infty} kJ _ {\ alpha} (kx) J _ {\ alpha} (k) \, dk = \ delta (x-1)}

После замены переменных получается уравнение замыкания:

∫ 0 ∞ Икс J α (ux) J α (vx) dx знак равно 1 u δ (u - v) {\ displaystyle \ int _ {0} ^ {\ infty} xJ _ {\ alpha} (ux) J _ {\ alpha } (vx) \, dx = {\ frac {1} {u}} \ delta (uv)}{\ displaystyle \ int _ {0} ^ {\ infty} xJ _ {\ alpha} (ux) J _ {\ alpha} (vx) \, dx = {\ frac {1} {u}} \ delta (uv)}

для α>−1/2. Преобразование Ханкеля может выразить довольно произвольную функцию как интеграл от функций Бесселя разного масштаба. Для сферических функций Бесселя соотношение ортогональности:

∫ 0 ∞ x 2 j α (ux) j α (vx) dx = π 2 u 2 δ (u - v) {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2} j _ {\ alpha} (ux) j _ {\ alpha} (vx) \, dx = {\ frac {\ pi} {2u ^ {2}}} \ delta (uv)}{\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2} j _ {\ alpha} (ux) j _ {\ alpha} (vx) \, dx = {\ frac {\ pi} {2u ^ {2}}} \ дельта (УФ)}

для α>−1.

Другое важное свойство уравнений Бесселя, которое следует из тождества Абеля, включает вронскиан решений:

A α (x) d B α dx - d A α dx В α (Икс) знак равно С α Икс {\ Displaystyle A _ {\ alpha} (x) {\ frac {dB _ {\ alpha}} {dx}} - {\ frac {dA _ {\ alpha}} {dx}} B _ {\ alpha} (x) = {\ frac {C _ {\ alpha}} {x}}}{\ displaystyle A _ {\ alpha} (x) {\ frac {dB _ {\ alpha}} {dx}} - {\ frac {dA _ {\ alpha}} {dx}} B _ {\ alpha} (x) = {\ frac {C _ {\ alpha}} {x}}}

где A α и B α - любые два решения уравнения Бесселя, и C α является константой, не зависящей от x (которая зависит от α и от конкретных рассматриваемых функций Бесселя). В частности,

J α (x) d Y α dx - d J α dx Y α (x) = 2 π x {\ displaystyle J _ {\ alpha} (x) {\ frac {dY _ {\ alpha}} {dx}} - {\ frac {dJ _ {\ alpha}} {dx}} Y _ {\ alpha} (x) = {\ frac {2} {\ pi x}}}{\ displaystyle J _ {\ alpha} (x) {\ frac {dY _ {\ alpha}} {dx}} - {\ frac {dJ _ {\ alpha}} {dx}} Y _ {\ alpha} (x) = {\ frac {2} {\ pi x}}}

и

I α (Икс) d К α dx - d я α dx К α (Икс) = - 1 Икс, {\ Displaystyle I _ {\ alpha} (x) {\ гидроразрыва {dK _ {\ alpha}} {dx}} - {\ frac {dI _ {\ alpha}} {dx}} K _ {\ alpha} (x) = - {\ frac {1} {x}},}{\ Displaystyle I _ {\ альфа} (х) {\ гидроразрыва {dK_ {\ alpha}} {dx}} - {\ frac {dI _ {\ alpha}} {dx}} K _ {\ alpha} (x) = - {\ frac {1} {x}},}

для α>−1.

При α>−1 четная целая функция рода 1, xJ α (x), имеет только действительные нули. Пусть

0 < j α, 1 < j α, 2 < ⋯ < j α, n < ⋯ {\displaystyle 0{\ displaystyle 0 <j _ {\ alpha, 1} <j _ {\ alpha, 2} <\ cdots <j _ {\ alpha, n} <\ cdots}

- все его положительные нули, тогда

J α (z) = (z 2) α Γ (α + 1) ∏ n = 1 ∞ (1 - z 2 j α, n 2) {\ displaystyle J _ {\ alpha} (z) = {\ frac {\ left ({\ frac {z} {2}} \ right) ^ {\ alpha}} {\ Gamma (\ alpha +1)}} \ prod _ { n = 1} ^ {\ infty} \ left (1 - {\ frac {z ^ {2}} {j _ {\ alpha, n} ^ {2}}} \ right)}{\ displaystyle J _ {\ alpha} (z) = {\ frac { \ left ({\ frac {z} {2}} \ right) ^ {\ alpha}} {\ Gamma (\ alpha +1)}} \ prod _ {n = 1} ^ {\ infty} \ left (1 - {\ frac {z ^ {2}} {j _ {\ alpha, n} ^ {2}}} \ right)}

(Есть большой число других известных интегралов и тождеств, которые здесь не воспроизводятся, но которые можно найти в справочных материалах.)

Повторяющиеся отношения

Функции J α, Y α, H. αи H. αудовлетворяют рекуррентным соотношениям

2 α x Z α (x) = Z α - 1 (x) + Z α + 1 ( х) {\ displaystyle {\ frac {2 \ alpha} {x}} Z _ {\ alpha} (x) = Z _ {\ alpha -1} (x) + Z _ {\ alpha +1} (x)}{\ displaystyle {\ frac {2 \ alpha} {x}} Z_ { \ альфа} (х) = Z _ {\ альфа -1} (х) + Z _ {\ альфа +1} (х)}

и

2 d Z α (x) dx = Z α - 1 (x) - Z α + 1 (x), {\ displaystyle 2 {\ frac {dZ _ {\ alpha} (x)} {dx} } = Z _ {\ alpha -1} (x) -Z _ {\ alpha +1} (x),}{\ displaystyle 2 {\ frac {dZ _ {\ alpha} (x)} {dx}} = Z _ {\ alpha -1} (x) -Z _ {\ альфа +1} (х),}

. где Z обозначает J, Y, H или H. Эти две идентичности часто комбинируются, например добавлены или вычтены, чтобы получить различные другие отношения. Таким образом, например, можно вычислить функции Бесселя более высоких порядков (или более высокие производные) с учетом значений более низких порядков (или более низких производных). В частности, следует, что

(1 xddx) m [x α Z α (x)] = x α - m Z α - m (x), (1 xddx) m [Z α (x) x α] = (- 1) м Z α + м (х) х α + м. {\ displaystyle {\ begin {align} \ left ({\ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {m} \ left [x ^ {\ alpha} Z_ { \ alpha} (x) \ right] = x ^ {\ alpha -m} Z _ {\ alpha -m} (x), \\\ left ({\ frac {1} {x}} {\ frac {d } {dx}} \ right) ^ {m} \ left [{\ frac {Z _ {\ alpha} (x)} {x ^ {\ alpha}}} \ right] = (- 1) ^ {m} {\ frac {Z _ {\ alpha + m} (x)} {x ^ {\ alpha + m}}}. \ end {align}}}{\ displaystyle {\ begin {align} \ left ({\ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {m} \ left [x ^ {\ alpha} Z _ {\ alpha} (x) \ right] = x ^ {\ alpha -m} Z _ {\ alpha -m} (x), \\\ left ({\ frac {1} {x}} {\ frac {d} {dx}} \ right) ^ {m} \ left [{\ frac {Z _ {\ alpha} (x)} {x ^ {\ alpha}}} \ right] = (- 1) ^ { m} {\ frac {Z _ {\ alpha + m} (x)} {x ^ {\ alpha + m}}}. \ end {align}}}

Модифицированные функции Бесселя подчиняются аналогичным отношениям:

e (x 2) (T + 1 T) знак равно ∑ N = - ∞ ∞ I N (x) tn {\ displaystyle e ^ {\ left ({\ frac {x} {2}} \ right) \ left (t + {\ frac {1} {t}} \ right)} = \ sum _ {n = - \ infty} ^ {\ infty} I_ {n} (x) t ^ {n}}{\ displaystyle e ^ {\ left ({\ frac {x} {2}} \ right) \ left (t + {\ frac {1} {t }} \ right)} = \ sum _ {n = - \ infty} ^ {\ infty} I_ {n} (x) t ^ {n}}

и

ez cos ⁡ θ знак равно I 0 (z) + 2 ∑ n знак равно 1 ∞ I n (z) cos ⁡ n θ. {\ displaystyle e ^ {z \ cos \ theta} = I_ {0} (z) +2 \ sum _ {n = 1} ^ {\ infty} I_ {n} (z) \ cos n \ theta.}{\ displaystyle e ^ {z \ cos \ theta} = I_ {0} (z) +2 \ sum _ {n = 1} ^ { \ infty} I_ {n} (z) \ cos n \ theta.}

и

1 2 π ∫ 0 2 π ez cos ⁡ (m θ) + y cos ⁡ θ = I 0 (z) I 0 (y) + 2 ∑ n = 1 ∞ I n (z) I mn (у). {\ displaystyle {\ frac {1} {2 \ pi}} \ int _ {0} ^ {2 \ pi} e ^ {z \ cos (m \ theta) + y \ cos \ theta} = I_ {0} (z) I_ {0} (y) +2 \ sum _ {n = 1} ^ {\ infty} I_ {n} (z) I_ {mn} (y).}{\ displaystyle {\ frac {1} {2 \ pi}} \ int _ {0} ^ {2 \ pi} e ^ {z \ cos (m \ theta) + y \ cos \ theta} = I_ {0} (z) I_ {0} (y) +2 \ sum _ {n = 1} ^ {\ infty} I_ {n} (z) I_ {mn} (y).}

. Рекуррентное соотношение имеет вид

C α - 1 (x) - C α + 1 (x) = 2 α x C α (x), C α - 1 (x) + C α + 1 (x) = 2 d C α dx, {\ displaystyle {\ begin {align} C _ {\ alpha -1} (x) -C _ {\ alpha +1} (x) = {\ frac {2 \ alpha} {x}} C _ {\ alpha} (x), \\ C _ {\ alpha -1} (x) + C _ {\ alpha +1} (x) = 2 {\ frac {dC _ {\ alpha}} {dx}}, \ end {align}}}{\ displaystyle {\ begin {align} C _ {\ alpha -1} ( x) -C _ {\ alpha +1} (x) = {\ frac {2 \ alpha} {x}} C _ {\ alpha} (x), \\ C _ {\ alpha -1} (x) + C_ {\ alpha +1} (x) = 2 {\ frac {dC _ {\ alpha}} {dx}}, \ end {align}}}

, где C α обозначает I α или eK α. Эти рекуррентные соотношения полезны для задач дискретной диффузии.

Теорема умножения

Функции Бесселя подчиняются теореме умножения

λ - ν J ν (λ z) = ∑ n = 0 ∞ 1 n! ((1 - λ 2) Z 2) N J ν + N (Z), {\ Displaystyle \ lambda ^ {- \ nu} J _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ {\ infty} {\ frac {1} {n!}} \ left ({\ frac {\ left (1- \ lambda ^ {2} \ right) z} {2}} \ right) ^ {n} J_ {\ nu + n} (z),}{\ displaystyle \ lambda ^ {- \ nu} J _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ {\ infty} {\ frac {1} {n!}} \ left ({\ гидроразрыва {\ left (1- \ lambda ^ {2} \ right) z} {2}} \ right) ^ {n} J _ {\ nu + n} (z),}

где λ и ν могут быть взяты как произвольные комплексные числа. Для | λ - 1 | < 1, the above expression also holds if J is replaced by Y. The analogous identities for modified Bessel functions and |λ − 1| < 1 are

λ - ν I ν (λ z) знак равно ∑ n = 0 ∞ 1 n! ((λ 2 - 1) z 2) N я ν + N (z) {\ Displaystyle \ lambda ^ {- \ nu} I _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ { \ infty} {\ frac {1} {n!}} \ left ({\ frac {\ left (\ lambda ^ {2} -1 \ right) z} {2}} \ right) ^ {n} I_ { \ nu + n} (z)}{\ displaystyle \ lambda ^ {- \ nu} I _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ {\ infty} {\ frac {1} {n !}} \ left ({\ frac {\ left (\ lambda ^ {2} -1 \ right) z} {2}} \ right) ^ {n} I _ {\ nu + n} (z)}

и

λ - ν K ν (λ z) = ∑ n = 0 ∞ (- 1) nn! ((λ 2 - 1) z 2) n K ν + n (z). {\ displaystyle \ lambda ^ {- \ nu} K _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ {\ infty} {\ frac {(-1) ^ {n}} {n !}} \ left ({\ frac {\ left (\ lambda ^ {2} -1 \ right) z} {2}} \ right) ^ {n} K _ {\ nu + n} (z).}{\ displaystyle \ lambda ^ {- \ nu} K _ {\ nu} (\ lambda z) = \ sum _ {n = 0} ^ {\ infty} {\ frac {(-1) ^ {n}} {n!}} \ left ({\ frac {\ left (\ lambda ^ {2} -1 \ right) z} {2}} \ right) ^ {n} K _ {\ nu + n} (z).}
Нули функции Бесселя

Гипотеза Бурже

Сам Бессель первоначально доказал, что для неотрицательных целых чисел n уравнение J n (x) = 0 имеет бесконечное число решений по x. Однако, когда функции J n (x) нанесены на один и тот же график, кажется, что ни один из нулей не совпадает для разных значений n, за исключением нуля при x = 0. Это явление известно как Гипотеза Бурже в честь французского математика XIX века, изучавшего функции Бесселя. В частности, он заявляет, что для любых целых чисел n ≥ 0 и m ≥ 1 функции J n (x) и J n + m (x) не имеют общих нулей, кроме одного при x = 0. Гипотеза была доказана Карлом Людвигом Сигелем в 1929 году.

Численные подходы

Для численных исследований нулей функции Бесселя см. Gil, Segura Temme (2007) ошибка harvtxt: нет цели: CITEREFGilSeguraTemme2007 (help ), Kravanja et al. (1998) ошибка harvtxt: нет цели: CITEREFKravanjaRagosVrahatisZafiropoulos1998 (help ) и Moler (2004) ошибка harvtxt: нет цели: CITEREFMoler2004 (help ).

См. Также
Примечания
Ссылки
Внешние ссылки
Последняя правка сделана 2021-05-12 14:02:15
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте