Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие

.pdf
Скачиваний:
98
Добавлен:
22.10.2023
Размер:
13.99 Mб
Скачать

Ряды Фурье по полным системам функций {Ф„} обла­ дают следующим замечательным свойством:

Т е о р е м а 1. Если система функций {Ф„} полна

в области D, то ряд Фурье для всякой функции с инте­ грируемым квадратом в D можно почленно интегрировать независимо от того, сходится он или расходится, т. е. для любой области D' a D справедливо равенство

СО

 

 

 

\f(M)dx = £

Сп $

Фn(M)dT.

 

 

 

 

 

D ’

/1=1

D'

 

 

 

Д о к а з а т е л ь с т в о . Оценим

разность

6 n — ^ f d x

 

 

 

 

 

 

 

 

D’

— i ]

С* $ ф* dr.

 

 

 

 

 

1

 

D'

 

 

 

 

 

 

is,

$

\ [ -

Z Ckd)k\dx

 

 

 

 

 

D'

'

6 = 1

/

 

 

 

 

 

 

 

T

f - £С *Ф *

f -

Z

dx.

 

 

 

D ’

k= 1

 

 

A= 1

 

Для оценки последнего интеграла воспользуемся неравен­

ством

Коши — Буняковского:

 

 

 

 

 

 

 

/ - X С , Ф* ch = jj V 9 t - 2 c kф *

dx

^

 

 

П Т " 5

 

 

D

k=i

 

D

 

 

 

k= i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

1

 

"

 

\2

 

/

 

 

 

 

«

у

 

 

 

С*ф*

dx Л /

 

 

 

 

 

\

 

k =

 

1

 

'

 

 

 

 

 

 

D

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

th

 

 

 

Ч

/

\ p/2 л

-

У

c i ji

а д |j

 

 

 

D

о

 

 

 

 

 

 

D

 

6=1

 

 

Последний

интеграл ограничен,

а

разность

 

 

 

 

 

 

S р/2 rfT -

2

W

* l l a

 

 

 

 

 

 

 

D

 

 

*= i

 

 

 

 

 

 

при П-У со стремится

к

нулю

по условию полноты. Сле­

довательно, 8,г->0

при

п —у со,

ч.

т.

д.

функций краевой

Т е о р е м а 2.

Система

собственных

задачи (5)—(6) полна.

НО

Д о к а з а т е л ь с т в о .

Возьмем

произвольную,

непре­

рывную в D функцию f(M). Тогда, каково

бы ни было

число £ > 0 в

классе

А

(см. § 1)

найдется

такая

функ­

ция g-(M)*), что

 

 

 

 

 

 

 

 

J P ( f - g ) 2 d r < ~ .

 

(43)

 

b

 

 

 

 

 

 

 

Функция g(M)

представляется

рядом Фурье

по собствен-

ным функциям,

g{M) =

со

СкФк>равномерно сходящимся

 

в D. Следовательно,

 

k =

i

 

ех> 0 найдется

 

для всякого

такое

/V (ех), что

 

 

 

 

 

 

 

 

g - Ц скФк < е х

для

п > Л г(е1).

(44)

 

4 = 1

 

 

 

 

 

 

 

 

(

п

\2

Тогда

со-

Покажем, что

 

^

СкФк) с/т< е.

 

D

1

4 = 1

/

 

 

 

гласно достаточному признаку полноты система {Ф„} будет полной. Очевидно,

П \ 2 / п \ 2

\p[f — У,

Скфк\ с/т =

[ f - g + g -

£

СкФк) d x ^

) \

к = 1

/

D

\

 

 

4 = 1

1

 

 

 

< 2

$ р ( / - £ ) 2с/т+2

\ p l g -

С*Ф*) dx.

 

 

 

 

 

 

 

D

'

4 = 1

Используя

неравенства (43) и (44),

получим

 

 

Ц

/

11

\2

с/т < ^

+

2е]

 

р dx <

е,

 

/ -

^

СкФкI

 

 

 

 

к —1

 

 

 

 

 

 

если

взять ех ■

I

е , где

В = \

р dx.

 

 

 

 

 

 

2 (-й

 

 

 

 

 

 

Эта теорема вместе с предыдущей позволяет интегри­ ровать почленно ряды Фурье по собственным функциям краевой задачи (5) — (6) для всякой функции, интегри­ руемой с квадратом, не заботясь не только о равномер­ ной сходимости этих рядов, но даже вообще об их схо­ димости в каждой точке.

*) См. Т о л с т о в Г. П., Ряды Фурье, Физматгиз, 1960.

