О задаче идентификации функции источника специального вида в двумерном параболическом уравнении Текст научной статьи по специальности «Математика»

УДК 517.9

О задаче идентификации функции источника специального вида в двумерном параболическом уравнении

Игорь В. Фроленков* Екатерина Н. Кригер^

Институт математики, Сибирский федеральный университет, Свободный, 79, Красноярск, 660041,

Россия

Получена 18.07.2010, окончательный вариант 25.08.2010, принята к печати 10.09.2010 В 'работе доказаны существование, единственность и устойчивость по входным данным 'решения задачи идентификации функции источника специального вида в параболическом уравнении в случае данных Коши.

Ключевые слова: задачи идентификации коэффициентов, обратные задачи, уравнения в частных производных, метод слабой аппроксимации, устойчивость решения.

В работе рассмотрена задача идентификации функции источника специального вида в параболическом уравнении в случае данных Коши. Существование и единственность решения задачи доказаны в классах гладких ограниченных функций. Для исследования задачи применяется метод, позволяющий, используя условия переопределения, привести исходную обратную задачу к прямой задаче для нагруженного уравнения [1, 2]. Существование решения прямой задачи доказано методом слабой аппроксимации [3, 4].

Задача идентификации функции источника, зависящей от (t, x), рассмотрена в работе [5]. В [2] изучен случай, когда условия переопределения задаются на гладкой кривой, а неизвестная функция источника зависит лишь от временной переменной. Краевые задачи идентификации функции источника исследованы в работах [6, 7].

В полосе G[o,t] = {(t,x, z) | 0 ^ t ^ T, (x,z) £ R2} рассматривается задача Коши для параболического уравнения

ut = uxx + uzz + f (t, x, z) • A(t, x, z), t £ (0, T), (x, z) £ R2, (1)

с начальным условием

u(0, x, z)= uo(x, z), (2)

где наряду с функцией u(t, x, z) нужно определить также функцию A(t, x, z), причем

A(t, x, z) = Ai(t,x) + A2(t, z). (3)

Пусть заданы условия переопределения:

u(t, x, a) = y(t, x), (4)

* kspk_job@mail.ru te_katherina@mail.ru © Siberian Federal University. All rights reserved

и(г,в,г) = ф(г,г),

где а, в — некоторые фиксированные постоянные. Считаем выполненными условия согласования:

по(х, а) = р(0, х),

по (в,г) = ф(0,г),

ф,в) = Ф(г,а),

и условия на функцию /(г,х,г):

\/(г,в,г)\ > ¿1 > о, I/(г,х,а)\ > ¿2 > о, Vг е [о,т], V(х,г) е м2,

(5)

(6)

(7)

(8)

(9)

где ¿1, ¿2 — некоторые постоянные.

Приведем данную обратную задачу к вспомогательной прямой. Для этого в уравнении (1) сначала положим х = в и выразим из получившегося равенства Л1(г,в), затем положим г = а и выразим Л2(г, а), далее возьмём х = в и г = а. В полученное соотношение подставим выражения на функции Л1(г,в) и Л2(г,а) и выразим Л(г,х,г).

Х(г х г) = Щ(г,в,г) - Пгг(г,в,г) + п(Ь, х, а) - пхх(г, х, а)

/(г, в, г)

/(г, х, а)

щ(г,в,а) - Пхх(г, в, а) - пгг(г, в, а) Пхх(г, в, г) игг(г,х,а)

/(г, в, а)

/(г, в, г) /(г,х,а)

Затем, используя условия переопределения (4),(5), получим выражение на неизвестный коэффициент при функции /(г, х, г) в следующем виде:

Л(Ь, х, г) = д(г, х, г) -

Пхх(г,в,г) пхх(Ь, х, а)

/(г, в, г) /(г,х,а)

(10)

где

д(г х г) = , г) - Ф*г(г, г) + ^(г, х) - Фхх(Ъ, х)

/(г, в, г)

/ (г, х, а)

Фг(г, а) - фхх(Ъ, в) - Фгг(г, а)

/(г, в, а)

(11)

— известная функция.

Теперь подставим выражение (10) в уравнение (1) и перейдем к прямой задаче для уравнения

П = пхх + пгг + /(г, х, г)

,(г,в,г) пгг (г,х,а)

/(г, в, г) /(г,х,а)

