Скачиваний:
0
Добавлен:
17.02.2026
Размер:
713.98 Кб
Скачать

Рис. 1

41

1

J 2( )

 

xJ2( x)dx

 

 

.

(8.2)

 

2

0

 

 

 

 

Равенство (8.1) дает возможность находить коэффициенты особого разложения какой-нибудь функции f (x) по функциям Бесселя, которое имеет вид

 

 

f (x) amJ ( mx).

(8.3)

m 1

 

Здесь 0 1 2 3 4 — положительные нули функ-

ции J (x), занумерованные в порядке их возрастания. Предположим, что разложение (8.3) существует и что ряд в

правой части этого равенства сходится на отрезке 0,1 равно-

мерно. Тогда помножим обе части равенства (8.3) на xJ ( k x) и

в правой части проведем почленное интегрирование. Таким образом, получается соотношение

1

 

1

x f (x)J ( k

x)dx am xJ ( m x)J ( k x)dx . (8.4)

0

m 1

0

В силу равенства (8.1) в правой части (8.4) равны нулю все

те слагаемые, где m k . Поэтому равенство (8.4) перепишется в таком виде:

1

 

 

 

 

 

 

1

J 2

(

 

)

 

x f (x)J ( kx)dx

ak xJ2( k x)dx ak

k

 

 

 

 

.

 

2

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

Отсюда получается следующее выражение:

 

 

 

 

 

 

 

 

 

2

 

 

1

 

 

 

 

 

ak

 

 

 

 

 

x f (x)J ( k x)dx.

 

 

 

 

(8.5)

J

2

(

 

)

 

 

 

 

 

 

 

k

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разложение (8.3), в котором коэффициенты am определяют-

ся с помощью формулы (8.5), называется разложением функции f (x) в ряд Фурье–Бесселя. Подобные разложения встречаются при исследовании многих вопросов.

42

Предыдущие рассуждения были основаны на двух гипотезах: предполагалось, что разложение такого вида для данной функции существует и что оно равномерно сходящееся, или, по крайней мере, что почленное интегрирование ряда, какое требовалось по ходу дела, законно. Для произвольной функции f (x) справедливость этих гипотез непосредственно не видна. Поэтому предыдущих рассуждений недостаточно, чтобы установить справедливость разложения функций в ряд Фурье–Бесселя.

Ниже без доказательства приводятся теоремы о разложении функций в ряды Фурье по функциям Бесселя. Эти теоремы уточняют общую теорему Стеклова о разложении функции в ряд Фурье по собственным функциям для частного случая, когда собственными функциями являются функции Бесселя.

Теорема 8.1. Пусть функция f (x) непрерывна на отрезке [0,l] и имеет абсолютно интегрируемую на (0,l) производную.

 

 

 

k

 

Тогда ряд Фурье этой функции по функциям Бесселя J

 

 

x

 

 

 

 

l

 

 

 

 

( 0) сходится равномерно к f (x) на всяком отрезке [ , l ]

где 0 l . Если f (l) 0, то сходимость будет равномерной

2

на всяком отрезке [ ,l]. (Здесь k — положительные корни

уравнения J (x) 0).

Замечание. Утверждение теоремы 8.1 справедливо также в случае, когда k — положительные корни уравнения (5.7)

J (x) xJ (x) 0 ,

если дополнительно потребовать, чтобы и опустить ус-

 

 

 

ловие f (l) 0.

 

 

Теорема 8.2. Если

f (x) непрерывна и дважды дифферен-

цируема на отрезке [0,l]

и f (0) f (0) 0,

f (l) l f (l) 0 ,

43

 

 

 

 

 

 

k

 

 

 

то ее ряд Фурье по функциям Бесселя J

 

 

x

порядка

0

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

сходится равномерно к f (x) на

отрезке

[0,l]. (Здесь

k

положительные корни уравнения

J

 

(x) xJ (x) 0 ).

 

 

 

 

 

 

 

 

 

 

§ 9. Применение функций Бесселя при решении задач

Многие задачи, возникающие в различных областях науки и техники, математически формулируются как задачи нахождения функции, удовлетворяющей некоторому заданному дифференциальному уравнению с частными производными и дополнительным краевым условиям. Одним из методов решения таких задач является метод разделения переменных. Основная идея этого метода состоит в поиске решения задачи в виде суммы (конечной или бесконечной) специальных частных решений дифференциального уравнения. Всякое из частных решений должно иметь специальную структуру: оно должно представлять собой произведение функций, каждая из которых зависит только от одной независимой переменной. Для определения функций одной переменной при этом возникают обыкновенные дифференциальные уравнения. Если при формулировке задачи используются полярные или цилиндрические координаты, то одним из упомянутых обыкновенных дифференциальных уравнений оказывается уравнение Бесселя (1.1).

