Математические заметки
RUS  ENG    ЖУРНАЛЫ   ПЕРСОНАЛИИ   ОРГАНИЗАЦИИ   КОНФЕРЕНЦИИ   СЕМИНАРЫ   ВИДЕОТЕКА   ПАКЕТ AMSBIB  
Общая информация
Последний выпуск
Скоро в журнале
Архив
Импакт-фактор
Правила для авторов
Лицензионный договор
Загрузить рукопись

Поиск публикаций
Поиск ссылок

RSS
Последний выпуск
Текущие выпуски
Архивные выпуски
Что такое RSS



Матем. заметки:
Год:
Том:
Выпуск:
Страница:
Найти






Персональный вход:
Логин:
Пароль:
Запомнить пароль
Войти
Забыли пароль?
Регистрация


Математические заметки, 2023, том 113, выпуск 2, страницы 217–235
DOI: https://doi.org/10.4213/mzm13882
(Mi mzm13882)
 

Эта публикация цитируется в 4 научных статьях (всего в 4 статьях)

Об асимптотике решений двучленных дифференциальных уравнений

Н. Н. Конечнаяa, К. А. Мирзоевb, А. А. Шкаликовb

a Северный (Арктический) федеральный университет имени М. В. Ломоносова, г. Архангельск
b Московский государственный университет имени М. В. Ломоносова, Московский центр фундаментальной и прикладной математики
Список литературы:
Аннотация: В работе получены асимптотические формулы при $x\to\infty$ для фундаментальной системы решений уравнения вида
$$ l(y):=(-1)^n(p(x)y^{(n)})^{(n)}+q(x)y=\lambda y,\qquad x\in[1,\infty), $$
где локально суммируемая функция $p$ допускает представление
$$ p(x)=(1+r(x))^{-1},\qquad r\in L^1 [1,\infty), $$
а $q$ – обобщенная функция, представимая при некотором фиксированном $k$, $0\le k\le n$, в виде $q=\sigma^{(k)}$, где
\begin{alignat*}{2} \sigma&\in L^1[1,\infty),&\qquad &\text{если }k<n, \\ |\sigma|(1+|r|)(1+|\sigma|)&\in L^1[1,\,\infty), &\qquad &\text{если }k=n. \end{alignat*}
Аналогичные результаты получены для уравнений $l(y)=\lambda y$, коэффициенты $p(x)$ и $q(x)$ которых допускают при некотором фиксированном $k$, $0\le k\le n$, представление
$$ p(x)=x^{2n+\nu}(1+r(x))^{-1},\qquad q=\sigma^{(k)},\quad \sigma(x)=x^{k+\nu}(\beta+s(x)), $$
где функции $r$ и $s$ удовлетворяют некоторым условиям интегрального убывания. Получены также теоремы об индексах дефекта минимального симметрического оператора, порожденного дифференциальным выражением $l(y)$ (при условии вещественности функций $p$ и $q$), и теоремы о спектрах соответствующих самосопряженных расширений.
Библиография: 19 названий.
Ключевые слова: дифференциальные операторы с коэффициентами-распределениями, квазипроизводные, асимптотика решений дифференциальных уравнений, дефектные числа дифференциального оператора.
Финансовая поддержка Номер гранта
Российский научный фонд 20-11-20261
Проведенные в этой статье исследования поддержаны Российским научным фондом, грант № 20-11-20261.
Поступило: 02.11.2022
Исправленный вариант: 16.11.2022
Англоязычная версия:
Mathematical Notes, 2023, Volume 113, Issue 2, Pages 228–242
DOI: https://doi.org/10.1134/S0001434623010261
Реферативные базы данных:
Тип публикации: Статья
УДК: 517.928

1. Формулировка результатов

Цель настоящей работы – получить асимптотические формулы при $x\to\infty$ для некоторой фундаментальной системы решений дифференциальных уравнений вида

$$ \begin{equation} l(y): =(-1)^n(p(x)y^{(n)})^{(n)}+q(x)y=\lambda y, \qquad x\in[1,\infty), \end{equation} \tag{1.1} $$
где $\lambda $ – комплексный параметр. Всюду далее предполагается, что комплекснозначная функция $p$ такая, что
$$ \begin{equation} p,\frac{1}{p}\in L^1_{\mathrm{loc}}[1,\infty), \end{equation} \tag{1.2} $$
где $L^1_{\mathrm{loc}}[1,\infty)$ – пространство локально суммируемых функций, $q$ – обобщенная функция, представимая в виде $q=\sigma^{(k)}$, где $k$ – некоторое целое число, $0\leqslant k\leqslant n$, а $\sigma$ – комплекснозначная функция такая, что
$$ \begin{equation} \begin{alignedat}{2} \sigma&\in L^1_{\mathrm{loc}}[1,\infty),&\qquad &\text{если}\quad k<n,\qquad \\ \frac{\sigma^2}{p}&\in L^1_{\mathrm{loc}}[1,\infty),&\qquad &\text{если}\quad k=n, \end{alignedat} \end{equation} \tag{1.3} $$
причем производная $k$-го порядка от $\sigma$ понимается в смысле теории распределений.

Более точные условия на функции $p$ и $q$ будут указаны в формулировках нижеследующих теорем. Как следствие асимптотических формул для решений рассматриваемых уравнений мы получим результаты об индексах дефекта дифференциальных операторов, порожденных выражениями $l(y)$ на полуоси $(1,\infty)$.

Новизна полученных теорем в сравнении с известными классическими результатами состоит в том, что мы допускаем, чтобы функция $q$ была обобщенной. Кроме того, результаты новы и для числовых функций $q$, для которых функция $|q(x)|$ растет на некоторой подпоследовательности интервалов при $x\to\infty$ сколь угодно быстро, но рост гасится быстрой осцилляцией так, что первообразная функции $q$ порядка $\leqslant n$ является суммируемой функцией.

В формулировках результатов данной работы используется понятие квазипроизводных заданной “допустимой” функции $y$, связанных с выражением $l(y)$, т.е. с функциями $p$, $\sigma$ и числом $k$. При фиксированном $j$, $j=0,1,\dots,2n-2$, для функции $y$ по приведенной ниже процедуре определяется квазипроизводная $y^{[j+1]}$, если квазипроизводная $y^{[j]}$ уже определена и является локально абсолютно непрерывной функцией ($y^{[j]}\in AC_{\mathrm{loc}}[1,\infty)$). При этом из условий (1.2) и (1.3) будет следовать, что $y^{[j+1]}\in L^1_{\mathrm{loc}}[1,\infty)$. Итак, зафиксируем целое число $k$, $0\leqslant k\leqslant n$, и определим квазипроизводную $y^{[j]}$, $j=0,1,\dots,2n-1$. Эту квазипроизводную при $j=0,1,\dots,n-1$ и при любом $k$ определим равенством

$$ \begin{equation} y^{[j]}:=y^{(j)}, \end{equation} \tag{1.4} $$
а квазипроизводную $y^{[n]}$ – равенством
$$ \begin{equation} y^{[n]}:=\begin{cases} py^{(n)}, &\text{если }k<n, \\ py^{(n)}+(-1)^n\sigma y, &\text{если }k=n. \end{cases} \end{equation} \tag{1.5} $$
Далее, при $k<n-1$ определим квазипроизводную $y^{[n+j]}$, $j=1,2,\dots,n-1$, положив
$$ \begin{equation} y^{[n+j]}:=\begin{cases} (y^{[n]})^{(j)}, &j=1,2,\dots,n-k-1, \\ (y^{[n+j-1]})'+(-1)^{j+k}C_k^{j+k-n}\sigma y^{[j+k-n]}, &j=n-k,\dots,n-1, \end{cases} \end{equation} \tag{1.6} $$
(здесь и далее $C_k^{j}$ – биномиальные коэффициенты), а при $k=n-1$ или $k=n$ определим эту квазипроизводную, положив
$$ \begin{equation} y^{[n+j]}:=(y^{[n+j-1]})'+(-1)^{j+k}C_k^{j+k-n}\sigma y^{[j+k-n]},\qquad j=1,2,\dots,n-1. \end{equation} \tag{1.7} $$

Сформулируем теперь основные результаты. Всюду далее через $o(1)$ обозначается функция, являющаяся бесконечно малой при $x\to\infty$.

Теорема 1. Пусть комплекснозначная локально интегрируемая функция $p$ допускает представление $p=(1+r)^{-1}$, где $r\in L^1[1,\infty)$, и $q$ – комплекснозначная обобщенная функция такова, что при некотором целом $k$, $0\leqslant k\leqslant n,$ она представима в виде $q=\sigma^{(k)}$, где $\sigma$ – регулярная обобщенная функция (т.е. принадлежит $L^1_{\mathrm{loc}}[1,\infty)$), а производная понимается в смысле теории распределений. Пусть при этом

$$ \begin{equation} \begin{alignedat}{2} \sigma&\in L^1[1,\infty),&\qquad &\textit{если}\quad k<n, \\ \sigma(1+|\sigma|)(1+|r|)&\in L^1[1,\infty),&\qquad &\textit{если}\quad k=n. \end{alignedat} \end{equation} \tag{1.8} $$