+ о(г, х, г),

г е (о,т), (х,г) е м2, (12)

с начальным условием (2).

Здесь 0(г, х, г) = /(г, х, г) ■ д(г, х, г) — известная функция. Пусть

Ф1(г,х,г) = ф^г), Ф2(г,х,г)= /(г,х,г)

/(г, в, г)

/ (г, х, а)

п

Относительно входных данных предположим, что они достаточно гладкие, имеют все непрерывные производные, входящие в соотношение (13), и удовлетворяют ему.

д 3к1

дг дхк1

у (г, х)

+

д дк2

дг дгк2

^(г,г)

+

дк1 дк2

дхк1 дгк2

мо(х,г)

+

дк1 дк2

/(г, х, г)

< С,

дхк1 дгк2'

кьк2 =0,1,..., 6. (13)

Докажем разрешимость прямой задачи.

Для доказательства существования решения прямой задачи будем использовать метод слабой аппроксимации [3]. Расщепим уравнение (12) на четыре дробных шага и сделаем сдвиг по времени на 4 в следах неизвестных функций

«Г = 4 • (иХх + ), пг <г < (п + з)т;

«Т = 4 • иХх (* - 4, в, г) • Ф^г, х, г), (п + 1) т < 4 < (п + 2) т;

«Т = 4 • - 4,х, а • Ф2(г,х, г), (п + 2) т < г < (п + |) т;

= 4 • С(г, х, г), (п + 4) т < г < (п + 1)т, мт(0, х, г) = мо(х, г), п = 0,1, 2,..., (Ж - 1), Жт = Т.

(14)

(15)

(16)

(17)

(18)

Для решения мт расщепленной задачи (14)—(18) доказаны равномерные по т оценки:

дк1 дк2

дхк1 дгк2

м' (г, х, г)

< С, кь к2 =0, 1, .. ., 6, (г, х, г) € С[0,Т].

(19)

Из оценок (19) следует, что правые части уравнений (14)—(17) ограничены равномерно по т на любом временном шаге, а значит, и левые части уравнений будут равномерно ограничены по т.

(г, х, г)| ^ С, (г, х, г) € С[о,т].

Дифференцируя уравнения (14)—(17) по х и г, в силу (19) получим оценки

< С, ^ ,к2 =0,1,..., 4. (20)

дк1 дк2

Оценки (19), (20) гарантируют выполнение условий теоремы Арцела о компактности. В силу теоремы Арцела некоторая подпоследовательность иТк (г, х, г) последовательности мт (г,х, г) решений задачи (14)—(18) сходится вместе с производными по х и по г до четвертого порядка включительно к функции м(г, х, г) € С0,4'4 (С[о,т]), которая в силу теоремы сходимости метода слабой аппроксимации [3] является решением задачи (12), (2), причем

