Журнал вычислительной математики и математической физики, 2020, T. 60, № 9, стр. 1604-1619
Об определении источников с компактными носителями для одномерного уравнения теплопроводности
В. В. Соловьёв *
Национальный исследовательский ядерный университет “МИФИ”
115409 Москва, Каширское ш. 31, Россия
* E-mail: soloviev.vyacheslav@gmail.com
Поступила в редакцию 10.05.2019
После доработки 03.02.2020
Принята к публикации 09.04.2020
Аннотация
Исследуется обратная задача определения источника в одномерном уравнении теплопроводности для случая первой краевой задачи. В качестве “переопределения” (дополнительной информации о решении прямой задачи) задан след решения прямой задачи на отрезках прямой внутри интервала в финальный момент времени. Доказана справедливость альтернативы Фредгольма для этой задачи и получены достаточные условия существования и единственности решения этой обратной задачи. Рассмотрение обратной задачи проводится в классах гладких функций, производные которых удовлетворяют условию Гёльдера. Библ. 15.
1. ПОСТАНОВКА ОБРАТНОЙ ЗАДАЧИ
Пусть $T > 0$, $0 < \alpha < 1$, $l > 0$ – фиксированные числа, обозначим интервал оси $x$ через $\Omega = (0,l)$, отрезок $\bar {\Omega } = [0,l]$, и на плоскости с декартовыми координатами точек $(x,t)$, определим прямоугольник с верхней крышкой ${{\Omega }_{T}} = \Omega \times (0,T]$ и параболической границей $\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\partial } {{\Omega }_{T}} = (\{ 0\} \times [0,T]) \cup (\{ l\} \times [0,T]) \cup ([0,l] \times \{ 0\} )$. В замкнутом прямоугольнике ${{\bar {\Omega }}_{T}} = \bar {\Omega } \times [0,T]$ рассмотрим первую начально-краевую задачу для уравнения теплопроводности с источником специального вида, а именно, рассмотрим задачу определения функции $u\,:{{\bar {\Omega }}_{T}} \to \Re $ из условий
(1)
$(Lu)(x,t) = {{u}_{t}}(x,t) - {{u}_{{xx}}}(x,t) = {{f}_{k}}(x){{h}^{k}}(x,t) + g(x,t) = q(x,t) + g(x,t),\quad (x,t) \in {{\Omega }_{T}},$(2)
$u(0,t) = {{\mu }_{1}}(t),\quad u(l,t) = {{\mu }_{2}}(t),\quad t \in [0,T],\quad u(x,0) = \varphi (x),\quad x \in \bar {\Omega }.$Так как правая часть уравнения (1) содержит функции, вообще говоря, имеющие разрыв непрерывности по $x$ I рода в точках $x \in A \cup B$, при фиксированном переменном $t \in (0,T]$, то естественным шагом является расширение области определения оператора $L$, т.е. использование для уравнения (1) некоторого обобщенного решения. Будем включать в область определения оператора $L$ функции $u(x,t)$, имеющие у второй производной по $x$ при фиксированном $t \in (0,T]$ разрывы I рода при $x \in A \cup B$. Именно таким образом и будет пониматься далее уравнение (1). При рассмотрении решений задачи (1), (2) и изучении их свойств будем использовать теорию решений уравнений параболического типа в пространствах Гёльдера, подробное изложение которой приведено в монографиях [1], [2]. Определим необходимые для дальнейших построений пространства функций (определение используемых в данной работе пространств Гёльдера и норм в них см. [1, с. 16])
Неформально говоря, для функции $u$ справедливо включение $u \in {{U}_{0}}({{\Omega }_{T}})$, если она непрерывна на всем замкнутом прямоугольнике ${{\bar {\Omega }}_{T}}$, в каждом из прямоугольников ${{P}_{T}}$ и прямоугольников ${{Q}_{T}}$ она имеет две непрерывные производные по $x$ и одну по $t,0 < t \leqslant T$, имеет непрерывную первую производную по $x$ во всем прямоугольнике ${{\Omega }_{T}}$. Для принадлежности функции $u$ пространству $U({{\Omega }_{T}})$ дополнительно требуется, чтобы в каждом из прямоугольников, составляющих множества ${{\widetilde P}_{T}},\;{{\widetilde Q}_{T}}$, вторая производная по $x$ и первая по $t$ удовлетворяли условиям Гёльдера по $x$ и $t$ с показателями $\alpha $ и $\alpha {\text{/}}2$ соответственно. Отметим здесь, что функциональные пространства ${{U}_{0}}({{\Omega }_{T}})$ и $U({{\Omega }_{T}})$ не являются общепринятыми при рассмотрении задачи (1), (2), которую в дальнейшем будем называть прямой задачей (см., впрочем, [1, с. 262]). Тем не менее, как будет видно из дальнейшего изложения, при изучении обратной задачи для уравнения (1) удобно рассматривать его решение именно в таких классах функций. Это даст возможность доказать для этого уравнения теорему существования и единственности этой обратной задачи в этих классах функций. Предварительно перед рассмотрением обратной задачи для уравнения (1) докажем существование и единственность решения прямой задачи (1), (2) в классах функций $u \in {{U}_{0}}({{\Omega }_{T}})$ и $U({{\Omega }_{T}})$, что, возможно, не лишено интереса и само по себе. Сначала докажем для задачи (1), (2) справедливость принципа максимума для случая $u \in {{U}_{0}}({{\Omega }_{T}})$ в форме так называемого “слабого принципа максимума”. Отметим здесь, что в классах функций $u \in {{U}_{0}}({{\Omega }_{T}})$ его надо специально доказывать, так как решение уравнения (1) в этом классе функций не является, вообще говоря, классическим решением. Всюду далее, как это обычно принято при рассмотрении различного вида обобщенных решений, будем писать $(x,t) \in {{\Omega }_{T}}$, понимая при этом, что, вообще говоря, при $x \in A \cup B$ левая и правая части уравнения (1) не определены, так как в этих точках вторая производная по $x$ может иметь разрывы непрерывности.
Лемма 1 (принцип максимума). Пусть функция $u \in {{U}_{0}}({{\Omega }_{T}})$ и удовлетворяет, в указанном выше смысле, неравенству $(Lu)(x,t) \leqslant 0$ ($(Lu)(x,t) \geqslant 0$), $(x,t) \in {{\Omega }_{T}}$. Тогда функция $u(x,t)$ принимает значение своего положительного максимума (отрицательного минимума) на параболической границе прямоугольника ${{\Omega }_{T}}$ (если, конечно, у этой функции они есть).
Доказательство будем строить от противного. Сначала рассмотрим случай $(Lu)(x,t) \leqslant 0$. Пусть значения положительного максимума есть где-то на ${{\overline \Omega }_{T}}$, но на параболической границе они не достигаются. Пусть
Следствие 1. Задача (1), (2) в классе функций ${{U}_{0}}({{\Omega }_{T}})$ не может иметь двух различных решений.
Следствие 2. Пусть существует решение задачи (1), (2) в классе функций ${{U}_{0}}({{\Omega }_{T}})$ и выполнены неравенства $\varphi (x) \geqslant 0$, $x \in \Omega $, ${{\mu }_{1}}(t) \geqslant 0$, ${{\mu }_{2}}(t) \geqslant 0$, $t \in [0,T]$, $g(x,t) \geqslant 0$, $q(x,t) \geqslant 0$, $(x,t) \in {{\Omega }_{T}}$. Тогда для функции $u$, являющейся решением задачи (1), (2), будет справедливо неравенство $u(x,t) \geqslant 0$, $(x,t) \in {{\bar {\Omega }}_{T}}$. Более того, если при этом хотя бы одна из этих функций отлична от нуля хотя бы в одной точке, то для функции $u$ будет выполнено строгое неравенство $u(x,T) > 0$, $(x,T) \in {{\Omega }_{T}}$. Кроме того, если известно, что ${{f}_{k}} \geqslant 0$, ${{h}^{k}} \geqslant 0$, $g \geqslant 0$, $\varphi \geqslant 0$, ${{\mu }_{1}} \geqslant 0$, ${{\mu }_{2}} \geqslant 0$ и хотя бы в одной точке $(x{\text{*}},t{\text{*}}) \in {{\Omega }_{T}}$ справедливо, что $u(x{\text{*}},t{\text{*}}) = 0$, то для любых точек $(x,t) \in {{\bar {\Omega }}_{T}}$, $t \leqslant t{\text{*}}$ будет выполнено равенство $u(x,t) = 0$.
Следствие 3. Пусть существует решение задачи (1), (2) в классе функций ${{U}_{0}}({{\Omega }_{T}})$. Тогда для функции $u$ справедлива априорная оценка
Доказательства этих утверждений проводятся стандартными способами так же, как для классических решений, используя принцип максимума (доказанная выше лемма 1) (см., например, [2, с. 59]).
Докажем теперь теорему существования и единственности решения задачи (1), (2) в определенных выше классах функций $U({{\Omega }_{T}})$, ${{U}_{0}}({{\Omega }_{T}})$.
Терема 1. Пусть справедливы включения $\varphi \in C(\bar {\Omega })$, ${{\mu }_{1}},{{\mu }_{1}} \in C[0,T]$, ${{h}^{k}},g \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, ${{f}_{k}} \in {{C}^{\alpha }}[{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, выполнены условия согласования ${{\mu }_{1}}(0) = \varphi (0)$, ${{\mu }_{2}}(0) = \varphi (l)$. Тогда задача (1), (2) имеет единственное решение в классе функций ${{U}_{0}}({{\Omega }_{T}})$. Если дополнительно к этим условиям справедливы включения $h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, то это решение из класса $U({{\Omega }_{T}})$.
Доказательство. Единственность решения следует из следствия 1 леммы 1. Докажем сначала существование решения задачи (1), (2) в классе функций ${{U}_{0}}({{\Omega }_{T}})$. Будем строить решение задачи (1), (2) в виде суммы двух функций
В качестве функции ${{u}_{0}}$ возьмем решение следующей начально-краевой задачи:(4)
${{u}_{0}}(0,t) = {{\mu }_{1}}(t),\quad {{u}_{0}}(l,t) = {{\mu }_{2}}(t),\quad t \in [0,T],\quad {{u}_{0}}(x,0) = \varphi (x),\quad x \in \bar {\Omega }.$Ясно, что $(0,l) = {{\widetilde Q}^{{(1)}}} \cup \overline P \cup {{\widetilde Q}^{{(2)}}}$, ${{\Omega }_{T}} = \widetilde Q_{T}^{{(1)}} \cup {{\widetilde P}_{T}} \cup \widetilde Q_{T}^{{(2)}}$. Пусть далее заданы функции $f{\kern 1pt} :\;\Re \to \Re $, $f(x) = 0$, $x \notin [a,b]$, $f \in {{C}^{\alpha }}[a,b]$, $h \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$. Тогда задача (5), (6) при $N = 1$ примет следующий вид:
Покажем, что задача (7), (8) имеет решение в классе функций ${{U}_{0}}({{\Omega }_{T}})$ или, подробнее, в принятых для случая $N = 1$ обозначениях
(9)
$w(x,t) = \int\limits_0^t {d\tau } \int\limits_0^l {G(x,\xi ;t,\tau )f(\xi )h(\xi ,\tau )d\xi } .$Покажем, что определенная фoрмулой (9) функция $w \in {{U}_{0}}({{\Omega }_{T}})$ и удовлетворяет всем требованиям к решению задачи (7), (8). По определению функции Грина задачи (7), (8) она имеет вид (см. [1, с. 464])
где функция $Z$ – фундаментальное решение уравнения (1), определяемое по формулеПредставив выражение функции $g(x,\xi ;t,\tau )$ через тепловые потенциалы двойного слоя с носителями при $x = 0$, $x = l$ (см. [1, с. 462], [4, с. 735]) получим, что эта функция бесконечно дифференцируема по $x,\;t$ при $0 < x < l$, $\tau < t \leqslant T$. Таким образом, формула (9) с учетом того, что носитель функции $f$ лежит на отрезке $[a,b]$, примет вид
(10)
$w(x,t) = \int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {Z(x - \xi ;t - \tau )f(\xi )h(\xi ,\tau )d\xi } - \int\limits_0^t {d\tau } \int\limits_0^l {g(x,\xi ;t,\tau )f(\xi )h(\xi ,\tau )d\xi } .$Пользуясь свойствами тепловых потенциалов и абсолютной интегрируемостью фундаментального решения, из формулы (10) получаем, что для функции $w$ справедливы условия (8). В формуле (10) второе слагаемое не является несобственным интегралом, удовлетворяет уравнению теплопроводности и, в силу бесконечной дифференцируемости функции $g(x,\xi ;t,\tau )$ в ${{\Omega }_{T}}$, лежит в пространстве $C({{\bar {\Omega }}_{T}}) \cap C_{{x,t}}^{{2,1}}({{\Omega }_{T}})$. Более того, в силу этих же причин оно есть элемент пространства ${{C}^{{2 + \alpha ,1 + \alpha /2}}}({{\Omega }_{T}})$. Для изучения первого слагаемого воспользуемся следующими известными представлениями его производных (см. [1, с. 306], [2, с. 24]):
(11)
$\frac{\partial }{{\partial x}}\int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {Z(x - \xi ;t - \tau )f(\xi )h(\xi ,\tau )d\xi } = \int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {\frac{{\partial Z}}{{\partial x}}(x - \xi ;t - \tau )f(\xi )h(\xi ,\tau )d\xi } ,\quad 0 < t \leqslant T,\quad x \in \Re ,$(12)
$\begin{gathered} \frac{{{{\partial }^{2}}}}{{\partial {{x}^{2}}}}\int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {Z(x - \xi ;t - \tau )f(\xi )h(\xi ,\tau )d\xi } = \int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {\frac{{{{\partial }^{2}}Z}}{{\partial {{x}^{2}}}}(x - \xi ;t - \tau )[f(\xi )h(\xi ,\tau )} - \\ \, - f(x)h(x,\tau )]d\xi ,\quad 0 < t \leqslant T,\quad x \in \Re ,\quad x \ne a,\quad x \ne b, \\ \end{gathered} $(13)
$\begin{gathered} \frac{\partial }{{\partial t}}\int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {Z(x - \xi ;t - \tau )f(\xi )h(\xi ,\tau )d\xi } = \int\limits_0^t {d\tau } \int\limits_{ - \infty }^\infty {\frac{{{{\partial }^{2}}Z}}{{\partial {{x}^{2}}}}(x - \xi ;t - \tau )[f(\xi )h(\xi ,\tau )} - \\ \, - f(x)h(x,\tau )]d\xi + f(x)h(x,t),\quad 0 < t \leqslant T,\quad x \in \Re ,\quad x \ne a,\quad x \ne b. \\ \end{gathered} $Каждый из интегралов (11)–(13) является несобственным только при $a < x < b$. В силу известных оценок фундаментального решения и его производных (см. [2, с. 19]) при указанных ниже постоянных $\mu $ справедливы оценки
(14)
$\left| {{{Z}_{x}}(x - \xi ,t - \tau )} \right| \leqslant {{C}_{2}}{{(t - \tau )}^{{ - \mu }}}{\text{|}}x - \xi {{{\text{|}}}^{{2\mu - 1}}},\quad 1{\text{/}}2 < \mu < 1,$Из этих оценок следует равномерная сходимость интегралов в формулах (10), (11) при $0 < x < l$, в формулах (12), (13) при $a < x < b$ (при $0 < x < a$, $b < x < l$ в формулах (12), (13) интегралы не являются несобственными). Таким образом, показали, что функция $w(x,t)$, определенная по формуле (9), является непрерывной на ${{\bar {\Omega }}_{T}}$, ее производная по $x$, функция ${{u}_{x}}(x,t)$, непрерывна на ${{\Omega }_{T}}$, производные ${{u}_{t}}(x,t),{{u}_{{xx}}}(x,t)$ непрерывны на прямоугольниках $Q_{T}^{{(1)}},\;{{P}_{T}},\;Q_{T}^{{(2)}}$. В силу указанной выше равномерной сходимости интегралов и известных свойств фундаментального решения $Z$, непосредственной проверкой устанавливается, что функция $w$ удовлетворяет уравнению (7) при $(x,t) \in {{\Omega }_{T}}$, $x \ne a,b$. То есть доказано существование решения задачи (7), (8) в классе функций ${{U}_{0}}({{\Omega }_{T}})$, представимое формулой (9). Покажем далее, что если дополнительно известно, что ${{h}_{t}} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, то эта функция принадлежит пространству $U({{\Omega }_{T}})$, т.е. справедливо включение $w \in {{U}_{0}}({{\Omega }_{T}}) \cap {{C}^{{2 + \alpha ,1 + \alpha /2}}}({{\widetilde P}_{T}}) \cap {{C}^{{2 + \alpha ,1 + \alpha /2}}}({{\widetilde Q}_{T}})$. Для этого рассмотрим функцию $v\,:{{\bar {\Omega }}_{T}} \to \Re $, определенную по правилу
(15)
$v(x,t) = \int\limits_0^t {d\tau } \int\limits_0^l {G(x,\xi ;t,\tau )f(\xi ){{h}_{\tau }}(\xi ,\tau )d\xi } + \int\limits_0^l {G(x,\xi ;t,0)f(\xi )h(\xi ,0)d\xi } = {{v}_{1}}(x,t) + {{v}_{2}}(x,t),\quad (x,t) \in {{\bar {\Omega }}_{T}}.$Определим пространства функций
(17)
${{v}_{2}}(0,t) = 0,\quad {{v}_{2}}(l,t) = 0,\quad t \in [0,T],\quad {{v}_{2}}(x,0) = f(x)h(x,0),\quad x \in \bar {\Omega },\quad x \ne a,\quad x \ne b.$(19)
$v(0,t) = 0,\quad v(l,t) = 0,\quad t \in [0,T],\quad v(x,0) = f(x)h(x,0),\quad x \in \bar {\Omega },\quad x \ne a,\quad x \ne b.$(21)
${{w}^{{(n)}}}(0,t) = 0,\quad {{w}^{{(n)}}}(l,t) = 0,\quad t \in [0,T],\quad {{w}^{{(n)}}}(x,0) = 0,\quad x \in \bar {\Omega },$(23)
${{v}^{{(n)}}}(0,t) = 0,\quad {{v}^{{(n)}}}(l,t) = 0,\quad t \in [0,T],\quad {{v}^{{(n)}}}(x,0) = {{f}^{{(n)}}}(x)h(x,0),\quad x \in \bar {\Omega }.$(24)
${{w}^{{(n)}}}(x,t) = \int\limits_0^t {{{v}^{{(n)}}}(x,\tau )d\tau ,(x,t) \in {{{\bar {\Omega }}}_{T}}} .$(26)
$\begin{gathered} \left| {v({{x}_{2}},{{t}_{2}}) - v({{x}_{1}},{{t}_{1}})} \right| \leqslant C(\varepsilon ,\delta )({{\left| {{{x}_{2}} - {{x}_{1}}} \right|}^{\delta }} + {{\left| {{{t}_{2}} - {{t}_{1}}} \right|}^{{\delta /2}}})\left\| f \right\|\left( {\left\| {h(.,0)} \right\| + \left\| {{{h}_{t}}} \right\|} \right), \\ {{x}_{2}},{{x}_{1}} \in [a - \varepsilon ,b + \varepsilon ],\quad {{t}_{2}},{{t}_{1}} \geqslant \varepsilon > 0. \\ \end{gathered} $Для этого рассмотрим следующую последовательность начально-краевых задач для уравнения теплопроводности:
(28)
$v_{1}^{{(n)}}(0,t) = 0,\quad v_{1}^{{(n)}}(l,t) = 0,\quad t \in [0,T],\quad v_{1}^{{(n)}}(x,0) = 0,\quad x \in \bar {\Omega }.$(29)
$\left| {v_{1}^{{(n)}}({{x}_{2}},{{t}_{2}}) - v_{1}^{{(n)}}({{x}_{1}},{{t}_{1}})} \right| \leqslant C(\delta )({{\left| {{{x}_{2}} - {{x}_{1}}} \right|}^{\delta }} + {{\left| {{{t}_{2}} - {{t}_{1}}} \right|}^{{\delta /2}}})\left\| f \right\|\left\| {{{h}_{t}}} \right\|,\quad ({{x}_{2}},{{t}_{2}}),({{x}_{1}},{{t}_{1}}) \in {{\bar {\Omega }}_{T}}.$(30)
$\left| {{{v}_{1}}({{x}_{2}},{{t}_{2}}) - {{v}_{1}}({{x}_{1}},{{t}_{1}})} \right| \leqslant C(\delta )({{\left| {{{x}_{2}} - {{x}_{1}}} \right|}^{\delta }} + {{\left| {{{t}_{2}} - {{t}_{1}}} \right|}^{{\delta /2}}})\left\| f \right\|\left\| {{{h}_{t}}} \right\|,\quad ({{x}_{2}},{{t}_{2}}),({{x}_{1}},{{t}_{1}}) \in {{\bar {\Omega }}_{T}}.$Для произвольной функции $v{\kern 1pt} :\;{{\overline \Omega }_{T}} \to \Re $ и произвольного числа $0 < \delta < 1$ обозначим
(32)
$\left\| {v(.,T)} \right\|_{\omega }^{{(\delta )}} = \sup \{ \left| {v(x,T)} \right| + \left| {v({{x}_{2}},T) - v({{x}_{1}},T)} \right|{\text{/}}{{\left| {{{x}_{2}} - {{x}_{1}}} \right|}^{\delta }},\;{{x}_{2}} \ne {{x}_{1}},\;x,{{x}_{1}},{{x}_{2}} \in {{\overline P }^{{(k)}}},\;k = 1,...,N\} .$Следствие 1. Пусть $w(x,t)$ – решение следующей задачи в классе функций ${{U}_{0}}({{\Omega }_{T}}){\text{:}}$
(33)
$(Lw)(x,t) = {{w}_{t}}(x,t) - {{w}_{{xx}}}(x,t) = {{f}_{k}}(x){{h}^{k}}(x,t),\quad (x,t) \in {{\Omega }_{T}},$(35)
$\left\| {{{w}_{t}}(.,T)} \right\|{\text{|}}_{\omega }^{{(\delta )}} \leqslant C(\delta )\left\| {{{f}_{k}}} \right\|\left( {\left\| {h_{t}^{k}} \right\| + \left\| {{\text{|}}{{h}^{k}}(.,0)} \right\|} \right).$Доказательство. Оценка (35) следует из применения оценки (30) к каждому слагаемому в правой части уравнения (33) и их суммирования при фиксированном $\varepsilon = \min \{ {{a}_{1}}{\text{/}}2,(l - {{b}_{N}}){\text{/}}2,T{\text{/}}2\} > 0$. Следствие 1 доказано.
Для формулировки следующего следствия определим множество $\bar {\Omega }_{T}^{*}$ и соответствующее функциональное пространство для случая $N > 1$, аналогично случаю $N = 1$. Пусть
Следствие 2. Пусть $u(x,t)$ – решение задачи (33), (34) в классе функций ${{U}_{0}}({{\Omega }_{T}})$, выполнены условия ${{h}^{k}},h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, ${{f}_{k}} \in {{C}^{\alpha }}[{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$. Тогда функция $v = {{u}_{t}} \in V(\Omega _{T}^{*})$ удовлетворяет задаче
(36)
$(Lv)(x,t) = {{v}_{t}}(x,t) - {{v}_{{xx}}}(x,t) = {{f}_{k}}(x)h_{t}^{k}(x,t),\quad (x,t) \in {{\bar {\Omega }}_{T}},$(37)
$v(0,t) = 0,\quad v(l,t) = 0,\quad t \in [0,T],\quad v(x,0) = {{f}_{k}}(x){{h}^{k}}(x,0),\quad x \in \bar {\Omega },\quad x \notin A \cup B.$Перейдем к постановке обратной задачи для уравнения (1). Для этого будем предполагать, что в уравнении (1), кроме функции $u$, неизвестны еще функции ${{f}_{k}},k = 1,...,N$, но при этом задана дополнительная информация о функции $u$. Будем предполагать, что, кроме условий (2), известны также следы функции $u$ на отрезках ${{\bar {\omega }}^{{(k)}}},k = 1,...,N$ (такого рода информация при постановке обратной задачи называется обычно переопределением). Далее для удобства изложения будем обозначать $f = ({{f}_{1}},...,{{f}_{N}})$, $\chi = ({{\chi }_{1}},...,{{\chi }_{N}})$, $h(x,t) = ({{h}^{1}}(x,t),...,{{h}^{N}}(x,t))$. Соответствующие функциональные пространства (очевидно банаховы) обозначим
Таким образом, рассмотрим обратную задачу определения $N + 1$ функции $(u,f) \in U({{\Omega }_{T}}) \times {{F}^{\alpha }}$ из условий
(39)
$u(0,t) = {{\mu }_{1}}(t),\quad u(l,t) = {{\mu }_{2}}(t),\quad t \in [0,T],\quad u(x,0) = \varphi (x),\quad x \in \bar {\Omega },$Обратная задача для уравнения теплопроводности с источниками, имеющими компактный носитель, рассматривалась в работе [7] для случая $N = 2$ и независимости функций ${{h}^{k}}$ от $x$. В этой работе при другом переопределении изучался вопрос о единственности ее решения и были построены примеры неединственности решения обратной задачи. Другие постановки обратных задач определения источников в уравнении теплопроводности, близкие к задаче (38)–(40), а также историю изучения таких задач см. в работах [8]–[13]. В данной работе для обратной задачи (38)–(40) получены некоторые достаточные условия ее единственности (в том числе носящие глобальный характер), доказана альтернатива Фредгольма и получены некоторые достаточные условия ее однозначной разрешимости.
2. ЕДИНСТВЕННОСТЬ РЕШЕНИЯ ДЛЯ ОБРАТНОЙ ЗАДАЧИ (36)–(38)
Рассмотрим сначала важнейший для постановки любой обратной задачи вопрос о единственности ее решения. Сначала сформулируем и докажем теорему, дающую достаточные условия единственности решения обратной задачи (38)–(40), имеющие локальный характер.
Теорема 2. Пусть для функций $h$ выполнены условия ${{h}^{k}},h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, $\left| {{{h}^{k}}(x,T)} \right| \geqslant {{h}_{T}} > 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, где ${{h}_{T}} > 0$ – некоторая постоянная. Тогда существует число ${{\nu }_{0}} > 0$ такое, что для любых не пересекающихся отрезков ${{\bar {\omega }}_{k}} \subset {{\Omega }_{T}}$ таких, что $max\{ {{b}_{k}} - {{a}_{k}}\} < {{\nu }_{0}}$, обратная задача (38)–(40) не может иметь двух различных решений.
Доказательство. В силу линейности задачи (38)–(40) достаточно доказать, что имеет только тривиальное решение обратная задача определения функций $(u,f) \in U({{\Omega }_{T}}) \times {{F}^{\alpha }}$ из следующих условий:
Доказательство строим от противного. Пусть задача (41)–(43) имеет нетривиальное решение. Но тогда, в силу единственности решения прямой задачи (41), (42), обязательно хотя бы одна из функций ${{f}_{k}}$ отлична от тождественного нуля. Подставим в уравнение (41) $t = T$ и $x \in [{{a}_{k}},{{b}_{k}}]$. В этом случае из условий переопределения (43) получим для функции $u$ условия (по индексу ${{k}_{0}}$ здесь и далее суммирование не проводится)(44)
${{u}_{t}}(x,T) = {{f}_{{{{k}_{0}}}}}(x){{h}^{{{{k}_{0}}}}}(x,T),\quad x \in [{{a}_{{{{k}_{0}}}}},{{b}_{{{{k}_{0}}}}}],\quad {{k}_{0}} = 1,...,N.$(45)
$(Lv)(x,t) = {{v}_{t}}(x,t) - {{v}_{{xx}}}(x,t) = {{f}_{k}}(x)h_{t}^{k}(x,t),\quad (x,t) \in {{\Omega }_{T}},$(46)
$v(0,t) = 0,\quad v(l,t) = 0,\quad t \in [0,T],\quad v(x,0) = {{f}_{k}}(x){{h}^{k}}(x,0),\quad x \in \bar {\Omega },\quad x \notin A \cup B.$(47)
$v(x,T) = {{f}_{{{{k}_{0}}}}}(x){{h}^{{{{k}_{0}}}}}(x,T),\quad x \in [{{a}_{{{{k}_{0}}}}},{{b}_{{{{k}_{0}}}}}],\quad {{k}_{0}} = 1,...,N.$(48)
$v(x,T) = \int\limits_0^T {d\tau } \int\limits_0^l {G(x,\xi ;T,\tau ){{f}_{k}}(\xi )h_{\tau }^{k}(\xi ,\tau )d\xi } + \int\limits_0^l {G(x,\xi ;T,0){{f}_{k}}(\xi ){{h}^{k}}(\xi ,0)d\xi } .$(49)
$\left| {G(x,\xi ;T,\tau )} \right| \leqslant {{C}_{8}}{{(T - \tau )}^{{ - 3/4}}}{{\left| {x - \xi } \right|}^{{1/2}}}.$(50)
$\begin{gathered} {{f}_{{{{k}_{0}}}}}(x) = 1{\text{/}}{{h}^{{{{k}_{0}}}}}(x,T)\left( {\int\limits_0^T {d\tau } \int\limits_0^l {G(x,\xi ;T,\tau ){{f}_{k}}(\xi )h_{\tau }^{k}(\xi ,\tau )d\xi } + \int\limits_0^l {G(x,\xi ;T,0){{f}_{k}}(\xi ){{h}^{k}}(\xi ,0)d\xi } } \right), \\ x \in [{{a}_{{{{k}_{0}}}}},{{b}_{{{{k}_{0}}}}}],\quad {{k}_{0}} = 1,...,N. \\ \end{gathered} $(51)
$\left| {{{f}_{{{{k}_{0}}}}}(x)} \right| \leqslant C(T)\nu _{0}^{{3/2}}{\text{/}}{{h}_{T}}\left\| {{{f}_{k}}} \right\|\left( {\left\| {h_{\tau }^{k}} \right\| + \left\| {{{h}^{k}}(.,0)} \right\|} \right) \leqslant C(T)\nu _{0}^{{3/2}}{\text{/}}{{h}_{T}}\left\| f \right\|Nmax\left\{ {\left\| {h_{\tau }^{k}} \right\| + \left\| {{{h}^{k}}(.,0)} \right\|} \right\} = \lambda \left\| f \right\|.$Приведем достаточные условия единственности решения обратной задачи (38)–(40), имеющие глобальный характер и не предполагающие выполнения каких-либо условий малости.
Теорема 3. Пусть для функций $h$ выполнены условия ${{h}^{k}},h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, ${{h}^{{{{k}_{0}}}}}(x,t)h_{t}^{{{{k}_{0}}}}(x,t) \geqslant 0$, $(x,t) \in P_{T}^{{({{k}_{0}})}}$, ${{k}_{0}} = 1,...,N$. Тогда для единственности решения обратной задачи (38)–(40) необходимо и достаточно, чтобы носитель каждой из функций ${{h}^{k}}(.,T)$ содержал отрезок $[{{a}_{k}},{{b}_{k}}],k = 1,...,N$, т.е.
Замечание. Так как каждая из функций ${{f}_{k}}$ может быть отличной от нуля только на отрезке $[{{a}_{k}},{{b}_{k}}]$, то накладываемые на носители функций ${{h}^{k}}(.,T)$ ограничения сводятся к ограничениям на носитель следа этих функций на отрезках ${{\bar {\omega }}^{{(k)}}}$. Условия (52) теоремы требуют, чтобы замыкание множества точек ${{h}^{k}}(x,T) \ne 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, содержало отрезок $[{{a}_{k}},{{b}_{k}}],k = 1,...,N.$ Вне этих отрезков, естественно, никаких ограничений на поведение каждой из функций ${{h}^{k}}$ не накладывается. Кроме этого, справедливость неравенств ${{h}^{{{{k}_{0}}}}}(x,t)h_{t}^{{{{k}_{0}}}}(x,t) \geqslant 0$ можно также считать верной, без ограничения общности, при всех $(x,t) \in {{\Omega }_{T}}$, что и будет далее использоваться, не оговаривая это специально.Доказательство. Необходимость условий (52) для единственности решения обратной задачи (38)–(40) доказывается совершенно так же, как и доказательство аналогичного факта для задачи с финальным переопределением (на всей верхней крышке) (см. [14]). Докажем достаточность условия (52) для единственности решения. Рассуждаем от противного. Пусть условие (52) выполнено, но обратная задача (41)–(43) имеет нетривиальное решение, т.е. существует $N$ функций $f \in {{F}^{\alpha }}$, $f \ne 0$, таких, что решение прямой задачи (41), (42) удовлетворяет условиям (43). Сделаем сразу же очевидное замечание о характере правой части уравнения (41) в случае справедливости условия (43). Из принципа максимума следует, что правая часть этого уравнения обязательно должна менять знак на прямоугольнике ${{\Omega }_{T}}$. Для дальнейшего проведения доказательства представим каждую из функций ${{f}_{k}}$ в виде
(53)
${{f}_{k}}(x) = f_{k}^{ + }(x) - f_{k}^{ - }(x),\quad f_{k}^{ + }(x) = max\{ {{f}_{k}}(x),0\} ,\quad f_{k}^{ - }(x) = \left| {{{f}_{k}}(x)} \right| - f_{k}^{ + }(x),\quad x \in \Re .$(54)
${{h}^{k}}(x,t) = {{h}^{k}}(x,0) + \int\limits_0^t {h_{\tau }^{k}(x,\tau )d\tau } = {{({{h}^{k}}(x,0))}^{ + }} - {{({{h}^{k}}(x,0))}^{ - }} + \int\limits_0^t {({{{(h_{\tau }^{k}(x,\tau ))}}^{ + }} - ({{{(h_{\tau }^{k}(x,\tau ))}}^{ - }}))d\tau } .$Для преобразования формулы (54) понадобится
Лемма 2. Пусть для функции $p:{{\bar {\Omega }}_{T}} \to \Re $ справедливы включения $p,{{p}_{t}} \in C({{\bar {\Omega }}_{T}})$ и выполнены условия $p(x,t){{p}_{t}}(x,t) \geqslant 0$, $(x,t) \in {{\Omega }_{T}}$. Тогда для этой функции справедливы следующие свойства:
1) если при некотором $x \in \bar {\Omega }$ выполнено $p(x,T) = 0$, то $p(x,t) = 0$ при $0 \leqslant t \leqslant T$,
2) если при некотором $x \in \bar {\Omega }$ выполнено $p(x,T) > 0$, то $p(x,t) \geqslant 0$, ${{p}_{t}}(x,t) \geqslant 0$ при $0 \leqslant t \leqslant T$,
3) если при некотором $x \in \bar {\Omega }$ выполнено $p(x,T) < 0$, то $p(x,t) \leqslant 0$, ${{p}_{t}}(x,t) \leqslant 0$ при $0 \leqslant t \leqslant T$,
4) функции ${{p}^{ + }}(x,t)$, ${{p}^{ - }}(x,t)$ непрерывно дифференцируемы по $t$ на ${{\bar {\Omega }}_{T}}$ при этом
(55)
${{({{p}^{ + }}(x,t))}_{t}} = {{(p{{(x,t)}_{t}})}^{ + }},\quad {{({{p}^{ - }}(x,t))}_{t}} = {{(p{{(x,t)}_{t}})}^{ - }}.$Доказательство. Докажем, например, свойство 1). От противного. Пусть в некоторой точке $(x{\text{*}},T)$ будет выполнено равенство $p(x{\text{*}},T) = 0$, но при этом в некоторой точке $(x{\text{*}},t{\text{*}}) \in \{ x{\text{*}}\} \times [0,T)$ будет справедливо, например, неравенство $p(x{\text{*}},t{\text{*}}) > 0$. Рассмотрим множество чисел $G = \{ t:t{\text{*}} < t \leqslant T,p(x{\text{*}},t) = 0\} $. Оно, очевидно, замкнуто и, следовательно, его нижняя грань ему принадлежит. Пусть это есть число $\bar {t} > t{\text{*}}$. Из определения этого числа следует, что при $t{\text{*}} \leqslant t < \bar {t}$ справедливо неравенство $p(x{\text{*}},t) > 0$. Тогда, по теореме Лагранжа, существует такое число $\xi \in (t{\text{*}},\bar {t})$, что выполнено равенство
(56)
$p(x{\text{*}},\bar {t}) - p(x{\text{*}},t{\text{*}}) = {{p}_{t}}(x{\text{*}},\xi )(\bar {t} - t{\text{*}}) = - p(x{\text{*}},t{\text{*}}) < 0.$Используя лемму 2, представим уравнение (41) в следующем виде:
(57)
$\begin{gathered} (Lu)(x,t) = {{f}_{k}}(x){{h}^{k}}(x,t) = \left[ {{{f}_{k}}{{{(x)}}^{ + }}{{{({{h}^{k}}(x,0))}}^{ + }} + {{f}_{k}}{{{(x)}}^{ - }}{{{({{h}^{k}}(x,0))}}^{ - }} + {{f}_{k}}{{{(x)}}^{ + }}\int\limits_0^t {{{{(h_{\tau }^{k}(x,\tau ))}}^{ + }}} + } \right. \\ \left. { + \;{{f}_{k}}{{{(x)}}^{ - }}\int\limits_0^t {{{{(h_{\tau }^{k}(x,\tau ))}}^{ - }}d\tau } } \right] - \left[ {{{f}_{k}}{{{(x)}}^{ + }}{{{({{h}^{k}}(x,0))}}^{ - }} + {{f}_{k}}{{{(x)}}^{ - }}{{{({{h}^{k}}(x,0))}}^{ + }} + {{f}_{k}}{{{(x)}}^{ + }}\int\limits_0^t {{{{(h_{\tau }^{k}(x,\tau ))}}^{ - }}} } \right. + \\ \left. {\, + {{f}_{k}}{{{(x)}}^{ - }}\int\limits_0^t {{{{(h_{\tau }^{k}(x,\tau ))}}^{ + }}d\tau } } \right] = {{q}_{ + }}(x,t) - {{q}_{ - }}(x,t),(x,t) \in {{\Omega }_{T}}. \\ \end{gathered} $Определим функции ${{u}_{ + }},{{u}_{ - }} \in U({{\Omega }_{T}})$ из условий
(59)
${{u}_{ + }}(0,t) = 0,\quad {{u}_{ + }}(0,t) = 0,\quad t \in [0,T],\quad {{u}_{ + }}(x,0) = 0,\quad x \in \bar {\Omega },$(61)
${{u}_{ - }}(0,t) = 0,\quad {{u}_{ - }}(0,t) = 0,\quad t \in [0,T],\quad {{u}_{ - }}(x,0) = 0,\quad x \in \bar {\Omega }.$Как следует из формулы (57), для функций ${{q}_{ + }},{{q}_{ - }}$ справедливы включения ${{q}_{ + }},{{q}_{ - }},{{({{q}_{ + }})}_{t}}$, ${{({{q}_{ - }})}_{t}} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, а тогда по теореме 1 прямые задачи (58), (59), (60), (61) действительно имеют единственное решение в указанном классе. По следствию 2 теоремы 1, примененному к этим задачам, функции ${{v}_{ + }} = {{({{u}_{ + }})}_{t}}$, ${{v}_{ - }} = {{({{u}_{ - }})}_{t}} \in V(\Omega _{T}^{*})$ и удовлетворяют задачам
(63)
${{v}_{ + }}(0,t) = 0,\quad {{v}_{ + }}(l,t) = 0,\quad t \in [0,T],\quad {{v}_{ + }}(x,0) = {{q}_{ + }}(x,0),\quad x \in \bar {\Omega },\quad x \notin A \cup B.$(65)
${{v}_{ - }}(0,t) = 0,\quad {{v}_{ - }}(l,t) = 0,\quad t \in [0,T],\quad {{v}_{ - }}(x,0) = {{q}_{ - }}(x,0),\quad x \in \bar {\Omega },\quad x \notin A \cup B.$(66)
${{u}_{ + }}(x,T) = {{u}_{ - }}(x,T) = {{\psi }^{k}}(x),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(67)
${{v}_{ + }}(x,T) = \psi _{{xx}}^{k}(x) + {{q}_{ + }}(x,T),\quad {{v}_{ - }}(x,T) = \psi _{{xx}}^{k}(x) + {{q}_{ - }}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(68)
${{v}_{ + }}(x{\text{*}},T) - \psi _{{xx}}^{k}(x{\text{*}}) = {{q}_{ + }}(x{\text{*}}T),\quad {{v}_{ - }}(x{\text{*}},T) - \psi _{{xx}}^{k}(x{\text{*}}) = {{q}_{ - }}(x{\text{*}}T).$(69)
${{q}_{ + }}(x{\text{*}},T) = f_{1}^{ + }(x{\text{*}}){{({{h}^{1}}(x{\text{*}},0))}^{ + }} + {{f}_{1}}{{(x{\text{*}})}^{ + }}\int\limits_0^T {{{{(h_{\tau }^{1}(x{\text{*}},\tau ))}}^{ + }}d\tau } = f_{1}^{ + }(x{\text{*}}){{({{h}^{1}}(x{\text{*}},T))}^{ + }} > 0.$3. ТЕОРЕМЫ ФРЕДГОЛЬМА И ОБРАТНАЯ ЗАДАЧА (38)–(40)
Вопрос о разрешимости обратной задачи (38)–(40) тесно связан с разрешимостью однородной обратной задачи, соответствующей этой задаче. Однородной обратной задачей, соответствующей задаче (38)–(40), будем называть задачу (38)–(40) при $\varphi = 0$, ${{\mu }_{1}} = 0$, ${{\mu }_{2}} = 0$, $g = 0$, ${{\chi }_{k}} = 0$, $k = 1,...,N$, т.е. задачу определения функций $(u,f) \in U({{\Omega }_{T}}) \times {{F}^{\alpha }}$ из условий (41)–(43). Связь между задачами (38)–(40) и (41)–(43) устанавливает следующее утверждение, указывающее на некоторую аналогию между теорией Фредгольма, справедливой для уравнений II рода с компактным оператором в банаховом пространстве, и характером разрешимости обратной задачи (38)–(40).
Теорема 4. Пусть для функций $h$ справедливы включения ${{h}^{k}},h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, при этом $\left| {{{h}^{k}}(x,T)} \right| \geqslant {{h}_{T}} > 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, где ${{h}_{T}} > 0$ – некоторая постоянная. Тогда однородная обратная задача (41)–(43) может иметь только конечное число линейно независимых решений. Если дополнительно к этому выполнены условия единственности, сформулированные в теореме 2 или теореме 3, то для любых функций $\varphi \in C(\bar {\Omega })$, ${{\mu }_{1}},{{\mu }_{2}} \in C[0,T]$, удовлетворяющих условиям согласования ${{\mu }_{1}}(0) = \varphi (0)$, ${{\mu }_{2}}(0) = \varphi (l)$ и любых функций $g \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, ${{\theta }_{k}} \in {{C}^{\alpha }}[{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, неоднородная обратная задача (38)–(40) имеет единственное решение, пару функций $(u,f)$, удовлетворяющую условиям (38), (39) и условиям переопределения (40) с некоторыми, однозначно определенными функциями $\chi \in {{H}^{{2,\alpha }}}$, при этом для этих функций справедливы равенства ${{({{\chi }_{k}}(x))}_{{xx}}} = {{\theta }_{k}}(x)$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$.
Доказательство. Покажем сначала, что если выполнены условия $\left| {{{h}^{k}}(x,T)} \right| \geqslant {{h}_{T}} > 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, то функция $f$ удовлетворяет в банаховом пространстве ${{F}^{\alpha }}$ уравнению Фредгольма II рода с компактным оператором. Определим линейный оператор $A\,:{{F}^{\alpha }} \to {{F}^{\alpha }}$ по следующему правилу:
(70)
${{(A(f))}_{k}}(x) = v(x,T;f){\text{/}}{{h}^{k}}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(71)
${{f}_{k}}(x) = {{(A(f))}_{k}}(x) - [{{({{\chi }_{k}})}_{{xx}}}(x) - {{({{u}_{0}})}_{{xx}}}(x,T)]{\text{/}}{{h}^{k}}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(72)
${{p}_{k}}(x) = - [{{({{\chi }_{k}})}_{{xx}}}(x) - {{({{u}_{0}})}_{{xx}}}(x,T)]{\text{/}}{{h}^{k}}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$В условиях (76) для получившегося как следы решения прямой задачи (74), (75) на отрезках ${{\omega }_{k}}$ набора функций $\zeta = {{\zeta }_{k}}(x)$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, в соответствии с теоремой 1, справедливо включение $\zeta \in {{H}^{{2,\alpha }}}$. В соответствии со следствием 2 теоремы 1 для функции $v = {{u}_{t}}$ будет справедлива следующая задача:
(77)
$(Lv)(x,t) = {{v}_{t}}(x,t) - {{v}_{{xx}}}(x,t) = {{f}_{k}}(x)h_{t}^{k}(x,t),\quad (x,t) \in {{\Omega }_{T}},$(78)
$v(0,t) = 0,\quad v(l,t) = 0,\quad t \in [0,T],\quad v(x,0) = {{f}_{k}}(x){{h}^{k}}(x,0),\quad x \in \bar {\Omega },\quad x \notin A \cup B.$(79)
$v(x,T) = {{({{\zeta }_{k}}(x))}_{{xx}}} + {{f}_{k}}(x){{h}^{k}}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(80)
${{(f - Af)}_{k}}(x) = {{({{\zeta }_{k}})}_{{xx}}}(x){\text{/}}{{h}^{k}}(x,T),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(81)
${{\psi }^{k}}(x) = {{u}_{ + }}(x,T) = {{u}_{ - }}(x,T) + {{\zeta }_{k}}(x),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$(82)
${{v}_{ + }}(x{\text{*}},T) - \psi _{{xx}}^{k}(x{\text{*}}) = {{q}_{ + }}(x{\text{*}},T),\quad {{v}_{ - }}(x{\text{*}},T) - \psi _{{xx}}^{k}(x{\text{*}}) - {{({{\zeta }_{k}})}_{{xx}}}(x{\text{*}}) = {{q}_{ - }}(x{\text{*}},T).$4. УСЛОВИЯ ОДНОЗНАЧНОЙ РАЗРЕШИМОСТИ ДЛЯ ОБРАТНОЙ ЗАДАЧИ (38)–(40).
Рассмотрим предварительно простой пример обратной задачи (38)–(40) при $N = 1$, $h = 1$, $0 < a < b < l$. В этом случае требуется найти пару функций $(u,f) \in U({{\Omega }_{T}}) \times {{C}^{\alpha }}[a,b]$ из условий
В уравнении (83), в соответствии с предыдущими обозначениями, функция $h(x,t) \equiv 1$. Так как ${{h}_{t}}(x,t) \equiv 0$, то в соответствии с теоремой 4 обратная задача (83)–(85) имеет единственное решение, определяемое с помощью функции ${{\chi }_{{xx}}}$. Уравнение для определения функции $f$, в соответствии с рассмотрениями этой теоремы, имеет вид В уравнении (86) функция $v(x,t) = v(x,t;f)$ является решением прямой задачи(88)
$v(0,t) = 0,\quad v(l,t) = 0,\quad t \in [0,T],\quad v(x,0) = f(x),\quad x \in \bar {\Omega },\quad x \ne a,b.$Проведем аналогичные построения в общем случае $N > 1$ для обратной задачи (38)–(40). Будем предполагать выполненными для функций $h$ условия теоремы 4, т.е. $\left| {{{h}^{k}}(x,T)} \right| \geqslant {{h}_{T}} > 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, и условия единственности, сформулированные в теореме 2 или теореме 3. В этом случае определим оператор $B\,:{{F}^{\alpha }} \to {{F}^{\alpha }}$ для любых функций $\mu \in {{F}^{\alpha }}$по правилу
(93)
$(B\mu )(x) = ({{(E - A)}^{{ - 1}}}\mu )(x),\quad x \in [{{a}_{k}},{{b}_{k}}],\quad k = 1,...,N.$Будем говорить, что функции $\chi \in {{H}^{{2,\alpha }}}$ удовлетворяют условиям согласования, если для них выполнены равенства
(94)
$\chi ({{a}_{k}}) = u({{a}_{k}},T; - B({{\chi }_{{xx}}} + {{u}_{0}}(x,T))),\quad \chi ({{b}_{k}}) = u({{b}_{k}},T; - B({{\chi }_{{xx}}} + {{u}_{0}}(x,T))),\quad k = 1,...{\kern 1pt} ,N.$Роль условий согласования в этом случае проясняет следующее утверждение.
Лемма 3. Пусть для функций $h$ выполнены условия теоремы 4, справедливы включения $\varphi \in C(\bar {\Omega })$, ${{\mu }_{1}},{{\mu }_{1}} \in C[0,T]$, $g \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, ${{f}_{k}} \in {{C}^{\alpha }}[{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, выполнены условия согласования ${{\mu }_{1}}(0) = \varphi (0)$, ${{\mu }_{2}}(0) = \varphi (l)$. Тогда функции $\chi \in {{H}^{{2,\alpha }}}$ будут следами решений прямой задачи (1), (2) на отрезках ${{\bar {\omega }}_{k}},k = 1,...,N$, в том и только в том случае, когда для них будут выполнены условия согласования (94).
Доказательство леммы 3 проводится совершенно аналогично доказанному выше аналогичному утверждению для случая $N = 1$.
Для обратной задачи (38)–(40) получаем окончательно следующую теорему существования и единственности.
Теорема 5. Пусть для функций $h$ справедливы включения ${{h}^{k}},h_{t}^{k} \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, при этом $\left| {{{h}^{k}}(x,T)} \right| \geqslant {{h}_{T}} > 0$, $x \in [{{a}_{k}},{{b}_{k}}]$, $k = 1,...,N$, где ${{h}_{T}} > 0$ – некоторая постоянная, выполнены условия единственности, сформулированные в теореме 2 или теореме 3, то для любых функций $\varphi \in C(\bar {\Omega })$, ${{\mu }_{1}},{{\mu }_{2}} \in C[0,T]$, удовлетворяющих условиям согласования ${{\mu }_{1}}(0) = \varphi (0)$, ${{\mu }_{2}}(0) = \varphi (l)$, любых функций $g \in {{C}^{{\alpha ,\alpha /2}}}({{\bar {\Omega }}_{T}})$, $\chi \in {{H}^{{2,\alpha }}}$, удовлетворяющих условиям согласования (94), обратная задача (38)–(40) имеет единственное решение, пару функций $(u,f)$, удовлетворяющую условиям (38)–(40).
Доказательство теоремы 5 проводится совершенно аналогично случаю $N = 1$, рассмотренному в начале раздела.
В заключение автор выражает глубокую благодарность всем участникам семинара по обратным задачам на механико-математическом факультете МГУ под руководством В.А. Садовничего и А.И. Прилепко за интерес, проявленный к изложенным выше исследованиям.
Список литературы
Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967.
Фридман А. Уравнения с частными производными параболического типа. М.: Мир, 1968.
Ландис Е.М. Уравнения второго порядка эллиптического и параболического типов. М.: Наука, 1971.
Смирнов В.И. Курс высшей математики, т. 4. М.: ГИФМЛ, 1958.
Смирнов В.И. Курс высшей математики, т. 2. М.: ГИФМЛ, 1961.
Михайлов В.П. Дифференциальные уравнения в частных производных. М.: Наука, 1983.
Денисов А.М. Задачи определения неизвестного источника в параболическом и гиперболическом уравнениях // Ж. вычисл. матем. и матем. физ. 2015. Т. 55. № 5. С. 830–835.
Денисов А.М. Введение в теорию обратных задач. М.: Изд-во МГУ, 1994.
Лаврентьев М.М. Условно-корректные задачи для дифференциальных уравнений. Новосибирск: Изд-во НГУ, 1973.
Prilepko A.I., Orlovsky D.G., Vasin I.A. Methods for solving inverse problems in mathematical physics. New York–Basel: Marcel Dekker Inc. 2000.
Isakov V. Inverse problems for Partial Differential Equations. New York.: Springer, 2006.
Костин А.Б. Обратная задача восстановления источника в параболическом уравнении по условиям нелокального наблюдения // Матем. сб. 2013. Т. 204. № 10. С. 3–46.
Соловьёв В.В. О разрешимости обратной задачи определения источника с переопределением на верхней крышке для параболического уравнения // Дифференц. ур-ния. 1989. Т. 25. № 9. С. 1577–1583.
Соловьёв В.В., Ткаченко Д.С. Определение источника в уравнении теплопроводности для случая негладких краевых и начальных условий // Вестник национального исследовательского ядерного университета “МИФИ”, 2017. Т. 6. С. 1–9.
Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1981.
Дополнительные материалы отсутствуют.
Инструменты
Журнал вычислительной математики и математической физики