Научная статья на тему 'Обратные задачи для уравнения теплопроводности'

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

CC BY
1231
227
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ / ПЕРВАЯ НАЧАЛЬНО-ГРАНИЧНАЯ ЗАДАЧА / ОБРАТНЫЕ ЗАДАЧИ / СПЕКТРАЛЬНЫЙ МЕТОД / ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ / ЕДИНСТВЕННОСТЬ / СУЩЕСТВОВАНИЕ / УСТОЙЧИВОСТЬ

Аннотация научной статьи по математике, автор научной работы — Зайнуллов А. Р.

На основании формулы решения первой начально-граничной задачи для неоднородного уравнения теплопроводности изучены обратные задачи по отысканию начального условия и правой части. Методом спектрального анализа установлен критерий единственности решения обратной задачи по отысканию начального условия. Правая часть уравнения теплопроводности представлена в виде произведения двух функций, одна из них зависит от пространственной координаты, другая от времени. В одной задаче наряду с неизвестным решением ищется множитель правой части, зависящий от времени, а в другой множитель, зависящий от пространственной координаты. Для этих задач доказаны теоремы единственности, существования и устойчивости решения.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Inverse problems for the heat equation

The inverse problem of finding initial conditions and the right-hand side had been studied for the inhomogeneous heat equation on the basis of formulas for the solution of the first initial-boundary value problem. A criterion of uniqueness of solution of the inverse problem for finding the initial condition was found with Spectral analysis. The right side of the heat equation is represented as a product of two functions, one of which depends on the spatial coordinates and the other from time. In one task, along with an unknown solution is sought factor on the right side, depending on the time, and in another a factor that depends on the spatial coordinates. For these tasks, we prove uniqueness theorems, the existence and stability of solution.

Текст научной работы на тему «Обратные задачи для уравнения теплопроводности»

УДК 517.95

А.Р. Зайнуллов1

ОБРАТНЫЕ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

На основании формулы решения первой начально-граничной задачи для неоднородного уравнения теплопроводности изучены обратные задачи по отысканию начального условия и правой части. Методом спектрального анализа установлен критерий единственности решения обратной задачи по отысканию начального условия. Правая часть уравнения теплопроводности представлена в виде произведения двух функций, одна из них зависит от пространственной координаты, другая — от времени. В одной задаче наряду с неизвестным решением ищется множитель правой части, зависящий от времени, а в другой — множитель, зависящий от пространственной координаты. Для этих задач доказаны теоремы единственности, существования и устойчивости решения.

Ключевые слова: уравнение теплопроводности, первая начально-граничная задача, обратные задачи, спектральный метод, интегральное уравнение, единственность, существование, устойчивость.

1. Постановка задач

Рассмотрим уравнение теплопроводности

Ьи = щ - а ихх = Р(х,г) = / (х)д(г) (1.1)

в прямоугольной области

б = {(х,г) | о <х<1, о <г <т}

и следующую начально-граничную задачу.

Первая начально-граничная задача. Найти в области Б функцию и(х,г), удовлетворяющую условиям:

и(х,г) е с (Б) п СХ(Б) п сХ:1(Б), (1.2)

Ьи = р(х,г), (х,г) е Б, (1.3)

и(о,г) = и(1,г) = о, о < г < т, (1.4)

и(х, о) = ф(х), о < х < I, (1.5) где ф(х) и Р(х, г) — заданные функции. !© Зайнуллов А.Р., 2015

Зайнуллов Артур Рашитович (arturzayn@mail.ru), кафедра математического анализа, Стер-

литамакский филиал Башкирского государственного университета, 453103, Российская Федерация, г. Стерлитамак, пр. Ленина, 49.

На основе этой прямой задачи для уравнения (1.1) рассмотрим следующие обратные задачи.

Обратная задача 1. Найти функции u(x,t) и ф(х), удовлетворяющие условиям (1.2)—(1.5) и, кроме того, дополнительному условию

u(xo,t) = h(t), 0 <t0 < t < ti < T, (1.6)

где xo — заданная фиксированная точка отрезка (0,l), to и t1 — заданные действительные числа, F(x,t), h(t) — заданные функции.

Отметим, что указанная обратная задача приведена в книгах [1, с. 119] и [2, с. 248] и для уравнения (1) при F(x,t) = 0, ux(0,t) = ux(l,t) = 0, 0 < t < T доказана теорема ее единственности при xo = 0 и xo = l/n.

В данной работе установлен критерий единственности решения (1.2)—(1.6) для любой точки xo € (0,l).

Обратная задача 2. Найти функции g(t) и u(x,t), удовлетворяющие условиям (1.2)—(1.5) и, кроме того, дополнительному условию

u(x0,t) = h(t), 0 < t < T, (1.7)