«^х г) € (С[0'Т^, где

С^'Х.г'2 (С[0'Т]) Н «(г, х, г)

дк д к1 + к2

€ С (С[о'Т0 ,

дгк ' дхк1 дгк2

к = 0,1,...,/, ^

0,1,...,/ьк2 =0,1,...,/2

При этом справедливы следующие оценки при (г, х, г) е ]: дк1 дк2

дхк1 дгк2

ь(Ь, х, г)

< С, к1,к2 =0,1,..., 4.

(21)

Теперь докажем, что пара функций и(г, х, г), Л(г,х, г), где Л(г,х,г) определяется соотношением (10) и удовлетворяет условию (3), является решением обратной задачи. Так как и(г, х, г) — решение прямой задачи, то при подстановке и(г, х, г) в (1) и (2) получим тождества.

Из (10) и (12), с учетом условий (9), (13), (21), следует, что и(г,х,г), Л(г,х,г) удовлетворяют неравенству

4 4

ЕЕ

к1 = 0 к2 = 0

дк1 д)к2

дхк1 дгк2

/,(г, х, г)

22 к1 =0 к2 = 0

дк1 дк2

дхк1 дгк2

Л(г, х, г)

< С.

(22)

Докажем выполнение условия переопределения (4). Для этого положим в (12) г = а. Имеем:

иг(г, х, а) = пхх(г, х, а) + ихх(г, х, а) + /(г, х, а)

д(г, х, а) -

Пхх(г,в,а) игг (г,х,а)

/ (г, в, а) / (г,х,а)

где функция д(г,х,а) известна и задана соотношением (11). Преобразуем выражение следующим образом:

иг(г, х, а) = Пхх(г, х, а) + ^(г, х) - <Рхх(г, х) + (фхх(г, в) - Пхх(Ь, в, а)) Введём обозначение:

у1(г, х) = и(г, х, а) - ^(г, х). В силу (6) справедливо следующее равенство:

V1 (0,х) = и(0, х, а) - <р(0,х) = 0.

/(г, х, а) /(г, в, а).

(23)

Используя обозначения (23), получим

(г,х) = х(г,х) -

/ (г, х, а)

™ /(г, в, а)

v1(0, х) = и(0, х, а) - р(0, х) = 0.

Vlx(t,в),

(24)

Рассмотрим в области П[0,т] = {(г, у) | 0 ^ г ^ Т,у е М} задачу Коши для уравнения

vt(г, у) = vyy (г, у) - s(г, у) ■ Иуу (г, п)

с начальным условием

v(0,y) = г (у).

(25)

(26)

Предполагаем, что функции в(г, у), г(у) и все их производные непрерывны и ограничены в П[0,т], П — фиксированная вещественная постоянная. Доказана

Лемма 1. Если решение у) € С1,4 (П[0;Т]) задачи (25), (26) существует и удовлетворяет условию

Е

й=0

то оно единственно.

д^

^ С,

Здесь

С1,4 (П[о,т]) = |

дй I

ду^ ^^ € С (П[0,Т]) , к = 0,1,..., 4

В силу леммы 1 решение задачи (24) единственно. Очевидно, что решением является функция и1^, х) = 0.

Таким образом, из (23) следует выполнение условия (4).

Теперь докажем условие переопределения (5). Положим в (12) х = в. Получим:

в, г) = в, г) + (I, в, г) + /(I, в, г)

д(г в г) _ «жж(!,в,*0 (I, в, а)

/(1,в,г) /(1,в,а) ]

где функция д(1, в, г) известна и задана соотношением (11). Используя условия согласования, получим уравнение

/(I, в, г)

^(1, в, г) = (I, в, г) + г) _ ^(I, г) _ (I, в, а) _ ^(I, а)) • /(¿^д,). Обозначим

и2(1,,г) = и(1,в,г) _ ^(М). (27)

Тогда из (7) следует равенство

и2(0, х) = и(0, в, г) _ ^(0, г) = 0.

Используя обозначения (27), получим задачу

