/
Автор: Просветов Г.И.
Теги: физика задачи по физике
Текст
ПРЕДИСЛОВИЕ
ГЛАВА 1
ЗАДАЧИ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
i.l. ПРИМЕРЫ УРАВНЕНИЙ МАТЕМАТИЧЕСКОЙ
ФИЗИКИ
Теория дифференциальных уравнений с частными производными
юзникла из рассмотрения некоторых важных задач физики и механи-
си. Наиболее часто встречаются дифференциальные уравнения второ-
•о порядка, которые в основном и будут предметом нашего изучения.
Дифференциальным уравнением с частными производными второ-
10 порядка называется соотношение, связывающее независимые пере-
менные xi,..., хп, неизвестную функцию u(xi,..., хп) и ее частные
Производные до второго порядка включительно:
. ди ди д'2 и -_______
ф(жь.. жП) и,—(i,j = i,n)) = о. (1-1.1)
oxi дхп oxidxj '
Материал пособия охватывает традиционные разделы теории л
нейных дифференциальных уравнений с частными производными, вс
никающих из рассмотрения некоторых важных задач физики и мех
ники. В главе 1 приводятся примеры физических задач, приводят
к дифференциальным уравнениям с частными производными, а та
же дается постановка и классификация задач математической физик
Глава 2 посвящена исследованию уравнений гиперболического тиг
которые возникают, в основном, при моделировании колебателью
процессов. В этой главе изучаются важнейшие задачи для уравнен!
гиперболического типа и представлены методы их решения. В гла
3 приводятся решения основных задач для уравнений параболическ
го типа, которые наиболее часто встречаются при изучении проц<
сов теплопроводности и диффузии. Глава 4 содержит теорию решен
краевых задач для уравнений эллиптического типа, возникающих, К4 ''равнение (1.1.1) называется линейным относительно старших про-
правило, при моделировании различных стационарных процессов щ
описании стационарных полей. В главе 5 приводится теорема Кош:
Ковалевской, которая является фундаментальной теоремой о сущее
вовании решения задачи Коши для весьма общих систем уравненш
частными производными.
вводных, если оно имеет вид
А , ч 32и . ди ди , ,,
, Хп)^ + Р(Х!..........xn> U, = 0. (1.1.2)
''равнение (1.1.2) называется линейным, если оно имеет вид
А / \ ®2и A L I \^U I ! \
А ан(хь • > Хп)-л—я-1- 2L о>(®1> • • • > Хп)-х— + t>o(xi,..., х„)и
CJ X OXj i—1 (УХ{
= /(xi,...,xn). (1.1.3)
Приведем несколько примеров задач математической физики, при-
водящих к дифференциальным уравнениям с частными производными.
1- Сравнение колебаний струны. Допустим, что в положении
'авновесия струна совпадает с осью х, и для простоты предполо-
<им, что колебания происходят в одной плоскости, и все точки струны
5
движутся перпендикулярно к оси х. Тогда процесс колебаний мож,
описать одной функцией и(т, 4), которая характеризует вертикально
смещение струны в точке х в момент времени 4. Предположим, ч-
струна однородна, толщина ее постоянна, что она нерастяжима и j
сопротивляется изгибу. Ограничимся только ’’малыми” колебания!,
струны, когда (u®)2 << 1. Тогда функция и(т,4) удовлетворяет сл
дующему уравнению ([1, с. 23]; [2, с. 56]):
utt = а2ихх + /(г, 4),
(1.1.4
где а > 0 - константа, зависящая от физических свойств струны,
/(х,4) характеризует внешние силы. Более общее уравнение одномеи
ных колебаний имеет вид
p(x)utt = (k(x)ux)x - q(x)u + f(x, t)y (1.1.
где функции р(ж), k(x), q(x) определяются свойствами среды, в кот<
рой происходят колебания. Например, уравнение продольных колеб
ний стержня (пружины) имеет вид (1.1.5), где и(т,4) - смещение то,
ки х от положения равновесия в момент времени 4, р(х) - плотное!
стержня, к(х) — модуль Юнга (коэффициент упругости), д(х) = С
- плотность внешних сил, рассчитанная на единицу длины [;
с. 27].
2. Волновое уравнение. Многомерные колебательные процеса
в п-мерном пространстве переменных х — (г,,... ,хп) описываютв
так называемым волновым уравнением
utt = а2 Ди + f(x, t),
(1.1.1
где Ди := ^2 g-j _ оператор Лапласа. При n = 1 (1.1.6) совпала! де температура в точке х в момент времени 4, а
с уравнением колебаний струны (1.1.4). При п — 2 уравнение (1.1.1
описывает, например, двумерные поперечные колебания тонкой меЛодности в одномерном случае имеет вид
браны [1, с. 31]. Колебания в трехмерном пространстве описывают,
уравнением (1.1.6) при п = 3. Более общее волновое уравнение име<
вид
р(х)ии = div(fc(z)gradu) - д(х)и + /(ж, 4),
div(fc(z)gradu) := £
i=l OXj OXj
(1.1.8)
I функции р(т), А(т),д(х) описывают свойства среды.
I з Телеграфное уравнение. Важным частным случаем волно-
I уравнения является так называемое телеграфное уравнение. Про-
L дение электрического тока по проводу описывается следующей сис-
Смой уравнений:
ix + Cvt 4" Ov — О,
vx + Lit + Ri = 0,
I i - сила тока, v - напряжение, C, L,R,G - коэффициенты, xa-
[актеризуюшие емкость, индуктивность, активное сопротивление и
[течку соответственно (рассчитанные на единицу длины) [2, с. 88]. Из
1.1.8) получаем
wxx = aowtt + 2bowf + egw, (1.1.9)
де w = i или v, a0 = LG, 26o = RC + GL, cq = GR. Уравнение (1.1.9)
Называется телеграфным уравнением. Замена w = uexp(—bgt/ao) при-
водит (1.1.9) к более простому виду
utt = a2uxx + Ви, (1.1.10)
де а = ао1/2, В = ао2(Ь§ - аосо). Заметим, что В = 0 тогда и толь-
Но тогда, когда RC = GL (условие Хевисайда), что соответствует
ртниям без искажений.
4. Уравнение теплопроводности. Рассмотрим однородный стер-
кень, теплоизолированный с боков и достаточно тонкий, чтобы в лю-
бой момент времени температуру во всех точках фиксированного по-
перечного сечения можно было считать одинаковой. Процесс распро-
странения тепла в стержне описывается следующим дифференциаль-
ным уравнением [1, с. 180]:
ut = а2ихх + /(х,4),
(1.1.11)
О - кон-
такта, зависящая от физических свойств материала стержня, /(я, 4)
- плотность тепловых источников. Более общее уравнение теплопро-
р(х)щ = (fc(x)u®)x - q(x)u + /(г, 4).
I Многомерное уравнение теплопроводности в п - мерном пространстве
(1.1.‘ 1еРеменных х = (ц,..., г„) имеет вид
щ = а2Аи + /(ж, 4)
(1.1.12)
7
6
или в более общем случае
р(х)и( = div(fc(x)gradu) - g(x)u + /(х, t).
5. Уравнение Лапласа. Уравнение
Ди = О
(1.1.14
называется уравнением Лапласа. К нему сводится решение многи
I ! есь и(х>У) ~ неизвестная функция, a ajk = Ojk(x,y) - заданные не-
I ывные в некоторой области функции, не обращающиеся в нуль
(1.1.1; ’ ноВременно. Свойства решений уравнения (1.2.1) в основном опре-
ляются членами со старшими производными, т.е. коэффициентами
(х !/) Проведем классификацию уравнений вида (1.2.1) по старшим
роизводным.
I (Определение 1.2.1.
физических и математических задач. Например, рассмотрим устано "равнение (1.2.1) называется гиперболическим, если “12 - “п“22 > 0.
L Сравнение (1.2.1) называется параболическим, если а{2 — = 0-
Падение (1-2.1) называется эллиптическим, если а?2 “ °па22 < 0-
предположи] равнение \
Тип уравнения (1-2.1) определяется поточечно. Если во всех точ-
вившийся во времени температурный режим в некоторой области п
мерного пространства переменных х = (х,,..., х„) и i ,
что источники тепла отсутствуют. Тогда температура будет зав1
сеть только от х, и уравнение (1.1.12) переходит в уравнение Лапла< ах некоторой области уравнение (1.2.1) имеет один и тот же тип, то
(1.1.14).
Неоднородное уравнение Лапласа
одходящей заменой переменных (т.е. переходом к новой системе ко-
рдинат) можно добиться, чтобы уравнение имело наиболее простой
канонический) вид. Сделаем замену переменных
f = ¥>(ac,J/), г] = ф(х,у),
Ди = f(x)
(1.1.1!
называется уравнением Пуассона. Например, гравитационный (эле1
тростатический) потенциал и(х), порожденный телом с плотность де <р,ф - дважды непрерывно дифференцируемые функции, причем
массы (заряда) р(х), удовлетворяет уравнению (1.1.15) при /(х)
—4тгр(х). Более общее уравнение, описывающее стационарные проце<
сы, имеет вид
(1.2.2)
Рх Фх
| <Ру Фу
1реобразуя производные к новым переменным, вычисляем
(1.2.3)
0 0.
div(A:(x)gradw) — g(x)u = /(г).
Другие примеры физических задач, приводящих к уравнениям
частными производными, и вывод основных уравнений математичек
кой физики см. в [1-6].
1.2. КЛАССИФИКАЦИЯ УРАВНЕНИЙ
ВТОРОГО ПОРЯДКА
1. Классификация уравнений второго порядка с двумя не
зависимыми переменными. Рассмотрим уравнение второго поряд
ка с двумя независимыми переменными х, у, линейное относителы
старших производных
auUxx + Занижу + атиуу + F(x, у, и, их, иу) = 0.
(1.2.1
их = U(ipx 4- иуфх, иу = U(y>y 4- щфу,
ихх = UtfiPx 4- 2и^<рхФх + иут)Фх 4" "U^^Pix + иуф
XXI
иху = UK<PxFy + Щч&хФу + Фуфх) + ОупФхфу 4- и(<Рху +
иуу = о((<ру 4- 2и&РуФу 4- и^ф2 4- и(у>уу 4- и^фуу.
Тодставляя эти выражения в (1.2.1), получаем
“и“« 4- 2di2u^ 4- а22и^ 4- Ё(£, т;, и, и5, иу) = О,
де
“11 = Оц(р* + 2а.12Ч>х<Ру + аггЧ>2у,
ап = atiy>xi/>x 4- а12(<рхфу + Фх¥>у) + а22^уФу,
022 = “11 + 2О12фхФу + “22^-
(1.2.4)
(1.2.5)
8
9
Замечание 1.2.1. Так как а22 — аца22 = (ai2 — аца22)И2, rj
П = (рхфу — ^xVy / О, то отсюда следует инвариантность типа и;
замене переменных (1.2.2).
Выберем замену (1.2.2) так, чтобы уравнение (1.2.4) имело наиб
лее простой (канонический) вид. В (1.2.2) в нашем распоряжении д
функции ip и i/>. Покажем, что их можно выбрать так, чтобы выпо
нялось одно из следующих условий: 1) ац = а22 = 0; 2) ац = aj2 =
(или симметрично а22 = а12 — 0); 3) ац = а22, а12 = 0. Пусть д:
определенности ац 0. В противном случае, либо а22 0 (и тогд
меняя местами х и у, получим уравнение, в котором ац 0), ли(
ац = а22 = 0, aj2 0 (и тогда уравнение (1.2.1) уже имеет требуем:
вид). Итак, пусть ац у! 0. Рассмотрим обыкновенное дифференциал
ное уравнение
ац(-у^)2 — 2aJ2-y^ + а22 = 0, (1.2
ах ах
которое называется характеристическим уравнением для исходнс
уравнения (1.2.1). Характеристическое уравнение иногда записыва!
в симметричной форме:
“п№)2 - Zaudydx + a22(di)2 = 0. (1.2.6J
Решения характеристического уравнения называются характерно!
ками. Уравнение (1.2.6) равносильно двум уравнениям
dy _ Д12 ± Уа12 ~ ana22 , .
dx ац '
Согласно определению 1.2.1, знак подкоренного выражения определя
тип уравнения (1.2.1).
Лемма 1.2.1.(Если у>(х,у) = С является общим интегралом ypt
нения (1.2.0), то
anVx + 2ai2<pI<p!( + a22<p2 = 0. (1.2
Доказательство. Пусть (то, уо) - некоторая точка нашей обл
ти. Для определенности считаем, что у>у(хо, Уо) 0- Пусть ф(аго, уо)
Cq. Проведем через точку (го, Уо) интегральную кривую у = у(х, (
уравнения (1.2.6). Ясно, что у(хо,Со) = уо- Так как <р(х,у) = С
10
ельн'1- jr
dy
бший интеграл, то на интегральной кривой <рх + <ру — = 0; следова-
„ !^ = -—. Тогда
4>у
,dy.2 о dy , /„ ,'АЛ2 4>х . х
0 = ац(^) -2“12^ +a22 = (“1l(—) + 2a12—+ a22)|s=IZ(x,Co)
частности, полагая х = х0, получим (1.2.7) в точке (го, Уо) В силу
роизвольности точки (хо,уо) лемма доказана. □
Случай 1: ah — ana22 > 0 - гиперболический тип. В этом случае
ерез каждую точку данной области проходят две различные характе-
истики уравнения (1.2.6). Пусть <р(х,у) = С и ^(г,у) — С - общие
нтегралы уравнений (1.2.6±) соответственно. Эти интегралы опре-
еляют два семейства характеристик уравнения (1.2.6). Отметим, что
ерно (1.2.3), так как
>рх _ 012 + У°12 ~ а11а22 Фх
Vv “и ’ Фу
012 — уа12 ~ а11а22
Оц
Vy^ Фу'
!делаем в (1.2.1) замену (1.2.2), где у>,ф взяты из общих интегралов,
'огда из (1.2.5) и леммы 1.2.1 следует, что ац = а22 = 0. Поэтому
равнение (1.2.4) приводится к виду
и(у = Ф(Сл,и, «£,“»), (1.2.8)
де Ф = — F/(251г). Это и есть канонический вид уравнения (1.2.1) для
иперболического случая. Часто встречается и другой канонический
ид гиперболического уравнения. Положим £ = а + /3, у = а — /3, т.е.
' = (£ + у)/2, р = (£-г))/2. Тогда и( = (ua + u/j)/2, un = (ua-u^)/2,
<1 = (uaa — Дда)/4, и уравнение (1.2.8) принимает вид
Uaa - иРР = $i(a,P,U,Ua,Up). (1.2.9)
>то второй канонический вид гиперболического уравнения. Простей-
шим пРимером гиперболического уравнения является уравнение коле-
>ании струны utt = ихх или более общее уравнение (1.1.4).
Случай 2:|а22 - аца22 = 0 - параболический тип. В этом слу-
ае уравнения (1.2.6+) и (1.2.6-) совпадают, и мы получаем один об-
ций интеграл уравнения (1.2.6): <р(х,у) = С. Сделаем замену (1.2.2),
11
где ф(х,у) - любая гладкая функция, удовлетворяющая (1.2.3). Т<
по лемме 1.2.1 йц = 0. С другой стороны, так как <цз =
то в силу (1.2.5) йц = (^/опУ’х + у'йгг^)2, следовательно, уйщ
у/а&Ру = °- Тогда а12 = (y/aU<px + V'nJii'PsXx/anV’x + y/aHV's) = 1
уравнение (1.2.1) приводится к каноническому виду
“чч = ф(£. ’?>«,“{, «ч), (1-2
где Ф = —F/(2a22) Простейшим примером уравнения параболич<
го типа является уравнение теплопроводности щ = ихх.
Случай 31 ,«12 ~ апа22 < 0 ~ эллиптический тип. В этом сл]
правые части (1.2.6±) комплексны и различны. Пусть <р(х,у) =
комплексный общий интеграл уравнения (1.2.6+). Тогда <р(х,у)\
является общим интегралом уравнения (1.2.6 ). Перейдем к комш
ным переменным, полагая ( — <р(ж,г/), г/ = <р(х,у), т.е. сделав!
мену (1.2.2), где ф = ф. Повторяя рассуждения из случая 1, п
чим, что уравнение (1.2.1) приводится к виду = Ф(£, Т), и, и%
Но это уравнение еще не является каноническим для эллиптиче<
уравнения, так как — комплексные переменные. Перейдем к
щественным переменным (а, 0), полагая £ — а + i/З, т] = а -
Тогда а = (£ + г/)/2, /3 = (£ — ту)/(2г). Имеем = (иа — iu[
и, = (иа + «“4)/2, «£ч = (“аа + “де)/4. Следовательно, уравв
(1.2.1) приводится к каноническому виду
иаа + ирр = Ф1(а,/3,и,иа,ир). (1;
Простейшим примером эллиптического уравнения является двуме
д2и дги
уравнение Лапласа —-г + —= 0.
охг ауг
Замечание 1.2.2. Так как мы выходили в комплексную плоек
то наши рассуждения для случая 3 справедливы лишь тогда, к
а^(х,у) ~ аналитические функции. В общем случае приведение :
ионическому виду для эллиптического уравнения является техни»
более сложной задачей. Мы же для простоты ограничились слу
аналитических коэффициентов.
2 Классификация уравнений второго порядка со многими
зависимыми переменными. Рассмотрим уравнение
V О (г1,....^)^- + ^............................ £-) = 0. (1.2.12)
Z. Qxidxj oxi охп
Ы=1
[е нарушая общности, считаем, что = ац. Сделаем замену пере-
[енных ___
Ct = £t(zi, • •, in), k = I?n. (1.2.13)
% „
(бозначим am := Югда
ди " ди d2u " d2u " ди d2(k , .
chc, t=i Э& dxtdxj t d£kd& i=13ft dx,dxj
[одставляя (1.2.14) в (1.2.12), получаем
E аы(С1, • • •, £n) o, + Ffa,..., C„,U, .... ^-) - 0,
k,i=i d^n
де
Oti = E ‘hjOlikaji.
ij=l
ассмотрим квадратичную форму
E ^ijPiPj, (1.2.15)
i,j=l
це a°ij = a,j(x“,..., s°) в некоторой фиксированной точке (г°,..., х®).
имена
Pi = Е “itgt
1=1
риводит (1.2.15) к виду'
п п п
S й?1«1«|, где a°kl = Е а°^ац.
к,1=1 i,j—l
аким образом, коэффициенты главной части уравнения изменяются
алогично коэффициентам квадратичной формы. Как известно, выбо-
I 9 ™неиного преобразования можно привести квадратичную форму
) к диагональному виду.
£ a?)PiPj = Е(±9|2), т<п,
причем число положительных, отрицательных и нулевых коэффц
ентов при gf инвариантно относительно линейного преобразована
Определение 1.2.2. 1. Если т = п и все коэффициенты при
одного знака, то уравнение (1.2.12) называется эллиптическим в то
хо-
2. Если т = п и все знаки кроме одного одинаковы, то уравнв
(1.2.12) называется гиперболическим в точке хд.
3. Если т = п и число знаков ” + ” и ” — ” более одного, то ypai
ние (1.2.12) называется улътрагиперболическим в точке Хо-
4. Если т < п, то уравнение называется параболическим в то
х0- |
Например, уравнение Лапласа Ди = 0 является эллиптичеа
волновое уравнение ult = Ди - гиперболическим, а уравнение те
проводности и< = Ди - параболическим во всем пространстве.
1.3. ПОСТАНОВКА ЗАДАЧ МАТЕМАТИЧЕСКОЙ
ФИЗИКИ
Из курса обыкновенных дифференциальных уравнений изве<
что решение обыкновенного дифференциального уравнения опре;
ется неоднозначно. Так, решение уравнения n-го порядка у^ — F(
у1,..j/"-1)) зависит, вообще говоря, от п произвольных постоям
у = <р(х, Ср..., Сп). Если мы хотим выделить единственное реш
уравнения, то нужно задавать некоторые дополнительные условия
пример, начальные условия J/|I=IO = од,..., = «п-1
Для уравнений с частными производными решение также оп[
ляется неоднозначно, причем общее решение зависит, вообще гов
от некоторых произвольных функций. Например, рассмотрим сл;
двух независимых переменных (х, у):
уравнение -- — 0 имеет общее решение и — /(х), где /(
оу
произвольная функция от х;
д2и J
2) уравнение = 0 имеет общее решение и = f(x) + g(y)
f(x) и g(t/) - произвольные гладкие функции от х и у соответс
но.
14
того чтобы выделить единственное решение уравнения мате-
й (Ьизики, нужно задать некоторые дополнительные условия
атическои ч*
еизвестную функцию. Поэтому каждая задача математической фи-
1 ставится как задача решения некоторого уравнения с частными
оизводными при определенных дополнительных условиях, которые
гктуются ее физической постановкой. Рассмотрим несколько приме-
>в постановки задач математической физики.
I рассмотрим задачу о колебаниях однородной струны длины I (0 <
< /) закрепленной на концах. В этом случае нужно искать решение
>авнения
иц = а2ихх + f(x, t) (1.3.1)
области Q < x <1, t> 0 при условиях
U|X=O = W|X=/ = 0.
(1.3.2)
акие условия называются граничными или краевыми. Однако задания
(ведения струны на концах недостаточно для описания процесса ко-
‘баний. Необходимо также задать положение струны (и ее скорость)
начальный момент времени t = 0, т.е. наложить начальные условия
“|<=0 = “ф=0 = V’C®)-
(1.3.3)
аким образом, мы приходим к задаче (1.3.1)—(1.3.3) о поиске решения
i,i) уравнения (1.3.1) в области 0 < х < I, t > 0, удовлетворяющего
>аевым условиям (1.3.2) и начальным условиям (1.3.3). Такая задача
1зывается краевой задачей для уравнения колебаний струны. В ли-
ратуре используется также другой термин - смешанная задача, так
1к в ней участвуют как граничные, так и начальные условия.
Отметим, что вместо (1.3.2) можно задавать и другие краевые усло-
1я. Например, если концы струны не закреплены, а движутся по за-
1Нному закону, то краевые условия имеют вид U|I=0 = p(t), U|I=i =
t). где p(t),j/(t) - известные функции. Если, например, конец х = 0
ободно перемещается, то краевое условие в точке х = 0 записыва-
ем в виде Ux|i=q = 0. Если конец х = 0 упруго закреплен, то краевое
ловие имеет вид (цх + /ш). _0 = о. Встречаются и другие виды кра-
“х условий.
Рассмотрим колебания бесконечной струны (—оо < х < оо).
tMo СЛУчае кРаевые условия отсутствуют, но по-прежнему необхо-
адавать начальные условия. Мы приходим к следующей задаче:
15
найти решение и(х, t) уравнения (1.3.1) в области —оо < х < оо, t
удовлетворяющее начальным условиям (1.3.3). Такая задача наз
ется задачей Коши для уравнения колебаний струны.
Аналогичным образом ставятся задачи и в многомерном сд
(х 6 R"), а также для уравнений параболического типа (подро
см. ниже главы 2-3).
III. Приведем пример постановки задачи для уравнения эллипп
кого типа. Уравнения эллиптического типа описывают, как прав
стационарные процессы. Например, задача нахождения стациона
го теплового поля в некоторой области D переменных х = (ац, х:
при условии, что на границе S области D поддерживается фикс
ванная температура 9з(х), записывается в следующем виде:
Д« = 0, U|E = <p, Ди:=£)^2. (1
Задача (1.3.4) называется задачей Дирихле. Другие постановки з г, „
' . , , зшение задачи (1.3.6) имеет вид
математической физики см. в [1-6].
Корректно поставленные задачи. Задача математическо! «n(z,y) = —е '^sinnxshnp, shny := - (e"v — е пу).
зики, состоящая в определении решения дифференциального ура
ния с частными производными при некоторых дополнительных у<
виях, является математической моделью того или иного физичеса
процесса. Решение соответствующей задачи зависит, как мы вид
от некоторых функций, входящих в уравнение, в начальные и 1
ничные условия, которые мы будем называть исходными данни
При исследовании решения задачи математической физики центр]
ное место занимают следующие вопросы: 1) существование реше]
2) его единственность; 3) зависимость решения от ” малых” измена
исходных данных.
Если ’’малые” изменения исходных данных приводят к ’’мал!
изменениям решения, то будем говорить, что решение устойчиво.!
нечно, понятие ’’малое изменение” каждый раз нуждается в точ
определении. Например, для задачи Дирихле (1.3.4) решение будем
зывать устойчивым, если для любого е > 0 существует <5 = <5(е)
такое, что если |^j(x) — ф(т)| < <5 при всех х g Е, то |д(х) — «(т)|
при всех х g D (здесь и, й - решения задачи Дирихле при гранич
условиях р, <р соответственно).
Важным классом задач математической физики, который был вве-
н Адамаром, является класс корректно поставленных задач.
Определение 1.3.1. Задача называется корректно поставленной,
ее решение существует, единственно и устойчиво.
Пример Адамара. Приведем пример некорректно поставленной
Дачи. Рассмотрим задачу Коши для двумерного уравнения Лапласа:
«хх + «уу - °, -°° < х < оо, у > 0, 1
«|у=о = 0, «у|у=о = 0. / V )
здача (1.3.5) имеет единственное решение и = 0. Наряду с (1.3.5)
осмотрим задачу
«хх + «уу = 0, —оо < х < оо, у > 0,
«|у=о = 0, «у|5=о = е-'"'" sin пт, п > 1.
(1.3.6)
ак как | sinnx| < 1, то Iе ^”sinnx| —> 0 при п —> оо. Однако
шение и„(х, у) может неограниченно возрастать. В самом деле,
«п(г,3/) = enB"'/5isinnx(l - е-2п“).
2п '
ри каждом фиксированном у > 0
max |«п(т, у)| —> оо
-оо<х<оо 1
16
(С С МД
I YW YjvCl
ко
С X ЧФ -Д 6
C * 51 ГЛАВА 2
Уравнения гиперболического типа возникают в основном при
делировании колебательных процессов. В этой главе мы изучаем
новные задачи для уравнений гиперболического типа и представ!
методы их решения.
2.1. ЗАДАЧА КОШИ ДЛЯ УРАВНЕНИЯ КО.
1. Рассмотрим задачу Коши i t? t
“i|i=O = V’(i)-
faecb x, t - независимые переменные (сфизи^^кой^гочки зрен
.переменная. t. 7 Bpe^)j неизв<Ж
гганта. Bj^m предполагать, что ^\т)непр(
проста
ция, с
дифференцируема (тр е С1)^ (р(х) дважды непрерывно диффере:
Определение 2.1.1. Функция u(x,t) называется решением за<
\Коши (2.1.1)—(2.1.2), если u(x,t) G C2(D), и u(x,t) удовлетвр
^**^>(2.1.1) и (2.1.2). Здесь и в дальнейшем запись и & Ст(р} означ
что функция и имеет в области D непрерывные частные произвол
до m-го порядка включительно.
-----Найдем сначала общее решение уравнения (2.1.1). Характера
ческое уравнение для (2.1.1) имеет вид — = ±а (см. п. 1.2), слеп
тельно, уравнение (2.1.1) имеет два семейства характеристик т±
const. Замена переменных
— х + at, г) = х — at,
„ривоДиГуравнение (2.1Л) к виду и( =
этого уравнения суть и - /(£) + д(т)), т<
(2.1.1) имеет вид
= 0. Таи как общее решение
, то общее решение уравнения
и(т,г) = /(т4-а4)+з(х-а«), | #Д.З) я
де f и g - произвольные дважды непрерывно дифференцируемые
ia общее решение уравнения (2.1.1) вида (2.1.3). Предположим, что
решение задачи (2.1.1)-(2.1.2) существует. Тогда оно имеет вид (2.1.3)
зри (21'2)’получаем
, f(x)+g(x) = <p(x),
или
J f(x) + g(x) = ч>(х),
af'(x) - ag'(x) = tf>(®)
1 гх
/(®) - ff(®) = - / ^(s) ds
О, dXQ
Разрешая эту систему относительно f и д, находим
/(х) = | ф(х) + £ V'(s) ds, д{х) = i <р(х) V>(«) ds.
Подставляя эти выражения в (2.1.3), приходим к формуле
u(x,t) = i(<p(x + at) + ¥>(х - at)) + f*+a‘i/>(s)ds, (2.1.4)
Kt называете^ формулой Даламбера. Таким образом, мы показа-
тели решение задачи (2.1.1)—(2.1.2) существует, то оно дается
формулой (2.1.4)ДОтсюда сразу следует единственность решения зада-
1И (2.1.1)-(2.1.2). С другой стороны, непосредственной проверкой не-
рудно убедиться, что функция u(x,t), задаваемая^^рмулой (2.1.4),
действительно дает решение задачи (2Л.1)-(2.1.2). Таким образом, мы
(оказали следующее утверждение.
„Г”- 2-Ы. Пусть <р(х) g С2, ф(х) С С1. Тогда решение задо-
и оши (2.1.!)-(!>.1.2) существует, единственно и дается формулой
тламбера (2.1. J).
следзем теперь устойчивость решения задачи (2.1.1)-(2.1.2).
1гт^~*"^в'”еНИе Решение задачи (2.1.1)-(2.1.2) называется
i ели для любых е > 0, Т > 0 существует S = 6(е, Т) > О
19
такое, что если |<^(х) — <р(х)| < й, |^>(х) — V'(x)| — <5 ПРИ всех J
|u(x, t)—й(х,<)| <е при всех (x,t) G D. Здесь u(a:,t) - решение aaj
Коши при начальных данных <р, ф.
Покажем, что решение задачи Коши (2.1.1)-(2.1.2) устойчиво в <
ле определения 2.1.2. В самом деле, в силу (2.1.4) имеем
1 1 |
|u(x, t) — й(х, t)| < - |^(х + at) — ф{х + at)| + - |<р(х — at) — ф(х -
Ла * ® Gt
Для е > 0, Т > 0 выберем <5 = е/(1 + Т) и предположим, что |<р|
£(®)1 < 6, |^(г) - 0(г)| < <5 при всех из. (2.1.5) нах
|u(x,t)-u(x,t)| < <5/2+5/2+(6/2a)2af = (1+Т)й = s. Т&им обй
задача Коши (2.1.1)-(2.1.2) корректно поставлена.
В дальнейшем нам потребуется следующее утверждение.
Лемма 2.1.1}
1. Если <р(х) и ф(х) - нечетные функции, то u(0,t) = 0.
2. Если <р(х) и ф(х) - четные функции, то ux(0, t) = 0.
Доказательство. Пусть <р(х) = — <р(—х), ф(х) = -ф(—х). 1
из (2.1.4) при х = 0 имеем Q'
1 1
“(°> <) = + V’(-aO) + V'(s) ds
Аналогично доказывается и второе утверждение леммы.
at at
Рис. 2.1.1
Замечание 2.1.1/ Решения вида и = f(x + at) и и = д(х
называются бегущими волнами, причем а - скорость распростри
волны. Волна и = f(x + at) перемещается влево, а волна и = д(,х
V0U.
. 91 Й Таким образом, общее решение (2.1.3) уравнения
вправо (рис. i i
1.1) является^ммойбегщо^^ . д
ассмотрим задачу Койти для^неоднородного уравнения колеба-
[ИЙ струны:
иц
= а2ихх + f(x,t),
U|i=o = ¥’(г). =
ешение задачи (2.1.6) ищем в виде u(z,t) = u(x,t)+v(x,t), где u(z,t)
решение задачи Коши для однородного уравнения:
йи = а2йхх, -оо < х < оо, t > 0; й|(=0 = V>(x), “ф=0 = V’fa). (2-17)
v(i,t) - решение задачи для неоднородного уравнения с нулевыми
ачальными условиями:
Vtt = a2vxx + f(x,t), -оо < X < оо, t > 0; пц=0 = 0, vt|f=0 = 0. (2.1,8)*
етрудно убедиться, что решение задачи (2.1.8) имеет вид (У^)
v(х, t) = I w(x, t, т) dr, (2.1.9) ’
ie функция w(x,t, т) при каждом фиксированном т > 0 является
штением следующей вспомогательной задачи:
wlt = a2wxx, -оо < х < оо, t > т; W|t=T = 0, wt|(=T = f(x, т). (2.1.10)
:щение задачи (2.1.10) выписывается с помощью формулы Даламбе-
, w(x,t,T) = lf+;(‘-;>/(s,r)(is.
, 2а '
ак как решение задачи (2.1.7) также выписывается с помощью фор-
^лы аламбера, то мы приходим к следующему утверждению.
Теопема *> i о гг хлС'Т)
с УУ&пь с 0(D) и f имеет непрерывную
с те^Ю пР°изв°дную fx(x,t), и пусть удовлетворяют услови-
МЫ ?огда решение задачи Коши (2.1.6) существует,
инственно и дается формулой
, ^Х' ~ + а4) + v(x — at)) + 1- [ i/i(s)ds~\- -
I 2a Jx—at
21
, 1 r‘ . rz+^t-r) ,, . , я
+— / dr / , , f(s,r)ds. M
2a JO Jz-a(t-r) J
Замечание 2.1.2. Из формулы (2.1.11) видно, что решение и|
в некоторой фиксированной точке (х, t) зависит от исходных дар
заданных только в характеристическом треугольнике Д(х,
{(s,r) : х — a(t — т) < s < х + a(t — т), 0 < т < t} (рис. 2.1.2). j
гими словами, изменение функций вне треугольника Д(х,|
изменяет значения решения u(x,t) в данной точке (x,i). Физия
это связано с тем, что скорость распространения волн конечна. 1
3. Рассмотрим теперь задачу для полубесконечной струны:
u(j = а2ихх, х > 0, t > О,
“|х=о = 0. “|<=о = <р(я)> «ф=0 = ^(®)-
Для решения данной задачи применим та^ назывгшмьд) МуТЖ<
жения. Пусть V’(s) 6 С1, <р(х) 6 С2, <р(б) = ^”(0)=V(0) = Д.
должим у>(х) и ^(x) на всю ось -оо < х < оо нечетным образо!
рассмотрим функции
(*’
Ф(х) =
-V'(-x),
Согласно формуле Даламбера функция
u(x,t) = ^(Ф(* + at) + Ф(х - ai)) + f * $(s)ds
2 Zu •Jx—at,
является решением задачи Коши
кк = а2ихх, -оо < х < оо, t > 0; U|(=o = Ф(я), “t|t=o =
22
лемме 2.1.1 имеем U|x=o = Рассмотрим функцию (2.1.13) при
> о, t > 0. Тогда
U|1=o = Ф(г) = ««|«=0 = ОД = ^(®)> х-°-
КИМ образом, формула (2.1.13) дает решение задачи (21.12). При^оЛ
> 0, t > о запишем (2.1.13) в следующем виде: И
-(у>(х + at) + »»(х - at)) + V"(s) ds, x > at,
u(x, t) = ' . 1 rx+at
l^x + at) - <p(flt - x)) + - /at x V-(s) ds, x <
2 MM
'МСТИМ, что в области х > at влияние отраженной от конца х = 0
лны еще не сказывается, и решение имеет тот же вид, что и для
сконечной струны.
(Замечание 2.1.3. f Аналогичным образом решается задача с крае-
1м условием иф=о = 0. В этом случае, согласно лемме 2.1.1, нужно
рдолжать <р(х) и il>(x) четным образом.
. „ 1ЛЕБАНИЙ СТРУНЫ
ДлГ, чЩхуГ
• Рассмотрим следующую краевую задачу:
Дд — а иХХ1 0 < х < I, t > 0,
(2.2.1)
(2.2.2)
(2.2.3)
u|z=0 = U|z=l = 0,
“|«=0 = <Р(х), и,|,=0 = V’(s)^
Ьзначим D = {(х,«) ; о < х < I, t > 0}.
Функция u(x,t) называется решением задачи
| 3) 2' если и u(x,t) удовлетворяет (2.2.1)-
Ь нет и ИССЛедУем слеДУюЩую вспомогательную задачу. ByjieS ис-
Гия vf ииальиы« (т.е. не равные тождественно нулю) частные ре-
ия уравнения Го о г
Рпускаю! Z 1 ’ УдовлетвоРяк>Щие краевым условиям (2.2.2)
ие разделение переменных, т.е. представимые в виде
u(x,t) = y(x)'T(t). ЧМЙ1
23
0кЛOp- \
-\aaaa -Ь
Для нахождения таких решений подставим (2.2.4) в (2.2.1) и (2
силу (2.2.1) имеем Y(x)T(t) = a2Y"(x)T(t). Здесь и в дальней]
обозначает дифференцирование по х, а ” • ” - дифференцировав
„ Y"(x) T(t) m
Разделяя переменные, получаем . -г- = ^ГГ1. . . 1 ак как х и
Y(x) arT(t)
зависимые переменные, то равенство возможно лишь тогда, kq
дроби равны константе; эту константу обозначим —Л. Следов!
Y"(x) At)
У(т) а2Т(4)
Кроме того, краевые условия (2.2.2) дают
Из (2.2.5) и (2.2.6) выводим
f (t) + a2XT(t) = 0,
Y"(x) + XY(x) = 0,
y(o) = y(0 = o.
Таким образом, функция T(t) является решением обыкновещ
ференциального уравнения (2.2.7), а функция У(т) является
ем краевой задачи (2.2.8)-(2.2.9), которая называется краевой
Штпурма-Лиувилля. Нас интересуют нетривиальные решеня
(2.2.8)—(2.2.9), но они существуют лишь при некоторых значв
раметра А.
Определение 2.2.2. Те значения А, при которых краев!
(2.2.8)-(2.2.9) имеет нетривиальные решения, называются с<
ными значениями, а соответствующие нетривиальные решев
ваются собственными функциями. Множество собственных !
называется спектром.
) Ясно, что собственные функции определяются с точност!
стоянного множителя, так как если У(т) - нетривиальное
задачи (2.2.8)—(2.2.9) при А = А* (т.е. У(х) - собственная’
соответствующая собственному значению А*), то и функци
С - const, также является решением задачи (2.2.8)-(2.2.9) пр
А = А*.
24
ветвенные значения и собственные функции краевой за-
Найдем сПусть Д = р2. Общее решение уравнения (2.2.8)
^каждом фиксированном А имеет вид
sinpx , „
УГт) = А—-—Ь Bcospx . „
(> р->
частности, при р = 0 (2.2.10) дает У(хУ= Аг + В). Так
то в (2.2.10) в = 0, т.е.
, . л sin рх
= Ар'
(2.2.11)
рставляя (2.2.11) В краевое условие Y(l) = 0, получаем А^у^ = 0.
с как мы ищем нетривиальные решения, то А^О, и мы приходим
равнению на собственные значения:
sin pi _ 0
Р
, \ п1г
эни уравнения (2.2.12) суть рп = —, следовательно, собственные
чения задачи (2.2.8)-(2.2.9) имеют вид
(2.2.12)|
силу (2.2.11), суть
18
ответствующие им собственные функциг
У„(х) = вшух, п = 1,2,3,... (2.2.14)
эчностью до постоянного множителя). Так как нетривиальные ре-
ия задачи (2.2.8)-(2.2.9) существуют лишь при Л — А„ вида (2.2.13),
равнение (2.2.7) имеет смысл рассматривать лишь при А = Ап:
fn(t) + a2AnTn(t) = 0.
лее решение этого уравнения имеет вид
Tn(«) = A„sm—t + BnCOS^f> п_ 12,3,..., (2.2.15)
ЗА)1 по Пр°извольные постоянные Подставляя (2.2.14) 0^2.15)
лучаем все решения вспомогательной задачи:
i3-» = (A sin f I D йП7Г птг
/ f+ B„cos-pt)sinyI, n = 1,2,3,... (2.2.16)'
25
Решения вида (2.2.16) называются стоячими волнами.
Будем искать решение краевой задачи (2.2.1)-(2.2.3) в виде
стоячих волн: 0<Г iX'LOxJ Ъф'Х?
u(x,t) = > . (А ——t)sin— х.
n=l * ‘
Проведем пока формальные рассуждения, не заботясь о сход
ряда. По построению функция u(x,t) удовлетворяет уравнению (
и краевым условиям (2.2.2) при любых Ап и Вп. Выберем А„
так, чтобы и(х, t) удовлетворяла и начальным условиям (2.2.3). С
целью, подставляя (2.2.17) в (2.2.3), находим
оо „я- iA . оо ап7Т п7г
у(х) = Е ВП sin —х, ф(х) = > . —7— А„ sin —х.
n=I J
Пользуясь формулами для коэффициентов Фурье [7, гл. 1], вычи
_ 2 fl , . . птг , , 2 А . . . nir я
Вп = f „V>(x)sm—xdx, Ап =-------- Mx)sm—xdx. Д
1 jo I an-п I
Покажем теперь, что функция u(x,t), определяемая соотноше
(2.2.17)—(2.2.18), является решением задачи (2.2.1)-(2.2.3). Для
надо наложить некоторые условия на ip(x) и
Теорема 2.2.1; Пусть ip(x) е С3[0, Z], ф(х) V4 C^jOJ], 4
^>(1) = </'(0) = ¥>"(/) = ф(0) = V'G) = 0- Тогда функция и{х, t),
деляемая равенствами (2.2.17)-(2.2.18), является решением з
(2.2.1)-(2.2.3).
Доказательство. Из теории рядов Фурье [7, гл. 1] известны i
ющие факты. м
1. Система функций полна и ортогональна в Ц
2. Пусть f(x) е £г(0, /) и
“n = 7^,^sinTI<ia:’ = 7Z •^I)cosTldl
- коэффициенты Фурье функции f(x). Тогда
Кроме того, отметим, что
если числа7п таковы, что 52 1?п|2 < оо, то J2 |?п|/^ < оо.
П X. п
(
4 26'
L является следствием неравенства |7п|/п < (п 2 + |7’|)/2.
Покажем, что
£ п2(|Лп| + 1^1) < <*>
(2.2.21)*
интегрируя по частям интегралы в (2.2.18), вычисляем
птг . , 2l2 fl „ ,.птг
cos dx, А„ =-----------3-3 ф (х) sin —X dx.
I атг3пл Jo ' ' l
используя утверждение (2.2.20) и свойства коэффициентов
самом деле,
рЬ^(2.2.19), приходим к (2.2.21). Из (2.2.21) следует, что функция
с () определяемая равенствами (2.2.17)-(2.2.18), имеет в области D
1рерывные частные производные до второго порядка включительно
е u(x,t) £ С2(Р)), причем эти производные могут быть получены
пенным дифференцированием ряда (2.2.17). Очевидно, что u(x,t)
1яется решением уравнения (2.2.1) и удовлетворяет краевым услови-
(2.2.2). Покажем, что и(х,2) удовлетворяет и начальным условиям
2.3). В самом деле, из (2.2.17) при t = 0 имеем j
00 717Г 00 CLTVTI Т).7Г
u(x, 0) = £ Вп sin —х, u((x, 0) = £ Ап sin ^x.
n=l * n=l 1 ‘
шовательно, <. k ~ ^9* Tc
~ 7 L u(z, 0) sin — x dx, An=--- f u((x, 0) sin ^xdx.
• JU I ann Jo I
1внивая эти соотношения с (2.2.18), получаем при всех п > 1 :
(и(х, 0) - v(x)) sin ™х dx = 0, 0) - V-(x)) sin у x dx = 0.
илу полноты системы функций { sin в Г2(0, /) заключаем,
и'х’ 0) — <р(х), ц((х, 0) = ^(х), т.е. функция и(х,2) удовлетворяет
ильным условиям (2.2.3). Теорема 2.2.1 доказана. □
ргииДИНСТВеНН°СТЬ решения «Раевой задачи. Интеграл
Терема 2.2.2. £сли р=-------
™вуеТп‘ то оно единственно.
решение краевой задачи (2.2.1)-(2.2.3) су-
>енно.
•3). Обти=," “усть иЦМ) И ц2(х,1)- решения задачи (2.2.1)—
чнм u(x,t) - ujjx, t) _ u2(x,t). Тогда u(x,t) € C’(P)
U|I=0 = U|X=I = 0, U|1=o = Иф=0 = 0. (2.2.22f
27
Рассмотрим интеграл
K(t) = I Г (i(u,(l’t))2 + ‘))2) dx'
который называется интегралом энергии. Тогда
£(*) ~ <)) + их(х, t))uxt(x, t)\dx.
“ TiZ eZMZlA
Интегрируя по частям второе слагвеячое, получаем
E(t) = l‘ ut(x, t)(^ut((x, t) - u„(z, t)) dx + loOtrfs, t)ut(x, t)
В силу (2.2.22) это дает E(t) = 0 и, следовательно, E(t) = coni
как
Я(0) = 2 £ (?(u,(l’0))2 + 0))2) dx = °’
то E(t) = 0. Отсюда заключаем, что U|(x,t) = 0, ux(a:,t) =
u(i,i) = const. Снова используя (2.2.22), получаем u(x,t) = 0,
рема 2.2.2 доказана.
3. Метод Крылова. Рассмотрим краевую задачу для neoi
кого уравнения колебаний струны:
utl = а2ихх + f(x, t), 0 < х < I, t > 0, 1 J
“|i=o = “|z=l = 0, U|t=o = u(|t=0 = ф(х'). J
Пусть функция f(x,t) непрерывна в D, дважды непрерывно
ренцируема по х и /(0, t) = /(!,t) = 0. Пусть функции <р(х)
удовлетворяют условиям теоремы 2.2.1. Разложим функцию /(
ряд Фурье по переменной х:
ОО 77 7Г 2 fl T1/7V
f(x, t) = E /„(t) sin —x, f„(t) = -10 f(x, t) sin — x <fc.
Будем искать решение задачи (2.2.23) в виде ряда по собст!
функциям задачи Штурма-Лиувилля (2.2.8)-(2.2.9):
, . “ птг
“(г10 ~ Е “n(0 Sln ~ГХ-
п=1 1
Подставляя (2.2.24) в (2.2.23) и приравнивая коэффициенты пр!
циях sin ™х, п = 1,2,..., получаем соотношения для опрв
коэффициентов un(t) :
^п(0 + (^y-)2Wn(t) = /п(0> un(0) = Вп, йп(0) = Ап,
28
д и Дп определяются по (2.2.18). Задача Коши (2.2.25)
1е числа „ пешение, которое может быть найдено, например,
«ррт единственн и
Хм вариации произвольных постоянных:
/77? 7Г QYWT I [t - , \ CLTITV . .
u„(t) = А„ sin + В„ cos —< + — /0 Ш «« ~(* - т) dr.
Unl ’ ' (2.2.26)
дно убедиться, что функция u(x,t), определяемая равенствами
°2 24) и (2.2.26), является решением задачи (2.2.23).
2.3. ЗАДАЧА ГУРСА
^Задача Гурса заключается в решении гиперболического уравнения
данными на характеристиках. Поэтому эта задача называется так-
е задачей на характеристиках. Нам удобнее будет рассматривать
>т канонический вид гиперболического уравнения, для которого ха-
1ктеристики параллельны осям координат. При этом для простоты
раничимся линейным уравнением. Рассмотрим следующую задачу
грса: , \
(ȣ)
их„ = а(х, у)их + Ъ(х, у)иу + с(х, у\и 4- f(x, у), (г, у) G П,
= = V>(x).
(2.3.1)
«|х=хо = ¥>(!/), и|ц=ю = V'(z). (2.3.2)
[есь т, у - независимые переменные, и(х,у) — неизвестная функция,
= {(*.У) хо х < Xi, уд < у < г/]} - прямоугольник. Уравнение
3.1) имеет два семейства характеристик х = const и у = const.
1ким образом, условия (2.3.2) представляют собой условия на харак-
ристиках х = х0 и у = у0.
Определение 2.3.1. Функция и(х,у) называется решением задачи
(2.3.2), если и(х,у) определена и непрерывна вместе со своими
’2ИзВ2)ДНЫМИ Ux,Uv,Uzv в прямоугольнике П и удовлетворяет (2.3.1)
и и МЯ ^Усть Функции a,b,c,f непрерывны в, П, а Функ- ,
существует и единственно.
/твой системе”"0 ^ведем задачУ Гурса (2.3.1)-(2.3.2) к эквива-
е и(г> О задачиИ/2ТчГ^аЛЬНЫХ ^’Равнений- Предположим, что реше-
1)-(2.3.2) существует. Положим v = их, w = uv.
29
Тогда
vy = av + bw + си + f, wx = av + bw + cu + f,'o'uy = w. (2
Интегрируя (2.3.3) и учитывая условия (2.3.2), получаем
v(i,j) = ^'(г) + I (av + bw + си + f)(x,y) dy,
w(x,y) = ¥>'(i/) + (av + bw + cu + f)(£,y)d£, (2
u(x, у) = ф(х) + [" w(x, y) dy.
'1/0
Таким образом, если и - решение задачи (2.3.1)-(2.3.2), то тр
функций и, v,w является решением системы (2.3.4). Справедля
обратное утверждение. В самом деле, пусть тройка и, v,w явл1
решением системы (2.3.4). Дифференцируя (2.3.4), получаем, что
ют место равенства (2.3.3) и, кроме того,
их(х, у) = ^'(х) 4- [У wx(x, ту) dr] =
JVo
•ф'(х) + [y(av + bw + cu + f)(x, у) dy = v(x, у).
'1/0
Таким образом, их = и. Вместе с (2.3.3) это дает, что функция и
является решением уравнения (2.3.1). Далее, из (2.3.4) при х =
У = Уо вычисляем
и(х,уа) = i/>(x), ,
п \1Ъ 2 4 Ч I
и(х0, у) = ^(хо) + / ui(x0, y)dy = V'(io) +
Jyo
т.е. u(x,y) удовлетворяет условиям (2.3.2). Таким образом, за
Гурса (2.3.1)-(2.3.2) равносильна системе (2.3.4).
2. Будем решать систему (2.3.4) методом последовательных пр1
жений. Положим
vo(x,y) = /(х) + l^f(x,y)dy, w0(x,y) = ip'(y) + Д /(€, v) 4
ио(х,у) = ф(х),
vn+i(x, у) = [\avn + bwn + cun)(a:, y) dy,
w„+i(t, y) = F(avn + bwn + cu„)(e, y) d^,
JlQ
Un+i(x,y) = l^wn(x,y)dy.
¥’(У)>
Jyo
tf.'UL^-^-
n!
нстанты M > 0 и К > 1 так, чтобы |u0|, |v0|, |w0| < M,
Выберем . Используя (2.3.5), методом математической индукции
получаем оценки:
S *** <2М>'
в силу оценок (2.3.6) ряды
и(т,у)= £“п(г>»). v(®.»)= £>-(*, у)> w(x,y) = Е^х,»)
П=°(^>) х "=0 (2-3-7)
сходятся абсолютно и равномерно в П (так как они мажорируются
,, г,„ (ц + У1 - х0 - Уо)" ч
сходящимся числовым рядом М 2_,q ^1 7» причем
|u(a,l/)|,|w(*I«)l>l«(a:>v)l Afexp(K(xi + yi -х0-у0)). (2.3.8)
Непосредственной проверкой убеждаемся, что функции и, v, w, постро-
енные по формулам (2.3.7), являются решением системы (2.3.4).
3. Докажем единственность. Пусть тройки (u,v,w) и (й, й,ш) явля-
ется решениями системы (2.3.4). Тогда функции и* = и—й, v* — v—v,
iu* = w — w являются решениями однородной системы
п*(;с>у) = Л (ап* + bw* + cu*)(x,r/) dr],
JVo
w’(x, у) = / (ап* + bw* + cu*)(£, у) d£,
JXq
u'(x,y) = [y w*(x,Tf) dr],
JVo
Гак как функции u*,v*,w* непрерывные П, то существует константа
1 > 0 такая, что |и*|, |и*|, |ш*| < М\. Повторяя предыдущие рассуж-
(ения, по индукции получаем оценку
lu‘(®.V)|,|w*(x,!,)|1|u«(a;,J,)| < MlKn (х + у - хо ~ Уо)"
1РИс” °° ЭТ° “*(Х’V) = W’(l’V) = W'(l’= °'
едуем теперь устойчивость решения задачи Гурса.
пРеделение 2 3 9 о.. п
>ьии, если для гешение задачи Гурса называется устойчи-
ели 1Сл(‘') (,л Л-1?„^?ГО £ > 0 существует б = 5(e) > 0 такое, что
(у) - ^)(у)| < 6) |v>M(x) _ ^)(i)| „ = 011 х 6
/’к)
n
31
[xo.ij, У 6 [1Л)>»1Ь то |u(x,y) - й(х,у)| < e, |ux(x,у) - йх(х,у}
|и„(х,у) - йу(х, у)| < е при всех (х,у) е П. Здесь й(х,у) - ре]
задачи Гурса при исходных данных ср, ip.
Покажем, что решение задачи (2.3.1)-(2.3.2) устойчиво в с
определения 2.3.2. В самом деле, обозначим и' = и — й, н* =
w*. — w — w, и>’ = <р — ib' = ib — ib, f' == Ол В силу (2.3.8), '
/ |в'(х, з/)|, |ш*(х, у)|, |u*(x, у)\ < М' exp(K(xi + yi - ха - дЗ
где М' = max(max|V'*(x)|,max|V’*’,(x)|,max|<p*,'(y)|) < <5. Bi
<5 = еехр(-К(х1 + у\ - хд — уо)) и учитывая, что М’ < <5,
чаем |и*(х, у)|, |н*(х,у)|, |w*(x, у)| < е. Таким образом, задачу
(2.3.1)—(2.3.2) корректно поставлена. Q VbC A Cui/ Д|"XJ J
Переформулируем теперь полученные выше результаты для
го канонического вида гиперболического уравнения. Рассмотри]
чу Гурса
J»- uII-titi+a(x,i)uI+b(x,t)u1+c(x,t)u = (x,t) е Д(хо><о)>1
u(x, х — хо + to) = ф(х), “(т, -х + хо + to) = V'(a:)>
где Д(хо,1о) = {(х, t)) : t - to + х0 < х < —t + to + хо, 0 < t <
характеристический треугольник (рис. 2.4.1). Уравнение (2.3.9)
два семейства характеристик х + t = const и х — t = const
образом, условия (2.3.10) представляют собой условия на хара
тиках Д : t = x-x0+t0 и I2 : t = -x + io+to- Задачу (2.3.9)-
можно свести к задаче вида (2.3.1)-(2.3.2) заменой переменны!
рот осей координат):
. г । „ , t + х хо — tg t — X ]
x = (-T] + to, t = C + tj-xo О C = —, r) = —2~ +
В самом деле, обозначим й({, rf) = и(£—y+to,i+r]-xo) = u(x,t)
Ux = (й( - й,)/2, ut = (й{ + йп)/2, ихх - ии = -й^, следом
уравнение (2.3.9) принимает вид
Чч = “(?> т')й + */)>
где a=(a + b)/2, Ь=-(а-Ь)/2, с = с, f = Характера!
и 12 уравнения (2.3.9) переходят в характеристики уравнения
32
( соответственно, а условия (2.3.10) принимают вид
= to " .(^ = ^(ч), где £(£) = ¥>(f) > V'(’)) = V’(-’7+io+to)-
(€, to) - W ’ TeoDeMa является следствием теоремы 2.3.1.
[оэтому следующая н
2 3.2. Пусть функции a,b,c,f непрерывны в Д(хоЛо), а
ф непрерывно дифференцируемы и <р(х0) = ф(х0'). Тогда
задачи Гурса (2.3.9)-(2.3.10) существует и единственно.
2.4. МЕТОД РИМАНА
Рассмотрим следующую задачу Коши:
- и- + а(т, «)“х + ь(®, г)“< + cfo t)u = / (г, t), -оо < т < оо, t > 0,
(2-4.1)
U|<=0 = ¥’(х), “I|t=o = 0(*)- (2.4.2)
бозначим D = {(т,4) : < x < oo, t > 0}. Будем предполагать,
го a,be Cl(D), C,f e C(P), <p e C2(R), ф G C’(R).
Определение 2.4.1. Функция u(x,t) называется решением задачи
.4.1)-(2.4.2), если u(x,t) G С2(Д) и удовлетворяет (2.4.1) и (2.4.2).
Вывод формулы Римана. Обозначим
и — u„ — uel + аих + but + cu, C'v = vxx — о(1 — (ан)х — + со.
t M(x0,t0)
Рис. 2.4.1
Фиксируем точку \ ,
>льник ДСт t i г, ° И РассмотРим характеристический тре-
; (Рие 2 41Л в{(1,<) : ‘ ~ < х < -t + 10 + 10, 0 <
с° + «о,0). Гра„ица ^ШИНами в точках Wo,to), Р(х0 - i0,0) и
' треугольника А состоит из трех отрезков:
33
оД = А U Ig U 1з, где 11 = МР с уравнением t — х — хо 4- to, 1г
с уравнением t = -х + х0 +tg и I3 = PQ с уравнением t = 0.
Предположим, что решение задачи (2.4.1)-(2.4.2) существуй
значим его и(х, t). Пусть v(x, t) g C2(D) - некоторая функод
как
Jfr v£u — u£’v = (yux — uvx + auv)x — (yut — uvt — buv)t,
то в силу формулы Грина [8, гл. 16] получаем
u£*v) dxdt = fg^(vut~uvt—buv) dx+(yux—uvx+auv) dt (
(с обходом <ЭД против часовой стрелки). Преобразуем интен
границе. Имеем
На Ii: t = х — xg + tg, dt = dx. Обозначим
ai(x) := u(x,x - xg + tg), @i(x) := v(x,x - x0 +10). j
Тогда ai(x) = (ux + 0'i(x) = (vx + Vi)|t=I_Io+to, c;
тельно,
(tifu; + ux) — u(yt + vz) + (a - b}uv)(x, X - Xg + tg) (
= (<3i(®)ai(i) - + (a - ь)(г, x - xq + to)ai(xj/3i(x)
= 1м(/31(1)“1(1)) - “i(a:)(2/3l(I) - (a - i)(x,x - Xo + to)0i(x]
Наложим первое условие на функцию v(x,t), а именно, пот;
2/3] (х) - (а - Ь)(х, х - х0 + to)/?i(x) = О.^о^} >
Решая это обыкновенное дифференциальное уравнение, вычи<
1 гХ
v(x, х - Хд + to) = exp (- / (a - (>)(£, f - х0 + to) <%).
2 Jx0
Тогда
=u(P)v(P)-u(M)v(M).
2. На Ig: t = —x + Xo + to, dt = — dx. Обозначим
a2(®) ~ «(x, -x + io + to), Дг(х) := v(x, -x + x0 + to)-
34
> = (их - т)|«=-х+хо+*» ^х~> = (Wl ~ Vi)|i=-x+xo+io. следова-
Гогда ot2\x) '
кльно,
у. _ jM („(U1 - Ul) - u{yt - vx~) - (a + b)uu)(x, -x + x0 + to) dx
- JM Pi^a^x) + a2(x)P'2(x) - (a + 6)(x, -x + x0 + ta}a2(x)fa(xy) dx
; -|"(/32(x)«2(i)) + !q a2(x')(2/3'2(x')-(a + b'){x,-x + x0 + t0')l32(x)')dx.
[аложим второе условие на функцию v(x,t), а именно, потребуем,
тобы
2/%(x) ~(a + b)(x< ~x + xo + to)ft(i) = 0.
ешая это дифференциальное уравнение, вычисляем
v(x, —х
'огда
+ х0 + t0) = exp (i Г(а + 6)(f, -f + x0 + to) <*£)
Z Jxo
(2.4.6)
(2-4.7)
Д = “(<ЭМ<2) - н(Л4)ы(М).
3. На 7з: t = 0, dt = 0 и, следовательно,
Д = L°-i° (V(x< °MZ) “ v»(®> 0)v(®) - Ь(®, 0)v(a:, 0)p(x)) dx. (2.4.8)
аложим третье условие на функцию v(x,t), а именно, потребуем,
гобы
£’и = 0, (1,<)еД. (2.4.9)
□отношения (2.4.9), (2.4.4), (2.4.6) - это задача Гурса. Ее решение -
ЦИЯ п' сУществует и единственна. Она называется функцией
висит л ™стим’ что "(М) = v(x,t;x0,t0), т.е. функция Римана
4 31 И J Т0ЧКИ Подставляя (2.4.5), (2.4.7), (2.4.8), (2.4.9) в
разрешая относительно u(x0,t0), получаем
о) 2^(Io + fc|)w(Io + to,O) + ^(xo-to)v(xo-to,O))
+1 [г°+1°
2 Ло-(о Wx, 0№(х) _ и((г| 0)^ _ 6(Х| 0^(Х) 0) у>(х)) dx
_ /'<0 rxo+io-l
(2.4.10)
CV((2)
^(^1 f’u°
Формула (2.4.10) называется формулой Римана. Таким обрад
доказали, что если решение задачи (2.4.1)-(2.4.2) существует
дается формулой (2.4.10). Отсюда сразу следует единственна
шения задачи (2.4.1)-(2.4.2). Можно показать [3, гл. 5], что ф
и, определяемая формулой (2.4.10), действительно является ре:
х задачи (2.4.1)-(2.4.2). Отметим, что существование решения
(2.4.1)-(2.4.2) можно также доказать независимо, используя ме(
ложенный в п. 2.3.
В некоторых частных случаях функция Римана может быть
лена явно. Приведем два примера.
Пример 2.4.1. Рассмотрим задачу Коши для неоднородной
нения колебания струны: а
«хх - -оо < х < оо, t > 0, 1
«|1=0 = (P(x), ut|i=o = '0(i)- J
Это частный случай задачи (2.4.1)-(2.4.2), когда а = Ь = с = О
видно, что в этом случае v(x,t) = 1, и формула Римана (2.4.10
нимает вид
“(М) = |(¥’(i + t) + cp(x-t)) + | r+til>(s)ds-^l‘ dr J f(s
I £ JU JX— t+T . (
Сравните (2.4.12) c (2.1.11)!
Пример 2.4.2. Рассмотрим задачу Коши (2.4.1)-(2.4.2) для
нения с постоянными коэффицинтами, т.е. в случае когда а,
константы. Отметим, что здесь, в частности, содержится теп"
ное уравнение (см. п. 1.1). Не нарушая общности, можно счита’
= Ъ = 0 (этого можно добиться заменой u(x, t) = й(х, t) ехр(-
bt/2)). Таким образом, рассмотрим следующую задачу Коши: j
Нхх - ии + си = f(x, t), —оо < х < оо, );'‘с
“|i=o = «i|i=o = V’(i)- J 1
Для функции Римана имеем задачу Гурса:
vxx — Utt + СП = 0, П||=х-хо+1о = 1> H|t=-x+x0+fo = 1.
Будем искать решение задачи (2.4.14) в виде
v(x,i) = w(z), где z = ^(t - to)2 - (® - то)2-
36
fl/-) точка (x,t) лежит на характеристиках I, или /2, то
“о Tz>оя5z
z3
z
Ut-----z Z Z, Z~
Додста^я в (2.4.14), получаем
ш"(г) + - cw(z) = О, w(o) = l. (2.4.15)
рамена ( = з/^г. »(<) = wW в С2'415) дает
»"(e)+^+y(e) = o, y(o) = i. fc.b.'iG)
Р ледовате л ьно,
y«) = W, где J0(C):=l_(|)2 +А_(|/_ А_(|)» + ...
функция Бесселя нулевого порядка [1, с. 636]. Таким образом, функ-
ия Римана имеет вид
v(x,t) = J0(^/-c((t - t0)2 - (х - io)2))-
2.5. ЗАДАЧА КОШИ ДЛЯ ВОЛНОВОГО УРАВНЕНИЯ.
МЕТОД СПУСКА АДАМАРА
смотрим задачу Коши для волнового уравнения
utt = Ди, х е R3, t > 0,
(2.5.1)
(2.5.2)
Десь х =
ИЗВестна* функция, Дц :=”Л &и
— г/
Н||=о — ф(т), Ut|(=0 = V’(l).
( 1,:с2>гз) е R , t > о независимые переменные, u(x, t) -
_ _ оператор Лапласа. Обозначим
-{(x.t): xeR3, t>0}
37
Определение 2.5.1. Функция u(x,t) называется решением зд
(2.5.1)-(2.5.2), если u(x,t) G C2(D) и u(x,t) удовлетворяет (2 J
(2.5.2).
Зафиксируем точку = (г?, х%, х®) € R3 и обозначим
г = ||х - х®|| =
, Е(^-4)2
V=i
- расстояние между точками х и хд. Пусть К(х°,а) := {х 6 R3;
а} - шар с центром в точке хо радиуса а, и дК(х°,а) := {х g
г = а} - сфера (граница шара К(х°, а)).
Лемма 2.5.1. Д(-) = 0^<' .0.. с|> - Ч й.
„ dr Xk — xi д . 1. Xi -
Доказательство. Так как —— = ---------, то ——(-) --------
\ dxk г дхк г
И д2 .1. 1 3(х* — х£)2 .1.
тгИ-) = —ч + ~g , следовательно, Д(-) = 0.
дхк г ГЛ г5 г
Вывод формулы Кирхгофа. Зафиксируем точку (х®, to) Е R4 и
смотрим в четырехмерном пространстве переменных (x,t) ков
с вершиной в точке (х°, to), образующими to — t = г и основ!
K(x°,to) С R3, лежащим в гиперплоскости t = 0 (рис. 2.5:1). I
to — t 1 '
жим v(x, t) =--------1. Из леммы 2.5.1 следует, что функция v
Г я
удовлетворяет уравнению (2.5.1): Vtt = Кроме того, v(x,t)
внутри конуса G и v(x,t) = 0 на боковой поверхности G. На
Рис. 2.5.1
38
j > 0. Обозначим fi(z°, <5) = {(x, t) e R4 : r < 6, t >
ЗафиксиРУеМосью x = xo радиуса <5 и рассмотрим Gj = G\Q(x°,<5)
I) - цилиндр c ° /см рИС 2.5.1). Граница Ss области G3 состоит
конус без пил SifU Si s и s3 S, где Si,s = {(х,t) : t = 0, i < г <
»T₽eX ХТХнования конуса, S2,s = {(х, t) : t0 - t = г, 0 < t < t0 - <5}
“боковая поверхность конуса, S3,s = {(x,t) : г = <5, 0 < t < t0 - <5} -
поверхность цилиндра.
’‘’предположим, что решение u(x,t) задачи (2.5.1)-(2.5.2) существу-
, £« Ф.РИУЛУ Ооуроу™^^ f
О = / (u(vu — Дп) — v(“tt — Ди)) <^хЛ
V-CCU. xiOMAXu
[ = !g. i^UVt -vu,)" £ ^(uVxt" 7iJ) dxdt
= Js ((uvt — vui) cos(n, t) - zZ (“Их» — l'“iJ cos(n, x*)) ds
A- , <ГС
= fSi ((v, cos(n, t) - £ vIt cos(n, x*))u
3
—v(ut cos(n, t) — £ uXk cos(n, x*))) ds,
*=1
ib n - внешняя нормаль к поверхности, a cos(n,t) и cos(n, ц) -
эсинусы углов между внешней нормалью и соответствующими осями
эординат. Выберем направление I так, чтобы cos(Z, t) = cos(n, t) и
>s(/,xt) = — cos(n,xj). Тогда
, r . av du. ,
4(“az~“ai)ds = 0-
1. Ha 5i x • n = / ~ — j j . Л
,d' dl-------Ht' dS ~ dx' ~ 0» слеД°вательно,
(u« -voi)ds = /Sii(-*£+^)dx
РИЛ-.О =4,1(v2 + (7-W))dx.
(2.5.3)
u u имеем
Ь/’“ ’§>*- /вд,., (^ + ф - 1)«M) *. (2.5.4)
//
39
Ct
2. На S2,s направление I совпадает с образующей конуса ]
dv '
которой v = 0. Следовательно, и — = 0 на S2 s. Это дает I
W ГЗ. Ha S3ts: r = 6. Тогда
(f du j ,to— t v г du .
ls„ vaids = k dt L™ d,ds-
»u OUA л 1
Так как |—| < С, to tUVJCLC j
n 1
Inn / h— ds = 0.
i->0 Js„ dl I
Далее, на S3,4: I = -n = f := (x, - x?,x2 - x%,x3 - x§). fl
dv _ dv |
~di~dr и
отсюда
/s ulfids = /s u7Tds = /s (~^-2~)\ds
1 Л-4 г
= _Й2Л ^to~t'>dt/aK(x<‘,ii^'t'>ds= -4,r/o (to-t)u(x°,i)d«
где <^o Z c|-i£fvVu ~ Ч1Г S1' 1
JS = " ~ ‘) dt /gKW)(“(^> 0 - “(Л 0) d^-
В силу непрерывности^ Q существует!
такое, что если г <д, то |u(i, t) - u(zD,tj| < е при всех t € Я
Следовательно,
|Л| - i /о °<(<о ~dt Lk^ а ds l^to - vdt-
'Z'O
Таким образом, I
fesL (uS_ ds=-4,r № ~ <)u(i°'f) dt-
йз(2.5.3)-(2-5-6)“здУеТ’1ЧТ°
£(<0 - № ^dt = ia + (г " 1)V,(l)) dx- {2-5-7)
дифференцируя (2.5.7) дважды по t0, получаем
(А(*°. «о) =^( Гdr !„к^ (“+(? - ь
ИЛИ 13, 1 г V>(x)
u(x°,to) - 47r5Zo^x(xo,to) tQ ds+4irJaK^M to ds' ^2'5'8^
формула (2.5.8) называется формулой Кирхгофа. Таким образом, мы
юказали, что если решение задачи (2.5.1)-(2.5.2) существует, то оно
1ается формулой Кирхгофа. Отсюда сразу следует единственность ре-
цения задачи (2.5.1)-(2.5.2).
^/Теорема 2.5.1. Пусть <р(х) 6 C3(R3), ф(х) 6 C2(R3). Тогда ре-
иение задачи (2.5.1)-(2.5.2) существует, единственно и дается фор-
мулой (2.5.8).
Доказательство. Достаточно доказать, что функция u(x,t), опре-
;еляемая формулой (2.5.8), удовлетворяет (2.5.1)-(2.5.2). Для это
>ассмотрим функцию ~ '
I; С ) = JL / Л£) ds ,2 5 9)
во ЛТ) достаточно гладкая функция. Сделаем в (2.5.9) замену пере-
енных xk =xi + t^k, k = U (т.е. & = Если x пробегает
^тел^нУ’<0) ’ ТО С пробегает сферу дК(0,1), и dsx = t20ds(. Сле-
з(2.5 ю) с W(x0’<o) = £/a*(o.i)/(l° + i^ds- (2-51°)
^(г), и ДУеТ’ ЧТ° ^нкция w(a;0» ^о) имеет ту же гладкость, что
w(to=o = 0.
'РУя (2.5.10), вычисляем
(У _ 1 г
3t0 4тг /дк(о,1) + ^оС) ds -j-
(2.5.11)
41
В частности, (2.5.12) дает
gU = /(x°). I
Отметим, что в (2.5.12) = cos(n,&). Поэтому, применяя Л
д I
лу Остроградского и учитывая (2.5.10) и соотношение —— = <
% “
получаем
dw = w(x°,to) to Г (Г —(?L(r0,f
dt0 to + 4тг Мод) д^дх+t°^^
C^OV- to +4%4(0,1)(^ат1( + °е))^
у_с v« w(x°,t0) 1 г ,а2/ „ ,
^сг-77- + ^/к(1»,1о)(Е^(х))^.
Таким образом,
-V = w(x°,t0) 1
СлсЦиуиЛЛ to +47Tt04’
где
(2
«- u«< i,3w)d” <2 й'*»'"'
Используя (2.5.14), вычисляем I
d2w _ 1 dw w 1 1 д
dt% t0 dt0 to 4?rt§ + 4-zrto dt0
_ 1J4 w 1 1 £Q_J_________
t0 t0 + 4Trt0 to Airto + 4nto dt0 4-rrt0 Oto
следовательно,
d2w _ 1 r . A d2 f
dt2o ~ 4%t0 (*?, dx2(l)) ds'
В частности, это дает (как и при доказательстве (2.5.11))
d2w _
dtT[[»=° - °-
(2
(2
42
— °LLfx)ds.
47rto JdK(xo,to) дх^ '
_ з d2w
(2.5.17)
(2.5.18)
„V. (2 5.10) дважды по x°k, получаем
1ифференпиРУя 1
d2v> _ r ^(x° + to?) ds = -1-1
^62 - 4?r »4 dx% 47rtn Jd
(месте c (2.5.15) это дает
d2w
dtg
функция w удовлетворяет волновому уравнению. Кроме того, ис-
ользуя (2.5.17), вычисляем
д2 .dw _ А д2 dw
эё0(дГ0'~
dw ,
е. функция — также удовлетворяет волновому уравнению. X/
Формула (2.5.8) имеет вид
и(х°, t0) = Wi(x°, t0) + w2(x°, t0),
Oto
не функции W] и w2 имеют вид (2.5.9) с. f — ip и f — ip соответ-
гвенно. Поэтому из (2.5.11), (2.5.13) и (2.5.16)-(2.5.18) следует, что
ункция и, определяемая формулой (2.5.8), удовлетворяет (2.5.1) и
'.5.2). Теорема 2.5.1 доказана. □
Исследуем теперь устойчивость решения задачи (2.5.1)-(2.5.2).
пределение 2.5.2. Решение задачи (2.5.1)-(2.5.2) называется
йчивылс, если для любых е > 0 и Т > 0 существует 6 = 6(е,Т)
х°е, что если |у,(х)-^(х)| < 6> Wx)_^)| < S, < <5
R3’ Т0 = при всех х^з, 0< t < Т.
1 решение задачи Коши при начальных данных ф и ip.
Релеления РрШение задачи (2.5.1)-(2.5.2) устойчиво в смысле
• • В самом деле, запишем (2.5.8) в виде
“(A t0) = -Lfli г t
! ^tp 4тг /ак(од) <P(I° + *о£) ds) + ip(x° + t0£) ds
= tp , Tv? ’ r
4tt -wc(o,i)
4^L(o,i)*’(l° + to?)<is
ГА 43 1
+ 4тг Ьк(о,1) £ дхк + to^ '^kds' t
Для е > О, Т > О выберем 6 = е/(4Т +1) и предположим, что IJ
£(т)| < 5, |^(х) - 4>(х)\ < 8, < <5 при всех х 1
Тогда дХк дХк
|u(x°, t0) - й(х°, to)| < /ак(од) |^(х° + to^) - ^(х° + toOl J
V +4^ /вк(0,1) l*>^° + to& ~ ^х° + to^lds
+ 4^ /att(o.i)1 “ 0^®° + ds ~ 6(iT + 1) =|
Таким образом, задача Коши (2.5.1)—(2.5.2) корректно поставлен^
Метод спуска Адамара. 1. Рассмотрим задачу Коши для
мерного волнового уравнения:
utt = Ди, х 6 В.2, t > 0, (2.5.1
^=0 = ^(1), Ui|e=0 = ^(x). (2.55
Здесь х = (xj, х3) € В2, t > 0 - независимые переменные, u(xi,»
д2и д2и
- неизвестная функция, Au := g-j 4- yj _ двумерный оператор Л
ласа. Решение задачи (2.5.19)-(2.5.20) может быть получено метода
аналогичным вышеизложенному. Однако удобнее получить репЛ
задачи (2.5.19)—(2.5.20) непосредственно из формулы Кирхгофа (2»
(т.е. считать задачу (2.5.19)-(2.5.20) частным случаем задачи (2.51
(2.5.2)).
Для этого рассмотрим задачу (2.5.1)-(2.5.2) и предположим,
функции и V не зависят от х3. Покажем, что в этом случае^)®
ция u(x°, t0), определяемая формулой (2.5.8), не будет зависеть от*
т.е. будет являться решением задачи (2.5.19)-(2.5.20). Обозначим!
y(xi-x?)2 + (xi-:i:?)2, а = <r(x?,x§,t0) = {(xi,x2) : р < to}
центром в точке (х?,!^) радиуса t0. Так как дК^х0,^ = S+(jS'il
S* - полусферы х3 = xg ± \/t)j - р* (рис. 2.5.2), то из (2.5.9) полу»
^°>^ = r/s+Zr2ds+r-/s f-?ds
47Г JS+ tQ 4тг Js- tQ
44
= 3- / /(xb z2,4 + yi§T7) J^x,
47Tto Ja cos(n, гз)
dx]dx2
cos(n, 3?з)
ds
п
/72 2
Пусть f не зависит от 13. Так как cos(n, Z3) = —-—— (рис. 2.5.3),
ТО t0
f-^^dxidx2.
w(x°.tn) = — [ —
1 J 27Г Л(х?^,«о)
Рис. 2.5.3
от ®з- тТ B Частности' вытекает, что функция w(x°, t0) не зависит
V и 0 не (^°6Разом, если в задаче (2.5.1)-(2.5.2) начальные данные
формула К ави< ят от *3, то и решение не будет зависить от Хз, и
Рхгофа (2.5.8) принимает вид
и(г1. ®°, tn) = -1A f *2) j_
2?г at0 А(»?,х",1о) dlldl2
45
Формула (2.5.21) называется формулой Пуассона. Она дает pJ
задачи (2.5.19)-(2.5.20).
2. Сделаем еще один шаг в методе спуска Адамара. Предпо™
что в (2.5.21) функции <р и не зависят от х2. Тогда
С С
где h2 = t% - (ti - г?)2 (рис. 2.5.4). Так как
fi°+h dx2
^~h 7^7
®2-Х2 |й+Л
= arcsm —— Ц_Л =
J=lG-ty^dxi-
Рис. 2.5.4
Аналогично преобразуется и второй интеграл в (2.5.21). Таким»
разом, (2.5.21) принимает вид
«(®1. fo) = | Мх° + to) + ~ to)) + | ^(з>1) dx\. (2.$
Формула (2.5.22) - это формула Даламбера (см. п. 2.1), дающая P0t
ние задачи Коши для уравнения колебания струны:
I ' ание 2.5.1. Уравнения колебаний для п = 3, п = 2 и
3аМ^зЫваются уравнениями сферических, цилиндрических и плос-
” = 1 лн соответственно. Формулы (2.5.8), (2.5.21) и (2.5.22) дают воз-
l“Ct В°сть выяснить физическую картину распространения волн. От-
М°*Н, «то при п = 3 решение u(i°,t0) зависит от начальных дан-
только на границе основания характеристического конуса (т.е. на
"лере дК(х°, to) )• ПРИ п = 2 и n = 1 решение и(х°, t0) зависит от на-
чальных данных на всем основании характеристического конуса (т.е.
в круге t°'> ”ЛИ На отрезке 1®1 ~ *0- х° + *»]) Другими слова-
ми при П = 3 начальное возмущение, локализованное в пространст-
ве, вызывает в каждой точке х° действие, локализованное во време-
ни (принцип Гюйгенса), т.е. волна имеет как передний, так и задний
фронт. Для п = 2 влияние локализованного начального возмущения
не локализовано во времени (принцип Гюйгенса не имеет место), т.е.
волна имеет передний фронт, но не имеет заднего фронта - точка х°
выйдя из положения равновесия, будет колебаться бесконечно долго.
Отметим, что задачу для п = 2 можно рассматривать как простран-
ственную задачу, когда начальные возмущения заданы в бесконечном
цилиндре и не зависят от третьей координаты.