Электронный научный журнал
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,737

КОМПЬЮТЕРНЫЙ АНАЛИЗ УСТОЙЧИВОСТИ РЕШЕНИЙ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ НА ОСНОВЕ РЕКУРРЕНТНЫХ ПРЕОБРАЗОВАНИЙ РАЗНОСТНЫХ СХЕМ

Катрич С.А. 1
1 Таганрогский институт имени А.П. Чехова (филиал) ФГБОУ ВПО «Ростовского государственного экономического университета (РИНХ)»
Подход к компьютерному анализу устойчивости решений задачи Коши для систем ОДУ в нормальной форме строится на основе рекуррентных преобразований разностных схем, преимущественно используется метод Эйлера-Коши. Условия устойчивости даны в форме оценок асимптотического поведения сходящихся бесконечных произведений. На практике выполняется аппроксимация бесконечных произведений частичными, которые циклически реализуются программно с целью компьютерной проверки устойчивости. Характер поведения частичных произведений обусловливает критерии компьютерного анализа: их ограниченность соответствует устойчивости, неограниченный рост – неустойчивости, стремление к нулю с ростом независимой переменной – асимптотической устойчивости. При компьютерной реализации анализа в качестве разностных схем помимо метода Эйлера-Коши используются методы Рунге-Кутта и Адамса четвертого порядка. Предложенный подход и данные условия устойчивости применяются для компьютерной оценки устойчивости решений систем ОДУ общего вида. Приводятся результаты численных экспериментов и компьютерного моделирования, даны их аналитические интерпретации.
обыкновенные дифференциальные уравнения
задача Коши
устойчивость
асимптотическая устойчивость
разностные методы
бесконечные произведения
компьютерный анализ устойчивости
численные эксперименты
1. Арнольд В.И. Математические методы классической механики. – М.: Наука, 1989. – 472 с.
2. Демидович Б.П. Лекции по математической теории устойчивости. – СПб.: Изд-во «Лань», 2008. – 480 с.
3. Зубов И.В. Методы анализа динамики управляемых систем. М.: Физматлит, 2003. – 224 с.
4. Катрич С. А. Разработка и исследование программного моделирования устойчивости решений нелинейных дифференциальных уравнений на основе разностных методов: Автореф.дис. канд. техн. наук. – Таганрог, 2006. – 20 с.
5. Ромм Я.Е. Мультипликативные критерии устойчивости на основе разностных решений обыкновенных дифференциальных уравнений // Кибернетика и системный анализ. – 2006. – № 1. – С. 127 – 142.
6. Ромм Я.Е. Моделирование устойчивости по Ляпунову на основе преобразований разностных схем решений обыкновенных дифференциальных уравнений // Известия РАН. Математическоемоделирование. – 2008. – Т. 20. – №12. – С. 105 – 118.
7. I.V. Astashova. Application of dinamical systems to the study of asymptotic properties of solutions to nonlinear higher-order differential equations // Journal of Mathematical Sciences. – 2005. – Vol. 126. – № 5. – P. 1361 – 1390.
8. B. Cannas, F. Pisano.A Piecewise Linear Approximation Method for the Evaluation of Lyapunov Exponents of Polynomial Nonlinear Systems // Chaos and Complex Systems. – Springer-Verlag Berlin Heidelberg. – 2013. – P. 439 – 447.
9. Maria Alice Bertolim, Alain Jacquemard. Time switched differential equations and the Euler polynomials // Fondazione Annali di Matematica Puraed Applicata and Springer-Verlag Berlin Heidelberg 2013.
10. K. Wright. Adaptive methods for piecewise polynomial collocation for ordinary differential equations // BIT Numerical Mathematics, 2007, 47. – P. 197 – 212.

Постановка задачи

Традиционно анализ устойчивости по Ляпунову – предмет качественной теории дифференциальных уравнений. Теоретические методы сложны, трудоемки и, как правило, не допускают компьютерной реализации. В работе предлагается подход к компьютерному анализу устойчивости на основе разностной схемы Эйлера-Коши решений систем обыкновенных дифференциальных уравнений (ОДУ) в нормальной форме, который сводит вопрос об устойчивости к исследованию асимптотического поведения бесконечных произведений, построенных на основе метода Эйлера-Коши. С целью компьютерной реализации анализа устойчивости выполняется циклическое накопление частичных произведений, значения которых при фиксированном шаге разностной схемы пробегают достаточно большой промежуток на полуоси. При данном подходе компьютерный анализ устойчивости сводится к заданию правых частей систем ОДУ и вычислению частичных произведений, попутно получается приближенное решение системы.

Описание общей схемы компьютерного анализа устойчивости

