Это не официальный сайт wikipedia.org 01.01.2023

Функции Бесселя — Википедия

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

(перенаправлено с «Дифференциальное уравнение Бесселя»)

Фу́нкции Бе́сселя в математике — семейство функций, являющихся каноническими решениями дифференциального уравнения Бесселя:

x 2 d 2 y d x 2 + x d y d x + ( x 2 α 2 ) y = 0 ,

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

График функций Бесселя первого рода

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

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

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

ПримененияПравить

Уравнение Бесселя возникает во время нахождения решений уравнения Лапласа и уравнения Гельмгольца в цилиндрических и сферических координатах. Поэтому функции Бесселя применяются при решении многих задач о распространении волн, статических потенциалах и т. п., например:

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

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

Функция Бесселя является обобщением функции синуса. Ее можно трактовать как колебание струны с переменной толщиной, переменным натяжением (или одновременно обоими условиями); колебаниями в среде с переменными свойствами; колебаниями дисковой мембраны и т. д.

ОпределенияПравить

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

Функции Бесселя первого родаПравить

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

J α ( x ) = m = 0 ( 1 ) m m ! Γ ( m + α + 1 ) ( x 2 ) 2 m + α .  

Здесь Γ ( z )   — это гамма-функция Эйлера, обобщение факториала на нецелые значения. График функции Бесселя похож на синусоиду, колебания которой затухают пропорционально 1 x  , хотя на самом деле нули функции расположены не периодично (однако расстояние между двумя последовательными нулями стремится к π   при x  )[1].

Ниже приведены графики J α ( x )   для α = 0 , 1 , 2  :

Если α   не является целым числом, функции J α ( x )   и J α ( x )   линейно независимы и, следовательно, являются решениями уравнения. Но если α   целое, то верно следующее соотношение:

J α ( x ) = ( 1 ) α J α ( x ) .  

Оно означает, что в этом случае функции линейно зависимы. Тогда вторым решением уравнения станет функция Бесселя второго рода (см. ниже).

Интегралы БесселяПравить

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

J α ( x ) = 1 π 0 π cos ( α τ x sin τ ) d τ .  

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

J α ( x ) = 1 2 π π π e i ( α τ x sin τ ) d τ .  

Для нахождения интегрального представления функции Бесселя в случае нецелых α   необходимо учесть, что имеется разрез вдоль оси абсцисс. Это вызвано тем, что подынтегральное выражение более не является 2 π  -периодическим. Таким образом, контур интегрирования разбивается на 3 участка: луч от   до 1  , где φ = π  , окружность единичного радиуса и луч от 1   до +   при φ = π  . Проделав несложные математические преобразования, можно получить следующее интегральное представление:

J α ( x ) = 1 2 π π π e i ( x s i n ( φ ) α φ ) d φ sin ( α π ) π 1 e 1 2 x ( r 1 r ) r α + 1 d r .  

Нетрудно убедиться, что при целых α   это выражение переходит в предыдущую формулу.

Функции НейманаПравить

Функции Неймана — решения Y α ( x )   уравнения Бесселя, бесконечные в точке x = 0  .

Эта функция связана с J α ( x )   следующим соотношением:

Y α ( x ) = J α ( x ) cos ( α π ) J α ( x ) sin ( α π ) ,  

где в случае целого α   берётся предел по α  , вычисляемый, например, с помощью правила Лопиталя.

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

y ( x ) = C 1 J α ( x ) + C 2 Y α ( x ) .  

Ниже приведён график Y α ( x )   для α = 0 , 1 , 2  :

В ряде книг функции Неймана обозначаются N α ( x )  .

Сферические функции БесселяПравить

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

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

x 2 d 2 y d x 2 + 2 x d y d x + ( x 2 n ( n + 1 ) ) y = 0.  

Два линейно-независимых решения называются сферическими функциями Бесселя jn и yn, и связаны с обычными функциями Бесселя Jn и Неймана Yn с помощью[2]

j n ( x ) = π 2 x J n + 1 2 ( x ) , y n ( x ) = π 2 x Y n + 1 2 ( x ) = ( 1 ) n + 1 π 2 x J n 1 2 ( x ) .  

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

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

j n ( x ) = ( x ) n ( 1 x d d x ) n sin x x , y n ( x ) = ( x ) n ( 1 x d d x ) n cos x x .  

Несколько первых сферических функций Бесселя[4]:

j 0 ( x ) = sin x x , j 1 ( x ) = sin x x 2 cos x x , j 2 ( x ) = ( 3 x 2 1 ) sin x x 3 cos x x 2 , j 3 ( x ) = ( 15 x 3 6 x ) sin x x ( 15 x 2 1 ) cos x x  

и Неймана[5]:

y 0 ( x ) = j 1 ( x ) = cos x x , y 1 ( x ) = j 2 ( x ) = cos x x 2 sin x x , y 2 ( x ) = j 3 ( x ) = ( 3 x 2 + 1 ) cos x x 3 sin x x 2 , y 3 ( x ) = j 4 ( x ) = ( 15 x 3 + 6 x ) cos x x ( 15 x 2 1 ) sin x x .  

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

Производящие функции сферических функций Бесселя[6]:

1 z cos ( z 2 2 z t ) = n = 0 t n n ! j n 1 ( z ) , 1 z sin ( z 2 2 z t ) = n = 0 t n n ! y n 1 ( z ) .  

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

В следующих формулах fn может быть заменено на jn, yn, h(1)
n, h(2)
n, где h(1)
n и h(2)
n — сферические функции Ханкеля, для n = 0, ±1, ±2, ...[7]:

( 1 z d d z ) m ( z n + 1 f n ( z ) ) = z n m + 1 f n m ( z ) , ( 1 z d d z ) m ( z n f n ( z ) ) = ( 1 ) m z n m f n + m ( z ) .  

