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

Учебник по УМФ первая часть

.pdf
Скачиваний:
102
Добавлен:
16.04.2015
Размер:
968.33 Кб
Скачать

Если мы теперь продифференцируем уравнение (79) по x, а уравнение (80) по t, а

затем исключим из полученных уравнений смешанную производную 2i , то получим

xt

следующее линейное дифференциальное уравнение второго порядка относительно v:

2v

= LC

2v

+(RC +GL)

v

+GRv

x2

t2

t

 

 

 

Аналогичным образом получим дифференциальное уравнение относительно i :

2i

= LC

2i

+(RC +GL)

i

+GRi

x2

t2

t

 

 

 

(81)

(82)

В результате получим, что напряжение v и сила тока i удовлетворяют одному и тому же дифференциальномууравнению:

2u

= a

2u

+ 2b

u

+cu ,

(83)

x2

t

2

t

 

 

 

 

где a = LC, 2b = RC +GL, c = GR .

Это уравнение называют телеграфным уравнением. Нетрудно заметить, что при

с = 0 оно по форме совпадает с уравнением колебания струны (или стержня) с учетом процесса затухания.

§ 11. Общая схема метода разделения переменных для одномерных гиперболических уравнений

В § 5 мы познакомились с методом разделения переменных, который является

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

однородной струны или однородного стержня, а именно:

 

 

u

 

2u

 

,

(84)

Lu(x)

 

k(x)

 

q(x)u = ρ(x)

t2

0 x l

 

 

x

x

 

 

 

 

где k(x), q(x) и x(ρ) – непрерывные на отрезке 0 ≤ x l, положительные функции (k >0, ρ>0, q ≥0). Для этого уравнения сформулируем начальные условия

u(x,0) = f (x)

t ≥ 0

(85)

 

ut (x,0) = g(x)

и граничные условия (однородные)

u(0,t) = 0

(86)

u(l,t) = 0

 

В соответствии со знакомым нам методом разделения переменных будем искать нетривиальное решение уравнения (84) в виде произведения

u(x,t) = X (x)T (t)

(87)

Подставляя предполагаемый вид решения (87) в уравнение (84) и пользуясь граничными условиями, после операции разделения переменныхполучим два уравнения

L X (x) +λρX (x) = 0

(88)

T ''(t) + a2λT (t) = 0

(89)

где λ – неизвестный пока параметр.

В силуод нородности граничных условий (86) мы можем записать однородные граничные условия и для уравнения (88):

X (0) = 0 и X (l) = 0

(90)

Для определения функции X (x) мы получаем известную уже задачу ШтурмаЛиувилля, в которой надо найти такие значения параметра λ, при которых существует

нетривиальное решение уравнения (88), удовлетворяющее граничным условиям (90), а также найти эти решения. Удовлетворяющие условиям задачи значения параметра λ

называются собственными значениями или собственными числами, а соответствующие им нетривиальные решения – собственными функциями задачи Штурма-Лиувилля.

В силу однородности уравнения (88) и граничных условий (90) собственные

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

1

ρ (x) Xk2 (x) dx =1 ,

(91)

0

 

 

котрое представляет собой условие нормировки собственных функций.

Собственные функции, удовлетворяющие условию (91), называют нормированными. Сформулируем некоторые общие свойства собственныхфункций и собственныхзначений

задачи Штурма-Лиувилля.

1.Существует счетное множество собственных значений λ1, λ2,…, λn,…, которым соответствуют собственные функции X1(x), X2(x),…, Xn(x),…

2.Все собственные числа при q ≥ 0 – положительны.

3.Собственные функции ортогональны междусобой с весом ρ(x) на отрезке 0≤ x l,

т.е.

1

ρ (x) Xn (x) Xm (x) dx = 0

(92)

0

 

 

4.Всякому собственному числу соответствует с точностью до постоянного множителя только одна собственная функция.

5.Произвольная дважды непрерывно дифференцируемая функция F(x), удовлетворяющая граничным условиям F(0) = F(l) = 0 , разлагается в равномерно и

абсолютно сходящийся ряд по собственным функциям {Xn(x)}:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(x) =