Пусть требуется исследовать устойчивость в смысле Ляпунова решения задачи Коши для нелинейной системы ОДУ в нормальной форме

, , (1)

где , – искомое решение, – вектор начальных данных, – заданная вектор-функция от переменных: независимой переменной и зависимых переменных , , …, .

Предполагается, что для системы (1) выполнены все условия существования и единственности решения на всей полупрямой . Предполагается, что эти же условия выполнены для всех решений с начальными условиями , если только

, (2)

где – некоторое постоянное число. Здесь и в дальнейшем под знаком понимается каноническая норма вектора.

Пусть – пространство точек , или,. Пусть множество всех точек и в -мерном пространстве , где – всевозможные решения, получаемые из (2), при .

На множестве предполагаются выполненными условия:

1) функция определена и непрерывно дифференцируема по на при и всех из (2);

2) выполнено условие Липшица

, ; (3)

3) существует константа , такая что

. (4)

Определение устойчивости по Ляпунову заимствуется из [2] с некоторыми упрощениями, допустимыми в рассматриваемых условиях. Решение устойчиво (справа), если для любого сколь угодно малого числа существует , , такое что влечет . Решение асимптотически устойчиво (справа), если оно устойчиво и найдется , , такое что влечет . Всюду ниже рассматривается устойчивость справа (слева аналогично), которая для краткости называется просто устойчивостью.

Базовая схема излагаемого способа компьютерного анализа устойчивости строится на основе разностных методов приближенного решения ОДУ, ниже для этой цели используется метод Эйлера-Коши. Существенной особенностью при этом является выбор шага численного интегрирования : предполагается, что для каждого произвольного , каково бы ни было ,

, , (5)

изменение переменной рассматривается как изменение правой границы промежутка . Иными словами, на любом промежутке , но шаг остается равномерным внутри промежутка.

В [4] на основе тейлоровских разложений невозмущенного решения и возмущенного решения задачи Коши (1) показано, что для разности между соответственными компонентами возмущенного и невозмущенного решений для любого из (5) имеет место равенство

,(6)

где величины суть и

,(7)

где здесь и всюду ниже всегда предполагается, что , .

Соотношение (6) является базовым для формулируемых в рамках предложенной схемы условий устойчивости. Первое слагаемое правой части представляет собой главную часть роста возмущения, поскольку не содержит множителя . В [4] на основании условий (3), (4) доказывается, что при, что равносильно , для любого из (5). Предельный переход в равенстве (6) влечет

,, (8)

где

. (9)

Аналогичные соотношения на основе метода Эйлера приводятся в [5, 6].

Пусть , где из (9), из (7), . Непосредственно на основании (8) формулируются условия устойчивости в виде следующей теоремы.

Теорема 1. В условиях 1) – 3) для устойчивости по Ляпунову решения задачи (1) необходимо и достаточно, чтобы существовало , , такое, что одновременно для всех решений при ограничении выполняется условие

, . (10)

Для асимптотической устойчивости необходимо и достаточно, чтобы выполнялось (10) и нашлось , , такое, что влечет

. (11)

Условия (10), (11) реализуются программно. Для этого выполняется циклическое накопление частичных произведений , , поведение этих произведений будет определять характер устойчивости решения: если при неограниченном росте будет наблюдаться ограниченность произведений, это будет означать устойчивость, стремление к нулю – асимптотическую устойчивость, неограниченность – неустойчивость. Такое моделирование составляет существо предложенного компьютерного анализа устойчивости.

На практике бесконечные произведения не могут быть вычислены точно. Моделирующая их поведение программа с необходимостью остановится на их приближении конечным числом сомножителей. Возникает вопрос, как такое приближение отразится на достоверности оценки устойчивости – необходимо исследовать поведение левых частей соотношений (10), (11) в зависимости от количества сомножителей. Решение этого вопроса дано в [4]. Утверждение теоремы 1 переносится на случай частичных произведений с некоторыми оговорками, в частности, в формулировку теоремы добавляется, помимо требования существования , требование существования номера (шага ), начиная с которого выполняется условие , . Из предельного перехода в неравенстве с учетом (10) следует, что это условие является достаточным для устойчивости решения задачи (1).Как показано в[4],при естественных ограничениях замена предельных значений в выражении условий (10), (11) на их конечные приближения сохраняет достоверность оценки устойчивости.

В [4] показана возможность компьютерного моделирования устойчивости при помощи циклического накопления произведений , приближаемых как по методу Эйлера-Коши, так и на основе методов Рунге-Кутта и Адамса, аналогичный анализ на основе метода Эйлера выполняется в [5, 6].

Численный эксперимент по моделированию устойчивости

Пусть при  рассматривается система

, , , (12)

общее решение которой ,,где ,,.