Тогда уравнение (1.1) при $\lambda\ne 0$ имеет фундаментальную систему решений $\{y_j\}_{j=1}^{2n}$ вида

$$ \begin{equation} y_j^{[s-1]}(x)=e^{z_jx}z_j^{s-1}(1+o(1)),\qquad s=1,2,\dots,2n, \end{equation} \tag{1.9} $$
где $z_j$ – различные корни степени $2n$ из $(-1)^n\lambda$.

Замечание 1. Полагая $s=1$ в формулах (1.9) и учитывая равенство (1.4), получаем, что уравнение (1.10) имеет фундаментальную систему решений вида

$$ \begin{equation*} y_j=e^{z_jx}(1+o(1)),\qquad j=1,2,\dots,2n, \end{equation*} \notag $$
и эти равенства допускают дифференцирование в классическом смысле до $(n-1)$-го порядка включительно. Более того, из формул (1.5)(1.7) следует, что это утверждение остается справедливым, если при $j\geqslant n$ вместо производных рассматриваются квазипроизводные.

Теорема 2. Пусть выполнены условия теоремы 1, но условия (1.8) заменены условиями

$$ \begin{equation*} \begin{alignedat}{2} x^{2n-1-k}\sigma&\in L^1[1,\infty),&\qquad &\textit{если}\quad k<n, \\ x^{n-1}\sigma(1+x^n|\sigma|)(1+|r|)&\in L^1[1,\infty),&\qquad &\textit{если}\quad k=n. \end{alignedat} \end{equation*} \notag $$
Тогда при $\lambda=0$ уравнение (1.1) имеет фундаментальную систему решений $\{y_j\}_{j=1}^{2n}$ вида
$$ \begin{equation*} y_j^{[s]}(x)=\begin{cases} \dfrac{x^{j-1-s}}{(j-1-s)!}(1+o(1)), &\textit{если }s=0,1,\dots,j-1, \\ x^{j-1-s}o(1), &\textit{если }s=j,j+1,\dots,2n-1. \end{cases} \end{equation*} \notag $$

Отметим, что в случае $k=0$ асимптотические формулы (1.9) для решений уравнения (1.1) хорошо известны и являются частным случаем более общего результата (см., например, [1; гл. VII, § 22, теорема 8] или [2; гл. 1, § 1.9, теорема 1.9.1]). В случае $k=0$ утверждение теоремы 2 также хорошо известно и представляет собой классический результат (см., например, [3; гл. X, § 17, теорема 17.1]). При $k> 1$ утверждения теорем 1 и 2 являются новыми.

Утверждения теорем 1 и 2 позволяют получить классическую асимптотику фундаментальной системы решений уравнения (1.1) для потенциалов $q$, которые могут расти сколь угодно быстро на некоторой последовательности интервалов, уходящих в бесконечность, при условии их подходящей осцилляции. Например, утверждения теорем 1 и 2 выполняются для коэффициентов вида $p(x)=1$ и $q=\sigma^{(k)}$, $1\leqslant k\leqslant n$, где

$$ \begin{equation*} \sigma(x)=\frac{\sin\varphi(x)}{1+x^2}\,. \end{equation*} \notag $$
Выбором функции $\varphi$ можно обеспечить произвольный рост функции $q$ на любой заданной последовательности непересекающихся отрезков.

Приведенные две теоремы дают асимптотику некоторой фундаментальной системы решений уравнений, которые можно рассматривать как возмущения канонического уравнения $y^{(2n)}-\lambda y=0$. Ниже приводим результат, где в качестве канонического уравнения рассматривается уравнение

$$ \begin{equation*} \alpha(x^{2n+\nu}y^{(n)})^{(n)}+\beta x^\nu y-\lambda y=0, \end{equation*} \notag $$
где $\alpha\ne 0$, $\beta$ – произвольные комплексные числа и $\nu\geqslant 0$. При $\nu=0$ это уравнение Эйлера.

Теорема 3. Пусть существуют числа $\nu \geqslant 0$ и целое число $k$, $0\leqslant k\leqslant n$, такие, что функция $p\in L^1_{\mathrm{loc}}[1,\infty)$ и обобщенная функция $q$, фигурирующие в уравнении (1.1), допускают представления

$$ \begin{equation} p(x)=\frac{\alpha x^{2n+\nu}}{1+r(x)}\,,\qquad q=\sigma^{(k)},\qquad \sigma(x)=x^{k+\nu}(\beta+s(x)), \end{equation} \tag{1.10} $$
где $\alpha\ne 0$ и $\beta$ – комплексные числа, а $r$ и $s$ – комплекснозначные функции. Пусть при этом функции $r$ и $s$ таковы, что
$$ \begin{equation} \frac{(\ln x)^{m-1}}{x}(1+|s(x)|)(|r(x)|+|s(x)|+|r(x)s(x)|)\in L^1[1,\infty), \end{equation} \tag{1.11} $$
где $m$ – наибольшее из чисел, равных кратностям корней многочлена от переменной $z$
$$ \begin{equation*} \begin{gathered} \, \mathscr F_{2n}(z,\nu,k)=\alpha\prod_{j=0}^{n-1} \biggl[\biggl(z+\frac{\nu}{2}\biggr)^2 -\biggl(\frac{\nu+1}{2}+j\biggr)^2\biggr] +(-1)^n\beta\prod_{j=1}^k(j+\nu)-\lambda(\nu), \\ \lambda(\nu)=\begin{cases} 0, &\textit{если }\nu>0, \\ (-1)^n\lambda, &\textit{если }\nu=0. \end{cases} \end{gathered} \end{equation*} \notag $$
Пусть $z_1$ – корень многочлена $\mathscr F_{2n}(z,\nu,k)$ кратности $l_1\leqslant m$. Тогда уравнение (1.1) имеет подсистему фундаментальных решений $\{y_j\}_{j=1}^{l_1}$ вида
$$ \begin{equation} \begin{alignedat}{2} y_j^{[s-1]}(x) &=x^{z_1-s+1/2}(\ln x)^{j-1}(c_s+o(1)), &\qquad &s=1,2,\dots,n, \\ y_j^{[n+s-1]}(x) &=x^{z_1+n+\nu-s+1/2}(\ln x)^{j-1}(c_{n+s}+o(1)), &\qquad &s=1,2,\dots,n, \end{alignedat} \end{equation} \tag{1.12} $$
где $c_s$ – некоторые комплексные числа и $c_1\ne 0$. Такую же асимптотику имеет и другая подсистема фундаментальных решений $\{y_j\}$, $j=l_1+1,\dots,l_1+l_2$, отвечающая корню $z_2$ многочлена $\mathscr F_{2n}(z,\nu,k)$ кратности $l_2$ и т.д.

Замечание 2. Многочлен $\mathscr F_{2n}(z,\nu,k)$ является характеристическим многочленом некоторой постоянной матрицы размерности $2n$. Структура этой матрицы такова, что геометрическая кратность любого собственного значения этой матрицы равна $1$, причем первая компонента соответствующего собственного вектора отлична от нуля. Числа $c_s$ в формулировке теоремы 3 являются координатами собственного вектора, соответствующего собственному значению $z_1$. Это будет ясно из приведенного ниже доказательства теоремы 3 в п. 5.

Как и утверждения теорем 1 и 2, утверждение теоремы 3 также справедливо для некоторых сильно осциллирующих функций $q(x)$. Это, очевидно, можно обеспечить подходящим подбором функций $r(x)$ и $s(x)$.

Аналог теоремы 3 для дифференциальных уравнений вида

$$ \begin{equation*} l_{2n}(y):=\sum^n_{j=0}(p_{j}(x)y^{(j)})^{(j)}=\lambda y \end{equation*} \notag $$
с вещественными аналитическими в некоторой окрестности бесконечности коэффициентами
$$ \begin{equation*} p_{j}(z):=z^{2j+\nu}\biggl[a_j+\sum_{i=1}^{+\infty}a_i^{(j)}z^{-i}\biggr],\quad j=0,1,\dots,n,\qquad |z|\geqslant x_0>1, \end{equation*} \notag $$
где $a_n\ne 0$, а $\nu\geqslant 0$ – целое число, впервые был получен Орловым в работе [4], и тем самым был найден класс линейных симметрических дифференциальных операторов с вещественными аналитическими коэффициентами, дефектные числа которых определяются как число корней полиномов $\mathscr F_{2n}(z,\nu)$, лежащих в левой полуплоскости (эти корни конструируются явно по числам $a_0,a_1,\dots,a_n$).

Позднее Пэрис и Вуд в работе [5] заново открыли теорему Орлова для частного случая $p_j(x)=a_jx^{2j+\nu}$, $j=0,1,\dots,n$, где $\nu$ – неотрицательное целое число, и подробно изучли полином $\mathscr F_{2n}(z,\nu)$. В частности, они установили, что существуют положительные числа $K$ и $\nu$ такие, что индекс дефекта минимального замкнутого положительного симметрического оператора, порожденного выражением

