
Учебник по УМФ первая часть
.pdf
Если мы теперь продифференцируем уравнение (79) по x, а уравнение (80) по t, а
затем исключим из полученных уравнений смешанную производную ∂2i , то получим
∂x∂t
следующее линейное дифференциальное уравнение второго порядка относительно 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 |
|
|