УДК 519.999
ОБРАТНАЯ ЗАДАЧА С ИНТЕГРАЛЬНЫМ УСЛОВИЕМ ПЕРЕОПРЕДЕЛЕНИЯ ДЛЯ ГИПЕРБОЛИЧЕСКОГО
УРАВНЕНИЯ
© 2014 А.Е.Савенкова1
В статье рассмотрена обратная задача определения правой части гиперболического уравнения с интегральным условием переопределения. Доказана теорема о существовании обобщенного решения.
Ключевые слова: гиперболическое уравнение, обратная задача, интегральное условие переопределения.
Введение
Обратные задачи определения правой части дифференциального уравнения возникают при математическом моделировании некоторых физических процессов в том случае, когда помимо решения уравнения требуется восстановить действие внешних источников. Обратные задачи для уравнений с частными производными различных типов исследовались во многих работах [1-3]. Отметим здесь некоторые из них, близкие к тематике данной статьи [4-6].
В обратных задачах вместе с начальными и граничными условиями, характерными для той или иной прямой задачи, задается дополнительная информация, необходимость которой обусловлена наличием неизвестных коэффициентов или правой части уравнения. Дополнительная информация, которая называется условием переопределения, может быть представлена в различных формах. Например, если известно значение искомого решения в определенный момент времени, то это дополнительное условие называют финальным переопределением. Однако часто такая информация поступает в усредненном виде, и тогда ее удобно представить как интеграл от искомого решения. Именно такое условие переопределения рассматривается в статье.
1. Постановка задачи
Рассмотрим в области Qт = (0,I) х (0, Т) уравнение
ии(х, Ь) — ихх(х, Ь) + с(х, Ь)и(х, Ь) = р(х)Н(х, Ь) (1.1)
хСавенкова Алеся Евгеньевна ([email protected]), кафедра математической физики Самарского государственного университета, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
и поставим для него следующую задачу: найти пару функций (и,р), удовлетворяющих уравнению (1.1), начальным данным
и(х, 0)=0,щ (х, 0) = 0, (1.2)
граничным условиям
и(0,г) = 0,и(1,г) = 0 (1.3)
и условию переопределения
I Н (г)и(х,г)аг = 5(х). (1.4) Jо
Прежде чем сформулировать основной результат, проделаем некоторые преобразования и докажем вспомогательные утверждения, нужные для обоснования разрешимости поставленной задачи. Обозначим
а(х) = I н(г)Н(х,г)аг. Jо
Пусть и £ С2^т)П С 1(Ят), Р £ С[0,1] и выполняются равенства (1.1)-(1.4). В этом случае пару функций (и,р) будем называть классическим решением задачи. Будем предполагать выполненными следующие условия:
Н(г) £ С2 [0, Т], Н(Т) = Н'(Т) = 0, ст(х) > сто > 0 Ух £ [0,1], (А)
6 £ С2 [0,1], ¿(0) = 6(1) =0. (В)
Умножим уравнение (1.1) на Н(г) и проинтегрируем по г от 0 до Т. После элементарных преобразований получим
1 Гт
р(х) = —-—- I" I г(х,г)и(х,г)А - 5''(х)], (1.5)
ст(х)1.]о
где г(х,г) = Н''(г) + с(х,г)Н(г).
Пусть теперь (и,р) удовлетворяет уравнению (1.1), и выполняются условия (1.2), (1.3), (1.5). Интегрируя (1.1), умноженное на Н(г), от 0 до Т, получим, принимая во внимание условия Н(Т) = Н (Т) = 0, равенство
Г т а2 гт гт
г(х,г)и(х,г)аг --¡-2 Н(г)и(х,г)аг = р(х) Н(г)Н(х,г)аг, ио ах ./о и о
откуда в силу (1.5)
¿2 Гт
Н(г)и(х,г)А - б''(х) = 0. (1.6)
а х ,] о
Из граничных условий (1.3) следует, что
[ Н(г)и(0,г)аг = 0, [ Н(г)и(1,г)аг = 0. (1.7)
оо
Уравнение (1.6) можно записать так:
а2
а2
( Н(г)и(х,г)аг - 5(х)) = 0.
ах2 о
Интегрируя его с учетом условий (1.7), получим
i H(t)u(x, t)dt = S(x), J 0
т. е. выполняется условие (1.4). Таким образом доказана
Лемма 1. Если выполняются условия (A) и (B), то задачи (1.1)—(1.4) и (1.1)—(1.3), (1.5) эквивалентны.
Доказанная лемма позволила разработать подход к обоснованию разрешимости поставленной задачи, который заключается в следующем.
Сначала мы докажем существование и единственность обобщенного решения задачи (1.1)—(1.3), (1.5), а затем покажем, что при выполнении некоторых условий на входные данные оно будет решением задачи (1.1)—(1.4). Обозначим
W1(Qt) = {v : v е W}(Qt), v(x, T) = 0}.
Введем понятие обобщенного решения задачи (1.1)—(1.3), (1.5). Для этого сначала выведем тождество, следуя известной процедуре [12, с. 93]: умножим (1.1) на v е W^Qt) и проинтегрируем по Qt. Получим
/'T /Л /Л z'T /Л
/ / ( —utvt+uxvx+cu)dxdt = / ф(x)v(x, 0)dx+ / p(x)h(x,t)v(x,t)dxdt. (1.8)
0 0 0 0 0
Определение. Обобщенным решением задачи (1.1)—(1.3), (1.5) будем называть пару функций (u,p) таких что u е W^Qt),Р е L2(0,l), u(x, 0) = 0, выполняется условие (1.5), и для всех v е W2i(QT) справедливо тождество (1.8).
Рассмотрим теперь задачу (1.1)-(1.3), предполагая, что p(x) известна. Обобщенное решение u(x,t) е W^Qt) этой задачи существует, единственно, и справедлива оценка [12, с. 209].
II u Ww^Qtc II P(x) Ilb2(0,0 •WHlzQt) (1.9)
Опуская доказательства единственности, которое проводится так же, как в [12, с. 209], остановимся на процедуре доказательства существования решения для того, чтобы уточнить значение постоянной c, входящей в (1.9). Пусть
c(x,t) е C (Qt ), h(x, t) е C (Qt ). (1.10)
Пусть функции wk(x) е C2[0,l],wk(0) = wk(l) = 0,wk(x) линейно независимы и образуют полную в W21 (0,1) систему. Будем искать решение задачи (1.1)—(1.3) в виде
m
um(x,t) = ^2 ck(t)wk(x)
k=1
из соотношений
d W
p(x)h(x t)wj
/ \Uimwj(х) + (х) + cumWj(х)в,х] = (х)&. (1-11)
■10 ■! 0
Напомним, что сейчас мы считаем функцию р(х) известной. Соотношение (1.11) вместе с условиями ск (0) = ск (0) = 0 представляют собой задачу Коши для систем линейных обыкновенных дифференциальных уравнений второго порядка, которая в силу условий (1-10) однозначно разрешима, что полностью определяет последовательность приближенных решений {пт(х,1)}.
Получим теперь оценку, нужную как для обоснования возможности перехода к пределу при т —> то, так и для нахождения условий разрешимости обратной
задачи. Умножим (1.11) на Cj (t), просуммируем по j от 1 до m и проинтегрируем по t от 0 до т, получим
[ [ \utt (x,t)ur[i(x,t) + um(x,t)uZt(x,t)+ cum(x,t)urm(x,t)dxdt] =
Jo Jo
f-T rl
p(x)h(x, t)u?(x, t)dxdt.
ио ио
Интегрируя по частям в левой части последнего равенства, получим
/■I
и"Ч 1 )) + (их
11 \(иТ(х,т))2 + №(х,т))\1х =
= [ ( р(х)Н(х,г)ит(х,г)йх& - [ [ сит(х,г)ит(х)йх&. (1.12)
о о о о
Заметим, что в силу граничных условий имеет место представление
ит(х,г) = (
о
из которого следует неравенство
(ит(х,г))2 < ¿[ (ит(х,г))2йх. (1.13)
о
Так как начальные условия однородны, то аналогично получаем, что (ит(х,т))2 < т[ (Ут(х,г))2йь,
Jо
откуда
f (um(x,T))2dx < т[ f (u™(x,t))2dxdt. (1.14)
Jo Jo Jo
Оценим правую часть (1.12): | i i p(x)h(x,t)urmi(x,t)dxdt\ ^ — f f p2(x)dxdt +--- ( i (u'ln(x,t))2)dxdt,
Jo Jo 2 Jo Jo 2 Jo Jo
где ho = max \h(x,t)\, Qt
if i cum(x,t))u™(x,t)dxdtl^ Co f f ((um(x,t))2 + (u™(x,t))2 )dxdt,
Jo Jo 2 Jo Jo
где Co = max\(x,t)\. С учетом полученных оценок и неравенств (1.13), (1.14) из
Qt
(1.12) следует неравенство
i [(v™(x,t))2 + (um(x,T))2 + (um(x,T))2]dx < ho i i p2(x)dxdt+
Jo Jo Jo
+C1 f f [(um(x,t))2 + (um(x,t))2 + (um(x,t))2]dxdt,
Jo Jo
r, Col2,
C1 = max{ho + Co + },
применив к которому лемму Гронуолла получим
i [(um(x,T))2 + (um(x,T))2 + (um(x,T))2]dx < hoeClT i i p2(x)dxdt.
Jo Jo Jo
Интегрируя полученное неравенство по т от 0 до Т, приходим к неравенству
II
которое влечет за собой
II um \\2wi(QT) < ¡hoc1(eClT — 1)|| p \\12{оЛ,
II u IlWi(QT) < lhoCi(eC1 T — 1)|| p I\lmy (1.15)
Заметим, что \ h(x,t) \\l2(qt ) = (Jo Jo h2 (x,t)dxdt)2 ^ hoVTl. Но тогда (1.9) можно записать следующим образом:
\\ u IIW1QT) < choV¥l\\ p \\L2(o,l).
Следовательно, c = \J~T[Ci(eClT — 1). Обозначим choVTl = C3.
Операторное уравнение. Рассмотрим соотношение (1.5). Введем оператор
1 Г
A(p) =-тт / r(x,t)u(x,t)dt, (1.16)
a(x) Jo
где u(x,t) — решение прямой задачи (1.1)—(1.3) с заданной функцией p(x).
Лемма 2. Оператор A, определяемый формулой (1.16), действует из L2(0,l) в L2(0,l).
Доказательство. Рассмотрим
fi 1 fT
\ A(p) \\L2(oi) = r(x,t)u(x,t)dt)2dx.
Jo a(x) Jo
Так как r(x,t) = H (t) + cH(t), то в силу условий (A) и (1.10) r(x,t) G C(QT).
Тогда 3 Ro > 0 такое, что max | r(x,t) Ro. Применив неравенство Коши-Буня-
QT
ковского, получим для почти всех x G (0, l)
T T 2 f T r T I- T
( / r(x,t)u(x,t)dt) ^ r2(x,t)dt u2(x,t)dt ^ TR^ u2(x,t)dt.
Jo Jo Jo Jo
Так как по условию a(x) ^ ao > 0, то
TR2 T i TR2
\ A(p) \\L2(o I) < —2Г V2(x,t)dxdt < tro4M2l2(o l),
2 ao o o ao 2
откуда
II A(p) II
L2(o,l) ^ c4\\pIIL2(o,l), (1.17)
RpThp аол/h
Запишем соотношение (1.5) в виде операторного уравнения:
где c4 = RoTh c1(eClT — 1). Таким образом, A(p) G L2(0,l). Лемма доказана.
р - АМ = --Х-. (1.18)
Найдем условия, при выполнении которых оператор А(р) сжимающий. В силу неравенства (1.17) это свойство будет иметь место, если С4 < 1. Так как
RoTho т С4 =-гт с\(е - 1),
ао\I
то А(р) сжимающий, если
К0ТН0е1(еС1Т - 1) < . (1.19)
Таким образом, если справедливо условие (1.19), то операторное уравнение (1.18) однозначно разрешимо, и его решение р(х) £ Ь2(0,1). Доказано следующее утверждение.
Теорема 1. Если выполнены условия (А), (В) и справедливо неравенство (1.19), то задача (1.1)—(1.3), (1.5) однозначно разрешима.
Теперь приступим к доказательству разрешимости задачи (1.1)—(1.4). Для этого нам потребуется показать, что найденное решение задачи (1.1)—(1.3), (1.5) принадлежит не только пространству W2(Qт), но и W2(Qт)■
Лемма 3. Если выполняются условия теоремы 1, и, кроме того, ^ £ ),
то решение задачи (1.1)-(1.3), (1.5) принадлежит пространству W2(QT). Доказательство. Продифференцируем (1.11) по
/ {иш(х,1)^ (х) + итг(х,1)тп (х) + сит(х, (х) + е^пт(х^)ю(х)^йх Jо
Г1
= / р(х)Н^х,1)тз(х)йх.
Jо
ае
з
проинтегрируем по £ от 0 до т :
Умножим полученное равенство на Сз (£) и просуммируем по ] от 1 до т, а затем
/ / (ит,иа + + еити™ + елитит)йхйЬ = [ ( рН1и11йхйЬ. (1.20)
ио Jо Jо Jо
Проинтегрируем слагаемые, стоящие в левой части последнего равенства
гт р I 1 /• I 1 /• I
,т„,тл—П — 1 I„.т/п ^Х\2 1 /п\\2,
£ I ититйхль = 11 (итх т))2з.х - 2 ^ ытх о))2йх, ! [ иттитгйхж = \ ( (итт(х,т))2йх,
Jо Jо 2 и о
[ [ еититйхйь = -1 [ [ ен(ит)2йхйг +1 [ е(ит(х,т))2йх. и о Jо 2 Jо Jо 2 Jо
Тогда из (1.20) получим 1-1
и"Нг (х, т)) + (ихг(х,1)) + (еи
2 [(ит(х,т))2 + (ит(х,т))2 + (еит(х,т))2]йх =
= 1 [«т(х,о))2йх +1 [ [ а(ит(х,г))2йхйг+ 2 ■! 2 ■! о Jо
+ [ [ егит(х,г)(ит(х,г))2йхйг + ( [ р(х)н^х,г)итйхйг. (1.21) ио Jо Jо Jо
оо
,-Т !■ I ,-Т !■ I
(х)
о о о о
Рассмотрим (1.11), положив в нем £ = 0, умножим на ет(0), а затем просуммируем по ] от 1 до т:
I [(ит(х, 0))2 + ит(х, 0)ит(х, 0) + еит(х, 0)ии(х, 0)]йх
Jо
= [ р(х)к^х, 0)и%(х, 0)3х. (1.22)
Jо
Учитывая, что игт^(х, 0) = и'Ча(х, 0) = 0, получим, применив неравенство Коши к правой части (1.22), неравенство
I (и%(х, 0))Чх < 2 I р2ь1(х, 0)3х + 11 ит(х, 0))23х, откуда следует оценка
[ (и%(х, 0))23х < [ р2к2(х, 0)3х < р1, (1.23)
ио Jо
где р1 = к!\\р\\Ь2(01), к! = тах\к^(х,Ь)\.
Продолжим оценку. Оценим правую часть равенства (1.21) по абсолютной величине. Второе слагаемое оценим с помощью (1.14), к третьему и четвертому слагаемым применим неравенство Коши "с е".
\1 [ I <н(и?(х,ь))23х31 кръ 2 и о ио
\ [ ( вгит(х,г)(ит(х,г))23х3ь [ ( (ит(х,г))23х3ь+е I [ (ит(х,г))23х3г,
и о Jо 2е ■! 0 и 0 2 ■) 0 и 0
\ I I р(х)к^х,1)ит(х,Ь)3х3Ь 1/ I р2(х)к'2(х,1)3х31+—( / (и^(х,Ь))23х3Ь, Jо Jо 2 е ■! 0 -¡0 2 ■! 0 -¡0
с2 рт I
(ит(х,г))23х3г < р2,
2— 0 0
р2(х)к2(х,1)3х3Ь ^ рз,
¡■т Г1
22
где
2— 0 0
с5 = тах \ сг(х, Ь) \, Ят
2
с5 + С0е . сТ .
р2 = —2——1к0С1(вС1Т — 1)р1, 1
рз =2—р1.
Выберем е так, чтобы 2е < 1. Тогда (■I
1(ии 1 )) + uxt(x, 1 )ии
т0 / [(ит(х,Т))2 + и™(х,т)ит(х,т) + сипт(х,г)ии(х,т)]3х < р1 + р2 + рз,
0
0
где
т0 = 1 — е.
Отсюда
\\ ит и^ < ра, \\ ит \\ь,
где ра = р1 + р2 + рз.
Оценим и™. В силу полученных оценок можем записать тождество:
/ / (иЦV + и'тух + сиу)3х3Ь = / рку3х3Ь. и0 Jо Jо Jо
Представим у(х,Ь) = Ф(1)У(х), где Ф(Ь) £ С(0,Т),У £ (0,1), тогда для почти всех Ь £ [0, Т]
/•'___ г1
/ (и^V + и^Ух + сиЪ)йх = рНУв,х, ■10 ■! 0
откуда
г1 г1
/ иттухв,х = (рн—ит — сит)ъйх.
00
о
Так как функция У произвольная из пространства (0,1), то последнее равенство означает, что существует обобщенная производная
ихх = сит + ит — рН.
Но тогда полученное решение и из Ш2^т) принадлежит Ш2^т )■ Поэтому уравнение (1.1) выполняется для почти всех (х,Ь) £ Qт, и выполняются условия (1.2),
(1.3), (1.5).
Покажем, что выполняется условие (1.4). Пусть р*(х) £ Ь2(0,1) является решением операторного уравнения (1.18), и*(х,Ь) — обобщенное решение из пространства^! ) прямой задачи (1.1)-(1.3) с функцией р* в правой части уравнения (1.1). Положим 5*(х) = ^ Н(Ь)и*(х,Ь)&. В силу свойств решения
о
и* (х,Ь) и Н(Ь) 6*(х) £Ш1 (0,1). Повторяя рассуждения, приведенные при доказательстве леммы 1, получим
1 Г т
р* (х) = —1—г( I г(х,г)и*(х,г)(й — 5*" (х)). (1.24)
°(х) Уо
Но, с другой стороны, р* является решением (1.18):
1 Гт
р* (х) = —-—-( I г(х,г)и*(х,гш — 5 (х)). (1.25)
°(х) ./о
Тогда, вычитая из (1.24) равенство (1.25), получим
5* (х) = 5 (х),х £ (0,1).
Так как 5*(0) = 5(0) = 0,5*(1) = 5(1) = 0, то 5*(х) = 5(х), следовательно, и* удовлетворяет условию (1.4), стало быть, является решением задачи (1.1)—(1.4), понимаемой как функция, принадлежащая пространству Ш2^т), удовлетворяющая условиям (1.2)—(1.4) в смысле равенства функций в Ь2, и почти всюду уравнению (1.1).
Литература
[1] Cannon J.R., Lin Y. Determination of a control parameter in a parabolic partial differential equation // J. Austral. Math. Soc. Ser. B, 1991. Vol. 33. № 2. P. 149--163.
[2] Cannon J.R., Lin Y. Determination of a parameter p(t) in some quasi-linear parabolic differential equations // Inverse Problems, 1998. V. 4. № 1. P. 35-45.
[3] Камынин В.Л. Об обратной задаче определения младшего коэффициента в параболическом уравнении при условии интегрального наблюдения // Мат. заметки. 2013. Т. 94. Вып. 2.
[4] Сафиуллова Р.Р. Обратная задача с неизвестным составным внешним воздействием при составном переопределении // Мат. заметки ЯГУ. 2006. Т. 13. Вып. 2.
[5] Павлов С.С. Обратная задача восстановления внешнего воздействия в многомерном волновом уравнении с интегральным переопределением // Мат. заметки ЯГУ. 2011. Т. 18. Вып. 1.
[6] Павлов C.C. Нелинейные обратные задачи для многомерных гиперболических уравнений с интегральным переопределением // Мат. заметки ЯГУ. 2011. Т. 18. Вып. 2.
[7] Камынин В.Л. Об обратной задаче определения правой части в параболическом уравнении с условием интегрального переопределения // Мат. заметки. 2005. Т. 77. Вып. 4.
[8] Камынин В.Л. Об обратной задаче определения старшего коэффициента в параболическом уравнении // Мат. заметки. 2008. Т. 84. Вып. 1.
[9] Прилепко А.И., Костин А.Б. Об обратных задачах определения коэффициента в параболическом уравнении. II // Сибирский мат. журнал. 1993. Т. 34. № 5.
[10] Амиров А.Х. К вопросу о разрешимости обратных задач // Сибирский математический журнал. 1987. Т. XXVIII. № 6.
[11] Денисов А.М. Обратная задача для гиперболического уравнения с нелокальным краевым условием, содержащим запаздывающий аргумент // Труды Института математики и механики УрО РАН. 2012. Т. 18. № 1.
[12] Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973. 407 с.
References
[1] Cannon J.R., Lin Y. Determination of a control parameter in a parabolic partial differential equation // J. Austral. Math. Soc. Ser. B. 1991. V. 33. № 2. P. 149-163.
[2] Cannon J.R., Lin Y. Determination of a parameter p(t) in some quasi-linear parabolic differential equations // Inverse Problems. 1998. V. 4. № 1. P. 35-45.
[3] Kamynin V.L. On inverse problem of determining the lower-order coefficient in parabolic equations with integral observation // Matematicheskie Zametki. 2013. V. 94. № 2. P. 207-217.
[4] Saffullova R.R. Inverse problem with unknown complex external influence at complex overdetermination // Matematicheskie Zametki YaGU. 2006. V. 13. № 2. P. 79-94.
[5] Pavlov S.S. Inverse problem of recovering the external influence for many-dimensional wave equation with integral overdetermination condition // Matematicheskie Zametki YaGU. 2011. V. 18. № 1. P. 81-92.
[6] Pavlov S.S. Nolinear inverse problems for many-dimensional hyperbolic equations with integral overdetermination // Matematicheskie Zametki YaGU. 2011. V. 18. № 2. P. 128-153.
[7] Kamynin V.L. On the inverse problem of determining the right-hand side in a parabolic equation under an integral overdetermination condition // Matematicheskie Zametki. 2005. V. 77. № 4. P. 522-534.
[8] Kamynin V.L. On the inverse problem of determining the leading coefficient in parabolic equation // Matematicheskie Zametki. 2008. V. 84. № 1. P. 48-58.
[9] Prilepko A.I., Kostin A.B. On inverse problems of determining a coefficient in a parabolic equation. II // Sibirsky Matematichesky Zhurnal. 1993. V. 34. № 5. P. 147-162.
[10] Amirov A.Kh. Solvability of the inverse problems // Sibirsky Matematichesky Zhurnal. 1987. V. XXVIII. № 6. P. 3-11.
[11] Denisov A.M. Inverse problem for a hyperbolic equation with nonlocal boundary condition containing a delay argument // Trudy Instituta Matematiki i Mekhaniki UrO RAN. 2012. V. 18. № 1. P. 139-146.
[12] Ladyzhenskaya O.A. Boundary problems of mathematical physics. M.: Nauka, 1973. 407 p.
Поступила в редакцию 19/77/2014; в окончательном варианте — 19/77/2014.
INVERSE PROBLEM WITH INTEGRAL OVERDETERMINATION CONDITION FOR A HYPERBOLIC EQUATION
© 2014 A.E. Savenkova2
In the paper, we study an inverse problem for a hyperbolic equation with integral overdetermination condition. The existence of a generalized solution is proved.
Key words: hyperbolic equation, inverse problem, integral condition of overdetermination.
Paper received 19/77/2014. Paper accepted 19/77/2014.
2Savenkova Alesya Evgen'evna ([email protected]), the Dept. of Equations of Mathematical Physics, Samara State University, Samara, 443011, Russian Federation.