$$ \begin{equation*} l_6(y)=-(x^{6+\nu}y^{(3)})^{(3)}+Kx^\nu y \end{equation*} \notag $$
равен $(5,5)$, тем самым уточнив результат Кауффмана [6] о том, что для этого оператора индекс дефекта не равен $(3,3)$. Позже результаты работы [5] вошли в книгу [7], в которой также были рассмотрены дифференциальные операторы, порожденные выражениями нечетного порядка частного вида, и исследованы соответствующие им полиномы.

В работе Мирзоева [8] исследуется задача об индексе дефекта минимального замкнутого симметрического оператора, порожденного симметрическим дифференциальным выражением произвольного (четного или нечетного) порядка. Изучается также характер спектра соответствующих самосопряженных расширений. Полученные в [8] результаты обобщают результаты работы [4]. Теорема 3 при $k=0$ является частным случаем утверждения из работы [8].

Отметим, что формулировки теорем 13 анонсированы ранее в работе авторов [9]. Но их доказательства проведены только в случае $n=1$. Ознакомление с ними в этом частном случае поможет читателю понять сущность дела в общем случае.

Работа [9] послужила началом исследований по нахождению асимптотик решений уравнений на полуоси с коэффициентами-распределениями при $x\to\infty$. Основная идея состояла в конструкции различных регуляризаций соответствующих уравнений с помощью приведения этих уравнений к эквивалентным системам первого порядка. Впоследствии в работе [10] была проведена регуляризация уравнения

$$ \begin{equation*} y^{(n)}+(a_1+p_1(x))y^{(n-1)} +(a_2+p_2'(x))y^{(n-2)}+\dotsb+(a_n+p_n'(x))y=\lambda y, \end{equation*} \notag $$
где $a_1,a_2,\dots,a_n,\lambda$ – комплексные числа,
$$ \begin{equation*} |p_1|+(1+|p_2-p_1|)\sum_{j=2}^n|p_j|\in L_{\mathrm{loc}}^1[1,\infty ), \end{equation*} \notag $$
а производные функций $p_j$ понимаются в смысле теории распределений. Было установлено, что при некоторых условиях убывания функций $p_j$ на бесконечности главный член асимптотики при $x\to\infty$ фундаментальной системы решений этого уравнения определяется, как и в классическом случае, по корням многочлена
$$ \begin{equation*} Q(z)=z^n+a_1z^{n-1}+\dotsb+a_n-\lambda. \end{equation*} \notag $$
В работе [11] анонсированы аналоги теорем 13 для дифференциального уравнения нечетного порядка вида
$$ \begin{equation*} i^{2n+1}\{(qy^{(n+1)})^{(n)}+(qy^{(n)})^{(n+1)}\}+py=\lambda y. \end{equation*} \notag $$
Общий подход к регуляризации дифференциальных уравнений высших порядков с коэффициентами-распределениями был предложен в работе Мирзоева и Шкаликова [12]. Но вопросы об асимптотиках решений таких уравнений при $x\to\infty$ в [12] не изучались.

Коротко скажем о структуре работы. В п. 2 мы проводим регуляризацию выражения $l(y)$, в частности, обосновываем целесообразность определения квазипроизводных допустимых функций формулами (1.4)(1.7). В пп. 3, 4 и 5 приведены доказательства теорем 1, 2 и 3 соответственно. Пункт 6 посвящен приложениям этих теорем к исследованию индексов дефекта соответствующих симметрических операторов (с вещественными функциями $p$ и $\sigma$) и изучению спектров самосопряженных расширений этих операторов.

2. Построение матрицы, согласованной с выражением $l(y)$. Квазипроизводные и квазидифференциальное выражение

Зафиксируем целое число $k$, $0\leqslant k\leqslant n$, и определим квадратную матрицу $\mathbf{F}_{2n}^k$ с элементы $f_{ij}^k(=:f_{ij})$, $i,j=1,2,\dots,2n$, размерности $2n$, полагая

  • i) при любом $k$
    $$ \begin{equation*} f_{j,j+1}:=1,\quad j=1,2,\dots,n-1,n+1,\dots,2n-1,\qquad f_{n,n+1}:=p^{-1}; \end{equation*} \notag $$
  • ii) при $k<n$
    $$ \begin{equation*} f_{2n-1+j-k,j}:=(-1)^{n+j}C_k^{j-1}\sigma,\qquad j=1,2,\dots,k+1, \end{equation*} \notag $$
    а при $k=n$
    $$ \begin{equation*} \begin{gathered} \, f_{n,1}:=(-1)^{n+1}\sigma p^{-1},\qquad f_{2n,n+1}:=-\sigma p^{-1},\qquad f_{2n,1}:=(-1)^n\sigma^2 p^{-1}, \\ f_{n-1+j,j}:=(-1)^{n+j}C_n^{j-1}\sigma,\quad j=2,3,\dots,n. \end{gathered} \end{equation*} \notag $$
Все неуказанные в этом списке элементы матрицы $\mathbf{F}_{2n}^k$ равны нулю.

Таким образом, матрица $\mathbf{F}_{2n}^k$ при $k<n$ имеет вид

$$ \begin{equation*} \mathbf{F}_{2n}^k=\begin{pmatrix} 0 &1 & & & & & &0 \\ & &\ddots & & & &\cdot^{\displaystyle{\displaystyle\cdot}^{\displaystyle\cdot}} & \\ & & &1 & &0 & & \\ & & & &p^{-1} & & & \\ & & & & &1 & & \\ & &D_k & & & &\ddots & \\ & & & & & & &1 \\ & & & & & & &0 \end{pmatrix}, \end{equation*} \notag $$
где элемент $p^{-1}$ находится на пересечении $n$-й строки и $(n+1)$-го столбца, диагональная матрица $D_k=\operatorname{diag}(f_{2n-k,1},f_{2n+1-k,2},\dots,f_{2n,k+1})$ размерности $k+1$ (см. равенство ii) при $k<n$) занимает левый нижний блок, а при $k=n$
$$ \begin{equation*} \mathbf{F}_{2n}^n=\begin{pmatrix} 0 &1 & & & & & &0 \\ & &\ddots & & & &\cdot^{\displaystyle{\displaystyle\cdot}^{\displaystyle\cdot}} & \\ & & &1 & &0 & & \\ (-1)^{n+1}\sigma p^{-1} & & & &p^{-1} & & \\ & & & & &1 & & \\ & &D_n & & & &\ddots & \\ & & & & & & &1 \\ (-1)^n\sigma^2p^{-1} & & & &-\sigma p^{-1} & & &0 \end{pmatrix}, \end{equation*} \notag $$
где диагональная матрица $D_n=\operatorname{diag}(f_{n+1,2},f_{n+2,3},\dots,f_{2n-1,n})$ размерности $n-1$ (см. равенства ii) при $k=n$) заполняет блок $\{f_{ij}\}$ с номерами строк $i=n+1,n+2, \dots,2n-1$ и столбцов $j=2,3,\dots,n$ матрицы $\mathbf{F}_{2n}^n$. На всех незаполненных местах матриц $\mathbf{F}_{2n}^k$ находятся нули.

Пусть $g\in L^1_{\mathrm{loc}}[1,\infty)$. Символом $\mathbf g$ обозначим вектор-столбец $\mathbf g:=(0,0,\dots,0,g)^T$ ($T$ – знак транспонирования) размерности $2n$, и рассмотрим систему линейных дифференциальных уравнений первого порядка

$$ \begin{equation} \mathbf y'=\mathbf{F}_{2n}^k\mathbf y+\mathbf g, \end{equation} \tag{2.1} $$
где $\mathbf y$ – $2n$-компонентная неизвестная вектор-функция. Из условий (1.2) и (1.3) легко извлечь, что все элементы матрицы $\mathbf{F}_{2n}^k$ являются локально интегрируемыми на $[1,\infty)$, а координаты $\mathbf g$ также являются таковыми. Таким образом, выполняются условия, обеспечивающие справедливость теоремы существования и единственности решения задачи Коши для уравнения (2.1), поставленной в произвольной точке полуоси $[1,\infty)$ (см., например, [1; гл. V, § 16, теорема 1]). В частности, компоненты любого решения $\mathbf y$ системы (2.1), которые в дальнейшем мы будем обозначать символами $y^{[0]},y^{[1]},\dots,y^{[2n-1]}$, принадлежат пространству $AC_{\mathrm{loc}}[1,\infty)$. С другой стороны, структура матрицы $\mathbf{F}_{2n}^k$ позволяет последовательно выразить компоненты $y^{[1]},\dots,y^{[2n-1]}$ вектор-столбца $\mathbf y$ через его первую компоненту $y^{[0]}$, а именно, из первых $2n-1$ уравнений системы (2.1) находим, что они удовлетворяют равенствам
$$ \begin{equation} y^{[j+1]}=f_{j+1,j+2}^{-1}(y^{[j]})',\qquad j=0,1,\dots,2n-2-k, \end{equation} \tag{2.2} $$
и
$$ \begin{equation} y^{[2n-k+j]}=f_{2n-k+j,2n-k+j+1}^{-1} \{(y^{[2n-k+j-1]})'-f_{2n-k+j,j+1}y^{[j]}\},\qquad j=0,1,\dots,k-1. \end{equation} \tag{2.3} $$

Таким образом, эти формулы позволяют выразить компоненту $y^{[j+1]}$ при фиксированном $j\in\{0,1,\dots,2n-2\}$ через $y^{[s]}$, где $s\leqslant j$. Следовательно, все компоненты $y^{[j]}$ в итоге выражаются через $y^{[0]}:=y$. Функция $y^{[j]}$, $j=0,1,\dots,2n-1$, определяемая рекуррентными формулами (2.2) и (2.3), называется квазипроизводной порядка $j$ допустимой функции $y$, ассоциированной с матрицей $\mathbf{F}_{2n}^k$. Несложные вычисления с учетом равенств i) и ii) показывают, что для них справедливы формулы (1.4)(1.7).