Ш

§ 5. Решение неоднородных краевых задач методом Фурье

Знание системы собственных функций {Ф„} и соответ­ ствующих им собственных значений {Хп} позволяет решать и неоднородные краевые задачи. Рассмотрим некоторые из них.

I. Пусть требуется найти решение задачи

 

L[u] + /(M,

i) — puit

(соответственно put),

(45)

и(М,

0) —0,

щ (М, 0) = 0,

(46)

 

 

 

(47)

непрерывное в замкнутой области

В =

е

D;

и

принадлежащее классу

А при всяком фиксированном зна­

чении t >

0.

решение

ы(Л4,

/)

принадлежит

Так

как искомое

классу А, согласно теореме Стеклова оно может быть

представлено в виде

ряда

Фурье

по собственным функ­

циям {Ф„} соответствующей

однородной задачи

(5) —(6):

и(М,

t)=

У) а д )Ф „ (М ),

(48)

где

1

П~ 1

 

 

с

 

t) ф„ (Р) dx.

(49)

ч'я (t) =

 

J ри (Р,

 

b

 

 

 

Выражая рФ„ под знаком интеграла (49) из уравнения (5), получим

_ АМФп, и] _

- 1

? ф

^«Ф я|р

К II

II2 ^ l wJ aT-

 

 

D

Выражая L [и] из уравнения (45), получим

W = М Ф Л 5 S pW/‘tlJ' ' + Я„ЦФЛ2 \ f®11dx’ (5°)

Первое слагаемое в правой части формулы (50) равно

Wn/kn. Второе слагаемое представляет известную функ­ цию, обозначим ее через f„(t)/Kn. Таким образом,