Fn Xn (x) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

причем

 

 

 

 

 

1

 

1

 

 

(x) ρ (x) dx,

(93)

 

 

 

F

=

 

 

 

 

F(x) X

 

 

 

 

 

 

Xn

 

 

 

 

 

 

 

n

 

 

 

2

0

 

n

 

 

 

Xn

 

2

1

Xn2 (x) ρ (x) dx

 

 

 

 

 

 

 

 

 

 

где

 

= ∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

Это свойство называюттеоремой разложимости Стеклова.

Доказательство свойств 1 и 5 обычно основывается на теории интегральных уравнений, которая выходит за пределы настоящего курса. Мы остановимся лишь на доказательствахсвойств 2, 3 и 4.

Обратимся к доказательству свойства 3. Его удобно доказать, воспользовавшись одной из формул Грина. Мы докажем эти формулы в главе VI. Здесь нам понадобится формула Грина для дифференциального оператора, представляющего левую часть

уравнения(84), а именно

 

 

u

 

 

Lu =

 

k(x)

x

 

q(x)u

 

 

x

 

 

Возьмем две произвольные

 

функции

 

Xn (x) Xm (x), дважды непрерывно

дифференцируемые на интервале a < x < b и имеющие непрерывную первую производную на отрезке a x b . Рассмотрим теперь выражение

uL[v]- vL[u]=u(kv' ) v(ku' )=[k(uv' ) k(vu' )]'

(94)

Интегрируя это равенство поx от a до b, мы и получим нужную нам формулуГрина

 

b

(uL[v]- vL[u])dx= k[(uv' ) (vu' )]

 

b

(95)

 

a

 

 

a

 

 

 

 

 

 

 

 

Возьмем теперь две собственные функции Xn(x) и Xm(x), соответствующие собственным значениям λn и λm. Если мы положим в формуле (95) u = Xm(x), v = Xn(x), а в

качестве границ интервала возьмем значения 0 и l, то учитывая граничные условия (90) получим

l

(Xm L[ Xn ]- Xn L[ Xm ])dx = 0,

(96)

o

 

 

 

откуда, пользуясь уравнением (88), получаем

 

 

l