Рассмотрим задачу о малых поперечных колебаниях однородной круглой мембраны радиуса R, закрепленной по краям, под действием начального возмущения. Математически эта задача формулируется так.

Пусть D — круг радиуса R, т. е. D x (x1,x2): x R . Требуется найти функцию

u(t,x) C2 D (0, ) C1

 

[0, ) ,

u(t,k)

,

D

удовлетворяющую уравнению

44

 

 

u

a2

x

u ,

 

 

(x,t) D (0, ) ,

(9.1)

 

 

tt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

начальным условиям

 

 

 

 

 

 

 

 

 

 

 

u(0,x) u0(x) ,

ut (0,x) u1(x)

x

 

,

(9.2)

 

 

D

и однородным граничным условиям первого рода

 

 

 

u(t,x)

 

D 0 ,

 

 

 

 

 

t 0.

(9.3)

 

 

 

 

 

 

 

 

Здесь u

 

(x) C1

 

, u (x) C

 

, и выполняются условия согла-

0

D

D

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

сования

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0(x)

 

D 0,

 

u1(x)

 

D 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если перейти к полярным координатам (r, )

 

 

 

 

x (rcos , rsin ),

 

r 0,

0 2 ,

 

и к новой функции w(t, r, ) u[t, x(r, )], то смешанная задача (9.1)–(9.3) запишется следующим образом:

 

a2

 

 

 

1

 

1

 

 

 

 

w

w

 

 

w

 

w

 

,

(9.4)

 

r2

tt

 

 

rr

 

r r

 

 

 

 

 

0 r R ,

0 2 ,

t 0 ,

 

w(0, r, ) w0(r, ),

wt(0, r, ) w1(r, ),

(9.5)

 

0 r R,

0 2 ,

 

 

w(t, R, ) 0,

 

 

(9.6)

 

0 2 ,

t 0,

 

 

и при 0 r R, 0 2 ,

t 0 должно выполняться

 

 

w(t, r, )

 

,

w0(R, ) 0,

w1(R, ) 0,

 

 

 

 

 

 

 

w(t, r, 2 ) w(t, r, ).

 

Условие ограниченности w является естественным

следствием физической постановки задачи, а условие периодичности по является следствием условия однозначности решения.

Для решения задачи (9.4)–(9.6) применим метод разделения переменных. Ищем решение уравнения (9.4) в виде произведения двух функций

w(t, r, ) T(t) V(r, ) .

45

Подставив это представление в уравнение (9.4), получим

 

T V a2 T

 

1

 

 

 

1

 

 

 

V

 

 

 

V

 

V

.

r

 

 

 

 

 

 

 

rr

 

 

 

r

 

r2

 

Разделяя переменные, будем иметь

 

 

 

 

 

 

 

T

1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

V

 

 

V

 

 

V

 

.

(9.7)

 

 

 

 

r2

 

 

a2T V

rr

r r

 

 

 

 

 

 

 

В равенстве (9.7) левая часть зависит

 

только от t , а правая

часть — только от

r, . Следовательно, общая величина этих

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

T(t) и V(r, ):

T (t) a2T(t) 0,

 

 

 

t 0 ,

 

 

 

(9.8)

V

 

 

1

V

 

1

V r 0 ,

 

0 r R,

 

 

0 2 .

(9.9)

 

 

r2

 

 

 

rr

 

r r

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения (9.8) находится без труда.

 

Будем искать такие решения уравнения (9.9), которые явля-

ются 2 -периодические по

(т. е. V(r, ) V(r, 2 )), кото-

рые ограничены в круге (т. е.

 

V(r, )

 

) и которые удовлетво-

 

 

ряют краевому условию V(R, ) 0.

 

Эти решения также будем

искать методом разделения переменных, а именно положим

 

 

 

 

 

 

 

 

V(r, ) ( ) Z(r).

 

 

 

Подстановка этого представления в (9.9) дает

 

 

 

 

 

 

 

 

( ) 2 ( ) 0,

 

 

 

(9.10)

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.11)

 

 

 

 

 

 

 

Z (r)

 

r

 

Z (r)

 

 

r

2

Z(r) 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие 2 -периодичности функции

 

V(r, ) приводит к

тому,

что функция ( ) также должна быть

2 -периодической.

Это возможно лишь при n – целом. Поэтому общий интеграл уравнения (9.10) есть

~

~

n – целое.