Последнее уравнение системы (2.1) в принятых нами обозначениях записывается в виде

$$ \begin{equation} \tau(y):=(y^{[2n-1]})'-f_{2n,k+1}y^{[k]}-f_{2n,1}y=g. \end{equation} \tag{2.4} $$

Выражение $\tau(y)$ называется квазидифференциальным выражением, ассоциированным с матрицей $\mathbf{F}_{2n}^k$. В терминах элементов этой матрицы выражение $\tau(y)$ при $k<n$ принимает вид

$$ \begin{equation} \tau(y)=(y^{[2n-1]})'+(-1)^{n+k}\sigma y^{[k]}, \end{equation} \tag{2.5} $$
а при $k=n$
$$ \begin{equation} \tau (y)=(y^{[2n-1]})'+\sigma p^{-1}y^{[n]}+(-1)^{n+1}\sigma^2 p^{-1}y. \end{equation} \tag{2.6} $$

Естественная область определения $\mathscr D(\tau)$ выражения $\tau$ – это множество всех комплекснозначных функций $y$ таких, что $y^{[j]}\in AC_{\mathrm{loc}}[1,\,\infty)$, $j=0,1,\dots,2n-1$. Отметим, что если $y\in\mathscr{D}(\tau)$, то $\tau(y)\in L^1_{\mathrm{loc}}[1,\,\infty)$.

Квазидифференциальное уравнение (2.4) и система уравнений (2.1) равносильны в том смысле, что если все квазипроизводные функции $y$ до $(2n-1)$-го порядка принадлежат классу $AC_{\mathrm{loc}}[1,\infty)$ и $y$ удовлетворяет уравнению (2.4) почти всюду, то $\mathbf y$ является решением системы (2.1). Наоборот, если вектор-функция $\mathbf y$ является решением системы (2.1), то ее координаты принадлежат классу $AC_{\mathrm{loc}}[1,\infty)$, а первая координата $y=y^{[0]}$ удовлетворяет уравнению (2.4) почти всюду. Таким образом, условия (1.2) и (1.3) обеспечивают справедливость теоремы существования и единственности решения задачи Коши для уравнения (2.4), поставленной в произвольной точке полуоси $[1,\infty)$.

Лемма 1. Пусть $y\in\mathscr D(\tau)$. Тогда распределение $l(y)$, определенное в (1.1), является регулярной обобщенной функцией и для него справедливо равенство

$$ \begin{equation*} l(y)=(-1)^n\tau(y). \end{equation*} \notag $$

Доказательство. При $k=0$ лемма 1, очевидно, справедлива.

Рассмотрим случай $k=1$, и пусть $y\in\mathscr D(\tau)$. При $n=1$, согласно формулам (1.5) и (2.6), квазипроизводная $y^{[1]}$ и квазидифференциальное выражение $\tau(y)$ имеют вид

$$ \begin{equation*} y^{[1]}=py'-\sigma y,\qquad \tau(y)=(y^{[1]})'+\sigma p^{-1}y^{[1]}+\sigma^2 p^{-1}y. \end{equation*} \notag $$
Из условий (1.2) и (1.3) следует, что в выражении $\tau(y)$ все скобки можно раскрыть, и в результате получим справедливость леммы 1 (см. работы [12]–[16]).

При $n\geqslant 2$ квазипроизводные и квазидифференциальное выражение (см. формулы (1.4)(1.6) и (2.5)) принимают вид

$$ \begin{equation*} \begin{gathered} \, y^{[j]} =y^{(j)},\quad j=0,1,\dots,n-1, \qquad y^{[n+j]} =(py^{(n)})^{(j)},\quad j=0,1,\dots,n-2, \\ y^{[2n-1]} =(py^{(n)})^{(n-1)}+(-1)^n\sigma y, \qquad \tau y =(y^{[2n-1]})'+(-1)^{n+1}\sigma y', \end{gathered} \end{equation*} \notag $$
и справедливость леммы 1 также очевидна.

Рассмотрим теперь случай $k\geqslant 2$. Применяя индукцию по $r$, покажем, что при $r=1,2,\dots,k-1$, для любого $n\geqslant 2$ справедлива формула

$$ \begin{equation} {y^{[2n-k+r]} =(py^{(n)})^{(n-k+r)}+(-1)^n\sigma^{(r)}y +(-1)^n\sum_{m=1}^r (-1)^mC_{k+m-r-1}^m\sigma^{(r-m)}y^{(m)},} \end{equation} \tag{2.7} $$
где производные всюду понимаются в смысле теории распределений.

Покажем справедливость формулы (2.7) при $r=1$. Из формул (1.6) (при $k<n-1$, полагаем $j=n-k+1$) и (1.7) (при $k=n-1$, полагаем $j=2$ и при $k=n$, полагаем $j=1$) получаем, что справедливо равенство

$$ \begin{equation*} y^{[2n-k+1]}=(y^{[2n-k]})'+(-1)^{n+1}C_k^{1}\sigma y^{[1]}. \end{equation*} \notag $$
Далее, при $k=n$ воспользуемся равенством (1.5). При $k=n-1$ положим $j=1$ в равенстве (1.7) и воспользуемся равенством (1.5). При $k<n-1$ положим сначала $j=n-k$, а затем $j=n-k-1$ в равенстве (1.6). Вновь используя (1.5), находим
$$ \begin{equation*} y^{[2n-k]}=(py^{(n)})^{(n-k)}+(-1)^n\sigma y. \end{equation*} \notag $$
Используя это равенство и равенство $y^{[1]}=y'$ (см. (1.4)), получаем
$$ \begin{equation*} y^{[2n-k+1]}=(py^{(n)})^{(n-k+1)}+(-1)^n\sigma'y-(-1)^nC_{k-1}^1\sigma y', \end{equation*} \notag $$
т.е. формула (2.7) справедлива при $r=1$.

Предположим теперь, что формула (2.7) справедлива при $r=s$, где все еще $s<k-1$, и докажем, что эта формула справедлива при $r=s+1$. Действительно, используя формулы (1.6) (при $k<n-1$, полагаем $j=n-k+s+1$) и (1.7) (при $k=n-1$ полагаем $j=s+2$, и при $k=n$ полагаем $j=s+1$), а также формулу (2.7) при $r=s$ получаем, что

$$ \begin{equation*} \begin{aligned} \, I &:=(-1)^n\{y^{[2n-k+s+1]}-(py^{(n)})^{(n-k+s+1)} -(-1)^n\sigma^{(s+1)}y\} \\ &\,\,=\sigma^{(s)}y'+\sum_{m=1}^s(-1)^mC_{k+m-s-1}^m [\sigma^{(s+1-m)}y^{(m)}+\sigma^{(s-m)}y^{(m+1)}] \\ &\,\,\qquad +(-1)^{s+1}C_k^{s+1}\sigma y^{[s+1]}. \end{aligned} \end{equation*} \notag $$
Таким образом,
$$ \begin{equation*} \begin{aligned} \, I &=\sum_{m=1}^{s+1}(-1)^m(C_{k+m-s-1}^m-C_{k+m-s-2}^{m-1}) \sigma^{(s+1-m)}y^{(m)} \\ &=\sum_{m=1}^{s+1}(-1)^mC_{k+m-s-2}^m\sigma^{(s+1-m)}y^{(m)}, \end{aligned} \end{equation*} \notag $$
так как при $t=k+m-s-1$ справедливо равенство
$$ \begin{equation*} C_t^m-C_{t-1}^{m-1}=C_{t-1}^m. \end{equation*} \notag $$
Следовательно, формула (2.7) справедлива при всех значениях $r\leqslant k-1$. В частности, при $r=k-1$ эта формула принимает вид
$$ \begin{equation*} y^{[2n-1]}=(py^{(n)})^{(n-1)}+(-1)^n\sigma^{(k-1)}y +(-1)^n\sum_{m=1}^{k-1}(-1)^m\sigma^{(k-1-m)} y^{(m)}. \end{equation*} \notag $$

При $k\leqslant n-1$, используя формулу (2.6) и полученную выше формулу, находим, что

$$ \begin{equation*} \begin{aligned} \, &(-1)^n\tau(y)-(-1)^n(py^{(n)})^{(n)}-\sigma^{(k)}y \\ &\qquad=\sigma^{(k-1)}y'+\sum_{m=1}^{k-1}(-1)^m \{\sigma^{(k-m)}y^{(m)}+\sigma^{(k-m-1)}y^{(m+1)}\} +(-1)^k\sigma y^{(k)}=0, \end{aligned} \end{equation*} \notag $$
т.е. в случае $k\leqslant n-1$ утверждение леммы 1 также справедливо.