где x0 — заданная фиксированная точка отрезка [0, l], ф(x), h(t) и f (x) — заданные достаточно гладкие функции, при этом ф(xo) = h(0).

Обратная задача 2 при ф(x) = 0, ux(0,t) = ux(l,t) = 0, 0 ^ t ^ T изучена в [1, с. 123; 2, с. 250] и доказана теорема ее единственности и существования при f(xo) = 0.

Здесь получены новые достаточные условия единственности и существования решения обратной задачи 2.

Обратная задача 3. Найти функци f (x) и u(x,t), удовлетворяющие условиям (1.2)—(1.5) и, кроме того, дополнительному условию

u(x,t0) = ф(x), 0 < x < l, (1.8)

где to — заданная фиксированная точка отрезка (0, T], ф^), ф(x) и g(t) — заданные функции.

В данной статье с введением дополнительного условия (1.8) получены теоремы единственности, существования и устойчивости решения обратной задачи 3.

2. Критерий единственности решения обратной задачи 1

Решение прямой задачи (1.2)—(1.5) может быть получено методом разделения переменных [3, с. 200]. Решение уравнения (1.1), удовлетворяющее условию (1.4), будем искать в виде u(x,t) = X(x)T(t) как сумму решений неоднородного и однородного уравнений. Подставляя последнее в однородное уравнение (1.1), т. е. при F(x,t) = 0, получим относительно X(x) спектральную задачу:

X"(x) + XX(x)=0, 0 < x < l, X = const,, (2.1)

X(0) = X(l) = 0. (2.2)

Как известно, задача (2.1)-(2.2) — это задача Штурма-Лиувилля, которая имеет счетное множество собственных значений Xk = (nk/l)2, k € N, и соответствующую систему собственных функций Xk (x) = sin ^x, цк = nk/l. Эта система ортогональна и полна в пространстве L2 [0, l] и поэтому образует в нем ортогональный базис.

Решение неоднородного уравнения (1.1) по системе собственных функций Xk (x) сводится к решению дифференциального уравнения

Tk(t)+ ш\Тк (t) = Fk (t), 0 < t < Т.

Тогда решение задачи (1.2)-(1.5) определяется в виде суммы ряда

u(x,t) = ^2, Фк 1 sin Hk X + ^ Fk (t)sin Цк x, (2.3)

k=1 k=1

где

