CC BY
62
Математика. Физика
УДК 517.9
Б01: 10.17277/уе81шк.2015.01.рр.141-147
УСТОЙЧИВОСТЬ РЕШЕНИЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ПРОИЗВОДНЫМИ ДРОБНОГО ПОРЯДКА
В. А. Рязанцев
Кафедра «Высшая и прикладная математика», ФГБОУ ВПО «Пензенский государственный университет», г. Пенза; ryazantsevv@mail.ru
Ключевые слова: дробная производная; логарифмическая норма; операторная экспонента; преобразование Фурье; устойчивость по Ляпунову.
Аннотация: Представлен подход к проблеме устойчивости решений уравнений в частных производных дробного порядка в смысле Римана-Лиувилля. Предлагаемый подход основан на применении преобразования Фурье и исследовании эквивалентных дифференциальных уравнений в спектральной области. В результате получены достаточные критерии устойчивости решений уравнений в частных производных дробных порядков.
Вопрос устойчивости решений дифференциальных уравнений в частных производных представляет значительный интерес для исследователей ввиду большого числа приложений этих уравнений в различных областях науки и техники. Несмотря на большое число работ по данной тематике, указанная проблема по-прежнему далека от окончательного решения. Основным методом исследования устойчивости в настоящее время является прямой метод Ляпунова [1, 2]. Однако этот весьма эффективный и прекрасно разработанный метод не лишен ряда недостатков, в числе которых следует назвать в первую очередь громоздкость и неудобство с вычислительной точки зрения. Описываемый же в данной работе подход отличается простотой и легко может быть реализован в форме вычислительного алгоритма. Кроме того, данный подход позволяет анализировать устойчивость решений уравнений в частных производных дробных порядков, исследования по которым в настоящее время отсутствуют. Эту работу можно рассматривать как обобщение и развитие идей, ранее изложенных в работах [3 - 7]. Рассмотрим следующую задачу Коши:
Я, т п
ди = + В(0и; (1)
у=1£=1
и(^, х) = ио(х), (2)
где А( ],к ^) = а(^к ^); Б^) = Ьрд ^); ] = 1т; к = 1,П; р = й; ц = й; и(^ х) = (ы]_^, х), х),..., и1 (^ х))т; X = (, Х2,..., хп )т. Символом -х Б Х^и обо-
* По материалам доклада на конференции ММТТ-27 (см. Вестник ТГТУ, т. 20, № 4).
значена дробная производная функции и(/, х) в смысле Римана-Лиувилля порядка V у > 0 по переменной х^, вычисляемая по формуле [8]
[v j ] + 1 X
Dv j - , 1 „ d j Г t dt
- xk-r(i-{vj)j]+1 I(X-^*
где Г - гамма-функция; [V у ] - целая часть порядка дробной производной V у,
а {V у } - дробная часть V у, так что V у = [V у ] + {V у}.
Придав тривиальному решению уравнения (1) некоторое начальное возмущение ио( х), рассмотрим возмущенную задачу Коши (1) - (2). Поставим задачу
об определении условий на элементы матриц А(у,к-*(/), В(/) и Vу, при которых
тривиальное решение уравнения (1) устойчиво в смысле Ляпунова.
Исследовать устойчивость будем в банаховом пространстве вектор-функций
/(х) = С/1(х), / (х),..., // (х))т с нормой
112
Xk-ii_rKij fk (X1'-X")2 dX1-dX"
hit, X)W -
При каждом фиксированном значении ? > ^ норма вектор-функции и(/, х) определяется формулой
т/ р+то р+то .2 112
к=1 .)—ток^то\ик&xl...,хпА ^1...(¡х„ .
ди
Пусть задача (1) - (2) имеет решение, вместе со своей производной д? суммируемое с квадратом по пространственным переменным. При этом условии применим к (1) - (2) преобразование Фурье по хь..., хп
и(/, ю) = Зи(/, х) = 1/(2п)п/2 Г"... Г+"и(/, х)е-г(Ю1х1 +"+юпхп)^х1...скп;
•1— то то
и(Г,х) = 3-1и(/,ю) = 1/(2п)п/2 Г+'к Г+'и(/,ю)ег(Ю1х1 + . +юпхп)dю1...dюn,
—ТО * — то
в результате чего, имея в виду, что [9, с. 196]
DVku(t,X)) - (-jfflk)vj U(t
(3)
получим следующую задачу Коши для параметрического дифференциального уравнения в спектральной области:
dU dt
-0(t, ra)U;
U (i0,ffl) = Uo(ffl), (5)
где ю = (ю1,...,ю„) и ф(/,Ю)=^(t)+Z™=1 ZI=1 а{р;!р(*)(-™к)v 1, P=1,z,
q = Ц.
Введем в рассмотрение следующие нормы векторов:
IUI = |ui|2 + |U2|2 +...+\иг\2]12, Ulli = Uli+\U2\ +...+\Ul|.
Имеют место неравенства |U||2 < |U||2 + q=1 • < ЦиЦ^,
p * q
U12=zp=1 Ы2+2ZP, q=1 pp\■ Kl <1 zp=1 Ы2 < U2.
p * q
Следовательно, справедлива двусторонняя оценка
IUI < U 1 <V71|||. (6)
Покажем, что при любом фиксированном ю е Rn траектория U(t, ю) уравнения (4) остается внутри замкнутого шара B[0,||Uо(ю)||] в пространстве R1 с центром в нуле радиусом ||Uо(ю)||. Для доказательства предположим обратное: пусть при каком-нибудь значении ю = ю е Rn траектория U (t, ю) покидает пределы шара B[0,|\Uо(ю)||] в момент времени t = T. При этом предположении представим уравнение (4) следующим образом
dU ю ю
-= Ф1 (co)U + Ф 2(t,co)U, (7)
dt
где Ф1(ю) = Ф(ю) и Ф2(t, ю) = Ф(t - T, ю).
При t > T > to решение уравнения (7) можно представить следующим образом [10]
U (t, ю) = еф (5)(t-T )U (T, ю) + f' еф (5)(t - 5)Ф 2 (5, S)U (s, S)ds. (8)
j t
На основании (8) имеем оценку сверху
\\U (t, ю)|| < |e°1(5)(t-T U (T, 5)| + JT еф1(5)(t - 5)Ф 2 (5, 5)U (5,S)ds
(9)
Далее последовательно оценим сверху каждое из слагаемых, составляющих правую часть неравенства (9). Для первого слагаемого в силу (6) получаем
||еф1 (5№-тр(Т, 5)|| < ||еф1 -т)| р(т, 5)|| < еЛ(Ф(5))^-т)|р(Т, ю)||, (10)
где Л(Ф1 (5)) - логарифмическая норма оператора Ф1 (5), которая в пространстве с нормой р| ¡1 вычисляется по формуле [11, с. 22]
Л(Ф1(5)) = Л(Ф(Т, 5)) = sup
q
Re(q,q (T,5))+ZP=1,kp,q (T,5)l
p * q
(11)
где фр ц (t,ю), р = 1,1, ц = 1,1 - элементы матрицы Ф(^ю).
Пусть для всех юе Яп и t > tо величина Л(Ф(^ю)) ограничена сверху некоторым отрицательным числом -ст
Л(Ф^,ю))<-ст, ст> 0. (12)
Тогда из (10) вытекает, что
||еФ1(5)0-т )р (т, 5)11 < e-o(t-т )р (Т, 5)||. (13)
Для второго слагаемого в силу (6) выполняется следующая оценка
\*те°1(5)(г-5)Ф2(5,5)и(5,5^ < еф1(5)(г-5)Ф2(5,5)и(5,5)& < ||еФ1 (5)(/-^Ф2(55)и(5,5)|| <к <41\Т||еф1(5)(/-4)| ||Ф2(5,5)и5)||& <
* Г,
-a(i - s)
||02(s,5)U (5, ю)| ds.
(14)
Пусть все коэффициенты ар/^-*(/), Ьрд (/) / = 1, да, к = 1, и, р = 1,/, д = 1, /
непрерывны при / >Тогда из структуры оператора Ф2(5,5) следует, что для любого как угодно малого е > 0 найдется такое значение Д/ > 0, что при Т < / < Т + Д/(е,5) будет иметь место неравенство ||Ф2 (5,5)и(5, 5)|| < е|!(5, 5)||. Стало быть, из (14) следует, что
J{еф1 (»)(i-5)ф2 (s, 5)u(s, 5)ds < ел/7J
e
-o(f - s)
||U (s, 5)| ds.
Неравенства (9), (13), (15) служат обоснованием следующей оценки
\\и(/,< е-о(/-т)|и(Т,5)| + е4/Г е-о(/-5)|\и(5,5)||<Ь.
Введем в рассмотрение функцию 6(s) = е
= e-°(i - s)
||U (s, 5)|| . Тогда (16)
(15)
(16)
можно
переписать следующим образом
0(/) < 0(Т) + елГ/6(5)^5. (17)
Применяя к неравенству (17) неравенство Гронуолла-Беллмана и переходя к исходной нотации, получаем
и(/,5)|| < е(-°+е7/)(г -т)|и(Т,5)||. (18)
Вследствие произвольности выбора е> 0 и Т > /о из неравенства (18) имеем окончательно: |и(/,5)|| < |и0(5)||.
Тем самым, приходим к противоречию, в силу произвольности 5 е Яп доказывающему, что траектория и(/, 5) уравнения (4) не выходит за пределы замкнутого шара В[0,||и0(5)||]. Иными словами, |и(/,5)||<|и0(5)! для всех 5е Яп.
Проинтегрировав последнее неравенство с квадратом по ю и применив формулу Планшереля, приходим к оценке ||м(/, х)||2 < ||м0 (х)||2, из которой сразу же следует неравенство ||м(/,х)^ < ||м0(х)\2. Тем самым, доказана следующая теорема.
Теорема 1. Пусть выполнены следующие условия:
1) функции ард)(/), Ьрд (/), / = 1, да, к = 1, и, р = 1,/, д = 1, / непрерывны при / > /0;
2) при всяком фиксированном наборе значений (/,5) е [/0,<»)х Я1 логарифмическая норма Л(Ф(/, 5)), вычисляемая по формуле (11), удовлетворяет неравенству (12).
Тогда тривиальное решение уравнения (1) устойчиво в смысле Ляпунова.
t
В качестве простейшего примера применения теоремы 1 выведем достаточный критерий устойчивости тривиального решения задачи Коши для следующей системы
^ = ап(0 •-< Щ + ап(0 •-< Б^2 и2 + ¿ll(t)Ul,
и (19)
= а21(0 •-< Б^11 и1 + а22 (t) •-< Б^12 и 2 + ¿22 ^ Щ.
. д
В предположениях теоремы 1 применим к (19) преобразование Фурье по переменной х, в результате чего получим уравнение (4) с матрицей Ф(^ ю), определяемой формулой
) + ац^-гаУ11 а12«(- /-юР2 ^ ч а2l(t)(-/ю)21 ¿22(0 + а22^)(-/ю)22, Известно [4, с. 86], что
Ф(/, ю) =
(-/'ю^рц =|ю^рч [со8(^рц!2)-i рц!2Tsgnю],
где sgn ю обозначает функцию знака переменной ю.
Тогда матрицу Ф(^ю) можно представить следующим образом
Ф^, ю) = ((О + арц^^^^^^рц СОБ^рд/2)-iapq (ф^рц ¡ИП(™р^2Tsgnю),
где р = 1,2, ц = 1,2, ¿12 ^) =¿21^) = 0. Условие 2 теоремы 1 вкупе с формулой (11) приводит нас к следующей системе условий:
{¿11(0 + а^фр1 СО8(™П/2) + |a21(t) ^ 21 <-ст, (20)
[¿22 (t) + 022 (t)|22 CQs(nV2^2) + \an (t)| • |ю|'
< -CT,
где с - какое-нибудь положительное число, для всех ю е Я. Таким образом, справедлива теорема.
Теорема 2. Пусть выполнены условия:
1) функции ар ц (0, ¿п^), ¿22^), р = 1,2, ц = 1,2 непрерывны при t > tо;
2) при каждом ю е Я справедливы неравенства (20).
Тогда тривиальное решение системы (20) устойчиво в смысле Ляпунова.
Список литературы
1. Шестаков, А. А. Обобщенный прямой метод Ляпунова для систем с распределенными параметрами / А. А. Шестаков. - М. : Наука, 1990. - 320 с.
2. Сиразетдинов, Т. К. Устойчивость систем с распределенными параметрами / Т. К. Сиразетдинов. - Новосибирск : Наука, Сиб. отд-ние, 1987. - 231 с.
3. Бойков, И. В. Критерии устойчивости решений дифференциальных уравнений в частных производных параболического типа / И. В. Бойков, В. А. Рязан-цев // Журн. Средневолжского мат. общества. - 2012. - Т. 14, № 3. - С. 12 - 20.
4. Бойков, И. В. Устойчивость решений параболических уравнений с дробными производными / И. В. Бойков, В. А. Рязанцев // Изв. высших учеб. заведений. Поволжский регион. Физ.-мат. науки. - 2012. - № 4(24). - С. 84 - 100.
5. Бойков, И. В. О достаточных критериях устойчивости решений дифференциальных уравнений гиперболического типа / И. В. Бойков, В. А. Рязанцев //
Изв. высших учеб. заведений. Поволжский регион. Физ.-мат. науки. - 2013. -№ 2(26). - С. 33 - 49.
6. Бойков, И. В. Устойчивость по Тьюрингу динамических систем, описываемых уравнениями с дробными производными / И. В. Бойков, В. А. Рязанцев // Журн. Средневолжского мат. общества. - 2013. - Т. 15, № 4. - С. 15 - 24.
7. Бойков, И. В. Устойчивость решений параболических уравнений с дробными производными по координатным переменным / И. В. Бойков, В. А. Рязанцев // Аналитические и численные методы моделирования естественнонаучных и социальных проблем : сб. ст. VII Междунар. науч.-техн. конф., г. Пенза, 22 - 25 окт. 2012 г. / Пенз. гос. ун-т. - С. 27 - 31.
8. Самко, С. Г. Интегралы и производные дробного порядка и некоторые их приложения / С. Г. Самко, А. А. Килбас, О. И. Маричев. - Минск : Наука и техника, 1987. - 688 с.
9. Учайкин, В. В. Метод дробных производных / В. В. Учайкин. - Ульяновск : Артишок, 2008. - 512 с.
10. Далецкий, Ю. Л. Устойчивость решений дифференциальных уравнений в банаховом пространстве / Ю. Л. Далецкий, М. Г. Крейн. - М. : Наука, 1970. -526 с.
11. Бойков, И. В. Устойчивость решений дифференциальных уравнений / И. В. Бойков. - Пенза : Изд-во Пенз. гос. ун-та, 2008. - 244 с.
Stability of Solutions of Differential Equations with Fractional Orders
V. A. Ryazantsev
Department "Higher and Applied Mathematics", Penza State University, Penza; ryazantsevv@mail.ru
Keywords: Fourier transformation; fractional derivative; logarithmic norm; Lyapunov stability; operator exponent.
Abstract: The paper deals with one approach to the problem of stability of partial differential equations solutions with fractional derivatives of Riemann-Liouville type. The suggested approach is based on the use of spatial Fourier transformation and the analysis of equivalent differential equations in spectral space.
References
1. Shestakov A.A. Obobshchennyi pryamoi metod Lyapunova dlya sistem s raspredelennymi parametrami (Generalized direct Lyapunov method for systems with distributed parameters), Moscow: Nauka, 1990, 320 p.
2. Sirazetdinov T.K. Ustoichivost' sistem s raspredelennymi parametrami (Stability of systems with distributed parameters), Novosibirsk: Nauka, Sibirskoe otdelenie, 1987, 231 p.
3. Boikov I.V., Ryazantsev V.A. Zhurnal Srednevolzhskogo matematicheskogo obshchestva, 2012, vol. 14, no. 3, pp. 12-20.
4. Boikov I.V., Ryazantsev V.A. Izvestiya vysshikh uchebnykh zavedenii. Povolzhskii region. Fiziko-matematicheskie nauki, 2012, no. 4(24), pp. 84-100.
5. Boikov I.V., Ryazantsev V.A. Izvestiya vysshikh uchebnykh zavedenii. Povolzhskii region. Fiziko-matematicheskie nauki, 2013, no. 2(26), pp. 33-49.
6. Boikov I.V., Ryazantsev V.A. Zhurnal Srednevolzhskogo matematicheskogo obshchestva, 2013, vol. 15, no. 4, pp. 15-24.
7. Boikov I.V., Ryazantsev V.A. Analiticheskie i chislennye metody modelirovaniya estestvennonauchnykh i sotsial'nykh problem (Analytical and numerical methods for modeling natural and socially insured cial problems), Collection of articles VII International Scientific and Technical Conference, Penza, 22-25 October 2012, Penza, pp. 27-31.
8. Samko S.G., Kilbas A.A., Marichev O.I. Integraly i proizvodnye drobnogo poryadka i nekotorye ikh prilozheniya (Integrals and derivatives of fractional order and some of their applications), Minsk: Nauka i tekhnika, 1987, 688 p.
9. Uchaikin V.V. Metod drobnykh proizvodnykh (The method of fractional derivatives), Ul'yanovsk: Artishok, 2008, 512 p.
10. Daletskii Yu.L. Krein M.G. Ustoichivost' reshenii differentsial'nykh uravnenii v banakhovom prostranstve (Stability of solutions of differential equations in Banach space), Moscow: Nauka, 1970, 526 p.
11. Boikov I.V. Ustoichivost' reshenii differentsial'nykh uravnenii (Stability of solutions of differential equations), Penza: Izdatel'stvo Penzenskogo gosudarstvennogo universiteta, 2008, 244 p.
Immunität der Lösungen der Differentialgleichungen mit den Ableitungen der Bruchordnung
Zusammenfassung: Im Artikel wird ein Herangehen an das Problem der Immunität der Lösungen der Gleichungen in den privaten Ableitungen der Bruchordnung im Sinne von Riemann-Loiouville beschrieben. Das angebotene Herangehen ist auf der Anwendung der Transformation von Fourier und der Forschung der äquivalenten Differentialgleichungen in dem spektralen Gebiet gegründet. Es sind die ausreichenden Kriterien der Immunität der Lösungen der Gleichungen in der privaten Ableitungen der Bruchordnungen erhalten.
Stabilité des solutions des équations différentielles avec les dérivées de l'ordre fractionnaire
Résumé: Dans l'article est décrite une approche envers le problème de la stabilité des solutions des équations dans les dérivés partielles de l'ordre fractionnaire du type Riemann-Liouville. L'approche proposée est fondée sur l'emploi des transformations de Fourier et sur l'étude des équations différentielles dans le domaine spectral. Sont reçues les critères exactes de la stabilité des solutions des équations dans les dérivés partielles de l'ordre fractionnaire.
Автор: Рязанцев Владимир Андреевич - ассистент кафедры «Высшая и прикладная математика», ФГБОУ ВПО «Пензенский государственный университет», г. Пенза.
Рецензент: Бойков Илья Владимирович - доктор физико-математических наук, профессор, заведующий кафедрой «Высшая и прикладная математика», ФГБОУ ВПО «Пензенский государственный университет», г. Пенза.