При $k=n$, используя формулу (2.6) и опять полученную выше формулу, находим

$$ \begin{equation*} \begin{aligned} \, &(-1)^n\tau y-(-1)^n(py^{(n)})^{(n)}-\sigma^{(n)}y \\ &\quad=\sigma^{(n-1)}y'+\sum_{m=1}^{n-1}(-1)^m \{\sigma^{(n-m)}y^{(m)}+\sigma^{(n-m-1)}y^{(m+1)}\} +(-1)^n\sigma p^{-1}y^{[n]}-\sigma^2 p^{-1}y \\ &\quad=0, \end{aligned} \end{equation*} \notag $$
что доказывает справедливость леммы 1 в случае $k=n$.

3. Доказательство теоремы 1

Исходя из леммы 1, дифференциальное выражение $l(y)$ отождествим с квазидифференциальным выражением $(-1)^n\tau(y)$ и перепишем дифференциальное уравнение (1.1) в виде

$$ \begin{equation} \tau(y)=(-1)^n\lambda y. \end{equation} \tag{3.1} $$
Это уравнение совпадает с уравнением (2.4), где $g=(-1)^n\lambda y$, а уравнение (2.4), как отмечено выше, равносильно системе (2.1) с той же матрицей $\mathbf{F}_{2n}^k$ и соответствующим вектором $\mathbf g:=(0,0,\dots,0,(-1)^n\lambda y)^T$. Полученную при этом систему запишем в виде
$$ \begin{equation} \mathbf{y}'=(\mathbf\Lambda+\mathbf R(x))\mathbf y. \end{equation} \tag{3.2} $$
Здесь $\mathbf y$ – вектор-столбец c элементами $\mathbf y:=(y,y^{[1]},\dots,y^{[2n-1]})^T$, элементы $\lambda_{ij}$ квадратной матрицы $\boldsymbol\Lambda $ размерности $2n$ определяются равенствами $\lambda_{2n,1}:=(-1)^n\lambda$, $\lambda_{j,j+1}:=1$ при $j=1,2,\dots,2n-1$ и $\lambda_{ij}:=0$ для всех остальных значений $i$ и $j$, т.е.
$$ \begin{equation*} \boldsymbol\Lambda=\begin{pmatrix} 0 &1 & & & & & &0 \\ & &\ddots & & & &\cdot^{\displaystyle{\displaystyle\cdot}^{\displaystyle\cdot}} & & \\ & & &1 & &0 & & \\ & & & &1 & & & \\ & & & & &1 & & \\ & & & & & &\ddots & \\ & & & & & & &1 \\ (-1)^n\lambda & & & & & & &0 \end{pmatrix}, \end{equation*} \notag $$
а элементы $r_{ij}$ матрицы $\mathbf R(x)$ совпадают с соответствующими элементами матрицы $\mathbf{F}_{2n}^k$, за исключением элементов вида $r_{j,j+1}$, $j=1,2,\dots,2n-1$, которые определены равенствами
$$ \begin{equation*} r_{j,j+1}=0\qquad \text{при}\quad j=1,\dots,n-1,n+1,\dots,2n-1,\quad r_{n,n+1}=r. \end{equation*} \notag $$
Таким образом, уравнение (3.1), а следовательно, и уравнение (1.1), в указанном выше смысле равносильно системе (3.2).

Если $\lambda\ne 0$, то матрица $\mathbf\Lambda$ имеет различные собственные числа

$$ \begin{equation*} z_j=\alpha\omega_j,\qquad j=1,2,\dots,2n, \end{equation*} \notag $$
где $\alpha$ – фиксированный корень степени $2n$ из числа $(-1)^n\lambda$, а $\omega_1,\omega_2,\dots,\omega_{2n}$ – все корни $2n$-й степени из $1$ (т.е. корни уравнения $\omega^{2n}=1$).

Из условий (1.8) теоремы 1 следует, что все элементы матрицы $\mathbf R(x)$ принадлежат пространству $L^1[1,\infty)$. Тогда, применяя утверждение задачи 29 из книги [17; гл. III] (см. также [2; теорема 1.8.1]), получаем, что система (3.2) имеет решения $\mathbf y_j$, представимые в виде

$$ \begin{equation*} \mathbf y_j=e^{z_jx}\mathbf p_j(1+o(1)),\qquad j=1,2,\dots,2n, \end{equation*} \notag $$
где $\mathbf p_j=(1,z_j,z_j^2,\dots,z_j^{2n-1})^T$ – собственный вектор матрицы $\mathbf\Lambda$, соответствующий собственному числу $z_j$.

Учитывая далее связь между решениями уравнения (3.1) и решениями системы (3.2), получаем, что уравнение (3.1) (а тогда и (1.1)) имеет фундаментальную систему решений $y_j$, $j=1,2,\dots,2n$, для которой при $x\to+\infty$ справедливы асимптотические формулы (1.9). Теорема доказана.

4. Доказательство теоремы 2

Мы уже показали, что уравнения (1.1) и (3.1) эквивалентны системе (3.2). При $\lambda=0$ матрица $\mathbf\Lambda$ в (3.2) имеет жорданову форму, содержащую одну жорданову клетку, и ей отвечает $2n$-кратное собственное значение $0$. Заметим, что в рассматриваемом случае система уравнений $\mathbf y'=\boldsymbol\Lambda\mathbf y$ имеет $2n$ линейно независимых решений

$$ \begin{equation*} \mathbf y_j=\biggl(\frac{x^{j-1}}{(j-1)!}\,, \frac{x^{j-2}}{(j-2)!}\,,\dots,1,0,0\dots,0\biggr)^T,\qquad j =1,2,\dots,2n. \end{equation*} \notag $$

Обозначим символом $\mathbf Q$ диагональную матрицу вида

$$ \begin{equation*} \mathbf Q:=\operatorname{diag}(1,x,x^2,\dots,x^{2n-2},x^{2n-1}), \end{equation*} \notag $$
и положим $\widetilde{\mathbf R}:=\mathbf Q\mathbf R\mathbf Q^{-1}$, где матрица $\mathbf R$ определена выше. Элементы $\widetilde r_{ij}$, $i,j=1,2,\dots,2n$, матрицы $\widetilde{\mathbf R}$ определяются следующим образом:

если $k<n$, то

$$ \begin{equation*} \widetilde r_{2n-1+j-k,j}=(-1)^{n+j}C_k^{j-1}x^{2n-k-1}\sigma,\quad j=1,2,\dots,k+1,\qquad \widetilde r_{n,n+1}=x^{-1}r; \end{equation*} \notag $$

если $k=n$, то

$$ \begin{equation*} \begin{gathered} \, \widetilde r_{n,1}=(-1)^{n+1}x^{n-1}\sigma(1+r),\qquad \widetilde r_{n,n+1}=x^{-1}r,\qquad \widetilde r_{2n,n+1}=-x^{n-1}\sigma(1+r), \\ \widetilde r_{n-1+j,j}=(-1)^{n+j}C_n^{j-1}x^{n-1}\sigma,\quad j=2,\dots,n,\qquad \widetilde r_{2n,1}=(-1)^nx^{2n-1}\sigma^2(1+r). \end{gathered} \end{equation*} \notag $$
Все неуказанные элементы матрицы $\widetilde{\mathbf R}$ – нули.

Согласно теореме 1.10.1 из главы 1, § 1.10 книги [2] уравнение (3.2) имеет решения $\mathbf y_j$, $j=1,2,\dots,2n$, вида

$$ \begin{equation} \begin{aligned} \, \nonumber \mathbf y_{j}=\biggl(&\frac{x^{j-1}}{(j-1)!}+o(x^{j-1}), \frac{x^{j-2}}{(j-2)!}+o(x^{j-2}),\dots, \\ &\qquad 1+o(1),o(x^{-1}),o(x^{-2}),\dots,o(x^{-2n+j})\biggr)^T \end{aligned} \end{equation} \tag{4.1} $$
при условии, что все элементы матрицы $\widetilde{\mathbf R}$ принадлежат пространству $L^1[1,\infty)$. Но это условие гарантируется условиями теоремы 2. Поэтому из представлений (4.1) получаем утверждение теоремы. Этим заканчивается доказательство.

5. Доказательство теоремы 3

Запишем уравнение (3.1) в виде системы

$$ \begin{equation} \mathbf y'=(\mathbf{F}_{2n}^k+\boldsymbol\Lambda_0)\mathbf y, \end{equation} \tag{5.1} $$
где матрица $\mathbf{F}_{2n}^k$ определена ранее, а элементы квадратной матрицы $\boldsymbol\Lambda_0=(\lambda_{ij}^0)$ размерности $2n$ определяются равенствами $\lambda_{2n,1}^0:=(-1)^n\lambda$ и $\lambda_{ij}^0:=0$ для всех остальных значений $i$ и $j$.

Обозначив символом $\mathbf D$ диагональную матрицу с элементами