(Xm (x) Xn (x) ρ(x)d x = 0

 

(λn λm )

(97)

 

o

 

 

Таким образом, если λn λm, тонулю равен интеграл

 

 

l

(Xm (x) Xn (x) ρ(x)d x = 0

(98)

 

 

o

 

 

что и выражает условие ортогональности собственных функций Xn(x) и Xm(x) с весом ρ(x).

Докажем теперь свойство 4.

Как мы знаем, всякая собственная функция, будучи решением обыкновенного дифференциального уравнения второго порядка, определяется однозначно по заданному

значению самой функции и её первой производной в некоторой точке, например при x=0. Допустим существование двух функций X1 и X2, отвечающиходномуи томуже значению

λ и обращающихся в соответствии с граничными условиями в нуль при x=0. Тогда, если мы возьмем функцию

X

12

(x) =

X1/

(0)

X

2

(x),

(99)

X2/

(0)

 

 

 

 

 

то убедимся, что она удовлетворяет тому же уравнению второго порядка (88) и тем же начальным условиям, что и функция X1

 

/

(0)

 

 

X12 (0) =

X1

X2 (0) = 0,

(100)

/

(0)

 

X2

 

 

dX12 (0) =

X1/ (0)

X /

(0) = X / (0) .

(101)

X2/ (0)

 

dx

 

2

 

1

 

Тем самым доказано, что X12 (x) = X1(x) и что

 

 

 

 

X

12

(x) =

A X

2

(x),

где А=

X1/ (0)

 

X2/ (0)

 

 

 

 

 

 

 

 

 

Отметим, что при определении функции X12 негласно предполагалось, что X2/ (0) 0 . Это предположение вполне оправдано, так как решение линейного уравнения (88)дляэтой функции при начальныхусловиях

X2 (0) = 0, X2/ (0) = 0

тождественно равно нулю и не может быть собственной функцией. В силулинейности и однородности уравнения (88)

Наконец, докажем свойство 2.

Пусть Xn – нормированная собственная функция, соответствующая собственному значению λn , такчто в соответствии с (88) мы можем написать равенство

L[ Xn (x)] = −λn ρ(x)Xn (x)

(102)

Умножая обе части этого равенства на Xn иинтегрируяпо x от 0 до l, получаем

l

Xn2

l

Xn (x)

λn

(x)ρ(x)d x= −∫

0

 

0

 

тогда с учетом нормировки функции Xn получим

l

 

d

 

 

λn = −∫

Xn (x)

k( x ) d Xn d x

 

0

 

dx

d x

L[ Xn (x)]d x ,

l

+ ∫q(x) Xn2 (x) d x .

0

Интегрируя последнее выражение по частям и пользуясь граничными условиями, получаем

 

 

 

l

 

l

λn = − Xn (x) k Xn(x)

 

l0

+ ∫k(x) [Xn(x)]2

d x + ∫q(x) Xn2 (x) d x=

 

 

 

 

0

 

0

l

 

l

 

= ∫k(x) [Xn(x)]2 d x + ∫q(x) Xn2 (x) d x,

0

 

 

 

0

 

откуда и следует, что λn>0, посколькупо условию k (x) >0 и q(x) ≥0.

Закончив обсуждение свойств собственных функций уравнения (88), обратимся к уравнению (89). Решение этого уравнения известно и для каждого λn оно будет иметь вид

Tn (t) = An cos(

λn

t) + B nsin(

λn

t)

(103)

где Ак и Вк – произвольные постоянные. Тогда каждая функция

u n (x,t) = Xn (x)Tn (t) =[An cos(λn t) + B n sin (λn t)]Xn (x)

будет решением уравнения (84), удовлетворяющим граничным условиям (86). Чтобы получить общее решение уравнения (84) запишем суммуряда

u (x,t) = [An cos (λn t) + B n sin (λn t)]Xn (x)

n=1

Чтобы эта сумма ряда была решением уравнения (84), нужно, чтобы сам ряд и ряды,

полученные двукратным почленным дифференцированием, сходились равномерно по x и t. Для определения коэффициентов Ak и Bk воспользуемся начальными условиями (85) и теоремой разложимости Стеклова. Для этого запишем начальные условия в виде

u(x,0) = An Xn (x) = f (x)

n

ut (x,0) = Bn λn Xn (x) = g (x)

n

Умножим теперь обе части этих равенств на ρ(x)Xk (x) и проинтегрируем по x от 0 до l.

Тогда воспользовавшись условием ортогональности (92) и нормировки (91) мы получим выражения

l

 

 

1

 

l

 

An = ρ(x) f (x)Xn (x)dx и

Bn =

 

 

ρ(x)g(x)Xndx

(104)

 

 

 

λ

0

 

 

n

0

 

§12. Задача Гурса

Вэтом параграфе мы коротко познакомимся с задачей определения решения уравнения в частных производных по данным на характеристиках. Такую задачу

называют задачей Гурса.

Рассмотрим задачу о колебании полубесконечной струны в случае, когда a2 =1. В этом случае уравнение будет иметь вид

utt = uxx

(105)

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

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

поскольку для переменной t в начальных условиях (при t = 0) задаются две функции u(x,0) и ux (x,0) , а для переменной x в граничном условии (при x=0) задается только одна

функция u(0,t) .

Наряду с постановкой дополнительных условий на прямых x=0 и t = 0 на фазовой плоскости (x,t) эти условия можно ставить и на кривых в фазовой плоскости (рис. ).

Например, граничные условия можно задавать на некоторой линии С1: x = f1(t) , а начальные условия на некоторой линии С2: x = f2 (t) . В качестве таких кривых могут

выступать характеристики уравнения в частныхпроизводных.

Для некоторых уравнений координатные прямые также могут являться характеристиками. Рассмотрим, например, простейшую задачу, в которой требуется решить уравнение

uxy = f (x, y)

(106)

с дополнительными условиями, заданными на прямых x=0 и y = 0, являющихся характеристиками уравнения (105):

u(x,0) =ϕ1

(x)

(107)

u(0, y) =ϕ2 (y)

 

Интегрируя уравнение (105) последовательно уравнение (105) поx и y, получим

x

uy (x, y)=uy (0, y)+f (ξ, y)dξ,

0

y x

u(x, y)=u(x,0)+u(0, y) u(0,0)+ ∫ ∫ f (ξ,η)dξdη

0 0

или

y x

u(x, y) =ϕ1(x) +ϕ2 (y) ϕ1(0)+∫ ∫ f (ξ,η)dξdη (108)

0 0

Таким образом, для простейшего уравнения, не содержащего первых производных

ux, uy и искомой функции, решение представляется в аналитической форме (108). Рассмотрим теперь более общий вид линейного гиперболического уравнения

uxy = a(x, y)ux +b(x, y)uy +c(x, y)u + f (x, y)

(109)

при дополнительныхусловияхна характеристиках x=0 и y = 0

u(x,0) =ϕ1

(x)

(110)

u(0, y) =ϕ2 (y)

 

где φ1 и φ2 удовлетворяют требованием дифференцируемости, а коэффициенты a, b и c – непрерывные функции x и y.

Последовательно интегрируя уравнение (109) по x и y, получим для функции u(x, y) интегродифференциальное уравнение

 

y x

ξ

η

 

1

 

2

1

 

y x

 

 

 

∫∫

 

 

∫∫

 

(111)

u(x, y)

a(ξ,η)u

+b(ξ,η)u

+c(xξ,η)u dξdη +ϕ

(x) +ϕ

 

(y) ϕ

(0)+

 

f (ξ,η)dξdη

 

0 0

 

 

 

 

 

 

 

 

0 0

 

 

Это уравнение решается приближенно с помощью метода последовательных приближений, что выходит за пределы настоящего курса.

§ 13. Теорема единственности решения краевых задач для одномерного волнового уравнения

При решении граничных задач следует убедиться в существовании решения и его

единственности. Существование решения либо устанавливается в процессе изложения метода решения, либо доказывается теорема существования. Единственность устанавливается доказательством соответствующей теоремы единственности.

В настоящем параграфе приводится доказательство теоремы единственности для одномерного волнового уравнения достаточно общего вида при заданных начальных

условиях и граничных условиях первого, второго и третьего рода. Эта теорема формулируется следующим образом.

Существует только одна функция u(x, t), определенная в области 0 x l, при t > 0, которая удовлетворяет уравнению

 

2u

 

 

u

+ F(x,t) ,

( ρ(x) > 0,

k(x) > 0)

(112)

ρ(x)

t2

=

 

k(x)

 

 

 

 

x

x

 

 

 

 

 

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

 

 

 

 

 

 

u(x,0) =ϕ(x),

ut (x,0) =ψ (x)

 

(113)

u(0,t) = µ1(t),

u(l,t) = µ2 (t)

(114)

если выполнены следующие условия:

1)функция u(x, t) вместе со своими первыми и вторыми производными непрерывна на отрезке 0 x l, при t > 0,

2)коэффициенты ρ (x) и k (x) непрерывны на отрезке 0 x l .

