RUDN Journal of MlPh 2018 Vol. 26 No. 2 103-118
Вестник РУДН. Серия МИФ http://joumals.rudn.ru/miph
Математика
УДК 517.95
DOI: 10.22363/2312-9735-2018-26-2-103-118
Разрешимость линейной обратной задачи для эволюционного уравнения с суперустойчивой полугруппой
И. В. Тихонов*, Ву Нгуен Шон Тунг1"
* Кафедра математической физики МГУ имени М. В. Ломоносова Ленинские горы, ГСП-1, Москва, Россия, 119991 ^ Кафедра математического анализа Московский педагогический государственный университет ул. Краснопрудная, д. 14, Москва, Россия, 107140
Для эволюционного уравнения в банаховом пространстве изучается линейная обратная задача о нахождении «источника». Требуется восстановить неизвестное неоднородное слагаемое при помощи дополнительного нелокального условия, выраженного в виде интеграла Римана—Стильтьеса. Основное предположение связано с суперустойчивостью (квазинильпотентностью) эволюционной полугруппы. Точнее, предполагается, что эволюционная полугруппа, ассоциированная с абстрактным дифференциальным уравнением, имеет бесконечный отрицательный экспоненциальный тип. Без других ограничений установлена теорема об однозначной разрешимости обратной задачи. Показано, что решение представи-мо сходящимся рядом Неймана. Предъявлены точные условия, при которых бесконечный ряд обращается в конечную сумму. Здесь алгоритм вычисления решения становится финитным. Разобраны модельные примеры, в том числе — важный пример обратной задачи с финальным переопределением. Перечисленные результаты могут найти применение в специальных разделах математической физики, связанных с теорией упругости и задачами линейного переноса. Как принято, наше исследование проходит «в случае общего положения» — при выборе комплексного поля скаляров, но основные факты справедливы также и в вещественном случае. Созданная теория допускает перенос на нелокальные задачи для эволюционных уравнений, когда для нахождения решения вместо традиционного начального условия используют специальные усреднения по времени.
Ключевые слова: эволюционное уравнение, обратная задача, суперустойчивая полугруппа, операторное уравнение, ряд Неймана, теорема существования и единственности решения
1. Введение
В банаховом пространстве Е при фиксированном значении Т > 0 рассмотрим абстрактную задачу Коши для эволюционного уравнения специального вида:
(и'(1) = Аи(1) + 0 < 1 < Т,
| и(0) = и0.
Предполагаем, что А — линейный замкнутый оператор в Е с плотной областью определения В(А) с Е, причём А порождает в Е полугруппу и(Ь) класса С0 (подробнее см. [1-5]). Элемент и0 задан в В(А). Скалярная функция <р(1) ф 0 непрерывна на отрезке [0, Т]. Считаем также, что <¿>(1) имеет ограниченную вариацию на [0, Т], т.е. (р е ВУ[0,Т].
Статья поступила в редакцию 12 января 2018 г.
Тогда при любом выборе элемента д € Е стандартная формула
и(г) = и (г)ио + ! и (г - з)ф)д ¿8, о < г < т, (2)
0
даёт классическое решение задачи (1), такое, что и € С1 ([0,Т];#) и и(Ь) € В (А) при всех £ € [0,Т] (см. [6]). Формула (2) однозначно восстанавливает нужное решение: иные решения, отличные от (2), в задаче Коши (1) отсутствуют (см. [4, с. 105-106]).
Допустим теперь, что элемент д € Е неизвестен. Для его нахождения добавим к (1) специальное нелокальное условие (переопределение) вида
т
/ и(1) ё^(^) = и1. (3)
Элемент ui задан в D(A). Интеграл в условии (3) есть стандартный векторный интеграл Римана—Стильтьеса (см. [1, гл. III]). Скалярная функция ß £ BV[0,Т] предполагается непрерывной справа в любой точке t £ [0,Т) и отличной от тождественной константы на [0,Т]. В частности, всегда
ß(0) = ß(0+) = Дт+^) (4)
Отсутствие скачка у функции ß(t) при t = 0 не есть ограничение общности, поскольку такой (возможный) скачок можно легко убрать из (3) с помощью заданного начального условия и(0) = и0.
Итак, требуется подобрать элемент д £ Е, при котором векторная функция u(t) из формулы (2) удовлетворяет дополнительному условию (3). Поставленная задача (1), (3) для нахождения пары ( u(t), д) называется линейной обратной задачей (с нелокальным условием). Общая теория подобных обратных задач была разработана в [6,7] в рамках известного направления [8]. Наше настоящее исследование (кратко изложенное в заметке [9]) проходит при дополнительном предположении, что полугруппа U(t), порождённая оператором A, является квазинильпотентной или суперустойчивой. Этот вырожденный класс специальных полугрупп привлёк внимание специалистов в последние десятилетия [10-15].
2. Квазинильпотентные операторы и полугруппы
Напомним нужные сведения из спектральной теории [1,2]. Исходную норму в пространстве Е обозначим через || ■ ||, а единичный оператор — через I. Как известно, линейный ограниченный оператор B: Е ^ Е называют квазинильпотентным или вольтерровым, если его спектральный радиус равен нулю:
r(B) = lim d || Bfc || = 0. (5)
fc^+ro v
Соотношение (5) эквивалентно тому, что o(B) = {0}, т.е. спектр квазинильпотент-ного оператора B состоит лишь из нуля. При этом резольвентный ряд
го
(Л I - B)-1 = £ — B» (6)
п=0
сходится по операторной норме всюду на множестве Л е С \ {0}. Использование разложений типа (6) без всяких ограничений на значение Л = 0 составляет главное преимущество квазинильпотентных операторов.
При проверке квазинильпотентности непосредственное обращение к формуле (5) часто бывает затруднительным из-за сложной структуры изучаемого оператора В. Здесь может пригодиться такой критерий (ср. с [16, с. 14-16]).
Лемма 1. Пусть В: Е ^ Е — линейный ограниченный оператор. Следующие утверждения эквивалентны
1) оператор В квазинильпотентен;
2) для Уе > 0 существует эквивалентная норма || ■ ||(е) в Е, такая,, что
\\В||(Е) < е. (7)
Доказательство. Пусть В квазинильпотентен, т. е. г(В) = 0. При любом фиксированном выборе £ > 0 имеем
lim d \\e-kBkу = lim e~1 f/ \\Bk|| = e~1 r(B) = 0.
В частности, последовательность \\e kВk\\ будет заведомо ограниченной, и найдётся число М > 0, такое, что Це"kBk\\ ^ М при всех k G N U {0}. Положим теперь
\\/\\(Е) = sup це-fcBkf \\, f G Е. (8)
fceNujo}
Свойства нормы для (8) проверяются непосредственно. Поскольку
\\/\\ = IIе"kßkf IIlfc-o < \\/H(s) < sup Це"fcBk\\\\/\\ < M\\f \\, f G E,
fceNujo}
то норма (8) эквивалентна исходной норме \\ ■ \\. При этом
\\Bf \\(Е) = sup Це-fcBk+1f \\ = sup £ Це"(fc+1)Bk+1f \\ = fceNujo} fceNujo}
= £ sup \\£-mBmf \\ < £ sup Це" mBmf \\ = £ \\f \\(e), f G E,
rneN rneNUjo}
т.е. для нормы (8) выполнена требуемая оценка (7). Итак, из 1) следует 2).
Обратная импликация почти очевидна. Действительно, пусть выполнено 2). В силу оценки (7) спектр о(В) принадлежит множеству De = {Л G C : |А| ^ е}. Поскольку значение е > 0 является произвольным, то о(В) = {0}. Тем самым оператор В квазинильпотентен. В результате 1) и 2) эквивалентны. Лемма доказана. □
Более тонкие обстоятельства, связанные с квазинильпотентными операторами, для нас сейчас несущественны. Для полноты картины отметим два сравнительно недавних исследования по теме [17,18], а также ряд работ [19-21], где указаны специальные условия, при которых обобщённые интегральные операторы являются квазинильпотентными.
Перейдём к полугруппам. Напомним, что экспоненциальным типом операторной полугруппы U(t) класса Co называется величина
„o S lim УШ».
+ t
Указанный предел заведомо существует со значением в промежутке [-то, +то).
В случае, когда ш0 = —то, полугруппа U(t) называется квазинильпотентной или ещё суперустойчивой. Мы будем использовать в основном второй термин, говоря именно про «суперустойчивость» полугруппы по аналогии с английским "superstability", принятым в современной западной литературе (см. [10-15]). Данное название вполне оправдано, поскольку требование Шо = —то эквивалентно тому, что при любом выборе числа а > 0 найдётся константа М = Ма ^ 1, для которой
||UС0Ц < Ме~at, t^ 0. (9)
Оценка (9) допускает дополнительную коррекцию при подходящем изменении нормы.
Лемма 2. Пусть U (t) — суперустойчивая полугруппа класса С0. Тогда при любом выборе числа а > 0 найдётся эквивалентная норма || ■ ||а в пространстве Е, такая, что
||U(Ща "ai, 0. (10)
Доказательство. Воспользуемся известным приёмом [1, с. 376]. Зафиксируем число а > 0 и запишем оценку (9) с соответствующей константой М = Ма ^ 1. Положим
||/|U = sup ||e^U(r)f ||, feE. (11)
т ^0
Учитывая (9), получаем, что ||/|| < ||/||а ^ М||/|| для любого f e Е. Свойства нормы для || ■ ||а проверяются непосредственно. Итак, формула (11) действительно корректно определяет эквивалентную норму в Е. Кроме того, если t ^ 0, то
||U(t) f ||а = sup ||е ^U(т + t)f || = sup е||е^+t)U(т + t)f || =
т^0 т^0
= е sup || е asU (s)f || < е sup || easU (s)f || = е ||/||а, feE,
откуда и следует нужная оценка (10). Лемма доказана. □
Суперустойчивая полугруппа U(t) называется также квазинильпотентной, поскольку все её операторы U(t), взятые при t> 0, являются квазинильпотентными в пространстве Е. При этом спектр о(А) порождающего оператора А будет пустым, а резольвентное множество р(А) заполнит всю комплексную плоскость. В частности, заведомо существует ограниченный обратный оператор А~1: Е ^ Е, что несколько упрощает формулы при изучении поставленной задачи (1), (3) (см. [6]).
Особым примером суперустойчивой (квазинильпотентной) полугруппы U(t) служит нильпотентная полугруппа, когда
U(t) = 0, V > to > 0, (12)
с некоторым фиксированным значением to > 0. Нильпотентные полугруппы естественно возникают при рассмотрении процессов простого переноса в той или иной ограниченной области (см. [1, с. 552] или [5, с. 363-364]).
3. Формулировка основных результатов
Вернёмся к рассмотрению обратной задачи (1), (3). Из теории, развитой в [6], следует такое утверждение (см. [6, с. 112]).
Теорема 1. Пусть линейный замкнутый оператор А порождает суперустойчивую полугруппу и(£) класса Со. Пусть е С[0,Т] П ВУ[0,Т] и ц е ВУ[0,Т] с выполненным условием (4). Пусть также
т
¡3 = у ^)ф(*) = 0. (13)
о
Тогда обратная задача (1), (3) имеет и притом единственное решение (и(!),д) при любом выборе элементов ио,и1 е Е(А).
Метод доказательства теоремы 1, предложенный в [6], носил неконструктивный характер и базировался на теории Гельфанда коммутативных банаховых алгебр, точнее, на специальной версии этой теории, разработанной в трактате [1] для описания спектров интегралов от полугрупп. Дополнительный анализ показывает, что можно избежать обращения к таким сложным концепциям и получить теорему 1 элементарным конструктивным методом.
В предположениях теоремы 1 рассмотрим операторное уравнение на неизвестный элемент д е Е. Согласно схеме [6] уравнение допускает запись
!3д - Вд = f (14)
со значением /3 = 0 из формулы (13), линейным ограниченным оператором
т т 4
В = J <р(0)и ф фф + у фф у и (г - в)ё^(в) (15)
о оо
и заданным элементом f е Е вида
í = А ^ У и(г)ио ф(*) - . (16)
Здесь и в дальнейшем все операторные интегралы типа (15) понимаются в сильной операторной топологии — со сходимостью на произвольном элементе Н е Е. Используя структуру выражения (15) и экспоненциальную оценку (10) с подходящим выбором значения а > 0, можно оценить оператор В в эквивалентной норме (11) и показать его квазинильпотентность. Наша цель — установить прямым методом следующий принципиальный факт.
Теорема 2. Пусть выполнены предположения теоремы 1. Тогда оператор В из формулы (15) является квазинильпотентным в пространстве Е. При этом решение однозначно разрешимого операторного уравнения (14) со значением Р = 0 представимо рядом Неймана
£
1 Вк/, (17)
к=о м
сходящимся по норме пространства Е. Если элемент / е Е взят в виде (16), а комплексное число Р = 0 — в виде (13), то элемент д из формулы (17) даёт второй компонент решения (и(£),д) поставленной обратной задачи (1), (3). При этом первый компонент решения — функция и(£) — выражается формулой (2).
Теорема 2 есть главный результат нашей работы. Она обосновывает итерационный алгоритм для поиска решения обратной задачи (1), (3) при помощи ряда Неймана (17). Алгоритм приобретает законченный вид в специальном случае ниль-потентной полугруппы и (£). Тогда, при некоторых дополнительных ограничениях на функции можно доказать такое утверждение.
Теорема 3. Пусть выполнены предположения теоремы 1. Пусть, дополнительно, полугруппа и(£) является нильпотентной, удовлетворяя соотношению (12) с фиксированным .значением ¿° > 0. Предположим также, что функция является постоянной на некотором промежутке [0, а] С [0, Т), а функция является постоянной на некотором промежутке [Т — Ь,Т] С (0, Т], причём
а + Ь >Т. (18)
Тогда оператор В из формулы (15) является нильпотентным и
Вк = 0, к е М, к > -. (19)
а+Ь—Т ' '
Соответственно, решение операторного уравнения (14) представимо в виде конечной суммы
N0 1 г ,
рк+1 к=0 м
замещающей бесконечный ряд Неймана (17).
а + Ь — Т
— 1, (20)
Здесь через [с] обозначен потолок числа с е К, т. е. наименьшее целое число, большее или равное с (см. [22, с. 88]). Теоремы 2 и 3 различны по характеру, и доказывать их удобно раздельно.
4. Доказательство основной теоремы 2
По поводу теории интеграла Стильтьеса см. [1,23]. Без всяких пояснений используем стандартные оценки типа
< тех ||ь(Щ ■ Уаг{/^)}
верные при любом выборе V е С([0,Т],Е^) и ц е ВУ[0, Т]. Особо отметим следующий полезный факт.
Лемма 3. Пусть у е С[0,Т] П ВУ [0,Т]. Тогда для любого заданного числа 5 > 0 можно подобрать такое значение т = т((р, 5) > 0, что при всех £ 1, ¿2 е [0,Т], удовлетворяющих условию 0 < ¿2 —t 1 < т, выполняется соотношение
Уаг{у( 8 )}
Доказательство. Определим функцию
№ = УахМ«)}
¿2 5 < -.
^ 2
0 < г < т,
0
монотонно возрастающую на [0,Т]. Поскольку по условию е С[0,Т] П ВУ[0,Т], то также и ф е С[0,Т] П ВУ[0,Т] (см. [23, с. 210]). По теореме Кантора функция 'ф(р) будет равномерно непрерывной на [0,Т], т.е. для любого 6 > 0 найдётся такое значение т > 0, что при всех ¿1, ¿2 е [0, Т], удовлетворяющих условию 0 < ¿2 — ^ ^ т, выполнено соотношение
- ^1)| = Ф(^) - Ф(Ь) = УахМ«)}
<2 5 < -.
^ 2
Но это и утверждается в (21). Лемма доказана. □
Оценим теперь отдельные слагаемые в основном операторе В из формулы (15). Все предположения теоремы 1, а, значит, и теоремы 2, относительно полугруппы и(¿) и функций считаем выполненными. Напомним, что операторные
интегралы рассматриваются (и оцениваются) в сильной операторной топологии — на произвольном элементе Н е Е, который для краткости исключаем из записей.
Лемма 4. Пусть оператор В1 определён формулой
т
Вх = у и (¿)ф(£). (22)
о
Тогда для любого 5 > 0 найдётся такое а1 > 0, что при всех а > а1 верна оценка,
11^1 ||а <6. (23)
Здесь || ■ ||а — эквивалентная норма, заданная в пространстве Е по формуле (11).
Доказательство. Зафиксируем произвольное 5 > 0. В силу предположения (4) найдётся такое 7 е (0, Т], что
7 6
Уаг{^)} о < 2. (24)
При этом, если 7 = Т, то оценка (23) выполняется при любом значении а > 0, так как, согласно (10), имеем
||В1||а < тах е~а1 ■ Уаг{^)} = Уаг{^)}
о<4<Т
т
о
7 6 Х о2
Итак, при 7 = Т значение «1 > 0 можно выбирать произвольно. Если же 7 е (0, Т), запишем оператор В1 в виде
7 Т 7 Т
В1 = у и (г) фф +1 и (г) фф = у и (г) фф + и (7) У и (г - 7) фф. (25)
о 7 о 7
Оценивая (25) с учётом (10) и (24), получаем
7
||В1||а < тах е~аг ■ Уаг{^)} ' + е"а7 ■ тах е"а{1~7) ■ Уаг{^)}
т
7
о
т § т
7 2 .....о
7
= Уаг{^)} + е~а7 ■ Уаг{^)} <- + е~а7 ■ Уаг{^)} . (26)
Выберем теперь а\ > 0 так, чтобы при всех а > а\ действовало соотношение
т §
е~а1 • Уагмт т < 2. (27)
Комбинируя (26) и (27), получаем нужную оценку (23). Лемма доказана. □
Для оценки второго слагаемого в операторе (15) требуется следующий предварительный результат.
Лемма 5. Пусть операторная функция Р(¿) определена формулой
г
р ф = I и (г — з) ¿<р(з), о < г < т. (28)
0
Тогда для любого 5 > 0 найдётся такое а2 > 0, что при всех а > а2 верна оценка,
||Р(¿)||а <6, о (29)
Здесь || • ||а — эквивалентная норма, заданная в пространстве Е по формуле (11).
Доказательство. Зафиксируем произвольное 5 > 0. Выберем значение т > 0 как в лемме 3, чтобы соотношение (21) выполнялось при всех £ 1, ¿2 € [0, Т], удовлетворяющих условию 0 < ¿2 — £ 1 < т. При а > 0 и 0 т, используя (10) и (21), получаем
IIP(i)IU < max е~a(t~s) ■ Var{(( S)} = Var{(( s)}
0<s<i
0
< о < 8. о 2
Тем самым оценка (29) выполнена при любом а > 0 и 0 т. Возможно,
конечно, что т = Т, и тогда все доказано с произвольным выбором значения а2 > 0. При 0 < т <Т иг <1 ^ Т перепишем Р (¿) в виде
t-T
P (i)=y и (t - s)d(p(s)+ J и (t - s)d((s), T<t ^T. (30)
0 t-T
Пусть снова a > 0. Оценивая (30) с учётом (10) и (21), получаем
t — T
IIP(i)IU < max e~a(t~s) ■ Var{(( s)} + max e~a(t~s) ■ Var{(( s)}
O^s^i—t 0 t—T ^s^t
T
0
e~aT ■ Var{(( s)} +Var{(( s)} "aT ■ Var{(( s)} + -. (31)
т s
02
Выберем теперь а2 > 0 так, чтобы при всех а > а2 действовало соотношение
т
е~ат • Var{^( ^ )} о <2. (32)
Комбинируя (31) и (32), получаем оценку (29) при т <1 ^ Т .В итоге с указанным выбором т > 0 и а2 > 0 нужный результат обоснован на всем отрезке [0, Т]. Лемма доказана. □
Отметим, что, если ф £ С1 [0, Г], то доказательство леммы 5 сильно упрощается. Действительно, обозначив К = max |w'(s)|, для функции (28) имеем оценку
I|F (i)IU ^ f е~а( 1 ~s) ( s )l d S < К j e ~a(t ~s) ds <К <s, 0 ^t^T,
OO
верную при a > a2 = К/5 с произвольным фиксированным ё > 0.
Приступим к доказательству самой теоремы 2. Оценим основной оператор В из формулы (15). Заметим, что В = ф(0)В1 + В2, где оператор В1 определён формулой (22), а оператор В2 имеет вид
T T
В2 = J dn(t) J U(t - s)d(p(s) = J F(t)dn(t) (33)
O O O
с функцией F(t) из формулы (28).
Возьмём какое-то значение ё > 0. Принимая во внимание предыдущие леммы 4 и 5 с их выбором a1 и a2, положим a0 = max(a1,a2). Тогда при всех a > a0 в эквивалентной норме || ■ ||a, заданной по формуле (11), получим, что
I ВЦ, < |w(0)| -S + S- Var{^.(t)}|T = (|w(0)| + Var{/Y.(t)|
о
Зафиксируем произвольное е > 0. Ясно, что при достаточно малых ё > 0 имеем оценку ||-В||а ^ £, верную при всех а > ао с соответствующим значением ао > 0. По лемме 1 это гарантирует квазинильпотентность оператора В. Но тогда операторное уравнение (14) при любом выборе числа /3 = 0 будет однозначно разрешимым, и его решение д е Е представимо рядом Неймана (17).
То, что уравнение (14) с оператором (15) и заданным элементом (16) эквивалентно исходной обратной задаче (1), (3), следует из прежних результатов [6]. Теорема 2 полностью доказана.
5. Доказательство теоремы 3
Пусть функции ф((р), удовлетворяют условиям теоремы 3, включая специальное требование (18). Тогда для операторов В1, В2 из формул (22), (33) возможны представления
т т
В1 = J и (г) ф(£) = и (а + Ь - Т )и (Т - Ь) I и (г - а) фф,
а а
т т-ь т т-ь
В2 = ! ^ц(Ъ) ! и (Ь - 8)й<р(8) = и (а + Ь -Т) ! ¿¡¿(г) ! и (Ь -а + Т -Ь- 8)й<р(8).
Следовательно, основной оператор В из формулы (15) выражается в виде
В = ф(0)В1 + В2 = U(a + b - T)Q, (34)
где
Т Т Т-Ь
я = <р(0)и (Т - Ь) I и (г - а)<^) + I ! и (г - а + Т - Ь - з)<(р(з). (35)
а а 0
Оператор и (а + Ь - Т) коммутирует с оператором Q. По предположению теоремы 3 полугруппа и(^ является нильпотентной в пространстве Е и удовлетворяет соотношению (12) с фиксированным значением ¿о > 0. Поэтому
Вк = и (к (а + Ь - Т))Як = 0, к е М, к(а + Ь - Т) > ¿о,
что и даёт формулу (19). Итак, нильпотентность оператора В установлена. В силу (19) ряд Неймана (17) переходит в конечную сумму (20). Теорема доказана.
Как видим, ключевую роль в доказательстве теоремы 3 сыграло специальное разложение (34) с оператором Q из формулы (35). Корректность конструкции обусловлена требованием (18). Отметим примеры, когда такая конструкция сильно упрощается.
6. Модельные примеры
Всюду в данном пункте считаем, что полугруппа и(^ является нильпотентной в пространстве Е, удовлетворяя условию (12) с фиксированным значением ¿о > 0.
Пусть ^>(1) = 1 на [0,Т]. Этот особый случай часто возникает на практике. Поскольку функция является постоянной на любом промежутке [Т - Ь,Т] С (0,Т], то для применения теоремы 3 достаточно требовать, чтобы функция была постоянной на промежутке [0, а] С [0,Т) с неким (хоть каким-то) значением а е (0,Т). В ответ можно всегда подобрать такое Ь е (0,Т), что требование (18) окажется выполненным. Впрочем, здесь нет смысла прибегать к теореме 3, поскольку прямое исследование гораздо нагляднее.
Действительно, зафиксируем значение а е (0,Т) и рассмотрим обратную задачу
(и'(г) = Аи(г) + д, 0 < г < т,
(36)
и(0) = и0, ! и(Ь)<^(Ь) = щ.
Стандартное требование (13) переходит в условие
Р = ц(Т) - ц(а) = 0. (37)
Оператор (15) приобретает вид
т т
В = ! и (г) фф = и (а) ! и (г - а) фф. (38)
а а
Поскольку и(ка) = 0 при к > ¿0/а, то Вк = 0 при любом натуральном к > 10/а, т.е. оператор (38) гарантированно является нильпотентным. Формула (17) для неизвестного элемента д е Е сводится к конечной сумме
N 1
£ 7^вк * ^ -
3к+ к=0 м
а
- 1, (39)
с константой 3 = 0 из формулы (37), оператором В из формулы (38), и элементом
/ = А (а) j и(Ь - а)ио ф(^) - щ| . (40)
В случае, когда а ^ ¿0 (а такое вполне возможно), формулы (38)—(40) радикально упрощаются, и элемент д выражается одним слагаемым
9 =--1 / ч Ащ. (41)
МТ) - Ма)
Любопытно, что такой ответ (41) зависит не от выбора той или иной конкретной функции а лишь от значения 3 = ^(Т) — ^(а) = 0.
Подведём итог: именно соотношения (37)—(41) и надо использовать для нахождения решения обратной задачи (36) с оператором А, порождающим нильпотентную полугруппу и( ).
Специальный интерес представляет частный случай обратной задачи (36), когда есть функция скачков, локализованных в промежутке [а, Т] С (0, Т]. Общее нелокальное условие (3) заменяется многоточечным условием, и обратная задача принимает вид
'и'(г) = Аи(г) + д, 0
и(0) = ио, ^аки(тк ) = щ. (42)
к
Здесь {Тк } — конечный или счётный набор различных точек из [ а, Т], занумерованных натуральным индексом к. Все числовые коэффициенты ак считаются отличными от нуля и такими, что ^ |ак| < то. Требование (13) переходит в условие
3 = ^ак = 0. (43)
к
Элемент д £ Е в обратной задаче (42) выражается прежней конечной суммой (39), где
В = ^аки(тк), / = А( ^аки(Тк)ио -щ), (44)
к к
а значение 3 вычисляется через (43). Окончательное упрощение формул (43), (44) происходит в случае обратной задачи с финальным переопределением.
Действительно, пусть = 1 на [0,Т], а есть функция единичного скачка в финальной точке £ = Т. Тогда исходная обратная задача (1), (3) принимает вид
IV (¿) = Аи(£) + д, 0
^ (45)
I и(0) = и0, и(Т) = щ.
Это и есть модельная обратная задача с финальным переопределением (см. [8]). Ввиду особой важности примера (45) разберём его более подробно. Следуя общим формулам (13), (15), (16), получаем
р = 1, В = и(Т), / = А(и(Т)ио -щ). (46)
Поскольку Вк = и(кТ) = 0 при любом натуральном к ^ Ьо /Т, то общая разрешающая формула (17) вновь сводится к конечной сумме, имеющей в данном случае особенно простой вид
N0 .
д = £ и (кТ)/, N0 = Т
к=0
- 1, (47)
с элементом f из (46). Возможно, конечно, что Т ^ ¿о, и тогда N0 = \tо/Т\ — 1 = 0. При этом также и(Т) = 0 и / = — Аиг. В результате общий шаблон (47) даёт ответ
д = —Аиг.
С другой стороны, возможно, что Т < ^ • Тогда, принимая во внимание представление (46) для элемента f и проводя элементарные преобразования, имеем
N0
д = £ и (кТ)/ = / + и (Т)/ + и (2Т)/ + ... + и (ЩТ)/ =
к=0
= А (и (Т)ио + и (2Т)ио + и (3Т)ио + ... + и ((N0 + 1)Т) ио) — — А (иг + и (Т )и\ + и (2Т)щ + ... + и ЩТ) иг) =
N0
= —А иг + А ^ и (кТ) (ио — иг), N0 =
к=1
Т
1.
Слагаемое и ((Щ0 + 1) Т) и0 пропадает, ибо (Щ0 + 1)Т > 10 и и ((Щ0 + 1)Т) = 0.
В итоге формулу (47) для второго компонента д решения обратной задачи (45) можно записать в виде
{—Аиг, 0 < ^ < Т,
Л (48)
—Аиг +А^и (кТ)(и0 — иг), ^ >Т> 0,
к=г
где Щ0 = \ ^/Т\ — 1. Выражение (48) удобно на практике. Для последующих применений данного результата зафиксируем точное утверждение.
Теорема 4- Пусть линейный замкнутый оператор А порождает нильпотент-ную полугруппу и(1) класса С0, удовлетворяющую условию (12) с фиксированным значением, 10 > 0. Тогда при любом выборе и0,иг е В (А) элемент д е Е, разрешающий обратную задачу (45), находится по формуле (48).
При этом функция и(1) восстанавливается по формуле
г
и(Ъ) = и(Ъ)и0 +[и(з)д<з, 0 < Ь < Т, (49)
непосредственно следующей из (2) в случае, когда <p(t) = 1. Формула (49) используется также в обратных задачах (36) и (42).
7. Заключение
Итак, теория линейной обратной задачи (1), (3) с оператором А, порождающим суперустойчивую полугруппу, оказывается в основном завершённой. Полученные результаты представляют интерес для некоторых специальных разделов математической физики, связанных с теорией упругости и задачами линейного переноса. Например, на основе формул, установленных нами в п. 6, можно разработать конкретные алгоритмы решения обратных задач типа (36) для многомерного уравнения переноса, взятого в среде без столкновений, но с возможным поглощением частиц. Этот материал полезно опубликовать отдельно.
Литература
1. Хилле Э., Филлипс Р. Функциональный анализ и полугруппы. — М.: ИЛ, 1962.
2. Данфорд Н, Шварц Д. Линейные операторы. Т. 1. Общая теория. — М.: ИЛ, 1962.
3. Крейн С. Г. Линейные дифференциальные уравнения в банаховом пространстве. — М.: Наука, 1967.
4. Pazy A. Semigroups of Linear Operators and Applications to Partial Differential Equations. — N.Y.: Springer Verlag, 1983.
5. Engel K.-J., Nagel R. One-Parameter Semigroups for Linear Evolution Equations. — N.Y.: Springer, 2000.
6. Тихонов И. В., Эйдельман Ю. С. Вопросы корректности прямых и обратных задач для эволюционного уравнения специального вида // Математические заметки. — 1994. — Т. 56, вып. 2. — С. 99-113.
7. Прилепко А. И., Тихонов И. В. Восстановление неоднородного слагаемого в абстрактном эволюционном уравнении // Известия РАН. Сер. матем. — 1994. — Т. 58, вып. 2. — С. 167-188.
8. Prilepko A. I., Orlovsky D. G., Vasin I. A. Methods for Solving Inverse Problems in Mathematical Physics. — N.Y.: Basel, 2000.
9. Тихонов И. В., Тунг В. Н. Ш. Метод решения обратной задачи для эволюционного уравнения с суперустойчивой полугруппой // Дифференциальные уравнения и процессы управления. — 2017. — Вып. 2. — С. 51-58.
10. Balakrishnan A. V. On Superstability of Semigroups // Systems modelling and optimization. Proceedings of the 18th IFIP Conference on System Modelling and Optimization. CRC Research Notes in Mathematics / Ed. by M. P. Polis et al. — 1999. — Pp. 12-19.
11. Balakrishnan A. V. Smart Structures and Super Stability // Evolution Equations and Their Applications in Physical and Life Sciences. Lecture Notes in Pure and Applied Mathematics / Ed. by G. Lumer, L. Weis. — Vol. 215. — 2001. — Pp. 43-53.
12. Balakrishnan A. V. Superstability of Systems // Applied Mathematics and Computation. — 2005. — Vol. 164, issue 2. — Pp. 321-326.
13. Creutz D., Mazo M., Preda C. Superstability and Finite Time Extinction for Co-Semigroups // arXiv:0907.4812. Submitted. — 2013. — Pp. 1-12.
14. Chen J.-H., Lu W.-Y. Perturbation of Nilpotent Semigroups and Application to Heat Exchanger Equations // Applied Mathematics Letters. — 2011. — Vol. 24. — Pp. 1698-1701.
15. Kmit I., Lyul'ko N. Perturbations of Superstable Linear Hyperbolic Systems // arXiv:1605.04703. Submitted. — 2017. — Pp. 1-26.
16. Приближенное решение операторных уравнений / М. А. Красносельский, Г. М. Вайникко, П. П. Забрейко и др. — М.: Наука, 1969.
17. Malinen J., Nevanlinna O., Zemanek J. Microspectral Analysis of Quasinilpotent Operators // arXiv:1211.4790v1. Submitted. — 2012. — Pp. 1-25.
18. Eskandari R., Mirzapour F. Hyperinvariant Subspaces and Quasinilpotent Operators // Bulletin of the Iranian Mathematical Society. — 2015. — Vol. 41, issue 4. — Pp. 805813.
19. Забрейко П. П. О спектральном радиусе интегральных операторов Вольтерра // Литовский математический сборник. — 1967. — Т. 7, вып. 2. — С. 281-287.
20. Жуковский Е. С. К теории уравнений Вольтерра // Дифференциальные уравнения. — 1989. — Т. 25, вып. 9. — С. 1599-1605.
21. Сумин В. И., Чернов А. В. Операторы в пространствах измеримых функций: вольтерровость и квазинильпотентность // Дифференциальные уравнения. — 1998. — Т. 34, вып. 10. — С. 1402-1411.
22. Грехэм Р., Кнут Д., Паташник О. Конкретная математика. — М.: Мир, 1998.
23. Натансон И. П. Теория функций вещественной переменной. — М.: Наука, 1974.
UDC 517.95
DOI: 10.22363/2312-9735-2018-26-2-103-118
The Solvability of the Inverse Problem for the Evolution Equation with a Superstable Semigroup
I. V. Tikhonov*, Vu Nguyen Son Tung^
* Faculty of Computational Mathematics and Cybernetics Lomonosov Moscow State University Leninskyie Gori, Moscow, 119991, Russian Federation ^ Mathematical Analysis Department Moscow State University of Education 14 Krasnoprudnaya Str., Moscow, 107140, Russian Federation
For the evolution equation in a Banach space, the linear inverse source problem is studied. It is required to recover an unknown nonhomogeneous term by means of an additional nonlocal condition written in the form of a Riemann-Stieltjes integral. The main assumption is related to the superstability (quasinilpotency) of the evolution semigroup. More precisely, it is assumed that the evolutionary semigroup associated with the abstract differential equation has an infinite negative exponential type. Without other restrictions, a theorem on the solvability of the inverse problem is obtained. It is shown that the solution can be represented by a convergent Neumann series. Exact conditions under which an infinite series becomes a finite sum are found. Here, the algorithm for calculating the solution becomes finite. Model examples are considered, including an important example of the inverse problem with final overdetermination. The above results can be applicated in special parts of mathematical physics related to the theory of elasticity and the linear transport theory. As is customary, our research takes place in the general case with a choice of the complex scalar field, but the main facts are also true in the real case. The created theory allows transfer to nonlocal problems for evolution equations, when instead of the traditional initial condition special time averaging is used to find the solution.
Key words and phrases: evolution equation, inverse problem, superstable semigroup, operator equation, Neumann series, existence and uniqueness theorem of the solution
References
1. E. Hille, R. Phillips, Functional Analysis and Semigroups, IL, Moscow, 1962, in Russian.
2. N. Dunford, J. Schwartz, Linear Operators. P. 1. General Theory, IL, Moscow, 1962, in Russian.
3. S. G. Krein, Linear Differential Equations in Banach Space, Nauka, Moscow, 1967, in Russian.
4. A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, Springer Verlag, N.Y., 1983.
5. K.-J. Engel, R. Nagel, One-Parameter Semigroups for Linear Evolution Equations, Springer, N.Y., 2000.
6. I. V. Tikhonov, Yu. S. Eidelman, Problems on Correctness of Ordinary and Inverse Problems for Evolutionary Equations of a Special Form, Math. Notes 56 (1994) 830-839, in Russian.
7. A. I. Prilepko, I. V. Tikhonov, Recovery of the Nonhomogeneous Term in an Abstract Evolution Equation, Russian Acad. Sci. Izv. Math. 58 (1994) 167-188, in Russian.
8. A. I. Prilepko, D. G. Orlovsky, I. A. Vasin, Methods for Solving Inverse Problems in Mathematical Physics, Basel, N.Y., 2000.
9. I. V. Tikhonov, V. N. S. Tung, A Method of Solving the Inverse Problem for the Evolution Equation with a Superstable Semigroup, Journal Differential Equations and Control Processes (2) (2017) 51-58, in Russian.
10. A. V. Balakrishnan, On Superstability of Semigroups, in: M. P. Polis, et al. (Eds.), Systems Modelling and Optimization, Proceedings of the 18th IFIP Conference on System Modelling and Optimization, CRC Research Notes in Mathematics, Chapman and Hall, 1999, pp. 12-19.
11. A. V. Balakrishnan, Smart Structures and Super Stability, in: G. Lumer, L. Weis (Eds.), Evolution Equations and Their Applications in Physical and Life Sciences. Lecture Notes in Pure and Applied Mathematics, Vol. 215, 2001, pp. 43-53.
12. A. V. Balakrishnan, Superstability of Systems, Applied Mathematics and Computation 164 (2005) 321-326.
13. D. Creutz, M. Mazo, C. Preda, Superstability and Finite Time Extinction for Co-Semigroups, arXiv:0907.4812. Submitted (2013) 1-12.
14. J.-H. Chen, W.-Y. Lu, Perturbation of Nilpotent Semigroups and Application to Heat Exchanger Equations, Applied Mathematics Letters 24 (2011) 1698-1701.
15. I. Kmit, N. Lyul'ko, Perturbations of Superstable Linear Hyperbolic Systems, arXiv:1605.04703. Submitted (2017) 1-26.
16. M. A. Krasnoselsky, G. M. Vainikko, P. P. Zabreiko, Y. B. Rutitsky, V. Y. Stetsenko, Approximate Solution of Operator Equations, Nauka, Moscow, 1969, in Russian.
17. J. Malinen, O. Nevanlinna, J. Zemanek, Microspectral analysis of quasinilpotent operators, arXiv:1211.4790v1. Submitted (2012) 1-25.
18. R. Eskandari, F. Mirzapour, Hyperinvariant Subspaces and Quasinilpotent Operators, Bulletin of the Iranian Mathematical Society 41 (2015) 805-813.
19. P. P. Zabreiko, On the Spectral Radius of Volterra Integral Operators, Lithuanian Math. Collection 7 (1967) 281-287, in Russian.
20. E. S. Zhukovskii, On the Theory of Volterra Equations, Differential Equations 25 (1989) 1132-1137, in Russian.
21. V. I. Sumin, A. V. Chernov, Operators in the Spaces of Measurable Functions: the Volterra Property and Quasinilpotency, Differential Equations 34 (1998) 1403-1411, in Russian.
22. R. L. Graham, D. E. Knuth, O. Patashnik, Concrete Mathematics: A Foundation for Computer Science, Mir, Moscow, 1998, in Russian.
23. I. P. Natanson, Theory of Functions of a Real Variable, Nauka, Moscow, 1974, in Russian.
© Тихонов И. В., Ву Нгуен Шон Тунг, 2018
Для цитирования:
Тихонов И. В., Ву Нгуен Шон Тунг Разрешимость линейной обратной задачи для эволюционного уравнения с суперустойчивой полугруппой // Вестник Российского университета дружбы народов. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 2. — С. 103-118. — БОТ: 10.22363/2312-9735-2018-26-2-103-118.
For citation:
Tikhonov I.V., Vu Nguyen Son Tung The Solvability of the Inverse Problem for the Evolution Equation with a Superstable Semigroup, RUDN Journal of Mathematics, Information Sciences and Physics 26 (2) (2018) 103-118. DOI: 10.22363/2312-9735-201826-2-103-118. In Russian.
Сведения об авторах:
Тихонов Иван Владимирович — доктор физико-математических наук, профессор кафедры математической физики факультета ВМК МГУ имени М. В. Ломоносова (e-mail: [email protected], тел.: +7 9162963261)
Ву Нгуен Шон Тунг — аспирант кафедры математического анализа Московского педагогического государственного университета (e-mail: [email protected], тел.: +7 9169735556)
Information about the authors:
Tikhonov I. V. — Doctor of Physical and Mathematical Sciences, professor of Department of Mathematical Physics, Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University (e-mail: [email protected], phone: +7 9162963261) Vu Nguyen Son Tung — Postgraduate Student of Mathematical Analysis Department, Moscow State University of Education (e-mail: [email protected], phone: +7 9169735556)