$$ \begin{equation*} d_i(x):=x^{-i+1/2},\quad d_{i+n}(x):=x^{n+\nu-i+1/2},\qquad i=1,\dots,n, \end{equation*} \notag $$
сделаем замену $\mathbf y=\mathbf D\mathbf u$ и перепишем уравнение (6.4) в виде
$$ \begin{equation} \mathbf u'=\mathbf A(x)\mathbf u, \end{equation} \tag{5.2} $$
где
$$ \begin{equation*} \mathbf A(x)=\mathbf D^{-1}(\mathbf{F}_{2n}^k+\mathbf\Lambda_0)\mathbf D -\mathbf D^{-1}\mathbf D'. \end{equation*} \notag $$
С учетом представлений (1.1) функций $p$ и $\sigma$ находим, что
$$ \begin{equation*} \mathbf A(x)=x^{-1}(\mathbf A_0+\mathbf B(x)). \end{equation*} \notag $$
Простые вычисления показывают, что элементы $a_{ij}$ числовой матрицы $\mathbf A_0$ и элементы $b_{ij}(x)$ матрицы-функции $\mathbf B(x)$ определяются следующими равенствами:
  • i) диагональные и первые наддиагональные элементы этих матриц при любом $0\leqslant k\leqslant n$ определяются равенствами
    $$ \begin{equation*} \begin{gathered} \, a_{i,i+1}=1,\qquad i=1,\dots,n-1,n+1,\dots,{2n-1}, \\ a_{n,n+1}=\frac{1}{\alpha}\,,\qquad b_{n,n+1}(x)=\frac{1}{\alpha}\,r(x), \qquad a_{i,i}=i-\frac{1}{2}\,,\quad i=1,2\dots,n, \\ a_{i+1+n,i+1+n}=i+\frac{1}{2}-n-\nu,\quad i=0, 1,2\dots,n-1; \end{gathered} \end{equation*} \notag $$
  • ii) при $k=0$
    $$ \begin{equation*} \begin{alignedat}{2} a_{2n,1} &=(-1)^{n+1}\beta+(-1)^n\lambda,\quad b_{2n,1}(x)=(-1)^{n+1}s(x),&\qquad &\text{если}\quad\nu=0, \\ a_{2n,1} &=(-1)^{n+1}\beta,\quad b_{2n,1}(x)=(-1)^{n+1}s(x)+\frac{(-1)^n\lambda}{x^\nu}\,,&\qquad &\text{если}\quad\nu>0; \end{alignedat} \end{equation*} \notag $$
    при $0<k<n$
    $$ \begin{equation*} \begin{aligned} \, a_{2n-1+i-k,i}&=(-1)^{n+i}C_k^{i-1}\beta, \\ b_{2n-1+i-k,i}(x)&=(-1)^{n+i}C_k^{i-1}s(x), \end{aligned}\qquad i=1,2,\dots,k+1, \end{equation*} \notag $$
    а также
    $$ \begin{equation*} \begin{alignedat}{2} a_{2n,1} &=(-1)^n\lambda,\quad b_{2n,1}(x)=0, &\qquad &\text{если}\quad\nu=0, \\ a_{2n,1} &= 0,\quad b_{2n,1}(x)=\frac{(-1)^n\lambda}{x^\nu}\,, &\qquad &\text{если}\quad\nu>0; \end{alignedat} \end{equation*} \notag $$
    при $k=n$
    $$ \begin{equation*} \begin{gathered} \, a_{n+i-1,i} =(-1)^{n+i}C_n^{i-1}\beta, \quad b_{n+i-1,i}(x) =(-1)^{n+i}C_n^{i-1}s(x),\qquad i=2,3, \dots,n, \\ a_{n,1} =(-1)^{n+1}\frac{\beta}{\alpha}\,, \qquad b_{n,1} =(-1)^{n+1}\frac{1}{\alpha}(s(x)+\beta r(x)+r(x)s(x)), \end{gathered} \end{equation*} \notag $$
    а также
    $$ \begin{equation*} \begin{gathered} \, a_{2n,n+1} =-\frac{\beta}{\alpha}\,,\qquad b_{2n,n+1}=-\frac{1}{\alpha}(s(x)+\beta r(x)+r(x)s(x)), \\ \begin{aligned} \, a_{2n,1} &=(-1)^n\frac{1}{\alpha}\,\beta^2+(-1)^n\lambda, \\ b_{2n,1}(x) &=(-1)^n\frac{1}{\alpha}(r(x)(\beta+s(x))^2+(2\beta+s(x))s(x)), \end{aligned}\qquad \text{если}\quad \nu=0, \\ \begin{aligned} \, a_{2n,1} &=(-1)^n\frac{1}{\alpha}\,\beta^2, \\ b_{2n,1}(x) &=(-1)^n\frac{1}{\alpha}(r(x)(\beta+s(x))^2+(2\beta+s(x))s(x)) +\frac{(-1)^n\lambda}{x^\nu}\,, \end{aligned}\quad \text{если}\quad \nu>0. \end{gathered} \end{equation*} \notag $$
Все неуказанные элементы в матрицах $\mathbf A_0$ и $\mathbf B(x)$ равны нулю.

Сделав замену $x=e^t$ в системе (5.2), получим систему уравнений

$$ \begin{equation} \mathbf v'(t)=(\mathbf A_0+\mathbf R(t))\mathbf v(t), \end{equation} \tag{5.3} $$
где $\mathbf v(t)=\mathbf u(e^t)$, $\mathbf R(t)=\mathbf B(e^t)$. Вычисления показывают, что характеристический многочлен матрицы $\mathbf A_0$ имеет вид (здесь $\mathbf 1$ – единичная матрица)
$$ \begin{equation*} \det(\mathbf A_0-z\boldsymbol 1) =\frac{1}{\alpha}\,\mathscr F_{2n}(z,\nu,k). \end{equation*} \notag $$

Структура матрицы $\mathbf A_0$ такова, что $a_{j,j+1}\ne 0$ при $1\leqslant j\leqslant 2n-1$ и $a_{jl}=0$ при $2\leqslant j+1<l\leqslant 2n$ (см. формулы для них, приведенные выше). Поэтому собственный вектор матрицы $\mathbf A_0$, соответствующий какому-либо собственному значению, однозначно определяется заданием своей первой координаты, которая должна быть отлична от нуля. Таким образом, каждому собственному значению матрицы $\mathbf A_0$ соответствует только один собственный вектор, т.е. геометрическая кратность любого собственного значения матрицы $\mathbf A_0$ равна единице. Иными словами, каждому собственному значению матрицы $\mathbf A_0$ соответствует только одна жорданова клетка в ее канонической форме. Поэтому максимум размерности жордановых клеток в канонической форме матрицы $\mathbf A_0$ совпадает с наибольшей алгебраической кратностью собственного значения этой матрицы.

Если матрица $\mathbf A_0$ имеет простое собственное число $z_1$, то система уравнений $\mathbf v'(t)=\mathbf A_0\mathbf v(t)$ имеет решение $e^{z_1t}\mathbf c$, где $\mathbf c$ – собственный вектор, соответствующий $z_1$. Из приведенных выше равенств для элементов $b_{ij}(x)$ матрицы $\mathbf B(x)$ и условия (1.11) теоремы 3 при $m=1$ следует, что все элементы матрицы $\mathbf R(t)$ принадлежат пространству $L^1[0,\infty)$. Применяя далее утверждение задачи 29 из книги [17; гл. III] (см. также [2; теорема 1.8.1]), получаем, что система (5.3) имеет решение $\mathbf v$, представимое при $t\to\infty$ в виде

$$ \begin{equation*} \mathbf v(t)=e^{z_1t}(\mathbf c+o(1)). \end{equation*} \notag $$
Остается заметить, что решения $\mathbf y(x)$ и $\mathbf v(t)$ связаны равенством
$$ \begin{equation*} \mathbf y=\mathbf D\mathbf v(\ln x). \end{equation*} \notag $$
Следовательно, уравнение (3.1) имеет решение $y$, представимое при $x\to\infty$ в виде
$$ \begin{equation*} \begin{alignedat}{2} y^{[s-1]}(x) &=x^{z_1-s+1/2}(c_s+o(1)), &\qquad s &=1,2,\dots,n, \\ y^{[n+s-1]}(x) &=x^{z_1+n+\nu-s +1/2}(c_{n+s}+o(1)), &\qquad s &=1,2,\dots,n, \end{alignedat} \end{equation*} \notag $$
где $c_1,c_2,\dots,c_{2n}$ – координаты собственного вектора $\mathbf c$ матрицы $\mathbf A_0$, причем $c_1\ne 0$.

Пусть теперь $z_1$ – собственное число матрицы $\mathbf A_0$ кратности $l_1$, $2\leqslant l_1\leqslant m$. Этому числу, как мы уже отмечали выше, с точностью до постоянного множителя соответствует только один собственный вектор $\mathbf c$. Тогда система уравнений $\mathbf v'(t)=\mathbf A_0\mathbf v(t)$ имеет решения, представимые в виде

