SSSR, 1973, vol. 210, no. 1, pp. 15-17 (in Russian). 11. Markushevich A. I. The theory of analytic functions : in 2 vol. Moscow, Nauka, 1968, vol. 2, 624 p. (in Russian).
12. Salimov R. B., Shabalin P. L. The regularizing factor method for solving a homogeneous Hilbert problem with an infinite index. Russian Math. [Izvestiya VUZ. Matematika], 2001, vol. 45, iss. 4, pp. 74-77.
УДК 517.95; 517.984
О КЛАССИЧЕСКОМ РЕШЕНИИ ОДНОЙ СМЕШАННОЙ ЗАДАЧИ ДЛЯ ВОЛНОВОГО УРАВНЕНИЯ
А. П. Хромов
Доктор физико-математических наук, профессор кафедры дифференциальных уравнений, Саратовский государственный университет им. Н. Г. Чернышевского, [email protected]
В статье методом Фурье дается классическое решение смешанной задачи для волнового уравнения с комплексным потенциалом при минимальных условиях гладкости начальных данных. Используется резольвентный подход, состоящий в привлечении в формальном решении метода Коши - Пуанкаре интегрирования резольвенты соответствующей спектральной задачи по спектральному параметру, не требующий никакой информации о собственных и присоединенных функциях и использующий лишь главную часть асимптотики собственных значений. Существенно используется прием А. Н. Крылова об ускорении сходимости рядов Фурье. Граничные условия таковы, что спектральная задача допускает кратный спектр и бесконечное множество присоединенных функций, что создает дополнительные трудности при анализе формального решения.
Ключевые слова: формальное решение, спектральная задача, резольвента, классическое решение.
ВВЕДЕНИЕ
Настоящая публикация приурочена к 150-летию со дня рождения выдающихся отечественных ученых В. А. Стеклова (1884-1926) и А. Н. Крылова (1883-1945), внесших весомый вклад в решение смешанных задач методом Фурье.
Обоснование метода Фурье в задачах математической физики традиционно опирается на доказательство равномерной сходимости ряда, представляющего формальное решение задачи, и рядов, полученных из него почленным дифференцированием нужное число раз. В. А. Стеклов, впервые давший строгое обоснование метода Фурье, придерживался этой точки зрения [1, с. 224], которая сделала метод Фурье очень популярным. Было проведено большое количество исследований и достигнуты значительные успехи.
Информация обзорного характера содержится, в частности, в книгах И. Г. Петровского [2], В. И. Смирнова [3], О. А. Ладыженской [4], В. А. Ильина [5] (см. также [6]), В. А. Чернятина [7]. Недостатком такого подхода является требование завышенной гладкости на начальные данные. Выход из этого положения намечен А. Н. Крыловым в его исследованиях по ускорению сходимости рядов Фурье и им подобных [8]. Суть его приема состояла в том, что изучение вопроса о дифференцировании ряда Фурье решается путем разбиения его на два ряда, один из которых точно суммируется (в этом случае не надо прибегать к почленному дифференцированию), а второй ряд сходится настолько быстро, что его можно почленно дифференцировать. На ряде прикладных задач им были успешно преодолены трудности, связанные с невозможностью почленного дифференцирования. Приведем его слова: «Изложенный прием усиления быстроты сходимости рядов Фурье и нахождения производных от функций, ими представляемых, может служить для доказательства или проверки того, что представляемая рядом функция действительно удовлетворяет тому дифференциальному уравнению, как решение коего она найдена, хотя бы самый ряд и нельзя было дифференцировать почленно требуемое число раз.» [8, с. 227].
В. А. Чернятин [7] приемом А. Н. Крылова с применением асимптотик для собственных значений и собственных функций успешно исследовал ряд задач методом Фурье и значительно ослабил условия гладкости, а в ряде случаев эти условия гладкости стали минимально возможными. Переход от формального решения к новому виду, вытекающему из исследований А. Н. Крылова и В. А. Чернятина, есть качественно новый шаг в методе Фурье, позволяющий с исчерпывающей полнотой исследовать краевые задачи методом Фурье и ставящий много новых вопросов в теории функций.
В настоящей статье дается дальнейшее развитие метода А. Н. Крылова и В. А. Чернятина путем привлечения метода Коши - Пуанкаре контурного интегрирования резольвенты оператора, порождаемого спектральной задачей из метода Фурье. Его преимущество сказывается особенно в случае несамосопряженной задачи из метода Фурье, в том числе и когда участвуют присоединенные функции в любом количестве. Такое положение и наблюдается в настоящей статье.
Мы будем изучать задачу
д2 д2 u(x,t) „ ч „ г
—= J2 - q(x)u(x,t), x е [0,1], t е ю), (1)
dt2 дх2
u(0, t) = uX(0, t) - uX(1, t) - au(1, t) = 0, (2)
u(x, 0) = p(x), u'(x, 0) = 0, (3)
где q(x) е C[0,1] и комплекснозначна, a — комплексное число, условие u'(x, 0) = 0 берется для простоты изложения.
Мы привлекаем еще спектральную задачу для оператора L0:
Loy = -y''(x), y(0) = y'(0) - y'(1) = 0.
Этот оператор рассматривался В. А. Ильиным (см. [5,6]) и замечателен тем, что все его собственные значения двукратны, причем для каждого собственного значения оператора Lo имеется одна собственная и одна присоединенная функции. Оператор L0 участвует в формировании эталонной смешанной задачи (см. [9]) при использовании приема А. Н. Крылова по усилению скорости сходимости рядов, подобных рядам Фурье, т.е. L0 выполняет одну из важнейших функций в методе Фурье. В результате получаем классическое решение задачи (1)-(3) при минимальных требованиях на p(x): p(x) е C2 [0,1], комплекснозначна, причем
р(0) = р'(0) - р'(1) - ар(1) = р''(0) = 0. (4)
1. ПРЕОБРАЗОВАНИЕ ФОРМАЛЬНОГО РЕШЕНИЯ
Метод Фурье связан со спектральной задачей для оператора L:
Ly = -y''(x) + q(x)y(x), y(0) = y'(0) - y'(1) - ay(1) = 0.
Теорема 1. Собственные значения An = рП (A = p2, Re ^ 0) оператора L образуют две серии {A'n}, {АП} с асимптотиками
A'n = p'n2, p'n =2nn + АП = рП2> рП = 2nn +
где n = n0,n0 + 1,..., e'n и е'П есть O(1/n). Если e'n = е'П, то они простые, если e'n = е'П, то двукратные.
Это хорошо известный факт (см., например, [10, с. 74-75]).
Замечание. В [10] приведена более точная асимптотическая формула для собственных значений, но она нам не потребуется.
Обозначим Yn = {р| |р - 2пп| = 5}, где 5 > 0 и достаточно мало, а n ^ n0 и n0 таково, что при всех n ^ n0 внутрь Yn попадают p'n и рП и других корней квадратных из собственных значений нет ни внутри, ни на границе Yn. Пусть 7n — образ Yn в А-плоскости (А = р2, Reр ^ 0). Обозначим через RA резольвенту оператора L, т.е. RA = (L-AE)-1, где E — единичный оператор и А — спектральный параметр. Формальное решение по методу Фурье возьмем в виде (см. также [11,12]):
u(x,t) = - -—т (RAp)cos рЫА - ^^ -—: /(rap)cos рЫА, (5)
I л I n>no ~
|A|=r ^ 0 7„
где r > 0 фиксировано и взято таким, что все собственные значения по модулю меньше r, имеют номера, меньше n0; на контуре |А| = r нет собственных значений. Таким образом, получаем, что в формальном решении u(x, t) не фигурируют ни собственные значения, ни собственные функции.
Проводим теперь дальнейшее преобразование ряда (5) с использованием эталонной задачи (см. п. 3).
Лемма 1. Пусть до не является собственным значением оператора L и таково, что |до | > r и до не находится ни внутри, ни на границе замкнутого контура 7 ни при каком п ^ n0, где 7 либо |А| = r, либо 7n. Тогда
I (RAy)cos ptdA = 1 --(RAg) cos ptdA, (6)
I I A - До
Y Y
где g = (L - доE
Доказательство. Так как y(x) G DL (область определения оператора L), то
g(x) = (L - AE)p(x) + (A - до)p(x).
Отсюда R\g = y(x) + (A — до)R\у и поэтому имеет место (6). □
Наше преобразование формального решения с учетом эталонной задачи дается следующей теоремой.
Теорема 2. Для формального решения u(x,t) имеет место формула
u(x, t) = ио(x, t) + ui(x, t) + u2(x, t), (7)
где
1 r 1 _ 1 r 1
uо(x,t) = -— т-(R° g)cos ptdA - V] — --(R° g) cos ptdA =
2m J A - до ^ 2ni J A - до
I л I n>U0 ~
| A|=r ^ 0 Yn
1
2n~i
( \ / + J/
\[A|=r n^n0Y„
1 f 1
1 (RA g)cos ptdA,
A - до
u1 (x, t) =--I --[RA g - R°A g] cos ptdA,
2ni J A - до
|A|=r
11
u2(x,t) = -— У] -- [R\g - RAg]cosptdA,
2ni ^ A - до J
n>n0 Yn
где R°A = (L0 - AE)-1. Считаем еще, что ранее приведенные требования на до выполняются и для оператора Lid, при этом все собственные значения оператора Lo вне |A| = r попадают внутрь контуров 7n при n ^ по.
В дальнейшем будем использовать вместо (5) именно представление (7).
2. ВСПОМОГАТЕЛЬНЫЕ УТВЕРЖДЕНИЯ
Рассмотрим уравнение
У" - q(x)y + p2y = 0. (8)
Пусть z1 (x,p), z2(x,p) — решения уравнения (8) с начальными условиями:
z1 (0,p) = 1, zi (0,p) = 0, z2 (0,p)=0, z2 (0, p) = 1.
Если q(x) = 0, то эти решения есть z°(x, p), z2(x,p). Таким образом,
о/ \ о/ \ sin px
z-г (x, p) = cos px, Zo (x, p) = -.
p
Теорема 3. Для RAf имеет место формула
raf = v1 (x p)(f, z1) + v2(x, p)(f, z2) + Mpf,
где
Mf = M (x, С, p)f (C) dC, M (x, С, p) =
Zi(x,p) Z2(x, p) zi (С,Р) z2 (С,Р)
vi(x, p) = 1 . (z2(1, p) + az2(l, p))z2(x, p), V2(x, p) = - (z1(l, p) + azi(l, p))z2(x, p)
A(p)
A(p)
A(p) = z2(0, p) - z2(1,p) - az2(1,p), (f,g) = y f (С)g(C) dC-
0
Доказательство. Общее решение неоднородного уравнения
y'' - q(x)y + p2У = f (x)
(9)
есть
у(х, р) = (ж, р) + С2^2(х, р) + Мр/(х), (10)
где с1, с2 — произвольные постоянные. Если у = , то у удовлетворяет (10) и граничным условиям
Получаем отсюда
где A(p) = «2(^2). Далее,
ui(y) = y(0) = 0, U2(y) = y'(0) - y'(1) - ay(1) = 0.
ci =0, C2 A(p) + U2(Mpf) =0,
U2 (Mpf) = -
zi(C,p) z2 (C,p) zi(1, p) + azi(1, p) z2 (1, p) + az2 (1, p)
f (C) dC-
Отсюда
«2(Мр/) = -(/, ¿1) (4(1, р) + 0^2(1, р)) + (/, ¿2) (4(1, р) + а*1 (1, р)). Таким образом, получаем:
Дл/ = [(/' ¿1) (4(1, р) + °^2(1, р)) - (/, ¿2) (¿К1, р) + р))] ¿2(х, р) + Мр/-А(р)
Теорема 4. Для Я°х/ имеет место формула
□
где
Raf = v0(x, p)(f, z0) + v0(x, p)(f, z0) + M0f,
v0(x, p) = (1, p)z2 (x, P), ^iO^ p) = Zo'(1, p)z2 (x P),
A0 (p) A0 (p)
A0(p) = 1 - cos p, Mf =
zi (x, p) (x,p)
zi (C,p) z20(C,p)
f (C) dC-
Доказательство такое же, как и доказательство теоремы 3 (здесь u0(y) = y(0), (у) = у'(0) -- У'(1)).
Теорема 5. В полосе | Im p| ^ h (h > 0 любое) имеют место асимптотические формулы:
zi(x, p) = cospx + О ( - ) , z'(x, p) = -psinpx + O(1), p1
p
sin px 1
Z2 (x, p) =--+ О —z
p p2
z2(x, p) = cos px + О
где оценки O(...) равномерны по x е [0,1].
x
1
1
x
1
p
Эти оценки содержатся в [10, с. 59].
Теорема 6. Для zj(x,p) (j = 1, 2) имеют место формулы:
x
z1 (x, p) = cos px + J K1(x,£)cos p£d£, (11)
о
x
z2 (x,p) = ^ +/ K^x.O^di, (12)
pp
о
где Kj(x,£) (j = 1, 2) непрерывно дифференцируемы по x и £, K2(x, 0) = 0.
Замечание. Формулы (11), (12) хорошо известны как формулы операторов преобразования (см. [13, с. 17, 33]).
Лемма 2. При p G Yn имеют место асимптотические формулы:
vj (x,p) = v0{j)(x1p) + O(pj-2), (13)
vj (x,p) = v0(j)(x,p) + O(pj-1), (14)
(•) dj
где j = 0,1, 2, v(j)(x, p) = -— vs(x, p) и оценки O(...) равномерны по x.
dxj
Доказательство. Имеем:
A(p) = z2(0, p) - z2(1,p) + az2(1,p) = 1 - cos p + O ^= Ao(p) + O ^^ . Отсюда при p G Yn имеем:
A(p) Ао (p)-
Поэтому
A(p)(z'2 (1,p)+ az2 (1'p))=A^ ^ (1'p)+ O ( p ) , A(p)(z1 (1,p) + az1 (1,p)) = A^)z°'(1,p) + OQ) ■
Тогда
v1(x, p) = ( Щ4 + o(^ (z$(x, p) + o(= ^(x, p) + o(,
,0/ 2—
Ас (р) ' ~ ) ' ~ \Р2 ' ~ \р2
V2(х,р) = £А(;(рр) + ^ (*с(х,р)+ °(р2)) = *(х,р) + р
и тем самым (13) и (14) при з =0 установлены.
Используя теперь вместо г2(х,р) асимптотики для г%\х,р) (з = 1, 2) придем аналогично к (13) и (14) при з = 1, 2. □
Лемма 3. Если д(х) е С[0,1] и р е 7п, то
(д,г1) = (д1(С)с08 сое 2пп^) - (д1(С)в1п ят2^^, (15)
(д, - г?) = -—-— [(д2 (С) с0й дС, эш 2пп^) + (д2 (С) вт д^ те 2пп^)], (1б)
2пп + д
где д1 (С) = д(С) + } ^(т, С)д(т) ¿г, д2 (С) = д(№(С, С) - } Щ (т, С)д(т) ¿г, р = 2пп + д.
1
1
Доказательство. Имеем:
5
(#,¿1) = (& сое р£ К1 (С,т) сое рт^т) = (У1 (С), сое рС) 0
и отсюда следует (15). Формула (16) получается аналогично. Лемма 4. Если #(х) е С[0,1] и р е 7п, то
(g,z2) =
1
[(gs(C) cos sin 2n^C) + (дз(C) sin cos 2n^C)],
2пп + р
- ¿0) = Т^-1-42 [(#4(С)сой соЭ2ПпС) - (#4(С)81П вт 2ПП^)] ,
(2пп + р)2
гае #з(С) = #(С) + / К(т, С)#(т) ¿т, #4(С) = С) + } ^2,5(т, С)#(т) ¿т.
5 5
Доказательство. Формула (17) получается так же как и (15). Далее, так как К2(С, 0) = 0, то имеем:
z2 - Z0 = — cos pxK2(x, x) +— cos ptK2 t(x, т) dr.
□
(17)
(18)
Отсюда приходим к (18).
□
Лемма 5. Обозначим через ^(x) функции cos x или sinx. Пусть f(x) е L2[0,1] и f(x,p) = = f (x)^(px), где p е y0 и e«(p) = (f (x, p),^(2nnx)). Тогда справедлива оценка
«2 1
Е n(р)к с
Vi.,
n2 '
(19)
где постоянная C не зависит от ni, n2 и р е y0.
Доказательство. По неравенству Коши - Буняковского имеем:
«2 1 «2 1 £ n|A'х £ n2
n = «l \ n = «l
По неравенству Бесселя получаем:
«2
Е ien (p)i ^ с ii/(x,p)y2,
E ien (p)i
(20)
(21)
где || ■ || — норма в Х2[0,1]. Но ||/(х, р)|| ^ С||/1|, где С не зависит от р е 70. Тем самым из (20) и (21) получаем (19). □
3. ИССЛЕДОВАНИЕ «о (х,£)
Лемма 6. Для м0(х,^) имеет место формула
it0(x,í) = -J + Е / /(RaPi)cospídA,
^ | A|=r 7™ '
(22)
где pi = R°0g.
x
n=«l
n=«l
n=«l
Доказательство. Действительно, по лемме 1, взяв RA вместо RA и у1 вместо у, имеем:
I (RAAy1)cos ptdA = I --R°Ag1 cos ptdA,
J J A - до
Y Y
где g1 = (L1 - доE)у. Но g1 = (L0 - доE)R°0g = g. Значит, имеет место (6) c заменой RA на RA и у на у1. □
Лемма 7. Собственные значения оператора L0 есть
An = p02 = 4п2 п2 (п = 0,1, 2,...). При этом A0 — простое собственное значение, а остальные A0 двукратные.
Доказательство. Все p0, отличные от нуля, есть корни Ао(p) = 1 -cosp. Но 1 -cosp = 2sin2(p/2). Отсюда p0 = 2nn и все они двукратны при п = 0. Докажем, что 0 есть простое собственное значение оператора L0. Действительно, если у" = 0, то y(x) = c0+c1x. Найдем c0 и c1 из граничных условий. Из у(0) = 0 следует c0 = 0, тогда y(x) = c1x. Отсюда y'(x) = c1, т. е. y(x) = x есть собственная функция для собственного значения A0 = 0. Покажем, что присоединенных функций для A0 нет. В противном
x
случае существует y(x) такое, что L0y = x, т. е. -у'' = x. Отсюда y(x) = c0 + c1x - f (x - t)t dt. Для
о
1
y(x) должны выполняться условия: у(0) = 0, у'(0) = у'(1). Но тогда c0 = 0, c1 = c1 - J tdt, чего
быть не может. □
Лемма 8. Для RA f имеет место формула
RA f = T-^—f f (i)cos p(1 - i) dj ^ +i ^^ii f (J) dj. (23)
1 - cos p p p
0 0
Утверждение леммы следует из теоремы 4. Теорема 7. Имеем место формула
те
n
n=
u0(x, t) = u0(x, t), (24)
где u0(x,t) = 2x(у1, 1), u0o(x,t) = an sin p0ox cos p0ot + bnx cos p°nx cos p0ot - bnt sin p0ox sin p0ot (n = 1, 2,...), an = 4((1 - £)у1 (J), sinp°nJ), bn = 4(у1®, cos p°nJ).
Утверждение теоремы следует из (22) и (23) по теореме вычетов.
Лемма 9. Для функций u°0(x,t) справедливы формулы:
u°n (x,t) = 2(vn (x + t)+ vn (x - t)) (n = 0,1,...), (25)
где v0(x) = u0(x, t) (u0(x, t) не зависит от t), vn(x) = an sinp0x + bnx cos p0x.
Доказательство. Для n = 0 равенство (25) очевидно. Пусть n ^ 1. Тогда, используя формулы для умножения тригонометрических функций, имеем:
u0o(x, t) = 2 {an(sin p°0(x + t) + sin p°0(x - t)) + bnx(cos p°0(x + t) + cos p°0(x - t))-
-bnt(cos p°0(x - t) - cos p°0(x + t))} = - {[an sinp°0(x + t) + bn(x + t) cos p°0(x + t)] +
2
+ [an sinp0(x - t)) + bn(x - t) cos p0(x - t)]} = ~(vn(x + t) + vn(x - t)).
□
Обозначим
те
F(x) = ^^ vn(x), F1 (x) = ^^ an sin p°nx, F2(x) = ^^ bn cos p°nx. 01
Лемма 10. Функции F(х), ^(х), F2(х) непрерывно дифференцируемы на оси (-го, го), ^(х) = -^1 (-х), (1 + х) = (х), ^г(х) = Т2(-х), х) = F2(х). При этом ^(х) = р 1 (х) при
х е [0,1].
Доказательство. Интегрируя два раза по частям формулы для ап и Ьп с учетом краевых условий на р 1(х), имеем:
_ ап 7 _ ап
«П = -02, = ,
рп рп
где через ап обозначаем различные числа ап, лишь бы ^ |ап|2 < го. Отсюда следует гладкость ^(х), ^1(х), ^2(х), если воспользоваться неравенством Коши-Буняковского. Свойства симметрии и периодичности (х) и ^2(х) очевидны. То, что F(х) = р1 (х) при х е [0,1], следует из разложения функции р1(х) (как функции из области определения оператора ) по собственным и присоединенным функциям оператора £0 в силу регулярности краевых условий. □
Лемма 11. Если х е [0,1], а т — любое целое число, то
F (т + х) = тр1 (1 — х) + (т + 1)р1 (х). (26)
Доказательство. Имеем F(х) = F1 (х) + xF2(х) при х е (—го, го). Пусть теперь х е [0,1]. Тогда
F(1 - х) = Fl (1 - х) + (1 - х^2 (1 - х) = (х) + (1 - х^(х) = = -F (х) + F2(x) = -р1 (х) + F2 (х).
Но F(1 - х) = р1(1 - х). Значит, F2 (х) = р1 (х) + р1 (1 - х). Поэтому
F (т + х) = F1(m + х) + (т + х^2 (т + х) = F1 (х) + (т + х^2(х) = = F(х) + mF2(х) = р1(х) + т(р1(х) + р1 (1 - х)) = тр1 (1 - х) + (т + 1)р1(х). □
Следствие. F(х) е С2(-го, го).
Доказательство. Имеем: -р'/(х) - р0р1 (х) = - р''(х) + д(х)р(х) - р0р(х). Так как р1 (0) = р(0) = р''(0) = 0, то и р1'(0) = 0. Поэтому
F'' (т + 0) = тр1'(1) + (т + 1)р1' (0) = тр1' (1), F''(т - 0) = (т - 1 + 1 - 0) = (т - 1)р1'(0) = тр1'(1) = тр1'(1),
т.е. F''(m + 0) = F''(m - 0).
Если же х = т, то F''(х) непрерывны в силу (26). □
Теорема 8. Функция и0(х,£) есть классическое решение эталонной задачи (1)-(3), т.е. когда д(х) = 0 и вместо р(х) берется р1 (х).
Доказательство. Имеем:
ио(х,г) = 2^(х + £) + F(х - £)]. (27)
Поэтому по следствию леммы 11
д2 и0 (х,£) д2и0 (х,£) д£2 дх2
Выполнение начальных и граничных условий следует из (27) и законности почленного дифференцирования ряда (24) по х. □ Результаты этого параграфа получены совместно с В. В. Корневым.
4. ИССЛЕДОВАНИЕ «2 (х,£)
По теоремам 3 и 4 с учетом того, что Мр/ и Мр0/ есть целые функции по А для и2(х,^), получаем следующее представление:
м2(х,£) = ^^ ап(х, £),
те^тео
где
ап(х,£) = I 2рС°8 ^ (х,р)(у,^1) + -У2(х,р)(у,^2) - (х,р)(у,^°)) - ^0(х,р)(у,г2)] ф. (28) 2п У р2 - ро
7™
(дз(С) с0й дС, вт^п^ (д4(С )соэ дС, соя2ппС),
(дз(С) вт дС, те 2п^С), -(д4(С) вт дС, вШ2"^).
Тогда по леммам 3 и 4 имеем:
(д,г1 )= Рп(д) + рп (д), (д,^2 ) =
1
(д,г1 - гС) =
1
2пп + д
[вп (д)+ вп (д)] ,
2пп + д
(д, г2 - гС) =
[вп (д)+ вп(д)] ,
(2пп + д)2
[вп (д)+ вп (д)]
Тогда, используя также лемму 2 и оценки
^(х,р) = С(р3-1), ^ (х,р) = 0(р),
получим из представления (29) оценку
(29)
Лемма 12. Ряды аЩх (х,Ь), Х^ аЩ^ (х,Ь) (з = 0,1, 2) сходятся абсолютно и равномерно по х е [0,1] и Ь е [—Т, Т] при любом фиксированном Т > 0.
Доказательство. Имеем:
J(х, р) = Л(х, р) + J2 (х, р), где J(х, р) — выражение в квадратных скобках в (28),
^(х,р) = (У1(х,р) - (х,р))(д,х1) + (У2(х,р) - у2(х, р))(д, г2), J2(x,р) = ^с(х,р)(д,^1 - г?) + (х,р)(д,г2 - г?).
Обозначим через вп(д) любой из функционалов:
(д1(С) с°э ДС, гов 2п^С), - (д1(С) вт дС, эш 2п^С), (д2 (С) С0Э ДС, с0й 2п^С), - (д2(С) эш дС, вт 2п^С),
(30)
J2j) (х,р) = 0(п-2) [|вп(д)| + |вп(д)|] + 0(п-7-1)
|2пп + д|
[|вп(д)| + |вп (д)|] +
+0(п-1)
|2пп + д|
[|вп (д)| + |вп(д)|]+ 0(п )
|2пп + д|2
[|вп (д)| + |вп(д)|] = 0(п-2 Рп (д)), (31)
где вп(д) = Х^ 1вп(д)| и оценки 0(...) равномерны по х е [0,1] и д е 7С. Отсюда 1
^ I (х,Ь)| =
0(п-3 )вп (д) Ж
70
Если 3 = 0,1, то отсюда |а^ (х,Ь)| = 0(п 2) и тем самым ряды Х^ |а„ Х (х, Ь)| сходятся равномерно по х е [0,1] и ь е [-Т,Т].
Далее из (31) следует оценка
к2 Х^(х,Ь)| = 0 (п-1Рп(д)) Ш,
70
и поэтому по лемме б
п2 Г {
Е ^(х,Ь)| = 0 (
п=п1 70 V
Е I 1^д1 = 0
А
п2 1 I
Е I ,
где оценка 0(...) равномерна по х е [0,1] и Ь е [-Т,Т]. Тем самым утверждение леммы получено для рядов Х^ ап(х,Ь).
Аналогично (даже проще) исследуются ряды Х^ аЩ(х,Ь)
□
1
1
1
1
п = п\
п = п\
5. КЛАССИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ (1)-(3)
Теорема 9. Формальное решение м(х,£) задачи (1)-(3) есть классическое решение при р(х) е е С2[0, —] и выполнении условий (4).
Доказательство. Для формального решения (5) в силу формулы (9) и аналитичности Мр/ по лемме 1 имеем формулу
и(м) = -2—/ Х-"р0(Дл^)с°8рЫА2—т/
|\|—„ те^тео ^
|Л|=г
У1 (х, р)(& ¿1) + У2(х, р)(& ¿2) А - ро
с°Б рЫА. (32)
Отсюда в силу оценок (31) по леммам 2-4 получаем, что ряд (32) и ряды, получающиеся из него почленным дифференцированием один раз по х и I, сходятся абсолютно и равномерно (см. также доказательство леммы 12). Таким образом, в этом случае процедура ускорения сходимости рядов не требуется. Поэтому и(х, £) удовлетворяет начальным и граничным условиям.
Далее, в силу теоремы 2 и леммы 12 формальное решение м(х,£) дважды непрерывно дифференцируемо. Докажем, что м(х,£) удовлетворяет уравнению (1). Обозначим через М оператор д2 д2
-. Тогда по теореме
Ми0 (х,£) = 0.
М=
д£2 дх2
Далее,
М«1 (х, £) = -
2—
|Л|=г
—-— М ((Длу)(х) С°Б р£) ¿А = М А - ро 2—
(
\
(ДЛр)(х) с°б р£ ¿А
= [ М ((Длр)(х)с°8р*) ¿А = ^
2— ^ 2— J
|Л|=г |Л|=г
\Л|=г
(Дл р)(х) с°б р£ ¿А.
/
Теперь имеем:
1 Г 1
М«2(х,^) = - — Е т-(х, р)(у, ¿1) + ^2(х, р)(у, ¿2)-
2— ^ 7 А - ро -^о(х, р)(у, ¿о) - г>°(х, р)(у, ¿о)] с°в р£) ¿А =
1 Г 1
— Е / А __ М([^ 1 (х,р)(у, ¿1)+ ^2 (х, р)(0, ¿2) С°Б р^^А.
2— ^^ У А - ро
те^тео ~
Но
Значит,
М (г^ (х, р) С°Б р£) = -г>"(х, р) С°Б р£ - р2(х, р) С°Б р£ = -^(х)^ (х, р) С°Б р£.
Ми2(х, £) = ^^ V [(ДЛр)(х) с°б рЫА, 2—г ^—' '
те^тео
т. е. Ми = -д(х)и(х, £).
Работа выполнена при финансовой поддержке РФФИ (проект № 13-01-00238).
Библиографический список
1. Стеклов В. А. Основные задачи математической физики. М. : Наука, 1983. 432 с.
2. Петровский И. Г. Лекции об уравнениях с частными производными. М. ; Л. : ГИТТЛ, 1953. 360 с.
3. Смирнов В. И. Курс высшей математики : в 4 т. М. : Гостехиздат, 1953. Т. 4. 204 с.
4. Ладыженская О. А. Смешанная задача для гипер-
□
болического уравнения. М. : Гостехиздат, 1953. 282 с.
5. Ильин В. А. Избранные труды : в 2 т. М. : Изд-во ООО «Макс-пресс», 2008. Т. 1. 727 с.
6. Ильин В. А. О разрешимости смешанных задач для гиперболических и параболических уравнений // УМН. 1960. Т. 15, вып. 2. С. 97-154.
7. Чернятин В. А. Обоснование метода Фурье в сме-
1
шанной задаче для уравнений в частных производных. М. : Изд-во Моск. ун-та, 1991. 112 с.
8. Крылов А. Н. О некоторых дифференциальных уравнениях математической физики, имеющих приложения в технических вопросах. М. ; Л. : ГИТТЛ, 1950. 368 с.
9. Бурлуцкая М. Ш., Хромов А. П. Смешанные задачи для гиперболических уравнений первого порядка с инволюцией // Докл. АН. 2011. Т. 441, № 2. С. 151-154.
10. Наймарк М. А. Линейные дифференциальные операторы. М. : Наука, 1969. 528 с.
11. Расулов М. Л. Метод контурного интеграла. М. : Наука, 1964. 462 с.
12. Вагабов А. И. Введение в спектральную теорию дифференциальных операторов. Ростов н/Д : Изд-во Рост. ун-та, 1994. 160 с.
13. Марченко В. А. Операторы Штурма - Лиувилля и их приложения. Киев : Наук. думка, 1977. 332 с.
About the Classical Solution of the Mixed Problem for the Wave Equation
A. P. Khromov
Saratov State University, 83, Astrakhanskaya str., 410012, Saratov, Russia, [email protected]
The classic solution of the mixed problem for a wave equation with a complex potential and minimal smoothness of initial data is established by the Fourier method. The resolvent approach consists of constructing formal solution with the help of the Cauchy-Poincare method of integrating the resolvent of the corresponding spectral problem over spectral parameter. The method requires no information about eigen and associated functions and uses only the main part of eigenvalues asymptotics. Krylov's idea of accelerating the convergence of Fourier series is essentially employed. The boundary conditions of the mixed problem can produce multiple spectrum and infinite number of associated functions in the spectral problem, thus making more difficult the analysis of the formal solution.
Key words: formal solution, spectral problem, resolvent, classical solution.
This work was supported by the Russian Foundation for Basic Research (project no. 13-01-00238).
References
1. Steklov V. A. Osnovnye zadachi matematicheskoi fiziki [The main tasks of mathematical physics]. Moscow, Nauka, 1983. 432 p. (in Russian).
2. Petrovsky I. G, Lectures on partial differential equations. Dover Publ. Inc., 1992, 245 p. (Rus. ed. : Petrovskii I. G. Lektsii ob uravneniiakh s chastnymi proizvodnymi. Moscow, GITTL, 1953, 360 p.).
3. Smirnov V. I. Kurs vysshei matematiki [A Course of Higher Mathematics : in 5 vol., vol. 4]. Moscow, Gostekhizdat, 1953, 204 p. (in Russian).
4. Ladyzhenskaya O. A. Smeshannaia zadacha dlia giperbolicheskogo uravneniia [Mixed problem for a hyperbolic equation]. Moscow, Gostekhizdat, 1953, 282 p. (in Russian).
5. Il'in V. A. Izbrannye trudy [Selected works : in 2 vol.]. Moscow, OOO «Maks-press», 2008, vol. 1, 727 p. (in Russian).
6. Il'in V. A. The solvability of mixed problems for hyperbolic and parabolic equations. Rus. Math. Surv., 1960, vol. 15, iss. 1, pp. 85-142.
7. Chernyatin V. A. Obosnovanie metoda Fur'e v sme-shannoi zadache dlya uravnenii v chastnykh proiz-vodnykh [Justification of the Fourier method in a mixed problem for partial differential equations]. Moscow, Moscow Univ. Press, 1991, 112 p. (in Russian).
8. Krylov A. N. O nekotorykh differentsial'nykh uravneniiakh matematicheskoi fiziki, imeiushchikh prilozheniia v tekhnicheskikh voprosakh [On some differential equations of mathematical physics with applications in technical matters]. Leningrad, GITTL, 1950, 368 p. (in Russian).
9. Burlutskaya M. Sh., Khromov A. P. Initial-boundary value problems for first-order hyperbolic equations with involution. Doklady Math., 2011, vol. 84, no. 3, pp. 783786.
10. Naymark M. A. Lineinye differentsial'nye operatory [Linear differential operators]. Moscow, Nauka, 1969, 528 p. (in Russian).
11. Rasulov M. L. Metod konturnogo integrala [The method of the contour integral]. Moscow, Nauka, 1964, 462 p. (in Russian).
12. Vagabov A. I. Vvedenie v spektral'nuiu teoriiu differentsial'nykh operatorov [Introduction to the spectral theory of differential operators]. Rostov-on-Don, Rostov Univ. Press, 1994, 106 p. (in Russian).
13. Marchenko V. A. Sturm - Liouville Operators and Applications. Kiev, Naukova Dumka, 1977, 332 p. (in Russian).