CC BY
8
2017 Математика и механика № 49
УДК 517. 977.5
Б01 10.17223/19988621/49/2
Г.Ф. Гулиев, Ю.С. Гасымов, Х.Т. Тагиев, Т.М. Гусейнова
ОБ ОБРАТНОЙ ЗАДАЧЕ НАХОЖДЕНИЯ ПРАВОЙ ЧАСТИ ВОЛНОВОГО УРАВНЕНИЯ С НЕЛОКАЛЬНЫМ УСЛОВИЕМ
Рассматривается задача определения правой части волнового уравнения с нелокальным условием. Эта задача приводит к задаче минимизации некоторого функционала, построенного с помощью дополнительной информации. Для новой задачи выводится необходимое и достаточное условие оптимальности.
Ключевые слова: обратная задача, волновое уравнение, нелокальные условия, условие оптимальности.
В последнее время обратные задачи для дифференциальных уравнений интенсивно изучаются. Отметим, что такие задачи возникают в самых разнообразных областях математики, геофизики, сейсмологии, астрономии, экологии и т.д. [1]. В данной работе рассматривается подход к решению одной обратной задачи для волнового уравнения. Поиск неизвестной правой части уравнения сводится к задаче минимизации функционала, построенного с помощью дополнительной информации. В результате получаем градиент функционала и условие оптимальности.
1. Постановка задачи
Для цилиндра 0>Т =Пх(0,Т) рассматривается краевая задача
д2м
--Дм =-Э(х, 0, (х, /) е О, ; (1)
д/2
дм( х, 0)
м(х,0) = ф0(х),-= Ф1(х), хеП ; (2)
д/
дм ду
= IК(х, у)м(у,/)ёу , (х,/) е Бт . (3)
Здесь Пей" - ограниченная область с гладкой границей дП; Бт = дПх(0,Т) -боковая поверхность цилиндра 0>т; V - внешняя нормаль к границе дП; Ф0 (х) е (П), ф2 (х) е Ь2 (П); К(х, у) е Ь2 (П х П) - заданная функции, а $(х,/) е Ъ2(0,Т) - неизвестная функция. Для того чтобы определить $(х,/), воспользуемся дополнительной информацией
м(х, Т) = я(х), х еП , (4)
где я (х) е ¿2(П) - заданная функция.
Б
Приводим эту задачу к задаче оптимального управления, т.е. на решениях задачи (1) - (3) минимизируем функционал
1 2
Л(9) = -|(( х, т; 9) - я (х))), (5)
2 о
где и(х,Т;9) является решением задачи (1) - (3), которое соответствует функции 9(х,t). Функцию 9(х,t) назовем управлением. Если мы найдем управление 9(х, t), которое доставляет функционалу (5) нулевое значение, тогда дополнительное условие (4) выполняется.
Отметим, что при каждом фиксированном управлении 9(х, t) е Ь2 (вТ) краевая
задача (1) - (3) имеет единственное обобщенное решение из W21 (вТ) [2].
2. О разрешимости задачи (1) - (3), (5)
Теперь рассмотрим следующую задачу: при каких условиях
М Л(9) = 0? (6)
9е^2 (вт)
Пусть у0 (х) - заданная функция из Ь2 (О), такая, что
|у0(х)и(х,Т;9)Сх = 0, У9е 12(вТ). (7)
О
Мы хотим выяснить, будет ли отсюда следовать у 0 (х) = 0 ? Введем функцию W(х, t) как решение задачи:
д2W
--ДW = | К(5, х^(5, t)сЬ , (х, t) е вт ; (8)
дt зо
W(х,Т) = 0, Ш(X,Т) = у0(х), хеО ; (9)
дt
Ж
дv
= 0, (х, 0 е 5Т . (10)
Как и в работе [2], можно показать, что (8) - (10) имеет единственное обобщенное решение из класса W21 (вт) и это решение обладает свойствами
W(х, t) е С([0, Т], W21 (О)), ^д^ *) е С([0, Т], Ь2 (О)).
В силу определения обобщенного решения задачи (1) - (3) имеем: при t = 0 выполняется условие и(х,0) = ф0(х) и интегральное тождество:
т ( ди дп ^ т
||1---- + УиУп I СхС-| | п(х,t)| К(х,у) и(у,t) dydsdt-
0 ОЧ 3* 3* ) 0 дО О
т
-1ф1(х)п(х,0)Сх =| |9(х,*)п(х,*)СхС* (11)
О 0 О
для любой функций п е W21 (вТ), п(х, Т) = 0 , где Уи =
( ди ди ди ^ Чдх1 'дх2 ""'дх„ )
5.
Т
В силу определения обобщенного решения задачи (8) - (10) имеем: при / = Т выполняется условие Ж (х, Т) = 0 и интегральное тождество:
Т ( дЖ дФ ^ Т
III---+ УЖУФ I <ыг-I I Ф(х,/)1 К(х,у)Ж(у,/) dydsdt-
0 ПЧ д/ д/ ) 0 дП П
-I дЖ(х,0) Ф(х,0)<х = ^^х)ф(х,Т)<х</, (12)
п д/ п
для любой функции Ф е Ж2: (0>Т).
Теперь в (11) за функцию п возьмем Ж, а в (12) за функцию Ф возьмем м, из (11) вычтем (12), тогда имеем
дЖ(х 0)
I-!—ф0(х)<х - х)Ж(х,0)<х -
п д/ П
Т
-119(х,/)Ж(х,/)<х</ + I у0(х)м(х,Т;9)<х</ = 0 .
0 П П
Если учесть условия (7), то получим
дЖ(х 0) Т
I-х)<х - I ф^х)Ж(х,0)<х -119(х,/)Ж(х,/)<х</ = 0 . (13)
П д/ П 0 П
Если соотношение (13) записать для произвольных 9^х,/) и 92(х,/), то из полученных двух равенств следует, что
Т
11 (9 (х, /) - 92 (х, /))Ж(х,/)<х</ = 0,, 92 е ¿2 (бТ).
0П
Отсюда, в свою очередь, следует, что Ж = 0 в Qт . Значит, в силу (8) у0 (х) = 0 в
Таким образом, в силу теоремы Хана - Банаха [3] получаем, что
М ,/0(3) = 0.
&еЬ2^т)
Если образ Ь2(^т) при отображении 9^-м(х,Т;9) замкнут в Ь2^т), то возможно существует такой элемент 90 (х, /) е Ь2 (QT), что
М J0(&) = = 0.
»еL2(Qт) 0 0 0
В задаче (1) - (3), (5) минимизирующий элемент 9(х, /) е ¿2®т), вообще говоря, не единственный. Теперь рассмотрим задачу минимизации функционала
Ja (9) = J0(9) +—Т I (9(х,/) -ю(х, /))2<х</ (14)
2 0 П
в выпуклом замкнутом множестве Па< е ¿2®т) при ограничениях (1) - (3), где ю(х,/) е Ь2^т) - заданная функция, а > 0 - заданное число.
Тогда, в силу известных результатов [5], в задаче (1) - (3), (14) существует единственный минимизирующий элемент.
3. Вычисление дифференциала функционала (14) и необходимое условие оптимальности
Теперь покажем, что функционал (14) дифференцируем в Ь2(вТ). Берем два допустимых управления 9,9 + 59е ПаС . Соответствующее решение задачи (1) -(3) обозначим через и(х, ^ 9) и и(х, ^ 9 + 59).
Пусть 5и(х, t) = и(х, ^ 9 + 59) - и(х, ^ 9). Ясно что 5и(х, t) является обобщенным решением краевой задачи
-д-5и-Д5и = 59(х,t), (х,t) е О, ; (15)
З* 2
5и (х,0) = 0, З5идх,0) = 0, х еО ; (16)
З5и
з =| К(х, у)5и (у, t)Су, (х, t) е ^ , (17)
дv
т.е. для любой функций п е W21(Qr), п(х, Т) = 0 выполняется интегральное тождество
Т
т i д5и dn ] т
| JI---- + V5u Vn I dxdt- J | n(x,t)J K(x,y)5u(y,t) dydsdt =
0 QV dt dt ) 0 3Q Q
т
= J J 59(x, t)n(x, t)dxdt.
0Q
(18)
Пусть у - обобщенное решение из W2, (QT) сопряженной задачи
д 2у
-Ду = J K(5, x)y(§, t)ds , (x,t) £ QT ; (19)
dt2 9q
у (x, T) = 0, x,T) = u (x, T; 9) - g (x), x £Q ; (20)
dt
dy dv
= 0, (x, t) £ ST . (21)
То есть для любой функции Ф £ W2, (QT) выполняется интегральное тождество
t i дУ дФ ^ T
JJI —--+ VyVФ I dxdt - J J Ф (x, t) J K (x, y) у (y, t) dydsdt -
0 QV dt dt ) 0 3Q Q
- j ду(x,0) ф(x,0)dx- J (u(x,T;9) -g(x)^(x,T)dx = 0. (22)
q dt q
Поскольку смешанная задача (19) - (21) является линейной относительно у (x, t), то эта задача в пространстве W2, (QT) имеет единственное решение[2]. Теперь вычислим приращение функционал (14). Ясно, что
5J— (9) = Ja (9 + 59) - Ja (9) = J (u( x, T; 9) - g (x))5u( x, T )dx +
Q
T 1 T
+aJ J(9(x,t)-ra(x,t))59(x,t)dxdt + - J(5u(x,T))2dx + —J J(59)2dxdt. (23)
0 Q 2 Q 2 0 q
S
т
Если в (18) положить п(х,/) = у(х,/), а в (22) Ф(х,/) = 5м(х,/) и вычесть полученные соотношения, то имеем
т
I (м(х, Т; 9) - я(х))5м(х, Т)<х =! 159(х, /)у(х, /)<х</.
П 0 П
Тогда, учитывая это равенство в (23), получим
т т
Ыа(9) = аI I (9(х,/)-ю(х,/))59(х,-1189(х,/)у(х,/)<х< + Я , (24)
0 П 0 П
1 т
где Я = -1 (5м (х, Т))2 <х + -11 (59)2 <х</.
2 П 2 0 П
Теперь оценим остаточный член Я, входящий в (24). Покажем, что
т 2
|Я| < с11 (59)2<х</ . (25)
0П
Для этого покажем, что
НЖ^т) < с1 К^). (26)
Здесь и в дальнейшем через с будем обозначать различные постоянные, не зависящие от оцениваемых величин и от допустимых управлений. Применяя метод Галеркина, из (15) - (17) получим
г д25uN (х, /) д5uN (х, /) ™ д5uN (х, /) д25uN (х, /) п д/2 д/ П i=1 дхг- дхг- д/
- I Шм^IК(х,у)5мы(у,/)= I59д5UNЛ.
дП д/ П П д/
N
Здесь 5м (х,/) - приближения Галеркина, т.е. 5м (х,/) = 2 С^ (/)фк (х), фк (х) -
к=1
базис в Ж21(П).
Интегрируя по / от 0 до /, получаем
I
д 5м .
+ У5м
д/ ) )
<х =
/
= 21 I д5мМ(х/) IК(х,у)5м1'(у,/)<у<5</ + 2II 59д5мд(^-/)<< (27)
0 дП д/ П 0 П д/
Преобразуем интеграл по боковой поверхности Бт следующим образом:
д5uN (х, /)
I I-К(х, у)5м (у, = i1 + i2 + ^
0 дП д/ П
^ ^ ^ т^, ч д5uN(у,/) , , , где г1 =-] 15м (х, /) ] К (х, у)-——1 ,
дП 0 П д/
^ = I 5uN (х, /) IК(х, y)5мN (у, /)<у<5,
П П
П
3 = - | 5un(х,0) | K(х, у) 5uN(y,0)dyds = 0 .
ЗО
Пользуясь неравенством | < c | (|W| + |VW|)Сх [4] и затем неравенством
зо
Коши - Буняковского, получаем
д5N (у, t)
| | 5и (х, t) | K (х, у)
0 зо о 3*
< c11 ((5uN(х,t))2 +|v5uN(х,t)| )СхС* + с|1
0 О ' ' 0 О
СуС^С* ' ( д5uN(у, t)
д*
СуС* ; (28)
| 5uN (х, 01К (х, у^и1* (у, ^уСз
дО
< сЧ/((5uN(х,t))2 +|v5uN(х,*)|2|Сх 12 ( | (5uN(х,t))2dx-2 . (29)
Введем обозначение
1 (t) ((5и )2 +
о
(55^ 1 ч )
+ У5им )Сх.
Ясно, что
| (5uN (х, *))2 Сх < | (*)С* ,
О 0
где
(*) = |
о
12
д
+ У5и'
'(х, * )|2
Сх.
Из (29) следует
ц < с(?ы (*))21 | (*)С* I .
При условиях на данные задачи и учитывая (28) и (30), из (27) имеем
(30)
ZN (*) < с| ZN (*) С* + | ZN (*)С* + с(ZN (*))2 I | ZN (*)С* |2 + с| | (59)2 СхС*.
0 0 Ч 0 ) 00
Обозначим тах ZN (5) = 1 (*), тогда
0<5<*
(*) < (с + 2*^ + (*) + с11 (59)2СхС*.
0
Определяя из условия (с + 2*1 )*1 = 2 , получим
(*) < с|| (59)2СхС*.
0
о
о
Продолжая этот процесс, за конечное числа шагов получаем неравенство:
д5 N V 2
П((5м^2 +1^ + |У5м^ )<х < с|5^|L,(Qт), V/ е [0,Т],
V /
из которого следует
||5uN| |21 < с| |59|| 2Г (31)
ii ) 11 ХХТ2(°Т)
В силу (31), из последовательности {5мN} можно выбрать подпоследовательность, сходящуюся слабо в Ж2: (От) к некоторому элементу 5м (х, /) е Ж2, (QT).
Поскольку норма в гильбертовом пространстве слабо полунепрерывна снизу, то отсюда следует, что для предельной функции 5м (х, /) справедлива оценка (26) [7].
По теореме выложения Ж2: (QT) с Т2 (П) [4] получаем, что
Nx,Т)Ё2(П) < с| 15м(^0|). (32)
Сопоставляя соотношения (26) и (32), имеем
1М^Т)||Т2(П) < с||59(x,/^).
Отсюда и из выражения Я следует справедливость оценки (25).
Тогда из равенства (24) и из оценки (25) следует, что функционал Ja (9) дифференцируем в L2(QT) и его дифференциал и градиент определяются выражениями
т
Л(9),59^ ) = 11 (а(9-ю)-у(х,/;9))59(х,/)<х</ (33)
2 т) 0 П
и Jа (9) = а(9-ю) - у (х, /; 9) (34)
соответственно.
Можно показать, что отображение 9^ Ja (9) из Па< ^ L2 (QT) является непрерывным.
Из (34) следует
J'a (9 + 59) - J'a (9) = а59-5у( х, /). (35)
Из выражения (35) имеем
т 1
иа (9 + 59) - Ja (9)||т ю < с( II ((59)2 +(5у)2)<х</)2. (36)
2 ^ ' 0 П
Ясно, что 5у(х, /) = у(х, /; 9 + 59) -у(х, /; 9) является решением задачи:
д 25у
-Д5у = I К(§, х)5у(§,/<, (х,/) е Qт ; (37)
д/2 дП
5у(х,Т) = 0, д5у(х,Т) =5м(х,Т), хеП ; (38)
д/
д5у
= 0, (х,/) е Бт . (39)
Б
Т
Как и для задачи (15) - (17), относительно 5у(x, t) имеем оценку
IN ti(Br ) * С115м ( x, T С □)• (40)
По теореме выложения W (QT ) с L2 (Q) [4] получаем, что
IN^T)IL2( Q) * clMx,0|W21(QT)• (41)
Сопоставляя соотношения (26), (40) и (41), имеем
J(9 + 59)-J(9)||l2(Qt) *C|K2(QT)• (42)
Из этого неравенства следует, что при ||59|| ^ 0
11 L2 (QT )
jа (9+59) - jа (9)| L2(Qt ) ^о.
Следовательно, отображение 9^ J'a(9) из Uad ^L2(QT) является непрерывным.
Теорема. Пусть выполняются вышеналоженные условия на данные задачи (1) - (4). Тогда для оптимальности управления 9, = 9,(x, t) е Uad в задаче (1) - (3), (14) необходимо, чтобы выполнялось неравенство
T
Ц (а(9, -ю) - у(x, t;9,))(9-9,)dxdt > 0 (43)
0 Q
при всех 9 е Uad .
Доказательство. Согласно доказанным утверждениям, функционал Ja (9) непрерывно дифференцируем по Фреще на L2(QT) и его дифференциал в точке 9(x,t) е Uad определяется равенством (33). В силу теоремы [6, с.28] на элементе 9, е Uad необходимо и достаточно выполнение неравенства
J (9), 9-9,}> 0 при всех 9е Uad .
Отсюда и из (34) следует справедливость равенства (43). Теорема доказана.
ЛИТЕРАТУРА
1. Кабанихин С.И. Обратные и некорректные задачи. Новосибирск: Сибирское научное издательство, 2009. 457 с.
2. Кожанов А.И., Пулькина Л.С. О разрешимости краевых задач с нелокальным граничным условием интегрального вида для многомерных гиперболических уравнений // Диффе-ренц. уравнения. 2006. Т. 42. № 9. С. 1166-1179.
3. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1981. 543 с.
4. Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973. 408 с.
5. Лионс Ж.Л. Оптимальное управление системами, описываемыми уравнениями с частными производными. М.: Мир, 1972. 414 с.
6. Васильев Ф.П. Методы решения экстремальных задач. М.: Наука, 1981. 400 с.
7. Tagiev H.T. Gradient of the functional in the optimal control problem with non-local conditions for the wave equation // Proceedings of the Institute of Mathematics and Mechanics NAS of Azerbaijan. V. XXXVII(XLV). 2012. P. 139-148.
Статья поступила 16.03.2017 г.
Guliyev H.F., Gasimov Y.S., Tagiyev H.T., Huseynova T.M.(2017) ON THE INVERSE PROBLEM OF FINDING THE RIGHT-HAND SIDE OF WAVE EQUATION WITH NONLOCAL CONDITION. Tomsk State University Journal of Mathematics and Mechanics. 49. pp. 16-25
DOI 10.17223/19988621/49/2
Recently, inverse problems for the differential equations have been intensively studied. Such problems arise in the various fields of mathematics, geophysics, seismology, astronomy, ecology, etc. In this paper, we propose an approach to solving the inverse problem for the wave equation. The search for the unknown right-hand side of the equation is reduced to the problem of minimizing the functional constructed using additional information. The gradient of the functional is calculated and the optimality condition is derived. In the cylinder QT =Qx(0,T), consider the problem
d2U
—--Au = 9(x,t), (x, t) e QT , (1)
dt
u( x,0) = 90( x), du (x,0) = x), x eQ , (2)
dt
du dv
= J K(x,y)u(y,t)dy , (x,t) e ST , (3)
ST Q
where Qe Rn is a bounded domain with a smooth boundary dQ , ST = 3Q x (0,T) is the laterial surface of QT, v is an outward normal to dQ, q>0(x) e ^(Q), x) e L2(Q), K(x,y) e L2(QxQ) are given functions, and 9(x,t) e L2(QT) is the unknown function. To determine 9(x, t), we use the following additional information:
u (x,T) = g (x), x eQ , where g (x) e L2(Q) is a given function. (4)
The problem is reduced to the following problem: minimize the functional
1 2 J0(9) = - J(u( x,T; 9) - g (x) )2dx (5)
2 Q
subject to (1)-(3), where u (x,T; 9) is a solution of problem (1)-(3) corresponding to 9( x, t) which is called a control. The solvability of problem (1)-(3), (5) is proved. Consider the functional
a t 2
J a (9) = J 0(9) + -J J(9( x, t) - rn( x, t) )2dxdt. (6)
2 0 Q
Then, the differential of this functional is calculated and the following theorem is proved:
Theorem. Under the considered conditions, for the optimality of the control 9, =9,(x,t)eUad in the problem (1)-(3), (6) it is necessary that the inequality
T
J J (a(9, -ro) - x,t; 9, ))(9 - 9, )dxdt > 0 (7)
0Q
is fulfilled for all 9eUad .
Keywords: inverse problem, wave equation, nonlocal conditions, optimality condition.
GULIYEV Hamlet Farman (Doctor of Physics and Mathematics, Baku State University, Baku, Azerbaijan). E-mail: hkuliyev@rambler.ru
GASIMOV Yusif Soltan (Doctor of Physics and Mathematics, Azerbaijan University, Institute of
Mathematics and Mechanics of Azerbaijan National Academy of Sciences, Baku, Azerbaijan)
E-mail: yusif.gasimov@au.edu.az
TAGIYEV Hikmet Tahir (Candidate of Physics and Mathematics, Baku State University, Baku,
Azerbaijan). E-mail: hitagiyev@gmail.com
GUSEYNOVA Tunzale Maharram (Candidate of Physics and Mathematics, Azerbaijan State
Pedagogical University, Baku, Azerbaijan). E-mail: tunzalemustafayeva@mail.ru
REFERENCES
1. Kabanikhin S.I. (2009) Obratnye i nekorrektnye zadachi [Inverse and ill-posed problems] Novosibirsk: Sibirskoe nauchnoe izdatel'stvo.
2. Kozhanov A.I., Pul'kina L.S. (2006) On the solvability of boundary value problems with a nonlocal boundary condition of integral form for multidimensional hyperbolic equations. Diff. Equat. 42. pp. 1233-1246. https://doi.org/10.1134/S0012266106090023
3. Kolmogorov A.N., Fomin S.V. (1981) Elementy teorii funktsiy i funktsional'nogo analiza [Elements of the theory of functions and functional analysis]. Moscow: Nauka.
4. Ladyzhenskaya O.A. (1973) Kraevye zadachi matematicheskoy fiziki [Boundary value problems of mathematical physics]. Moscow: Nauka.
5. Lions J.L. (1968) Controle optimal de systèmes gouvernés par des équations dérivées partielles. Paris, Dunod, Gauthier-Villars.
6. Vasilyev F.P. (1981) Metody resheniya ekstremal'nykh zadach [Methods for solving extreme problems]. Moscow: Nauka.
7. Tagiev H.T. (2012) Gradient of the functional in the optimal control problem with non-local conditions for the wave equation. Proceedings of the Institute of Mathematics and Mechanics NAS of Azerbaijan. XXXVII(XLV). pp. 139-148.
