Математические заметки СВФУ Июль—сентябрь, 2021. Том 28, № 3
УДК 517.9
НАЧАЛЬНЫЕ ЗАДАЧИ ДЛЯ НЕКОТОРЫХ КЛАССОВ ЛИНЕЙНЫХ ЭВОЛЮЦИОННЫХ УРАВНЕНИЙ С НЕСКОЛЬКИМИ ДРОБНЫМИ ПРОИЗВОДНЫМИ В. Е. Федоров, К. В. Бойко, Т. Д. Фуонг
Аннотация. Исследованы вопросы однозначной разрешимости начальных задач для линейных неоднородных уравнений общего вида с несколькими дробными производными Герасимова — Капуто в банаховых пространствах. Рассмотрена задача Коши для разрешенного относительно старшей дробной производной уравнения, содержащего ограниченные операторы при младших производных, решение представлено с помощью интегралов типа Данфорда — Тейлора. Полученный результат позволил исследовать начальную задачу для линейного неоднородного уравнения с вырожденным оператором при старшей дробной производной при условии, что относительно этого оператора 0-ограниченным является оператор при второй по величине порядка производной. Абстрактные результаты использованы при изучении одного класса начально-краевых задач для уравнений с несколькими производными Герасимова — Капуто по времени и с многочленами от самосопряженного эллиптического дифференциального по пространственным переменным оператора.
Б01: 10.255877SVFU.2021.75.46.006
Ключевые слова: дифференциальное уравнение дробного порядка, дробная производная Герасимова — Капуто, вырожденное эволюционное уравнение, задача Ко-ши, начально-краевая задача.
1. Введение
За последние несколько десятилетий произошел резкий рост интереса исследователей к дробным дифференциальным уравнениям, в первую очередь из-за их возрастающего значения в моделировании различных явлений, возникающих в физике, химии, математической биологии, технике [1,2]. Подробнее о дробных дифференциальных уравнениях и тесно связанных с ними интегро-дифференциальных уравнениях Вольтерра см. в монографиях [3-10].
Задачи для различных классов уравнений с несколькими дробными производными исследовались многими авторами, в частности изучались начально-краевые задачи для телеграфных [11], диффузионных уравнений [12,13] такого
Работа поддержана Приказом 211 Правительства Российской федерации, договор 02.A03.21.0011, и Российским фондом фундаментальных исследований, грант 21—51—54003.
© 2021 Федоров В. Е., Бойко К. В., Фуонг X. Д.
вида, уравнения с запаздыванием [14,15], различные уравнения в локально выпуклых (или просто в банаховых) пространствах с приложениями к уравнениям в частных производных [16-20].
В этой статье мы исследуем вопросы однозначной разрешимости для линейных неоднородных уравнений общего вида с несколькими дробными производными Герасимова — Капуто в банаховых пространствах. Рассмотрена задача Коши для разрешенного относительно старшей дробной производной уравнения, содержащего ограниченные операторы при младших производных, решение представлено с помощью интегралов типа Данфорда — Тейлора. Полученный результат позволил исследовать начальную задачу для линейного неоднородного уравнения с вырожденным оператором при старшей дробной производной при условии, что относительно этого оператора 0-ограниченным является оператор при второй по величине порядка производной. Абстрактные результаты использованы при изучении одного класса начально-краевых задач для уравнений с несколькими производными Герасимова — Капуто по времени и с многочленами от самосопряженного эллиптического дифференциального по пространственным переменным оператора.
Ранее условие относительной ограниченности использовалось, как правило, для пары операторов при старшей производной и при самой функции [13-21]. Отметим также касающиеся вырожденных дробных эволюционных уравнений работы [24-27], близкие к данной работе по объектам и методам исследования. В данной работе вторым оператором в относительно ограниченной паре является оператор при второй по величине порядка дробной производной. Мы продемонстрировали успешность предложенного подхода.
2. Однородное уравнение
Пусть ^ — банахово пространство, ^(^) — пространство линейных ограниченных на ^ операторов. Рассмотрим линейное однородное уравнение дробного порядка
п
= £ ва Лкг(1) (1)
к=1
с заданными начальными условиями
2(г)(0)= г1, I = 0,1,..., то - 1. (2)
Здесь 0 < а1 < а2 < ■ ■ ■ < ап < а, т = [а] — наименьшее целое число, не превосходимое числом а, тк = [ак], к = 1, 2,...,п. Решением задачи (1), (2) будем называть функцию 2 <Е Ст_1(М+; для которой В^кАкг <Е к = 1,2, и выполняются равенства (1) при всех I € М+ :=
М+ и {0} и (2).
Преобразование Лапласа функции к : М+ ^ ^ будем обозначать через к или Ьар[к], если выражение для к слишком громоздкое.
Сначала при фиксированном I € {0,1,..., т — 1} рассмотрим задачу
•г(0(0) = г1, г(к)(0) = 0,к = {0,1,...,т — 1}\{1} , (3)
для уравнения (1). В предположении, что к решению г такой задачи применимо преобразование Лапласа, рассмотрим два случая:
(1) I < тк - 1, тогда (А) = Аак £(А) - Аак-г-1гг;
(2) I > тк - 1, тогда Щ^г(А) = Аак г (А).
Обозначим щ = тш{к € {1, 2,..., п} : I < тк — 1} и из равенств (1), (3) получим следующее:
Aaz(A) - А'
a-,-1z, = Aak Akz(A) - ^ Aak-l-1Akz,,
k= 1
k=ni
2(A) = AOT - ^ Aak A J Aa-,-1Z - ^ Aak-1-1A J z,.
V k=1 ) V k=nl )
Применяя обратное преобразование Лапласа, получаем равенство
Аа/ - Y^ Aak A,
k=1
k
A
i-i-1
I - E Aak- Ak z,eAt dA,
k=ni /
где можно взять 7 = {А € С : Ие А = а} при достаточно большом а € К. Таким образом, получили вид решения. Уточним некоторые детали. Обозначим
' 1
: k = 1,2
..., n j"
го := (2An) а-°
Y1 := {A G C : |A| = ro, arg A G (-n, n)} , 72 := {A G C : arg A = n, A G [-r0, -<»)} ,
3
73 := {A G C : arg A = -n, A G (-ro, -ro]} , Y := У Yk,
k=1
J (Xai-J2XakAk) E A^-'-Ufc ) extdX, t> 0,
Y \ k=1 / \ k=ni /
при l = 0, 1, .. ., m - 1.
Лемма 1. Пусть Ak G L(Z), k = 1, 2,..., n, z, G Z при некотором l G {0,1,..., m-1}. Тогда функция z(t) = Z,(i)z, является единственным решением задачи (1), (3).
i
Доказательство. Для всех А g 7 имеем |А| > г0 > (2п ||Дк||) , к = 1,2,...,n, поэтому |A|ak-a||Ak||l(z) < 1/(2n),
Y,Aak-aAk
k=1
<
L (Z)
Aa / Aak Ak
k=1
<
L (Z)
|A|
E
j=o
Ea°
k=1
Ak
<
L (Z)
|A|a
-1
1
1
2
j
1
2
поэтому существует обратный оператор ( Аа1 — Аак Ak ) Е L (Z).
k=1
Пусть R > го,
4
Гд = [J Гг,я, Г1,д = Yi, Г2,я = {А Е C : | А| = R, arg А Е (п, —п)} ,
i=i
Гэ,я = {А Е C : arg А = п, А Е [—го, — R]} ,
Г4,я = {А Е C : arg А = —п, А Е [—R, —го]} ,
замкнутый контур Гд обходится по часовой стрелке. Введем в рассмотрение также контуры
гб,я = {А Е C : arg А = п,А Е (—R, —го)} , Гв,я = {А Е C : arg А = —п, А Е (—го, —R)} ,
тогда y = Г5,д U Г6,я U Гд \ Г2,д.
При t > 0 определяющий оператор Z;(t) интеграл сходится, перепишем его в следующем виде:
l\I+ 2т
- ni — 1
АаI — ^ Аак Ak ^А'
ak—1—1AkeAt ¿А
k=1
k=1
ni -1
АаI — ^ Аак A k Y Аак—1—1AkeAt ¿А.
k=1
k=1
При t Е [0,1], А Е y имеем следующие неравенства:
ni -1
АаI — ^ Аак Aj А'
ak—1—1Ak eAt
k=1
k=1
ni -1
< 2^ IА |
k=1
L (Z)
ak— а—l —1
IA ||
sro <
2Anero
- |A|1+a '
Здесь использован тот факт, что при к € {1, 2,..., п; — 1} выполняется тк — 1 < I, поэтому ак < I. Отсюда при j € {0,1,..., т — 1}
\ -1,
ni -1
Аа1 — Аак Afc Y Аак—1—1+jA
„At
k e
k=1
k=1
<
2Anero
— |А|2+а—m '
L (Z)
при этом 2 + а — m = 1+ а — (m — 1) > 1. Следовательно,
// n \ 1 ni —1
Аа1 — Y Аак AH Y Аак —1—1+j Ak eAt ¿А
Y V k=1 / k=1
L(Z)
< 2Anero | I — I + I +
Гя Г2,я Гб,я Гб,я
ds
|А|
2+а—m
-1
при К ^ го. В силу полученных оценок интегралы сходятся равномерно по 4 € [0,1], поэтому при j = 0,1,..., I
¿-з 1 г( « V1 "р-1
21 № := 77—м 7 + / А"7 " Е АЧ Е <*а, 4 > 0,
(I — j)! 2П
7 4
при j = I + 1,1 + 2,. .., т — 1
к=1
к=1
1
- «г — 1
Аа1 — ^3 Аак АН —1—1+3АкеА* ¿А, 0
к=1
к=1
.(3)(0) = 0 при j = I, .(1)(0) = I. Таким образом, функция удовлетворя-
ет условиям (3).
Далее при I = 0,1,..., т — 1 оценим интегралы, определяющие .(£), по контуру 72:
2А
«г —1
2пг I ^
к=1
+ю
< ^е^ / -
п У г
Го
¿г
а+1+1 — а„
< С1е
— го*
так как а + I + 1 — ап > 1, г0 > 1по определению. На контуре 73 получаем аналогичные неравенства. Для контура 71 имеем
2А П"г — 1
к=1
п
< еГ0*с°8 ^ < 2пС2еГо*.
Отсюда следует экспоненциальная ограниченность оператор-функции . (4), а значит, к ней можно применять преобразование Лапласа: при Ие м > г0
сю
.г(м) = I е-^ 0
(^"Е^^ ) (А"-'-1/ - /е^-^^Ш
7 V к=1 / V к=пг / 0
= ¿/¡Т^А (^"Е*"^) £ А—'-1^ ¿А
с учетом неравенства
к=1
к=пг
1
М — А
/ V к=«г
Следовательно,
ад = 1
к=1
++
^ )
С
- Таг
2пг 7 7 7 У м — А
Гя Г2,я Гб,я Гв,я
ГТ Аа/-£А^
к=1
— 1
х А"-1-1/ - X dA
V k=n; /
( n \ —1 / n N
- X A J Ma-1-1/ - X Mafc-1-1Ak
V k=1 / \ k=n; /
-1
n; —1
-X A J Ak + X A J M—1—1
k=1 / V k=1 k=1 )
/ n \—1 n; —1
= М-1—1/ + -X) a J 1—1Ak •
\ к=1 / к=1
Выражение получено по интегральной формуле Коши для интеграла по Гд при К > когда точка ц лежит внутри этого контура. При этом интегралы по при ] = 2, 5, 6 стремятся к нулю при К ^ го. Тогда
= - X Ак) и-1-1/- X мак-г-1Ак) - ма-г-1/
\ к=1 / \ к=и /
/и и \ / п \ —1
= ( - X Ак + X Ак - X Ак )
к=1 к=1 ) V к=1
ч^ / 1 — 1 т \ ^ 1— 1 л 1 —1—1т-
х м / - М k AU - М /
k=n;
1
хм"к A j ма/ -x) мак a j ма—1—1/ - x Мак—1—1Ak
Х^Г AkZi
k=1
(м).
k=1 \ k=1 / \ k=n;
nnn
- X Мак—1—1Ak = XМакAkZ,(m)- X Мак—1—1Ak = Lap
k = n; k=1 k = n;
Поэтому
n
DtaZ,(i) = XAkZ,(t) при t> 0.
k=1
В то же время при таком к0, что mk0 - 1 < l,
DT°Z,(M) = Mafc0 - XMak—(м—1—1/- X Mak—a—1—1Ak)
V k=1 / \ k=n; /
/ n \—1 n; —1
= мако—1—1/ + /-Xмак—"AJ X мак—a—1—1+afco Ak, k=1 k=1
+ I ) £ !iak-a-l-1+ak° Аке^ dti.
2пг „
Y \ k=1 / k=1
Поэтому при ako < l а в случае ako = l Пусть теперь mj0 — 1 > l, тогда
lim D"'0 Zi(t) = 0, t^o+ t
lim D"'0 Z (t) = I.
D"'0 Z(/) = /"'o I — £ /"k—"Ak м—1—1I — £ Mak—a—1—1Ak
- //"k0
k=1 k=ni
/ n \ —1 ni — 1
"'o —1 —1r = I r _ V^/("k — a ,
I = I — ^ —"Ak £ —a—1—1+ako A
k
k=1 k=1
1 / n \ 1 ni — 1
пг Y V k=1 J k=1
Имеем ak — а — l — 1 + ak0 < —1, так как ak0 — a< 0 и ak — l < 0 при k = 1, 2,..., n; — 1, поэтому lim D"Z;(t) = 0. Таким образом, в силу ограниченности
операторов Ak, k = 1, 2,..., n, Z; (-)z; удовлетворяет уравнению (1) и при t = 0. Тем самым Z;(-)z; — решение задачи (1), (3).
Пусть zi(t) и ^(i) — два решения задачи (1), (3) на R+. Зафиксируем T > 0, тогда y(t) = z1(t) — z2(t) — решение задачи Коши z(l) (0) = 0, l = 0,1,... ,m — 1, для уравнения (1) на отрезке [0,T]. Доопределим функцию y нулем при t Е (T, +го). Полученная функция ограничена и также является решением этой задачи Коши для уравнения (1), кроме, может быть, t = T. Подействуем преобразованием Лапласа на обе части равенства
n
D"y(t) = £ D"k Aky(t)
k=1
и получим
D"y(/) = /"£(/) = Lap
Поэтому
£Dt"k Aky
k=1
(/) = £ mu"'Ak£(/).
k=1
/"I — £ Aj = 0
k=1
1 ni -1
1
k
и при Ие^ > 0 имеем у(^) = 0. Отсюда — г2(4) = у(4) = 0 при 4 € [0,Т). Так как Т > 0 можно выбрать сколь угодно большим, г^) = г2(4) при всех
I > 0. □
Из леммы 1 и линейности уравнения (1) сразу получим следующее утверждение.
Теорема 1. Пусть А € ^ (^), к = 1, 2, ...,п, гг € ^, I = 0,1,...,т — 1. Тогда функция
т— 1
г (4) = X ЭДг,
г=о
является единственным решением задачи (1), (2).
3. Неоднородное уравнение
Рассмотрим неоднородное уравнение
п
^+ /(4), I € [0,Т]. (4)
к=1
Обозначим
Лемма 2. Пусть А € ^ ), к = 1,2,...,п, / € С([0,Т]; ^). Тогда функция
t
г/ (¿) = | Я(4 — (а) (5)
о
ением задачи
(1) (0) = 0, I = 0,1,..., т — 1, (6)
о
является единственным решением задачи
г
для уравнения (4).
Доказательство. Имеем при а > 1, I = 0,1,..., т — 2, А € 7
А1 ( Аа/ — X Аак Ак )
^ I _
к=1 / при этом а — I > 1, поэтому для всех 4 > 0
е
2
<
|А|а—1:
^ (^)
/е Не А ^^ " С'2еГ0*'
7
Я (1)(0)=0, г/р) (4) ^ У Я(р) (4 — 5)/(а) ¿з, р = 0, 1,...,т — 1.
При I = т — 1 имеем
!■ е4 Не А !■ е4 Не А
П^-^Ьс*) <с] у ^^ ^
Г 71
сю сю
Г "теА \ <Т{т-а)Г-гп, ] = 1,2,
— I —т+ I — I — т+ I — 4 ' то т 1
у |А|а—т+1 — у ^а—т+1 — у ра—т+1
Yj г0 о
поэтому для всех 4 > 0
||Я(т—(*) < СГ—теГо4, в частности, ||Я(т—1)(4)|^(^) = 0(4а—т) при 4 ^ 0+. Таким образом,
4
7 J а — т + 1
при 4 ^ 0+ и условия Коши (6) для функции г/ вида (5) выполняются.
Положим /(4) = 0 при 4 > Т, тогда г = Я * /, г = Как при доказательстве леммы 1, нетрудно показать, что
ад= — Е Мак А^ , (7)
так как для А € 7
1п
М _ А .
^ \ к=1
С
<
- |А|а+1
Следовательно,
¿^г(м) = (м) = (v/ — /(м)
к=1
= /(м)+£МакАл Ма1 — еМакАк /(м) = Дм)+еАкг(м). (8)
к=1 V к=1 ) к=1
Подействовав обратным преобразованием Лапласа на обе части этого равенства, получим при I € (0,Т] равенство (4), так как / € С([0,Т]; X). Обозначим
—1
а Л еА ¿А, 4> 0,
тогда при к = 1, 2,..., п
1 Г е4 Не А ¿в
1
—1
В силу (8)
ЯГ г (¿) = | .ак (4 — я)/(5) 0
поэтому
а — ап
при 4 ^ 0+. В силу ограниченности операторов , к = 1, 2,... ,п, равенство (4) выполняется и при 4 = 0 в предельном смысле.
Доказательство единственности такое же, как для леммы 1. □ Теперь можно сформулировать общий результат.
Теорема 2. Пусть € ^), к = 1, 2,..., п, / € С([0,Т]; ^), € ^, I = 0,1,..., т — 1. Тогда функция
т— 1
т—1 л
= ^ .(ф, + .(4 — 5)/(я) ¿я 1=0 ^
является единственным решением задачи (2), (4).
Замечание 1. Заметим, что при I = 0,1,..., т — 1
п
ад = (¿) — £ (¿)Ак
к=пг
Действительно, из (7) следует, что при Ие м > г0
Ьар
(м)
к=пг
п \ —1 ( п \
ма1 — £макЛ) ма—1 — £ мак—= .(м).
к=1 / \ к=пг
4. Вырожденное уравнение
Пусть X, & — банаховы пространства, ^ (X; &) — пространство линейных непрерывных операторов, действующих из X в &, С 1(Х; &) — множество всех линейных замкнутых операторов с областью определения, плотной в X, действующих в &. Будем предполагать, что
п € М, Ь,МьМ2,...,Мп— 1 € ^ (X; &), кег Ь = {0}, Мп € СЧ^; &),
Ом„ — область определения оператора Мп, на которой задана норма графика
НЬм„ := II ■ 11х + ||Мп -|к.
Обозначим
рь(М«) := {м € С : (мЬ — М«) —1 € ^; X)},
Я£(М„) := (МЬ - Мп)-1Ь, Ь^(М„) := Ь(МЬ - М„)-1. Оператор Мп называется (Ь, ст)-ограниченным, если
За > 0 V € С (|м| > а) ^ (м € РЬ(М„)). В случае (Ь, <г)-ограниченности оператора Мп определим проекторы
7 7
где 7 := {м € С : |м| = г > а} (см. [28, с. 89, 90]). Положим
X0 := кег Р, X1 := 1т Р, := кег д, &1 := 1т Q.
Обозначим для краткости Р° := I — Р, д° := I — д, через Ьг (Мк,г) — сужение оператора Ь (Мк) на Xг (ДМп г := ПХг при к = п), г = 0,1, к = 1, 2,..., п. При этом известно (см. [28, с. 90, 91]), что ЬР = дЬ, МпРх = дМпх для ж € Дм„, поэтому
Мп,1 € &(X1; &1), Мп,° € С1(Х0; Ьг € &(Xг; ), г = 0,1.
Кроме того, в этой ситуации существуют операторы М- € &X°), Ь-1 € &(&1; X1). Оператор Мп будем называть (Ь, 0)-ограниченным, если Ь° — нулевой оператор.
Рассмотрим начальную задачу
х(1)(0) = XI, I = 0,1,... ,тп — 1, (Рж)(1)(0) = XI, I = тп,тп + 1,.. .,т — 1, (9)
для линейного неоднородного уравнения дробного порядка
п
Я^Ьх^) = £ ^ Мк х(*) + дф, (10)
к=1
которое называется вырожденным в случае кегЬ = {0}. Предполагается, что, как и прежде, 0 < а1 < а2 < ■ ■ ■ < ап < а, т = [а], тк = |~ак], к = 1, 2,..., п, д € ([0,Т]; &).
Решением задачи (9), (10) будем называть функцию х : [0,Т] ^ , для которой х € Ст"-1([0,Т]; X), Д^Ьх, Мк х € С ([0,Т]; &), к = 1, 2 ...,п, выполняются равенства (10) при всех 4 € [0,Т] и (9).
В определении решения используется тот факт, что при условии (Ь, 0)-огра-ниченности оператора Мп гладкость функции Рх такая же, как у функции Ьх, поскольку Рх = Ь-1Ьх, Ьх = ЬРх. При этом условие Д^Ьх € С([0,Т]; X) в силу определения используемой здесь дробной производной Герасимова — Ка-путо означает, в частности, что Ьх € Ст-1([0, Т]; X). Тем самым все начальные условия (9) на функции из класса решений имеют смысл.
Теорема 3. Пусть Ь, Мк € ^(X; У), к = 1, 2,..., п — 1, Мп € ; У)
(Ь, 0)-ограничен, МкР = дМк, к = 1, 2,..., п - 1, д € С([0, Т]; У), € X при 1 = 0,1,..., тп — 1, ж; € X1 при 1 = тп, тп + 1,..., т — 1. Тогда существует единственное решение задачи (9), (10).
Доказательство. Нетрудно показать, что из условий МкР = фМк при к =1, 2,..., п — 1 сразу следует, что
Мк,г € ^(Xr; Уг), г = 0,1, к = 1, 2,..., п — 1.
Подействуем на (10) оператором Ь-1^ € ^(У1; X1) и получим уравнение
п
^ад) = ^ ^ ь-1Мк^(4)+ь-^), к=1
где = Рж(£). Если же аналогичным образом используем оператор М— ¿^о € ^(У0; X0), то получится уравнение
п- 1
= — х ^ Мп-0Мк^(4) — м—ододф, (11)
к=1
= Р0ж(4). При этом использованы равенства
МкРо = Мк — дМк = £0Мк, к =1, 2,..., п — 1. (12)
Уравнения (11) и (12) снабжены начальными условиями
«(;)(0) = Рж;, 1 = 0,1,...,т — 1, (13)
ад(;)(0) = Р0жг, 1 = 0,1,...,тп — 1. (14)
По теореме 2 каждая из задач (11), (13) и (12), (14) имеет единственное решение, при этом оно имеет вид соответственно
г
V
т —1 „
(¿) = ^ ХгдС^Ржг + ад — я)Ь—^(я) &
7_П
;=0 0
и
тп —1 г
^(¿) = — X ад^х* — / ад — а)М—0^0д(«) &
7_П
1=0
где для 4 > 0
ад) / (а^-ХА-^Мм
2П „ ,
у ^ к=1
да-11 — ^ да,—1Ь—1Мк, Л еАг ¿А, 1 = 0,1,..., т — 1,
к=пг /
-1
-1
^ Ь-1ММ ) еА ¿л,
^ V к=1 /
™ V к=1 /
X (ла--;-11 — £ лак-г-1М-°Мм) еА ¿л, 1 = 0,1,...,тп — 1,
V к=(п-1)г )
Х0(г) := [ (ха"1 - £ Х^М-^МьЛ еА* ¿А,
контуры 71, 7° строятся так же, как контур 7 для невырожденного случая, но с использованием в качестве радиуса окружности константы
год := ( 2п тах | , \\Ь1 : к = 1,2,...
(Г 1 И 1 и 1 \
2(п - 1) тах | 2{п — 1)' : Л = 1, 2, .. ., П - 1 И
соответственно. □
Замечание 2. Из доказательства теоремы 3 следует, что задача Коши
х(;) (0) = х;, 1 = 0,1,..., т — 1,
для уравнения (10) разрешима лишь при выполнении дополнительных условий
Р°хт„+г = ад(тп+г) (0), 1 = 0, 1,.. ., т — тп — 1,
связывающих между собой начальные данные задачи.
Замечание 3. Нетрудно заметить, что если оператор Мп (Ь, 0)-ограничен, условия (Рх)(;)(0) = х; € X1 при 1 = тп,тп + 1,...,т — 1 эквивалентны условиям (Ьх)(;) (0) = = Ьх; € &1.
5. Приложение к начально-краевым задачам
Пусть заданы многочлены
Рр(л) = Е^л*, д^(л) = 5, ¿=° ¿=°
с^-, ¿к € С, j = 0,1,... ,р, к =1, 2,..., п, ср = 0, О С М^ — ограниченная область с гладкой границей дО,
Ы<2г 1 2 ... Л
д
(В1п)(з) = X Ь1чсС°°(Щ, 1 = 1, 2,..., г,
|д|<П 1 2 " ' 0
9 = (91, 92,''' , 90) € N0, Ы = 91 +----+ 90, операторный пучок Л, Вь В2,... , Рг
регулярно эллиптичен [29]. Пусть оператор Л1 € С1(Р2(О)) с областью определения
= Н{2Вг} (О) := {V € Я2г(О) : В«(в) =0, I = 1, 2,... , г, в € дО}
действует согласно равенству Л1и = Ли. Предположим, что Л1 — самосопряженный оператор, тогда спектр ст(Л1) оператора Л1 действительный и дискретный [29]. Пусть, кроме того, спектр ст(Л1) ограничен справа и не содержит нуля, : к € М} — ортонормированная в Р2(О) система собственных функций оператора Л1, занумерованных по невозрастанию соответствующих собственных значений {Лк : к € М} с учетом их кратности. Рассмотрим начально-краевую задачу
—г(з,0)=щ(з), 1 = 0,1,...,т-1, вбО, (15)
дг
ВгЛки(в,4) = 0, к = 0,1,.",« - 1, I = 1, 2,.", г, (в, 4) € дО х [0,Т ], (16)
П
ДаРр(Л)и(я,*)=^ дк(Л)и(М)+ Л(я,*), (5,4) € О х [0,Т], (17)
к=1
где 0 < а1 < а2 < ■ ■ ■ < ап < а, к : О х [0, Т] ^ М. Возьмем X = {V € Я2гр(О) : РгЛк ф) = 0, к = 0,1, - . ,р - 1,1 = 1, 2,... ,г, в € дО}, У = Р2(О), Р = Рр(Л) € & (X; У), Мк = ф^Л) € & (X; У), к = 1, 2,.", п.
Пусть Рр(Лк) = 0 для всех к € М, тогда существует обратный оператор Р-1 € & (У; X) и задача (15)-(17) представима в виде задачи (2), (4), где % = X, = Р-1Мк € &(%), к = 1,2, .",п, ,гг = иг(-), 1 = 0,1,...,™ - 1, f (4) = Р-1к(-,£). По теореме 2 существует единственное решение задачи (15)-(17) при любых и € X, I = 0,1,..., т — 1, и к € С([0, Т]; Р2(О)) (в таком случае Р-1к € С([0,Т]; X)).
Пример 1. Возьмем
Р2 (Л) = Л2, Ф2(Л)= а, ^2(Л) = Ьо + &1Л, й =1,
д2и
0=(0,тг), г= 1, Аи=—, В1=1, «1 = 1/4, а2 = 4/3, а = 5/2.
дв2
Тогда т = 3 и задача (15)-(17) имеет вид
= а01/4и(з,1) + (ъ0 + Ьх^ (я,*) € (0,тг) х Ж+,
и(0,1) = и(тг, I) = 1) = ^—(тг, I) = 0, í е м+,
дв2 дв2
ди д2 и
и(в,0) =и0 (я), -тгг М) =И1(в), —(в,0) =и2(я), 5 6(0,7Г). д£ д42
Теперь рассмотрим вырожденный случай. Предположим, что Рр(Ак) = 0
при некотором к € М, тогда при условии, что многочлены Рр и Я не имеют
общих корней на множестве (А^}, оператор Мп (Ь, 0)-ограничен (см. [30]), при
этом проекторы имеют вид
Р = £ , Я = £ ,
РР(Ак )=0 Рр(Ак )=0
где > — скалярное произведение в Ь2(0). Начальные условия с учетом замечания 3 зададим в виде
д1 и
тгп-(я,0) = гф), г = 0,1,...,то„-1, в е П,
, ^ (18) д'Р (Л)и
--(в, 0) = г = то„,то„ + 1,...,то-1, в е О,
дг
тогда задача (16)-(18) представима в виде (9), (10) с выбранными выше пространствами X, ^ и операторами Ь, М. Из теоремы 3 следует однозначная разрешимость задачи (16)-(18) при любых начальных данных и; € X, I = 0,1,..., тп — 1, у; € Ь2(0), I = тп, тп + 1,..., т — 1, таких, что (у;, ^> = 0 при Рр(Ак) = 0, и Н € С([0,Т]; ¿2(П)). Пример 2. Пусть
Р2(А) = А(А + 9), п = 2, Я2(А) = а, Я2(А) = 1 + А, й = 1,
д2и
0=(0,тг), г = 1, Аи=—, В1=1, «1 = 1/4, а2 = 4/3, а = 5/2.
дв2
Тогда т = 3, т2 = 2 и задача (16)-(18) имеет вид
= aZ)t1/4ii(s,i) + Z?t4/3 (bo u(s>*)> (М) е (0, тг) х R+,
d2u. . d2u, . —
u 0, t) = и тг, t) = —5-(0, t) = — (тг, i) = 0, i G R+, ds2 ds2
д2 / д4 д2 \
u(s,0) = «o(s), —(s,0) =wi(s), ^ f — +9^2 )u(s,0) = j/2(s), sG (0,тг).
ЛИТЕРАТУРА
1. Учайкин В. В. Метод дробных производных. Ульяновск: Изд-во Артишок, 2008.
2. Tarasov V. E. Fractional dynamics: Applications of fractional calculus to dynamics of particles, fields and media. New York: Springer, 2011.
3. Самко С. Г., Килбас А. А., Маричев О. И. Интегралы и производные дробного порядка и некоторые их приложения. Минск: Наука и техника, 1987.
4. Pruss J. Evolutionary integral equations and applications. Basel: Birkhauser-Verl., 1993.
5. Podlubny I. Fractional differential equations. Boston: Acad. Press, 1999.
6. Kiryakova V. Generalized fractional calculus and applications. Harlow: Longman Scientific & Technical, 1994; copublished in New York: John Wiley & Sons, Inc.
7. Bajlekova E. G. Fractional evolution equations in Banach spaces: PhD thesis. Eindhoven: Eindhoven Univ. Technol., Univ. Press Facilities, 2001.
8. Псху А. В. Уравнения в частных производных дробного порядка. М.: Наука, 2005.
9. Kilbas A. A., Srivastava H. M., Trujillo J. J. Theory and applications of fractional differential equations. Amsterdam; Boston; Heidelberg: Elsevier Sci. Publ., 2006.
10. Kostic M. Abstract Volterra integro-differential equations. Boca Raton: CRC Press, 2015.
11. Jiang H., Liu F., Turner I., Burrage K. Analitical solutions for the multi-term time-space Caputo—Riesz fractional advection-diffussion equations on a finite domain //J. Math. Anal.
Appl. 2012. V. 389, N 2. P. 1117-1127.
12. Liu F., Meerschaert M. M., McGough R. J., Zhuang P., Liu Q. Numerical methods for solving the multi-term time-fractional wave-diffussion equation // Fract. Calc. Appl. Anal. 2013.
V. 16, N 1. P. 9-25.
13. Alvarez-Pardo E., Lizama C. Mild solutions for multi-term time-fractional differential equations with nonlocal initial conditions // Electron. J. Differ. Equ. 2014. V. 2014, N 39. P. 1-10.
14. Singh V., Pandey D. N. Existence results for multi-term time-fractional impulsive differential equations with fractional order boundary conditions // Malaya J. Math. 2017. V. 5, N 4. P. 625-635.
15. Singh V., Pandey D. N. Mild solutions for multi-term time-fractional impulsive differential systems // Nonlinear Dyn. Syst. Theory. 2018. V. 18, N 3. P. 307-318.
16. Глушак А. В. Задача типа Коши для абстрактного дифференциального уравнения с дробными производными // Мат. заметки. 2005. Т. 77, № 1. C. 28-41.
17. Lizama C., Prado H. Fractional relaxation equations on Banach spaces // Appl. Math. Lett. 2010. V. 23. P. 137-142.
18. Karczewska A., Lizama C. Solutions to stochastic fractional oscillation equations // Appl. Math. Lett. 2010. V. 23. P. 1361-1366.
19. Li C.-G., Kostic M., Li M. Abstract multi-term fractional differential equations // Kragujevac J. Math. 2014. V. 38, N 1. P. 51-71.
20. Fedorov V. E., Kostic M. On a class of abstract degenerate multi-term fractional differential equations in locally convex spaces // Euras. Math. J. 2018. V. 9, N 3. P. 33-57.
21. Плеханова М. В. Задачи стартового управления для эволюционных уравнений дробного порядка // Челяб. физ.-мат. журн. 2016. Т. 1, № 3. С. 15-36.
22. Федоров В. Е., Плеханова М. В., Нажимов Р. Р. Линейные вырожденные эволюционные уравнения с дробной производной Римана — Лиувилля // Сиб. мат. журн. 2018. Т. 59, № 1. С. 171-184.
23. Байбулатова Г. Д. Задачи стартового управления для одного класса вырожденных уравнений с младшими дробными производными // Челяб. физ.-мат. журн. 2020. Т. 5, вып. 3. С. 271-284.
24. Федоров В. Е., Авилович А. С. Задача типа Коши для вырожденного уравнения с производной Римана — Лиувилля в секториальном случае // Сиб. мат. журн. 2019. Т. 60, № 2. С. 461-477.
25. Федоров В. Е., Гордиевских Д. М., Балеану Д., Таш К. Критерий приближенной управляемости одного класса вырожденных распределенных систем с производной Римана — Лиувилля // Мат. заметки СВФУ. 2019. Т. 26, № 2. С. 41-59.
26. Федоров В. Е., Нагуманова А. В. Линейные обратные задачи для вырожденного эволюционного уравнения с производной Герасимова — Капуто в секториальном случае // Мат. заметки СВФУ. 2020. Т. 27, № 2. С. 54-76.
27. Федоров В. Е., Фуонг Т. Д., Киен Б. Т., Бойко К. В., Ижбердеева Е. М. Один класс полулинейных уравнений распределенного порядка в банаховых пространствах // Челяб. физ.-мат. журн. 2020. Т. 5, № 3. С. 342-351.
28. Sviridyuk G. A., Fedorov V. E. Linear Sobolev type equations and degenerate semigroups of operators. Utrecht; Boston: VSP, 2003.
29. Трибель Х. Теория интерполяции. Функциональные пространства. Дифференциальные операторы. М.: Мир, 1980.
30. Федоров В. Е. Сильно голоморфные группы линейных уравнений соболевского типа в локально выпуклых пространствах // Дифференц. уравнения. 2004. Т. 40, № 5. С. 702— 712.
Поступила в редакцию 25 .марта 2021 г. После доработки 25 .марта 2021 г. Принята к публикации 26 августа 2021 г.
Федоров Владимир Евгеньевич
Челябинский государственный университет,
кафедра математического анализа,
ул. Бр. Кашириных, 129, Челябинск 454001;
Южно-Уральский государственный университет,
лаборатория функциональных материалов,
пр. Ленина, 76, Челябинск 454080
Бойко Ксения Владимировна Челябинский государственный университет, кафедра математического анализа, ул. Бр. Кашириных, 129, Челябинск 454021 kvboyko@mail. ru
Та Дуй Фуонг
Институт математики Вьетнамской академии наук и технологий, отделение численного анализа и научного вычисления, Ханой, Вьетнам [email protected]
Математические заметки СВФУ Июль—сентябрь, 2021. Том 28, № 3
UDC 517.9
INITIAL VALUE PROBLEMS FOR SOME CLASSES OF LINEAR EVOLUTION EQUATIONS WITH SEVERAL FRACTIONAL DERIVATIVES V. E. Fedorov, K. V. Boyko, and T. D. Phuong
Abstract: The problems of unique solvability of initial problems for linear inhomoge-neous equations of a general form with several Gerasimov—Caputo fractional derivatives in Banach spaces are investigated. The Cauchy problem is considered for an equation solved with respect to the highest fractional derivative containing bounded operators at the lowest derivatives. The solution is presented with the use of Dunford—Taylor type integrals. The obtained result allowed us to study an initial problem for a linear inhomogeneous equation with a degenerate operator at the highest fractional derivative, provided that, with respect to this operator, the operator at the second largest derivative is 0-bounded. Abstract results are applied to the study of a class of initial-boundary value problems for equations with several Gerasimov—Caputo time derivatives and with polynomials with respect to a self-adjoint elliptic differential operator in space variables.
DOI: 10.25587/SVFU.2021.75.46.006
Keywords: fractional order differential equation, Gerasimov—Caputo fractional derivative, degenerate evolution equation, Cauchy problem, initial-boundary value problem.
REFERENCES
1. Uchaykin V. V., Method of Fractional Derivatives [in Russian], Artishok, Ul'yanovsk (2008).
2. Tarasov V. E., Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media, Springer, New York (2011).
3. Samko S. G., Kilbas A. A., and Marichev O. I., Fractional Integrals and Derivatives, Theory and Applications, Gordon and Breach Sci. Publ., Yveron (1993).
4. Pruss J., Evolutionary Integral Equations and Applications, Birkhauser-Verl., Basel (1993).
5. Podlubny I., Fractional Differential Equations, Acad. Press, Boston (1999).
6. Kiryakova V., Generalized Fractional Calculus and Applications, Longman Scientific & Technical, Harlow (1994); copublished in John Wiley & Sons Inc., New York.
7. Bajlekova E. G., Fractional Evolution Equations in Banach Spaces, PhD Thes., Eindhoven Univ. Technol., Univ. Press Facilities, Eindhoven (2001).
8. Pskhu A. V., Partial Differential Equations of Fractional Order [in Russian], Nauka, Moscow (2005).
9. Kilbas A. A., Srivastava H. M., and Trujillo J. J., Theory and Applications of Fractional Differential Equations, Elsevier Sci. Publ., Amsterdam; Boston; Heidelberg (2006).
10. KostiC M., Abstract Volterra Integro-Differential Equations, CRC Press, Boca Raton (2015).
11. Jiang H., Liu F., Turner I., Burrage K., "Analitical solutions for the multi-term time-space Caputo—Riesz fractional advection-diffussion equations on a finite domain," J. Math. Anal. Appl., 389, No. 2, 1117-1127 (2012).
12. Liu F., Meerschaert M. M., McGough R. J., Zhuang P., and Liu Q., "Numerical methods for solving the multi-term time-fractional wave-diffussion equation," Fract. Calc. Appl. Anal., 16, No. 1, 9-25 (2013).
© 2021 V. E. Fedorov, K. V. Boyko, T. D. Phuong
13. Alvarez-Pardo E. and Lizama C., "Mild solutions for multi-term time-fractional differential equations with nonlocal initial conditions," Electron. J. Differ. Equ., 2014, No. 39, 1—10 (2014).
14. Singh V. and Pandey D. N., "Existence results for multi-term time-fractional impulsive differential equations with fractional order boundary conditions," Malaya J. Math., 5, No. 4, 625-635 (2017).
15. Singh V. and Pandey D. N., "Mild solutions for multi-term time-fractional impulsive differential systems," Nonlinear Dyn. Syst. Theory, 18, No. 3, 307-318 (2018).
16. Glushak A. V., "A Cauchy-type problem for an abstract differential equation with fractional derivatives," Math. Notes, 77, No. 1, 26-38 (2005).
17. Lizama C. and Prado H., "Fractional relaxation equations on Banach spaces," Appl. Math. Lett., 23, 137-142 (2010).
18. Karczewska A. and Lizama C., "Solutions to stochastic fractional oscillation equations," Appl. Math. Lett., 2, 1361-1366 (2010).
19. Li C.-G., KostiC M., and Li M., "Abstract multi-term fractional differential equations," Kragu-jevac J. Math., 38, No. 1, 51-71 (2014).
20. Fedorov V. E. and KostiC M., "On a class of abstract degenerate multi-term fractional differential equations in locally convex spaces," Euras. Math. J., 9, No. 3, 33-57 (2018).
21. Plehanova M. V., "Start control problems for fractional order evolution equations [in Russian]," Chelyab. Fiz.-Mat. Zhurn., 1, No. 3, 15-36 (2016).
22. Fedorov V. E., Plekhanova M. V., and Nazhimov R. R., "Degenerate linear evolution equations with the Riemann-Liouville fractional derivative," Sib. Math. J., 59, No. 1, 136-146 (2018).
23. Baybulatova G. D., "Start control problem for a class of degenerate equations with lower order fractional derivatives [in Russian]," Chelyab. Fiz.-Mat. Zhurn., 5, No. 3 271-284 (2020).
24. Fedorov V. E. and Avilovich A. S., "A Cauchy type problem for a degenerate equation with the Riemann-Liouville derivative in the sectorial case," Sib. Math. J., 60, No. 2, 359-372 (2019).
25. Fedorov V. E., Gordievskih D. M., Baleanu D., and Tash K., "Criterion of the approximate controllability of a class of degenerate distributed systems with the Riemann-Liouville derivative [in Russian]," Mat. Zametki SVFU, 26, No. 2, 41-59 (2019).
26. Fedorov V. E. and Nagumanova A. V., "Linear inverse problems for degenerate evolution equations with the Gerasimov-Caputo derivatives in sectorial case [in Russian]," Mat. Zametki SVFU, 27, No. 2, 54-76 (2020).
27. Fedorov V. E., Phuong T. D., Kien B. T, Boyko K. V., and Izhberdeeva E. M., "A class of distributed order semilinear equations in Banach spaces [in Russian]," Chelyab. Fiz.-Mat. Zhurn., 5, No. 3, 342-351 (2020).
28. Sviridyuk G. A. and Fedorov V. E., Linear Sobolev Type Equations and Degenerate Semigroups of Operators, VSP, Utrecht; Boston (2003).
29. Triebel H., Interpolation Theory, Function Spaces, Differential Operators, North-Holland Publ. Co., Amsterdam; New York; Oxford (1978).
30. Fedorov V. E., "Strongly holomorphic groups of linear equations of sobolev type in locally
convex spaces," Differ. Equations, 40, No. 5, 753—765 (2004).
Submitted March 25, 2021 Revised March 25, 2021 Accepted August 26, 2021
Vladimir E. Fedorov, Chelyabinsk State University, Mathematical Analysis Department,
129 Brothers Kashirin Street, Chelyabinsk 454021, Russia;
South Ural State University,
Laboratory of Functional Materials,
76 Lenin Avenue, Chelyabinsk 454080, Russia
Kseniya V. Boyko, Chelyabinsk State University, Mathematical Analysis Department,
129 Brothers Kashirin Street, Chelyabinsk 454021, Russia [email protected] Ta D. Phuong,
Institute of Mathematics of the Vietnamese Academy of Science and Technology, Department of Numerical Analysis and Scientific Computing Hanoi, Vietnam [email protected]