Доклады Российской академии наук. Математика, информатика, процессы управления, 2021, T. 501, № 1, стр. 11-15
ЕДИНСТВЕННОСТЬ ВЕРОЯТНОСТНОГО РЕШЕНИЯ УРАВНЕНИЯ КОЛМОГОРОВА С МАТРИЦЕЙ ДИФФУЗИИ, УДОВЛЕТВОРЯЮЩЕЙ УСЛОВИЮ ДИНИ
В. И. Богачев 1, 2, 3, 4, *, С. В. Шапошников 1, 2, 4, **
1 Московский государственный университет
имени М.В. Ломоносова
Москва, Россия
2 Национальный исследовательский университет “Высшая школа экономики
Москва, Россия
3 Православный Свято-Тихоновский
гуманитарный университет
Москва, Россия
4 Московский центр фундаментальной
и прикладной математики
Москва, Россия
* E-mail: vibogach@mail.ru
** E-mail: starticle@mail.ru
Поступила в редакцию 25.06.2021
После доработки 25.06.2021
Принята к публикации 22.07.2021
Аннотация
В этом сообщении мы исследуем стационарное уравнение Колмогорова и доказываем, что в случае, когда матрица диффузии удовлетворяет условию Дини, а коэффициент сноса локально интегрируем в степени выше размерности, отношение двух вероятностных решений входит в класс Соболева, а при наличии функции Ляпунова или глобальной интегрируемости коэффициентов относительно решения вероятностное решение единственно.
В этом сообщении мы исследуем стационарное уравнение Колмогорова
(1)
${{\partial }_{{{{x}_{i}}}}}{{\partial }_{{{{x}_{j}}}}}({{a}^{{ij}}}\varrho ) - {{\partial }_{{{{x}_{i}}}}}({{b}^{i}}\varrho ) = 0$Предположим, что коэффициенты удовлетворяют следующим условиям.
$({{{\mathbf{H}}}_{{\mathbf{a}}}})$ Матрица $A(x) = ({{a}^{{ij}}}(x{{))}_{{i,j \leqslant d}}}$ симметрична и положительно определена, причем для всякого шара $B \subset {{\mathbb{R}}^{d}}$ есть такая возрастающая непрерывная функция ${{\omega }_{B}}$ на $[0, + \infty )$, что $\omega (0) = 0$ и
Кроме того, для каждого шара $B$ есть такое число ${{\nu }_{B}} > 0$, что
где I – единичный оператор.$({{{\mathbf{H}}}_{{\mathbf{b}}}})$ Для всякого шара $B$ есть такое число $p = p(B) > d$, что ${\text{|}}b{\text{|}} \in {{L}^{p}}(B)$.
Матричный след обозначим через tr и положим
где $\nabla u\, = \,{{({{\partial }_{{{{x}_{i}}}}}u)}_{{i \leqslant d}}}$, ${{\nabla }^{2}}u\, = \,{{({{\partial }_{{{{x}_{i}}}}}{{\partial }_{{{{x}_{j}}}}}u)}_{{i,j \leqslant d}}}$. Функция $\varrho $ ∈ ∈ $L_{{{\text{loc}}}}^{1}({{\mathbb{R}}^{d}})$ является решением уравнения (1), еслиРешение $\varrho $ называется вероятностным, если $\varrho \geqslant 0$ и интеграл от $\varrho $ равен 1. Через ${{W}^{{p,1}}}(B)$ и ${{W}^{{p,2}}}(B)$ обозначим классы Соболева на открытом шаре B, состоящие из функций, лежащих в ${{L}^{p}}(B)$ вместе с обобщенными частными производными первого и второго порядка соответственно, наделенные их стандартными соболевскими нормами. Через $W_{0}^{{p,1}}(B)$ обозначим замыкание $C_{0}^{\infty }(B)$ в ${{W}^{{p,1}}}(B)$.
Проблеме единственности вероятностного решения посвящены, в частности, работы [1–3]. В них предполагается, что функции ${{a}^{{ij}}}$ входят в локальный класс Соболева $W_{{{\text{loc}}}}^{{p,1}}({{\mathbb{R}}^{d}})$ с некоторым $p > d$. Это условие существенно используется при обосновании единственности. Благодаря ему удается применить неравенство Харнака для решений дивергентных уравнений и получить положительность непрерывной версии вероятностного решения. Доказательство единственности в цитированных работах основано на рассмотрении функции ${v}$, равной отношению двух вероятностных решений $\sigma $ и $\varrho $, и анализе дивергентного эллиптического уравнения
Пусть $B$ – открытый единичный шар в ${{\mathbb{R}}^{d}}$, $\overline B $ – его замыкание и $p = p(B) > d$ – показатель интегрируемости из условия ${\text{(}}{{{\mathbf{H}}}_{{\mathbf{b}}}}{\text{)}}$. Следующее утверждение известно в случае ограниченных коэффициентов и в случае нулевого сноса (см. [9, 10]), для полноты изложения приведем доказательство.
Лемма 1. Для всякой функции $f \in {{L}^{p}}(B)$ в ${{W}^{{p,2}}}(B) \cap W_{0}^{{p,1}}(B) \cap {{C}^{1}}(\overline B )$ существует решение уравнения $Lu = f$.
Доказательство. Пусть ${{L}_{A}}u = {\text{tr(}}A{{\nabla }^{2}}u)$. Тогда в силу [9, лемма 9.17] для всякого u ∈ ∈ ${{W}^{{p,2}}}(B) \cap W_{0}^{{p,1}}(B)$ верна оценка
При этом верна оценка (см. [11, глава 5])
По принципу максимума (см. [9, теорема 9.1]) величина ${{\left\| u \right\|}_{{{{L}^{\infty }}(B)}}}$ оценивается сверху через ${{C}_{3}}{{\left\| {Lu} \right\|}_{{{{L}^{p}}(B)}}}$. Значит, для всех $u \in {{W}^{{p,2}}}(B) \cap W_{0}^{{p,1}}(B)$ имеем
В силу [9, теорема 9.15] существует единственная функция u из пространства ${{W}^{{p,2}}}(B) \cap W_{0}^{{p,1}}(B)$, удовлетворяющая уравнению ${{L}_{A}}u = f$. Теперь стандартным методом продолжения по параметру выводятся существование и единственность решения уравнения $Lu = f$. Из теоремы вложения получаем $u \in {{C}^{1}}(\overline B )$.
В следующем утверждении мы доказываем единственность решения задачи Дирихле на $B$ для дважды дивергентного уравнения
(2)
${{\partial }_{{{{x}_{i}}}}}{{\partial }_{{{{x}_{j}}}}}({{a}^{{ij}}}u) - {{\partial }_{{{{x}_{i}}}}}({{b}^{i}}u) = 0,\quad u{{{\text{|}}}_{{\partial B}}} = g,$Отметим, что это определение решения отличается от определения в [8] и не включает в интегральное тождество граничные слагаемые.
Лемма 2. Решение задачи Дирихле (2) единственно.
Доказательство. Покажем, что в случае g = 0 задача Дирихле имеет лишь нулевое решение. Пусть $f \in C_{0}^{\infty }(B)$. По лемме 1 существует решение ${v} \in {{W}^{{p,2}}}(B) \cap W_{0}^{{p,1}}(B)$ уравнения
Положим ${{r}_{k}} = 1 - {{k}^{{ - 1}}}$. Рассмотрим последовательность таких функций ${{\zeta }_{k}} \in C_{0}^{\infty }(B)$, что ${{\zeta }_{k}}(x)$ = 1 при ${\text{|}}x{\text{|}} < {{r}_{k}}$, $0 \leqslant {{\zeta }_{k}} \leqslant 1$, ${\text{|}}\nabla {{\zeta }_{k}}{\text{|}} \leqslant {{C}_{1}}k$ и $\left\| {{{\nabla }^{2}}{{\zeta }_{k}}} \right\| \leqslant {{C}_{1}}{{k}^{2}}$. Имеем
Так как ${v} \in {{C}^{1}}(\overline B )$ и ${v} = 0$ на $\partial B$, то ${\text{|}}{v}(x){\text{|}} \leqslant {{C}_{2}}{{k}^{{ - 1}}}$ при ${{r}_{k}} < {\text{|}}x{\text{|}} < 1$. Заметим, что объем множества $\{ x:{{r}_{k}} < {\text{|}}x{\text{|}} < 1\} $ оценивается через ${{C}_{3}}{{k}^{{ - 1}}}$. Кроме того, ${\text{|}}\nabla {{\zeta }_{k}}(x){\text{|}} = 0$ и ${{\nabla }^{2}}{{\zeta }_{k}}(x) = 0$ при ${\text{|}}x{\text{|}} < {{r}_{k}}$. Следовательно,
Таким образом,
В дальнейших рассуждениях важную роль будет играть процедура приближения решения уравнения (1) последовательностью решений аналогичных уравнений, но с гладкими матрица-ми A. Пусть $\varrho $ – непрерывное решение урав-нения (1) на ${{\mathbb{R}}^{d}}$. Пусть $\psi \in C_{0}^{\infty }({{\mathbb{R}}^{d}})$, $\psi \geqslant 0$, $\psi (x) = 0$ при ${\text{|}}x{\text{|}} > 1$ и ${{\left\| \psi \right\|}_{{{{L}^{1}}({{\mathbb{R}}^{d}})}}} = 1$. Положим ${{\psi }_{k}}(x) = {{k}^{d}}\psi (kx)$, где $k \in \mathbb{N}$. Пусть
Обозначим через $B{\kern 1pt} '$ шар $2B$. Для всех $x,y \in B$ имеем
Кроме того, равномерно на $\overline B $ отображения Ak сходятся к A, а функции ${{\varrho }_{k}}$ сходятся к $\varrho $. Поскольку уравнение
Замечание 1. Из доказательства [6, теоре- ма 1.8] можно усмотреть, что
Лемма 3. Функции ${{u}_{k}}$ равномерно на $\overline B $ сходятся к $\varrho $.
Доказательство. В силу [6, теорема 1.8] (см. замечание 1) функции ${{u}_{k}}$ равномерно ограничены и равностепенно непрерывны на $\overline B $. Следовательно, есть подпоследовательность $\{ {{u}_{{{{k}_{j}}}}}\} $, которая сходится равномерно на $\overline B $ к функции $u \in C(\overline B )$. Ясно, что $u = \varrho $ на $\partial B$. Пусть $\varphi \in C_{0}^{\infty }(B)$. Имеем
Полагая ${{k}_{j}} \to \infty $, получаем аналогичное равенство с $A$ и $u$ вместо ${{A}_{k}}$ и ${{u}_{k}}$. Итак, функция $u$ является решением задачи Дирихле (2) на $B$ с граничным условием $u{{{\text{|}}}_{{\partial B}}} = \varrho $. Другим решением этой задачи является функция $\varrho $. По лемме 2 имеем $u = \varrho $. Так как предел для всякой сходящейся подпоследовательности $\{ {{u}_{{{{k}_{j}}}}}\} $ равен $\varrho $, то вся последовательность $\{ {{u}_{k}}\} $ сходится к $\varrho $.
Пусть $\varrho $ и $\sigma $ – вероятностные решения уравнения (1). Согласно [4, теорема 3.1] функции $\varrho $ и $\sigma $ непрерывны на ${{\mathbb{R}}^{d}}$. Кроме того, для всякого шара $B$ есть такое число $C(B) > 0$, что $\varrho (x) > C(B)$ и $\sigma (x) > C(B)$ при $x \in B$. Это следует из неравенства Харнака. Если выполнено условие $({{H}_{a}})$, а $b$ – локально ограниченное векторное поле, то нужное неравенство Харнака доказано в [4, следствие 3.6]. В готовящейся к публикации работе условие локальной ограниченности $b$ заменено на $({{H}_{b}})$.
Положим $v = \frac{\sigma }{\varrho }$. Пусть $f \in {{C}^{2}}([0, + \infty ))$ и $f{\kern 1pt} '{\kern 1pt} ' \geqslant 0$.
Теорема 1. Верно включение $v \in W_{{{\text{loc}}}}^{{2,1}}({{\mathbb{R}}^{d}})$, причем для всякой неотрицательной функции $\psi \in C_{0}^{\infty }({{\mathbb{R}}^{d}})$ выполнено неравенство
(3)
$\int \,\varrho {{\left| {\sqrt A \nabla v} \right|}^{2}}f{\kern 1pt} '{\kern 1pt} '({v})\psi dx \leqslant \int \,\varrho f({v})L\psi dx.$Доказательство. Пусть $B = B({{x}_{0}},r)$ – шар, содержащий носитель функции $\psi $ и $B{\kern 1pt} ' = B({{x}_{0}},r + 1)$. Как отмечено выше, существует такое число $C(B) > 0$, что $\sigma > C(B)$ и $\varrho > C(B)$ на $B$. Пусть $\{ {{u}_{k}}\} $ и $\{ {{w}_{k}}\} $ – построенные перед леммой 3 последовательности решений уравнения
(4)
$\int \,{{u}_{k}}{{\left| {\sqrt {{{A}_{k}}} \nabla {{{v}}_{k}}} \right|}^{2}}f{\kern 1pt} '{\kern 1pt} '({{{v}}_{k}})\psi dx = \int \,{{u}_{k}}f({{{v}}_{k}})L\psi dx.$Пусть $f(t) = {{t}^{2}}$. Тогда
Заметим, что правая часть оценивается некоторой константой $C(\psi )$, которая не зависит от $k$. Пусть $Q$ – шар, причем $\overline Q \subset B$ и $\psi = 1$ на $Q$. Тогда
Следовательно, существует подпоследовательность $\{ {{v}_{{{{k}_{j}}}}}\} $, слабо сходящаяся в ${{W}^{{2,1}}}(Q)$ к некоторой функции $v$. Таким образом, $v \in {{W}^{{2,1}}}(Q)$. Из-за равномерной сходимости вся последовательность $\{ {{v}_{k}}\} $ сходится слабо в ${{W}^{{2,1}}}(Q)$ к $v$. Пусть теперь $\psi $ – произвольная неотрицательная функция из $C_{0}^{\infty }(B)$. Можно считать, что $Q$ содержит носитель $\psi $. Функции $f{\kern 1pt} '{\kern 1pt} '({{v}_{k}})$ равномерно на $B$ сходятся к $f{\kern 1pt} '{\kern 1pt} '(v)$. Так как
Используя слабую сходимость ${{{v}}_{k}}$ в ${{W}^{{2,1}}}(Q)$ и равномерную сходимость ${{A}_{k}}$, ${{u}_{k}}$, ${{v}_{k}}$, $f{\kern 1pt} '{\kern 1pt} '({{{v}}_{k}})$ и $f({{v}_{k}})$, с учетом (4) при $k \to \infty $ получаем (3).
Замечание 2. Из доказательства видно, что теорема верна для произвольных положительных решений на области $\Omega \subset {{\mathbb{R}}^{d}}$, при условии, что в определяющее равенство подставляются функции $\psi $ с носителем в $\Omega $. Таким образом, для всяких положительных решений $\varrho $ и $\sigma $ уравнения Колмогорова на области $\Omega $ их отношение $\sigma {\text{/}}\varrho $ лежит в $W_{{{\text{loc}}}}^{{2,1}}(\Omega )$, хотя сами функции $\varrho $ и $\sigma $ могут не иметь соболевских производных.
Теперь, используя неравенство (3) и дословно повторяя рассуждения из [3, теорема 4.1.6], приходим к нашему основному результату.
Теорема 2. Пусть $\varrho $ – вероятностное решение уравнения (1) и выполнено одно из следующих условий:
(i) ${{(1 + \left| x \right|)}^{{ - 2}}}{\text{|}}{{a}^{{ij}}}(x){\text{|}},(1 + \left| x \right|{{)}^{{ - 1}}}{\text{|}}{{b}^{i}}(x){\text{|}} \in {{L}^{1}}(\varrho dx)$,
(ii) найдется функция $V \in {{C}^{2}}({{\mathbb{R}}^{d}})$ с $\mathop {\lim }\limits_{|x| \to \infty } V(x)\, = \, + \,\infty $ и $LV \leqslant {{C}_{1}} + {{C}_{2}}V$.
Тогда $\varrho $ – единственное вероятностное решение.
Например, это верно, если коэффициенты оцениваются через $C + C{\text{|}}x{\text{|}}$.
Замечание 3. Условие Дини является значительно более общим по сравнению с включением ${{a}^{{ij}}} \in W_{{{\text{loc}}}}^{{p,1}}$, $p > d$, поскольку последнее влечет гёльдеровость решения. Однако нам неизвестно, оптимально ли данное условие единственности. В частности, неизвестен ответ на следующий вопрос. Пусть $\nu \cdot I \leqslant A(x) \leqslant {{\nu }^{{ - 1}}} \cdot I$ для всех $x$, ${{a}^{{ij}}}$ – непрерывные функции и $b(x) = - x$. Может ли стационарное уравнение Колмогорова с такими коэффициентами иметь несколько вероятностных решений? О неединственности см. [12].
Список литературы
Богачёв В.И., Рёкнер М., Штаннат В. // Матем. сб. 2002. Т. 197. № 7. С. 3–36.
Bogachev V.I., Röckner M., Shaposhnikov S.V. // J. Math. Sci. (New York). 2011. V. 176. № 6. P. 759–773.
Bogachev V.I., Krylov N.V., Röckner M., Shaposhnikov S.V. Fokker–Planck–Kolmogorov equations. Amer. Math. Soc., Providence, Rhode Island, 2015.
Bogachev V.I., Shaposhnikov S.V. // Annali Matem. 2017. V. 196. P. 1609–1635.
Dong H., Kim S. // Comm. Partial Differ. Equ. 2017. V. 42. P. 417–435.
Dong H., Escauriaza L., Kim S. // Math. Ann. 2018. V. 370. P. 447–489.
Bauman P. // Ark. Mat. 1984. V. 22. P. 153–173.
Escauriaza L., Montaner S. // Rend. Lincei Mat. Appl. 2017. V. 28. P. 49–63.
Гилбарг Д., Трудингер Н. Эллиптические дифференциальные уравнения с частными производными второго порядка. М.: Наука, 1989.
Chiarenza F., Frasca M., Longo P. // Trans. Amer. Math. Soc. 1993. V. 336. P. 841–853.
Мазья В.Г. Пространства С.Л. Соболева. Л.: Изд-во ЛГУ, 1985.
Красовицкий Т.И. // Докл. РАН. 2019. Т. 487. № 4. С. 361–364.
Дополнительные материалы отсутствуют.
Инструменты
Доклады Российской академии наук. Математика, информатика, процессы управления