Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
имеет вид
y(x) = C1Jν(x) + C2Yν(x),
где C1 и C2 — произвольные постоянные.
На рис. 4.1 и 4.2 представлены графики функций Бесселя J0(x), J1(x) и Неймана Y0(x), Y1(x) для действительного аргумента. Следует отметить, что функции Неймана при x → 0+ не ограничены.
Рис. 4.1
Рис. 4.2
3. Для функций Бесселя (Неймана) справедливы следующие
формулы дифференцирования:
dxd xνJν(x) = xνJν−1(x);
dxd x−νJν(x) = −x−νJν+1(x)
и вытекающие из них рекуррентные соотношения
Jν−1(x) + Jν+1(x) = 2νx−1Jν(x);
Jν−1(x) − Jν+1(x) = 2Jν0 (x).
26
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
4. Уравнение
αJ |
(μ) + βJ0 |
(μ) = 0, α |
≥ |
0, β |
≥ |
0, α + β > 0, |
ν |
ν |
|
|
|
при ν > −1 имеет бесконечное множество действительных корней
μn (μ1 < μ2 < ∙ ∙ ∙ < μn < . . . ), причем μn → +∞ при n → ∞. 5. Если на линейном пространстве функций, интегрируемых на
[0, a], ввести скалярное умножение согласно правилу
|
|
|
|
|
|
|
|
(f, g) = Z0a xf(x) g(x) dx, |
|
|
|
||||||||||||
то система функций |
Jν |
|
μn |
x |
|
∞ |
будет ортогональной отно- |
||||||||||||||||
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
a |
|
n=1 |
|
|
|
|
|
|
|
|||
сительно этого |
скалярного умножения, т. е. |
|
|
|
|
|
|||||||||||||||||
|
|
|
n |
|
|
|
|
|
|
o |
|
|
|
|
|
|
|
|
|||||
a |
|
|
|
μ |
|
|
μ |
|
|
|
|
|
|
|
|
|
|
|
|||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Z |
xJν |
|
|
|
k |
x |
Jν |
|
n |
x |
dx = 0, если |
μk 6= μn, |
|
||||||||||
|
|
|
a |
|
a |
|
|||||||||||||||||
при этом квадрат нормы |
|
J |
2 |
функции J |
|
μn |
x |
определяется |
|||||||||||||||
|
ν a |
||||||||||||||||||||||
по формуле |
|
|
k |
|
|
|
|
k |
|
νk |
|
|
− μn2 |
|
|
|
|||||||
|
k |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
Jν |
|
|
2 = |
a2 |
J02 |
(μ ) + 1 |
|
ν2 |
J2(μ ) |
. |
(4.3) |
||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
ν |
|
|
n |
|
|
|
|
|
n |
n |
|
|
|
4.2. Краевая задача на собственные значения и собственные функции для уравнения Бесселя
Пусть ν > 0. Рассмотрим следующую краевую задачу на собственные значения и собственные функции:
|
d |
dX(x) |
|
ν2 |
|
|||
− |
|
x |
|
|
+ |
|
X(x) = λxX(x), 0 < x < a; |
(4.4) |
dx |
dx |
x |
||||||
αX(a) + |
βX0(a) = 0; |
(4.5) |
||||||
|X(0)| < |
∞. |
|
|
|
(4.6) |
|||
Здесь условие (4.6) ограниченности решения при x = 0 связано с тем, что точка x = 0 является особой точкой уравнения (4.4).
27
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Вводя новую переменную ξ = √ |
|
x и обозначая X |
√ξλ |
= |
|||||||
λ |
|||||||||||
= y(ξ), уравнение (4.4) можно привести к уравнению Бесселя |
|||||||||||
ν-го порядка: |
+ ξdξ + ξ2 − ν2 y = 0, |
|
|
||||||||
ξ2 dξ2 |
|
|
|||||||||
|
d2y |
|
dy |
|
|
||||||
поэтому общее решение уравнения (4.4) имеет вид |
|
|
|||||||||
X(x) = AJν √ |
|
x + BYν √ |
|
x , |
|
|
|||||
λ |
λ |
|
|
||||||||
где A и B — произвольные постоянные.
Так как эта функция должна удовлетворять условию (4.6), то
константа В = 0, а константу A примем равной единице, т. е. |
|
X(x) = Jν √λx . |
(4.7) |
Далее, подставляя (4.7) в граничное условие (4.5), получаем уравнение для определения собственных значений λ:
αJν √λa + β√λJν0 √λa = 0.
Oбозначая μ = √λa > 0, приходим к уравнению |
|
αaJν(μ) + βμJν0 (μ) = 0. |
(4.8) |
Пусть μn > 0 — n-й корень уравнения (4.8), тогда собственные значения и собственные функции задачи (4.4)–(4.6) примут вид
λn = |
μ |
|
2 |
, n N; Xn(x) = Jν |
μ |
|
n |
|
n |
x . |
|||
a |
|
a |
||||
Квадрат нормы собственной функции определяется по формуле (4.3).
Замечание. Если α = ν = 0, то существует собственное значение λ0 = 0, которому соответствует собственная функция
X0(x) ≡ 1.
28
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
5. РЕШЕНИЕ КРАЕВЫХ ЗАДАЧ ДЛЯ УРАВНЕНИЯ ЛАПЛАСА В ЦИЛИНДРЕ
5.1. Примеры решения задач
Задача 1 (смешанная краевая задача). Решить задачу
|
u(r, ϕ, z) = 0, |
0 ≤ r < a, |
0 ≤ ϕ ≤ 2π, |
0 < z < h; (5.1) |
||||
(αu + ur) r=a = 0; |
|
|
|
|
(5.2) |
|||
|
|
= f1(r, ϕ), |
|
|
|
(r, ϕ). |
|
|
u |
|
u |
|
|
= f2 |
(5.3) |
||
|
z=0 |
|
|
|
z=h |
|
|
|
Решение. В соответствии с общей схемой метода Фурье поставим задачу на отыскание собственных значений и собственных
∂2
функций оператора Lϕ = −∂ ϕ2 :
Lϕ[Φ(ϕ)] = νΦ(ϕ), |
0 < ϕ < 2π; |
(5.4) |
Φ(ϕ + 2π) = Φ(ϕ). |
|
(5.5) |
Здесь условие периодичности (5.5) связано с тем, что решение u(r, ϕ, z) краевой задачи (5.1)–(5.3) должно быть периодично по ϕ с периодом, равным 2π.
Напомним, что собственные значения и собственные функции задачи (5.4), (5.5) имеют вид
ν0 = 0, Φ(1)0 (ϕ) ≡ 1;
νn = n2, Φ(1)n (ϕ) = cos nϕ, Φ(2)n (ϕ) = sin nϕ, n N.
Будем искать решение u(r, ϕ, z) краевой задачи (5.1)–(5.3) в форме разложения в ряд Фурье по системе собственных функций оператора Lϕ:
u(r, ϕ, z) = v0(1)(r, z)Φ(1)0 (ϕ)+
|
∞ |
|
|
|
|
|
X |
(r, z)Φ(1) |
(ϕ) + v(2) |
(r, z)Φ(2) |
(ϕ), (5.6) |
+ |
v(1) |
||||
|
n |
n |
n |
n |
|
n=1
где v0(1)(r, z) и vn(k)(r, z), k = 1, 2, — коэффициенты Фурье функции u(r, ϕ, z), которые необходимо определить так, чтобы решение
29
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
u(r, ϕ, z), записанное в виде ряда (5.6), удовлетворяло граничным условиям (5.2) и (5.3).
Подставляя ряд (5.6) в уравнение (5.1) и учитывая равенство (5.4), приходим к следующему уравнению относительно искомых коэффициентов Фурье:
1 ∂ |
r |
∂vn(k) |
||
|
|
|
|
|
r ∂r |
∂r |
|||
0 ≤ r < a,
− n2 v(k)
r2 n
0 < z < h;
|
∂2vn(k) |
|
+ |
∂z2 = 0, |
(5.7) |
k = 1, 2, n ≥ k − 1.
При этом, согласно граничному условию (5.2), функции vn(k), k = 1, 2, n ≥ k − 1, должны удовлетворять следующему условию:
|
[αvn(k) + (vn(k))r] r=a = 0, |
(5.8) |
||
а в силу граничных условий (5.3) |
|
|
|
|
v0(1)(r, 0) Φ0(1)(ϕ)+ |
|
|
|
|
|
∞ |
|
|
|
|
X |
|
|
|
+ |
vn(1)(r, 0) Φn(1)(ϕ) + vn(2)(r, 0) |
Φn(2)(ϕ) = f1(r, ϕ); |
|
|
|
n=1 |
|
|
|
v0(1)(r, h) Φ0(1)(ϕ)+ |
|
|
|
|
|
∞ |
|
|
|
|
X |
|
|
|
+ |
vn(1)(r, h) Φn(1)(ϕ) + vn(2)(r, h) Φn(2)(ϕ) = f2(r, ϕ). |
|
||
n=1
Записанные здесь соотношения представляют собой разложения функций f1(r, ϕ) и f2(r, ϕ) в ряды Фурье по системе собственных
функций оператора L |
ϕ |
. Коэффициенты Фурье v(k) |
(r, 0) и v(k) |
(r, h) |
|||||||
|
|
|
|
|
|
|
|
n |
n |
|
|
этих разложений могут быть вычислены по формулам |
|
||||||||||
|
|
|
1 |
|
|
2π |
|
|
|
|
|
v(k)(r, 0) = |
|
|
|
|
f |
|
(r, ϕ) Φ(k)(ϕ)dϕ = b(k)(r); |
|
|||
|
|
|
|
|
|
1 |
|
||||
|
|
|
|
|
|
|
|||||
n |
k |
Φn(k) |
k |
2 |
Z |
n |
n |
|
|||
|
|
|
|
|
|||||||
|
|
|
|
|
0 |
|
|
|
|
||
|
|
|
|
1 |
|
|
2π |
|
|
|
|
vn(k)(r, h) = |
|
|
|
|
|
f2(r, ϕ) Φn(k)(ϕ)dϕ = cn(k)(r). |
|
||||
|
|
|
|
|
|
|
|||||
|
Φn(k) |
|
|
||||||||
|
k |
2 |
Z |
|
|
|
|
||||
|
|
|
|
k |
0 |
|
|
|
|
||
30 |
|
|
|
|
|
|
|
|
|
|
|