УДК 517.956.32
ЗАДАЧИ КОШИ И ГУРСА ДЛЯ УРАВНЕНИЯ 3-ГО ПОРЯДКА
В.В. КарачиК
Рассматриваются задачи Коши и Гурса для гиперболического уравнения 3-го порядка. Доказана теорема существования функции Римана и на основе этого построены решения задач Коши и Гурса.
Ключевые слова: задача Коши; задача Гурса; гиперболическое уравнение 3-го порядка; функция Римана.
Введение
Рассмотрим уравнение третьего порядка следующего вида
L(D)u(x) ° Dpu + Y ai(x) — + c(x)u = f(x), xe G, (1)
i=i dxi
где x = (x1, x2, x3), Dx3 = Э3 / dx1dx2dx3, f e C(G), а необходимая гладкость коэффициентов и область G будут конкретизированы ниже. Очевидно, что характеристиками уравнения (1) являются плоскости Xj = const. Введем следующие обозначения. Пусть s(x) = 0 - некоторая поверхность в М3 класса C4 и x0 e М3. Область, ограниченную поверхностью s(x) = 0 и плоскостями xt = x°, i = 1,2,3, обозначим через G, а пересечение G с этими плоскостями через St. Пусть также Tj = s n Si. Вектор, получающийся из вектора x отбрасыванием i -ой компоненты, обозначим x(i). Исследуем для уравнения (1) следующие задачи.
Задача Коши. Найти функцию u(x), такую, что ue C2(G) и D^ue C(G), удовлетворяющую уравнению (1) и следующим условиям
Э ku
ЭГ
k
= jk(5), 5es, k = 0,2, (2)
где 1 e М3, || 11| = 1, l(x)e C3(s).
Рассмотрим область V = {хе М3 : х° < х1 < х1,1 = 1,з| . Пусть 5Дхк) = V = хк} , где
з
к = 0,1. В соответствии с предыдущими обозначениями = (х°), о = ^Б1 (х1).
i=1
Задача Гурса. Найти функцию u(x), такую, что ue Cl(D) и D£.2) ue C(D), D^ue C(D), удовлетворяющую уравнению (1) и следующим условиям
u( x)IS (/) =j(j = 1,2,3. (3) Функцию u(x), если она существует, будем называть регулярным решением сформулированных задач. Уравнения, содержащие операторы вида , и их итерации принято называть в
литературе уравнениями Манжерона. Если at = 0, i = 1,3, c(x) = const, то решение задачи Гурса получено, как пример, в [1]. В данной постановке задача Гурса рассматривалась также в [2, 3]. В [2] исследовалось существование решения, а в [3] строилась функция Римана. Путем применения методики исследования задач Коши и Гурса для уравнений 2-го порядка гиперболического типа, изложенной в [4, 5], получим условия разрешимости сформулированных задач.
1. Функция Римана
Введем оператор L*(D), сопряженный с L( D):
1 Карачик Валерий Валентинович - доктор физико-математических наук, профессор, кафедра математического и функционального анализа, Южно-Уральский государственный университет. E-mail: karachik@susu.ru
s
3 Ъ
L(D)v(x) ° БР\(х) + Xa,(x)v(x))-c(x)v(x). (4)
Очевидно, что оператор L (D) определен на функциях v(x), имеющих следующую гладкость: vе C^D), D^ vе C(D) и Dpve C(D). Обозначим через xA и x^ первую и вторую компоненты
вектора x(i).
Определение. Назовем функцией Римана уравнения (1) функцию, являющуюся регулярным решением следующей задачи
L*( D)v(x) = 0, xe D,
(5)
1 = w.( i , = 13
lxi=x
где функция w (x^) удовлетворяет условиям
v( x)|x.=x = W(^nl i =1,3,
ч(^ + а(%=хч(^ = х«е 51(^
11 ° } (6)
ч(^^==15=k,Л х° < х < х■
Обозначим функцию Римана уравнения (1) через Я ■ Очевидно, что Я = Я( х, х1).
Теорема 1. Пусть выполнены условия с, а, да1 /дх1 е С(V). Тогда функция Римана уравнения (1) существует и единственна.
Доказательство. Пусть решение задачи (3)-(4) - функция Я(х, х1) существует. Тогда, применяя к уравнению £ (В) Я(х, х1) = 0 оператор I ■ IX , где хе V , получим
3 х
R( x, x1) - £ (R( Aj^=x, x1) - R( ^=^, x1)) - R( x, x1) +
i=1
£¡^)(a,(X^X^,x1)- a,(X)R(Xx=x,,x1))dX) -£ c(X)R(X,x1)dX = 0.
+Х | ^ (а- (Х)Я(Х х 1, х1) - а. (£)Я(£Х =„, х1) IIX,, - Г с(£)Я(£, х1) IX = 0.
Здесь и в дальнейшем под выражением вида а(Х)Я(ХХ = х, У) будем понимать следующее выражение а(х1,Х2,Х3)Я(х1,Х3;Л,У2,Уз) ■ Учитывая условия (6), налагаемые на граничные функции задачи (5), найдем
Я(х,х1) = ]Г ч(хд) - 2 + ]ГI"11) (а,(Х)х=(х(1))-а,(Х)Я(Х|£.=х1,х1)) Х - с(Х)Я(Х,х1) IX.
1=1 1=1 ^ х
Нетрудно видеть, что
ч(х(1)) + ^ ^(Х)1Х=х1 ч(х(1)) Х =}, хе , 1 =1,3
Поэтому
R(x, x1) = 1 - £ ¡^ a, (X) R(XX=x,, x1) dX,- ¡x c(X) R(X, x1) dX, xe D. (7)
i=1 x(,) X
Обозначая R(x, x1) через v(x), запишем уравнение (7) в виде
)|X
v(x) = 1 - £¡^) a,(X) v(X)|X,=x, dX(,) - ¡^ c(X) v(X) dX. (8)
1=1 ^)
Для решения уравнения (8) воспользуемся методом последовательных приближений. Нетрудно убедиться, что Vе С(Р) и у(х) можно записать в виде
кх) = X(V,(х)- vп_l(х))+ v°(х), (9)
п=1
где v°( х) = 1, а при п > 1
уп(х) = 1 - £¡^ а(ХК_1(Х)|£.ах{1) — ¡х с(Х)уп_1(Х) ах,
1=1 х1)
если ряд (9) равномерно сходится в области V . Аналогичный результат получим после применения метода нормированных функций [6] к задаче (5)-(6). Оценим члены ряда (9). Пусть
| с|,| а{ |< М, в области V и К = тах,{х1 — х0,1}. Для сокращения дальнейших записей введем
следующие обозначения | х — у ^ = х] — у] + хк — ук,| х — у |= х1 — у + х] — у] + хк — ук , где 1, 7, к
различные числа, принимающие значения 1, 2, 3. Так как при хе V
(х1 — х)(^ — х]) < К/21 х1 — х^, (х1 — х,)(х* — х])(хк — хк) < К2 /31 х1 — х|, (10)
то для хе V справедлива оценка
| х) — Уо(х) |< М((х1 — Л1)(х2 — х2) + (х2 — х2)(х1 — х,) +
+(х — х1)(х3 — х3) + (х — х1)(х2 — х2)(х3 — х3)) < 4/3МК2 | х1 — х|.
Предположим, что для хе V и некоторого п > 1 верно неравенство
| п.1(х) — ^(х)|< К/3(4КМ)п | X — х|2п—1,!. (11)
Подставляя в (11) вместо п значение п + 1 , после несложных преобразований находим | х) — ^(х) |< К/3(4КМ)пМ(31 х1 — х|2п+1,! +1 х1 — х|2п+2,!) <
< К/3(4КМ)пМ | х1 — х|2п+1,! (3+1 х1 — х| /(2п + 2)). Так как 3 /(2п + 2) < 1, при п > 1, то
| Кп+1 (х) — ^(х) |< К/3(4КМ)п+1 | х1 — х|2п+1,!, т.е. равенство (11) справедливо для любого значения п > 1 и хе V . Применяя оценку (11) к формуле (9), убеждаемся, что функция х) существует, непрерывна в V и удовлетворяет неравенству
| Их)|< 1 + 2 / 3>/МК3 (V4МК | хд — х|).
Докажем единственность решения уравнения (8). Пусть имеется два решения. Разность этих решений обозначим и-(х). Очевидно, что и-(х) должна удовлетворять следующему уравнению
к х)=— ¡х с(Х) КХ) ах—£ ¡^ а (X) п(Х)х=х. ах{1).
х 1=1 х(1)
Учитывая сделанные выше обозначения, найдем
( 1 3 x Л
| мI jx |w(X)|dX+Xf W£%=x Xi) . ^ 1 i=1 x(i) J
Из полученного неравенства при W = sup | w(x) | будем иметь | w(x) |< 4/3MWK2 | X - x|.
xeD
Используя найденное неравенство, из оценки для w(x), аналогично (11), можно получить
| w(x)|< W/3(4KM)nK| X1 - X|2n-1J, xe D. Очевидно, что если n ® ¥ , то правая часть найденного неравенства стремится к 0, а так как левая его часть не зависит от n, то w(x) = 0 для xe D . Что и требовалось установить.
Итак, если функция R( x, x1) существует, то ее можно найти из интегрального уравнения (7). Если теперь мы покажем, что решение уравнения (8) имеет следующую гладкость: ve CX(D), D2 ve C(D), Dpve C(D), и удовлетворяет задаче (5)-(6), то функция Римана уравнения (1) существует и единственна.
Сначала докажем, что функция у(х) обладает в V непрерывными производными первого порядка. Для этого исследуем дифференцируемость ряда (9), например по х1. Очевидно, что функции Уп (х) дифференцируемы в V и справедливо равенство
Эх1
гх2
( ул+1(х) - уп(х)) = -¡х2 аз(х)( Уп(х) - Уп-1(х)),Х2=Х Х -
'\хз=хз
гх2 ,*3
х2 Г 4 Г 4
1x2 ) х3
-|хз3 в2( х) ( Уп (х) - Уп-1 ( х))\хз=хз ¿Хз - {х; | с( х) + —^ х) ( Уп (х) - Уп-1 ( Х))\Х2=Х2 ¿ХзХ -
_Э_ Э х1
хз =Хз
| 21 з ^(х)—(Уп(х)- Уп-1 (х))|х2=Х2 ХХ Эх1 1хз=хз
(12)
хз =Хз
Пусть М1 = 8ир(М, Эа1 / Э х), тогда, используя (10), нетрудно убедиться, что
хе'в
и значит
V" (У1 (х) - Уо(х))
Э х1
Э
<
М1 (2(х2 - х2)(хз - хз)+ \х - х \а)),
—(у(х) - Уо(х)) Э 1
< М1(К+1) \ х - х\а).
(1з)
Э
Теперь оценим-(уп+1(х) - уп(х)). Воспользовавшись (10), найдем
Э 1 п 1 п
(уп+1(х) - Уп(х))
Э 1
< 2/зК(К+1)М1(4КМ1)п \ х - х\2п,! +
х2 схз
+М111
•'ЛЬ
х2 •>хз
Э ( )
^(Уп - Уп-1 ) Э 1
(14)
Отсюда нетрудно получить цепочку неравенств
< 2/зК(К + 1)М1 (4КМ1)п \ х -х\2п,! +
^ (уп+1(х) - Уп (х))
Э х1
+2/зК(К+1)М2(4КМ1)п-1^(\ х1 -х\(з) +хз -Хз)2"2,! ¿ХзХ
ётзёт2ё%зё%2 << 4/зК(К + 1)М1(4КМ1)п \ х1 -х\2п,! +
х2 Г 4 г 4 г 4
+
+М2Гх2 Гхз Г4 Г4
1 ¡х2 Jхз •'Х Х
Э , Л
— (уп-1 - Уп-2 )
Э х1
+м2 Гх2 Гхз Гх2 Г
1 Л х-, ^ Х -Х
Э , Л
— (уп-1 - Уп-2 ) Э 1
с^зе^е^е^.
х2 хз Х2 Чз
Если данный процесс повторить п - 2 раз и использовать неравенство (1з), то получим
^ (уп+1(х) - Уп(х))
Э х1
< 2/зпК(К + 1)М1(4КМ1)п \ ^ - х\2п,! +(К + 3)M1n+1 \ х1 - х^
Следовательно, если ряд (9) продифференцировать по 1 , то полученный таким образом ряд
будет равномерно сходиться в области V . Значит, по известной теореме сумма ряда (9) - дифференцируемая по х в области V функция. Аналогичное утверждение можно сделать и относительно переменных х2 и хз . Итак, уе С1 (V). Очевидно, что более этого, уе
Исследуем возможность применения операторов Е^2 и Е^р к функции у(х) в области V .
Легко видеть, что правая часть уравнения (8) допускает применение к ней операторов Е-2-1 при
i = 1,3, а также оператора Dp, и полученные в результате этого функции будут непрерывны в D . Поэтому DP v, DP ve C(D). Далее нетрудно подсчитать, что
а значит
г* f ^
DX2)v+aiv=j; z(aJv)xj -cv
' \ j=i j
D2) v + av 1 = 0, xe S,.
x(i) i |x,=X ' i
dX, (15)
Поскольку кроме этого ^ _ 1 для 7 = 1,3 (I _ 7) и Vе С(V), то функции Щ(х(1)) _ v(х) _ 1 удовлетворяют условиям (6). Для окончательного доказательства теоремы подействуем на равенство (15) оператором Бх . Получим £ (П) х) _ 0 . Итак, функция х), находимая из (8), удовлетворяет задаче (5)-(6). Теорема доказана.
Замечание 1. Если с, а{ е С1(Р), Э2а; /Эх2 е С(Т>), тогда функция Римана дважды дифференцируема и Яхх (х, х1) е С(V).
Исследуем теперь свойства функции и(х) _ Я(х0, х).
Теорема 2. Пусть коэффициенты уравнения (1) обладают гладкостью, требуемой в теореме 1. Тогда и(х) имеет следующую гладкость:
ие С\Р), П^ие С(Р), Прие С(V), а также удовлетворяет однородному уравнению (1) и условиям (6)
Пх2) Щ(х(1)) + а(х)\х._х0 Щ(х(1)) _ 0, е Б;(x0),
11 0 1 (16)
Щ(х(1))\х _х0 _1,5 _ k,7х < хз < 4,
\хэ _ х0
где Щ(х(1)) _ и(х)\х_х0.
Доказательство. Положим х_ х° и х1 _ х в уравнении (7). Его решение представляется рядом (9), который перепишем в виде
и(х) _ X ип(х,х0), хе V, (17)
п_0
где и0(х, х0) _ 1, а при п > 1
ип (х, х0) с(Х)ип_1(х,Х) ах-X || а(Х)ип_1( х,Х)х_х0 Х (18)
Пусть ип (х, х0) обладает гладкостью, требуемой в теореме 2 от функции и(х) . Тогда из формулы (18) сразу следует, что ип+1( х, х0) также обладает этой же гладкостью. Это видно из того, что порядок производных справа от ип (х, х0) не может быть выше порядка производных слева от ип+х(х,х0). Оценим члены ряда (17). Из (11) следует, что
(36Мк3 )п
\ип (х, х0)\<--,
п 9(2п -1)!
если \ х1 - х0 \< к, \ а1 \,\ с \< М . Аналогичные оценки имеют место и для Пхип, Пр^^п, Прип. Значит, сумма ряда (17) обладает гладкостью, требуемой от функции и(х). Положим в формуле
(18) х1 _ х° . Тогда, обозначая ип (х) _ 0 _ Щ, (х(Л, х0) будем иметь
\х;_ х1
Щ(х^,х0) аДХ)йп-1(х(1),Х)\Х_х0 Х (19)
Положим х, = х,, для , = 1. Тогда шп
(Хл, х0 )| =0 = 0 для п > 1. Значит, (х^Д = о = 1. Подей-
х=х, <х,=х,
ствуем на равенство (19) оператором Ох при , = 1. Если п > 1, будем иметь
или при n > 2
(xw, x0)=-Q a(X) (, X)x=xo dxk - a(X)Dxj (), X)ix=xo X), k x'j=Xj 0 ' '
D
j wn (x0.),x0)=-j|) a(X)Dxj wn-i( xw,X)|X=x° dXw.
j0\_ fXo-)
Ц)
Действуя на полученное равенство оператором D , найдем
Dg W (X(„ x°)=-{| a(X)DX2.; wn-i( X(„X)|ü=x0 X
„°ч_ fx(i)
Если n = 1, то будем иметь D(2) W (x( ~, x0) = -аДx) = ° и значит
(2°)
% W)'" ' —Ix=x°
Dxl wi(x(i), x0)=-a(x)|x.=x° w°(x( ), x0).
x( i) iv (i)
Отсюда, используя (19) и (2°), нетрудно получить
Dg Wn(x(i),x0) = -a(x)|x=xo Wn-i(x(i),x0). (21)
Переходя к пределу при nи учитывая, что lim wn(x( x°) = wt(x,)), убеждаемся в справедливости (i6).
Выпишем рекуррентное дифференциальное соотношение для функций Un (x, x0). Подсчет показывает, что если i, j, k - не равные между собой числа, принимающие значения i, 2, 3
D^Un (x, x°) =
= -| [^4 (X)'ix/=x? dX + ¡^k (X)xk=x0 dX + c(X d^) ]nF)Un-i(x,X)x=Xl + (22)
+ щ(X)Dx3)Un-i (x,X)x=x, dX(i) 1 - Г° c(X)Dx3)Un-i (x, X) dX - c(x)Un-i (x, x°)
•*x(i) i * x
Очевидно, что Un (x, x) = 0, если n > i. Положим n = i в (22). Тогда, так как
U°( x, x0 ) = i, полу-
чим
Dx3)Un+i( x, x0) = -£ Щ (x) DxUn (x, x0) - c( x)Un (x, x0)
i=i
при n = 0 . Далее из (2i) следует, что при X = xj
DP]Un(x,X) = -щ(x)Un-i(x,X) .
(23)
(24)
Предположим, что формула (23) верна при п = т — 1. Докажем ее справедливость и при п = т. Для этого, подставим значения из (23) и (24), взятые при п = т, в формулу (22), в которой п = т +1. Получим
43)ит+1(х, х0) = ]Г а, (х) (Х)х.=х0 X + ^Ч (Х)&=х0 X + |хх01) С(Х)' ^^т—1(х,Х)х=х1. +
\ 3 ^
+ fx°l)а,(X) Iа/(x)DxUm_i(x,X) + c(x)Um-i(x,X)
Jx(i) j
.j=i
dXrz
(i)
Xi = xk
+
(25)
(а (x)c(X)DxUm-i( x,X) + c( x)c(X) Um-i (x,X)) d£.
i=i
Обозначим коэффициент при а1 (х) через ^(х). Применяя равенство (18), получим
Ь1 (х)=(^Ч (£)|х=%+С°ч %%=х° %+ут_1( +
+1Д^ (Х)ОхРт.1(х,%%=х0 ЧХ(у) + с(%)х,%) = -(х,х°). Аналогично, обозначая коэффициент при с(х) через Ч(х) получим
Ч( х) = I |Т а (Ххит-1 (х,Х)|Х=х° ЧХ(1) + Г° с(Х)ит-1 (х,Х) чх=-^т (х, х°).
1=1 х(1) ь 1~х1 ■'х
Подставляя полученные значения коэффициентов Ь1(х) и Ч(х) в формулу (25) получим (23) при п = т. Значит, формула (24) верна для любого п е N . Для окончательного доказательства теоремы перейдем в формуле (23) к пределу при п ® ¥ . Тогда будем иметь С(П)и(х) = 0 . Теорема доказана.
Замечание 2. Если коэффициенты имеют гладкость с, а1 е С1(П), Э2а1 / Эх2 е С(V), тогда Яхх (х0, х) е С (V).
2. Задача Коши
Исследуем задачу Коши. Рассмотрим следующие тождества для функций и, Vе С3(С)
уС(П)и + и£*( D)v=
= ( их2х3 V + а1 uv)x + ( -их3 ^ + а2uv)^ + ( и^х2 + а3 uv)^ ,
=(
+ a1uv) + ^uv) + + auv
х2 х3
^х3 ^ °2 и ^ + (■~ихг + а3 uv)^ ,
= (-их2 + а1 uv) + (и^. + а,uv) + ( и^ v + аuv]
а21
12 +(их:х2 V + а3т)х3.
Складывая их, найдем
+
3(кСф)и + и€ ( ^) = (их2х3 v- их2 ^ + и%х3 + 3a1uv) х + ( и*!*3 V - их3 ^ + и\х3 + 3a2uv) +( и^х2 V - их! ^ + и\х2 + 3a3uv)
(26)
х2 V 1 2 12 12 х3
Из обозначений, сделанных в начале статьи, следует, что область С определяется точкой х° , если поверхность о(х) = 0 фиксирована. Поэтому С = С(х°). Вернемся к предыдущему тождеству (26) и заменим в нем х на %, вместо и%) подставим Я(%, х) и проинтегрируем по С(х). Обозначим интегралы от каждой скобки через ^ , 72 и 73 соответственно. Компоненты вектора внешней нормали к поверхности о(х) = ° обозначим п%%, 1 = 1,2,3 . Вычислим 7 . Имеем
71 = х)
=I
%3 К- 1/2 ( % КХ3 + иХ3 КХ2 ) + иК%2%3 + 3а1иК]% = %Х31/2 (% %+% Ч)+и +3а1и к
Используя свойства функции Римана (6), найдем
71 = 0 х) Ги%2%3 1/2 (% + и%3 % )+ и Ч%3 + 3а1и К
п% Чо.
п% Чо-
^х) |(иК)%2%3 - 3/2 [( и Я%2 )%3 + ( иЧ )%2 ] + 3( Я%2%3 + а1 Я) и|
Здесь о( х) - часть поверхности о( х) = °, находящаяся в области С(х). Обозначим Р(х) = Ту(х) о тк(х). Интеграл по 51(х) обозначим через 11 и вычислим его. Опять используя свойство функции Римана (6), найдем
1 = 2 к (х) ][(и К)%2 - 3(и )% + [(и ^)%3 - 3( и К%3 )
2* Э1( х)
]%3
%2
Чз =
= 2% - 3(иЯХ)]Х + Я)Хз - 3(иЯХз)]X} *-
-212 (иЯ)Х Х -21з (и% Х = = 2х){[(и% - 3(иЯХ)]X + [(и% - 3(иЯХз)]X } *-
-1/2и(Р2)Я(Р2,х) - 1/2и(рз)Я(P3,х) + и(х), где п'х, - компоненты вектора внешней нормали к кривой Ь1 (х) в плоскости Х = х1. Итак, имеет следующее значение
Л = ^х) [иХ2Хз 1/2( Х % + иХз КХ2 ) + и^Хз + за1иЯ
п^ (1а-
Пе } ёэ +
- 2 Ь х){[(и ^ - з(и ЯХ2 )] п^з + [(и % - з(и %)] 42 }
+1 / 2и(Р2ШР2,х) + 1/2и(Рз)Я(рз,х) - и(х). Аналогично вычисляются интегралы У2 и Уз. Складывая значения У2 и Уз, получим
и( х) = з и (^ВД, х) + з и( Р2)Я( Р2, х) + з и (рз)Я( рз, х) -
1 з
-1 I [
2¿1 Ь (х)
(Х Я- 2и)4 + (% Я- 2иЯ& )%
((Ь +
з [ 1
+х) ^ Я 2(
Я--1 их. Я& + и^к Ях I
+ и Яе е + за, и Я
)-
ХХк-
Х1(1а-|с.х)Я(Х, х) /(£) (1Х. (27)
1С (х)
Для того, чтобы полученная формула давала решение задачи Коши, необходимо выразить значена а( х) = 0 через данные Коши. Покажем, что это возможно, когда направление
ния ихк , ихЛ
1( х) не касается поверхности а( х) = 0 ни в одной точке. В случае, когда уравнение имеет второй порядок, это известно [5].
Пусть поверхность а(х) = 0 на некоторой своей части задается уравнением х1 = а(х2, хз) и
<ге О4 . Можно считать, что (рк = ((х2, хз) . Легко показать, что на поверхности а(х) = 0 верны равенства
ди Эх1
ди _ Э(0 Эи
х2
а2 + -— = аз + -— = —11— +12 . .з
Эх2 Эх2 Эх1 Эхз Эхз Эх1 Эх2 Эхз
в которых а1 = да/Эх1. Полученная система уравнений для определения Эи/Эх^а разрешима при любой правой части, если А ° (1, п) = 0 в произвольной точке поверхности а(х) = 0 . Далее из (28) нетрудно получить систему уравнений для определения Э2 и / Эх^-Эх^. Если неизвестные записать
в виде вект°ра X = ( их1х1 , их1х2 , их1хз , их2х2 , их2хз , ихзхз )|а , то матрица этой системы
А X = Ь, (29)
будет иметь вид
Эи = Э(0 ? Эи
хз
Эи Эи +1,
■ = (\,
(28)
' а2&з аз а2 0 1 01
2а2 0 1 0 0
А = аз2 1lа2 0 11 + 12а2 2аз 1за2 0 12 0 1з 1 0
1lаз 12аз 11+а 0 12 1з
1 А2 21112 211з 12 12 2121з 12 1з /
а вектор правой части записывается в форме
Ь =
(
Эх2Эх3
Э 2(
"х^23 Эх2
Э 2(
'^а22' Эх2
Эл -Л Э1, Эи Эл -Л Э1, Эй — 1 Э/у Эи
Эх2 ,-=1 Эх2 Эха ' Эх3 ,-=1 Эх3 Эх,' а ' 2 и7=1 1 Эх, Эху \а
Iй у
После элементарных преобразований определитель матрицы А примет вид detA = А4. Сформулируем основной результат для задачи Коши.
Теорема 3. Задача Коши (1)-(2) имеет единственное решение, если с, е С1(С), Э2а,. /Эх2 е С(С), (рк е С4-к(а) и векторное поле 1(х) не касается поверхности ае С4 ни в одной точке. Это решение можно найти по формуле (27), если Эи / Эх,\а определить из системы (28),
а Э2и/ЭхДхуа из системы (29).
Доказательство. Единственность решения задачи (1)-(2) следует из представления (27). Для доказательства существования решения сделаем в уравнении (1) замену переменных
V = и - и,
(х1 -а)их1а- (х1-а) и
Х1Х1\а •
где их\а и ихх\а находятся через ( из (28) и (29). Тогда, уравнение (1) примет вид C(D)v = Г1, а
условия (2) станут однородными. В силу замечаний к теоремам 1 и 2 имеем ихх е С (С). Ясно,
что Л е С(С(хи). Если теперь показать, что функция
к х) = -С (х)Я(Х, х) /;(Х) аХ
является решением однородной задачи (1)-(2), то теорема будет доказана. Нетрудно подсчитать, что
Эv
= 15..(х)ЯХ х) Л ^Х,=х, Х - Сх)х) Л Х аХ
Эх,. и, (х) Э2v
С(х) Эх,.
Э Я.
ЭхА х Я(Хх) «%=«аХ>+и < х) ЭЯХх) =х,х+
4=хк
+15к(х)Э^, х) ^1(Х)\Хк=«"¿(к) - [
Э 2 Я С( х) Эх,.Эхк
(Х, х) 1(Х) с!Х,
Эъ\
Эх1ЭхуЭхк
= Я( х, х) /¡( х)
Л гР ( х) Э Я
1Р -
Эх
(Х, х) /¡(Х)\хк=хках +
\Ху=х,
а поэтому
3 Э2р ЭЗ Р
+— 1 (Х, х) л (Х)\Х = х аХ(1) -1 Э Я (Х, х) л (Х) аХ,
5 (х)Эх,.Эхк ' х Ь(,) •'с(х)Эх,.Эх,Эх* 1
¿(^ = х) - -^(х)х) /¡(Х)Хк=хк аХ +
й-1 5(х}
( Э2 Я
Эхк Эх,
\ к 7
\
(Х, х) + а,. (х)Я(Х, х)
л (Х)\Х=х, аХ(,) - |С (х) ¿х (D) Я(Х, х) л (Х) аХ.
В силу свойств функции Римана L(D) v = Л1. Легко видеть, что = vx\а = v\а = 0 для любых, значений ,, к = 1,2,3, не обязательно различных. Теорема доказана.
3. Задача Гурса
Для исследования существования и единственности решения задачи Гурса необходимы следующие рассуждения. Пусть и(х) и ^х) обладают гладкостью, требуемой от решения задачи Гурса, и кроме того е С(Т>). Нетрудно убедиться, что
Пх3)(^) _ vD?u + иПх\ + X(Vx¡DX¡¡)и + и^ V).
1=1
Если теперь к обеим частям полученного тождества добавить следующее равенство
3
vC(D)u + и£*(ПК_ vDX')u + иПхз) V + X(а1иу)х ,
(30)
1=1
33 П^) + X (_ кС(П) и + и£*(П) V + £ (VxJ. и + v).
1=1 1=1
справедливое в V , то будем иметь
3
'.1 а 1=1
Рассмотрим аналог тождества (30)
П2 (иу) _ и + иП(2) V + их vx + их vx .
х(1'Г 7 х(1) х(1) хк х7 х7 хк
Положим в нем V _ vx и просуммируем по 1. Будем иметь
(31)
X П2) ) _ X П2) и + 3иПхз) V + 2Х их П2) V.
¿—I х(1) хг х1 х(,) х х1 х(1)
1=1 1=1 1=1 Вычтем полученное тождество из равенства (31). Тогда, после некоторых простых преобразований, получим следующее тождество в V
3 3
Пx3)(uv) -X ^х) _ vC(D)u + и£(П^-X и(^ V + V)
1=1 1=1
(32)
Выберем вместо х) функцию Римана уравнения (1), т.е. положим х) _ Я(х, х1). Это возможно в силу теоремы 1 и включения Пx3vе Сф). Тогда ¿( П) V _ 0. Затем положим х _Х и проинтегрируем последовательно по X (1 _ 1,2,3) в пределах от х0 до х{. Получи
м
и(х)Я(х,х)- и(х0)Я(х0,х)-X(и(х)Я(^_х?,х)- и(х0)Я(х^,_х,х))
1=1
-X /х
1=1 '
3
и(х) Ях (х\х1 _х, х) - X и(х) (х\х _х0, х) + и(х0) Я^ (х^р, х)
1=1,Г /
х _х
а£_
3
_|х0я(х,х)Сх(П)и(Х)ах-XÍХрi)(и(Х)(П2я(Х,х) + а (Х)Я(Х,х)))\х ^ ах№. (33)
Предположим, что
х(1))\х._х0 (х(7))\х1 _х? (хк);
У (0) (0) _Ук (0) _У0,
(34)
где 1, 7, к — различные числа, принимающие значения 1, 2, 3. Учитывая, что функция Я(х, х1) является решением задачи (5)—(6), получим
и(х) _ XX (х(1)) Я(_х0, х) - у (х) Я(х^°х0_х, х)) + у Я(х0, х) -
1=1
(
-XX /х
1=1 1
3
у(Х)Ях(й_х,х)- X ъ(х(7))Ях(х\х,_х0,х) 7=1,7_1 | 7 X
х1 _Х
х+
ÍХрi)Ь(Х(1))(Пх2)я(Х,х) + а (Х)Я(Х,х))х_0 ах{1) + /х /(Х)Я(Х,х)ах. (35)
Итак, сформулируем основной результат для задачи Гурса.
Теорема 4. Пусть f, c, at, Эа; / Эх; е C(D), (i = 1,3), а функции j (х(;)), удовлетворяющие условиям согласования (34), имеют гладкость E^je C(D), тогда регулярное решение задачи
Гурса существует, единственно и записывается в виде (35).
Доказательство. Очевидно, что Е^ие C(D), а поэтому единственность решения задачи Гурса в требуемом классе функций следует из представления (35). Докажем существование. Нетрудно видеть, что правая часть формулы (35) обладает гладкостью, необходимой для искомого решения. Проверим, что и(х), находимая из (35), удовлетворяет уравнению (1) и краевым условиям (3). Начнем с последнего. Если воспользоваться свойствами функции Римана (5)-(6) и условиями согласования (34), то нетрудно убедиться, что интегральные члены в формуле (35) при
xt = х0 обращаются в нуль, а сумма членов, не содержащих интегралы, будет равна j (x(i)). Значит, условия (3) выполнены. Обозначим
и0(х) = JX f(X)R(X,х)dx. (36)
Покажем, что L(Е)и0(х) = f (х) . Действительно,
43Ч(х) = f (х) + f (X) R(X, х)|£.=х. dX(i) + j; f (X) Е(3) R(X, х) dX,
i=1
= ПХЖХ,х^Х + R(x,х)Х
и значит, используя свойство функции Римана (16), получим
С(П)Щ( х) = /(х)+£0 f (Х)£>х (D)R(Х, х) Х
Опять вспоминая теорему 2, убеждаемся в верности доказываемого равенства.
Обозначим и1(х) = и(х) - и0(х). Покажем, что £(П)и1(х) = 0 . Очевидно, для и1(х) верна формула
(35), записанная при /(х) = 0 . Если подставим в эту формулу выражения
Ъ(^ = и(^=х? , у (х1) = и(х) ^=х0 ,
х1=х°
где числа 1, 7, к не равны между собой и принимают значения 1, 2, 3, то в силу выполнимости условий (3) получим формулу (33) без слагаемого
{(Х)£х(°)и(Х) Х
Далее, проводя рассуждения, обратные сделанным при выводе формулы (33), найдем
3
!х0(43) (иДхШХ,х))-££ХХ2) (Ц((Х,х)) + £ц(х)(4^(Х,х) + а (Х)R(Х,х)) )Х = 0.
х 1=1 1=1 Х
Положим х = х1, х° = х, R(х,х1) = у(х) и применим к обеим частям полученного равенства оператор Бх3"1. Будем иметь
Б?^- £К) + X(и (Б* V + а(хЙ' = 0.
1=1 1=1 х1
Вспоминая тождество (32) и рассматривая его при ^х) = R(х, х1), получим vC(D)u1 = 0 . Так как Кх) ° 0, то С(Б)и1 = 0 . Теорема полностью доказана.
Замечание 3. Функция и0(х) из равенства (36) является решением неоднородной задачи Гурса с однородными краевыми условиями.
Пример. Исследуем задачу Гурса при (х) = 0 для 1 = 1,2,3 и с(х) = 1
Б? и + 1и = f( х), и( х)^ (х0) = Ъ( х^), 1 = 1,2,3.
Функция Римана этого уравнения легко находится из формулы (7). Она имеет вид
и
Ri( XxX) = Z (-l)n (* -X)n,!( -&)n,!( *3 -X)n,!.
n=0
Решение задачи Гурса (35), которое с помощью метода нормированных систем функций [6] и операторов U(D) и V(D) из [6, теорема 4.13] можно также записать как решение задачи Cp при
Ь = (1,1,1), легко преобразуется к виду
з * ( 3 ^ *
U( *) = Z (j(ХИУ>(*) ) + Уо-ЦХХ01 z (j (t(i))-У (t )) + Уо RlCt x) dt + j*0 f (t)Rl(t, *)
где функции y (*) необходимо брать из (34). Найденное решение совпадает с полученным другим путем в [1]. Метод нормированных систем функций [6], примененный к уравнению Лапласа позволяет строить специальные полиномы [7, 8] и решение задачи Дирихле [9], а примененный к линейным ОДУ - новое представление решений задачи Коши [10].
Литература
1. Карачик, В.В. Разработка теории нормированных систем функций и их применения к решению начально-краевых задач для уравнений в частных производных: дисс. ... д-ра физ.-мат. наук / В.В. Карачик. - Ташкент, 2001. - 213 с.
2. Мюнц, Г. Интегральные уравнения / Г. Мюнц. - М.: ГТТИ, 1934. - 330 с.
3. Жегалов, В.И. Трехмерный аналог задачи Гурса / В.И. Жегалов // Некласические уравнения и уравнения смешанного типа. - Новосибирск: Ин-т матем. СО АН СССР, 1990. - C. 94-98.
4. Бицадзе, А.В. Уравнения математической физики / А.В. Бицадзе. - М.: Наука, 1976. -296 с.
5. Владимиров, В.С. Уравнения математической физики / В.С. Владимиров. - М.: Наука, 1976. - 436 с.
6. Карачик, В.В. Метод нормированных систем функций / В.В. Карачик. - Челябинск: Издательский центр ЮУрГУ, 2014. - 452 с.
7. Karachik, V.V. On some special polynomials / V.V. Karachik // Proceedings of the American Mathematical Society. - 2004. - Vol. 132, no. 4. - P. 1049-1058.
8. Карачик, В.В. Полиномиальные решения дифференциальных уравнений в частных производных с постоянными коэффициентами I / В.В. Карачик // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». - 2011. - Вып. 4. - № 10(227). - С. 4-17.
9. Карачик, В.В. Полиномиальные решения задачи Дирихле для бигармонического уравнения в шаре / В.В. Карачик, Н.А. Антропова // Дифференциальные уравнения. - 2013. - Т. 49, № 2. -C. 250-254.
10. Карачик, В.В. Метод построения решений линейных обыкновенных дифференциальных уравнений с постоянными коэффициентами / В.В. Карачик // ЖВМиМФ. - 2012. - Т. 52, № 2. -С. 237-252.
Поступила в редакцию 24 декабря 2014 г.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2015, vol. 7, no. 2, pp. 31-43
CAUCHY AND GOURSAT PROBLEMS FOR DIFFERENTIAL EQUATION OF THIRD ORDER
V.V. Karachik
Cauchy and Goursat problems for the hyperbolic equation of third order are considered. The theorem of the existence of the Riemann function is proved and on the basis of the aforementioned function solutions for the Cauchy and Goursat problems are introduced.
Keywords: Cauchy and Goursat problems; hyperbolic equation of third order; Riemann function.
References
1. Karachik V.V. Razrabotka teorii normirovannyh sistem funkcij i ih primeneniya k resheniyu nachal'no-kraevyh zadach dlya uravnenij v chastnyh proizvodnyh. Diss. dokt. fiz.-mat. nauk (The development of the theory of normed systems of functions and their application to the solution of the initial-boundary value problems for equations with differential derivatives. Dr. phys. and math. sci. diss.). Tashkent, 2001. 213 p.
2. Munz G. Integral'nye uravneniya (Integral equations). Мoscow, GTTI Publ., 1934. 330 p.
3. Zhegalov V.I. Trechmerniy analog zadachi Gursa (Three-dimensional analogue of the Goursat problem). Neklasicheskie uravneniya i uravneniya smeshannogo tipa (Non-classical equations and the equations of mixed type). Novosibirsk, Institut matematiki SO AN SSSR Publ., 1990. pp. 94-98.
4. Bitsadze A.V. Uravneniya matematicheskoy fiziki (The equations of mathematical physics). Moscow, Nauka Publ., 1976. 296 p. (in Russ.).
5. Vladimirov V.S. Uravneniya matematicheskoy fiziki (The equations of mathematical physics). Moscow, Nauka Publ., 1981. 512 p. (in Russ.).
6. Karachik V.V. Metod normirovannyh sistem funkcij (Method of normed systems of functions). Chelyabinsk: Izdatel'skij centr YuUrGU, 2014. 452 p.
7. Karachik V.V. On some special polynomials. Proceedings of the American Mathematical Society. 2004. Vol. 132, no. 4. pp. 1049-1058.
8. Karachik V.V. Polinomial'nye resheniya differencialnykh uravnenij v chastnykh proizvodnykh s postoyannymi koefficientami I [Polynomial solutions to partial differential equations with constant coefficients I]. Bulletin of South Ural State University. Series of "Mathematics. Mechanics. Physics". 2011. Issue 4. no. 10(227). pp. 4-17.
9. Karachik V.V., Antropova N.A. Polynomial solutions of the Dirichlet problem for the biharmonic equation in the ball. Differential Equations. 2013. Vol. 49, no. 2. pp. 251-256.
10. Karachik V.V. Method for constructing solutions of linear ordinary differential equations with constant coefficients. Computational Mathematics and Mathematical Physics. 2012. Vol. 52, no. 2. pp.219-234.
Received 24 December 2014
1 Karachik Valeriy Valentinovich is Dr. Sc. (Physics and Mathematics), Professor, Mathematical and Functional Analysis Department, South Ural State University. E-mail: karachik@susu.ru