УДК 517.968
Б01: 10.21779/2542-0321-2019-34-1-56-60 Т.И. Ибавов
Задача Коши для системы п дифференциальных уравнений с производными
дробного порядка Капуто
Дагестанский государственный университет; Россия, 367001, г. Махачкала, ул. М. Гаджиева, 43а; [email protected]
В работе исследована задача Коши для системы п линейных дифференциальных уравнений с дробной производной Капуто. При достаточном условии существования дробной производной производная Капуто представлена через дробную производную Римана-Лиувилля. С помощью преобразований Лапласа и метода Крамера получены выражения для образа Лапласа, разложение выражений для образа Лапласа по корням характеристического уравнения. Посредством известных соотношений для образа Лапласа выведены выражения для оригиналов функции. Результаты сформулированы в виде теоремы. Показано, что при п = 2 полученные решения переходят в известные ранее решения для системы двух дифференциальных уравнений с дробными производными.
Ключевые слова: дробная производная Капуто, преобразование Лапласа, задача
Коши.
Введение
Проникновение идей фрактальной геометрии [3] в естествознание сформировало новую концепцию - концепцию фракталов. Особый интерес в концепции фрактала представляет аналитический подход, основанный на использовании математического аппарата дробных дифференциальных уравнений. Многие динамические процессы во фрактальных структурах, как указано в работах [2; 3; 5-7; 9; 10], описываются дробными дифференциальными уравнениями.
Исследованиям задачи Коши для дифференциальных уравнений дробного порядка посвящено много работ. В монографии [1] приведено решение задача Коши для дифференциального уравнения с дробной производной Римана-Лиувилля, доказаны существование и единственность решения. В работах [2; 3] исследована задача Коши для осцилляторного уравнения с дробной производной Капуто. В [4] проведено качественное исследование динамической системы, описываемой двумя дифференциальными уравнениями с производными дробного порядка. В [5] исследована задача Коши для уравнений Эйри с дробной производной Капуто. Работа [6] посвящена исследованию задачи Коши для уравнения нечётного порядка с дробной производной. В [7; 8] исследована задача Коши для системы трех дифференциальных уравнений с дробной производной Капуто. Работы [9; 10] посвящены исследованию нелинейных динамических систем, описываемых дробными дифференциальными
уравнениями. В [9] на основе дробных дифференциальных уравнений разработана и исследована динамика электронов в ветвящихся системах газоразрядных каналов со спадающей концентрацией газа.
Нами исследуется задача Коши для системы п дифференциальных уравнений с дробными производными Капуто.
Задача Коши для системы п дифференциальных уравнений с дробной
производной Капуто
Рассмотрим в области в систему п дифференциальных уравнений с производ ными дробного порядка вида
50, х = лх , (1)
где 0 <а < 1, а1 при г = 1, 2, ... п, 1 = 1, 2, ... п - некоторые числа,
1 , и '(5)
5 0{и (,) =
Л =
г (1 - а )0 (, - 5) ( а 11 а 12 а 21 а 22
у а п1 а п 2
а
- дробная производная Капуто,
Л
а1п а 2 п
а пп у
X (,) =
( х1(,) ^
X 2^ ) Vхп (,).
Предположим, что Л Ф 0. Для системы (1) рассмотрим следующую задачу. Задача. Найти решение х (г) - системы уравнений (1), удовлетворяющее
начальным условиям х (0) = х 0).
Пусть в а с ЛС 2[0,Т] - область определения оператора 5 . Тогда для любого вектора X(г) е Эа справедливо равенство [2]
50 ,х1 (г) = В 0 ,х1 (г )■
х1 (0)
Г (1 - а ) г
г = 1, 2, ..., п.
Воспользовавшись равенствами (2), получим
В 0 ,х,(г) = £ а^х 1 +
х, (0)
г (1 - а ) ,
г = 1, ..., п.
(2)
(3)
Предполагая, что вектор х(г) еВа - решение системы (3), применим к системе
(3) преобразования Лапласа. Тогда относительно вектора х (р) = (х 1 (р ), X 2( р ), . ., х ( р ) У получим следующую систему п дифференциальных
уравнений:
р - а г у ~,(р)-£ а ц х; (р) = х,(0 у р г = 1, 2, ..., п.
' ]=1
•* 1
Запишем систему (4) в матричной форме:
в ■ х (р ) = х (0) • р
а - 1
(4)
(5)
1=1
где B =
( a
Р - a 11 - a 21
- a 12
а
Р - a 22
- a 1 n
- a 2 n
- a n1 - an 2
а
p - a
nn 1
Определитель матрицы В отличен от нуля, то есть ёй В Ф 0 . Тогда решение системы (4) будет иметь вид
, п п-1 , , ч п _,,, , п
хк(0)ра(п-1) + Е Ра(п-1) Е МП= + Е л1 3 I = 2 /, ]=1 3=1
i,3 ф к (6)
~к (Р) =
,1-а
an
-TrA -pa(n-1) + V (-1)lpa(n-l)Mn= , + y A1 j l=2 j= 1
к = 1, 2, ..., п,
где Мп - 1 - миноры порядков (п - I) диагональных элементов матрицы А, в ко-
i, 3, i = 3
торой вместо ]-го столбца стоит столбец свободных членов системы (5); ТгЛ - след матрицы А; л1 3 - алгебраические дополнения первой строки матрицы А.
Систему (6) можно представить в виде
x
( Р ) =
П (.-à)p
ХШт У (- 1Г 1
П ti,-Äj)
,*=1
- а
1 - Àmp
У M.'j ',= j
n - 1 ,, j = 1 n
+ y ui-y (- 1 )m+1
L-! 1 - a ( n -1 ) L-! V V
l = 2 p m = 1
П (д , - a j)
1 - Amp
, к = 1, 2, ..., n,
где i = 1, 2,..., п, Х1, Хг, ■■■, Хп - корни характеристического уравнения системы (5). В случае, когда Яер > 1, имеют место следующие соотношения [1]:
Ю Я 1 1
| е - гг р- 1 £а п (/а 2= —, Ы < 1, А а,п 1 - г 11
œ
= J e ptEa 11 At" 1 dt , Re p > 1
p (1 -X-p a) 0
a
(7)
p1 -a(1 -X-p-a) 0
= J e ptt aEa1 -a(Xta) dt , Re p > 1.
Воспользовавшись соотношениями (7) для системы (6), получим выражение для оригиналов функций Х1 ((), X2(íX ..., X ):
x
(t ) =
П (X,-X
i, j = 1 i * j
[ t
n - 1
x,
(0 )У П (X, -X j ) E a ,1 -a ( n - 1)(Amt a )
m = 1 i, j = 1 i * m
J = 1
j = 1
1
1
+
к
i = 1,2,..., n, (8)
+ 2 Г"(4+° 2 МП -+ 12 П -А; )Е, Л- а (п-оА а)
к —1 , ; — 1 ^ ? з т =! ,, j=1
' ^ , Ф т
,,; * к
да _к
где е „(2) — 2 -, « > 0, ^ > 0 - функция Миттага-Леффлера.
а 'р к = 0 Г(ак +
Следовательно, имеет место следующая теорема.
Теорема. Вектор X(1) е Ва будет решением системы (1), если ёй А Ф 0 и элементы XI (1), , = 1, 2,. ., п этого вектора представимы в виде (8).
Рассмотрим случаи, когда п = 2. При п = 2 система (8) примет вид
Х1(0 = () + ЕдА ))+ ) - Х(0)а22 +Г Х1(0) ЕЖ ) - ЕадАа )) , 2 2/ У + (а21^2 -
х2 (t) = ^ ((,1 ОМ* ) + Eа,1 (Я21а ))+ Q)a21 - х2(Q)a11 + Г- x2
2VГ2 + (a21 a12 " a11a22 )
+ x,(Q)a 21 - x 2(0)ап + r-x 2(q) ( ,1(^2 ,g ^jCV )), (9)
где л12 —у±л1у2 + (а 21 а12 - а11 а 22 ) - корни характеристического уравнения,
У — а11 + а 22 2 '
Решение (9) совпадает с решением, полученным в [4] для системы двух линейных дифференциальных уравнений с дробными производными.
Литература
1. Самко С.Г., Килбас А.А., Маричев О.И. Интегралы и производные дробного порядка и некоторые их приложения. - Минск: Наука и техника, 1987. - C. 498.
2. Нахушев А.М. Дробное исчисление и его применение. - М.: Физматлит, 2QQ3. - C. 272.
3. Бейбалаев В.Д. Решение начальной задачи для дифференциального уравнения «фрактального» осциллятора // Вестник Самарского гос. техн. ун-та. Сер.: Физ.-мат. науки. - 2QQ9. - Т. 1 (19).
4. Назаралиев М.А., Бейбалаев В.Д. Динамические системы, описываемые двумя дифференциальными уравнениями с производными дробного порядка // Владикавказский математический журнал. - 2Q13. - Т. 15, № 1. - С. 3Q-4Q.
5. Паровик Р.И. Задача Коши для обобщенного уравнения Эйри // Доклады АМАН. - 2Q14. - Т. 16, № 3. - С. 64-69.
6. Псху А.В. Задача Коши для уравнения нечётного порядка с дробной производной // Материалы Международной конференции «Современные методы теории краевых задач». МГУ им. М.В. Ломоносова. - М.: ООО «МАКС Пресс», - 2Q18. -C. 185.
7. Ибавов Т.И. Решение начальной задачи для одной системы трёх дифференциальных уравнений с производной дробного порядка // Вестник ДГУ. - 2Q18. -
Т. 33, вып. 1. - С. 78-84.
8. Ибавов Т.И. Задача Коши для одной системы трёх дифференциальных уравнений с производной дробного порядка // Сборник материалов IV международной научной конференции «Актуальные проблемы прикладной математики». -Нальчик, 2018. - С. 104.
9. Тренькин А.А., Карелин В.И., Алисултанов З.З., Бейбалаев В.Д., Рагимха-нов Г.Б. Динамика электронов в ветвящихся системах газоразрядных каналов со спадающей концентрацией газа // Нелинейный мир. - 2017.- Т. 15, № 3. - С. 24-31.
10. Аливердиев А.А., Бейбалаев В.Д., Мейланов Р.Р., Назаралиев М.А. Обобщённый нелинейный осциллятор Дуффинга // Тезисы докладов XIII международной научной конференции «Теория операторов, комплексный анализ и математическое моделирование». - Владикавказ, 2016. - C. 141-142.
Поступила в редакцию 26 февраля 2019 г.
UDC 517.968
DOI: 10.21779/2542-0321-2019-34-1-56-60
Cauchy problem for a system of n differential equations with the Caputo
fractional derivative
T.I. Ibavov
Dagestan State University; Russia, 367001, Makhachkala, M. Gadzhiev st., 43a; [email protected]
The paper investigates the Cauchy problem for system of n linear differential equations with the Caputo fractional derivative. At sufficient condition for the existence of fractional derivative, the Caputo derivative is represented by the fractional derivative of Rieman-Liouville. An expression for Laplace's image is obtained by means of Laplace transform and Kramer rules. Decomposition of expressions for the image of Laplace is obtained by the roots of characteristic equation. Using the known relations for the image of Laplace, expressions for originals of functions is obtained. The results are formulated as a theorem. It is shown that for n = 2 the obtained solutions pass into previously known solutions for a system of two differential equations with the Caputo fractional derivative.
Keywords: fractional derivative Caputo, Laplace transform, Cauchy problem.
Received 26 February, 2019