2(+ \ 2 /, \ /(1,в, г) 2 \ V2 (I, г) = (г, г)--——-—- • (I, а),

а ; 2 /(г,в,а) ^ ' (28)

и2(0, х) = 0.

Решение задачи (28) единственно (в силу леммы 1). И таким решением является функция и2(1, г) = 0.

Отсюда, учитывая (27), следует выполнение условия (5).

Итак, условия переопределения выполнены. Следовательно, пара функций и(1, х, г), А(1, х, г) является решением исходной задачи (1)—(5). Доказана

Теорема 1. Пусть выполняются условия (3), (9), (13). Тогда существует решение и(1, х, г), А(1, х, г) обратной задачи (1)-(8) в классе

2(Т) = {и(1,х,г),А(1,х,г) и € С^;4(С[о,тО,А(1,х,г) € С/Д'^С^т])},

удовлетворяющее соотношению (22).

Докажем единственность решения задачи (1)-(5) при выполнении условий (9), (13), (22). Пусть {и(г,х, г), Л(г,х, г)} и {и(г,х, г), Л(г,х, г)} — два классических решения задачи (1)-(5). Пара функций и(г, х, г), Л(г, х, г) — решение, определяемое соотношением (10) и удовлетворяющее условию (3), а пара функций и(г,х,г),Л(г,х,г) = Л1(г,х) + Л2(г,г) —некоторое другое решение задачи (1)-(5), удовлетворяющее условию (22). Тогда справедливы следующие соотношения:

иг(г, х, г) = ихх(г, х, г) + ихх(г, х, г) + /(г, х, г) ■ Л(г, х, г), (29)

иг(г, х, г) = Лхх(г, х, г) + ихх(г,х, г) + /(г,х, г) ■ Л(г,х, г), (30)

и(0, х, г) = и0(х, г), и(0, х, г) = и0(х, г),

и(г,х,а) = р(г,х), и(г,х,а) = р(г,х), (31)

и(г, в, г) = Ф(г, г), и(г, в, г) = Ф(г, г). Пусть

л(Ь, х, г) - и(г, х, г) = ш(Ь, х, г), Л(г, х, г) - Л(г, х, г) = 7(г, х, г) = 71(г, х) + 72(г, г). Пара функций т(г, х, г), 7(г, х, г) является решением задачи Коши:

тг(г, х, г) = тхх(г, х, г) + тгг(г, х, г) + /(г, х, г) ■ 7(г, х, г), (32)

ю(0, х, г) = 0,

т(г, х, а) = 0, (33)

■ю(Ь, в, г) = 0.

Полагаем в уравнении (32) х = в, г = а. Используя (33), выражаем коэффициент при функции /(г, х, г). Подставляя его выражение в (32), получим

/ шхх(г,в,г) (г,х,а)\

'Шг(г, х, г) = 'Шхх(г, х, г) + (г, х, г) - / (г, х,г) ■ —г- + —т~-г- , (34)

у / (г, в, г) / (г,х,а) у ю(0,х,г) = 0. (35)

Доказано, что т(г,х,г) = 0 в С^т], то есть и(г,х,г) = Л(г,х,г) в С[0,т]. Из (32) и (33) для 7(г,х,г) следует выполнение соотношения

/(г, х, г) ■ 7(г,х,г)=0. (36)

Рассматривая (36) в точках х = в и г = а, в С[0,т] получаем справедливость равенств:

/(г, в, г) ■Ыг,в)+ 12(г,г))=0, /(г, х, а) ■ (71 (г, х) + 72 (г, а)) =0,

/(г,в, а) ■ Ыг,в) + ъ(г,а)) =0.

В силу (9)

и(г,в)+ 12(г,г)=0, (37)

ц(г,х)+ 72(г,а) = 0, (38)

ц(г,в)+ 72(г,а)=0. (39)

Сложим уравнения (37) и (38).

71(1, х) + 72(1, г) + 71(г, в) + 72(1, а) = 0.

Из (39) следует, что

7(1,х, г) = 71(1, х) + 72 (I, г) = 0, А(г, х, г) _ Л(г, х, г) = 0, I € [0, Т],

А(г,х,г) = А(г,х, г), I € [0,Т].

то есть

а следовательно,

Доказана

Теорема 2. Решение и(!, х, г), А(!, х, г) задачи (1)-(9), удовлетворяющее (22), единственно в классе 2 (Т).

Из теорем 1 и 2 следует

Теорема 3. Пусть выполняются условия (6)-(9), (13). Тогда в классе 2(Т) существует единственное решение и(!, х, г), А(!, х, г) задачи (1)-(5), удовлетворяющее соотношению (22).

Докажем непрерывную зависимость решения задачи (1)-(5) от входных данных.

Рассмотрим в С[о,т] две задачи Коши:

= + + /г(!, х, г) • Аг(!, х, г), г = 1, 2, А4(!,х,,г) = А!(!,х) + А|(М), мг(0, х, г) = и0(х, г), иг(!, х, а) = у>г(!, х), и4(!,в,г) = ^(М),

где а, в — некоторые фиксированные постоянные.

Введем обозначения:

и = и1 _ и2, л = л1 + л2 = а1 _ а2 + а2 _ а2,

12 хо _ ^

^ = /1 _ /2, ио = ио _ «2, Ф = ^ _ ^ ф = ^ _ ^

(40)

Вычтем из задачи при г =1 задачу при г = 2. Согласно обозначениям (40) получим задачу:

и = + игг + ^ • А1 + /2 • Л, (41)

и(0, х, г) = ио(х,г), (42)

и(I, х, а) = Ф(!, х), (43)

и (!,в,г) = Ф(!,г), (44)

Полагая в уравнении (41) х = в, г = а и используя (43) и (44), получаем выражение для неизвестного коэффициента

д/, ^ в, г) (¿, х, а) Л(£, х, г) =--т,. п—^---¡^-Г +

f 2(t,e,z) f 2(t,x,a)

Ф^, z) - (t, z) - F(t, в, z) ■ A1(i, в, z) +

+ лО ^ ^ \ +