Фк = 2 j Fk (t) = fk j g(ra"(t-T>dr, (2.4)

2 íl nk fk = j J f (£)sin Ик£d£, pk = k G N.

Следуя [3, с. 210; 4, c. 264], можно показать справедливость следующего утверждения.

Теорема 1. Если ф(х) G C 1[0,j], ф(0) = ф(1) = 0, f(x) G C2[0,l], f(0) = f(l) = = 0, g(t) G C[0,T], то существует единственное решение задачи (1.2)-(1.5), и оно определяется в виде суммы ряда (2.3), где коэффициенты находятся по формулам (2.4).

Рассмотрим обратную задачу (1.2)-(1.6), т. е. задачу 1. Полагая в формуле (2.3) х = хо, с учетом условия (1.6) получим уравнение относительно неизвестной функции ф(х):

ж 0 „ i

Е 2 2 21 j Ф>(£) sin pkt sin pkхо =

k=i

ж

= h(t) -J2 Fk(t)sinpkхо = ho(t), to < t < ti. (2.5)

k=1

Теорема 2. Если sin пкх0 = 0 при всех k G N, где ж0 = х0/1, то решение интегрального уравнения (2.5) единственно в L2[0,l].

Доказательство проведем, следуя [1, с. 119]. В силу линейности уравнения (2.5) достаточно показать, что оно имеет только нулевое решение при ho(t) = 0. Положив в (2.5) h(t) =0 и F(х, t) = 0, имеем

ж „ i

V 2 ф(0 sin pk£d£e-^a2t sin пкх0 = 0. (2.6)

k=i j0

Рассмотрим в комплексной полуплоскости Re z ^ а, где постоянная а G (0, t0), функцию комплексной переменной

ж

Ф(z) = ^2фke-^* sinpkх0. (2.7)

k=i

Так как при Re z ^ а

2 2 2 2 \фк sinnkx0e-^a z| ^ Ce-»ka (2.8)

здесь C = const > 0, то в этой полуплоскости ряд, стоящий в правой части соотношения (2.8), сходится равномерно. Учитывая то, что каждый член этого ряда является аналитической функцией при Re z ^ а, и применяя теорему Вейерштрас-са [5, с. 68], получаем, что функция Ф^) является аналитической при Re z ^ а.

Поскольку в силу (2.7) Ф(-г) =0 на отрезке [ío, ti] действительной оси t из области аналитичности Ф(-г), но на основании теоремы единственности для аналитических функций следует Ф(-г) = 0 при Re z ^ а. Отсюда следует, что равенство (2.7) выполнено при всех t ^ to. Умножив равенство (2.7) на e(^ia) 1 и в полученном равенстве переходя к пределу при t ^ найдем

sin пХ0 / ф(£) sin = 0.

o

Последовательно повторяя аналогичные действия, получим

Í1

o

sin пкх0 <^>(£)sin ¡k= 0, к £ N. o

o

Отсюда, поскольку sin пкХо = 0, при всех к следует, что

Í ф(С) sin = 0, к = 1, 2,.... (2.9)

o

Из равенства (2.10), в силу полноты системы функций {sin¡nx}k'^i в пространстве L2 [0, l], следует, что ф(х) = 0 почти всюду на [0,1]. ■

Тогда из теорем 2 и 1 следует единственность решения обратной задачи 1. Пусть при некоторых Хо и к = p £ N нарушено условие теоремы 2, т. е.

sin прх0 = 0. (2.10)

Тогда обратная задача (1.2)-(1.6) при h(t) =0 и F(x,t) = 0 имеет ненулевое решение

up(x, t) = sin ¡pxe~(^pa) *, фр(х) = up(x, 0) = sin ¡¡px. (2.11)

Из уравнения (2.10) найдем значения

Хо n

Х0 = -o = -, n £ N, n<p, (2.12)

l p

при которых нарушается условие теоремы 2, т. е. нарушается единственность решения задачи 1.

Следовательно, установлен следующий критерий единственности решения обратной задачи (1.2)-(1.6).

Теорема 3. Условие sin пкх0 = 0 при всех к £ N необходимо и достаточно для единственности решения обратной задачи (1.2)-(1.6).

Из (2.12) следует, что когда Х принимает рациональные значения вида n, n < p, n £ N, p £ N, то нарушается единственность решения обратной задачи. При остальных значениях Х0 из (0,1), например, когда Х0 принимает иррациональные значения, условие теоремы 3 выполнено при всех к £ N, следовательно, для таких Х0 обратная задача может иметь не более одного решения.

3. Обратная задача 2

При условии существования функции д(г) решение задачи (1.2)—(1.5) определяется формулой (2.2). Полагая здесь х = хо, поменяв местами порядок интегрирования и суммирования, получим для искомой функции д(Ь) интегральное уравнение Вольтерра первого рода:

[ к(г,т)д(т)йт = н(г), о < г < т, (3.1)

о

с ядром

K(t, т) = J2 fne-(Mnа) (t-T) sin MnXQ (3.2)

= 1

и правой частью

то

X—^ 2 2t

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

h(t) = h(t) - 2_фке-*ка t sin Hkxq.

k=i

Теорема 4. Пусть ф(х) G C3[0,/], ф(0) = ф(1) = 0, ф"(0) = ф'(1) = 0, f (х) G C3[0,/], f (0) = f (/) = 0, f ''(0) = f''(/) = 0 и sinnnXQ = 0 при всех n G N, a;Q = x// G (0,1). Тогда, если f (xq) =0 и h(t) G C 1[0,T], h(0) = ф(х0), то уравнение (2.22) имеет единственное решение g(t) в классе функций C[0,Т].

Доказательство. Прежде всего найдем скорость убывания коэффициентов fn при n ^ то. В интеграле формулы (2.3), интегрируя по частям три раза, получим

1 2 Гl fП3)

fn = -— т / f "'(Ocos Hnd = - fV, (3.3)

Hn L J Q Hn

причем в силу неравенства Бесселя ряд из квадратов

(3)

1п ^ 0 при п ^ ж. Отсюда будем иметь

f(3) 2

сходится, поэтому

с

\fn\ < -3, (3.4)

¡П

где £n > 0 и £n ^ 0 при п ^ то.

В силу оценки (3.4) ряд (3.2) сходится равномерно и допускает почленное дифференцирование по t при 0 ^ т ^ t ^ T. Поэтому функция K't(t,r) непрерывна на указанном множестве. Дифференцируя уравнение (3.1) по t, имеем

K(t,t)g(t)+ ¡ Kt(t, т)д(т) = h'(t), 0 < t < T. (3.5)

J 0

Положив в (3.2) т = t, получим

K(t,t) = ^2 fn sin ¡nxo. (3.6)

n=1

Правая часть равенства (3.6) представляет собой разложение в ряд функции f (x) по системе sin¡nx в точке x = xo. По условию K(t,t) = f (xo) = 0. Поэтому уравнение (3.5) представляет собой интегральное уравнение Вольтерра второго рода с непрерывным ядром и непрерывной правой частью, следовательно, уравнение (3.1) имеет единственное решение g(t) G C[0, T]. I Теперь покажем, что условие

sin пих0 =0, Уп G N (3.7)

является существенным. Пусть для некоторых п = m и Xo G [0, /] выражение sin nmxo = 0. Тогда для функции f (x) = sin ¡mx и h(t) G C[0, T] существует ненулевое решение обратной задачи 2 (где y>(x) = 0)

Г t 2 2

w(x, t) = sin ¡mx g(s)e-^2ma2(t-s)ds. (3.8)

Jo

Как видим, в этом примере нарушается и условие f(xo) = 0 теоремы 4. Возникает вопрос о том, что если f (xo) = 0, то существует ли единственное решение обратной задачи (1.2)-(1.5), (1.7)?

n

Теорема 5. Пусть ф(х) £ C5[0,l], ф(0) = ф(1) = 0, ф'(0) = ф'(1) = 0, ф(4)(0) = = ф(4)(1) = 0, f(x) £ C5[0,l], f(0) = f(l) = 0, f''(0) = f'(l) = 0, f(4)(0) = f(4)(l) = 0 и sinnnxo = 0 при всех n £ N, xo = x/l £ (0,1). Тогда, если f(xo) = 0, f''(xo) = = 0 и h(t) £ C2[0,T], h'(0) = о?ф''(xo), то уравнение (3.1) имеет единственное решение g(t) в классе функций C[0, T].

Доказательство. Аналогично доказательству теоремы 4 имеем

J*(3) 1 О f l Á5)