В [2] аналитически доказано, что тривиальное решение системы (12) неустойчиво. Ниже этот же результат получается на основе численного моделирования (код программы полностью приводится в [4]). Программные модели представляют собой стандартные подпрограммы Delphi.

В табл.1 результаты численного эксперимента даны при использовании метода Рунге-Кутта 4-го порядка при начальных данных , и для шага .

Таблица 1

Численное моделирование неустойчивости тривиального решения системы (12)

11.00

101.00

501.00

901.00

991.00

1.109E+0001

1.010E+0002

5.009E+0002

9.009E+0002

9.909E+0002

Норма произведения на промежутке возросла от значения 11 до 990, что соответствует неустойчивости тривиального решения.

Замечание 1. Из равенства (8) следует, что Поэтому, где при , Иными словами, аналогично норме произведения  ведет себя отношение разности возмущенного и невозмущенного решений к величине возмущений начальных данных, взятых по норме. Этот факт подтверждается численным экспериментом, результаты которого полностью повторяют данные табл. 1 [4 – 6].

Пусть дана система ОДУ

,(13)

где, . Общее решение системы (13) имеет вид.

Решения системы (13) с начальными данными асимптотически устойчивы, при условии – устойчивы, если – неустойчивы [3].

В табл.2 результаты численного эксперимента даны при использовании метода Эйлера-Коши для начальных данных , и с шагом.

Таблица 2

Численное моделирование асимптотической устойчивости решения системы (13)

10.00

100.00

400.00

900.00

1000.00

1.387E+0000

4.320E-0001

6.829E-0002

2.002E-0002

1.371E-0002

Норма произведения на промежутке убывает к нулю, что согласно (11) интерпретируется как признак асимптотической устойчивости. Стремление нормы к нулю более выражено в случае (например, и ).

В табл.3 даны результаты численного эксперимента при использовании метода Рунге-Кутта 4-го порядка для начальных данных , и с шагом .

Таблица 3

Численное моделирование устойчивости решения системы (13)

10.00

100.00

400.00

900.00

1000.00

1.680E+0000

1.749E+0000

1.803E+0000

2.286E+0000

1.926E+0000

Норма произведения на промежутке изменяется, но в целом не превосходит значения 2.3, что согласно условию (10) является признаком устойчивости.

В табл.4 приведены результаты численного эксперимента при использовании метода Адамса 4-го порядка с шагом для начальных данных , и .

Таблица 4

Численное моделирование неустойчивости решения системы (13)

1.00

10.00

40.00

47.00

54.00

1.746E+0000

2.514E+0000

2.077E+0001

1.017E+0002

1.393E+0004

Норма произведения на промежутке монотонно возрастает с большой скоростью, что в продолжение процесса приводит к переполнению, это – нарушение условия (10) и интерпретируется как неустойчивость. Возрастание более быстро приводит к переполнению при.

Для численных экспериментов в табл.2 – 4 можно повторить замечание 1.

Изложенный подход к компьютерному анализу устойчивости отличается от известных [1, 2] построением на основе рекуррентных преобразований разностных схем. Отличия сохраняются относительно современного состояния исследований [7 – 9], включая численное моделирование [10], отличие от [5, 6] заключается в использовании разностных методов более высокого порядка.

Заключение

Предложена схема компьютерного анализа устойчивости для систем ОДУ в нормальной форме на основе рекуррентных преобразований разностных приближений решений задачи Коши. Разность между значениями возмущенного и невозмущенного решений выражается как бесконечное произведение, сомножители которого включают в себя разностные приближения и возмущения начальных данных. Асимптотическое поведение таких произведений однозначно определяет характер устойчивости решения. Необходимым и достаточным условием устойчивости является равномерная ограниченность бесконечных произведений на полуоси. Асимптотическая устойчивость имеет место при стремлении этих произведений к нулю с ростом независимой переменной. Развернутый численный и программный эксперимент представлен в [4].

Рецензенты:

Веселов Г.Е., д.т.н., доцент, директор института компьютерных технологий и информационной безопасности ЮФУ, г. Таганрог;

Карелин В.П., д.т.н., профессор, заведующий кафедрой прикладной математики и информационных технологий ТИУиЭ, г. Таганрог.


Библиографическая ссылка

Катрич С.А. КОМПЬЮТЕРНЫЙ АНАЛИЗ УСТОЙЧИВОСТИ РЕШЕНИЙ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ НА ОСНОВЕ РЕКУРРЕНТНЫХ ПРЕОБРАЗОВАНИЙ РАЗНОСТНЫХ СХЕМ // Современные проблемы науки и образования. – 2014. – № 6.;
URL: http://science-education.ru/ru/article/view?id=15863 (дата обращения: 26.08.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.252