CC BY
14
УДК 517.929.2
DOI: 10.14529/mmph170402
СРАВНЕНИЕ НЕСКОЛЬКИХ ПРИЗНАКОВ УСТОЙЧИВОСТИ ЛИНЕЙНЫХ РАЗНОСТНЫХ УРАВНЕНИЙ
1 2 Д.А. Комиссарова', М.М. Кипнис2
1 Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: komissarovada@susu.ru
2Южно-Уральский государственный гуманитарно-педагогический университет, г. Челябинск, Российская Федерация E-mail: mmkipnis@gmai.com
Сравниваются признаки асимптотической устойчивости линейных разностных уравнений авторов публикации и некоторые известные признаки других авторов. Достаточные условия авторов лучше, чем известные ранее признаки. Конкурирующими являются признаки авторов и китайских исследователей. Признаки Кипниса и Комисаровой содержат линейные ограничения на коэффициенты уравнения, а в работе китайских исследователей были найдены нелинейные условия асимптотической устойчивости. Существуют области в пространстве положительных коэффициентов уравнения, устойчивость в которых диагностируется с помощью признаков авторов, но не выявляется признаком китайских исследователей, и наоборот. Показаны области гарантированной устойчивости, которые выявляются различными признаками, на примере уравнения с двумя запаздываниями. Указаны классы разностных уравнений, в которых признаки Кипниса и Комиссаровой заведомо лучше. Доказаны соответствующие теоремы. Приведены примеры, иллюстрирующие возможности применения различных признаков.
Ключевые слова: устойчивость; разностные уравнения.
1. Основные известные результаты
Рассмотрим линейное разностное уравнение порядка k
k
Xn = Xn-1 - Iasxn-s , (!)
s=1
где as î Ш, as > 0 (1 < s < k ). Мы называем уравнение (1) асимптотически устойчивым, если все
его решения стремятся к нулю при n ® ¥ . Цель работы - сравнить силу нескольких признаков устойчивости уравнения (1).
Авторы настоящей статьи опубликовали следующий результат [1].
Теорема 1.1. Если as > 0 (1 < s < k ) и
k a
0 < I-< 1, (2)
s=1 2sin-
2(2s -1)
то уравнение (1) асимптотически устойчиво.
Теорема 1 является многомерным обобщением известного результата Левина и Мэя [2]. Из теоремы 1.1 вытекает следующий результат.
Теорема 1.2 ([1]). Если as > 0 (1 < s < k) и
k ж
0<Isas <p, (3)
s=1 2
то уравнение (1) асимптотически устойчиво.
- ra TT [Al D ITATAm.TV DAJAPTA
ra m Г4 UUP ТГЛА/Г
Теорема 1.2 усиливает результаты работ [3] и [4], в которых вместо рр в неравенстве (3) были числа 1 и 1 +1 соответственно, а также результат статьи Танга и Джианга [5] с числом -у вме-
сто Ж в неравенстве (3).
Теорема 1.3 ([1]). Если as > 0 (1 < s < к)
к Р
0 < У sas < 2к sin-, (4)
s 2(2к -1)' W
то уравнение (1) асимптотически устойчиво.
Признак устойчивости теоремы 1.3 сильнее признака теоремы 1.2, но обладает тем недостатком, что его правая часть зависит от порядка уравнения (1). В то же время признак теоремы 1.3 слабее признака теоремы 1.1. Теорема 1.1 в некотором смысле не может быть улучшена. Следующая теорема подтверждает это высказывание.
Теорема 1.4 ([1]). Пусть As > 0 (1 < s < к), причем существует такое s (1 < s < к), что
Р к a
As > 2sin-. Тогда найдется точка (a1v..,ак) такая, что as > 0 (1 < s < к), 0 < У^ < 1 и
2(2s - 1) s=1 As
уравнение (1) с коэффициентами alv..,ak неустойчиво.
Замечание. Неулучшаемость констант в условии (2) теоремы 1.1 не свидетельствует о необходимости условия (2) устойчивости уравнения (1). Как показано в работе [6], для некоторых уравнений симплекс устойчивости, определенный неравенством (2), заполняет лишь малую часть области устойчивости.
Р
Постоянная — в теореме 1.2 и правая часть неравенства (4) в теореме 1.3 также не могут
быть улучшены (см. [1]). Приведем теперь для сравнения результаты Кука, Дьери и Хартунга. Теорема 1.5 (Кук-Дьери, [3]). Если as > 0 (1 < s < к) и
0< ysas < 1, (5)
s=1
то уравнение (1) асимптотически устойчиво.
Теорема 1.6 (Дьери-Хартунг, [4]). Если as > 0 (1 < s < к) и
к1
0 < X sas < 1 + -, (6)
,=1 е то уравнение (1) асимптотически устойчиво.
1 р
Поскольку 1 < 1 + — < — , теорема 1.2 сильнее теорем 1.5 и 1.6.
е 2
2. Сравнение с работой Танга и Джианга
Результаты работы [5] Танга и Джианга конкурируют с ранее указанными в настоящей статье признаками.
Теорема 2.1 (Танг-Джианг, [5]). Уравнение (1) асимптотически устойчиво, если
1+т (1 -то/4), 1 <то < 2;
к
X sas < <
s=1
1+2m - 5mo2/4, 1/2 <m < i; (7)
3/2+m0 (2-m0)/4, о <m0 < 1/2,
где mo = X
к
as
s=1
Условие (7) выполняется, если 0 < XSas < _, поэтому теорема 2.1 сильнее теорем 1.5 и 1.6.
s=1 2
Сравним теперь результаты Танга-Джианга с результатами работы [1] авторов настоящей статьи.
Пример 1. Рассмотрим уравнение xn = xn-1 - 1,86xn-1 - 0,01xn-12. В этом случае k p
X sas = 1,98 > —. Поэтому признак теоремы 1.2 не может быть применен. Теорема 1.1 также не-
s=1 2
и
Комиссарова Д.А., Кипнис М.М.
Сравнение нескольких признаков устойчивости линейных разностных уравнений
применима, поскольку I-
5= 2 81п -
§
= 1,003 > 1. Рассмотрим признак Танга-Джианга. Соглас-
2(25 -1)
но (7), т0 = 1,87 и 1 + т0 (1 -т0/4) = 1,996 > 1,98. Условие теоремы 2.1 выполнено, данное уравнение асимптотически устойчиво. Таким образом, для примера 1 теорема 2.1 оказалась эффективнее теоремы 1.1.
Рассмотрим разностное уравнение с двумя запаздываниями
Х„ = хп-1 - ахп-т - Ьхп-к , (8)
где а, Ь е М, запаздывания т, к е N , причем т < к .
Рисунок показывает, как теоремы 1.1, 1.2 и 2.1 обнаруживают области устойчивости уравнения (8) при »7 = 1, к = 12 .
Рисунок показывает некоторые области в пространстве коэффициентов уравнения (8), устойчивость в которых гарантируется теоремой 1.1, но не теоремой 2.1, и наоборот. Далее мы намерены показать, что при а1 = 0 признак теоремы 1. 1 сильнее признака Танга-Джианга.
Лемма 2.1. Для любых натуральных чисел 5 > 5,
0
§
2 ( 250 - 1) 81
81П
§
81П-
2(250 -1)
2(25 -1) 25 -1
Доказательство. Ввиду выпуклости
функции 81п х на
ство 8Ш х >
81П а
верно неравен-
Г 1 §
при х е [0;а], 0 < а < — . Поэтому для всех 5 >
Области гарантированной устойчивости уравнения (8) при т = 1, к = 12 . Пунктирная линия - граница симплекса (3), штриховая линия - граница симплекса (2), сплошная линия -граница области Танга-Джианга
5>5
§
2 ( 250 - 1) 81
81П
§
81П-
>
2(250 -1)
§
2 ( 250 - 1) 81
81П
§
2(250 -1)
§
2(25 -1)
25 -1
2(25 -1) Лемма доказана.
Теорема 2.2. Если а1 = 0, то из неравенства (7) следует неравенство (2).
Доказательство. По лемме 2.1 имеем
а5 < ^ (^ 2-к
5=2 281И-
§
< (25 - 1)=т ПI (5 - 2)
5=2
5 = 2
2(25 -1)
Пусть 1 < т0 < 2 . Из (7) следует, что
5=3
2т +1(5 - 2)а5 < 1 + т |1 I.
(9)
Отсюда т + 4т0 - 4 < 0 . Тогда т0 < -2 + 2^2 » 0,83 < 1, противоречие. Остаются два случая.
Случай 1: ~<М0 < 1. Из (7) имеем
к 5
2т +1(5 -2)а5 < 1+2т -^т2.
5=3
Отсюда
а
5
х
а
а
5
k
X(s-2)as < 1 -4mo2. (10)
s=3 4
Ввиду (9) и (10) при любом m0 получаем
Zas 2 5 2 1 -s~P— <m +7-7mo <1.
s=2 2 sin-p--3 6
2(2s -1)
Случай 2: 0 < m0 £ 1. Из (7) имеем
2m0 + X(s -2)as < 2+2m-"4m02.
s=3 2 2 4
Отсюда
Ввиду (9) и (11) получаем
XX(s-2)as <3-2m0 -4m02. (11)
s=3 2 2 4
Ек а„ 1 2 1 1 2 1
-^— <т+1 -т —то =1 <1
—ТО • р 6 6
5-22sln-
2(25 -1)
при любом то > о . Теорема доказана.
Пример 2. Рассмотрим уравнение хп - хп-1 - о,9хп-1 - о,о5хп-3.
к р р
В этом случае X5а5 -1,95 . Так как 1,95 >— и 1,95 > 6sln— »1,854, теоремы 1.2 и 1.3 неприменимы для исследования устойчивости этого уравнения. Поскольку
к а
X-5-~ о,981 < 1, согласно теореме 1.1 данное уравнение асимптотически устойчиво.
5-1 2 sin-
2(25 -1)
Имеем то - о,95 и 1 + 2то -5то2/4 -1,77 < 1,95 , условие (7) не выполняется. Значит, с помощью признака Танга-Джианга нельзя сделать выводов об устойчивости.
Укажем класс уравнений, для которых признаки теорем 1.2 и 1.3 сильнее признака Танга-Джианга.
Теорема 2.3. Если а1 - а2 -... - а1о - о, то из неравенства (7) следует неравенство (3). Доказательство. Предположим, что 1 < то < 2 . Из (7) получаем
11т + X (5 -11)< 1+то -т].
Отсюда то + 4ото - 4 < о. Тогда то < -2о + 2лЯо1 » о,1 < 1, противоречие. Теперь предположим, что 1 < то £ 1. Из (7) получим
k 5
11m + X (s -11) as < 1+2m0 - 4 m0.
s=12 4
Отсюда 5m02 + 36m0 - 4 < 0 . Тогда m0 < 18 344 » 0,11 < 1, противоречие. Остается случай 0 < m0 £ 1. Из (7) имеем
k 3 11
1m0 + X (s -11) as < 2+2 m- 4 m02.
s=12 2 2 4
Комиссарова Д.А., Сравнение нескольких признаков устойчивости
Кипнис М.М. линейных разностных уравнений
Отсюда т2 + 42^о - 6 < 0 • Тогда ¡и0 <-21 + ^447 «0,142.
из (7) получаем
2 4 V 2 J
^эаэ <- + -4(-21+ Т447)(23-^447)« 1,556<-.
Поскольку функция -2 + -2щ -1Щ возрастает на f 0;-2
s=п 2 4' Теорема доказана.
Так как -2 < 2ksin ^((^—1) при любом натуральном к, получаем следствие:
Следствие теоремы 2.3. Если a1 = a2 = ... = a10 = 0, то из неравенства (7) следует неравенство (4).
Теорема 2.3 и следствие из нее фиксируют преимущество любой из теорем 1.1, 1.2, 1.3 перед признаком Танга-Джианга для разностных уравнений с большими запаздываниями. Следующий пример показывает, что теорема 2.3 в некотором смысле неулучшаема. Пример 3. Рассмотрим уравнение xn = xn-1 -1572,4-10-4Xn-10 -0,0007 • 10-4Xn-480 . Это уравнение асимптотически устойчиво, поскольку (см. (7)) щ = 1572,4007• 10-4 и 3 к
- + щ (2-щ )/4 « 15724,3893 • 10-4 > ^ sas = 15724,336 10-4. Но неравенство (3) не выполнено,
2 s=1
к Ж
так как £sas = 15724,336•Ю-4 >-.
s=1 2
Замечание. Признак Танга-Джианга и теорема 1.1 требуют для устойчивости уравнения (1), чтобы сумма коэффициентов, обозначенная посредством щ , была меньше 2. Но следующий пример показывает, что уравнение (1) может быть устойчивым и при щ > 2 . Пример 4. Рассмотрим уравнение
Xn = Xn-1 - (1,9Xn-1 + 0,8Xn-2 + 0,7Xn-3 + 0,6Xn-4 ). (12)
В данном случае щ = 3. Согласно результатам Энестрема-Какейа 1892 г. [7], уравнение Xn + ^к=хasXn-s = 0 устойчиво, если 1 > ax > a2 > ... > ak > 0 . Уравнение (12) удовлетворяет этому
условию, поэтому оно асимптотически устойчиво.
Работа выполнена при поддержке ФГБОУ ВО «Красноярский государственный педагогический университет им. В.П. Астафьева», договор № 16-746.
Литература
1. Kipnis, M.M. A note on explicit stability conditions for autonomous higher order difference equations / M.M. Kipnis, D.A. Komissarova // Journal of Difference Equations and Applications. -2007. - Vol. 13, № 5. - P. 457-461.
2. Levin, S.A note on difference-delay equations / S. Levin, R. May // Theoret. Popul. Biol. - 1976.
- Vol. 9, № 2. - P. 178-187.
3. Cooke, K.L. Numerical approximation of the solutions of delay differential equations on an infinite interval using piesewise constant arguments / K.L. Cooke, I. Gyori // Computers & Mathematics with Applications. - 1994. - Vol. 28. - Issue 1-3. - P. 81-92.
4. Gyori, I. Stability in delay perturbed differential and difference equations / I. Gyori, F. Hartung // Fields Institute Communications. - 2001. - Vol. 29. - P. 181-194.
5. Tang, X.H. Asymptotic behavior of Volterra difference equations / X.H. Tang, Z. Jiang // Journal of Difference Equations and Applications. - 2007. - Vol. 13, № 1. - P. 25-40.
6. Левицкая, И.С. Область устойчивости линейного разностного уравнения с двумя запаздываниями / И.С. Левицкая // Известия Челябинского научного центра. - 2004. - Вып. 3(24). -С.12-16.
7. Прасолов, В В. Многочлены / В.В. Прасолов. - М.: МЦ МНО, 2001. - 336 c.
Поступила в редакцию 4 июля 2017 г.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 4, pp. 13-18
DOI: 10.14529/mmph170402
COMPARISON OF SEVERAL STABILITY CONDITIONS FOR LINEAR DIFFERENCE EQUATIONS
1 2 D.A. Komissarova , M.M. Kipnis
1 South Ural State University, Chelyabinsk, Russian Federation E-mail: komissarovada@susu.ru
2 South Ural State Pedagogical University for Humanities, Chelyabinsk, Russian Federation E-mail: mmkipnis@gmai.com
The paper compares the asymptotic stability criteria for linear difference equations, and some other well-known criteria found by the authors of the publication. Sufficient conditions given by the authors are better than the known ones. Criteria by the authors compete with the criteria of the Chinese researchers. The criteria of Kipnis and Komisarova consist of the linear constraints on the coefficients of the equations, while in the work of the Chinese researchers nonlinear conditions for asymptotic stability were found. There are domains in the space of positive coefficients of an equation, stability of which is revealed by the authors' criteria, but is not detected by the Chinese researchers test, and vice versa. The areas of guaranteed stability, which are revealed by different conditions, are shown for the linear difference equation with two delays. The classes of the difference equations are given in which the Kipnis and Komissarova criteria are certainly better. The corresponding theorems are proved. Examples illustrating the possibilities of using various criteria are given.
Keywords: stability; difference equations.
References
1. Kipnis M.M., Komissarova D.A. A note on explicit stability conditions for autonomous higher order difference equations. J. Difference Equ. Appl., 2007, Vol. 13, no. 5, pp. 457-461. DOI: 10.1080/10236190601132933
2. Levin S., May R. A note on difference-delay equations. Theoret. Popul. Biol, 1976, Vol. 9, no. 2, pp. 178-187. DOI: 10.1016/0040-5809(76)90043-5
3. Cooke K.L., Gyori I. Numerical approximation of the solutions of delay differential equations on an infinite interval using piesewise constant arguments. Comp. Math. Appl., 1994, Vol. 28, Issue 1-3, pp. 81-92. DOI: 10.1016/0898-1221(94)00095-6
4. Gyori I., Hartung F. Stability in delay perturbed differential and difference equations. Fields Institute Communications, 2001, Vol. 29, pp. 181-194. http://www.ams.org/books/fic/029/14 DOI: 10.1090/fic/029
5. Tang X.H., Jiang Z. Asymptotic behavior of Volterra difference equations. J. Difference Equ. Appl., 2007, Vol. 13, no. 1, pp. 25-40. DOI: 10.1080/10236190601008810
6. Levitskaya I.S. Izvestiya Chelyabinskogo nauchnogo tsentra, 2004, no. 3(24), pp. 12-16. (in Russ). http://csc.ac.ru/news/2004_3/2004_3_1_3.zip
7. Prasolov V.V. Polynomials. Springer Berlin Heidelberg, 2004, 301 p. DOI: 10.1007/978-3-64203980-5
Received July 4, 2017