fn = - íf3- = 41 f (5)(0sin Hnd = fnr • Hn Hn l J 0 Hn

Отсюда получим

с

\fn\ < -5• (3.9)

иП

В силу оценки (3.9) ряд (3.2) сходится равномерно и допускает почленное дифференцирование по г дважды при 0 ^ т ^ г ^ Т, поэтому функция К'ь' (г,т) непрерывна на указанном множестве.

Из условия ](хо) = 0 следует, что К(г,г) = 0, тогда уравнение (3.5) примет

вид г

/ Кг(г,т)д(т)в,т = Н'(г), 0 < г < Т. (3.10)

о

Дифференцируя уравнение (3.10) по г, имеем

К(г,г)д(г) + ГЩ(г,т)д(т)в.т = к"(г), 0 <г < Т, (3.11)

где

то

к'(t t) = -\^

l J 0

то 2 íl

K'(t,t) = -У~] -¡ f (C)sin HnCd^(Pno)2 sin Hnxo. (3.12)

n=1 l Jo

Проинтегрировав два раза по частям интеграл в ряде (3.12), получим

то 2 I' i

K't(t, t) = o2V - f''(£)sin sin Hnxo =

n=1 l Jo

2a1 то fi

—¡~ У2 f' '(O sin Hn^d^ sin Hnxo = a2 f'(xo). (3.13)

l n=1Jo

Дифференцируя ряд (3.5) дважды по t, найдем

то 2¡1 Гl

K'(t,T) = -a4 Y -^ f"(0sin Hn^e-(^a)2(t-T) sin Hnxo• (3.14)

n=1 l o

Теперь, интегрируя по частям два раза в интеграле ряда (3.14) и положив т = t, имеем

то — С i

K'(t,t) = a4 Y, у/ f(4\Osin Hnd cos Hnxo = a4f(4)(xo). (3.15)

n=1 l Jo

Ряды (3.14) и (3.15) сходятся равномерно при 0 ^ т ^ t ^ T в силу оценки (3.9).

Правая часть равенства (3.13) представляет собой разложение в ряд по системе sinHnx функции a2f'' (x) в точке x = xo. Следовательно, K'(t,t) = a2f '(xo) = 0. Уравнение (3.11) представляет собой интегральное уравнение Вольтерра второго рода с непрерывным ядром и непрерывной правой частью, а значит, оно имеет единственное решение g(t) в классе непрерывных на [0, T] функций. ■

o

На основании доказанных выше теорем можно показать справедливость следующего утверждения.

Теорема 6. Пусть ф(х) е C2k+3[0,l], ф(0) = ф(1) = 0, ф"(0) = ф''(1) = 0, ф(4)(0) = = ф(4)(1) = 0,...,ф(2к+2)(0) = ф(2к+2)(1) = 0, f(х) е C2k+3[0,l], f(0) = f(l) = 0,

f ''(0) = f '(l) = 0, f (4)(0) = f (4)(l) = 0,..., f (2k+2)(0) = f (2k+2)(l) =0 и sinnnx0 = 0 при всех n е N, x0 = x/l е (0,1). Тогда, если f (xo) = f ''(x0) = ... = f 2k-2(x0) = 0, f (2k)(xo) =0 и h(t) е Ck+1[0,T], h(k)(0) = a2k p(2k)(xo), где к е No, то уравнение (3.1) имеет единственное решение g(t) в классе функций C[0, T].

Доказательство. Аналогично доказанным выше теоремам найдем представления для коэффициентов

f(3) f(5) f(2k+1)

г _ jn _ Jn _ _ f i\k '" Jn ---"IT - —T - " ' - (-1)

t/3 1/5 ^fV ^^ X

Гп Гп ^n

J (2k+3) k+1 Jn

2k+3 '

f(2k+3)(e)cosrnd -(-1) 1гп jo Г

Отсюда получим оценку

\fn\ < -fc • (3-i6)

гп

В силу оценки (3.16) ядро K(t,T) уравнения (3.1) имеет непрерывную производную K(k+1) (t,T) при 0 < t < t < T. Дифференцируя уравнение (3.1) по t k+1 раз, с учетом условий f (xo) - f ''(xo) - • • • - f (2k-2)(xo) - 0 имеем

Г1

K(k)(t,t)g(t)+ / K(k+1)(t,T)g(T)- h(k+1)(t), 0 < t < T, (3.17)

Jo

где

2 1

K(k)(t,t)- -a2krlV] 2 f(2k-2)(e)sin Mn£d£ sin ^xo-