$$ \begin{equation*} e^{z_1t}\mathbf c\qquad \text{и}\qquad e^{z_1t}t^l\mathbf c+O(e^{z_1t}t^{l-1}),\quad l=1,2,\dots,l_1-1. \end{equation*} \notag $$
Из вышеприведенных таблиц для элементов $b_{ij}(x)$ матрицы $\mathbf B(x)$ и условия (1.11) теперь следует, что все элементы матрицы $t^{m-1}\mathbf R(t)$ принадлежат пространству $L^1[0,\infty)$. Применив далее утверждение задачи 35 из книги [17; гл. III] (см. также [2; теорема 1.10.1]), легко показать, что уравнение (1.1) имеет подсистему фундаментальных решений $\{y_j\}_{j=1}^{l_1}$, представимых при $x\to\infty$ в виде (1.12). Этим завершается доказательство теоремы.

6. Индексы дефекта и спектр

В этом пункте, как прежде, мы предполагаем, что $k$ – фиксированное целое число, $0\leqslant k\leqslant n$, а функции $p$ и $\sigma$ удовлетворяют условиям (1.2) и (1.3). Помимо этих условий здесь будем полагать, что эти функции вещественные. Эти условия обеспечивают справедливость формулы Лагранжа для квазидифференциального выражения $\tau(y)$, определенного в (2.5) и (2.6). А именно, для любых функций $f,g\in\mathscr D(\tau)$, где $\mathscr D(\tau)$ – естественная область определения выражения $\tau(y)$ (см. выше п. 2), выполнено равенство

$$ \begin{equation} \int_\alpha^\beta\overline g\tau(f) -\int_\alpha^\beta f\overline{\tau(g)} =[fg](\beta)-[fg](\alpha),\qquad \alpha,\beta\in[1,\infty), \end{equation} \tag{6.1} $$
где
$$ \begin{equation*} [fg](x)=\sum_{j=0}^{2n-1}(-1)^{n+1-j}f^{[j]}(x)\overline{g^{[2n-1-j]}(x)}, \end{equation*} \notag $$
а квазипроизводные $y^{[j]}$, $j=0,1,\dots,2n-1$, определяются формулами (1.4)(1.7). Как мы отмечали выше, при $k=0$ определение квазипроизводных и квазидифференциального выражения в этой статье и в книге [1; глава 5, § 15, п. 2, формула (3)] совпадают и формула (6.1) в этом случае доказана (см. [1; глава 5, § 15, п. 3]). При $k\geqslant 1$ доказательство формулы (6.1) можно провести с помощью тех же рассуждений.

Квазидифференциальное выражение $(-1)^n\tau(y)$ и, следовательно, дифференциальное выражение $l(y)$ в представлении (1.1) порождает в пространстве $L^2[1,\infty)$ оператор $L_{\max}$ с областью определения $\mathscr D(L_{\max})$ по формуле

$$ \begin{equation*} L_{\max}y=(-1)^n\tau(y),\qquad \text{если}\quad y\in\mathscr D(L_{\max}), \end{equation*} \notag $$
где
$$ \begin{equation*} \mathscr D(L_{\max})=\{y\in\mathscr D(\tau)\mid y\in L^2[1,\infty)\qquad \text{и}\qquad \tau(y)\in L^2[1,\infty)\}, \end{equation*} \notag $$
Оператор $L_{\max}$ называется максимальным оператором, порожденным выражением $(-1)^n\tau(y)$ в пространстве $L^2[1,\infty)$. Рассмотрим сужение этого оператора на область
$$ \begin{equation*} \mathscr D(L_{\min}^0) =\{y\in\mathscr D(L_{\max})\mid\operatorname{supp}y\subset(1,\infty)\}. \end{equation*} \notag $$
Здесь запись $\operatorname{supp} y\subset(1,\infty)$ означает, что носители функций $y$ являются компактами в $(1,\infty)$, т.е. функции $y\in\mathscr D(L_{\min}^0)$ обращаются в нуль в окрестностях точек $1$ и $\infty$. Замыкание этого оператора в пространстве $L^2[1,\infty)$ и его область определения обозначим $L_{\min}$ и $\mathscr D(L_{\min})$ соответственно. Оператор $L_{\min}$ называется минимальным оператором, порожденным выражением $(-1)^n\tau(y)$ в пространстве $L^2[1,\infty)$. Из (6.1) легко извлечь, что оператор $L_{\min}$ является симметрическим оператором.

Заметим также, что если $f,g\in\mathscr D(L_{\max})$, то из формулы (6.1) следует, что пределы

$$ \begin{equation*} \lim_{\alpha\to 1}[fg](\alpha)(:=[fg](1))\qquad \text{и}\qquad \lim_{\beta\to\infty}[fg](\beta)(:=[fg](\infty)) \end{equation*} \notag $$
существуют и справедливо равенство
$$ \begin{equation} (L_{\max}f,g)-(f,L_{\max}g)=[fg](\infty)-[fg](1), \end{equation} \tag{6.2} $$
где $(\,\cdot\,,\,\cdot\,)$ – скалярное произведение в пространстве $L^2[1,\infty)$.

Отметим, что максимальный и минимальный операторы взаимно сопряжены, а именно,

$$ \begin{equation*} L^*_{\min}=L_{\max},\qquad L^*_{\max}=L_{\min}. \end{equation*} \notag $$

Определим дефектные подпространства $D^\pm$ оператора $L_{\min}$, полагая

$$ \begin{equation*} D^\pm=\{y\in\mathscr D(L_{\max})\mid L_{\max}y=\pm iy\} =\{y\in\mathscr D(\tau)\cap L^2[1,\infty)\mid(-1)^n\tau(y)=\pm iy\}. \end{equation*} \notag $$
Числа $n_\pm=\dim D^\pm$ называются дефектными числами этого оператора. Согласно общепринятой терминологии, условия (1.2) и (1.3) позволяют говорить, что левый конец луча $[1,\infty)$ является регулярным концом для рассматриваемого дифференциального выражения, а правый конец сингулярным. В этом случае известно [18; разд. 5.2] (см. также [19; добавление II, п. 126, теорема $4^*$]), что $n\leqslant n_\pm\leqslant 2n$ (в нашем случае с использованием квазипроизводных другого вида, порождаемых матрицами $\mathbf{F}_{2n}^k$, доказательство сохраняется). Кроме того, условие вещественности функций $p$ и $\sigma$, влечет равенство $n_+=n_-$. Следовательно,
$$ \begin{equation*} n_+=n_-=n+s,\qquad \text{где}\quad 0\leqslant s\leqslant n. \end{equation*} \notag $$
Число $s$ определяется поведением функции $q$ на бесконечности. Вопрос о том, какие значения являются недопустимыми для этого числа даже в случае регулярной функции $q$ при $2n\geqslant 6$ является открытым. Известно, что при $2n\leqslant 4$ число $s$ может принимать все возможные значения.

Равенство дефектных чисел позволяет утверждать, что в рассматриваемом случае оператор $L_{\min}$ допускает самосопряженные расширения и эти расширения описываются следующей теоремой.

Теорема 4. Пусть функции $p$ и $\sigma$ вещественны, выполнены условия (1.2), (1.3), а индексы соответствующего минимального оператора $L_{\min}$ равны $(n+s,n+s)$. Тогда область определения $\mathscr D(L)$ произвольного самосопряженного расширения оператора $L_{\min}$ есть совокупность всех функций $y\in\mathscr D(L_{\max})$, подчиненных условиям

$$ \begin{equation} [yw_j](\infty)-[yw_j](1)=0,\qquad j=1,2,\dots,n+s, \end{equation} \tag{6.3} $$
где $w_1,w_2,\dots,w_{n+s}$ – некоторые функции из $\mathscr D(L_{\max})$, линейно независимые по модулю $\mathscr D(L_{\min})$, такие, что
$$ \begin{equation} [w_jw_i](\infty)-[w_jw_i](1)=0,\qquad j,i=1,2,\dots,n+s. \end{equation} \tag{6.4} $$
Обратно, каковы бы ни были функции $w_1,w_2,\dots,w_{n+s}\in\mathscr D(L_{\max})$, линейно независимые по модулю $\mathscr D(L_{\min})$ и удовлетворяющие соотношениям (6.3), совокупность всех функций $y\in\mathscr D(L_{\max})$, подчиненных условиям (6.3), есть область определения некоторого самосопряженного расширения $L$ оператора $L_{\min}$.

Если помимо указанных условий выполнены условия (1.8), то индексы дефекта оператора $L_{\min}$ равны $(n,n)$, а условия (3.2) принимают упрощенную форму

$$ \begin{equation*} U_j(y):=\sum_{s=1}^{2n}a_{j,s}y^{[s-1]}(1)=0,\qquad j=1,2,\dots,n, \end{equation*} \notag $$
где числа $a_{j,s}$ таковы, что
$$ \begin{equation*} \sum_{l=1}^n(a_{j,l}\overline a_{s,2n-l +1} -a_{s,2n-l+1}\overline a_{s,l})=0,\qquad j,s=1,2,\dots,n. \end{equation*} \notag $$
При этом непрерывный спектр произвольного самосопряженного расширения $L$ совпадает с положительной полуосью $\mathbb R^+=[0,\infty)$, а на отрицательной полуоси могут находиться только конечное число собственных значений конечной кратности.

