2017 Математика и механика № 45
УДК 517.95
DOI 10.17223/19988621/45/4
Р.К. Тагиев, Р.А. Касумов
ОБ ОПТИМИЗАЦИОННОЙ ПОСТАНОВКЕ КОЭФФИЦИЕНТНОЙ ОБРАТНОЙ ЗАДАЧИ ДЛЯ ПАРАБОЛИЧЕСКОГО УРАВНЕНИЯ С ДОПОЛНИТЕЛЬНЫМ ИНТЕГРАЛЬНЫМ УСЛОВИЕМ
Рассматривается оптимизационная постановка коэффициентной обратной задачи для параболического уравнения с дополнительным интегральным условием. Исследованы вопросы корректности оптимизационной постановки обратной задачи. Доказана дифференцируемость целевого функционала и найдена формула для его градиента. Установлено необходимое условие оптимальности в виде вариационного неравенства.
Ключевые слова: оптимальное управление, параболическое уравнение, интегральное граничное условие, условие оптимальности.
Обратные задачи для уравнений с частными производными могут быть сведены к задачам оптимального управления соответствующими системами. Например, обратные задачи для уравнений тепломассообмена могут рассматривается как задачи оптимального управления тепловыми режимами технических объектов. При этом управляющие воздействия обычно входят в коэффициенты уравнений тепломассообмена или граничные условия для них. Эти воздействия должны быть определены таким образом, чтобы удовлетворить условия. Обычно эти критерия качества составляются на основе дополнительных информации и их называют функционалами невязки, или целевыми функционалами.
Одним из основных типов обратных задач для уравнений с частными производными являются задачи, в которых подлежат определению коэффициенты уравнений или величин, в них входящих, по некоторой дополнительной информации. Такие задачи называются коэффициентными обратными задачами для уравнений с частными производными.
В работе А.Н.Тихонова [1] предложена идея использования методов теории оптимального управления для решения обратных задач. В работах [2-7] и др. обратные задачи об определении коэффициентов соответствующих уравнений с частными производными сводились к задачам оптимизации для этих уравнений с управлениями в коэффициентах, т.е. исследовались оптимизационные постановки коэффициентных обратных задач. Во многих этих работах дополнительные условия, по которым подлежат определению коэффициенты уравнений, являются локальными. Оптимизационные постановки коэффициентных обратных задач с дополнительными нелокальными условиями мало изучены.
В данной работе рассматривается оптимизационная постановка одной коэффициентной обратной задачи для параболического уравнения с целевым функционалом, соответствующим дополнительному интегральному условию. Исследованы вопросы корректности оптимизационной постановки обратной задачи. Доказана дифференцируемость по Фреше целевого функционала и найдено выражение для его градиента. Установлено необходимое условие оптимальности в виде вариационного неравенства.
1. Постановка задачи
Пусть управляемый процесс описывается в QT = {(x,t)eR2:0<x<£,0<t<T} следующей начально-краевой задачей для линейного параболического уравнения
ut-(k(x,t)ux)x + q(x,t)u = f (x,t), (x,t)e QT ; (1)
M|t=0 =ф(x), 0 < x < I; (2)
Uxlx=0 = UxU = 0, 0 < t < T . (3)
Здесь t, T > 0 - заданные числа, f (x,t)e L2 (QT), ф(x)e W2, (0,1) - заданные функции, k(x,t), q(x,t) - неизвестные коэффициенты, u(x,t) = (k(x,t),q(x,t)) - управление, u = u (x, t ) = u (x, t; u) - решение задачи (1) - (3) - состояние процесса, соответствующее управлению u = u(x,t).
Введем множество допустимых управлений
V = {u(x,t) = (k(x,t),q(x,t)) e H = W1 (QT)x L2 (QT ) :0 < v < k(x,t) <
lkx (^t )l < i Ik (xt) < < q0 < q (^t) < qi п.в.на qt Ь (4)
где ^>v>0, ц1,ц2 >0, q1 >q0 >0 - заданные числа.
Поставим следующую коэффициентную обратную задачу типа оптимального управления: среди всех допустимых управлений и(x, t) = (k (x, t), q (x, t))eV найти управление и» (x, t) = (k„ (x, t), q„ (x, t)) e V, минимизирующее функционал
t i
J(u) = J IK (x,t)u (x,t; u)dx- E(t)
0 0
где K (x, t )e Lx (QT), E (t )e L2 (0, T) - заданные функции, причем |K (x, t)|<|a3, п.в. на QT , |3 = const > 0 . Эту задачу ниже будем называть задачей (1) - (5).
Обозначения используемых в работе функциональных пространств соответствуют принятым в [8, с. 12-17]. Ниже положительные постоянные, не зависящие от оцениваемых величин и допустимых управлений, обозначаем через Mi (i = 1,2,...).
Под решением краевой задачи (1) - (3), при каждом фиксированном управлении u = u(x,t) e V , будем понимать обобщенное решение из V21,0 (QT), т.е. функцию u = u (x, t) = u (x, t; u) e V21,0 (QT), удовлетворяющую для всех П = П(x, t) e W (QT), n(x,T) = 0 интегральному тождеству
i
J^-u^ + k(x,t)uxnx + q(x,t)un]dxdt = ^(x)n(x,0)dx + JJ f (x,t)ndxdt. (6) QT 0 Qt
При сделанных предположениях краевая задача (1) - (3) имеет единственное обобщенное решение u = u (x, t; u) из V21,0 (QT) при каждом фиксированном u = u(x,t) e V и справедлива оценка [8, с. 181-189]
К-(Qt) = IuQt = 0mTlu (x,t; u))L2 (0,i)+lluxllL2 (Qt) < M1 [llfllL2 (Qt )+Hl2 (0,/)] . (7)
2
dt , (5)
Более того, это решение принадлежит также пространству ^г22'1 (ОТ), удовлетворяет уравнению (1) при почти всех (х,t) е ОТ и справедлива оценка [8, с. 203-211]
ЬЬрць) < М1 [1М^'(о,/) +1 А\ь2(о )]. (8)
Используя оценки (8) теорем вложения [8, с.78; 9, с. 33] и рассуждая аналогично работе [10] можно показать, что для решения краевой задачи (1) - (3) справедлива также оценка
\\и\кЖ) +И16 (От ) * М 2 [||фЦ21(0,/) +1 Л\ь2 О)] . (9)
Из оценки (7) следует, что функционал (5) определен на V и принимает конечные значения. Отметим, что функционал (5) нелинеен и исследование его выпуклости весьма сложно.
Задача (1) - (5) тесно связана с коэффициентной обратной задачей, заключающейся в определении функций {и (х, /; и), к (х, t), q (х, t)}, удовлетворяющих условиям (1) - (4) и дополнительному условию
I
| К (х, ^и (х, ^ и)ск = Е (), 0 < t < Т . (10)
0
Функционал (5) является функционалом невязки в Ь2 (0, Т) соответствующей условию (10). Если в задаче (1) - (5) окажется, что существует управление и„( х, t ) = ( к„( х, t), q»( х, t ))е V, такое, что 3 (о„) = ^ {,/ (и): ие V } = 0, то это управление решает обратную задачу (1) - (4), (10).
Задача (1) - (5) является задачей оптимального управления для параболического уравнения с управлениями в коэффициентах. Такие задачи при других целевых функционалов исследованы в работах [10-12] и др.
2. Корректность постановки задачи
Теорема 1. Пусть выполнены условия, принятые в п. 1. Тогда множество оптимальных управлений задачи (1) - (5) К»={и„е V: 3 (о„) = ./„} не пусто,
слабо компактно в Н и любая минимизирующая последовательность {и„ } =
= {(к„ (х, t), qn (х, t))} с V функционала (5) слабо в Н сходится к множеству V,.
Доказательство. Множество V, определяемое равенством (4), выпукло, замкнуто и ограничено в рефлексивном банаховом пространстве Н , и поэтому оно слабо компактно в Н [13,с. 49]. Покажем, что функционал (5) слабо в Н непрерывен на множестве V . Пусть и(х, t) = (к(х,t),q(х,t))е V - некоторый
элемент, {и„ (х, t)} = {(к„ (х,t), q„ (х,t))} с V - произвольная последовательность, такая, что и„ ^ и слабо в Н , т.е.
к„ (х, tк (х, t) слабо в W21 (ОТ ); (11)
qn (х, t) ^ q (х, t) слабо в Ь2 (ОТ ). (12)
Из (11) и компактности вложения W21 (ОТ) ^ Ь (ОТ) [8, с. 75] следует, что
к„ (х, t) ^ к (х, t) сильно в Ь (ОТ), (13)
где г1 > 2 - произвольное конечное число.
Кроме того, в силу однозначной разрешимости краевой задачи (1) - (3), каждому управлению ип е V соответствует единственное решение ип = и (х, ип) из
^22,1 (От) задачи (1) - (3) и справедлива оценка
,2,1(0, )< М 3 (п = I2, ■■■)..
' п|Ц2,1(От
(14)
т.е. последовательность {ип} равномерно ограничена в пространстве ^22,1 (От).
Тогда из (14) и компактности вложения ^22,1 (От) ^ Ь(От) [9, с. 33] следует, что из последовательности {ип} можно извлечь подпоследовательность {ип }, такую , что
ип (х, t) ^ и (х, t) слабо в W22,1 (0Т) и сильно в Ь (От ), (15)
где г2 > 2 - произвольное конечное число, и (х, t) - некоторый элемент
W22,1 (От ).
Покажем, что и (х, ^ = и (х, ^ и), (х, ^е От, т.е. и (х, ^ является решением задачи (1) - (3), соответствующим управлению ие V . Ясно, что справедливы тождества
Я[-ип„П + кпп (x,t)иптхПх + Япт (x,)ип]^ =
От
I
= |ф(х)п(х,0)ах +Ц / (x,t )ПСХС1,
0 ОТ
Уп = п(х,t) е W21 (0Т), п(х,t) = 0. Используя соотношения (13) - (15) и неравенство Гельдера, имеем
(16)
Я кпп и t) иппхПхСхС - Я к ^ 1:) ихnxdxdt
О, От
ЯК (х1: )- к ^1: ^ иппх Пх^Л
От
\\кп - к\.....
Лк(x, 1:)[иппх -их~]ЦхСхС:
От
Л к (x, t)[иппх - их ^х^1
От
" " ^ 0,
Ь(От)« пп 1Ьб(ОтГ 2Шт при п ^го .
Кроме того, используя соотношения (12), (15), получаем
(17)
Ццп (х,t)ип цСх&-Цц(х,t)ипйхск
О, От
Я[Яп (х,О-Ц(х,ип^СхС < ц
От
Я Цпп ^t )ипп - и ]пСхС1
От
ип -и НпНт \ +
"п 1Ь2(ОТ Г 2(ОТ)
Ц[дп (х,t)- ц(х,t)]ипёхЛ
От
^ 0,
(18)
при п ^ГО .
Тогда, переходя к переделу при т ^го в равенстве (16) и учитывая соотношения (15), (17), (18), получаем, что функция и (х, t) удовлетворяет тождеству
(6). Отсюда и из включения и (х, ^е W22,1 (ОТ) следует, что и (х, ^ = и (х, ^ и) т.е. функция и (х, t) является решением задачи (1) - (3), соответствующим управлению и е V.
Используя единственность решения задачи (1) - (3), соответствующего управлению ие V, нетрудно показать, что соотношение (15) справедливо с функцией
и (х,t) = и (х,^ и) не только для под последовательности {и„ }, но и для всей последовательности {и}, т.е.
ип (х, t) = и„ (х,^ и) ^ и (х, t) = и (х,^ и)
слабо в W22,1 (ОТ ) и сильно в Ьг (ОТ).
(19)
Рассмотрим теперь функционал цели 3(и), определенный формулой (5). Используя (5), нетрудно убедится, что
3 (и„)-3 (и)| =
|К(х,t)и„ (х,t)dx-Е()
0
|К(х,t)[и„ (х,t)-и(х,t)]]
Ь2 (0,Г)
Ь2 (0,т;
| К (х, t) и (х, t) dx - Е ()
0
I
|К(х,t)и„ (х,t)dx-Е()
ь2(0,т;
Ь2(0,Т)
| К (х, t) и (х, t )dx - Е (t)
Ь2 (0,Т)]
(20)
Используя неравенство Коши - Буняковского и условие |К (х, t)| < ц3 п.в. на ОТ, имеем
| К (х, t )[и„ (х, t)-и (х, t )]dx
Ь2 (0,Т)
| |К (х, t)[и„ (х, t)- и (х, t)] с1х
2 II
dt I <
\т /1
<Ц1 |К2 (х,t)dx II ||и„ (х,t)-и(х,t)|2 dx
dt I <ц3 ||и„ - и|
Ь2 (От
I Т I
|К(х,t)и„ (х,t)(Лх-Е() =Ц |К(х,t)и„ (х,t)dx-Е()
Ь2 (0,Т)
0 0
dt !> <
гт гI
Н |К2 (x,t)с!х\и„2 (x,t)с!х + Е2 (t) ЛI <л/2 [Ц3 ||иЛь2(вт)+ +(0,т)] ;
^0 [0
I
|К(x,t)и(x,t)(^х-Е(t) <721>3 ||и||ь2(0,т)+11Е11ь2(0,Т)] .
Ь2 (0,Т)
Учитывая эти неравенства в (20), получаем оценку
VК)-3(и)1<^3 к -4ь2(дт(||ип\\ь2Шт)+|Нь2(дт)) + 211ЕИ2(0,т)] . Тогда, используя оценку (14) и соотношение (19) получаем, что 3 (ип ) ^ 3 (и) при п , т.е. целевой функционал (5) непрерывен на V в слабой топологии пространства Н . Кроме того, множество V является выпуклым замкнутым ограниченным множеством в гильбертовом пространстве Н . Тогда применяя результат из [13, с. 49] получаем, что задача (1) - (5) корректно поставлена в слабой топологии пространства Н , т.е. справедливы все утверждения теоремы 1. Теорема 1 доказана.
3. Дифференцируемость целевого функционала и необходимое условие оптимальности
Для задачи (1) - (5) введем сопряженное состояние у = у (х, t) = у (х, и) как решение задачи
Уt +(к (хt )у х) х- ч (хt )у =
= -2К (х, t)
|К(§,t)и (§,П Е()
0
у|. = 0, 0 < х< I;
(х,t) е дт ;
(21)
(22)
У х 1х=0 =У х\х=1 = 0, 0 <t < т. (23)
Под решением краевой задачи (21) - (23), соответствующим управлению ие V, будем понимать функцию у = у( х, t ) = у( х, V, и) из V1,0 (йт), удовлетворяющую для всех п = п(х, t (йт), п( х,0) = 0 интегральному тождеству
{{[УП + к (х, t )ух пх + ч (х, t )уn]aXdt =
От
■2 Ц К (х, t) | К (§, t) и (§, V, о) dх - Е ()
йт
ndхdt.
(24)
Используя очевидное неравенство (а + Ь)2 < 2(а2 + Ь2) и неравенство Коши -Буняковского, имеем
ЦК2(х,t) |К(§,t)и(§,V,и)^-е(t)
йт
2
ёхЛ <
: 2 Ш К2 (х, К2 (§, t)d^ и2 (§, Г, о)^ ёхЛ + Ц К2 (х, ^Е2
1йт ь
йт
< 2^2^"1 иII12(йт)+ 4Е112(0,т).
Отсюда следует, что правая часть уравнения (21) является элементом пространства Ь2 (йт). Тогда из результатов работы [8, с. 181-189] следует, что для
каждого заданного ие V задача (21) - (23) имеет единственное решение из V,1,0 (ОТ) и справедлива оценка
1М1^0 О) < М4 [И Ь2 (еТ ) +1 Иь2 (0,Т)].
Учитывая здесь оценки (7), получаем
Ы^ (От ) < М5 [II А\Ь2 (вт) +ИфЬ2 (0,.) +11Е11Ь2 (0,Т)] (25)
Кроме того, решение задачи (21) - (23) принадлежит также пространству W22,1 (От), удовлетворяет уравнению (21) при почти всех (х, t) е От и справедлива оценка [8, с. 203-211]
ЫЦ^О) <М6 [IIЛ\Ь2(От) +Пфк(0,/) +1ЕЬ2(0,т)] .
Используя эту оценку, теорему вложения [8, с. 78; 9, с. 33] и рассуждая аналогично работе [10] получаем оценку
1МЬ„(ОТ ) +1^х11ь6 (От) < М 7 [11 /1Ь2 (От ) + 11фк(0,/) +1 1Е1Ь2 (0,Т)] . (26)
Введем еще одну вспомогательную краевую задачу для определения функции ю = ю( х, t ) = ю( х, t; и) из условий:
-Юхх + Ю = ихУх , (x, t)е ОТ ; (27)
юх1х=0 =юх1х=^ = 0, 0 < х < £ ; (28)
ю^=0 = ю^=т = 0, 0 < t < Т. (29)
Под решением краевой задачи (27) - (29), при заданном иеV, будем понимать функцию ю = ю(х,t) = ю(х,V,и) из (От), удовлетворяющую интегральному тождеству
Ц(сохпх + +®n)dxdt = Цихухndxdt, (30)
От От
при любой функции п = п(х,t) е (От ).
Из вложений их е Ь6 (От), ух е Ь6 (От) следует, что ихух е Ь3 (От). Тогда из результатов работы [14, с. 200] следует, что краевая задача (27) - (29), при заданном ие V, имеет единственное обобщенное решение из (От) и справедлива оценка
Нк(От) < М8 ИихУх\\ьф(От) < М8 11их11ь6(От) ИУх|1ь2(От) .
Учитывая здесь оценки (8) и (26), получаем оценку
ИЦО) <М9 [ИфЦ(0,/) +11 /1Ь2(от)][11ф1ь2(0,/) +1 /1Ь2(От) +ИЕ11ь2(0,Т)] . (31)
Теорема 2. Пусть выполнены условия теоремы 1. Тогда функционал (5) непрерывно дифференцируем на V в норме пространства Н и его градиент в точке ие V определяется равенством
3'(и) = (ю(х,V,и),и (х,V, и)у(х,V,и)), (х,t)е От . (32)
Доказательство. Пусть ueV - некоторое управление, Au = (AK,Aq)еH -произвольное приращение управления и, такое, что u + Aue V. Пусть
Au (х, t) = u (x, t; u +Au) - u (x, t; u), u (x, t) = u (x, t; и). Из условий (1) - (3) следует, что Au является решением краевой задачи
Aut -((k + Ak)Aux )x +(q + Aq)Au = (Akux )x -Aqu , (x, t)e QT ; (33) Au|t=0 = 0 , 0 < x < £ ; (34)
Auxlx=o =AuX=l = 0, 0 <t < T . (35)
Можно показать, что для решения задачи (34) - (35) справедлива оценка [8, с. 181-189]:
llAullv2i,0 (qt ) < Mw [llAkuxl L2 (qt ) +llAquIL2 (qt )]. (36)
Используя ограниченность вложений
W1 (Qt L4 (QT ), W22,1 (Qt ) ^ Lm (QT ) и оценки (8), получаем следующие неравенства:
\\Akux\\L2(Qt) < INIL4(Qt) IluxllL4(Qt) < M11 llAkllw2(QT) ,
IML,(qt ) <INL,(qt )llul„(qt) < M12 IMl,Q) .
Учитывая эти неравенства в (36), получаем оценку
11А
u\l ,1,0
Vj-° (Qt
< M„ II Au
IlH •
Приращение функционала (5) имеет вид
AJ (u) = J (u + Au)-J (u) =
t | e e 1
= 2 J IK (|, t)u (|, t; u)d|-E (t) IK (x, t)Au (x, t)dxldt-
(37)
t t
l-| |K(x,t)AK(x,t)dx
dt.
(38)
С помощью решений краевых задач (21) - (23), (27) - (29) и (33) - (35) преобразуем приращение (38). Для решения краевой задачи (33) - (35) справедливо
Ц[Дигу + (к + Ак)Дихух +(ч + Ад)Дuу]dхdt = -Ц(Дкихух + Дquу)dхdt (39)
йт йт
Если в тождестве (24) положим п = Ди и полученное равенство вычтем из (39), то получим
21-
IK (|, t) u (|, t; u)d E (t)
L 0
1 1
| K (x, t )Au (x, t )dx 1t =
= II (ux^x Ak - Aux^xAk + uyAq - AuyAq) dxdt.
Qt
Учитывая это равенство в (38), имеем
AJ(u) = U(uxyxAk + uyAq)dxdt + R ,
(40)
Qt
2
где
R = | |K(x,t)Ди(x,t)dx
2
dt -Ц(Дмху x Дк + Аu\Аq)dxdt. (41)
Qt
Полагая в тождестве (30) п = Дк и учитывая полученное равенство в (40), имеем
AJ(и) = Ц(юАк + raxД^ +ratAkt + ы\\Дд)dxdt + R . (42)
Qt
Используя теоремы вложения, можно показать [10], что первое слагаемое в правой части (42) при заданном ие V определяет линейный ограниченный функционал от Ди е H .
Теперь проведем оценку остаточного члена R , определяемого равенством (41). Используя неравенства Коши - Буняковского и Гельдера, имеем
T (l Y1 \
RI 41 i K 2 (x,t)dx iДи 2 (x,t)dx dt +11 ДиЛь2 (QT )l\xlL4 (QT )MIl4 (QT ) +
0V0 7Vo У
+НДи1^2 (qt ) IML(qt ) ^(qt ) - ^ I ML (Qt ) + M12 I Ml^q ) xlL, (qt ) INW'Q ) + +\\Ди\\ь2(Ят ) IMl„(qt )M\L2(Qt )- M131 MlH.
Учитывая в (42) эту оценку, заключаем, что функционал (1) дифференцируем по Фреше на V и справедлива формула (32). Непрерывность отображения и ^ J'(и) доказывается как в работе [10]. Теорема 2 доказана.
Необходимое условие оптимальности в задаче (1) - (5) устанавливает Теорема 3. Пусть выполнены условия теоремы 1. Тогда для оптимальности управления и» =(к»,q»)e V в задаче (1) - (5) необходимо, чтобы выполнялось неравенство
£[[ю»(к - k») + ro»x (kx - k„x ) + ю»( (kt - k»t ) + u»\»(q - q»)]dxdt > 0 (43)
Qt
для любого u = (k,q)eV, где u» = u(x,t;и»), \»=\(x,t;и»), ю»=ю(x,t;и») -решения задач (1) - (3), (21) - (23), (27) - (29) соответственно при и = и».
Доказательство. Множество V , определяемое равенством (4), выпукло в H . Кроме того, согласно теореме 2, функционал (1) непрерывно дифференцируем на V . Тогда в силу теоремы 5 из [13, с. 28] на элементе и» е V» необходимо выполнение неравенства ( J '(и»), и-и») > 0 при всех ие V. Отсюда и из (32) следует справедливость неравенства (42). Теорема 3 доказана.
ЛИТЕРАТУРА
1. Тихонов А.Н. О решении некорректно поставленных задач и методе регуляризации // ДАН СССР. 1963. Т. 151. № 3. С. 501-504.
2. Гласко В.Б. Обратные задачи математической физики. М.: Изд-во МГУ, 1984. 112 с.
3. Искендеров А.Д. О вариационных постановках многомерных обратных задач математической физики // ДАН СССР. 1984. Т. 274. № 3. С. 531-533.
4. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач. М.: Наука, 1988. 288 с.
5. Karchevsky A.L. Properties the misfit functional for a nonlinear one-dimensional coefficient hyperbolic inverse problem // J. Inverse III-Posed. Probl. 1997. V. 5. No. 2. P. 139-165. DOI: https://doi.org/10.1515/jiip.1997.5.2.139.
6. Кабанихин С.И., Искаков К.Т. Обоснование метода наискорейшего спуска в интегральной постановке обратной задачи гиперболического уравнения // Сиб. матем. журн. 2001. Т. 42. № 3. С. 567-584.
7. Кабанихин С.И., Даирбаева Г. Обратная задача нахождения коэффициента уравнения теплопроводности // Международная конференция «Обратные некорректные задачи математической физики», посвященная 75-летию академика М.М. Лаврентьева, 20-25 августа 2007 г., Новосибирск, Россия. С. 1-5.
8. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967. 736 с.
9. ЛионсЖ-Л. Управление сингулярными распределенными системами. М., 1987. 368 с.
10. Тагиев Р.К. Оптимальное управление коэффициентами в параболических системах // Дифференц. уравнения. 2009. Т. 45. № 10. С. 1492-1501.
11. Тагиев Р.К. Оптимальное управление коэффициентами квазилинейного параболического уравнения // Автоматика и телемеханика. 2009. № 11. С. 55-69.
12. Тагиев Р.К. Задача оптимального управления для квазилинейного параболического уравнения с управлениями в коэффициентах и с фазовыми ограничениями // Дифференциальные уравнения. 2013. Т. 49. № 3. С. 380-392.
13. Васильев Ф.П. Методы решения экстремальных задач. М.: Наука, 1981. 400 с.
14. Ладыженская О.А., Уральцева Н.Н. Линейные и квазилинейные уравнения эллиптического типа. М.: Наука, 1973, 576 с.
Статья поступила 07.09.2016 г.
Tagiev R.K., Kasumov R.A. (2017) ON THE OPTIMIZATION FORMULATION OF THE COEFFICIENT INVERSE PROBLEM FOR A PARABOLIC EQUATION WITH AN ADDITIONAL INTEGRAL CONDITION. Tomsk State University Journal of Mathematics and Mechanics. 45. pp. 49-59
DOI 10.17223/19988621/45/4
Let a controlled process be described in QT = {(x,t)e R2 : 0 < x < ^,0 < t < T} by the following initial-boundary value problem for a linear parabolic equation:
ut -(k(x,t)ux) + q(x,t)u = f (x,t), (x,t)e QT ,
u| 0 = ^(x), 0 < x < I, uxlx=0 = ux\x =, = 0 < t < T . Here i, T > 0 are given numbers f (x,t)e L2 (QT), ^(x)eW21 (0,i) are given functions, k(x,t), q(x,t) are unknown coefficients, u(x,t) = (k(x,t),q(x,t)) is a control, u = u(x,t) = u(x,t;u) is the solution to the boundary value problem corresponding to the control u = u(x,t). Let us introduce a set of admissible controls
V = {u(x,t) = (k(x,t),q(x,t)) e H = W,1 (QT) x L2 (QT): 0 < v < k(x,t) < |kx (x,t)| < |1,|kt (x,t)| <|2,0 < q0 < q(x,t) < q1 п.в.на QT}, where |> v> 0 , |1, |2 > 0 , q1 > q0 > 0 are given numbers.
Let us state the following coefficient inverse problem of optimal control type: among all the admissible controls u(x,t)=(k(x,t),q(x,t))eV, find the controls u,(x,i)=(k,(x,i),q,(i,t))eV minimizing the functional
t i
J (u) = U K (x, t )u (x, t; u) dx - E (t)
0 0
where K(x,t),E(t) are known functions, u = u(x,t) is a control u = u(x,t) = u(x,t;u) is a
2
dt
generalized solution to the boundary value problem from V21,0 (QT) corresponding to the control u = u(x, t) is a given set.
In the present work, the optimization formulation of the coefficient inverse problem for a parabolic equation with an additional integral condition is considered. The questions of correctness of the optimization formulation of the inverse problem are investigated. The differentiability of the objective functional is proved and the formula for its gradient is found. A necessary condition of optimality is found in the form of a variational inequality.
Keywords: optimal control, parabolic equation, integral boundary condition, optimality condition
TAGIEV Rafig Kalandar (Dr. Math. Sciences, prof. Baku State University, Azerbaijan) E-mail: [email protected]
KASUMOV Atakhan Rashid oglih (Senior Lecturer, Lankaran State University, Azerbaijan)
E-mail: [email protected]
REFERENCES
1. Tikhonov A.N. (1963) O reshenii nekorrektno postavlennykh zadach i metode regulyarizatsii [On the solution of ill-posed problems and regularization method]. Dokl. USSR Academy of Sciences. 151(3). pp. 501-504.
2. Glasko V.B. (1984) Obratnye zadachi matematicheskoy fiziki [Inverse problems of mathematical physics]. Moscow: MSU Publ.
3. Iskenderov A.D. (1984) O variatsionnykh postanovkakh mnogomernykh obratnykh zadach matematicheskoy fiziki [On the variational formulations of multidimensional inverse problems of mathematical physics]. Dokl. USSR Academy of Sciences. 274 (3). pp. 531-533.
4. Alifanov O.M., Artyuhin E.A., Rumyantsev S.V. (1988) Ekstremal'nye metody resheniya nekorrektnykh zadach [Extreme methods for solving ill-posed problems]. Moscow: Nauka.
5. Karchevsky A.L. (1997) Properties of the misfit functional for a nonlinear one-dimensional coefficient hyperbolic inverse problem. J. Inverse Ill-Posed Probl. 5(2). pp. 139-165. DOI: https://doi.org/10.1515/jiip.1997.5.2.139
6. Kabanikhin S.I., Iskakov K.T. (2001) Justification of the Steepest Descent Method for the Integral Statement of an Inverse Problem for a Hyperbolic Equation. Sib. Mat. Zh. 42(3). pp. 478-494. DOI: 10.1023/A:1010471125870.
7. Kabanikhin S.I. Dairbaeva G. (2007) Obratnaya zadacha nakhozhdeniya koeffitsienta uravneniya teploprovodnosti [The inverse problem of finding a coefficient of the heat equation]. Proc. International Conference «Inverse ill-posed problems of mathematical physics" devoted to the 75th anniversary of academician M.M. Lavrentev, Novosibirsk, Russia. pp. 1-5.
8. Ladyzhenskaya O.A., Solonnikov V.A., Uraltseva N.N. (1967) Lineynye i kvazilineynye uravneniya parabolicheskogo tipa [Linear and quasilinear equations of the parabolic type]. Moscow: Nauka.
9. Lions J.L. (1985) Control of Distributed Singular Systems. Paris: Gauthier-Villars.
10. Tagiyev R.K. (2009) Optimal control of coefficients in parabolic systems. Diff. equation. 45(10). pp. 1492-1501. DOI: 10.1134/S0012266109100164.
11. Tagiyev R.K. (2009) Optimal control of coefficients of a quasilinear parabolic equation. Automation and Remote Control. 70(11). DOI: 10.1134/S0005117909110058.
12. Tagiyev R.K. (2013) An optimal control problem for a quasilinear parabolic equation with controls in the coefficients and phase constraints. Diff. equation. 49(3). pp. 369-381. DOI: 10.1134/S0012266113030129.
13. Vasilyev F.P. (1981) Metody resheniya ekstremal'nykh zadach [Methods for solving extreme problems]. Moscow: Nauka.
14. Ladyzhenskaya O.A., Uraltseva N.N. (1973) Lineynye i kvazilineynye uravneniya ellipticheskogo tipa [Linear and quasilinear equations of the elliptic type]. Moscow: Nauka.