+

f 2(t,e,z)

Ф4(М) - Фжж(^ж) - F(t,x, a) ■ A1(t, x, a)

f 2(t, x, a)

Ф^, a) - Фяа!& в) - (t, a) + F(t, в, a) ■ A1(t, в, a) f 2M, a)

. (45)

Подставляя выражение (45) в уравнение (41), получаем задачу Коши для уравнения и4(г, х, г) = х, г) + (£, х, г) + ^(£, х, г) • А1^, х, г) +

+ f2(t,x,z) ■ -

Uxx(t, в, z) Uzz(t,x, a)

f 2(t,^z) f 2(t, x, a)

+

Ф^, z) - (t, z) - F(t, в, z) ■ A1(t, в, z) +

+ rO / , r> \ +

+

f 2(t,в,z)

Ф^,х) - ФХХ(М) - F(t,x, a) ■ A1(t, x, a)

f 2(t, x, a)

ФеМ - Ф**М) - (t, a) + F (t,в, a) ■ A1 M,a)4

f2 (t, в, a)

(46)

с начальным условием (42).

Введем следующие обозначения:

4 4

||D1(t,x,z)||1 = ЕЕ SUP SUP

d kl + fc2

dxkl dzk2

A(£,x,z)

44

||d2(.

x, z)|U =

||Ds(t,.

||D4(t,:

E E, sup

k1=0 fc2 = 0 (x>z 1 6

sup sup

^^ ^^ zeR

d fcl + fc2

i=0 j=0 16

dxkl dzk2

D2(x, z)

ЕЕ sup вир

i=0 j=0

d^dzJ'

dt^xJ'

Ds(^,z)

■D4(e,x)

(47)

Функции ^, г = 1, 2,..., 4, из (47) и их производные, входящие в (47), ограничены и непрерывны в .

Для всех £ € (0, Т] согласно обозначениям (47) доказана равномерная по т оценка

l|U||1 < C -g|U0||2 + ||F||1 + ||Ф||з + ||Ф||^. Из (13) и (45) следует оценка

l|U||1 + ||Л|| 1 < C ■ (1IU0II2 + ||FII1 + ||Ф||з + ||Ф||4

з

4

то есть

Цп1 - п2||1 + 11Л1 - Л2||1 <

< С ■ (Ш - п0||2 + ||/ - /21|1 + и^1 - ||з + И/ - Л). (48)

Теорема 4. При выполнении условий (6)-(9), (13) для 'решения и(г,х,г), Л(г,х,г) задачи (1)-(5) выполняется оценка (48) устойчивости по входным данным.

Список литературы

[1] Ю.Я.Белов, И.В.Фроленков, Некоторые задачи идентификации коэффициентов полулинейных параболических уравнений, Докл. РАН, 404(2005), №5, 583-585.

[2] Ю.Я.Белов, И.В.Фроленков, О задаче идентификации двух коэффициентов параболического полулинейного уравнения с условиями переопределения, заданными на гладкой кривой, Вычислительные технологии, 11(2006), №1, 46-54.

[3] Ю.Я.Белов, С.А.Кантор, Метод слабой аппроксимации, Красноярск, КрасГУ, 1999.

[4] Н.Н.Яненко, Метод дробных шагов решения многомерных задач математической физики, Новосибирск, 1967.

[5] О.А.Афиногенова, Ю.Я.Белов, И.В.Фроленков, О стабилизации решения задачи идентификации функции источника одномерного параболического уравнения, Докл. РАН, 424(2009), №4, 439-441.

[6] A.I.Prilepko, D.G.Orlovsky, I.A.Vasin, Methods for solving inverse problems in mathematical physics, New York: Marcel Dekker, 2000.

[7] О.Н.Черепанова, Т.Н.Шипина, Об одной задаче идентификации функции источника в параболическом уравнении, Журнал СФУ. Математика и физика, 2(2009), №3, 370-375.

An Identification Problem of the Source Function of the Special Form in Two-Dimensional Parabolic Equation

Igor V. Frolenkov Ekaterina N. Kriger

The existence, uniqueness and stability of solution by input data of the identification problem for parabolic equation with source function of the special form in the case of Cauchy's data has been proved in this article.

Keywords: problem of the identification of source function, inverse problem, equations in partial derivatives, method of weak approximation, solution stability.