Таким образом, при г -v-j-oo справедливы следующие асимптотические представления:
Так как функция ez не имеет нулей, то асимптотиче ские представления для функций Ганкеля можно записать в виде
|
|
|
|
|
|
|
(90) |
|
|
|
|
|
|
|
(91) |
7. |
Из |
асимптотических |
представлений непосредственно |
следует теорема |
3 § 4, |
а также утверждение: расстояние |
между двумя |
соседними нулями функций Бесселя Jv (г) |
(а также |
функций Неймана) |
стремится к л с неограни |
ченным ростом абсолютных величин нулей. |
|
На рис. |
35 |
приводятся |
графики |
функций |
Бесселя, |
а на рис. 36 — графики функций Неймана. Если восполь |
зоваться |
формулами (58) и (59) (§ 6), |
то легко |
получить |
следующие асимптотические представления функций Iv (2) |
и Kv (г) |
при |
больших |
| z | и Jarg 2 1< |
зт |
|
у — б: |
|
|
|
|
|
|
|
|
(92) |
|
|
|
|
|
|
|
(93) |
На рис. 37 и 38 приведены графики функций Kv(x)
иIv(x).
8.Одним из наиболее употребительных методов по лучения асимптотических представлений является метод перевала. Сущность этого метода состоит в том, что при
б о л ь ш и х з н а ч е н и я х п е р е м е н н о й х в е л и ч и н а и н т е г р а л а
/(*) = № (£) ex^ d l
с
определяется главным образом тем участком Си контура
интегрирования С, на котором |
| е*Ф (£>| = ех Re |
велик по |
сравнению со значениями этого |
модуля на остальной части |
контура С. При этом интеграл по участку Си оценивается тем легче, чем меньше этот участок и чем круче падает величина xRe(p(|). При применении метода перевала ста-
раются деформировать путь интегрирования С в наиболее
выгодный, в указанном выше смысле, контур С. По тео реме Коши такая деформация, если она не выводит за пределы области аналитичности функций ср (|) и if (|) и области существования интеграла, не меняет значения интеграла. В силу аналитичности функции ф (|) = и (а, (3) + + ш(а, р), | —а + ф, направление наибыстрейшего изме
нения функции |
и (а, Р) |
совпадает с направлением линии |
у (a, Р) = |
const. |
Контур |
Си должен содержать точку | 0 = |
= а 0 + /р |
0, в которой и (а, Р) достигает наибольшего зна |
чения (среди значений этой функции на С).
Нетрудно показать, что ср' (|0) = 0. Действительно, про
изводная от и (а, Р) вдоль линии С, |
взятая в точке | 0, |
равна нулю, ди |
■0, так как в точке |„ функция и (а, |
Р) |
|
ds ; = lo |
|
|
|
|
|
достигает максимального значения (вдоль С). |
|
|
Далее, |
I16 = 60 = 0 , |
поскольку в окрестности точки | = |
| 0 |
имеем v = const (вдоль С). |
Поэтому |
|
|
|
|
|
ди |
, -dv |
|
|
|
|
ф' ( У ; ' ds |
+ ldl 6 = 1. |
|
|
Точка | 0 для поверхности и = и(а, Р) |
является, |
очевидно, |
точкой перевала (седловой точкой). |
метода |
перевала |
Таким |
образом, |
при |
применении |
к асимптотической оценке |
интеграла |
ф (|) e'vtp |
d\ путь |
с
интегрирования С надо деформировать в путь С, прохо дящий через точку | 0, в которой ф' (|0) = 0, и в окрест ности этой точки совпадающий с линией у (а, Р) = const =
= v (а0, Ро) *).
*) Заметим, |
что окрестность точки перевала f0 |
разбивается |
линией уровня и (а, (5) = и ( а 0, ро) на 2п секторов ( п ^ 2 , |
где (п— 1) — |
кратность |
нуля |
функции <р' (|) в точке | 0), над которыми |
поверх |
ность и = |
и { а , Р) |
находится попеременно то выше, то |
ниже |
своей |
касательной плоскости в точке (а0, р0, и0). Линия v ( a , |
р) = |
о ( |
а 0, р0) |
вокрестности точки £0 состоит из п линий, проходящих через точку | 0
внаправлении биссектрис упомянутых секторов. Одну из таких линий
и следует взять |
в |
качестве С. Если и = |
и ( а , Р) |
имеет |
несколько |
точек перевала |
| 0, |
то в качестве С надо |
выбрать |
линию |
наиболее |
крутого перевала. |
|
|
|
|
О ц е н к а и н т е г р а л а
$ |
1p(l)exf(^ d l |
|
С \ \ |
|
|
производится следующим образом. |
их приближенными |
Функции г|з (|) и / (|) |
заменяются |
значениями |
|
|
(1) ~ do) + ( £ - go)V |
do), / (gW |
(go) + ^ 2— Г (go), |
а интеграл по контуру Си заменяется интегралом *)
S ^ (So) exf{l°] ^ (|- Ы2Г(Ы =
Сц
■ /у \ v f/ t \ Г
= № ) е * м ы \ cn
т (1 — i o ) 2 /"" (lo ) |
« . |
e2 |
d|. |
Последний интеграл соответствующей заменой переменной интегрирования приводится к интегралу
о(х)
Т; - оIй е~ ° ' ^
вкотором а(х)->оо при х->оо. Он равен
Следовательно,
I *®е"* Л/ I * ®)^ ,ы+о (ггЬо)'
сп
Здесь указана лишь идея получения асимптотических
представлений интеграла $ г|з (£) е*А£) dl при х - ^ оо . Под-
с
робное обоснование метода перевала смотрите, например,
в книге А. Г. Свешникова и А. Н. |
Тихонова **). |
*) При этом пренебрегаем слагаемым (| —| 0) 'К (5о)> малым в срав |
нении с 1)3 (£о). |
А. Н., Теория функций |
**) С в е ш н и к о в А. Г., Т и х о н о в |
комплексной переменной, «Наука», 1967. |
|
§ 8 . Ф ун к ц и и Эйри
Ряд задач физики (например, задача о движении заря женной частицы в однородном электрическом поле и др.) приводит к уравнению
Произведем замену неизвестной функции по формулам
. . [ V xz (х) |
для |
х ^ О , |
у{х) = 1 г___ |
для |
х < 0 . |
\ \ r— xz(x) |
Для функции z (х) получим уравнение
z" + -- z' (х) ~ |
+ х ) г = 0. |
(95) |
Для построения общего решения этого уравнения про изведем в нем замену независимой переменной по фор мулам
|
|
v х3/2 |
|
|
|
для |
х 5= 0 , |
Ч |
2 |
|
|
2 |
для |
Х <0 . |
|
|
3“ (— х)з/'2 |
= -3 j * |3/-2 |
При этом уравнение (95) перейдет в уравнения |
d4 |
+ t |
dz |
1 / 9 |
+ 1 2 = 0 |
ДЛЯ |
x ^ O , |
dP |
dt |
|
|
d4 |
, 1 |
dz . |
|
|
1 / 9 |
|
x < 0. |
dP |
T ’ dt + |
|
|
P 2 = 0 |
Д Л Я |
Это — уравнения цилиндрических |
функций. |
Их общие |
решения можно записать в следующем виде: |
|
г(/) = С1/ _ 1/;,(0 + |
С2/ 1/з(0 |
для |
х 2 г 0 , |
2 ( 0 = D1J _ 1/3(/) + |
D2J 1/3 (if) |
для |
х < |
0 . |
Следовательно, общее решение уравнения (94) можно запи сать в виде
V x |
C J—1/з ( з х 3/'2 ) + C.J j /з |
для х ^ |
0 , |
У(х) = |
2 |
для х < |
0 . |
К|Х |
DiJ—]/з(“з’ I х |3/ij + D2J |
Если произвольные постоянные Си С2, |
Du D-2 взять рав- |
ными |
|
|
|
Сх = — С2 = Dx= Ds = 1/3 и Ct =C2 = D1= — D2 = 1/3,
п о л у ч и м ф ун кц и и Э й р и A i (х ) и B i (х):
I - 1/3 ( з' *3/'2) - |
Л/з ( |
з *3/г)] |
Для |
* Ss О, |
Л/ (л-) = • |
|
|
|
|
|
|
|
/ И |
|
|
|
|
|
для х < 0 ; |
^-м /з(|Ч *13/2)+ Л /з (|Ч * 1 3/2)] |
' j/jg |
X3/2 |
Л, • |
|
|
|
для х ^ О , |
3 |
|
м |
|
) + |
|
|
|
|
Bi (х) : |
|
|
|
|
|
|
|
VV a |
2 1X |
- |
Ji^ \ |
х |3/'2) |
для |
х < 0 . |
|
|
|
зг 1 |
|
|
Из представления |
функций Iv (t) |
и Jv (t) в виде |
обоб |
щенных степенных рядов следует, что |
|
|
Лг(х)-> гт~— ----- , |
хАо |
|
— ------. |
|
>^18 Г (2/3) |
|
у 18 Г (2/3) |
|
|
Применяя приведенные в § 7 рассуждения (с очевидными несущественными изменениями) к функциям /._v (/) qp 7V(/) и J_v(t)± Jv (t), нетрудно получить следующие асимпто тические представления функций Эйри:
Ai1(х) - |
~ |
х |
~ |
I/* е * |
Х' Н [ 1 + |
О (х- »/*)] |
для |
х-> + оо |
ЛЦлг) |
у |
я |
х |
Щ 5ш ( з | д г ^ + |
Т ) х |
|
|
|
|
|
|
|
X [ 1 + |
О (1 х |— 1/*)] |
для |
X —> — оо |
Bi (х) = |
T/ Т |
|
— д;3/'2 |
[1-1-0 (х~з/*)] для |
х~> + оо, |
^т= х ~ 1/4 е3 |
|
у Я |
|
|
|
|
|
|
Bi (х) = |
|
|
|
sin ( I |
| * Г + ~ ) х |
|
|
|
|
|
|
|
X [ \+ |
0 ( \ х \ - '11)] |
для |
х->— оо. |
Имеются таблицы функций Эйри. |
|
|
ЗА Д А Ч И
1.Найти температуру бесконечного круглого цилиндра, начальная
температура |
которого равна |
и (г, 0) = А |
1 |
а |
на поверхности |
его поддерживается температура, равная нулю. |
|
радиуса R длительное |
2. |
Цилиндрический |
однородный |
проводник |
время нагревался постоянным током силы 1. Исследовать процесс |
остывания проводника после выключения тока, если в течение всего |
процесса |
на |
поверхности |
проводника |
происходил |
теплообмен по |
закону Ньютона со средой нулевой температуры.