Допустим, что существует два решения рассматриваемой задачи:

 

u1(x,t) и u2 (x,t)

 

 

Рассмотрим разность

 

 

 

 

 

 

 

 

 

v(x,t) = u1(x,t) u2 (x,t)

 

Функция v(x,t) , очевидно, удовлетворяет однородному уравнению

 

 

2u

 

 

 

u

 

(115)

ρ(x)

t2

=

 

 

k(x)

 

 

x

x

 

 

 

 

 

 

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

 

 

 

v(x,0) = 0,

 

 

vt (x,0) = 0

(116)

v(0,t) = 0,

 

 

 

v(l,t) = 0

(117)

а также условию 1) теоремы.

 

 

 

 

 

 

 

 

 

Докажем что функция v(x,t) тождественно равна нулю. Для этого

рассмотрим

функцию

 

 

 

 

 

 

 

 

 

E(t) =

1 l

 

 

 

2

+ ρ (vt )

2

(118)

2 k(vx )

 

d x

 

0

 

 

 

 

 

 

 

 

и покажем, что она не зависит от t. В задаче о колебании струны эта функция

представляет собой полную энергию струны в момент времени t. Продифференцировав E(t) по t получим

dE(t)

 

l

 

 

=

(kv v

+ ρ v v ) d x

