Разрешимость линейной обратной задачи для эволюционного уравнения с суперустойчивой полугруппой
- Авторы: Тихонов И.В.1, Ву Нгуен Ш.Т.2
-
Учреждения:
- МГУ имени М. В. Ломоносова
- Московский педагогический государственный университет
- Выпуск: Том 26, № 2 (2018)
- Страницы: 103-118
- Раздел: Математика
- URL: https://journals.rudn.ru/miph/article/view/18365
- DOI: https://doi.org/10.22363/2312-9735-2018-26-2-103-118
Цитировать
Полный текст
Аннотация
Для эволюционного уравнения в банаховом пространстве изучается линейная обратная задача о нахождении «источника». Требуется восстановить неизвестное неоднородное слагаемое при помощи дополнительного нелокального условия, выраженного в виде интеграла Римана-Стильтьеса. Основное предположение связано с суперустойчивостью (квазинильпотентностью) эволюционной полугруппы. Точнее, предполагается, что эволюционная полугруппа, ассоциированная с абстрактным дифференциальным уравнением, имеет бесконечный отрицательный экспоненциальный тип. Без других ограничений установлена теорема об однозначной разрешимости обратной задачи. Показано, что решение представимо сходящимся рядом Неймана. Предъявлены точные условия, при которых бесконечный ряд обращается в конечную сумму. Здесь алгоритм вычисления решения становится финитным. Разобраны модельные примеры, в том числе - важный пример обратной задачи с финальным переопределением. Перечисленные результаты могут найти применение в специальных разделах математической физики, связанных с теорией упругости и задачами линейного переноса. Как принято, наше исследование проходит «в случае общего положения» - при выборе комплексного поля скаляров, но основные факты справедливы также и в вещественном случае. Созданная теория допускает перенос на нелокальные задачи для эволюционных уравнений, когда для нахождения решения вместо традиционного начального условия используют специальные усреднения по времени.
Полный текст
1. Введение В банаховом пространстве E при фиксированном значении T > 0 рассмотрим абстрактную задачу Коши для эволюционного уравнения специального вида: uў(t) = Au(t) + f(t)g,0 Ј t Ј T, u(0) = u0. (1) Предполагаем, что A - линейный замкнутый оператор в E с плотной областью определения D(A) М E, причём A порождает в E полугруппу U(t) класса C0 (подробнее см. [1-5]). Элемент u0 задан в D(A). Скалярная функция f(t)в‰ў0 непрерывна на отрезке [0,T]. Считаем также, что f(t) имеет ограниченную вариацию на [0,T], т. е. f О BV [0,T]. Тогда при любом выборе элемента g О E стандартная формула u(t) = U(t)u0 + т 0tU(t - s)f(s)gds,0 Ј t Ј T, (2) даёт классическое решение задачи (??), такое, что u О C1([0,T];E) и u(t) О D(A) при всех t О [0,T] (см. [6]). Формула (??) однозначно восстанавливает нужное решение: иные решения, отличные от (??), в задаче Коши (??) отсутствуют (см. [4, с. 105-106]). Допустим теперь, что элемент g О E неизвестен. Для его нахождения добавим к (??) специальное нелокальное условие (переопределение) вида т 0Tu(t)dm(t) = u 1. (3) Элемент u1 задан в D(A). Интеграл в условии (??) есть стандартный векторный интеграл Римана-Стильтьеса (см. [1, гл. III]). Скалярная функция m О BV [0,T] предполагается непрерывной справа в любой точке t О [0,T) и отличной от тождественной константы на [0,T]. В частности, всегда m(0) = m(0+) є limt®0+m(t). (4) Отсутствие скачка у функции m(t) при t = 0 не есть ограничение общности, поскольку такой (возможный) скачок можно легко убрать из (??) с помощью заданного начального условия u(0) = u0. Итак, требуется подобрать элемент g О E, при котором векторная функция u(t) из формулы (??) удовлетворяет дополнительному условию (??). Поставленная задача (??), (??) для нахождения пары (u(t),g) называется линейной обратной задачей (с нелокальным условием). Общая теория подобных обратных задач была разработана в [6, 7] в рамках известного направления [8]. Наше настоящее исследование (кратко изложенное в заметке [9]) проходит при дополнительном предположении, что полугруппа U(t), порождённая оператором A, является квазинильпотентной или суперустойчивой. Этот вырожденный класс специальных полугрупп привлёк внимание специалистов в последние десятилетия [10-15]. 2. Квазинильпотентные операторы и полугруппы Напомним нужные сведения из спектральной теории [1, 2]. Исходную норму в пространстве E обозначим через ∥Ч∥, а единичный оператор - через I. Как известно, линейный ограниченный оператор B: E ® E называют квазинильпотентным или вольтерровым, если его спектральный радиус равен нулю: r(B) є limk®+Ґв€ҐBk∥k = 0. (5) Соотношение (??) эквивалентно тому, что s(B) = {0}, т. е. спектр квазинильпотентного оператора B состоит лишь из нуля. При этом резольвентный ряд (lI - B)-1 = е n=0Ґ 1 ln+1Bn (6) сходится по операторной норме всюду на множестве l О в„‚ \{0}. Использование разложений типа (??) без всяких ограничений на значение l№0 составляет главное преимущество квазинильпотентных операторов. При проверке квазинильпотентности непосредственное обращение к формуле (??) часто бывает затруднительным из-за сложной структуры изучаемого оператора B. Здесь может пригодиться такой критерий (ср. с [16, с. 14-16]). Лемма 1. Пусть B: E ® E - линейный ограниченный оператор. Следующие утверждения эквивалентны 1) оператор B квазинильпотентен; 2) для "ε > 0 существует эквивалентная норма ∥Ч∥(ε) в E, такая, что ∥B∥(ε) Ј ε. (7) Доказательство. Пусть B квазинильпотентен, т. е. r(B) = 0. При любом фиксированном выборе ε > 0 имеем limk®+Ґв€Ґε-kBk∥k = lim k®+Ґε-1∥Bk∥k = ε-1r(B) = 0. В частности, последовательность ∥ε-kBk∥ будет заведомо ограниченной, и найдётся число M > 0, такое, что ∥ε-kBk∥Ј M при всех k О в„• И{0}. Положим теперь ∥f∥(ε) є supkОв„•И{0}∥ε-kBkf∥,f О E. (8) Свойства нормы для (??) проверяются непосредственно. Поскольку ∥f∥ = ε-kBkf k=0 Јв€Ґf∥(ε) Ј supkОв„•И{0}∥ε-kBk∥∥f∥Ј M∥f∥,f О E, то норма (??) эквивалентна исходной норме ∥Ч∥. При этом ∥Bf∥(ε) = supkОв„•И{0}∥ε-kBk+1f∥ = sup kОв„•И{0}ε∥ε-(k+1)Bk+1f∥ = = εsupmОℕ∥ε-mBmf∥Ј εsup mОв„•И{0}∥ε-mBmf∥ = ε∥f∥ (ε),f О E, т. е. для нормы (??) выполнена требуемая оценка (??). Итак, из 1) следует 2). Обратная импликация почти очевидна. Действительно, пусть выполнено 2). В силу оценки (??) спектр s(B) принадлежит множеству Dε = {l О в„‚ : |l|Ј ε}. Поскольку значение ε > 0 является произвольным, то s(B) = {0}. Тем самым оператор B квазинильпотентен. В результате 1) и 2) эквивалентны. Лемма доказана. Более тонкие обстоятельства, связанные с квазинильпотентными операторами, для нас сейчас несущественны. Для полноты картины отметим два сравнительно недавних исследования по теме [17, 18], а также ряд работ [19-21], где указаны специальные условия, при которых обобщённые интегральные операторы являются квазинильпотентными. Перейдём к полугруппам. Напомним, что экспоненциальным типом операторной полугруппы U(t) класса C0 называется величина w0 є limt®+Ґln∥U(t)∥ t . Указанный предел заведомо существует со значением в промежутке [-Ґ,+Ґ). В случае, когда w0 = -Ґ, полугруппа U(t) называется квазинильпотентной или ещё суперустойчивой. Мы будем использовать в основном второй термин, говоря именно про «суперустойчивость» полугруппы по аналогии с английским “superstability”, принятым в современной западной литературе (см. [10-15]). Данное название вполне оправдано, поскольку требование w0 = -Ґ эквивалентно тому, что при любом выборе числа a > 0 найдётся константа M = Ma і 1, для которой ∥U(t)∥Ј Me-at,t і 0. (9) Оценка (??) допускает дополнительную коррекцию при подходящем изменении нормы. Лемма 2. Пусть U(t) - суперустойчивая полугруппа класса C0. Тогда при любом выборе числа a > 0 найдётся эквивалентная норма ∥Ч∥a в пространстве E, такая, что ∥U(t)∥a Ј e-at,t і 0. (10) Доказательство. Воспользуемся известным приёмом [1, с. 376]. Зафиксируем число a > 0 и запишем оценку (??) с соответствующей константой M = Ma і 1. Положим ∥f∥a = suptі0∥eatU(t)f∥,f О E. (11) Учитывая (??), получаем, что ∥f∥Јв€Ґf∥a Ј M∥f∥ для любого f О E. Свойства нормы для ∥Ч∥a проверяются непосредственно. Итак, формула (??) действительно корректно определяет эквивалентную норму в E. Кроме того, если t і 0, то ∥U(t)f∥a = suptі0∥eatU(t + t)f∥ = sup tі0e-at∥ea(t+t)U(t + t)f∥ = = e-at sup sіt∥easU(s)f∥Ј e-at sup sі0∥easU(s)f∥ = e-at∥f∥ a,f О E, откуда и следует нужная оценка (??). Лемма доказана. Суперустойчивая полугруппа U(t) называется также квазинильпотентной, поскольку все её операторы U(t), взятые при t > 0, являются квазинильпотентными в пространстве E. При этом спектр s(A) порождающего оператора A будет пустым, а резольвентное множество r(A) заполнит всю комплексную плоскость. В частности, заведомо существует ограниченный обратный оператор A-1: E ® E, что несколько упрощает формулы при изучении поставленной задачи (??), (??) (см. [6]). Особым примером суперустойчивой (квазинильпотентной) полугруппы U(t) служит нильпотентная полугруппа, когда U(t) = 0,"t і t0 > 0, (12) с некоторым фиксированным значением t0 > 0. Нильпотентные полугруппы естественно возникают при рассмотрении процессов простого переноса в той или иной ограниченной области (см. [1, с. 552] или [5, с. 363-364]). 3. Формулировка основных результатов Вернёмся к рассмотрению обратной задачи (??), (??). Из теории, развитой в [6], следует такое утверждение (см. [6, с. 112]). Теорема 1. Пусть линейный замкнутый оператор A порождает суперустойчивую полугруппу U(t) класса C0. Пусть f О C[0,T] З BV [0,T] и m О BV [0,T] с выполненным условием (??). Пусть также b єт 0Tf(t)dm(t)№0. (13) Тогда обратная задача (??), (??) имеет и притом единственное решение (u(t),g) при любом выборе элементов u0,u1 О D(A). Метод доказательства теоремы 1, предложенный в [6], носил неконструктивный характер и базировался на теории Гельфанда коммутативных банаховых алгебр, точнее, на специальной версии этой теории, разработанной в трактате [1] для описания спектров интегралов от полугрупп. Дополнительный анализ показывает, что можно избежать обращения к таким сложным концепциям и получить теорему 1 элементарным конструктивным методом. В предположениях теоремы 1 рассмотрим операторное уравнение на неизвестный элемент g О E. Согласно схеме [6] уравнение допускает запись bg - Bg = f (14) со значением b№0 из формулы (??), линейным ограниченным оператором B = т 0Tf(0)U(t)dm(t) + т 0Tdm(t)т 0tU(t - s)df(s) (15) и заданным элементом f О E вида f = A т 0TU(t)u 0dm(t) - u1 . (16) Здесь и в дальнейшем все операторные интегралы типа (??) понимаются в сильной операторной топологии - со сходимостью на произвольном элементе h О E. Используя структуру выражения (??) и экспоненциальную оценку (??) с подходящим выбором значения a > 0, можно оценить оператор B в эквивалентной норме (??) и показать его квазинильпотентность. Наша цель - установить прямым методом следующий принципиальный факт. Теорема 2. Пусть выполнены предположения теоремы 1. Тогда оператор B из формулы (??) является квазинильпотентным в пространстве E. При этом решение однозначно разрешимого операторного уравнения (??) со значением b№0 представимо рядом Неймана g = е k=0Ґ 1 bk+1Bkf, (17) сходящимся по норме пространства E. Если элемент f О E взят в виде (??), а комплексное число b№0 - в виде (??), то элемент g из формулы (??) даёт второй компонент решения (u(t),g) поставленной обратной задачи (??), (??). При этом первый компонент решения - функция u(t) - выражается формулой (??). Теорема 2 есть главный результат нашей работы. Она обосновывает итерационный алгоритм для поиска решения обратной задачи (??), (??) при помощи ряда Неймана (??). Алгоритм приобретает законченный вид в специальном случае нильпотентной полугруппы U(t). Тогда, при некоторых дополнительных ограничениях на функции f(t),m(t), можно доказать такое утверждение. Теорема 3. Пусть выполнены предположения теоремы 1. Пусть, дополнительно, полугруппа U(t) является нильпотентной, удовлетворяя соотношению (??) с фиксированным значением t0 > 0. Предположим также, что функция m(t) является постоянной на некотором промежутке [0, a] М [0,T), а функция f(t) является постоянной на некотором промежутке [T - b,T] М (0,T], причём a + b > T. (18) Тогда оператор B из формулы (??) является нильпотентным и Bk = 0,k О в„•,k і t0 a + b - T. (19) Соответственно, решение операторного уравнения (??) представимо в виде конечной суммы g = е k=0N0 1 bk+1Bkf,N 0 є t0 a + b - T - 1, (20) замещающей бесконечный ряд Неймана (??). Здесь через вЊ€cвЊ‰ обозначен потолок числа c О в„ќ, т. е. наименьшее целое число, большее или равное c (см. [22, с. 88]). Теоремы 2 и 3 различны по характеру, и доказывать их удобно раздельно. 4. Доказательство основной теоремы 2 По поводу теории интеграла Стильтьеса см. [1, 23]. Без всяких пояснений используем стандартные оценки типа т 0Tv(t)dm(t) Ј max 0ЈtЈT∥v(t)в€ҐЧ Var{m(t)}0T, верные при любом выборе v О C([0,T],E) и m О BV [0,T]. Особо отметим следующий полезный факт. Лемма 3. Пусть f О C[0,T] З BV [0,T]. Тогда для любого заданного числа d > 0 можно подобрать такое значение t = t(f,d) > 0, что при всех t1,t2 О [0,T], удовлетворяющих условию 0 < t2 - t1 Ј t, выполняется соотношение Var{f(s)}t1t2 Ј d 2. (21) Доказательство. Определим функцию y(t) є Var{f(s)}0t,0 Ј t Ј T, монотонно возрастающую на [0,T]. Поскольку по условию f О C[0,T] З BV [0,T], то также и y О C[0,T] З BV [0,T] (см. [23, с. 210]). По теореме Кантора функция y(t) будет равномерно непрерывной на [0,T], т. е. для любого d > 0 найдётся такое значение t > 0, что при всех t1,t2 О [0,T], удовлетворяющих условию 0 < t2 - t1 Ј t, выполнено соотношение |y(t2) - y(t1)| = y(t2) - y(t1) = Var{f(s)}t1t2 Ј d 2. Но это и утверждается в (??). Лемма доказана. Оценим теперь отдельные слагаемые в основном операторе B из формулы (??). Все предположения теоремы 1, а, значит, и теоремы 2, относительно полугруппы U(t) и функций f(t), m(t) считаем выполненными. Напомним, что операторные интегралы рассматриваются (и оцениваются) в сильной операторной топологии - на произвольном элементе h О E, который для краткости исключаем из записей. Лемма 4. Пусть оператор B1 определён формулой B1 = т 0TU(t)dm(t). (22) Тогда для любого d > 0 найдётся такое a1 > 0, что при всех a > a1 верна оценка ∥B1∥a < d. (23) Здесь ∥Ч∥a - эквивалентная норма, заданная в пространстве E по формуле (??). Доказательство. Зафиксируем произвольное d > 0. В силу предположения (??) найдётся такое g О (0,T], что Var{m(t)}0g < d 2. (24) При этом, если g = T, то оценка (??) выполняется при любом значении a > 0, так как, согласно (??), имеем ∥B1∥a Ј max0ЈtЈTe-at Ч Var{m(t)} 0T = Var{m(t)} 0g < d 2 < d. Итак, при g = T значение a1 > 0 можно выбирать произвольно. Если же g О (0,T), запишем оператор B1 в виде B1 = т 0gU(t)dm(t) + т gTU(t)dm(t) = т 0gU(t)dm(t) + U(g)т gTU(t - g)dm(t). (25) Оценивая (??) с учётом (??) и (??), получаем ∥B1∥a Ј max0ЈtЈge-at Ч Var{m(t)} 0g + e-ag Ч max gЈtЈTe-a(t-g) Ч Var{m(t)} gT = = Var{m(t)}0g + e-ag Ч Var{m(t)} gT < d 2 + e-ag Ч Var{m(t)} 0T.(26) Выберем теперь a1 > 0 так, чтобы при всех a > a1 действовало соотношение e-ag Ч Var{m(t)} 0T < d 2. (27) Комбинируя (??) и (??), получаем нужную оценку (??). Лемма доказана. Для оценки второго слагаемого в операторе (??) требуется следующий предварительный результат. Лемма 5. Пусть операторная функция F(t) определена формулой F(t) = т 0tU(t - s)df(s),0 Ј t Ј T. (28) Тогда для любого d > 0 найдётся такое a2 > 0, что при всех a > a2 верна оценка ∥F(t)∥a < d,0 Ј t Ј T. (29) Здесь ∥Ч∥a - эквивалентная норма, заданная в пространстве E по формуле (??). Доказательство. Зафиксируем произвольное d > 0. Выберем значение t > 0 как в лемме 3, чтобы соотношение (??) выполнялось при всех t1,t2 О [0,T], удовлетворяющих условию 0 < t2 - t1 Ј t. При a > 0 и 0 Ј t Ј t, используя (??) и (??), получаем ∥F(t)∥a Ј max0ЈsЈte-a(t-s) Ч Var{f(s)} 0t = Var{f(s)} 0t Ј d 2 < d. Тем самым оценка (??) выполнена при любом a > 0 и 0 Ј t Ј t. Возможно, конечно, что t = T, и тогда все доказано с произвольным выбором значения a2 > 0. При 0 < t < T и t < t Ј T перепишем F(t) в виде F(t) = т 0t-tU(t - s)df(s) + т t-ttU(t - s)df(s),t < t Ј T. (30) Пусть снова a > 0. Оценивая (??) с учётом (??) и (??), получаем ∥F(t)∥a Ј max0ЈsЈt-te-a(t-s) Ч Var{f(s)} 0t-t + max t-tЈsЈte-a(t-s) Ч Var{f(s)} t-tt = = e-at Ч Var{f(s)} 0t-t + Var{f(s)} t-tt Ј e-at Ч Var{f(s)} 0T + d 2.(31) Выберем теперь a2 > 0 так, чтобы при всех a > a2 действовало соотношение e-at Ч Var{f(s)} 0T < d 2. (32) Комбинируя (??) и (??), получаем оценку (??) при t < t Ј T. В итоге с указанным выбором t > 0 и a2 > 0 нужный результат обоснован на всем отрезке [0,T]. Лемма доказана. Отметим, что, если f О C1[0,T], то доказательство леммы 5 сильно упрощается. Действительно, обозначив K є max0ЈsЈT|fў(s)|, для функции (??) имеем оценку ∥F(t)∥a Јт 0te-a(t-s)|fў(s)|ds Ј Kт 0te-a(t-s)ds < K a < d,0 Ј t Ј T, верную при a > a2 є K∕d с произвольным фиксированным d > 0. Приступим к доказательству самой теоремы 2. Оценим основной оператор B из формулы (??). Заметим, что B = f(0)B1 + B2, где оператор B1 определён формулой (??), а оператор B2 имеет вид B2 = т 0Tdm(t)т 0tU(t - s)df(s) = т 0TF(t)dm(t) (33) с функцией F(t) из формулы (??). Возьмём какое-то значение d > 0. Принимая во внимание предыдущие леммы 4 и 5 с их выбором a1 и a2, положим a0 = max(a1,a2). Тогда при всех a > a0 в эквивалентной норме ∥Ч∥a, заданной по формуле (??), получим, что ∥B∥a < |f(0)|Ч d + d Ч Var{m(t)}0T = |f(0)| + Var{m(t)} 0T d. Зафиксируем произвольное ε > 0. Ясно, что при достаточно малых d > 0 имеем оценку ∥B∥a Ј ε, верную при всех a > a0 c соответствующим значением a0 > 0. По лемме 1 это гарантирует квазинильпотентность оператора B. Но тогда операторное уравнение (??) при любом выборе числа b№0 будет однозначно разрешимым, и его решение g О E представимо рядом Неймана (??). То, что уравнение (??) с оператором (??) и заданным элементом (??) эквивалентно исходной обратной задаче (??), (??), следует из прежних результатов [6]. Теорема 2 полностью доказана. 5. Доказательство теоремы 3 Пусть функции f(t), m(t) удовлетворяют условиям теоремы 3, включая специальное требование (??). Тогда для операторов B1, B2 из формул (??), (??) возможны представления B1 = т aTU(t)dm(t) = U(a + b - T)U(T - b)т aTU(t - a)dm(t), B2 = т aTdm(t)т 0T-bU(t - s)df(s) = U(a + b - T)т aTdm(t)т 0T-bU(t - a + T - b - s)df(s). Следовательно, основной оператор B из формулы (??) выражается в виде B = f(0)B1 + B2 = U(a + b - T)Q, (34) где Q = f(0)U(T - b)т aTU(t - a)dm(t) + т aTdm(t)т 0T-bU(t - a + T - b - s)df(s). (35) Оператор U(a + b - T) коммутирует с оператором Q. По предположению теоремы 3 полугруппа U(t) является нильпотентной в пространстве E и удовлетворяет соотношению (??) с фиксированным значением t0 > 0. Поэтому Bk = U(k(a + b - T))Qk = 0,k О в„•,k(a + b - T) і t 0, что и даёт формулу (??). Итак, нильпотентность оператора B установлена. В силу (??) ряд Неймана (??) переходит в конечную сумму (??). Теорема доказана. Как видим, ключевую роль в доказательстве теоремы 3 сыграло специальное разложение (??) с оператором Q из формулы (??). Корректность конструкции обусловлена требованием (??). Отметим примеры, когда такая конструкция сильно упрощается. 6. Модельные примеры Всюду в данном пункте считаем, что полугруппа U(t) является нильпотентной в пространстве E, удовлетворяя условию (??) с фиксированным значением t0 > 0. Пусть f(t) є 1 на [0,T]. Этот особый случай часто возникает на практике. Поскольку функция f(t) является постоянной на любом промежутке [T - b,T] М (0,T], то для применения теоремы 3 достаточно требовать, чтобы функция m(t) была постоянной на промежутке [0, a] М [0,T) с неким (хоть каким-то) значением a О (0,T). В ответ можно всегда подобрать такое b О (0,T), что требование (??) окажется выполненным. Впрочем, здесь нет смысла прибегать к теореме 3, поскольку прямое исследование гораздо нагляднее. Действительно, зафиксируем значение a О (0,T) и рассмотрим обратную задачу uў(t) = Au(t) + g,0 Ј t Ј T, u(0) = u0,т aTu(t)dm(t) = u 1. (36) Стандартное требование (??) переходит в условие b є m(T) - m(a)№0. (37) Оператор (??) приобретает вид B = т aTU(t)dm(t) = U(a)т aTU(t - a)dm(t). (38) Поскольку U(ka) = 0 при k і t0∕a, то Bk = 0 при любом натуральном k і t0∕a, т. е. оператор (??) гарантированно является нильпотентным. Формула (??) для неизвестного элемента g О E сводится к конечной сумме g = е k=0N0 1 bk+1Bkf,N 0 є t0 a - 1, (39) с константой b№0 из формулы (??), оператором B из формулы (??), и элементом f = A U(a)т aTU(t - a)u 0dm(t) - u1 . (40) В случае, когда a і t0 (а такое вполне возможно), формулы (??)-(??) радикально упрощаются, и элемент g выражается одним слагаемым g = - 1 m(T) - m(a)Au1. (41) Любопытно, что такой ответ (??) зависит не от выбора той или иной конкретной функции m(t), а лишь от значения b є m(T) - m(a)№0. Подведём итог: именно соотношения (??)-(??) и надо использовать для нахождения решения обратной задачи (??) с оператором A, порождающим нильпотентную полугруппу U(t). Специальный интерес представляет частный случай обратной задачи (??), когда m(t) есть функция скачков, локализованных в промежутке [a,T] М (0,T]. Общее нелокальное условие (??) заменяется многоточечным условием, и обратная задача принимает вид uў(t) = Au(t) + g,0 Ј t Ј T, u(0) = u0,е kaku(tk) = u1. (42) Здесь {tk} - конечный или счётный набор различных точек из [a,T], занумерованных натуральным индексом k. Все числовые коэффициенты ak считаются отличными от нуля и такими, что е |ak| < Ґ. Требование (??) переходит в условие b єе kak№0. (43) Элемент g О E в обратной задаче (??) выражается прежней конечной суммой (??), где B = е kakU(tk),f = A е kakU(tk)u0 - u1 , (44) а значение b вычисляется через (??). Окончательное упрощение формул (??), (??) происходит в случае обратной задачи с финальным переопределением. Действительно, пусть f(t) є 1 на [0,T], а m(t) есть функция единичного скачка в финальной точке t = T. Тогда исходная обратная задача (??), (??) принимает вид uў(t) = Au(t) + g,0 Ј t Ј T, u(0) = u0,u(T) = u1. (45) Это и есть модельная обратная задача с финальным переопределением (см. [8]). Ввиду особой важности примера (??) разберём его более подробно. Следуя общим формулам (??), (??), (??), получаем b = 1,B = U(T),f = A(U(T)u0 - u1). (46) Поскольку Bk = U(kT) = 0 при любом натуральном k і t0∕T, то общая разрешающая формула (??) вновь сводится к конечной сумме, имеющей в данном случае особенно простой вид g = е k=0N0 U(kT)f,N0 є t0 T - 1, (47) с элементом f из (??). Возможно, конечно, что T і t0, и тогда N0 є вЊ€t0∕TвЊ‰- 1 = 0. При этом также U(T) = 0 и f = -Au1. В результате общий шаблон (??) даёт ответ g = -Au1. С другой стороны, возможно, что T < t0. Тогда, принимая во внимание представление (??) для элемента f и проводя элементарные преобразования, имеем g = е k=0N0 U(kT)f = f + U(T)f + U(2T)f + † + U(N0T)f = = A U(T)u0 + U(2T)u0 + U(3T)u0 + † + U N0 + 1Tu0 - - A u1 + U(T)u1 + U(2T)u1 + † + U N0Tu1 = = -Au1 + Aе k=1N0 U(kT) u0 - u1 ,N0 є t0 T - 1. Слагаемое U N0 + 1Tu0 пропадает, ибо N0 + 1T і t0 и U N0 + 1T = 0. В итоге формулу (??) для второго компонента g решения обратной задачи (??) можно записать в виде g = -Au1, 0 < t0 Ј T, -Au1 + Aе k=1N0 U(kT) u0 - u1 ,t0 > T > 0, (48) где N0 є вЊ€t0∕TвЊ‰- 1. Выражение (??) удобно на практике. Для последующих применений данного результата зафиксируем точное утверждение. Теорема 4. Пусть линейный замкнутый оператор A порождает нильпотентную полугруппу U(t) класса C0, удовлетворяющую условию (??) с фиксированным значением t0 > 0. Тогда при любом выборе u0,u1 О D(A) элемент g О E, разрешающий обратную задачу (??), находится по формуле (??). При этом функция u(t) восстанавливается по формуле u(t) = U(t)u0 + т 0tU(s)gds,0 Ј t Ј T, (49) непосредственно следующей из (??) в случае, когда f(t) є 1. Формула (??) используется также в обратных задачах (??) и (??). 7. Заключение Итак, теория линейной обратной задачи (??), (??) с оператором A, порождающим суперустойчивую полугруппу, оказывается в основном завершённой. Полученные результаты представляют интерес для некоторых специальных разделов математической физики, связанных с теорией упругости и задачами линейного переноса. Например, на основе формул, установленных нами в п. 6, можно разработать конкретные алгоритмы решения обратных задач типа (??) для многомерного уравнения переноса, взятого в среде без столкновений, но с возможным поглощением частиц. Этот материал полезно опубликовать отдельно.
Об авторах
Иван Владимирович Тихонов
МГУ имени М. В. Ломоносова
Автор, ответственный за переписку.
Email: ivtikh@mail.ru
доктор физико-математических наук, профессор кафедры математической физики факультета ВМК МГУ имени М. В. Ломоносова
Ленинские горы, ГСП-1, Москва, Россия, 119991Шон Тунг Ву Нгуен
Московский педагогический государственный университет
Email: vnsontung@mail.ru
аспирант кафедры математического анализа Московского педагогического государственного университета
ул. Краснопрудная, д. 14, Москва, Россия, 107140Список литературы
- Хилле Э., Филлипс Р. Функциональный анализ и полугруппы. — М.: ИЛ, 1962.
- Данфорд Н., Шварц Д. Линейные операторы. Т. 1. Общая теория. — М.: ИЛ, 1962.
- Крейн С. Г. Линейные дифференциальные уравнения в банаховом пространстве. — М.: Наука, 1967.
- Pazy A. Semigroups of Linear Operators and Applications to Partial Differential Equations. — N.Y.: Springer Verlag, 1983.
- Engel K.-J., Nagel R. One-Parameter Semigroups for Linear Evolution Equations. — N.Y.: Springer, 2000.
- Тихонов И. В., Эйдельман Ю. С. Вопросы корректности прямых и обратных задач для эволюционного уравнения специального вида // Математические заметки. — 1994. — Т. 56, вып. 2. — С. 99–113.
- Прилепко А. И., Тихонов И. В. Восстановление неоднородного слагаемого в абстрактном эволюционном уравнении // Известия РАН. Сер. матем. — 1994. — Т. 58, вып. 2. — С. 167–188.
- Prilepko A. I., Orlovsky D. G., Vasin I. A. Methods for Solving Inverse Problems in Mathematical Physics. — N.Y.: Basel, 2000.
- Тихонов И. В., Тунг В. Н. Ш. Метод решения обратной задачи для эволюционного уравнения с суперустойчивой полугруппой // Дифференциальные уравнения и процессы управления. — 2017. — Вып. 2. — С. 51–58.
- 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.
- 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.
- Balakrishnan A. V. Superstability of Systems // Applied Mathematics and Computation. — 2005. — Vol. 164, issue 2. — Pp. 321–326.
- Creutz D., Mazo M., Preda C. Superstability and Finite Time Extinction for C0– Semigroups // arXiv:0907.4812. Submitted. — 2013. — Pp. 1–12.
- 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.
- Kmit I., Lyul’ko N. Perturbations of Superstable Linear Hyperbolic Systems // arXiv:1605.04703. Submitted. — 2017. — Pp. 1–26.
- Приближенное решение операторных уравнений / М. А. Красносельский, Г. М. Вайникко, П. П. Забрейко и др. — М.: Наука, 1969.
- Malinen J., Nevanlinna O., Zem´anek J. Microspectral Analysis of Quasinilpotent Operators // arXiv:1211.4790v1. Submitted. — 2012. — Pp. 1–25.
- Eskandari R., Mirzapour F. Hyperinvariant Subspaces and Quasinilpotent Operators // Bulletin of the Iranian Mathematical Society. — 2015. — Vol. 41, issue 4. — Pp. 805– 813.
- Забрейко П. П. О спектральном радиусе интегральных операторов Вольтерра // Литовский математический сборник. — 1967. — Т. 7, вып. 2. — С. 281–287.
- Жуковский Е. С. К теории уравнений Вольтерра // Дифференциальные уравнения. — 1989. — Т. 25, вып. 9. — С. 1599–1605.
- Сумин В. И., Чернов А. В. Операторы в пространствах измеримых функций: вольтерровость и квазинильпотентность // Дифференциальные уравнения. — 1998. — Т. 34, вып. 10. — С. 1402–1411.
- Грехэм Р., Кнут Д., Паташник О. Конкретная математика. — М.: Мир, 1998.
- Натансон И. П. Теория функций вещественной переменной. — М.: Наука, 1974.