Материал: Решение краевых задач для уравнения Лапласа (96

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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) можно привести к уравнению Бесселя

ν-го порядка:

+ ξ+ ξ2 − ν2 y = 0,

 

 

ξ2 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

 

 

 

 

 

 

v(k)(r, 0) =

 

 

 

 

f

 

(r, ϕ) Φ(k)(ϕ)dϕ = b(k)(r);

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

n

k

Φn(k)

k

2

Z

n

n

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

vn(k)(r, h) =

 

 

 

 

 

f2(r, ϕ) Φn(k)(ϕ)dϕ = cn(k)(r).

 

 

 

 

 

 

 

 

 

Φn(k)

 

 

 

k

2

Z

 

 

 

 

 

 

 

 

k

0

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

Источник: https://studfile.net/preview/16723585/