^ ( 0

- V Я + U

112

Следовательно, Чгя (/) есть решение уравнения

+

=

= fn{t) (соответственно '¥,n + hny¥n = fn) с

дополнительными

условиями

 

 

 

 

 

 

(0) = |Г ад S

(Р-

0) ф « Y

) = °>

 

 

 

ь

 

 

 

 

 

(°) = r i s

^ ри< (Л

0) фп (Р) dx = 0.

 

 

Решение такой задачи имеет вид

 

 

 

 

I

 

 

 

 

 

х¥п (t) =

[ sin V K

(t - 9) fn(9) rf0

 

 

(соответственно 4f„ (^) =

^e

«

)fn (Q)dO\.

 

 

\

о

 

 

/

 

 

Подставляя полученные функции 4^(0 в формулу (48),

получим искомое решение в виде ряда Фурье по соб­

ственным функциям. Если,

в частности,

f(M,

t) —

= f(t)&(M,

Мо),

то

 

 

 

 

 

 

 

 

h {t) =

р Ь

S f (/) 6 (P ’ Щ

Фп {P) dT = Щ

/ (0

 

 

 

ф„ (Mn)

i

 

 

(t-Q)f (9) dB

 

 

 

^ sin V K

 

 

У м фг

 

 

 

 

 

 

 

 

^соответственно ¥„(/) =

 

^ e

h*(t

9)/(0)^0j.

 

Если требуется

решить задачу

 

 

 

 

 

 

L[«]-b/(M, 0 — pU/t

(ри,),

 

 

и (М, 0) = ср(М),

ut (M,

0) = фх (М);

(yi J +

?2wjs = 0,

то будем искать

решение в виде суммы двух функций

 

и(М,

t) = v(M,

t)-\-w(M, t),

 

 

являющихся решениями следующих задач:

 

 

 

 

 

v. L[v]--^pvtt

(pv,)t

 

 

 

v(M, 0) = ср (М),

vt (М) =

cpi (М);

 

(y i^ +

Ya»)s =

0;

 

w: L [да] + / (M,

i) = pwtt

(рда,),

 

 

w(M, 0 ) = wt (M, 0) = 0;

(Y i^ +

Y2^ )s

 

Каждую из этих задач мы уже умеем решать.

113

2. Пусть требуется найти решение задачи

 

L[u]+f(M,

i)=-putt

щ),

(51)

 

(viS + Y2«)s = P(A1, ty,

(52)

и(М,

0) = Ф(/И),

ut (M,

0)=q>1 (M),

(53)

непрерывное в

замкнутой

области

В = { М е 5 ;

/5г0}.

Мы рассмотрим следующий способ решения этой задачи. Среди_функций v(M, t), непрерывных в замкну­ той области В и имеющих в этой области непрерывные частные производные первого и второго порядков, возь­ мем какую-нибудь функцию v2 (M, t), которая удовлетво­ ряет заданным краевым условиям (52). Будем искать

функцию и(М, t) в

виде суммы u = v1 (M,

t) + w(M,

t),

где для функции w(M, t), непрерывной в

области

В,

задача ставится следующим образом:

 

 

 

 

 

L[w]+fi(M,

t) = pw/i

wt),

 

 

w(M, 0) = ф(Л1),

wt(M, 0) = фх (M);

fYl^

+ Y.^ )s = o,

где

fi(M,

t) = f (M, t) + L [vj -

pvltt,

 

 

 

Ф(М) = ф ( М ) - ц 1(М, 0),

ф! (M) =

ф! (Af) — vlt (M, 0).

Эту задачу мы уже рассмотрели в п.

1. Функцию

v1(M, t) или угадывают,

или же находят методом Дюа-

меля (см. гл. V).

 

 

 

 

 

 

 

П р и м е р

6. Требуется решить задачу

 

 

 

 

 

 

 

 

=

 

 

 

 

и(х, 0) = ср! (jc),

ut (X,

0) =

ф2 (х),

и (0, f) = M

0 .

u(l,

f) = M 0 -

(*)

В качестве

функции

(л:, t), удовлетворяющей

краевым усло­

виям (*), берем функцию *)

 

 

 

 

 

 

vi (х>0 —j— H-i (0+ ~2 11з (О-

 

 

*) Мы предполагаем,

что функции щ (t)

и

(0 дважды диффе­

ренцируемы.

 

 

 

 

 

 

 

 

114

Решение и{х, t) ищем

в

виде

суммы

и (х,

t) = vl (x,

i)-^-w(x, ().

Функция

w (х, t), очевидно,

будет решением следующей задачи:

 

 

 

 

rfVxx + ^ - f - W i

 

Ц) -

*

n'/(t) = wth

 

 

w (x,

 

0) =

cPl (x) +

 

 

Pi (0) -

*

p2 (0) = ф! (x),

 

wt {x,

 

0) = ф2 (*) + ^

 

И;( 0) - *

 

И>(0) = ф2 W .

 

 

 

 

 

 

 

w (0,

t) —w(l,

t) = 0.

 

 

 

Функцию

w(x,

t) ищем

также

 

в виде суммы w = R(x,

t)-\-Q(x, t),

где R (x,

t)

есть

решение однородной

краевой

задачи

 

R(x,

 

 

 

 

 

Rt (x,

a2Rxx = Rtt<

R (0, t) = R(l,

 

0) =

$г (х),

0) = ф3 (х);

0 = 0

и имеет вид (см.

пример

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R (х,

0 =

 

(Сп cos а У Хп t +

Dn sin a Y %п t) sin — x,

 

 

 

 

 

n — 1

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

.

 

n%2

 

 

 

 

2

 

 

.

и

 

 

 

 

 

>

 

(' .

 

 

 

 

 

 

^71 —

p

 

0

\

 

(») S111 ^

 

 

 

 

 

 

 

 

 

 

 

l

 

0

 

 

 

 

 

 

 

 

 

 

 

 

2

_

 

.

nn .

 

 

 

 

 

 

 

 

 

l*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sm T

5 dg,

 

a Q (x, 0

есть решение следующей

задачи:

 

 

 

 

 

 

 

 

 

 

 

 

o2Q,vx + /(^. t)=Qth

 

 

 

 

Q(x,

 

0) = Qt (x,

0) =

0;

 

Q (0,

t) = Q(l, 0 = 0 ,

где f (x,

t) = x- j i ^ [ ' ( t )

*-цУ (0 -

 

 

 

 

 

 

Согласно

п.

1 решение имеет вид

 

 

 

 

 

 

 

 

 

 

 

Q (*, 0 = ^

 

 

(0 sin

 

 

*•

 

 

 

 

 

 

 

 

 

 

 

Я~1

 

 

 

 

 

 

 

Функции

 

 

(0

 

вычисляются

по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W =■“

 

\

sin ^

^ -

®) /я (0) d0>

 

где

 

 

 

 

 

 

Г Ая О

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

In (S) =

у

Jj / (I,

6) sin ™

E ^ =

 

 

 

 

 

л2п [ ( - 1)" n_" (e ) - ( ii' (6)J,

О

115

З а м е ч а н и е

1. Если

краевые

условия

имеют вид

а ,М 0, / ) - М ( 0 ,

() =

Pi(0,

«аМЛ 0 +

О

!Ч('>,

то в качестве функции v1 (x, t) (удовлетворяющей

этим

краевым условиям) можно взять функцию

 

 

щ (х,

/) - D xy 2 (/) - С (х -

If iu (/),

 

 

где С - 1/(2o^ + fM2).

D ~= 1/(2а2/ +

(У2).

 

 

З а м е ч а н и е

2. Иногда легко найти функцию vx{x, t),

удовлетворяющую не только заданным неоднородным крае­ вым условиям, но также и заданному уравнению.

П р и м е р 7.

Требуется решить задачу

(54)

 

Ч2Щ-л" =

 

и (х, 0) = ф! (х), щ (х, 0) = ф2 (х);

(55)

и{0,

0 = 0. u(l, t) = A s m a t ( ^ l ^ n j i ' j

(56)

Среди функций вида F (х) sin Ш нетрудно найти решение урав­ нения (54) i\(x, /), удовлетворяющее краевым условиям (56). Действи­ тельно, подставляя такую функцию в уравнение (54) и деля обе части равенства на sin со(, получим уравнение для F (х)\

 

a2F" -р со2/7 =

0.

(57)

Из краевых условий (56) находим, что

 

 

F (0) =

0,

F(l) = A.

(58)

Решение задачи (57) — (58),

очевидно,

имеет вид

 

 

 

 

.

со

 

 

F (х) =

sin -а- X

 

А

 

 

 

 

 

 

sm соI

 

Следовательно,

 

 

 

а

 

 

 

.

со

 

 

 

 

 

 

 

 

 

 

sin — X

 

V, (х,

t) =

А --------- sin соt.

 

 

 

 

.

со

,

 

 

 

 

sm

а

I

 

Решение задачи (54) —(56) будем искать в виде

и (х, t) = vx (x, t)-\-w(x, t),

где w (х, t) является решением следующей задачи:

аъ1£)хх = wit,

щ(0, 0 = 0,

со (/, 0 = 0;

 

sm -со х

w (х, 0) = фх (х), wt (х, 0) =

фг (х) — А со---------= q2 (х).

 

sin а I

Это однородная краевая задача, которая решается методом раз­ деления переменных.

116

§ 6. Применение метода Фурье к решению краевых задач для уравнений эллиптического типа

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

П р и м е р 8. Найти функцию и (г, ср), гармоническую в круге DR

радиуса R, непрерывную в замкнутой области

и принимающую

на границе этой области (r = R) заданные значения f (ср), т. е.

 

Аи = 0,

 

(59)

 

u(R, <р)= /(ф).

 

(60)

В силу однозначности

искомого решения

и (г, ф) оно должно

быть периодическим по ф с периодом 2я,

т. е.

 

и (г,

ф + 2я ) ^ м ( г ,

ф).

(61)

Из непрерывности решения в замкнутой области D^ следует его

ограниченность в D

Среди функций вида Ф (г) ¥ (ф) ищем ограниченные в 6% и периодические по ф (с периодом 2я) решения уравнения (59). Запи­ сывая лапласиан в полярных координатах

 

 

 

~

(г“л) + "2 «(fcp = 0

(59t)

и разделяя переменные,

получим

 

 

 

 

 

 

г ~ ( г Ф ’) - Х Ф = 0,

(62)

 

 

 

 

¥ " + Л¥ = 0.

 

(63)

Из

условия

(61)

находим

 

 

 

 

 

 

 

¥ ( ф +

2л)=е ¥ ( ф).

(64)

При

Д ,<0

уравнение (63)

не имеет

решений, удовлетворяющих

условию (64). Следовательно, Х ^ О . Для

^ > 0

находим

 

 

¥

(ф) = A sin УХ ф + В cos УХ ср.

Из условия

(64)

находим,

что (,'Х 2я = 2яп.

Отсюда Хп — п2, где

п — произвольное целое

неотрицательное число.

Таким образом, соб­

ственные значения задачи (63)—(64) суть

 

 

 

 

 

Хп = п2

(л= 0,1,2,...),

 

а им соответствующие собственные функции суть

1, sin ф, cos ф....... sin «ф, cos «ф, ...

117

При ^,= 0 общим решением уравнения (63) будет

^(сГ) = Л0ф + й0.

Лишь при Л0 = 0 оно будет удовлетворять условию (64). Таким образом, А —0 соответствует собственная функция Ч|'п( ф ) = 1.

Обратимся к уравнению (62). При Х — п2 имеем

/•2ф"_|_гф '_ П2ф_0.

Общее решение этого уравнения имеет вид

 

Ф» (')=<:„/■» + -£*.

(я > 0),

(65)

Фо(0 = С0 +

П0 1п

(л = 0).

(66)

В силу ограниченности

искомого

решения в формулах

(65) и

(66) надо положить £)„ = 0 (я = 0, 1,...).

Таким образом, ограниченными решениями уравнения (59^ вида Ф (г) W (ф), удовлетворяющими условию (61), будут функции

un = rn (Ап cos яф +

й„ sin яф).

 

Решение задачи (59)—(61) представится в виде ряда

 

СО

 

 

 

и (г, ф) = 2 ]

гп п cos Яф+ Б,, sin Яф).

(67)

п=о

 

 

Коэффициенты Ап и Вп находим из условия (60):

 

СО

 

 

 

/ (ф) — 2

(Л« C0S +

Вп sin ПФ) Rn>

(68)

гс —0

 

 

 

пользуясь ортогональностью собственных функций на отрезке [0, 2л]

с весом р = 1:

 

 

 

 

 

 

 

Aq=

1

f

 

 

 

1

/ (t) cos пГ dl,

\ f ©

^

 

л п = лRn ^

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

Вп =

nRn

\

/ (С)sin

=

1. 2....).

 

 

о

 

 

 

 

 

З а м е ч а н и е

1.

Ряд

(67) с коэффициентами, вычисляемыми

по формулам

(69),

нетрудно

просуммировать. Однако мы не будем

здесь этого делать, так как в гл. VII задача (59)—(61) будет решена

другим методом, позволяющим

получить результат в конечном виде.

З а м е ч а н и е

2.

Решение

краевой задачи (59)—(61) для внеш­

ности круга представляется

рядом

 

 

 

 

 

со

1

 

 

 

U(о ф) =

V

 

 

 

7.

гн (Ап cos nq> + Bn sin Яф),

 

 

 

 

4яш

 

 

 

 

 

гс = 0

 

 

 

118

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

С

О

и = ^

( C n ^ + ^ - ' j i A n cos тр-{-Вп sin rup) + Aa + B0 lnr,

коэффициенты которого (Ап, б 0, СпА,и СпВп, DnAn и DnBn) опре­ деляются из краевых условий *)

ф)=fj (ф). ф)=Мф)-

П р и м е р 9. Решить краевую задачу

Дм =

0,

 

 

(70)

и (0, у) — и (I,

у) ~ 0,

(71)

и(х, 0 ) = f l (x),

и(х,

b) = f2(x)

(72)

в прямоугольной области {O sgxsg/;

0

 

 

В классе функций вида Ф (х) ¥ (у) ищем решения уравнения (70), удовлетворяющие лишь однородным краевым условиям (71). Подстав­ ляя такую функцию в уравнение (70) и разделяя переменные,

ПОЛуЧИМ

 

ф//

т //

 

 

 

-W + - T - 0-

 

Чтобы это

равенство

было тождеством, необходимо, чтобы Ф"/Ф =

= — К, Чг"/Дг = Я, где

Я —постоянное число.

и ¥ (у):

Таким образом, получаем уравнения для функций Ф (х)

 

 

ф " +

1 ф = 0 ,

( 7 3 )

 

 

 

=

(74)

Из условий (71) находим, что

 

 

 

 

ф (0) =

Ф (/)= 0 .

(75)

Задача (73),

(75) имеет лишь положительные собственные

значения

 

кп = лая2//2

(я = 1, 2,...).

 

TCft

Им отвечают собственные функции Фя (х) = sin —j— х **). Обра­

тимся к уравнению (74). При к — кп оно имеет общее решение вида

(</)= C„ ch V K y + Dn sh \ 'кпу.

Следовательно, решения уравнения (70), удовлетворяющие лишь краевым условиям (71), имеют вид ип (х, у) = Фп (х)Чп (у).

Решение задачи (70)—(72) представляется рядом

С О

 

и (х, у) = У ( сп ch VknV + Dn Sh | кпу) sin

х,

п= 1

 

*) Читателю рекомендуется написать соответствующие формулы для определения этих коэффициентов.

**) См. гл. IV, § 3, пример I.

119

Соседние файлы в папке книги из ГПНТБ