Анализ устойчивости решений систем обыкновенных дифференциальных уравнений (ОДУ) опирается на методы качественной теории [1, 2] и теории автоматического управления [3]. Компьютеризация системных исследований приводит к вопросу о применении в этой области средств вычислительной техники. Подход к анализу устойчивости на основе численного моделирования был предложен в [4, 5], исследование продолжено, его современное состояние отражено в [6, 7]. Известные методы, если они не относятся к специальным видам ОДУ, используют различные способы вычисления функций Ляпунова [8, 9]. В излагаемой ниже работе подход опирается на разностное (в общем случае приближенное) решение системы, с выполнением которого совмещается численное моделирование устойчивости. Предварительные результаты приведены для линейных и нелинейных ОДУ в [10, 11]. В некоторых разновидностях предлагаемые методы дают возможность аналитического применения, в целом трактуются как компьютерно-ориентированные. Их построение опирается на рекуррентные и аддитивные преобразования разностных методов, для анализа устойчивости используются также функции правой части ОДУ и аналитическое выражение производных этих функций. Цель работы – представить целостные способы анализа устойчивости, ориентированные на компьютерную реализацию, изложить обоснование и указать ограничения. Требуется показать их применимость в различных условиях для нелинейных, автономных и линейных систем, привести примеры и результаты численных экспериментов. Излагаемые признаки могут входить частью в основу сравнительно доступных средств анализа для автоматизированного контроля устойчивости.
В работе ставится цель построить компьютерно-ориентированный метод анализа устойчивости в смысле Ляпунова решения системы обыкновенных дифференциальных уравнений, который имеет следующие отличия от известных аналогов.
1. Метод инвариантен относительно класса дифференциальных уравнений, применим к нелинейным, автономным и линейным системам, при этом позволяет учитывать специфику системы.
2. Предложенный подход не использует построение функций Ляпунова, не опирается на первый и второй методы Ляпунова.
3. Метод включает алгоритмизацию, сравнительно простую программную реализацию и может использоваться для численного моделирования по ходу решения системы.
4. Отличия дают преимущество относительно известных способов компьютерного анализа устойчивости, позволяющее достигать инвариантности достоверной оценки устойчивости в условиях сходимости метода приближенного решения системы.
Цель работы включает математическое обоснование предложенных критериев, выполнение их алгоритмизации. Кроме того, требуется представить развернутый численный эксперимент, иллюстрирующий практическую реализуемость комплекса предложенных средств компьютерного анализа устойчивости.
Исходные соотношения. Рассматривается задача Коши для системы
(1)
где
,
,
Предполагается, что существует , такое, что на полупрямой
выполнены все условия существования и единственности для невозмущенного решения и для каждого его возмущения
,
, с начальным вектором в границах
. Применяются канонические согласованные нормы матрицы и вектора. По умолчанию
, в численных экспериментах используется эвклидова норма. Предполагается, что в области
функция
определена, непрерывна и непрерывно дифференцируема (в t0 – справа), компоненты функции удовлетворяют неравенству:
. (2)
Из (2) следует условие Липшица . В этих условиях определение устойчивости упрощено: решение
устойчиво, если
найдется
, такое, что
влечет
. Решение асимптотически устойчиво, если оно устойчиво и найдется
, такое, что из неравенства
следует
, если
. Метод Эйлера решения задачи (1)
, (3)
включая запись с остаточным членом, на произвольном отрезке , применяется в предположении, что значение
является произвольно фиксированным, при этом индекс i неограниченно растет одновременно с убыванием равномерного шага:
. (4)
В форме с остаточным членом метод (3) примет вид ,
, где i, h из (4),
– остаточный член формулы Тейлора для k-го компонента приближения:
, аналогично для возмущенного решения,
,
,
,
. Предварительно предполагается, что
, (5)
в дальнейшем ограничение (5) будет ослаблено или снято. Исходными являются соотношения ,
, или
,
. (6)
С учетом (2) при ограничении (5) выполняется неравенство:
. (7)
Отсюда
. (8)
Рекуррентное преобразование (6) влечет
,
, (9)
где h соответствует (4), из (6),
. Для дальнейшего приводятся некоторые утверждения из [10, 11].
Лемма 1. В рассматриваемых условиях .
Доказательство. Для
:
,
. Поэтому
,
,
. Из (9) с учетом (8)
, следовательно,
. Отсюда
, тем более,
, и
при
. Лемма доказана.
Следствие 1. В тех же условиях
. (10)
Лемма 2. В рассматриваемых условиях для устойчивости решения задачи (1) необходимо и достаточно, чтобы существовало , такое, что
, при ограничении
выполняется неравенство
. (11)
Для асимптотической устойчивости необходимо и достаточно, чтобы выполнялось предыдущее утверждение и существовало , такое, что
влечет
. (12)
На частичное произведение
при изменении i меняет одновременно все свои сомножители и h в каждом из них согласно (4).
Расширенныe условия устойчивости. Из (10)
. (13)
Следствие 2. При формулировка и утверждение леммы 2 дословно сохраняются при замене (11) на соотношение
, (14)
и (12) – на соотношение
. (15)
Частное в (14) выражает отношение возмущения решения именно к вызвавшему его возмущению начальных значений при всех их вариациях в границах выполняется
.
Предварительно предполагалось выполнение (5). Требуется показать, что обе леммы и следствия сохраняются при ослаблении или исключении этого ограничения.
Лемма 3. Пусть для решений задачи (1) рассматриваются
, L из (2), ti из (4) и произвольные значения
, включая случай
, где
, при этом в некоторой окрестности
выполнено
. В этих условиях
(16)
Доказательство. В достаточно малой окрестности с учетом
и (2) выполняется (7). Функции
и
из (1) определены и непрерывны на
, поэтому
и
, разности под знаками пределов конечны на любом отрезке
, включающем
вместе с рассматриваемой окрестностью. Так что
, кроме того,
вследствие (2). Отсюда
, где
. В частности, для
,
,
,
. Лемма доказана.
Лемма 4. В условиях леммы 3
(17)
Для доказательства достаточно заметить, что относительно отношения под знаком предела в (17) с точностью до обозначений сохраняются рассуждения и оценки доказательства леммы 3. Поэтому
,
,
, отсюда следует (17).
Условия леммы 3 исключают неустранимые особенности в и
.
Следствие 3. Соотношение (10) будет сохраняться на произвольном отрезке , если условия леммы 1 расширить до условий леммы 3 и, кроме того, допускать любое конечное число точек, в которых на этом отрезке нарушается ограничение (5).
Доказательство. С учетом (16), (17) при однократном нарушении (5) ни одно из соотношений доказательства леммы 1 в расширенных условиях не изменится, поэтому (10) будет сохраняться. Если на количество точек
, где
конечно, то они взаимно отделены. В достаточно малой окрестности каждой из них нарушение ограничения (5) однократно, следовательно, (10) сохраняется. Следствие доказано.
Из следствия 3, соотношения (10) и леммы 2 вытекает
Теорема 1. В исходных предположениях, с тем изменением, что на произвольном отрезке допускается любое конечное число нарушений ограничения (5), для устойчивости, а также для асимптотической устойчивости решения задачи (1) сохраняются необходимые и достаточные условия леммы 2, включая соотношения (11), (12).
Согласно (13) имеет место
Следствие 4. Теорема 1 распространяется на условия устойчивости следствия 2, включая соотношения (14), (15) при .
Условия устойчивости можно дополнительно расширить. Пусть решение системы (1) устойчиво. Пусть произвольно зафиксировано k, и рассматривается
. Пусть сохраняются все исходные условия, с тем исключением, что множество точек
, в которых нарушается (5),
, является произвольным. Вне данного множества
,
, для разности
согласно следствию 1 сохраняется (10). В каждой такой точке необходимо выполняется (11) и (14), при этом в
-окрестности
константа
не меняется. В точках
, где
, функция
равна нулю. По непрерывности этой функции каждое
можно заключить в столь малый интервал
, что
. В результате все точки отрезка
покрываются системой интервалов
, на каждом из которых выполнено (14). По лемме Бореля из системы
можно выбрать конечную подсистему, покрывающую
. На каждом интервале конечного покрытия, следовательно, во всех точках
, для выбранных k и t выполнено (14). Рассуждения сохраняются
с одной и той же
. Ввиду произвольности t отсюда следует, что в случае устойчивости (14) с необходимостью выполнено
,
. Такое обоснование необходимости опирается на (10), что априори предполагает непрерывную дифференцируемость
. В то же время выполнение (14) является достаточным условием устойчивости без этого предположения. Если
,
,
,
, то
выполнено
, лишь только
. Отсюда
, где
, влечет
, и (14) – достаточное условие устойчивости без отмеченного ограничения. Из изложенного вытекает
Теорема 2. Пусть сохраняются все исходные условия, включая непрерывную дифференцируемость из (1), со следующим изменением: на произвольном отрезке
,
допускается любое множество точек
, в которых
. Тогда для устойчивости решения задачи (1) необходимо, чтобы существовало такое
, что
, при ограничении
,
верно неравенство
. (18)
Эти же условия являются достаточными условиями устойчивости без предположения о непрерывной дифференцируемости из (1). Для асимптотической устойчивости необходимо и достаточно, чтобы в соответственных условиях устойчивости существовало
, такое, что
влечет
. (19)
Остается доказать утверждение относительно (19). Поскольку , то (19) выполняется тогда и только тогда, когда
, поэтому утверждение равносильно определению асимптотической устойчивости. Теорема доказана.
Следствие 5. Пусть система (1) имеет точку покоя. Теорема 2 дает необходимые и достаточные условия ее устойчивости при , (18) переходит в соотношение
. (20)
Необходимые и достаточные условия асимптотической устойчивости точки покоя получаются из этой теоремы при переходе (19) в соотношение
. (21)
С учетом (2) имеет место
Следствие 6. В условиях теоремы 2 для устойчивости решения задачи (1) необходимо, чтобы существовало , такое, что
, при ограничении
выполняется неравенство
. (22)
Для асимптотической устойчивости необходимо, чтобы решение было устойчиво и существовало , такое, что
влечет
. (23)
В самом деле, согласно (2) ,
, и в случае устойчивости, с учетом (18)
,
. Отсюда
,
. Аналогично доказывается (23).
Замечание 1. В условиях теоремы 2 в случае устойчивости соотношение (22) влечет равномерную на полуоси оценку правой части (1): , и
выполнено
,
, лишь только
, где
. Таким образом, необходимым условием устойчивости является соотношение
. (24)
Замечание 2. Для асимптотической устойчивости необходимо, чтобы решение было устойчиво и существовало , такое, что
влечет
. (25)
Утверждение следует из того, что тогда и только тогда, когда
. Необходимость (25) следует также из условия Липшица.
Теорема 2, следствия и замечания применимы к нелинейным, автономным и линейным системам.
Условия устойчивости в случае автономной системы. Пусть для задачи Коши
(26)
где ,
,
, сохраняются все условия теоремы 2. В семейство решений (26) вместе с
всегда входит сдвиг
при любом
. Если в (18) положить
, то применительно к (26), в силу единственности решения,
. При таких начальных значениях в случае устойчивости условие
теоремы 2, или,
, влечет
. (27)
Таким образом, (27) – необходимое условие устойчивости решения задачи (26).
Поскольку выполнены условия теоремы о среднем, то
,
,
,
.
Из (27)
. (28)
В (28) меняется в зависимости от t, но неравенство переносится на
. Если предположить, что это не так и в условиях устойчивости
найдутся
,
, такие, что
, то с учетом
получится противоречие. В самом деле,
. В силу непрерывности
последнее неравенство сохраняется на некотором интервале, включающем
:
,
. Вследствие устойчивости на этом интервале выполняется (27). В частности,
, и
, или,
,
. Таким образом,
, и вместе с тем
, что невозможно. Следовательно, предположение неверно, и
, или
,
. Отсюда и из (24)
. В результате имеет место
Лемма 5. Если в условиях теоремы 2 решение задачи (26) устойчиво, то существует
,
, где
из (24),
из (27), такое, что
. (29)
В рассматриваемых ограничениях соотношение (29) – необходимое условие устойчивости решения автономной системы. Ниже определяются достаточные условия.
Лемма 6. Пусть выполнены условия теоремы 2 применительно к задаче (26). Тогда для устойчивости решения этой задачи необходимо и достаточно, чтобы , такие, что
, имеет место соотношение
. (30)
Условие (30) является достаточным для устойчивости решения (26), если выполнено (29).
Доказательство. Пусть . Можно считать, что
, иначе все сдвиги совпадали бы с невозмущенным решением, что для автономной системы невозможно [12]. Необходимость доказывается следующим образом. Если предположить, что решение устойчиво, но (30) не выполняется, то, как бы велико ни было c0 и как бы мало ни было
, найдутся
, Δ,
, и
, такие, что
. Тогда
можно выбрать c0 так, что
, и
. Но в силу устойчивости при достаточно малом
,
всегда выполнено
, что противоречит предположению. Достаточность получается с применением теоремы о среднем:
,
. Поэтому
,
и
.
Отсюда
,
тем более,
,
или,
.
Отсюда и из (29) . Таким образом, (30) влечет
,
, и, следовательно,
,
,
. В результате
выполнено
, лишь только
. Лемма доказана.
Для асимптотической устойчивости решения задачи (26) необходимо и достаточно соответственное выполнение условий леммы 6 и существование -окрестности
,
, в начальных условиях из которой выполняется соотношение
. (31)
Достаточность получается подстановкой (31) в правую часть (30) с учетом условий леммы, необходимость – из определения асимптотической устойчивости применительно к возмущению . Из изложенного вытекает
Теорема 3. Пусть выполнены все условия теоремы 2 применительно к задаче (26). Тогда для устойчивости решения этой задачи необходимо выполнение соотношений (27), (29) и (30). При этом если (29) выполнено, то выполнение (30) является достаточным условием устойчивости. Необходимым и достаточным условием асимптотической устойчивости решения задачи (26) является выполнение соответственных условий устойчивости при одновременном выполнении с ними соотношения (31).
Следствие 7. Из (27) следует известный факт равномерной непрерывности устойчивого решения автономной системы:
,
,
, где
из (27).
Из замечаний 1, 2 вытекает
Следствие 8. В условиях теоремы 3 для устойчивости решения задачи (26) необходимо, чтобы существовало , такое, что
, верно неравенство
. (32)
Для асимптотической устойчивости необходимо, чтобы решение было устойчиво и существовало , такое, что
выполнялось соотношение
. (33)
Замечание 3. В случае устойчивости решения задачи (26) соотношение (32) влечет равномерную непрерывность правой части на полуоси. Именно при
, выполнено
,
, лишь только
. Таким образом,
. (34)
При данных ограничениях (34) – необходимое условие устойчивости решения (26).
Условия устойчивости точки покоя автономной системы. Пусть система (26) преобразована к виду, при котором анализ устойчивости ее решения сводится к анализу устойчивости нулевого решения (точки покоя):
(35)
где ,
,
, не умаляя общности, можно считать
[1]. В дальнейшем формально предполагается, что
, если
, и, кроме того, соотношение
возможно только, если
. Применительно к анализу устойчивости нулевого решения задачи (35) вычитаемые в делимом (27), в соотношениях (30) – (34) становятся нулями. Ниже возмущение нулевого решения не будет отмечаться волной.
Следствие 9. Условия и утверждения теоремы 3 с точностью до обозначения повторяются для решения задачи (35). В частности, это относится к ограниченности в условиях устойчивости правой части
из (35). В случае устойчивости точки покоя задачи (35) соотношение (27) примет вид
,
,
и (35) перейдет в соотношение
,
.
Следствие 10. Применительно к (35) сохраняется следствие 8, при этом (32), (33) переходят в соотношения ,
, и
. Сохраняется также замечание 3, при этом (34) перейдет в соотношение
.
Согласно этим следствиям в случае устойчивости точки покоя системы (35) решение с начальными значениями из некоторой окрестности
и правая часть
ограничены. При некоторых условиях ограниченность V(t) может определяться знаком компонентов
, а
– знаком компонентов
. Покомпонентное сопоставление знаков производных приводит к достаточным условиям устойчивости, асимптотической устойчивости и неустойчивости, излагаемым ниже. По-прежнему относительно (35) сохраняются все условия теоремы 2.
Чтобы исключить неверный смысл последующих утверждений, дополнительно предполагается, что всюду, где используется знак нестрогого неравенства относительно ,
,
, решения системы (35) исключают каждый из случаев, когда
или
, или
.
Очевидны следующие утверждения.
1. Если , такое, что
,
верны неравенства
,
,
, то
ограничена, а точка покоя системы (35) устойчива.
В данных условиях функция не убывает, поэтому по модулю не возрастает, и
,
. Функция
не возрастает, поэтому
. Отсюда
, при выборе
, выполнено
,
.
2. Если , такое, что
,
верны неравенства
,
,
, то
ограничена, а точка покоя системы (35) устойчива.
Утверждение оправдывается по аналогии с пунктом 1.
3. Если для некоторых выполнены условия пункта 1, а для всех остальных
– условия пункта 2, то
ограничена, а точка покоя системы (35) устойчива.
При выполнено:
,
. Если
, то
, и
, при выборе
, верно
.
Ограниченность в общем случае – следствие устойчивости согласно лемме 5, ниже это свойство отдельно не оговаривается. Из пунктов 1–3 вытекает
Предложение 1. Если , такое, что
,
выполняется любая из пар неравенств
,
или
,
, то точка покоя системы (35) устойчива.
Теорема 4. Пусть выполнены условия предложения 1. Если , такое, что
,
выполняется пара неравенств
,
или
,
, и при этом
, то точка покоя системы (35) асимптотически устойчива.
Доказательство. Согласно предложению 1 точка покоя устойчива. Пусть при условии произвольно зафиксировано
, выбрано
, при котором выполняется
,
, и пусть
. Предположение, что
, окажется в противоречии с устойчивостью. В этом предположении
, такое, что, начиная с некоторого
, будет верно неравенство:
. Тогда
, где
, и
, если
. Отсюда
при
, что означает неустойчивость нулевого решения вопреки условию. Предположение неверно, и в рассматриваемых условиях выполнено
. Поскольку
, и
тогда и только тогда, когда
, с учетом произвольности выбора
, соотношение
возможно только, если
. Поэтому точка покоя асимптотически устойчива. Аналогично рассматривается случай
,
.
Предложение 2. Пусть точка покоя системы (35) устойчива. Если , такое, что
,
выполнены неравенства
,
, то точка покоя асимптотически устойчива.
Доказательство. Пусть зафиксировано и выбрано
. Функция
убывает и ограничена снизу. Поэтому
. Если предположить, что
, возникнет противоречие. В самом деле,
,
,
,
, такие, что
. Отсюда
, где
. Следовательно,
вопреки устойчивости точки покоя. Предположение неверно, и
. С учетом произвольности выбора
и
соотношение
возможно только, если
. Отсюда следует асимптотическая устойчивость точки покоя.
Предложение 3. Пусть точка покоя системы (35) устойчива. Если ,
выполнено
,
, то точка покоя асимптотически устойчива.
Доказательство аналогично предыдущему.
Если для некоторых выполнены условия предложения 2, а для всех остальных
– условия предложения 3, то точка покоя системы (35) асимптотически устойчива. Отсюда вытекает
Теорема 5. Если точка покоя системы (35) устойчива и , такое, что
,
выполнено либо
,
, либо
,
, то точка покоя асимптотически устойчива.
Замечание 4. Согласно замечанию 2 соотношение (25) – необходимое условие асимптотической устойчивости. Для системы (35) оно примет вид: ,
выполняется
. Для доказательства теоремы 4 и для доказательства предложения 2 достаточно выполнения самого этого соотношения и факта устойчивости, без использования условий, в силу которых это соотношение имеет место, в частности неравенств из условий предложений 1, 2.
На основании замечания 4 имеет место
Теорема 6. Если точка покоя системы (35) устойчива, то для ее асимптотической устойчивости необходимо и достаточно, чтобы в некоторой окрестности выполнялось
.
Непосредственно из теорем 4, 5 вытекает
Теорема 7. Если , такое, что
,
выполняется любая из троек неравенств
,
,
или
,
,
, то точка покоя системы (35) асимптотически устойчива.
В ограничениях начального предположения все знаки строгих неравенств в условиях теорем 4, 5, 7 можно заменить знаками соответственных нестрогих неравенств с сохранением утверждений. На пункты 4, 5 аналогичное утверждение не распространяется.
4. Если имеет место устойчивость, теорема 6 дает необходимые и достаточные условия асимптотической устойчивости. Кроме того, в условиях постоянства знака компонентов решения достаточные условия теорем 4, 5, 7 становятся необходимыми.
Предложение 4. Пусть , такое, что
,
задача (35) имеет решения со свойством постоянства знаков компонентов,
, то есть либо
, либо
. Пусть
из рассматриваемой окрестности начальных данных
. Тогда для асимптотической устойчивости точки покоя необходимо и достаточно выполнение одной из пар неравенств
,
или
,
.
Доказательство. Достаточность условий доказывается аналогично тому, как доказывалась достаточность условий теоремы 4. Пусть теперь точка покоя системы (35) асимптотически устойчива, тогда решение V(t) ограничено. Если предположить, что при некотором
условия теоремы 6 нарушены, например в виде
,
, то возникнет противоречие. Функция vk положительна, не убывает и ограничена, следовательно,
. Но значение
противоречит асимптотической устойчивости, в силу которой
. Предположение неверно, и в рассматриваемых условиях необходимо
,
. Если нарушение условий возникло в форме
,
, то в силу устойчивости функция vk(t) отрицательна, ограничена на полуоси, при этом не возрастает. Поэтому
. Такое значение
противоречит асимптотической устойчивости. Отсюда необходимо
,
. Случаи нарушения условий
,
рассматриваются аналогично.
С использованием ограниченности согласно (29) аналогично доказывается
Предложение 5. Пусть точка покоя системы (35) устойчива. Тогда для ее асимптотической устойчивости необходимо и достаточно выполнение одной из пар неравенств ,
или
,
.
5. Следующие условия достаточны для неустойчивости.
Предложение 6. Если для задачи (35) :
,
,
, такие, что неравенства
,
верны
, при этом либо
, либо
, то точка покоя неустойчива. Аналогично, если
,
.
Доказательство. Пусть ,
, где k из условий предложения. Если функция vk(t) ограничена сверху, то
, при этом
с учетом возрастания vk(t). Если в этом случае предположить, что точка покоя устойчива и
, то согласно предложению 5 и теореме 6 необходимо
. Тогда
, что противоречит
. Если же предположить, что точка покоя устойчива и
, то в силу леммы 5 функция uk(t) ограничена, не убывает, следовательно, имеет предел
. Отсюда, как и раньше,
, что противоречит устойчивости. Остается предположить, что vk(t) не ограничена сверху, но это также означает, что точка покоя неустойчива. Случай
,
рассматривается аналогично.
Нетрудно видеть, что к тому же результату приводит выполнение соотношений ,
или
,
при сохранении условий относительно
.
С помощью аналогичных рассуждений доказывается
Теорема 8. Точка покоя системы (35) неустойчива, если :
,
,
, такие, что
верна хоть одна из пар следующих неравенств при сохранении условий предложения 6 относительно
:
,
;
,
;
,
;
,
. В тех же условиях решение неустойчиво, если верна хоть одна из пар следующих неравенств:
,
;
,
;
,
;
,
.
Элементарные примеры и программная реализация. Нулевое решение уравнения неустойчиво:
, если
, то
,
; если
, то
,
(теорема 8). Нулевое решение уравнения
асимптотически устойчиво:
, если
, то
,
; если
, то
,
(теорема 7). Уравнение
неустойчиво:
(теорема 8), уравнение
асимптотически устойчиво:
(теорема 7). Для системы
знаки
и
,
и
покомпонентно противоположны:
,
. Однако нельзя сделать вывод об асимптотической устойчивости – не выполнены исходные условия, которые предполагают теоремы 4, 6, 8, в частности не выполнено (2). В общем случае аналитически идентифицировать асимптотическое поведение знака решения, его первой и второй производной далеко не всегда возможно. Однако это не составляет затруднения при применении компьютера: достаточно по ходу решения выводить приближения vk(t), uk(t),
. Например, для уравнения
следующая программа (Delphi) даст искомую информацию:
Program sgnnorm;{$APPTYPE CONSOLE} uses SysUtils;
const h = 0.000001;eps =0.001; tt=1000000; var t,v: extended; k: longint;
function u (t,v:extended):extended; begin u:=exp(v)-1; end;
function u1 (t,v:extended):extended; begin u1:=exp(v)*(exp(v)-1); end;
begin k:=0; v:= eps; t:=0; while t <= 6 do begin v:=v+h*u (t,v); t:=t+h; k:=k+1; if k = tt then begin writeln ('t=',t:4,' ');
writeln; writeln (' ','znaki':32,' '); writeln; writeln (' ',v,' ',u (t,v),' ',u1 (t,v)); writeln; writeln; k:=0 end; end;
readln end.
В качестве возмущения нулевого решения с отсчетом tt шагов выводится приближение по методу Эйлера решения v (v), , правой части
(u(t,v)), ее производной
(u1(t,v)). Их знаки положительны на всем интервале приближения, что означает признак неустойчивости. Интервал взят малым, иначе возникнет переполнение, что дает дополнительный признак: заведомо нарушается необходимое условие устойчивости (20). Если по данной программе приближенно решать уравнение
, интервал можно произвольно продлить, знаки решения, первой и второй производной будут чередоваться,
, все эти величины будут убывать к нулю, означая асимптотическую устойчивость. При варьируемых начальных значениях чередование сохраняется. В продолжение процесса обнулятся значащие цифры выводимых значений в расширенном формате. Согласно (21) и теореме 6 стремление решения, а также производной к нулю дает дополнительный признак асимптотической устойчивости. По аналогичной программе можно выводить отношение решения к начальному значению (v/eps) для проверки необходимых и достаточных условий (20), (21).
В общем случае для исследования системы (35) программа с очевидностью модифицируется. В результате она выполнит аналогичные действия для каждого уравнения в отдельности, при этом производная от априори аналитически определяется по формуле
, или,
[13]. Значения
и
задаются подпрограммами-функциями, их приближения выводятся одновременно с компонентами приближенного решения
. Использование компонента
вектор-функции
из (35) сходно с применением производной в силу системы [1, 3]. Принципиальное отличие функции Ляпунова, использующей известное построение, в том, что в общем случае она не является вектор-функцией, ее значение принадлежит R как функции n переменных. Известное видоизменение – вектор-функция Ляпунова [3] – имеет число компонентов меньше n и использует вспомогательную систему. Рассматриваемое применение
ограничено условиями существования и непрерывности
,
, что не необходимо для применения функции Ляпунова. Ее построение не детерминировано формулой и оставляет возможность оценки устойчивости с помощью выбора переменных, причем в более общем случае. Однако для численного моделирования устойчивости существенна инвариантность задания
по формуле производной сложной функции.
Признаки устойчивости линейных систем. Как частный случай автономной системы (35) рассматривается однородная система
(36)
с матрицей вещественных коэффициентов n×n. К (36) применимы теоремы 4–8 и данные ранее утверждения относительно автономных и нелинейных систем. Выражение производных наглядно просто: . Одним из условий асимптотической устойчивости системы является противоположность знаков k-х компонентов AY и
при каждом
. Одно из условий неустойчивости – наличие при некотором k противоположных знаков k-х компонентов AY и
. Пусть, например, рассматривается система
. В этом случае
,
,
,
,
.
Из записи метода Эйлера для приближенного решения этой системы, , видно, что существует решение с положительными компонентами. При
знаки компонентов Y, AY и
совпадают, решение неустойчиво (теорема 8). Такой вывод можно повторить для любой положительной матрицы A, что согласуется с теоремой Перрона о наибольшем по модулю собственном значении положительной матрицы: оно вещественно и положительно [14]. Пусть рассматривается система
. Здесь
,
,
,
,
.
Второй компонент меняет знак при переходе от Y к AY и от AY к . Если, начиная с некоторого
, выполнено
, то и первый компонент меняет знак аналогично. В этом предположении знаки компонентов Y, AY и
попарно противоположны
, система асимптотически устойчива согласно теореме 7. Легко убедиться непосредственной подстановкой, что в семейство решений входит
, это решение соответствует предположению и иллюстрирует асимптотическую устойчивость. Приближенное решение системы проявило бы те же и другие свойства асимптотического поведения ненулевого решения, в частности признак асимптотической устойчивости:
, если
(теорема 6).
Ниже оговаривается, что для оценки устойчивости системы (36) достаточно выполнения обсуждаемых признаков для одного произвольного решения со всеми ненулевыми компонентами начальных значений.
Невыполнение достаточных условий теорем и предложений не означает неустойчивость. Так, система не отвечает условиям этих утверждений, но устойчива: собственные значения матрицы
мнимые и не кратные [14].
Как и в общем случае, для системы (36) нельзя предложить универсальный алгоритм аналитической идентификации знаков компонентов решения и производных. Однако компьютерная идентификация знаков выполняется относительно без затруднений. В программу вводится матрица A, программно вычисляется A2. Для компонентов AY и задаются подпрограммы-функции, соответственные каждой из строк матриц, значения функций выводятся по шагам приближенного решения одновременно с компонентами решения Y в случае возмущенных начальных значений. Такой процесс обладает естественным параллелизмом, в последовательной реализации A2 можно не вычислять: выводятся компоненты Y, AY, по ассоциативности –
.
Для дальнейшего исследования потребуется рассмотреть структуру общего решения системы (36), которое в форме Коши имеет вид [15]
, (37)
где – начальные значения, которые считаются произвольными,
– компоненты нормированной фундаментальной системы,
. Нормированность означает [15], что любое частное решение с единственностью получается из (37) при фиксировании начальных значений именно этого решения:
(включая
). Отсюда любое частное решение системы (36) в выражении из (37) в качестве элементов линейной комбинации содержит все n вектор-функций фундаментальной системы, каждую с ненулевым коэффициентом, при условии, что начальные значения решения имеют все ненулевые компоненты:
. Поэтому одно такое решение сохраняет свойства всех решений фундаментальной системы в смысле асимптотического поведения, определяющего характер устойчивости. Вектор-функция фундаментальной системы представляет собой решение с компонентами
,
, где каждому собственному числу λi кратности
матрицы A соответствуют полиномы
степени не выше
, в совокупности имеющие
произвольных коэффициентов, через которые выражаются остальные [15]. При
полином имеет нулевую степень и выражается в виде одного свободного коэффициента. В случае комплексного λi отделяются его вещественная и мнимая части, по отдельности они порождают вещественные компоненты
. Вхождение каждого
с ненулевым коэффициентом в (37) сохраняет асимптотический характер роста в зависимости от знака действительной части λi либо от кратности мнимой части, если действительная часть равна нулю. В любом случае асимптотическое поведение одного частного решения, начальные значения которого составляют только ненулевые компоненты, эквивалентно в смысле устойчивости поведению всей фундаментальной системы, соответственной n собственным значениям матрицы системы. Это коррелирует с известным свойством системы (36): для ее устойчивости необходима и достаточна ограниченность одновременно всех решений, для асимптотической устойчивости необходимо и достаточно стремление к нулю одновременно всех решений при
[1]. Однако по асимптотике какого именно из решений можно точно судить об устойчивости всей системы, ясно только на основании (37): по поведению любого частного решения, все компоненты начальных значений которого не равны нулю. На этой основе для анализа устойчивости системы (36) не требуется указывать окрестность возмущения начальных значений, как в следствии 5 в общем случае. Достаточно формулировать условия этого следствия для произвольного решения с ненулевыми компонентами начальных значений. Применительно к (36) не требуется проверять выполнение условий предложений и теорем в окрестности нулевых начальных значений: достаточно их проверить для произвольно выбранных ненулевых компонентов начального вектора. Выполнение условий влечет в асимптотике ограниченность, стремление к нулю или неограниченность решения, характерные для всей фундаментальной системы.
Пусть в (36) ,
. Компоненты
,
,
находятся из соотношений
,
,
,
. (38)
Согласно изложенному из теорем 4–8 применительно к (36) вытекает
Теорема 9. Пусть решение системы (36) выбрано произвольно при условии
и компоненты
,
,
определяются из (38). Если
выполняется любая из пар неравенств
,
или
,
, то система (36) устойчива. Если
выполняется пара неравенств
,
или
,
, и при этом
, то рассматриваемая система асимптотически устойчива. Если эта система устойчива и
выполнено либо
,
, либо
,
, то она асимптотически устойчива. Если система (36) устойчива, то для ее асимптотической устойчивости необходимо и достаточно, чтобы выполнялось
. Если
выполняется любая из троек неравенств
,
,
или
,
,
, то система (36) асимптотически устойчива.
Теорема 10. Cистема (36) неустойчива, если ,
, такие, что
верна хоть одна из пар неравенств:
,
;
,
;
,
;
,
; и при этом либо
, либо
. В тех же условиях решение неустойчиво, если верна любая из следующих пар неравенств:
,
;
,
;
,
;
,
.
Имеет место
Предложение 7. Пусть решение системы (36) выбрано при условии
. Тогда для устойчивости системы необходимо и достаточно, чтобы хоть одна производная решения произвольного порядка была ограничена на полуоси:
,
,
. Для асимптотической устойчивости необходимо и достаточно, чтобы для такой производной выполнялось соотношение
,
.
Доказательство. Очевидно,
,
,
и
.
Поэтому производная ограничена в том и только в том случае, если ограничено решение
. В рассматриваемых условиях ограниченность одного решения определяет ограниченность одновременно всех решений, что необходимо и достаточно для устойчивости системы (36) [1]. Для ее асимптотической устойчивости необходимо и достаточно стремление к нулю одновременно всех решений [1]. При рассматриваемом выборе решения это определяется тем, что
при
. Последнее соотношение выполняется тогда и только тогда, когда
,
. Предложение доказано.
Следствие 11. Предложение 7 сохраняется, если в качестве условия устойчивости требовать ограниченности на полуоси производных одновременно всех порядков, в качестве условия асимптотической устойчивости – их одновременного стремления к нулю при .
Замечание 5. В условиях предложения 7 в случае неасимптотической устойчивости системы (36) необходимо найдется , такое, что компонент производной
любого порядка циклически меняет знак с ростом t на полуоси.
Если предположить обратное, например, что
,
:
, то
не убывает. По следствию 11 решение и производные всех порядков ограничены в силу устойчивости. Отсюда
,
. Аналогично в случае
. В любом случае
, такое, что знак
, а вместе с ним и знак
постоянны
. Для
можно повторить рассуждение, в результате
, такое, что знак
постоянен
. По индукции
, такое, что знак
постоянен
, и
, такое, что знак
постоянен
. Не умаляя общности, можно считать, что
. Тогда
компоненты
и
имеют постоянные знаки. Если постоянные знаки
,
совпадают (знак
постоянен), то по теореме 10 решение неустойчиво. Если эти знаки противоположны (
), то по теореме 9 решение асимптотически устойчиво. В обоих случаях это противоречит условию, по которому решение устойчиво, но не асимптотически. Следовательно, предположение неверно. С учетом того, что рассуждения распространялись на случай
при
, утверждение замечания доказано.
С учетом окрестности нулевых начальных значений сходное утверждение доказывается относительно решения системы (35).
Замечание 6. В случае неасимптотической устойчивости решения системы (35)
, такое, что
, при котором компоненты
,
циклически меняют знак с ростом t на полуоси.
Для компьютерного анализа устойчивости системы (36) целесообразно применять еще один способ численного моделирования. Он основан на том, что система (36) устойчива тогда и только тогда, когда ,
,
, для асимптотической устойчивости необходимо и достаточно, чтобы
, где h из (4) [10, 11]. Вычисление степени
заменяется вычислением
:
(39)
При достаточно малом h реализация (39) сводится к умножению матрицы на себя до момента достижения заданной границы изменения t. Следующая программа совмещает анализ устойчивости системы (36) на основе теорем 8, 9 и соотношений (39).
Program sgnlinconst;{$APPTYPE CONSOLE} uses SysUtils;
const n=3; h=1.1e-14; x0=0; xn=1500; h1=1.1e-4;TT=150000;
type matr=array[1..n,1..n] of extended; vect=array[1..n] of extended;
const A: matr= ((-1, -0.09, 0.077),
(0.087, -0.9, 0.005),
(-0.034, 0.034, -0.2));
var a1,c: matr; y1,y2,y3,z1,z2,z3: extended; s0,s1,x: extended; I,j,l,k,k0,kk: integer;
function ay1 (const a:matr; var y1,y2,y3: extended): extended;
begin ay1:=a[1,1]*y1 +a[1,2]*y2+a[1,3]*y3; end;
function ay2 (const a:matr; var y1,y2,y3: extended): extended;
begin ay2:=a[2,1]*y1 +a[2,2]*y2+a[2,3]*y3; end;
function ay3 (const a:matr; var y1,y2,y3: extended): extended;
begin ay3:=a[3,1]*y1 +a[3,2]*y2+a[3,3]*y3; end;
procedure ummatr (var a1,c: matr);
var s1: extended; I,j,l : integer;
begin for I := 1 to n do for j := 1 to n do
begin s1:=0; for l:= 1 to n do s1:= s1+a1[I,l]*a1[l,j]; c[I,j]:= s1 end; end;
begin kk:=0; for i:=1 to n do for j:=1 to n do begin a1[I,j]:=a[I,j]*h; if i=j then a1[I,j]:=a1[I,j]+1 end;
x:=x0; y1:=2;y2:=2;y3:=2; while x<=xn do begin kk:=kk+1;
y1:=y1+ay1(a,y1,y2,y3)*h1; y2:=y2+ay2(a,y1,y2,y3)*h1; y3:=y3+ay3(a,y1,y2,y3)*h1;
if kk=TT then begin z1:=ay1(a,y1,y2,y3);z2:=ay2(a,y1,y2,y3);z3:=ay3(a,y1,y2,y3);
write (' ':4,y1:4,' ':4,ay1(a,y1,y2,y3):4,' ':4,ay1(a,z1,z2,z3):4);writeln;
write (' ':4,y2:4,' ':4,ay2(a,y1,y2,y3):4,' ':4,ay2(a,z1,z2,z3):4);writeln;
write (' ':4,y3:4,' ':4,ay3(a,y1,y2,y3):4,' ':4,ay3(a,z1,z2,z3):4);writeln;
writeln; writeln; kk:=0; end; x:=x+h1; end; writeln; writeln;
k:=0; x:=x0; while abs(x) <= 1e9{3} do begin ummatr (a1,c); k:=k+1;x:=h*exp((k+1)*ln(2));
for i:=1 to n do for j:=1 to n do a1[I,J]:=c[I,j]; s0:=0; for i:=1 to n do for j:=1 to n do
s0:=s0+sqr(a1[I,j]); s0:= sqrt(s0); write (' ':2, s0:2,' ':8); end; writeln; writeln; writeln (' ':2, 'шаг =',h,' ');
writeln (' ':2, 'параметр k=',k,' ');writeln (' ':2, 'параметр x=',x:2,' '); readln end.
Результат выполнения программы:
Y
6.8E-0003 -1.3E-0003 2.8E-0004
1.3E-0003 -2.8E-0004 5.7E-0005
7.2E-0002 -1.4E-0002 3.0E-0003
2.4E-0004 -4.9E-0005 9.9E-0006
4.8E-0005 -9.8E-0006 2.0E-0006
2.5E-0003 -5.2E-0004 1.0E-0004
………………………………………………………
4.4E-0132 -8.9E-0133 1.8E-0133
8.8E-0133 -1.8E-0133 3.6E-0134
4.7E-0131 -9.4E-0132 1.9E-0132
1.7E + 0000 1.7E + 0000 1.7E + 0000 1.7E + 0000
1.7E + 0000 1.7E + 0000 1.7E + 0000 1.7E + 0000
……………………………………………………………......
1.7E + 0000 1.7E + 0000 1.6E + 0000 1.5E + 0000
1.3E + 0001 1.1E + 0000 8.0E-0001 5.4E-0001
……………………………………………………………......
1.9E-0009 3.7E-0018 1.4E-0035 1.9E-0070
3.5E-0140 1.2E-0279 1.4E-0558 2.0E-1116
……………………………………………………………......
0.0E + 0000 0.0E + 0000 0.0E + 0000 0.0E + 0000
0.0E + 0000 0.0E + 0000 0.0E + 0000 0.0E + 0000
Система (36) с матрицей A, 3×3, заданной в разделе описания констант, решается методом Эйлера с шагом на отрезке
. По ходу решения с начальными значениями
,
выводятся компоненты
,
,
. Компоненты убывают до
, по следствию 11 это признак асимптотической устойчивости. Знаки компонентов чередуются на всем отрезке решения:
,
,
,
. По теореме 8 это также признак асимптотической устойчивости. Соотношение (39) реализуется с шагом
путем умножения матрицы
на себя до момента
(закомментировано
на случай переполнения от роста нормы при неустойчивости). Норма убывает от
до
при количестве умножений матрицы на себя
. Согласно (39) это дополнительный признак асимптотической устойчивости. Признаки указывают правильный результат, поскольку матрица имеет диагональное преобладание отрицательных элементов.
Численный эксперимент. Используются эвклидовы нормы матрицы и вектора. Пусть в (36) в качестве матрицы коэффициентов поочередно рассматриваются
,
. (40)
Если в программе sgnlinconst ввести изменения, соответственные размерности и конкретным данным, выводить левую часть (20) в обозначении при
,
,
, то в случае A1 получится
: 1.02×100, …, 5.16×10-17, …, 3.72×10-29, …, 8.07×10-86, …, 1.28×10-104, …, 2.67×10-134. Согласно (21) это признак асимптотической устойчивости (соответствует диагональному преобладанию отрицательных элементов). Для этой же матрицы получится признак асимптотической устойчивости, соответственный теореме 9:
,
,
,
. Кроме того, компоненты
,
,
,
убывают от 6.1×10-11 до 2.6×10-142, по следствию 11 это также указывает на асимптотическую устойчивость. Реализация (39) даст для A1 еще один признак асимптотической устойчивости
: 2.0×100, 2.0×100, … , 1.9×100, 1.8×100,…, 2.3×10-36, 5.2×10-72,…, 0, 0. В случае A2 преобладают положительные диагональные элементы, рост нормы соответствует неустойчивости,
: 5.75×100, …, 1.01×1020, …, 3.68×10100, …, 1.85×10491, …, 4.30×10656, …, 8.62×10772, …, 4.18×10908. Другой признак неустойчивости:
;
;
;
. Компоненты растут до 8.8×10893 , и далее, до переполнения, нарушается необходимое условие устойчивости следствия 11. Наконец,
: 2.0×100, 2.0×100, … , 2.0×100, 2.1×100, 2.2×100, …, 1.1×10240, 1.3×10480, далее следует переполнение. Для обеих матриц (40) признаки указывают правильные результаты. Можно рассмотреть пограничные состояния устойчивости линейных уравнений, взяв в (36) матрицу Фробениуса с нижней строкой, соответственной характеристическому полиному с кратными (или напротив) мнимыми корнями. Пусть, например,
,
.
Последняя строка – коэффициенты характеристических полиномов, за вычетом единичных при старших степенях, взятые с обратным знаком и расположенные в порядке возрастания индексов. Характеристический полином A3 имеет кратные мнимые корни: или
, в этом случае система неустойчива [14]. Программа sgnlinconst для A3,
,
,
, даст
: 3.35×101, …, 3.37×1067, …, 3.71×10146, …, 1.76×10241, …, 1.32×101276 , 2.82×101290. Ввиду нарушения условия (20) это признак неустойчивости. То же получается в соответствии с теоремой 10:
,
,
,
. Компоненты
,
,
растут до 2.9×101285, нарушая необходимое условие устойчивости следствия 11. Еще один признак неустойчивости,
: 2.4×100, 2.4×100, …, 2.4×100, 2.5×100, …, 9.9×1002, …, 3.5×10344, 6.0×10688, далее следует переполнение, что нарушает (39). Для матрицы A4 характеристический полином имеет вид
или
, мнимые корни не кратны, действительные части остальных корней отрицательны, что соответствует неасимптотической устойчивости [14]. Программа с теми же входными параметрами даст
: 1.25×101 , …, 2.25×101 , …, 2.36×101 , …, 2.21×101 , …, 1.79×101 , …, 1.55×101, …, 1.95×101, …, 1.93×101. Это признак устойчивости согласно (20). Изменение компонентов
,
,
,
ограничено по модулю, например при
соответственно от 1.4×10-8, 4.9×10-9, 5.6×10-8 до 9.0×10-9, 2.2×10-8, 3.6×10-8, аналогично при
, что указывает на устойчивость по следствию 11. В тройке
,
,
знаки циклически меняются с ростом t. Точнее, при каждом
всегда имеется пара противоположных знаков соседних компонентов тройки, для смежной пары один из знаков повторяется. При этом комбинация знаков меняется с ростом t и циклически повторяется, что по замечанию 5 иллюстрирует именно неасимптотическую устойчивость системы. На том же отрезке
: 2.4×100, …, 2.4×100, …, 7.3×100, …, 1.0×101, …, 1.2×101, …, 9.4×100, …, 6.7×100, …, 1.1×101, …, 7.0×100. Согласно (39) признак указывает на устойчивость. В результате признаки соответствуют неустойчивости системы (36) с матрицей A3 и ее неасимптотической устойчивости в случае A4. Пусть теперь рассматривается система
. Известно, что ее нулевое решение устойчиво, все остальные неустойчивы. Программа, реализующая (20), аналогична приведенной выше sgnrost, отличается размерностью и функциями правой части. Признак устойчивости нулевого решения получается, в частности, при
,
,
, шаг метода Эйлера
. В прежнем обозначении выводится
: 1.4, 1.4, … , 1.4. Признак сохраняется при вариации начальных значений. То, что устойчивость нулевого решения не асимптотическая, иллюстрируется, в частности, в случае
: компоненты
и
циклически меняют знаки с ростом
в соответствии с замечанием 6, при этом
не выходит из отмеченных выше границ. Неустойчивость ненулевых решений на основе соотношения (18) проявляется, начиная с начальных значений
, при возмущении
. В обозначении
получится
: 1.41, 1.41, …, 1.42, …, 1.43, …, 1.447. Если
, то
: 1.4, 1.5, …, 1.6, …, 1.9, …, 2, …, 2.1, …, 3, …, 3.3. Если
, то
: 2, …, 4, …, 12, …, 89, …, 91, …, 180, …, 251, …, 300, …, 312. Признак сохраняется при вариации возмущений начальных значений. Для ненулевого решения признак в соответствии с теоремой 7 непосредственно не получится, поскольку теорема дана для нулевого решения системы (35). Однако если численно смоделировать замены переменных
, сводящие анализ устойчивости ненулевого решения
к анализу нулевого решения
, то требуемый признак проявится при
и
,
. На всем отрезке
знаки второго компонента решения и производной будут совпадать:
, при этом
монотонно растет. Согласно эксперименту при возмущениях любых ненулевых начальных значений
неизменно растет, указывая на неустойчивость ненулевых решений.
Инвариантность признаков (18)–(21) для нелинейных уравнений далее иллюстрирует пример с использованием аналитических преобразований и оценок. Из (10) . Лемма 2 сохранится при замене (11) на соотношение
,
и (12) – на соотношение
[11]. Согласно (16), (17) в условиях теоремы 1 при достаточно малом h под знаком логарифма исключены неустранимые особенности. С учетом
в тех же условиях лемма 2 сохранится при замене (11) на соотношение
и (12) – на соотношение
. На этой основе в рамках теоремы 1 ограничение (5) снимается с интегральных аналогов
как пределов интегральных сумм на
, элементами которых являются
из (6), (10). Пусть задача (1) записана в виде
(41)
где ,
,
, и пусть она имеет нулевое решение
,
,
, устойчивость которого требуется исследовать. В дальнейшем для задачи (41) предполагаются выполненными условия теоремы 1 относительно ограничения (5), и любое ненулевое решение (41) интерпретируется как возмущение решения
. При замене
на интеграл справедлива
Лемма 7 [11]. Лемма 2 сохраняется для из (41), если соотношения (11) заменить на
,
и (12) – на
,
.
В точках устранимых особенностей, , дроби
доопределяются:
, и ограничение (5) снимается с условий двух теорем из [11].
Теорема 11. Если :
, выполнено
,
, где
, то решение
устойчиво. В частности,
,
. Если
,
,
:
, где
, то решение
неустойчиво. В частности,
,
.
Теорема 12. Если для (41) выполнены условия устойчивости теоремы 1 и ,
,
:
,
, где
, то решение
асимптотически устойчиво. В частности,
,
,
,
.
Условиям теоремы 12 удовлетворяет система
(42)
где ,
;
,
;
,
– вещественные числа
,
. Для (42)
, где
, при этом
, или,
. Кроме того,
, где
,
, и
. Решение
асимптотически устойчиво. На основе (20), (21),
для значений n = 2,
,
,
,
модифицированная программа sgnrost,
,
,
, шаг метода Эйлера
, даст соответственный признак
: 5.33×10-1, …, 6.17×10-3, …, 9.66×10-10, …, 9.24×10-14, …, 4.92×10-16, …, 2.56×10-16.
Кроме того, ,
,
,
, для правой части в данном приближении выполнено
(на конце отрезка
).
Если в (42) множитель заменить на –1, то полученная система снова будет удовлетворять условиям теоремы 12, можно положить
, откуда следует асимптотическая устойчивость решения
. Программа даст
: 3.99×10-5, …, 7.65×10-12, …, 8.26×10-21, …, 2.31×10-38, …, 1.91×10-54, …, 1.42×10-63, …, 1.45×10-66.
Как и в предыдущем случае, ,
,
,
, на конце отрезка
.
Если заменить на +1 то
, где
,
нулевое решение окажется неустойчивым. Если предположить обратное – что устойчиво, то
,
,
: при условии
выполнено
, так что
,
. С учетом
при достаточно малом ε верно
. Отсюда
.
В силу , начиная с некоторого
, выполняется
. С точностью до устранимых особенностей
, где
. Поскольку
, то
,
, и
. Функция
ограничена
, отсюда
. По теореме 10 решение
неустойчиво вопреки предположению. Предположение неверно, и решение неустойчиво. Программа даст признак
: 1.75×103, …, 6.05×1022, …, 8.43×1040, …, 1.20×1050, …, 1.17×1066, …, 4.96×1074, …, 8.42×1077 .
Кроме того, ,
(
,
),
. Все признаки сохраняются в единичной окрестности нулевых начальных значений.
Система (42) не является автономной, поэтому для оценки устойчивости нельзя сослаться на теоремы 4–8 и предложения 1–6. Тем не менее в рассмотренных трех видоизменениях этой системы экспериментальное поведение знаков компонентов решения и производной, стремление производной к нулю на полуоси соответствуют предложению 4, теоремам 6 и 8. Можно допустить, что при некоторых ограничениях предложения 1–6 теоремы 4–8 переносятся на общий случай системы (41) и на ее частный случай (42). В самом деле, при доказательстве предложений и теорем в случае автономной системы нигде не использовалась специфика системы (35), а применялись только предельные свойства монотонной функции одной переменной в зависимости от знака ее производной. Эти свойства были нужны, чтобы получить стремление к нулю компонентов правой части:
. Как только эти соотношения выполнялись, из них автоматически следовало стремление к нулю решения:
. Отсюда рассматриваемые рассуждения переносятся на случай конкретной разновидности системы (41), если для нее выполняется
и то ограничение, что из
следует
. Нетрудно видеть, что для системы (42)
, а соотношения знаков
и
при достаточно большом t1 тривиально выводятся аналитически. Из этих соотношений непосредственно следует устойчивость точки покоя первых двух рассмотренных разновидностей системы (42) и ее неустойчивость в случае третьей разновидности. Однако доказать асимптотическую устойчивость точки покоя двух первых разновидностей системы с помощью аналитического вывода соотношения
уже затруднительно. Более того, принципиальную трудность представляет показать, что из этого соотношения следует
: правая часть (42) может стремиться к нулю за счет множителей – степенных функций с отрицательными показателями. Однако можно отметить, что экспериментальное выполнение
согласуется с оценками на основе теоремы 12 – из оценок следует
, если
, что возможно только, если
при
. Однако анализ на основе теоремы 12 уже влечет вывод об асимптотической устойчивости и не требует иных доказательств. Остается заключить, что представленные методы не перекрывают друг друга, они целесообразны в соответствии с классами задач и возможностями численного моделирования устойчивости.
Относительно ограничений параметров в ходе численного эксперимента можно отметить следующее. Для линейных систем изложенные признаки применимы в границах , для нелинейных –
, h из (4). Границы t, в которых численные признаки достоверны, сокращаются с ростом h и размерности n, увеличиваются с точностью приближенных методов.
Заключение
Изложены компьютерно-ориентированные методы анализа устойчивости на основе приближенного решения дифференциальных систем, применимые к нелинейным, автономным и линейным системам. Методы включают необходимые и достаточные условия устойчивости с учетом класса задач, используются для построения численных признаков и для аналитического исследования. Предлагается совокупность способов одновременного компьютерного контроля устойчивости по ходу решения системы в реальном времени.