Курсовая работа: Беселеві функції

Курсова робота

"Беселеві функції"


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