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

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

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

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

Добавлен: 16.10.2024

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

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

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

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

Возьмем

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

непре­

рывную в 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