n=1 l ^

Проинтегрировав в последнем интеграле по частям дважды, получим

2 Г1

K(k)(t,t)- a2kV - f (2k)(C)sin rnîdÇ sin rnxo -

n=1 l Jo

-a2k œ Г1

-"^E/ f(2k)(^)sin MnCde sin rnxo - a2k f(2k)(xo)- (3.18)

l n=1Jo

Правая часть равенства (3.18) представляет собой разложение в ряд функции a2kf(2k)(x) в точке x - xo. По условию K(k)(t,t) - a2kf(2k)(xo) - 0. Таким образом, уравнение (3.17) представляет собой интегральное уравнение Вольтерра второго рода такого же типа, что и уравнения (3.5) и (3.11), а значит, оно имеет единственное решение в классе функций C[0, T]. I

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Отметим, что в теоремах 5 и 6 условие sin nnxo - 0 при всех n G N существенно, так как в противном случае обратная задача 2 имеет ненулевое решение (3.8) при f (x) - sinrmx, <£>(x) = 0 и любой функции h(t) G C[0,T].

4. Обратная задача 3 при g(t) = 1

Пусть u(x,t) и f (x) — решение задачи (1.2)-(1.5), (1.8). Следуя работам [6-8], рассмотрим интегралы

2 fl

Uk(t) - 2 u(x,t)sinrkxdx, (4.1)

l Jo

2 íl

fk = i j f (ж) sin Hkxdx, (4.2)

l J 0

где Hk = пк/l, к £ N.

На основании (4.1) введем вспомогательную функцию

А-е

2 [ е

ик,е (t) = i J u(x,t)sin Hkxdx, (4.3)

где £ > 0 — достаточно малое число.

Дифференцируя равенство (4.3) один раз по t и учитывая уравнение (1.3), получим

2 í'l—е 2 í'l—е

u'k e(t) = — ut(x,t)sin Hk xdx =— a2 uxx (x,t)sin Hkxdx + fk. (4.4)

l е l е

Проинтегрируем дважды по частям интеграл в правой части равенства (4.4) и, переходя в полученном равенстве к пределу при £ ^ 0 с учетом граничных условий (1.4), заключаем, что uk (t) удовлетворяет дифференциальному уравнению

uk (t) + (Hk a)2uk (t) = fk. (4.5)

Общее решение уравнения (4.5) при к £ N определяется по формуле

uk(t) = f + Ck e—[an )2t, (4.6)

где Ck — произвольные постоянные.

Для нахождения коэффициентов Ck и fk воспользуемся граничными условиями (1.5), (1.8) и формулой (4.1):

2 f1 2 f1 uk (0) = у u(x, 0) sin Hkxdx = — <£>(x) sin Hk xdx = фk, (4.7)

l 0 l 0

2 fl 2 ¡'l uk (to) = 7 u(x,to)sin Hk xdx = — »(x)sin Hk xdx = фk. (4.8)

l 0 l 0

Теперь, удовлетворив функцию (4.6) условиям (4.7) и (4.8), получим систему относительно fk и Ck :

4 + Ck = Фk, 4 + Ck e—»k10 = фk. (4.9)

Hk Hk

Из системы (4.9) найдем

Ck = • (4.10)

fk = HkФk - Hk ф —»1 , к £ N. (4.11)

1 — e »kt0

Подставив (4.10), (4.11) в (4.6), найдем окончательный вид функций

ЛЛ Ф — Фk)(1 - (412) uk(t) = фk +--7-27-, 0 < t < T, к £ N. (4.12)

1 — e »kt0

Пусть теперь ф^) = 0 и »(x) = 0 на [0,1]. Тогда фk = фk =0 и из (4.11), (4.12) следует, что uk (t) = 0 на сегменте [0, T] и fk =0 при всех к £ N. Отсюда в силу (4.1) и (4.2) имеем

