|
Задача Наймарка для обыкновенного дифференциального уравнения дробного порядка
Л. Х. Гадзова Институт прикладной математики и автоматизации – филиал Федерального государственного бюджетного научного учреждения "Федеральный научный центр "Кабардино-Балкарский научный центр Российской академии наук", г. Нальчик
Аннотация:
Рассматривается задача для обыкновенного дифференциального уравнения дробного порядка, где краевые условия задаются в форме линейных функционалов. Это позволяет охватить достаточно широкий класс линейных локальных и нелокальных условий. Дробная производная понимается в смысле Герасимова–Капуто. Получено необходимое и достаточное условие однозначной разрешимости задачи. В терминах специальных функций найдено представление решения исследуемой задачи. Доказана теорема существования и единственности решения.
Библиография: 27 названий.
Ключевые слова:
дробная производная Герасимова–Капуто, задача Наймарка,
производная дробного порядка, уравнение дробного порядка,
функционал, функция Миттаг-Леффлера.
Поступило: 14.02.2023 Исправленный вариант: 07.03.2023
1. Введение Рассмотрим уравнение
$$
\begin{equation}
\partial_{0x}^{\alpha} u(x)+\lambda u(x)=f(x), \qquad x\in(0,1),
\end{equation}
\tag{1.1}
$$
где $n-1<\alpha\leqslant n$, $n\in {\mathbb N}$, $\lambda \in {\mathbb R}$, $\partial_{0x}^{\alpha}{u(x)}$ – производная Герасимова–Капуто [1; c. 11]. Дробное исчисление привлекает внимание исследователей во многих прикладных областях. Многие явления в механике жидкости, вязкоупругости, динамики населения и в других науках можно описать моделями с помощью математического аппарата из теории дробного исчисления. В связи с этим дифференциальные уравнения дробного порядка в настоящее время интенсивно развиваются. Обширный обзор литературы по дробному исчислению и его применению можно найти в монографиях [1]–[12]. В работе [13] Барретт построил в терминах функции Миттаг-Леффлера решение начальной задачи для обыкновенного дифференциального уравнения дробного порядка
$$
\begin{equation*}
\begin{gathered} \, D_{ax}^{\alpha}u(x)=\lambda u(x)+v(x), \qquad a<x<b, \\ \lim_{x\to a}D_{ax}^{\alpha-i}u(x)=b_i, \quad i=1,2,\dots, n,\qquad \text{где} \quad n-1<\alpha \leqslant n, \quad n \in {\mathbb N}, \end{gathered}
\end{equation*}
\notag
$$
где $D_{ax}^{\alpha}$ – оператор дробного интегро-дифференцирования Римана–Лиувилля. С течением времени уравнения такого вида стали называть уравнениями Барретта [1; c. 98], [9; c. 137], [14; c. 192]. Существенный вклад в изучение теории дробных дифференциальных уравнений внесли авторы работ [15], [16] и [17]. Задаче Коши в нелокальной постановке посвящена работа [15]. В работе [16] была исследована краевая задача для дифференциального оператора дробного порядка типа Штурма–Лиувилля, а работа [17] посвящена видоизмененной задаче Коши для дифференциального уравнения дробного порядка. В работе [18] исследована задача Коши в классической постановке для обыкновенного дифференциального уравнения дробного порядка. Начальная задача для линейного обыкновенного дифференциального уравнения дробного порядка решена в работе [19]. Отметим также работы [9] и [14], в которых можно найти изложение результатов и библиографию работ, связанных с задачей Коши для линейных обыкновенных дифференциальных уравнений дробного порядка. В работах [20] и [21] исследованы задачи с обобщенными краевыми условиями для линейных обыкновенных дифференциальных уравнений. В работе [22; c. 16] Наймарк предложил задачу определения функции $y$, удовлетворяющую условиям
$$
\begin{equation*}
L(y)=0,\qquad U_{\nu}(y)=0, \quad \nu=1,2,\dots,n,
\end{equation*}
\notag
$$
где $L(y)$ – дифференциальное выражение $n$-го порядка, $U_{1}(y),U_{2}(y),\dots,U_{n}(y)$ – линейно независимые линейные непрерывные функционалы. Такие задачи Наймарк называет обобщенными краевыми задачами. В [23; c. 247] приводятся различные виды обобщения краевых условий. Обобщенная краевая задача для уравнения второго порядка с дробной производной Римана–Лиувилля и уравнения с оператором дробного дискретно-распределенного дифференцирования была исследована в работах [24] и [25]. В работе [26] методом функции Грина получено представление решения для уравнения (1.1) при $1<\alpha\leqslant 2$. В настоящей работе для уравнения (1.1) сформулирована и решена задача с краевыми условиями в форме линейных функционалов, которые позволяют охватить достаточно широкий класс линейных локальных и нелокальных условий. В терминах специальных функций найдено представление решения. Получено необходимое и достаточное условие разрешимости исследуемой задачи. Доказана теорема существования и единственности решения.
2. Вспомогательные сведения Оператор $D_{sx}^{\alpha}$ дробного интегро-дифференцирования в смысле Римана–Лиувилля порядка $\alpha$ определяется следующим образом [1; c. 9]:
$$
\begin{equation*}
\begin{gathered} \, D_{sx}^{\alpha}u(x)= \frac{\operatorname{sign}(x-s)}{\Gamma(-\alpha)} \int_{s}^{x}\frac{u(t)\,dt}{|x-t|^{\alpha+1}}\,, \qquad \alpha<0, \\ D_{sx}^{\alpha}u(x)=u(x),\qquad \alpha=0, \\ D_{sx}^{\alpha}u(x)=\operatorname{sign}^{n}(x-s)\, \frac{d^n}{dx^n}D_{sx}^{\alpha-n}u(x), \qquad n-1<\alpha \leqslant n, \quad n \in {\mathbb N}, \end{gathered}
\end{equation*}
\notag
$$
где $\Gamma(z)$ – гамма-функция Эйлера. Производная Герасимова–Капуто $\partial_{0x}^{\alpha}{u(x)}$, определяется с помощью равенства [14], [27]
$$
\begin{equation}
\partial_{sx}^{\alpha}u(x)=\operatorname{sign}^{n}(x-s) D_{sx}^{\alpha-n}u^{(n)}(x), \qquad n-1<\alpha \leqslant n, \quad n \in {\mathbb N}.
\end{equation}
\tag{2.1}
$$
Формула дробного дифференцирования свертки Лапласа двух функций
$$
\begin{equation}
D_{0x}^{\alpha}(f\ast g)(x)=(D_{0x}^{\alpha}g\ast f)(x)+ f(x)\lim_{x \to 0}D_{0x}^{\alpha-1}g(x), \qquad 0<\alpha\leqslant1,
\end{equation}
\tag{2.2}
$$
где
$$
\begin{equation*}
(g\ast h)(x)=\int_{0}^{x}g(t)h(x-t)\,dt.
\end{equation*}
\notag
$$
Для произвольного $\alpha \in {\mathbb R}$ имеет место закон композиции операторов дробного интегро-дифференцирования [11; c. 15]
$$
\begin{equation}
D_{0x}^{\alpha}D_{0x}^{\beta}g(t)=D_{0x}^{\alpha+\beta}g(t), \qquad \beta\leqslant 0.
\end{equation}
\tag{2.3}
$$
Также для оператора дробного интегро-дифференцирования справедлива формула дробного интегрирования по частям [11; c. 15]
$$
\begin{equation}
\int_{a}^{b}g(x)D_{ax}^{\alpha}h(x)\,dx= \int_{a}^{b}h(x)D_{bx}^{\alpha}g(x)\,dx, \qquad \alpha\leqslant0.
\end{equation}
\tag{2.4}
$$
Функция (см., например, [9], [11])
$$
\begin{equation*}
E_{\delta,\mu}(x)= \sum_{j=0}^{\infty}\frac{x^j}{\Gamma(\delta j+\mu)}
\end{equation*}
\notag
$$
называется функцией Миттаг-Леффлера. Для функции Миттаг-Леффлера справедлива формула дробного интегрирования и дифференцирования
$$
\begin{equation}
D_{0x}^{\beta}x^{\mu-1}E_{\alpha,\mu}(\lambda x^{\alpha})= x^{\mu-\beta-1}E_{\alpha,\mu -\beta}(\lambda x^{\alpha}),
\end{equation}
\tag{2.5}
$$
где $\beta \in {\mathbb R}$, $\mu > 0$, если $\beta \notin {\mathbb N}\cup \{0\}$, и $\mu \in {\mathbb R}$, если $\beta \in {\mathbb N}\cup \{0\}$. Решение задачи Коши для уравнения (1.1) имеет вид (см., например, [11])
$$
\begin{equation}
u(x)=\int_{0}^{x}f(t)\mathcal{W}_{\alpha}(x-t)\,dt+ \mathcal{W}(x)\overline{u},
\end{equation}
\tag{2.6}
$$
где
$$
\begin{equation}
\begin{gathered} \, \nonumber \mathcal{W}(x)=(\mathcal{W}_{n}(x),\mathcal{W}_{n-1}(x),\dots, \mathcal{W}_{1}(x)),\qquad \overline{u}=\begin{pmatrix} u(0)\\ u'(0)\\ \cdots\\ u^{(n-1)}(0) \end{pmatrix}, \\ \mathcal{W}_{\mu}(x)=\begin{cases} x^{\mu-1}E_{\alpha,\mu}(-\lambda x^{\alpha}), & \text{если}\ x >0, \\ 0, & \text{если}\ x\leqslant 0. \end{cases} \end{gathered}
\end{equation}
\tag{2.7}
$$
Приведем утверждение из [26]. Лемма 1. Пусть
$$
\begin{equation*}
\frac{\partial^{i}}{\partial x^{i}}K(x,t)\in C([0,1]\times[0,1]), \qquad i=0,1,\dots,p,
\end{equation*}
\notag
$$
$\ell$ – линейный ограниченный функционал в пространстве $C^{p}[0,1]$. Тогда справедливо соотношение
$$
\begin{equation}
\ell\biggl[\int_{0}^{1} K(x,t)\,dt\biggr]= \int_{0}^{1}\ell[K(x,t)]\,dt.
\end{equation}
\tag{2.8}
$$
Заметим, что в формуле (2.8) функционал $\ell$ применяется к $K(x,t)$ как к функции переменной $x$. Доказательство. В силу условий, наложенных на $K(x,t)$, можем записать
$$
\begin{equation*}
\int_{0}^{1} K(x,t)\,dt=\lim_{m\to \infty} \frac{1}{m}\sum_{i=0}^{m}K\biggl(x,\frac{i}{m}\biggr).
\end{equation*}
\notag
$$
Причем $K(x,i/m)\in C^{p}[0,1]$ для любого $i$ и $\ell[K(x,t)]$ как функция переменной $t$ является непрерывной на отрезке $[0,1]$, поэтому
$$
\begin{equation*}
\begin{aligned} \, \ell\biggl[\int_{0}^{1} K(x,t)\,dt\biggr]&= \ell\biggl[\lim_{m\to \infty}\,\sum_{i=0}^{m} K\biggl(x,\frac{i}{m}\biggr)\frac{1}{m}\biggr]= \lim_{ m\to \infty}\ell\biggl[\,\sum_{i=0}^{m} K\biggl(x,\frac{i}{m}\biggr)\frac{1}{m}\biggr] \\ &=\lim_{ m\to \infty}\sum_{i=0}^{m}\ell \biggl[K\biggl(x,\frac{i}{m}\biggr)\biggr] \frac{1}{m}=\int_{0}^{1}\ell[ K(x,t)]\,dt. \end{aligned}
\end{equation*}
\notag
$$
Лемма доказана.
3. Постановка задачи и представление решений Регулярным решением уравнения (1.1) назовем функцию $u=u(x)$ из класса $AC^{n-1}[0,1]$ (функций, имеющих абсолютно непрерывные производные до порядка $n-1$ на отрезке $[0,1]$), удовлетворяющую этому уравнению во всех точках интервала $(0,1)$. Задача. Найти регулярное решение уравнения (1.1) в интервале $(0,1)$, удовлетворяющее условиям
$$
\begin{equation}
\ell[u]=b,
\end{equation}
\tag{3.1}
$$
где
$$
\begin{equation*}
b=\begin{pmatrix} b_1\\ b_2\\ \cdots\\ b_n \end{pmatrix}
\end{equation*}
\notag
$$
– заданный вектор-столбец, $b_{k}\in \mathbb{R}$,
$$
\begin{equation*}
\ell=\begin{pmatrix} \ell_1\\ \ell_2\\ \cdots\\ \ell_n \end{pmatrix}
\end{equation*}
\notag
$$
– вектор-столбец, $\ell_{k}$ – заданные линейные ограниченные функционалы с областью определения $D(\ell_{k})\subset AC^{n-1}[0,1]$, $k=1,\dots,n$. Теорема 1. Пусть $n\geqslant 2$,
$$
\begin{equation*}
D(\ell_{k})=C^{n-2}[0,1], \qquad k=1,\dots,n,
\end{equation*}
\notag
$$
функция $f(x)\in C(0,1)$ представима в виде
$$
\begin{equation}
f(x)=D_{0x}^{\alpha-n}g(x), \qquad g(x)\in L(0,1),
\end{equation}
\tag{3.2}
$$
и выполняется условие
$$
\begin{equation}
\det A\ne 0,
\end{equation}
\tag{3.3}
$$
где $A=\|\ell_i[\mathcal{W}_j(x)]\|$ – квадратная матрица, $i,j=1,\dots,n$. Тогда существует единственное регулярное решение задачи (3.1) для уравнения (1.1), определяемое равенством
$$
\begin{equation}
u(x)=\int_{0}^{1} f(t)G(x,t)\,dt+\mathcal{W}(x)A^{-1}b,
\end{equation}
\tag{3.4}
$$
где
$$
\begin{equation*}
G(x,t)=\mathcal{W}_{\alpha}(x-t)- \mathcal{W}(x)A^{-1}\ell[\mathcal{W}_{\alpha}(x-t)].
\end{equation*}
\notag
$$
Доказательство. Применяем условие (3.1) к представлению (2.6) и находим неизвестные константы $\overline{u}$. С учетом леммы 1 получаем систему линейных уравнений следующего вида:
$$
\begin{equation*}
A\overline{u}=b-\mathcal{F},
\end{equation*}
\notag
$$
где
$$
\begin{equation*}
\mathcal{F}=\int_{0}^{1}f(t)\ell[\mathcal{W}_{\alpha}(x-t)]\,dt.
\end{equation*}
\notag
$$
Если $\det A\ne 0$, то существует обратная матрица $A^{-1}$, и можем написать
$$
\begin{equation}
\overline{u}=A^{-1}b-A^{-1}\mathcal{F}.
\end{equation}
\tag{3.5}
$$
Подставляя найденное значение для $\overline{u}$ в уравнение (2.6), после несложных преобразований получаем (3.4). Покажем, что построенное решение является регулярным решением уравнения (1.1). В силу формулы дифференцирования функции Миттаг-Леффлера (2.5) имеем
$$
\begin{equation*}
(\partial_{0x}^{\alpha}+\lambda)\mathcal{W}(x)=0.
\end{equation*}
\notag
$$
Поэтому с учетом (2.6) и (2.7) остается показать, что
$$
\begin{equation*}
\partial_{0x}^{\alpha}(f\ast\mathcal{W}_{\alpha}(x))+ \lambda(f\ast\mathcal{W}_{\alpha}(x))=f(x).
\end{equation*}
\notag
$$
Учитывая представление (3.2) функции $f(x)$ и определение дробной производной Герасимова–Капуто (2.1), получим
$$
\begin{equation*}
\partial_{0x}^{\alpha}(f\ast\mathcal{W}_{\alpha}(x))= D_{0x}^{\alpha-n}\,\frac{d^n}{dx^n} (D_{0x}^{\alpha-n}g\ast\mathcal{W}_{\alpha}(x)).
\end{equation*}
\notag
$$
Используя свойства свертки (2.2), закон композиции операторов дробного интегро-дифференцирования (2.3) и формулу дробного интегрирования по частям (2.4), будем иметь
$$
\begin{equation*}
\begin{aligned} \, D_{0x}^{\alpha-n}\,\frac{d^n}{dx^n}(D_{0x}^{\alpha-n} g\ast\mathcal{W}_{\alpha}(x))&=D_{0x}^{\alpha-n}\,\frac{d^n}{dx^n} (g\ast\mathcal{W}_{n}(x)) =D_{0x}^{\alpha-n}\,\frac{d}{dx} (g\ast\mathcal{W}_{1}(x)) \\ &=D_{0x}^{\alpha-n}(g\ast\mathcal{W}_{1}(0)+ g\ast\mathcal{W}_{0}(0)) \\ &=D_{0x}^{\alpha-n}g(x)- D_{0x}^{\alpha-n}g\ast \lambda \mathcal{W}_{\alpha}(x) =f(x)-\lambda(f\ast \mathcal{W}_{\alpha}(x)). \end{aligned}
\end{equation*}
\notag
$$
С учетом формулы (2.8) покажем, что равенство (3.4) удовлетворяет краевым условиям (3.1):
$$
\begin{equation*}
\ell[f\ast \mathcal{W}_{\alpha}(x)]- AA^{-1}\ell[f\ast \mathcal{W}_{\alpha}(x)]+AA^{-1}b=b.
\end{equation*}
\notag
$$
Потому как $A=\ell[\mathcal{W}(x)]$, $AA^{-1}=E$ – единичная матрица.
4. Случай нарушения условия разрешимости (3.3) Покажем, что если условие разрешимости (3.3) нарушается, то решение задачи, вообще говоря, неединственно. Более точно, однородная задача в этом случае имеет нетривиальное решение. Рассмотрим функцию
$$
\begin{equation*}
\widetilde{u}(x)=\mathcal{W}(x)C,
\end{equation*}
\notag
$$
где $\mathcal{W}(x)$ определяется равенством (2.7),
$$
\begin{equation*}
C=\begin{pmatrix} C_1\\ C_2\\ \cdots\\ C_n \end{pmatrix}
\end{equation*}
\notag
$$
– произвольные постоянные. Из (2.6) следует, что функция $\widetilde{u}(x)$ является решением уравнения
$$
\begin{equation}
\partial_{0x}^{\alpha}\widetilde{u}(x)+\lambda\widetilde{u}(x)=0.
\end{equation}
\tag{4.1}
$$
Так как $\det A=0$, найдется такой ненулевой вектор $\widetilde{C}$, что
$$
\begin{equation*}
\ell[\widetilde{u}(x)]=\ell[\mathcal{W}(x)]\widetilde{C}= A\widetilde{C}=0.
\end{equation*}
\notag
$$
Таким образом, получено, что функция
$$
\begin{equation*}
\widetilde{u}(x)=\mathcal{W}(x)\widetilde{C}
\end{equation*}
\notag
$$
является решением однородного уравнения (4.1) и удовлетворяет нулевым условиям
$$
\begin{equation*}
\ell[\widetilde{u}(x)]=0.
\end{equation*}
\notag
$$
|
|
|
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
|
|
|
1. |
А. М. Нахушев, Дробное исчисление и его применение, Физматлит, М., 2003 |
2. |
А. А. Килбас, Теория и приложения дифференциальных уравнений дробного порядка (курс лекций), Методологическая школа-конференция “Математическая физика и нанотехнологии”, Самара, 2009 |
3. |
М. М. Джрбашян, Интегральные преобразования и представления функций в комплексной области, Наука, М., 1966 |
4. |
K. B. Oldham, J. Spanier, The Fractional Calculus, Academic Press, New York–London, 1974 |
5. |
R. L. Bagley, P. J. Torvik, “Fractional calculus in the transient analysis of viscoelastically damped structures”, AIAA Journal, 23:6 (1985), 918–925 |
6. |
С. Г. Самко, А. А. Килбас, О. И. Маричев, Интегралы и производные дробного порядка и некоторые их приложения, Наука и техника, Минск, 1987 |
7. |
K. S. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley, New York, 1993 |
8. |
I. Podlubny, “Fractional Differential Equations”, Academic Press,, Math. in Sci. and Eng., 198, San Diego, CA, 1999 |
9. |
A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Math. Stud., 204, Elsevier, Amsterdam, 2006 |
10. |
R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, River Edge, NJ, 2000 |
11. |
А. В. Псху, Уравнения в частных производных дробного порядка, Наука, М., 2005 |
12. |
В. В. Учайкин, Метод дробных производных, Изд.“Артишок”, Ульяновск, 2008 |
13. |
J. H. Barrett, “Differential equations of non-integer order”, Canad. J. Math., 6:4 (1954), 529–541 |
14. |
R. Gorenflo, A. A. Kilbas, F. Mainardi, S. V. Rogosin, Mittag-Leffler Functions, Related Topics and Applications, Springer Monogr. in Math., Springer, Heidelberg, 2014 |
15. |
М. М. Джрбашян, А. Б. Нерсесян, “Дробные производные и задача Коши для дифференциальных уравнений дробного порядка”, Изв. АН Армянской ССР. Матем., 3:1 (1968), 3–29 |
16. |
М. М. Джрбашян, “Краевая задача для дифференциального оператора дробного порядка типа Штурма–Лиувилля”, Изв. АН Армянской ССР, 5:2 (1970), 71–96 |
17. |
I. Ozturk, “On the theory of fractional differential equation”, Reports of Adyghe (Circassian) International Academy of Sciences, 3:2 (1998), 35–39 |
18. |
N. Hayek, J. Trujillo, M. Rivero, B. Bonilla, J. C. Moreno, “An extension of Picard–Lindelöff theorem to fractional differential equations”, Appl. Anal., 70:3–4 (1999), 347–361 |
19. |
А. В. Псху, “Начальная задача для линейного обыкновенного дифференциального уравнения дробного порядка”, Матем. сб., 202:4 (2011), 111–122 |
20. |
J. Tamarkin, “Some general problems of the theory of ordinary linear differential equations and expansion of an arbitrary function in series of fundamental functions”, Math. Z., 27:1 (1928), 1–54 |
21. |
А. С. Смогоржевский, “Функции Грина линейных диференциальных систем в области одного измерения”, Матем. сб., 7 (49):1 (1940), 179–196 |
22. |
М. А. Наймарк, Линейные дифференциальные операторы, Наука, М., 1969 |
23. |
Э. Камке, Справочник по обыкновенным дифференциальным уравнениям, Наука, М., 1971 |
24. |
Л. Х. Гадзова, “Обобщенная краевая задача для дифференциального уравнения второго порядка с дробной производной”, Доклады Адыгской (Черкесской) Международной АН, 21:4 (2021), 10–14 |
25. |
L. Kh. Gadzova, “Generalized boundary value problem for a linear ordinary differential equation with a fractional discretely distributed differentiation operator”, Bull. Karaganda Univ. Math. Ser., 106:2 (2022), 108–116 |
26. |
Л. Х. Гадзова, “Обобщённая краевая задача для обыкновенного дифференциального уравнения дробного порядка”, Челяб. физ.-матем. журн., 7:1 (2022), 20–29 |
27. |
О. Г. Новоженова, Биография и научные труды Алексея Никифоровича Герасимова. О линейных операторах, упруго-вязкости, элевтерозе и дробных производных, Перо, М., 2018 |
Образец цитирования:
Л. Х. Гадзова, “Задача Наймарка для обыкновенного дифференциального уравнения дробного порядка”, Матем. заметки, 114:2 (2023), 195–202; Math. Notes, 114:2 (2023), 159–164
Образцы ссылок на эту страницу:
https://www.mathnet.ru/rus/mzm14008https://doi.org/10.4213/mzm14008 https://www.mathnet.ru/rus/mzm/v114/i2/p195
|
Статистика просмотров: |
Страница аннотации: | 209 | PDF полного текста: | 61 | HTML русской версии: | 148 | Список литературы: | 37 | Первая страница: | 16 |
|