dt

 

x xt

t tt

 

 

 

 

 

0

 

 

Интегрируя первое слагаемое по частям, получим

l

l

 

kvxvxt dx =[kvxvxt ]l0 vt (kvx )x dx

(119)

0

0

 

В силу условий (116) и (117) первое слагаемое правой части равно нулю, следовательно

dE(t)

l

l

= [ρ vtvtt vt (kvx )x ]d x = vt [ρ vtt (kvx )x ]d x = 0 ,

dt

0

0

 

т.е. E(t) = const . Тогда учитывая начальные условия, получаем

E(t) = const=E(0) =

1 l

 

2

+ ρ (vt )

2

 

(120)

2

k(vx )

 

 

t=0 d x = 0

 

 

0

 

 

 

 

 

 

а тогда и формула (118)принимает вид:

1 l k(vx )2 + ρ (vt )2 d x = 0 2 0

откуда, учитывая положительность ρ (x) и k (x), заключаем, что

 

 

 

vx (x,t) = 0

и vt (x,t) = 0

 

 

 

 

 

(121)

Это, в свою очередь, означает, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v(x,t) = const = C0 ,

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v(x,0) = C0 = 0 ,

 

 

 

 

 

 

 

А тем самым доказано, что

 

v(x,t) 0 . Следовательно,

 

если существуют две функции

u1(x,t) и u2 (x,t), удовлетворяющие всем условиям теоремы, то

 

 

 

 

 

 

 

 

 

u1(x,t) u2 (x,t)

 

 

 

 

 

 

(122)

Единственность решения задачи с граничными условиями второго рода

доказывается аналогично. Введенная в рассмотрение функция

v = u1 u2

в этой задаче

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vx (0,t) = 0,

 

vx (l,t) = 0 ,

 

 

 

 

(123)

выполнение которыхтакже приведет кобращению в нуль

 

первого слагаемого в формуле

(119). Дальнейшее доказательство проводится также каки для первой краевой задачи.

Для третьей краевой задачи

также рассматриваются два решения u1

и u2. Тогда для

функции v = u1 u2 граничные условия будут однородными:

 

 

 

 

 

 

v(0,t) h1v(0,t) = 0

(h1 0)

 

 

 

(124)

 

 

v(l,t) + h1v(0,t) = 0

(h1 0)

 

 

 

 

 

 

 

 

 

Теперь представим первое слагаемое в формуле в виде

 

 

 

 

 

 

 

[kvxvxt ]l = − k

h2v2 (l,t) + h1v2

 

(l,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

2 t

 

 

 

 

 

 

 

 

 

 

и проинтегрируем d E в пределахот 0 до t. В результате получим

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

l

 

 

(kvx )x

]d x dt

 

 

 

 

 

 

 

E(t) E(t) =