В случае $k=0$ доказательство первой части теоремы 4 имеется в книге [1; § 18, пп. 1–3], а также в несколько ином виде в [19; добавление 2, пп. 125, 127]. В случае $k\geqslant 1$ доказательство принципиально не меняется, аналогичные рассуждения проходят, если квазипроизводные вводить согласно нашим определениям. Доказательство второй части теоремы 4 при выполнении дополнительных условий (1.8) имеется в работе авторов [9].

В случае вещественных функций $p$ и $\sigma$, починенных условиям теоремы 3, определение максимального и минимального операторов сохраняется. В этом случае мы имеем следующий результат.

Теорема 5. Пусть выполнены условия теоремы 3, причем функции $p$ и $\sigma$ вещественны. Обозначим через $\kappa=\kappa(\nu)$ число корней полинома $\mathscr F_{2n}(z,\nu,k)$, лежащих в открытой левой полуплоскости $\operatorname{Re}z<0$. Тогда максимальное число линейно-независимых решений уравнения (1.1), принадлежащих пространству $L^2[1,+\infty)$, при $\nu>0$ равно $\kappa\geqslant n$ (в частности, индексы дефекта оператора $L_{\min}$ равны $(\kappa,\kappa)$. При этом спектр произвольного самосопряженного расширения $L_{\min}$ дискретен.

При $\nu=0$ число $\kappa(0)$ невещественных значений $\lambda$ равно $n$.

Доказательство. Пусть $\nu>0$. Из формул (1.12) следует, что фундаментальная система решений уравнения (1.1) состоит из функций вида

$$ \begin{equation*} y (x)=x^{z-1/2}(\ln x)^s(c+o(1)), \end{equation*} \notag $$
где $s\geqslant 0$ – целое число, $c\ne 0$, а $z$ – корень уравнения $\mathscr F_{2n}(z,\nu,k)=0$. Функции такого вида принадлежат пространству $L^2[1,\infty)$ тогда и только тогда, когда $\operatorname{Re}z<0$. Следовательно, индексы дефекта оператора $L_{\min}$ равны $(\kappa,\kappa)$. В случае $k=0$ в книге [18; разд. 5.2] доказано, что $\kappa\geqslant n$, а в книге [5] доказана дискретность спектра произвольного самосопряженного расширения минимального оператора. Более того, показано, что резольвента таких расширений является мероморфной функцией от $\lambda$ с ядром Гильберта–Шмидта, если $\lambda$ не совпадает с собственным значением. При $k\geqslant 1$ аргументы доказательств сохраняются.

Пусть $\nu=0$. Так же, как в случае $\nu>0$, получаем, что индексы дефекта $n_\pm$ минимального оператора $L_{\min}$ совпадают с числом корней в открытой левой полуплоскости $\operatorname{Re}z<0$ полинома от $z$

$$ \begin{equation*} \mathscr F_{2n}(z,0,k,\lambda) =\alpha\prod_{j=0}^{n-1}\biggl[z^2-\biggl(\frac{1}{2}+j\biggr)^2\biggr] +(-1)^nk!\beta-(-1)^n\lambda \end{equation*} \notag $$
при $\lambda=\pm i$, т.е. числом к корней в полуплоскости $\operatorname{Re}z<0$ уравнений
$$ \begin{equation} \mathscr F_{2n}(z,0,k,\pm i)=0. \end{equation} \tag{6.5} $$
Воспользуемся равенством
$$ \begin{equation*} \mathscr F_{2n}(-\overline z,0,k,\lambda) =\mathscr F_{2n}(z,0,k,\overline\lambda\,). \end{equation*} \notag $$
Число корней уравнений (6.5) в открытой правой полуплоскости $\leqslant 2n- n_\pm$, поэтому из последних равенств получаем $n_++n_-\leqslant 2n$. С другой стороны известно [18], что для дефектных чисел $n_+$ и $n_-$ выполняются неравенства $n\leqslant n_+$, $n_-\leqslant 2n$. Следовательно, что $n_+=n_-=n$. Этим заканчивается доказательство теоремы 5.

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. М. А. Наймарк, Линейные дифференциальные операторы, Наука, М., 1969  mathscinet  zmath
2. M. S. P. Eastham, The Asymptotic Solution of Linear Differential Systems, Clarendon Press, Oxford, 1989  mathscinet
3. Ф. Хартман, Обыкновенные дифференциальные уравнения, Мир, М., 1970  mathscinet  zmath
4. C. А. Орлов, “Об индексе дефекта линейных дифференциальных операторов”, Докл. АН СССР, 92:3 (1953), 483–486  mathscinet
5. R. B. Paris, A. D. Wood, “On the $\mathscr L_2(I)$ nature of solutions of $n$-th order symmetric differential operator and McLeod's conjecture”, Proc. Roy. Soc. Edinburg Sect. A, 90 (1981), 209–236  crossref  mathscinet
6. R. M. Kauffman, “On the limit-$n$ classification of ordinary differential operators with positive coefficients”, Proc. London Math. Soc. (3), 3:35 (1977), 496–526  crossref  mathscinet
7. R. B. Paris A. D. Wood, Asymptotics of High Order Differential Equations, Pitman Res. Notes in Math. Ser., 129, Longman Sci. & Tech., Harlow, 1986  mathscinet
8. К. А. Мирзоев, “О теореме Орлова об индексе дефекта дифференциальных операторов”, Докл. АН, 380:5 (2001), 591–595  mathnet  mathscinet
9. Н. Н. Конечная, К. А. Мирзоев, А. А. Шкаликов, “Об асимптотике решений двучленных дифференциальных уравнений с сингулярными коэффициентами”, Матем. заметки, 104:2 (2018), 231–242  mathnet  crossref  mathscinet
10. Н. Н. Конечная, К. А. Мирзоев, “Главный член асимптотики решений линейных дифференциальных уравнений с коэффициентами-распределениями первого порядка”, Матем. заметки, 106:1 (2019), 74–83  mathnet  crossref  mathscinet
11. Н. Н. Конечная, К. А. Мирзоев, “Об асимптотике решений линейных дифференциальных уравнений нечетного порядка”, Вестн. Моск. ун-та. Сер. 1. Матем., мех., 2020, № 1, 23–28  mathnet  mathscinet  zmath
12. К. А. Мирзоев, А. А. Шкаликов, “Дифференциальные операторы четного порядка с коэффициентами-распределениями”, Матем. заметки, 99:5 (2016), 788–793  mathnet  crossref  mathscinet
13. А. М. Савчук, А. А. Шкаликов, “Операторы Штурма–Лиувилля с сингулярными потенциалами”, Матем. заметки, 66:6 (1999), 897–912  mathnet  crossref  mathscinet  zmath
14. А. М. Савчук, А. А. Шкаликов, “Операторы Штурма–Лиувилля с потенциалами-распределениями”, Тр. ММО, 64, М., 2003, 159–212
15. J. Eckhardt, F. Gesztesy, R. Nichols, G. Teschl, “Weyl–Titchmarsh theory for Sturm–Liouville operators with distributional potentials”, Opuscula Math., 33 (2013), 467–563  crossref  mathscinet
16. К. А. Мирзоев, “Операторы Штурма–Лиувилля”, Тр. ММО, 75, МЦНМО, М., 2014, 335–359  mathnet
17. Э. А. Коддингтон, Н. Левинсон, Теория обыкновенных дифференциальных уравнений, ИЛ, М., 1958  mathscinet  zmath
18. W. N. Everitt, L. Marcus, Boundary Value Problems and Symplectic Algebra for Ordinary Differential and Quasi-Differential Operators, Amer. Math. Soc., Providence, RI, 1999  mathscinet
19. Н. И. Ахиезер, И. М. Глазман, Теория линейных операторов в гильбертовом пространстве, Т. 2, Вища школа, Харьков, 1978  mathscinet

Образец цитирования: Н. Н. Конечная, К. А. Мирзоев, А. А. Шкаликов, “Об асимптотике решений двучленных дифференциальных уравнений”, Матем. заметки, 113:2 (2023), 217–235; Math. Notes, 113:2 (2023), 228–242
Цитирование в формате AMSBIB
\RBibitem{KonMirShk23}
\by Н.~Н.~Конечная, К.~А.~Мирзоев, А.~А.~Шкаликов
\paper Об асимптотике решений двучленных дифференциальных уравнений
\jour Матем. заметки
\yr 2023
\vol 113
\issue 2
\pages 217--235
\mathnet{http://mi.mathnet.ru/mzm13882}
\crossref{https://doi.org/10.4213/mzm13882}
\transl
\jour Math. Notes
\yr 2023
\vol 113
\issue 2
\pages 228--242
\crossref{https://doi.org/10.1134/S0001434623010261}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85149996584}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/mzm13882
  • https://doi.org/10.4213/mzm13882
  • https://www.mathnet.ru/rus/mzm/v113/i2/p217
  • Эта публикация цитируется в следующих 4 статьяx:
    Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Математические заметки Mathematical Notes
    Статистика просмотров:
    Страница аннотации:292
    PDF полного текста:31
    HTML русской версии:195
    Список литературы:45
    Первая страница:20
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2024