Файл: Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 156

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Таким образом, при г -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.Одним из наиболее употребительных методов по­ лучения асимптотических представлений является метод перевала. Сущность этого метода состоит в том, что при

11*

323


324

б о л ь ш и х з н а ч е н и я х п е р е м е н н о й х в е л и ч и н а и н т е г р а л а

/(*) = № (£) ex^ d l

с

определяется главным образом тем участком Си контура

интегрирования С, на котором

| е*Ф (£>| = ех Re

велик по

сравнению со значениями этого

модуля на остальной части

контура С. При этом интеграл по участку Си оценивается тем легче, чем меньше этот участок и чем круче падает величина xRe(p(|). При применении метода перевала ста-

325

раются деформировать путь интегрирования С в наиболее

выгодный, в указанном выше смысле, контур С. По тео­ реме Коши такая деформация, если она не выводит за пределы области аналитичности функций ср (|) и 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,

то в качестве С надо

выбрать

линию

наиболее

крутого перевала.

 

 

 

 

326


О ц е н к а и н т е г р а л а

$

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.

 

327


§ 8 . Ф ун к ц и и Эйри

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

У"~ху = 0.

(94)

Произведем замену неизвестной функции по формулам

. . [ 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,

3 2 8


п о л у ч и м ф ун кц и и Э й р и 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. Исследовать процесс

остывания проводника после выключения тока, если в течение всего

процесса

на

поверхности

проводника

происходил

теплообмен по

закону Ньютона со средой нулевой температуры.

329