fl fl

/ u(x,t)sin Hkxdx = 0, / f (x) sin Hk xdx = 0, к £ N.

00

В силу полноты системы {sin¡kx}+W в пространстве L2[0, l] из последних равенств следует, что u(x,t) =0 и f (x) = 0 почти всюду на [0,1] при любом t £ [0, Т]. Поскольку в силу (1.2) функции u(x,t) и f (x) непрерывны соответственно на D и (0, l), то u(x, t) = 0 в D и f (x) = 0 на [0, l].

Таким образом, доказано следующее утверждение. Теорема 7. Если существует решение задачи (1.2)-(1.5), (1.8), то оно единственно.

Теперь при определенных условиях на функции ф(x) и ф(x) покажем, что функции

u(x, t) = ик (t) sin ¡kx, (4.13)

k=i

w

f (x) = fk sin ¡kx (4.14)

k=i

удовлетворяют условию (1.2), где uk(t) и fk — определяются формулами (4.11), (4.12).

Лемма 1. При любых к £ N и t £ [0, Т] справедливы оценки:

u(t)| < Ki(ф| + |); \fk| < К2к2(ф| + |), к £ N, u'k(t)| < К3к2(1фк| + фk|), к £ N, где Ki — здесь и далее положительные постоянные, не зависящие от ф(x) и ф(x).

Справедливость данных оценок непосредственно следует из формул (4.11), (4.12).

Лемма 2. Если функции ф(x), ф(x) £ С3[0, l], ф''(0) = ф''(1) = ф''(0) = ф''(1) = 0, то справедливы равенства

Ф(3) ф(3)

фи = , Фи = , (4.15)

¡k ¡k

где фk, фk — коэффициенты разложения <^>(3)(x) и ф^^) в ряд по системе функ-

i+w lk=Q,

ций {cos ¡kx}+=°Q, при этом справедливы оценки

ефз))2 < 2 ||ф(з)|ц2[о,г]^ (Фkз))2 < 2 ^(3)nw (4.16) k=i l k=i l

Формально из (4.13) почленным дифференцированием составим ряды:

w

uxx(x,t) = ^J(-¡2k )uk (t)sin ¡kx, (4.17)

k=i

w

ut(t) = ^^ u'k (t) sin ¡kx. (4.18)

k=i

Ряды (4.13), (4.14), (4.17) и (4.18) при любом (x,t) £ D в силу лемм 1 и 2 мажорируются сходящимся числовым рядом

w

К4£ ^^ + ф^^)- (4.19)

k=i к

Тогда ряды (4.13), (4.14), (4.17) и (4.18) на основании признака Вейерштрасса сходятся абсолютно и равномерно на D. Тогда функции u(x,t) и f (x), определенные рядами (4.13) и (4.14), удовлетворяют условию (1.2).

Следовательно, нами доказана

Теорема 8. Если функции ф(х) и Ф(х) удовлетворяют условиям леммы 2, то существует единственное решение задачи (1.2)-(1.5), (1.8), и оно определяется рядами (4.13) и (4.14).

Теорема 9. Для решения (4.13) и (4.14) задачи (1.2)-(1.5), (1.8) справедливы оценки:

