Беселеві функції
Курсова робота
"Беселеві функції"
1. Беселеві функції з будь-яким
індексом
Рівняння Лапласа в циліндричних
координатах
Щоб пояснити походження Беселевих
функцій, розглянемо рівняння Лапласа в просторі:
. (1)
Якщо перейти до циліндричних
координат по формулах:
, ,
,
те рівняння (1) прикмет наступний
вид:
. (2)
:
,
Нехай є рішення згаданого виду. Підставляючи його в (2),
одержимо:
,
звідки (після ділення на )
.
Записавши це у вигляді:
,
знайдемо, що ліва частина не
залежить від , права не
залежить від , ; отже, загальна величина цих
виражень є деяка постійна .
Звідси:
; ;
; ;
.
В останній рівності ліва частина
не залежить від , права не
залежить від ; отже,
загальна величина цих виражень є деяка постійна . Звідси:
, ;
, .
Таким чином, , , повинні
задовольняти лінійним диференціальним рівнянням другого порядку:
,
(3)
, ,
з яких друге й третє є найпростіші
лінійні рівняння з постійними коефіцієнтами, а перше є лінійним рівнянням зі
змінними коефіцієнтами нового виду.
Обернено, якщо , , задовольняють
рівнянням (3), тобто рішення
рівняння (2). Справді, підставляючи в ліву частину (2) і ділячи потім на , одержимо:
.
Таким чином, загальний вид всіх
трьох рішень рівняння (2), які є добутком трьох функцій, кожна з яких залежить
від одного аргументу, є , де
, , –
будь-які рішення рівнянь (3) при будь-якому виборі чисел , .
Перше з рівнянь (3) у випадку , називається рівнянням Беселя. Думаючи в цьому
випадку , позначаючи
незалежну змінну буквою (замість
), а невідому функцію –
буквою (замість ), знайдемо, що рівняння Беселя
має вигляд:
. (4)
Це лінійне диференціальне рівняння
другого порядку зі змінними коефіцієнтами відіграє більшу роль у додатках
математики. Функції, йому задовольняючі, називаються Беселевими, або
циліндричними, функціями.
Беселеві функції першого роду
Будемо шукати рішення рівняння
Беселя (4) у вигляді ряду:
.
Тоді
,
,
,
.
Отже, приходимо до вимоги
або до нескінченної системи
рівнянь
,
яка розпадається на дві системи:
Перша з них задовольниться, якщо
взяти … У другій системі можна взяти довільно; тоді … однозначно визначаються (якщо не є цілим негативним числом).
Взявши
,
знайдемо послідовно:
,
,
,
і як рішення рівняння (4) одержимо
ряд:
Цей ряд, що формально задовольняє
рівнянню (4), сходиться для всіх позитивних значень і, отже, є рішенням рівняння (4) в області (у випадку цілого в області ).
Функція
(5)
називається бесселевой функцією
першого роду з індексом .
Вона є одним з рішень рівняння Беселя (4). У випадку цілого ненегативного
індексу одержимо:
, (5`)
і, зокрема,
. (5``)
Загальне рішення рівняння Беселя
У випадку нецілого індексу функції і є
рішеннями рівняння (4). Ці рішення лінійно незалежні, тому що початкові члени
рядів, що зображують ці функції, мають коефіцієнти, відмінні від нуля, і
містять різні ступені .
Таким чином, у випадку нецілого індексу загальне рішення рівняння Беселя є:
. (6)
Якщо (ціле негативне число), то функція, обумовлена
формулою (5) (з огляду на, що дорівнює
нулю для …), приймає вид:
(5```)
або, після заміни індексу
підсумовування на ,
, (7)
звідки видно, що задовольняє разом з рівнянню Беселя
.
Але формула (6) у випадку цілого вже не дає загального рішення
рівняння (4).
Думаючи
(
– не ціле) (8)
і доповнюючи це визначення для (ціле число) формулою:
, (8`)
одержимо функцію , що задовольняє рівнянню Беселя (4) і у
всіх випадках лінійно незалежну від (у випадку , де – ціле). Функція називається беселевою функцією другого роду з
індексом . Загальне рішення
рівняння Беселя (4) можна записати у всіх випадках у вигляді:
. (9)
2. Формули приведення для Беселевих
функцій
Маємо:
; ;
, ;
.
Отже,
. (10)
Таким чином, операція (що складається в
диференціюванні з наступним множенням на ), застосована до , підвищує в цьому вираженні індекс на одиницю й міняє знак. Застосовуючи цю
операцію раз, де – будь-яке натуральне число,
одержуємо:
. (10`)
Маємо:
;
Отже,
. (11)
Таким чином, операція , застосована до , знижує в цьому вираженні індекс на одиницю. Застосовуючи цю
операцію раз, одержуємо:
. (11`)
З виведених формул можна одержати
деякі наслідки. Використовуючи (10), одержимо:
; ;
.
Звідси, зокрема, треба, що . Використовуючи (11), одержимо:
; ;
.
По членне додавання й вирахування
отриманих рівностей дає:
, (12)
. (13)
Формула (13) дозволяє виразити всі
Беселеві функції із цілими індексами через , .
Дійсно, з (13) знаходимо (думаючи ):
, (13`)
звідки послідовно одержуємо:
,
, …………………
3. Беселеві функції з
напівцілим індексом
Беселеві функції, загалом кажучи,
є новими трансцендентними функціями, що не виражаються через елементарні
функції. Виключення становлять Беселеві функції з індексом , де – ціле. Ці функції можуть бути виражені через
елементарні функції.
Маємо:
,
,
отже,
.
Але , значить:
. (14)
Далі
,
,
отже,
.
Але , тому
. (15)
За допомогою (10') знаходимо:
,
а з огляду на (14)
,
отже, при цілому позитивному
. (14`)
За допомогою (11') знаходимо:
,
але в силу (15)
,
і, отже, при цілому позитивному
. (15`)
4. Інтегральне подання Беселевих
функцій із цілим індексом
Виробляюча функція системи функцій
Розглянемо систему функцій (з будь-якою загальною областю визначення),
пронумерованих за допомогою всіх цілих чисел:
Складемо ряд
,
де – комплексна змінна. Припустимо, що при кожному (приналежному області визначення
розглянутих функцій) цей ряд має кільце збіжності, що містить усередині себе
одиничну окружність .
Зокрема, це кільце може являти собою повну площину комплексної змінної без
крапок 0 і?.
Функція
(16)
(де x лежить в області визначення
функцій системи , – усередині кільця збіжності, що
відповідає розглянутому значенню ) називається виробляючою функцією системи .
Обернено, нехай задана функція , де пробігає деяку множину, перебуває усередині деякого кільця, що залежить від
, із центром 0 і утримуючого
усередині себе одиничну окружність. Тоді, якщо при кожному аналітичне відносно усередині відповідного кільця, тобто виробляюча функція деякої системи функцій. Справді, розклавши при
кожному функцію в ряд Лорана по ступенях :
,
знайдемо, що система коефіцієнтів цього ряду буде шуканою системою
.
Формули для коефіцієнтів ряду
Лорана дозволяють виразити функції розглянутої системи через виробляючу функцію.
Застосовуючи ці формули й перетворюючи потім інтеграл уздовж одиничної
окружності в простий
інтеграл, одержимо:
. (17)
Виробляюча функція системи Беселевих
функцій із цілими індексами
Покажемо, що для системи Беселевих
функцій першого роду із цілими індексами (…)
виробляюча функція є:
.
Маємо:
, ,
звідки після по членного
перемножування цих рівностей знайдемо:
(тому що в передостанній
внутрішній сумі й були зв'язані залежністю , то ми могли покласти , одержавши підсумовування по
одному індексі ). В останній
внутрішній сумі підсумовування виробляється по всіх цілих , для яких , отже, при це буде ; при це буде . Таким чином, у всіх випадках внутрішня сума є в силу формул (5`) і (5```).
Отже,
, (18)
але це й доводить, що є виробляюча функція для системи
.
Виведемо деякі наслідки з формули
(18). Думаючи в ній ,
одержимо:
,
звідки після поділу дійсної й
мнимої частини (з огляду на, що )
(18`)
(18``)
Заміняючи в (18`) і (18``) на , знайдемо:
, (18```)
. (18````)
Інтегральне подання Jn(x)
Тому що, по доведеному, при маємо , те по формулі (17) одержуємо (використовуючи в
перетвореннях формули Ейлера):
де прийнято в увагу, що є парна функція від є непарна функція від . Отже, доведено, що для
будь-якого цілого числа
. (19)
Формула (19) дає подання Беселевих
функцій із цілим індексом у вигляді певного інтеграла, що залежить від
параметра . Ця формула
називається інтегральним поданням Беселя для , права частина формули називається інтегралом
Беселя. Зокрема, при знайдемо:
. (19`)
5. Ряди Фур'є-Беселя
Розглянемо на якому-небудь
інтервалі (кінцевому або
нескінченному) два диференціальних рівняння
, , (20)
де й –
безперервні функції на .
Нехай і – ненульові рішення цих рівнянь. Множення
на й на й наступне вирахування дають
.
Нехай і належать
і , тоді після інтегрування в межах від до одержимо
. (21)
Якщо й –
сусідні нулі рішення , то
між і зберігає постійний знак, нехай, наприклад, на (, )
(у противному випадку варто замінити на ), тоді , (рівність
нулю виключено, тому що –
ненульове рішення диференціального рівняння другого порядку). Якщо на , то повинна, принаймні, раз звертатися в нуль між і , тому що інакше збереже постійний знак на (,).
Нехай, наприклад, на (,) (у противному випадку заміняємо на ), і тоді з (21) одержимо протиріччя, тому що ліва
частина ≤0, а права >0. У такий спосіб доведена теорема порівняння
Штурму: якщо P(x)<Q(x) на розглянутому інтервалі I і якщо y і z – ненульові
рішення рівнянь (20), те між кожними двома сусідніми нулями y(x) перебуває
принаймні один нуль z(x).
З теореми порівняння Штурму
випливають нижченаведені наслідки. Якщо на , то кожне ненульове рішення рівняння може мати на не більше одного нуля (це легко бачити,
якщо покласти й взяти ). Якщо на (де ), то для всяких двох сусідніх нулів і ()
кожного ненульового рішення рівняння маємо (це легко бачити, якщо покласти , взяти й помітити, що нулями будуть тільки числа виду , ціле).
Якщо на (де ), то для всяких двох сусідніх нулів кожного
ненульового рішення рівняння маємо
(це легко бачити, якщо
покласти й взяти ). Із сказаного випливає, що
якщо на , те для всяких двох сусідніх нулів і ()
кожного ненульового рішення рівняння маємо .
Викладене показує, що якщо безперервно на й перевищує деяке позитивне число поблизу
+∞, те кожне ненульове рішення рівняння має на нескінченно багато нулів. Якщо ще поблизу не звертається в нуль, то ці нулі утворять
нескінченну зростаючу послідовність , що має межею +∞, а якщо, крім того, , де , те .
Розглянемо рівняння Беселя
на інтервалі . Підстановка приводить до рівняння
.
Очевидно, і мають
ті самі нулі. Тому що , де – ціла функція, то не має нулів на при досить малому , і тому що при , те при кожному нулі на утворять нескінченну зростаючу послідовність
причому .
Якщо , то задовольнить рівнянню
на інтервалі (0, +∞).
Підстановка приводить до
рівняння
і, отже, задовольняє цьому рівнянню. Таким чином, при
будь-яких позитивних і маємо
, де ,
, де ,
звідки
,
отже,
, де . (22)
Нехай тепер . Розкладання по ступенях починається зі члена, що містить , розкладання по ступенях починається зі члена, що містить , тому що коефіцієнт при дорівнює нулю, що легко бачити, виходячи
з формули (5). Отже, з (22) при одержимо
,
тобто
, (23)
звідки видно, що якщо і є різними нулями функції , те
. (23`)
Цим доведено, що при система функцій
на інтервалі є ортогональної щодо ваги .
Переходячи до межі при в співвідношенні
і використовуючи правило Лопиталя,
одержимо при всякому
, (24)
отже, якщо є нулем функції , те
. (24`)
Таким чином, при кожному всякій безперервній функції на , що задовольняє вимозі
,
поставлений у відповідність ряд
Фур'є-Беселя
, (25)
коефіцієнти якого визначаються
формулами
. (25`)
Можна довести, що система функцій на , ортогональна щодо ваги , замкнута. Зокрема, якщо ряд Фур'є-Беселя (25)
рівномірно сходиться до його безперервної функції, що породжує.
Можна показати, що якщо й безперервна на й функція, то ряд Фур'є-Беселя цієї функції
сходиться до неї при .
6. Асимптотичне подання Беселевих
функцій із цілим індексом для більших значень аргументу
Нехай – позитивна функція й – яка-небудь функція для досить більших значень . Запис
при
означає, що найдуться такі числа й M, що при маємо .
Подібний запис уживається й в
інших аналогічних випадках. Наприклад, якщо – позитивна функція й – яка-небудь функція, визначені для досить малих
позитивних значень , то запис
при
означає, що найдуться такі числа й , що на .
Допоміжна лема
Якщо двічі безупинно диференцюєма на , то для функції
має місце асимптотичне подання
при .
Доведемо цю лему. Заміняючи на , одержимо:
.(26)
Розглянемо інтеграл, що фігурує в
правої частини формули (20). Заміняючи на ,
знайдемо:
,
але, замінивши на , одержимо:
.
Якщо позитивно, убуває й прагнути до нуля при , то й , а отже, і є при
, тому
при ,
звідки
при .
Отже, одержуємо асимптотичне
подання:
при . (27)
Розглянемо тепер інтеграл, що
фігурує в другому складати^ся правої частини формули (20). Маємо:
,
.
Очевидно, двічі безупинно на , але існують і ,
тому стає безупинно диференцуєма
на . Інтегрування вроздріб
дає:
,
де перший доданок правої частини є при , а інтеграл у другому мажорирується інтегралом, що
складається при нижній межі
,
який сходиться, тому що
при ;
отже, другий доданок є теж при .
Отже, маємо:
при . (28)
З (26), (27), (28) одержуємо
шукане асимптотичне подання:
при . (29)
Із цієї формули, переходячи до
сполучених величин, знайдемо ще:
при . (29')
Формули (29) і (29`) вірні й для
функцій .
Висновок асимптотичної формули для
Jn(x)
Заміняючи на , одержимо:
(з огляду на, що є парна функція від , а є непарна функція від ). Підстановка дає:
,
де є, мабуть, поліном n-й ступеня (поліном Чебишева),
тому що з формули Муавра видно, що є поліном n-й ступеня відносно . Але
і, заміняючи в першому із цих
інтегралів на , одержимо:
Тому що й на
мають похідні всіх
порядків, то до двох останніх інтегралів застосовні формули (29) і (29`), і ми
одержуємо:
;
але ; ,
отже,
.
Отже, маємо шукане асимптотичне
подання беселевої функції першого роду із цілим індексом для більших значень
аргументу:
при . (30)
Ця формула показує, що з точністю складається до
порядку, що, є загасаючою
гармонікою із хвилею постійної довжини й амплітудою, що убуває обернено
пропорційно квадратному кореню з абсциси.
Зокрема,
при ; (30`)
при . (30'')
Графіки цих функцій зображені ні
малюнках 1 і 2.
Розглянемо кілька прикладів
рішення рівняння Беселя.
1. Знайти рішення рівняння Беселя
при
,
задовольняючим початковим умовам
при , і .
Рішення.
На підставі формули (5') знаходимо
одне приватне рішення:
.
2. Знайти одне з рішень рівняння:
, .
Рішення.
Зробимо заміну
.
При одержимо:
.
При будемо шукати рішення у вигляді узагальненого
статечного ряду:
.
Рівняння на має вигляд ;
, ,
, , тому
,
, .
Рисунок 1 – Графік функції y=J0 (x)
Рисунок 2 – Графік функції y=J1 (x)
Висновок
Розглянуті усі рішення рівнянь,
які можуть бути представлені у вигляді добутку трьох функцій. Складені графіки
функцій.
Список літератури
1. Пискунов Н.С.
Диференціальне й інтегральне вирахування, навчальний посібник для вузів. – К.,
2003
2. Романовський П.
І. «Ряди Фур'є. Теорія поля. Аналітичні й спеціальні функції. Перетворення
Лапласа», навчальний посібник для вузів. – К., 2004
3. Самарський А.А.,
Гулін А.В. Чисельні методи. – К., 2003
4. Синіцин
О.К., Навроцкий А.А. Алгоритми обчислювальної математики. – К., 2003