( ) An cosn Bn sinn ,

46

Уравнение (9.11) имеет теперь вид

 

 

1

 

 

 

n

2

 

 

Z

 

 

 

 

 

 

(9.11’)

(r)

 

Z (r)

 

r

2

Z(r) 0.

 

 

r

 

 

 

 

 

Введем новую независимую переменную x r и новую функцию y(x) y( r) Z(r). Непосредственной подстановкой убеждаемся в том, что функция y(x) удовлетворяет уравнению Бесселя:

 

1

 

 

 

n

2

 

 

 

 

 

 

 

 

(9.12)

y (x)

 

y (x)

1

 

r

2

y(x) 0.

 

x

 

 

 

 

 

Общее решение уравнения Бесселя (9.12) имеет вид

 

y(x) c1Jn(x) c2Yn(x).

 

 

Условие ограниченности функции

V(r, ) требует,

чтобы

решения уравнения Бесселя также выбирались ограниченными.

Поэтому y(x) c1Jn(x). Отсюда Z(r) c1Jn( r).

Краевое условие V(R, ) 0 приводит к характеристическому уравнению для определения числа :

V(R, ) 0,

Z(R) 0,

Jn( R) 0.

Следовательно, R (kn) , (k 1, 2, ) – k -й положительный

корень функции Бесселя первого рода целого индекса n. Для числа получаем

(n)(n) k k R

2

, (n 0,1, 2, ; k 1, 2, ).

 

 

 

 

(n)

 

 

Таким образом, функции Z

nk

(r) J

 

k

r

являются ограни-

R

 

 

n

 

 

 

 

 

 

 

 

 

ченными решениями уравнения (9.11), обращающимися в нуль при r R.

47

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

 

 

 

 

 

a (n)

 

 

 

 

a (n)

 

 

 

 

 

(n)

 

 

w (t, r, ) A

 

 

cos

 

k

t

 

B

sin

 

 

k

t cosn J

 

 

k

 

r

 

 

 

 

R

 

 

R

 

R

nk

 

 

nk

 

 

 

 

 

nk

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a (n)

 

 

 

 

 

a (n)

 

 

 

 

(n)

 

 

 

 

C

 

cos

 

 

k

 

t D

 

 

sin

k

t

sinn J

 

k

 

 

r

,

 

 

 

 

R

 

 

 

R

R

 

 

 

 

 

nk

 

 

 

 

nk

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n 0,1, 2, ;

k 1, 2, )

 

 

 

 

 

 

 

 

 

являются ограниченными решениями уравнения (9.4), удовлетворяющие условиям гладкости и граничному условию.

Решение задачи (9.4)–(9.6) ищется в виде формального ряда

w(t, r, ) wnk (t, r, ).

 

 

 

 

 

 

n 0 k 1

 

 

 

Ank , Bnk , Cnk , Dnk оп-

Неизвестные коэффициенты разложения

ределяются в результате разложения функций w0(r, )

и w1(r, )

по системе функций

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n)

 

 

 

 

 

(n)

 

 

 

 

 

 

k

 

 

 

 

k

 

 

 

 

J

 

 

r

cosn , J

 

 

r

sinn .

 

R

R

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

48

Список литературы

1.Кузьмин Р.О. Бесселевы функции. Л.-М.: ГРОЛ, 1935.

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

3.Коренев Б.Г. Введение в теорию бесселевых функций. М.: Нау-

ка, 1971.

3.Владимиров В.С. Уравнения математической физики. М.: Нау-

ка, 1981.

4.Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1977.

5.Уроев В.М. Уравнения математической физики. М.: ИФ «Яу-

за», 1998.

49

СОДЕРЖАНИЕ

 

Введение . . . . . . . . . . . . . . . . .

3

§ 1.

Дифференциальное уравнение Бесселя . . . .

4

§ 2.

Решение уравнения Бесселя . . . . . . . . .

5

§ 3. Линейные зависимости между функциями

17

 

Бесселя . . . . . . . . . . . . . . . . . .

 

§ 4.

Функции Бесселя, индекс которых равен целому

21

 

числу с половиной . . . . . . . . . . . . .

 

§ 5.

Ортогональность функций Бесселя . . . . . .

24

§ 6. Поведение функций Бесселя при больших значе-

 

 

ниях аргумента . . . . . . . . . . . . . .

29

§ 7.

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

30

§ 8.

О разложении функций в ряды Фурье–Бесселя .

40

§ 9. Применение функций Бесселя при решении задач

44

 

Список литературы . . . . . . . . . . . . .

49

50

Соседние файлы в папке Вспомогательные