Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

COMPUTER STABILITY ANALYSIS OF SOLUTIONS OF ORDINARY DIFFERENTIAL EQUATIONS ON THE BASIS OF RECURRENT CONVERSION OF DIFFERENCE SCHEMES

Katrich S.A. 1
1 Taganrog Institute of A.P. Chekhov (a branch) FSBEE HVT «RostovStateEconomicUniversity (RINC)»
The approach to computer analysis of solution stability of Cauchy problem for systems of ODE in normal form is based on the recurrence of difference schemes mainly the method of Euler-Cauchy. The stability conditions are given in the form of estimates of the asymptotic behavior of convergent infinite products. In practice, the approximation of infinite products is performed by partial ones that cyclically are implemented in software with the aim of checking computer stability. The behavior of the partial products determines criteria of computer analysis: their limitations corresponds to stability, unlimited growth – instability, aspiration for zero with the increase of the independent variable – asymptotic stability. When computer realizes analysis as the method of finite difference schemes in addition to the Euler-Cauchy, methods of fourth order of Runge-Kutta and Adams are used. The proposed approach and data stability conditions apply for computer evaluation of stability solutions of ODE systems of general form. The results of numerical experiments and computer modeling and their analytical interpretation are given.
ordinary differential equations
Cauchy problem
stability
asymptotic stability
difference methods
infinite products
computer analysis of stability
numerical experiments

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

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

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

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

, , (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].

Рецензенты:

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

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