СвойстваПравить

ОртогональностьПравить

Пусть μ 1 , μ 2   — нули функции Бесселя J α ( x )  . Тогда[1]:

0 1 x J α ( μ 1 x ) J α ( μ 2 x ) d x = { 0 ; μ 1 μ 2 1 2 ( J α ( μ 1 ) ) 2 ; μ 1 = μ 2  .

АсимптотикаПравить

Для функций Бесселя первого и второго рода известны асимптотические формулы. При малых аргументах ( 0 < x α + 1 )   и неотрицательных α   они выглядят так[8]:

J α ( x ) 1 Γ ( α + 1 ) ( x 2 ) α ,  
Y α ( x ) { 2 π [ ln ( x / 2 ) + γ ] ; α = 0 Γ ( α ) π ( 2 x ) α ; α > 0  ,

где γ   — постоянная Эйлера — Маскерони (0,5772…), а Γ   — гамма-функция Эйлера. Для больших аргументов ( x | α 2 1 / 4 |  ) формулы выглядят так:

J α ( x ) 2 π x cos ( x α π 2 π 4 ) ,  
Y α ( x ) 2 π x sin ( x α π 2 π 4 ) .  

Использование следующего члена асимптотического разложения позволяет значительно уточнить результат. Для функции Бесселя нулевого порядка он выглядит следующим образом:

J 0 2 π x cos ( x π 4 ) + 1 4 x 2 π x sin ( x π 4 ) .  

Гипергеометрический рядПравить

Функции Бесселя могут быть выражены через гипергеометрическую функцию:

J α ( z ) = ( z / 2 ) α Γ ( α + 1 ) 0 F 1 ( α + 1 ; z 2 / 4 ) .  

Таким образом, при целых α   функция Бесселя однозначная аналитическая, а при нецелых — многозначная аналитическая.

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

Существует представление для функций Бесселя первого рода и целого порядка через коэффициенты ряда Лорана функции определённого вида, а именно

e z 2 ( w 1 w ) = n = + J n ( z ) w n .  

СоотношенияПравить

Формула Якоби — Ангера и связанные с нейПравить

Получается из выражения для производящей функции при a = 1  , w = e i ϕ  [9]:

e i z sin ϕ = J 0 ( z ) + 2 n = 1 J 2 n ( z ) cos ( 2 n ϕ ) + 2 i n = 1 J 2 n 1 ( z ) sin ( 2 n 1 ) ϕ .  

При a = 1  , t = i e i ϕ  [9]:

e i z cos ϕ = J 0 ( z ) + 2 n = 1 i n J n ( z ) cos ( n ϕ ) .  

Рекуррентные соотношенияПравить

Для функций Бесселя существует ряд рекуррентных соотношений. Приведём здесь некоторые из них:

J α + 1 = α x J α J α ( x ) ;  
J α + 1 ( x ) + J α 1 ( x ) = 2 α x J α ( x ) ;  
J α + 1 ( x ) J α 1 ( x ) = 2 J α ( x )  [10].

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

Для любого целого n и комплексных z 1  , z 2   выполняется[11]

J n ( z 1 + z 2 ) = k = J k ( z 1 ) J n k ( z 2 ) .  

Интегральные выраженияПравить

Для любых a   и b   (в том числе комплексных) выполняется[12]

0 e a t J n ( b t ) d t = b n a 2 + b 2 ( a 2 + b 2 + a ) n .  

Частным случаем последней формулы является выражение

0 e a t J 0 ( b t ) d t = 1 a 2 + b 2 .  

См. такжеПравить

ПримечанияПравить

  1. 1 2 Зубов В. И. . Функции Бесселя. — М.: МФТИ, 2007. Архивная копия от 24 июня 2016 на Wayback Machine
  2. Abramowitz and Stegun, p. 437, 10.1.1 Архивная копия от 2 сентября 2006 на Wayback Machine.
  3. Abramowitz and Stegun, p. 439, 10.1.25, 10.1.26 Архивная копия от 21 декабря 2009 на Wayback Machine.
  4. Abramowitz and Stegun, p. 438, 10.1.11 Архивная копия от 30 апреля 2009 на Wayback Machine.
  5. Abramowitz and Stegun, p. 438, 10.1.12 Архивная копия от 30 апреля 2009 на Wayback Machine.
  6. Abramowitz and Stegun, p. 439, 10.1.39 Архивная копия от 21 декабря 2009 на Wayback Machine.
  7. Abramowitz and Stegun, p. 439, 10.1.23, 10.1.24 Архивная копия от 22 декабря 2019 на Wayback Machine.
  8. Arfken G. B., Hans J. W. . Mathematical Methods for Physicists. 6th ed. — San Diego: Harcourt, 2005. — ISBN 0-12-059876-0.
  9. 1 2 Бейтмен, Эрдейи, 1974, с. 15.
  10. В. С. Гаврилов и др. Функции Бесселя в задачах математической физики Архивная копия от 26 ноября 2019 на Wayback Machine, стр. 7
  11. Лаврентьев, Шабат, 1973, с. 670.
  12. Лаврентьев, Шабат, 1973, с. 671.

ЛитератураПравить

  • Ватсон Г. . Теория бесселевых функций. — М.: ИЛ, 1949.
  • Бейтмен Г., Эрдейи А. . Функции Бесселя, функции параболического цилиндра, ортогональные многочлены // Высшие трансцендентные функции. Т. 2. 2-е изд / Пер. с англ. Н. Я. Виленкина. — М.: Наука, 1974. — 296 с.
  • Лаврентьев М. А., Шабат Б. В. . Методы теории функций комплексного переменного. — М.: Наука, 1973. — 736 с.