∫∫vt [ρvtt

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

{h2

 

2

(l,t) v

2

 

 

2

(0,t) v

2

(0,0)

 

 

2

v

 

 

(l,0)

+ h1 v

 

 

},

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда в силу уравнения для v и граничныхусловий следует, что

 

 

 

 

 

E(t)= k h

v2

(l,t) h v2 (0,t) 0 ,

 

 

 

(125)

 

 

 

 

 

2 2

 

 

 

 

1

 

 

 

 

 

 

 

но в силу неотрицательности подынтегральной функции

 

E(t) должно быть больше или

равно нулю. Из чего следует, что E(t) 0, а следовательно и

 

 

 

(126)

 

 

 

 

 

v(x,t) 0

 

 

 

 

 

 

 

Тем самым единственность третьей краевой задачи доказана.

Гл а в а III. Двумерные и трехмерные задачи для волнового уравнения

Вдвумерном и трехмерном случаях волновое уравнение можно записать следующим образом

u a2

u = F(M ,t) ,

(1)

tt

 

 

где М – точка на плоскости или в пространстве. Двумерное волновое уравнение обычно связывают с задачей о колебании мембраны. К трехмерному волновому уравнению приводятся задачи о течении жидкости, скорость которой имеет потенциал, о распространения звука в газе, о распространении электромагнитных полей в

непроводящей среде и задачи теории упругости. Эти задачи мы рассмотрим в последующих главах.

Во всех этих задачах волновое уравнение описывает процесс распространения волн. В одномерном случае мы в этом убедились при р ассмотрении формулы Даламбера. Такие волны называют плоскими волнами. Однако и в трехмерном пространстве в случае

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

плоскими. Двумерные волны называются цилиндрическими волнами, а трехмерные –

сферическими волнами.

Выведенную нами ранее формулу Даламбера можно обобщить соответствующим образом на двухмерный и трехмерный случаи. Начнем с трехмерного случая.

§ 1. Волны в трехмерном пространстве

Сферически симметричная задача

Сначала рассмотрим для однородного уравнения

u = a2

u ,

(2)

tt

 

 

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

искатьрешения уравнения(2)

u(M ,t) = u(r,t) ,

где r – расстояние междуточками М и М0. В этом случае уравнение (2) после записи его в сферической системе координат можно свести к одномерному уравнению для функции

v = ru

v =

1 v

(3)

rr

a tt

 

Причем, если функция u(r,t) ограничена

при r = 0 , то функция

v = ru при r = 0

обращается в нуль. В результате задача Коши для уравнения (2) с начальными условиями

u(r,0) =ϕ(r) и u t(r,0) =ψ (r)

(4)

сводится к задаче о колебаниях полуограниченной струны (0 r < ∞) с закрепленным концом в точке r = 0 :

vrr = 1a vtt

v(r,0) = rϕ(r)

(5)

vt(r,0)

= rψ (r)

 

v(0,t) = 0

Эта задача была нами решена в §4 гл. II, поэтому, не умаляя общности, мы можем записать общее решение уравнения (3) в более удобном для нас виде

v(r,t) = f1(t ar ) + f2 (t + ar ) ,

где f1 и f2 – произвольные дважды дифференцируемые функции. Тогда для функции u(r,t) будем иметь

 

 

 

u(r,t) =

1

f

(t

r

) +

1

f

2

(t +

r

)

(6)

 

 

 

 

r

 

r

 

 

 

 

 

 

 

1

 

a

 

 

 

a

 

 

 

 

Слагаемые в правой части (6) представляют собой частные решения уравнения (2)

u

1

(r

,t) = 1

f (t

r

)

и

 

u

2

(r,t) = 1

f

2

(t +

r

)

a

 

 

 

 

r

1

 

 

 

 

 

 

 

 

r

 

 

 

a

и являются сферическими

волнами;

u1(r,t)

есть

расходящаяся сферическая волна, а

u 2 (r,t) – сходящаяся сферическая волна. В отличие от плоских волн, сферическая волна

убывает обратно пропорционально расстоянию от центра.

Учитывая теперь нулевое граничное условие v(0,t) = 0, получим

f1(t) + f2 (t) = 0 или f2 (t)

Тогда решение (6) примет вид

u(r,t) = 1r f (t + ar ) + 1r

= − f1(t) = f (t)

f (t ar )

и при r = 0 , воспользовавшись формулой Лагранжа, можем записать u(0,t) = 2r f '(t)

Формула Пуассона

Теперь решим однородное волновое уравнение

utt = a2 (uxx +uyy +uzz )

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

u(x, y, z,0) =ϕ (x, y.z) u t(x, y, z,0) =ψ (x, y.z)

(7)

(8)

(9)

(10)

Будем предполагать, что φ(x,y,z) непрерывна вместе со своими производными до третьего порядка, а ψ (x,y,z) – до второго порядка включительно во всем пространстве.

Покажем сначала, что интеграл

u(x, y, z,t) =

1

 

w(ξ,η,ς) dσr ,

(11)

4πa ∫∫

 

r

 

 

 

Sat