\\п(х,г)\\Ь2[0А < К( (\\ф\\ь2[0,1] + \\ф\\ь2[0,1]) , (4.20)

(х)\\ь2[0,1] < К (\\ф\\^22[о,,] + \\'Ф\\ш2[0,1]) , (4.21)

ЫхМоф) < К (\\ф\\С1[о,1] + \\Ф\\с1[0,1]) , (4.22)

\\1 (х)\\о2 [0,1] < К (\\ф\\с3[0,г] + \\Ф\\С3[0,«0 • (4.23)

Доказательство. Поскольку система Хк(х) ортогональна в Ь2[0,1], то из формулы (4.13) на основании леммы 1 получим

1 (1 \

\Кх,*)\\|2[о,г] = 2 Е^ Ю < 2К? 2 Е(^1 + Ф2) < к=1 V к=1 )

( 1 1 \ ) < 2К? 2+ 2ЕФ2 < К2 (\\ф\\=2М + НФНич) •

V к=1 к=1 )

Отсюда следует справедливость оценки (4.20). Аналогично из формулы (4.14) получим

1 (1 \

(х)\\2Ь2[0,1] = 2Е/к2 < 2К2 -^(фк + Фк) • (4.24)

к=1 V 2=1 )

По условию фк и фк можно представить в следующем виде:

ф(2> Ф(2)

Фк = -ФкТ~, Фк = -ф22г, (4.25)

¡к ¡к

где

2 {'I 2 {'I

Ф(2) = — ф"(х)^тНкхг1х, ф(2) = — ф''(х)вт¡кхдх.

1 ио 1 ио

Ю 1 J0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Тогда из оценки (4.24) имеем

(х)\\1т^ К (I)4 - £ (ф(2) + Ф (2>) = 2К| (4)4 X

* (2 £ Ф2)+ф(2^ < к (\\ф'\\1а10Л+\\Ф"\\1т) <

< К6 (\\ф\\ж|[о,г] + \\Ф\\ж22[0,г]

Из оценки (4.26) непосредственно следует оценка (4.21).

Пусть (х,Ю) произвольная точка Б. Тогда на основании леммы 1 имеем

/ +ТО \

\п(х,г)\ ^ \ик(Ю)\ < К1 ]Т(\фк\ + \Фк\) • (4.27)

к=1 к=1

Поскольку

фк = , Фк = , (4.28)

¡ к ¡ к

где

2 fi 2 fi

ф^ = 2 ф'(x)cOS Hkxdx, ф^Р = 2 ф'(х) cos Hkxdx,

l J 0 l J 0

то из соотношений (4.27) и (4.28) будем иметь

\u(x,t)\ < п f I (ф ^+\ф (1}о) <

<кl (f1 г й^Г+Ш'^У

= к^3 (\Ф{1Чх)\\Ь2т + \\ф(1Чх)\\Ь2т) .

Отсюда с учетом неравенств:

\ \ ф\\ь2[0,1] \ф\\с[0,г], \ \ ф\\ь2[0,1] \ф\\с[о,1],

\ \ ф' \ \ Ь2[0,1] \ ф' \ \ С[0,1], \ \ф' \ \ Ь2[0,1] < V/\ \ ф' \ \ с[0,1]

получим

\ и(х,г)\ < к (\ \ ф'\\С[0,г] + \\ ф'\\С[0,1]) = к (\ \ ф \\ С1[0,1] + \\ ф\\С1[0,ч)

из которой уже следует (4.22).

Аналогично из формулы (4.14) на основании леммы 1 имеем

\ f (x) \ = £ \ fk К K2J2 k2( \ Ф к \ + \ фк \).

'Л:

к=1 к=1

Из последнего соотношения с учетом (4.25) и (4.15) получим

I 1

/ 7 \ 3 1

\ /(х) \ < К2( Л Е 1 ( \ ф(3) \ + \ ф (3) \ ) <

^ ' к = 1

< К2((\\ ф(3) \\ Ы0,1] + \\ф(3) \\ ь2[0,1]) <

< К (\ \ ф(3) \\ с[0,г] + \\ ф(3) \\ С[0,ч) < к ( \\ ф \\ С3[0,,] + \\ ф\\С3[0,1]) . (4.29)

Из (4.29) следует справедливость оценки (4.23).■

5. Обратная задача 3 при д(1) ф 1

Рассуждая аналогично п. 4, введем функцию (4.1) и для нее получим уравнение

ик (ь) + (цка)2п к (г) = !кд(ь). (5.1)

Общее решение уравнения (5.1) при к € N определяется по формуле

и к (г) = Скв-(^ а)Ч + ¡ид к (г), (5.2)

здесь

С1

дк (г)= д(в)в-(^к a)2(t-s)ds. (5.3) 0

Удовлетворяя функцию (5.2) граничным условиям (4.7) и (4.8), найдем неизвестные постоянные С и Д:

Ck — фк, fk

при условии, когда при всех к £ N

1

Як (to)

фк - Фке-(1Лка)Но

Як (to)—0.

Подставляя (5.4) в (5.1), построим в явном виде функции

Uk

(t) — фкв-(^а)2г +

Як(t) Як (to)

фк - фкe-{»kа)Н

(5.4)

(5.5)

(5.6)

Теперь аналогично выше исходя из равенств (5.6), (4.1) и (4.2) на основании полноты системы {sin Икх}к^о доказывается единственность решения обратной задачи 3 при произвольной непрерывной функции g(t) и выполнении условий (5.5) при всех к £ N.

Если при некоторых t0 и к — p выражение gp(t0) — 0, то однородная задача 3 (где ф(х) — ф(х) = 0) имеет ненулевое решение

u(x,t) — fpgp(t)sin црх, f (x) — fp sin ppx,

здесь fp — 0 — произвольная постоянная.

Следовательно, нами установлен критерий единственности решения задачи 3. Теорема 10. Если существует решение задачи 3, то оно единственно только тогда, когда при всех к £ N выполнены условия (5.5).

Лемма 3. Если g(y) непрерывна на [0, T] и \g(y)\ ^ т — const > 0, то существует постоянная C0 > 0 такая, что при всех к £ N

\Як (to^ > <C2 ■

Доказательство. На основании теоремы о среднем из (5.3) имеем

(5.7)

Як (to) — я(0 а) (to-s)ds — я(0

Отсюда получим оценку снизу

1 — e-(na)2to (Ик a)2 ''

е £ [0,to]■

\Як(to)\ >{—) m(1 - e-(»ia)2to)к-2, \na J

из которой уже следует (5.7).■

Решение в этом случае строится в виде суммы рядов (4.13) и (4.14), где только коэффициенты (Ь) и Д определяются формулами (5.4) и (5.6). Теорема 11. Если функции ф(х) и ф(х) удовлетворяют условиям леммы 2, а функция д(у) условиям леммы 3, то существует единственное решение задачи 3, которое определяется рядами (4.13) и (4.14), где коэффициенты находятся по формулам (5.4) и (5.6).

Отметим, что для решения (4.13) и (4.14) задачи 3 при д(Ь) ф 1 справедливы оценки (4.20)-(4.23), установленные в теореме 9, но только с другими постоянными г = 5, 8.

o

2

Литература

[1] Денисов А.М. Введение в теорию обратных задач. М.: Изд-во МГУ, 1994, 208 с.

[2] Кабанихин С.И. Обратные и некорректные задачи. Новосибирск: Сибирское научное издательство, 2009, 457 с.

[3] Тихонов А.Н., Самарский А.Н. Уравнения математической физики. 3-е изд. М.: Физматлит, 1966. 724 с.

[4] Сабитов К.Б. Уравнения математической физики. М.: Физматлит, 2013, 352 с.

[5] Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. М.: Наука, 1973, 749 с.

[6] Сабитов К.Б.,Сафин Э.М. Обратная задача для уравнения параболо-гиперболиче-ского типа в прямоугольной области // ДАН. 2009. Т. 429. № 4. С. 451-454.

[7] Сабитов К.Б., Сафин Э.М. Обратная задача для уравнения смешанного параболо-гиперболического типа в прямоугольной области // Известия вузов. Сер.: Математика. 2010. № 4. С. 55-62.

[8] Сабитов К.Б., Сафин Э.М. Обратная задача для уравнения смешанного параболо-гиперболического типа // Мат. заметки. 2010. Т. 87. № 6. С. 907-918.

References

[1] Denisov A.M. Introduction to the theory of inverse problems. M., Izd-vo MGU, 1994, 208 p. [in Russian].

[2] Kabanikhin S.I. Inverse and ill-posed problems. Novosibirsk, Sibirskoe nauchnoe izdatelstvo, 2009, 457 p. [in Russian].

[3] Tikhonov A.N., Samaraskii A.N. Equations of mathematical physics. M., Izd-vo Fizmatlit, 1966, 724 p. (3rd edition) [in Russian].

[4] Sabitov K.B. Equations of mathematical physics. M., FIZMATLIT, 2013, 352 p. [in Russian].

[5] Lavrentev M.A., Shabat B.V. Methods of the theory of complex variable. M., Nauka, 1973, 749 p. [in Russian].

[6] Sabitov K.B., Safin E.M. Inverse problem for a parabolic-hyperbolic type for rectangular area. DAN [Proceedings of the Russian Academy of Sciences], 2009, Vol. 429, no. 4, pp. 451-454 [in Russian].

[7] Sabitov K.B., Safin E.M. Inverse problem for a mixed parabolic-hyperbolic type for a rectangular area. Izvestiia vuzov. Matematika [News of Higher Educational Institutions. Mathematics], 2010, no. 4, pp. 55-62 [in Russian].

[8] Sabitov K.B., Safin E.M. Inverse problem for a mixed parabolic-hyperbolic type. Mat. zametki [Mathematical notes], 2010, Vol. 87, no. 6, pp. 907-918 [in Russian].

A.R. Zaynullov2

INVERSE PROBLEMS FOR THE HEAT EQUATION

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

The inverse problem of finding initial conditions and the right-hand side had been studied for the inhomogeneous heat equation on the basis of formulas for the solution of the first initial-boundary value problem. A criterion of uniqueness of solution of the inverse problem for finding the initial condition was found with Spectral analysis. The right side of the heat equation is represented as a product of two functions, one of which depends on the spatial coordinates and the other from time. In one task, along with an unknown solution is sought factor on the right side, depending on the time, and in another — a factor that depends on the spatial coordinates. For these tasks, we prove uniqueness theorems, the existence and stability of solution.

Key words: heat equation, first initial-boundary value problem, inverse problems, spectral method, integral equation, uniqueness, existence, stability.

Статья поступила в редакцию 20/777/2015. The article received 20/777/2015.

2 Zaynullov Artur Rashitovich (arturzayn@mail.ru), Department of Mathematical Analysis, Sterlitamak Branch of Bashkir State University, 49, Prospekt Lenina Street, Sterlitamak, 453103, Russian Federation.

i Надоели баннеры? Вы всегда можете отключить рекламу.