Электронная библиотека Веда
Цели библиотеки
Скачать бесплатно
Доставка литературы
Доставка диссертаций
Размещение литературы
Контактные данные
Я ищу:
Библиотечный каталог российских и украинских диссертаций

Вы находитесь:
Диссертационные работы России
Технические науки
Математическое моделирование, численные методы и комплексы программ

Диссертационная работа:

Катрич Сергей Анатольевич. Разработка и исследование схем программного моделирования устойчивости решений нелинейных дифференциальных уравнений на основе разностных методов : дис. ... канд. техн. наук : 05.13.18 Таганрог, 2006 220 с. РГБ ОД, 61:07-5/1212

смотреть содержание
смотреть введение
Содержание к работе:

ВВЕДЕНИЕ 4

ГЛАВА L РАЗРАБОТКА И ИССЛЕДОВАНИЕ ПРОГРАММИРУЕМЫХ
МУЛЬТИПЛИКАТИВНЫХ КРИТЕРИЕВ УСТОЙЧИВОСТИ РЕШЕНИЙ
СИСТЕМ НЕЛИНЕЙНЫХ ОДУ.
33

  1. Построение базовой схемы анализа устойчивости решений нелинейных ОДУ 34

  2. Синтез и обоснование мультипликативных критериев устойчивое і и решений нелинейных ОДУ 39

  3. Программное моделирование мультипликативных критериев устойчивости 44

1А Обоснование мультипликативных критериев устойчивости с конечным
числом сомножителей для программного моделирования 48

  1. Разностные оценки асимптотического поведения подущенного решения задачи Коши 52

  2. Критерии устойчивости для систем нелинейных ОДУ 62

  3. Об отличительных особенностях критериев устойчивое ж її случае систем линейных ОДУ 69

  4. Выводы 73

ГЛАВА 2. ИССЛЕДОВАНИЕ ДОСТОВЕРНОСТИ МУЛЬТИПЛИКАТИВНЫХ
КРИТЕРИЕВ УСТОЙЧИВОСТИ В ЗАВИСИМОСТИ ОТ ПОГРЕШНОСТИ
РАЗНОСТНЫХ СХЕМ.
75

  1. Перенос критериев устойчивости на разностные приближения решений ОДУ по методу Эйлера-Коши 76

  2. Достоверность мультипликативных критериев устойчивости на основе метода Эйлера-Коши в случае конечного числа сомножителей 79

  3. О построении мультипликативных критериев устойчивости на основе методов Эйлера, Руше-Кутта и Адамса 84

  4. Программное моделирование мультипликативных критериев устойчивости для представленных разностных схем 85

  5. Достоверность мультипликативных критериев устойчивости па основе метода Эйлера-Коши в случае систем нелинейных ОДУ 93

  6. Оценка накопления погрешности метода Эйлера-Коши в условиях устойчивости 96

  7. Выводы 99

ГЛАВА 3. ПРОГРАММНОЕ МОДЕЛИРОВАНИЕ ПО ПРИМЕНЕНИЮ
МУЛЬТИПЛИКАТИВНЫХ КРИТЕРИЕВ УСТОЙЧИВОСТИ И ЧИСЛЕННЫЙ
ЭКСПЕРИМЕНТ ПО ОЦЕНКЕ ИХ ДОСТОВЕРНОСТИ
102

3.1. Построение основной программной модели и реализация компьютерного
анализа устойчивости решений систем нелинейных ОДУ 102

  1. Инвариантность программной модели анализа устойчивости решений систем нелинейных ОДУ относительно размерности системы 123

  2. Моделирование мультипликативных критериев устойчивости и оценка их достоверности в зависимости от величины шага разностной схемы 126

  3. Моделирование мультипликативных критериев устойчивости и оценка их достоверности в зависимости от величины возмущений начальных данных 133

3.5. Выводы 139

ЗАКЛЮЧЕНИЕ 141

ЛИТЕРАТУРА 145

ПРИЛОЖЕНИЕ 1 152

ПРИЛОЖЕНИЕ 2 157

ПРИЛОЖЕНИЕ 3 162

ПРИЛОЖЕНИЕ 4 166

ПРИЛОЖЕНИЕ 5 176

ПРИЛОЖЕНИЕ 6 196

ПРИЛОЖЕНИЕ 7 217

Введение к работе:

Акгуалмюсіь проблемы. Понятие устойчивости движения, устойчивости динамических процессов зародилось и получило значительное развитие главным образом в механике задолго до появления классических и современных теорий устойчивости и управления.

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

В настоящее время число понятий устойчивости настолько велико, что этот термин справедливо считается перегруженным. Например, под устойчивостью системы автоматического регулирования обычно понимшог свойство системы возвращаться к первоначальному состоянию после прекращения внешнего воздействия [1]. Полагая, что система автоматического регулирования описывается системой обыкновенных дифференциальных уравнений (ОДУ) приходим к понятию устойчивости движения (решения) в смысле Ляпунова. Однако и в математической теории устойчивости данный термин обозначает весьма обширные качества и свойства решений ОДУ. Помимо устойчивости по Ляпунову различают устойчивость по Лагранжу, Пуассону (Пуанкаре), орбитальную устойчивость, экспоненциальную устойчивость, устойчивость порядка т по Бирктофу [2] и некоторые др. (определения некоторых понятий устойчивости приводятся ниже).

Кроме тою, в численных методах существует понятие устойчивости численных схем, в классической механике говорят об устойчивости равновесия, Солнечной системы, п бурно развивающейся отрасли современной математики, теории особенностей дифференцируемых отображений, помимо просто устойчивости гладких отображений различают устойчивость ннфиннтезн-мальиую, деформационную, дифференцируемую, лаграижеву, лежандрову,

5 тоїюлоіическую и некоторые др, [3],

Непосредственно возникновение .математической теории устойчивости движения связано с именем A.M. Ляпунова. Основы ею теории были разработаны более 100 лет назад («Общая задача об устойчивости движения», 1892). С начала 30-х годов теория устойчивости Ляпунова получила интенсивное развитие вследствие появления новых задач науки и техники, В настоящее время теория устойчивости по Ляпунову является общепринятой и применяется во многих областях естествознания.

Предлагаемая диссертационная работа посвящена именно анализу устойчивости по Ляпунову (ниже устойчивости).

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

Исследование устойчивости - предмет качественной теории дифференциальных уравнений. Ляпуновым предложены два метода для анализа устойчивости.

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

развит в фундаментальных трудах [4 - 16] и в трудах многих других специалистов по нелинейному анализу [17].

Второй (качественный, прямой) метод Ляпунова в основном связан с введением вспомогательных оценивающих функции, получивших его имя. Этот метод - результат синтеза идей качественной теории дифференциальных уравнений, идей А. Пуанкаре и A.M. Ляпунова [17].

Второй метод стал основным методом классической и современной теории устойчивости, а также качественной теории дифференциальных уравнений. Он эффективен для разнообразных квазилинейных и нелинейных систем.

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

Практикуемые в настоящий момент подходы к автоматизации анализа устойчивости опираются на вычисление функций Ляпунова [18 - 20] и применение схем символьной обработки [21]. Эти подходы продолжают быть связанными с традиционными методами математического анализа устоичивосіи на основе первого и второго методов Ляпунова [2, 22 - 25] или, в случае систем линейных ОДУ с постоянными коэффициентами, с исследованием корней характеристическою полинома матрицы коэффициентов на основе аппарата линейной алгебры [26, 27], Подходы сложны, трудоемки, предлаїаемьіе решения зависят от вида правой части и начальных данных, по крайней мере, в случае нелинейных ОДУ- Такая зависимость обусловлена тем, что не сущее г-вует единообразного алгоритма построения функций Ляпунова для нелинейных ОДУ общего вида.

7 Тем самым известные подходы к компьютерному анализу устойчивости

на данный момент не отвечают качествам технологичности.

Отсюда актуальна задача разработки метода компьютерного анализа устойчивости решений нелинейных ОДУ общего вида на основе разностных приближений, который не использовал бы аналитические решения и аппарат функций Ляпунова.

Существующие методы анализа устойчивости традиционно включаю г аішиз устойчивости на основе определения, первый метод Ляпунова, устойчивость линейных ОДУ, второй метод Ляпунова.

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

=F(tJ)y

0)

где Y = Y(t), Y = (у{(l),y2(/),.„,;'*('))- искомая вектор-функция, У0=(-Кі('о)^г('о)'"-^п('о))-веІСГОРнаііальньїх данных, F(i,Y)- (/Д^К),/,(^),...,/^(/,К)) - заданная векгор-функция от и+ 1 переменных; независимой переменной / и п зависимых переменных )^(/),

^,(/),...,^(0.

Функции /,, / = 1,2,.,.,/г, предполагаются определенными на некотором множестве S, элементами которого служат точки {^У^У^^-^У,, )

Ниже приводятся некоторые определения устойчивости и анализируются их различия.

Если F(t4Y) непрерывна на множестве S точек (/,Г) и S открыто

справа от /0, если решение К(/)=К(/;/0, KD) системы (1) существует па полупрямой /0 (t,Y(t))eS для всех />/0, то Y(i) называется устойчивым по Ляпунову (справа) [2], если выполнены следующие условия:

1) Найдется такое h>0, что каждое решение Г(/; /0, Г,) существует па

/э (t,Y(t))eS для всех />/0, если для начального вектора выполне-

8 но неравенство ' Yl - YQ ' < b.

2) Если дано є>0, то можно найти 5 = 5(c;F,rJ, 0<5<Л, іакое, что

Y}-Y0 <5 влечет Y{t\t^Yx)-Y{t\t^ Г0)]]<с для всей полупрямой

/,

Устойчивость по Ляпунову по сути дела является равномерной непрерывностью Y(i;t09 Y0) на [/0,о) по отношению к начальному вектору К0.

Иными словами, решение У(?;ґ00) устойчиво, если достаточно близкие к нему при любом /0 решения Г(/;/0, К,) целиком пофужаются в сколь угодно узкую є-трубку, построенную вокруг решения Y(t;t0,YQ).

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

ра, например, і, К = )\ + у ,

+ ...+

У»

Решение Y(i; ?0,17о) называется асимптотически устойчивым (справа) [2], если при выполнении 1) и 2) выполняется также условие:

3) Можно найти такое 50 =S0(/r»K0), 00<, что из Y}-Y0 <5;|

следует Yint^Y^-Yint^Ytl^QnpHt-xo.

Существуют разновидности [22, 28] данных определений, где б и б„ зависят от є, ?0.

В [29] впервые было введено определение устойчивости, равномерной по /0: в определении устойчивости положительное число 6 выбирается не зависящим от /0. В случае, если 6 не зависит и от /0 и от Y0, то решение Y{t)=K(/;f0, YQ) называется просто равномерно устойчивым (справа) [2].

Аналоїично определяется понятие равномерной асимптотической устойчивости.

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

Решение Y(t)=Y(t;t0i YG) системы (1), не являющееся устойчивым, на-

зывается неустойчивым. Не приводя определения неустойчивости на языке с, 5, следует отметить, что из отрицания определения устойчивости, неустойчивым следует считать решение К(/;/0, Y0) непродолжаемое при f->oo, i.e. не

существующее на [^,^), или такое, для которою в любой окрестности точки К0 найдется точка Г1? порождающая при t = tG решение К(/;/0, К,), непродолжаемое на [/и?о) [22].

Кроме тою, сели решение К(/;г0, К0) (а<г<со, а -число или а-~х)

системы (1) с непрерывной правой частью устойчиво при каком-нибудь /с є(д,со), то оно будет устойчиво при любом друюм t\ є(«,со) [22]. Таким

образом, можно ограничиться проверкой устойчивости решения, а также ею асимптотической устойчивости, лишь для некоторого заданного / = /0. Отсюда

также следует, что если решение Г(г;/00) (#

/ =/0, то оно является неустойчивым для любого друюі о t[ є(д,оо).

На этой основе в теоремах устойчивости tQi как правило, считается фиксированным [225 30 - 32],

Если в определении устойчивости (асимптотической устойчивости)

r(f;/05rt)-r(r;r05ru)J заменить на \yx(t;tQ, Y})-y}(t\t^ Г0)| +

+ \у2U; /D, I7, )-;s(/;/0, К0)| + ... + |^(/; f0, ^,)-^,(^,^)1 при некотором

m9 то получится соответствующее определение устойчивости (асимптотической устойчивости) решения системы (1) относительно части переменных ^JV-'>J'„ [28,33].

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

Равномерная устойчивость относительно части переменных изучалась в [34], дальнейшее развитие її приложение к механике - в [25, 35]. Современное сопояпие теории устойчивости и стабилизации движения по отношению к части переменных дано в [36]-

Рассмотренные выше понятия устойчивости связаны с изменением только начальных условий. В [37] было впервые введено понятие устойчивости при постоянно действующих возмущениях, которое затем было обобщенно в [38]. Данный вид устойчивости учитывает изменение самой правой части (1). Такой вид устойчивости приложим к большинству практических задач, в которых встречаются возмущения, действующие не только в начальный момент времени, но и во время движения.

Если устойчивость при постоянно действующих возмущениях является равномерной по /0, то ее называют тотальной.

Устойчивость и асимптотическая устойчивость по Ляпунову являются качественными характеристиками поведения систем па бесконечном интервале времени. Они не всегда могут характеризовать поведение систем. Большинство технических систем функционируют в течение конечного промежутка времени, и при этом представляет интерес не только факт устойчивости или неустойчивости, но и количественные оценки их поведения, а также приемлемость их моделей в реальных условиях.

Первые исследования устойчивости на конечном промежутке времени даны в [39 - 41]. В дальнейшем они получили развитие как теория практической устойчивости, например, в [42 - 45],

Пусть система (1) обладает свойством единственности решения Y(t;t09 Yv)> где /0є(й,со) и У0 принадлежит некоторой открытой области

действительного //-мерного векторного пространства.

Система (1) называется устойчивой по Лагранжу [22, 43, 46], если:

1) Каждое решение Y(t;t0, Г0), где /0 є(д,со), неограниченно продол
жаемо вправо, т.е. имеет смысл при /0

2) К(/;/00) оіраниченана [f0,co).

Движение называется устойчивым по Лагранжу, если его траектория вечно остается в ограниченной области фазового пространства [47].

Требование ограниченности решений приводит к разумному и естественному типу устойчивости, который находит свою область применения в за-

31 дачах небесной механики [25].

Если произвольное движение механической системы бесконечно часто

возвращается к своему начальному состоянию, то это свойство механических

систем называется устойчивостью по Пуассону (Пуанкаре) [2].

Понятие орбитальной устойчивости есть фактическая адаптация устойчивости в смысле Ляпунова к свойствам периодических решений автономных систем, поскольку из устойчивости решения по Ляпунову следует ею орбитальная устойчивость. Однако из орбитальной устойчивости решения, вообще говоря, не вытекает устойчивость его по Ляпунову, а тем более асимптотическая устойчивость [22].

Понятие экспоненциальной устойчивости формализует свойство асимптотически устойчивых решений, содержащих экспоненты с отрицательными показателями степени, заключающееся в стремлении к нулю на положительной полупрямой при аргументе, стремящемся к бесконечности, и имеет тог смысл, что стремление к нулю всех решений не является достаточным условием для асимптотической устойчивости.

Решение Y(t\t0, Y0) системы (1) называется экспоненциально устойчивым при f —>со [22, 32], если для каждого решения Y(t;t0, К,) этой системы в некоторой области ^„h = со, справедливо нера-

венство Y(ntQ,Y})-Y(t;i^Y,)\\(/>/ДгдеЛ/ иа-положительные постоянные, не зависящие от выбора решения Г(/;/0> У,).

Очевидно, что из экспоненциальной устойчивости следует асимптотическая- Однако из асимптотической устойчивости, вообще говоря, не следует экспоненциальная устойчивость за исключением случая однородных линейных систем с постоянными коэффициентами [22]. Отсюда понятие асимптотической устойчивости является более общим.

Из приведенных понятий устойчивости в смысле различных определений видно, что наиболее формально общим и математически абстрактным является определение устойчивости в смысле Ляпунова. Другие определения либо являются адаптацией устойчивости по Ляпунову в приложениях к кон-

12 кретиым задачам, как, например, устойчивость по Пуассону и Лагранжу применительно к задачам небесной механики, либо выражают свойства решений систем некоторого частного случая - орбитальная и экспоненциальная УСТОЙЧИВОСТИ.

В дальнейшем рассматривается только устойчивость по Ляпунову, которая для краткости называется устойчивостью.

Пусть дана действительная дифференциальная система (1), где вектор-функция F{t„Y) в некоторой области R = {a <со,Ке D) (а -число или символ -со, D - открытое множество действительного евклидова «-мерного пространства 9IJ) непрерывна по независимой переменной / и имеет непрерывные частные производные первого порядка по зависимым переменным >',(/)* >':(')» > уЛО- Тогда для каждой точки (/00)є Я справедлива локальная теорема существования и единственности решения К(/)=К(/;/ОЇ YQ) системы (і) с начальными условиями У(/0;/0, К0) = К0.

Пусть г|-П(0 (f0"'; то ) ~ решение системы (1) (невозмущенное движение), устойчивость которого требуется исследовать, причем //-окрестность этою решения такова, что Ulf{x]{t))cD при / є[ґ0,со), где

ии(ц(!))={10{\<Н<<ю}.

Положим X ~ К-Г|(/), т.е. X есть отклонение решения У отрешения г|(/). Так как r\(t) = F(tir\(t)), то для X получаем дифференциальное уравнение

Т^М). (2)

где Fl(l,X)= F(t,X +r|(/))-F(/,r|(0) и удовлетворяет условиям существования и единственности решения в некоторой области Z = {aя, [Х[<И}Л причем, очевидно, /^(/,0)=0. Следовательно, система (2) имеет тривиальное решение Х = 0, которое в пространстве ЭД* соответствует данному решению г) = п(0- Система (2) называется приведенной [22], по Ляпунову она называется системой уравнений возмущенного

движения [28].

Таким образом, исследование устойчивости решения г| = П(0 в пространстве 91" сводится к исследованию устойчивости тривиального решения (положения равновесия) X = 0 в пространстве 91^ -

Естественно переход от системы (1) к системе (2) не всегда возможен, так как для отого необходимо знать y\(t) в явном виде, и не всегда дает преимущество, потому что может случиться, что система (2) более сложна, чем система (1), Например, если система (1) автономна, то система (2), вообще говоря, не является автономной [25],

Следует отмстить, что компьютеризация анализа устойчивости основанная непосредственно на ее определении практически невозможна, хотя бы по причине тою, что в подавляющем большинстве случаев точное аналитическое решение системы (1) неизвестно.

Устойчивость линейных ОДУ. Пусть дана система линейных дифференциальных уравнений

(3)

= A(t)V+F(t)9

гдеУ =

У,(0

,A{t) =

4,(0-^(0 ^(О—^ДО

aJO-Mjt)

Однородная система вида

(4)

называется соответствующей неоднородной системе (3).

Решения систем линейных дифференциальных уравнений обладают тем качеством, что они либо все одновременно устойчивы, либо все неустойчивы, в отличие от случая нелинейных уравнений, некоторые решения которых могут быть устойчивы, а другие - неустойчивы [22],

Отсюда линейные системы называются [22] устойчивыми или неустойчивыми.

Исследование устойчивости решений системы (3) можно заменить ио следованием устойчивости решений соответствующей однородной системы (4), поскольку для устойчивости системы (3) с непрерывными коэффициентами при любом свободном члене F(t) необходимо и достаточно, чтобы было устойчиво тривиальное решение Г = 0 системы (4) [22],

Следовательно, система линейных дифференциальных уравнений с непрерывными коэффициентами устойчива, когда устойчиво хотя бы одно решение этой системы, и неустойчива, если неустойчиво некоторое решение её.

Данные свойства решений систем (3) справедливы также по отношению к равномерной и асимптотической устойчивостям [22].

Линейные системы обладают качеством, которого нет в случае нелинейных уравнений: существует взаимозависимость между устойчивостью и ограниченностью. Именно, система (4) с непрерывными коэффициентами устойчива (справа) тогда и только тогда, когда каждое решение У = Y(t) этой системы ограничено справа. Если линейная неоднородная система с непрерывными коэффициентами устойчива (справа) и одно из ее решений ограничено справа, то и все остальные ограничены справа. Если решения линейной неоднородной системы с непрерывными коэффициентами ограничены справа, то они устойчивы (справа) [2].

Для нелинейной системы из ограниченности сё решений не следует их устойчивость. Например, в случае системы

dt dy>

= ~Уг^Уу +У2 >

= Уірі +У2 .

общее решение которой имеет вид

\)\=cicos(c]t + c2), [v, =cisin(c-,/ + t'2),

где с\, с, произвольные постоянные, из очевидной ограниченности решений

не следует их устойчивость, поскольку нулевое решение рассматриваемой системы не асимптотически устойчиво, а все остальные неустойчивы [2].

15 Еще один отличительный аспект линейных уравнений заключается в

том, что система (4) с непрерывными коэффициентами асимптотически устойчива тогда и только тогда, когда все её решения Y = Y(i) стремятся к нулю

при / -> со, т.е. НтГ(/) = О,

Для нелинейной системы стремление к нулю всех решений не является достаточным условием асимптотической устойчивости тривиальною решения [22].

Пусть дана система

dt t dt t

СІУі Уі ,2

= —-' УхУг *

допускающая тривиальное решение yx=, y2 -0. Общее решение имеет вид

v\ =—, t

в частности, при /0 =1,

}'ХО = }\(1а)!е

-yfUuMt-U

уг{0=

Л('о)

Очевидно, >'i(/)^->0 и y2(t)-^ при г-^со. Однако для любою 5>0 при

>4У т

О J в

v,(/0) = S:, y2(t0) = 5 будет иметь место неравенство у

Следовательно, решение у,=0, у2=0 не является устойчивым и тем более

асимптотически устойчивым при t—>.

Наиболее простым и изученным является случай сне гем линейных дифференциальных уравнений с постоянными коэффициентами. Пусть в (4) Л -матрица постоянных вещественных коэффициентов. Тогда

(5)

Исследование устойчивости системы (5) сводится к выяснению знаков действительных частей собственных чисел (характеристических корней) чаї-рицы А. Точней, линейная однородная система с постоянными коэффициентами устойчива (асимптотически устойчива) тогда и только тогда, когда все характеристические корни ее матрицы имеют неположительные (отрицательные) вещественные части, причем характеристические корни, имеющие нулевые вещественные части, должны быть простыми [22],

Таким образом, согласно классической теории, для суждения об устойчивости линейной однородной системы дифференциальных уравнений с постоянными коэффициентами следует знать характер расположения корней характеристического уравнения на комплексной плоскости корней. Теоретически эта задача леї ко решаема, практическая же реализация этою в сущности простого подхода сопряжена с рядом трудностей, связанных с решением «векового» уравнения det(Л-ХЕ)= 0. Трудности эти обусловлены развертыванием характеристического многочлена, в процессе которого моїут оказаться возмущенными его коэффициенты. Вычисление корней многочлена на практике можеі повлечь неустойчивость вследствие возмущения коэффициентов [48]. Эти трудности усугубляются с ростом размерности системы. Непосредственный подход к решению алгебраической проблемы собственных значений, а, следовательно, и вопроса об устойчивости, применяют лишь при малой размерности матрицы А (я = 2,3); уже при п>А на первый план выходят специальные численные методы решения таких задач [49].

Критерий Гурвица дает необходимые и достаточные условия, при которых алгебраическое уравнение с действительными коэффициентами имеет корни лишь с отрицательными вещественными частями. Пусть

Р{\) = а}:+а„ .Г"1 +... + А0 (6)

характеристический многочлен матрицы А и д0>0, а„^0, и>1. Полином

вида (6), не имеющий, очевидно, нулевых корней, называется стандартным полиномом. Квадратная матрица порядка п

G =

a, я00 ...О

д^ а2 а} ,., О а5 я4 ау.,, О

0 0 0 ...а

л /

по їлавной диагонали которой расположены коэффициенты а,, я,,...,и л,

вправо по строке от этих элементов - коэффициенты с убывающими номерами, влево - с возрастающими, при этом полагается а,=0,если /<0 или />иэ

называется матрицей Гурвица. Для того чтобы действительные части всех корней характеристического многочлена (6) матрицы А были отрицательными, соответственно, система (5) была бы асимптотически устойчива, необходимо и достаточно, чтобы все главные диатнальные миноры матрицы G были бы положительны (условие Гурвица). Чтобы стандартный полином вида (6) имел нули, лежащие лишь в замкнутой левой полуплоскости Re?. <0, ішобхо-

димо и достаточно, чтобы все главные диагональные миноры ею матрицы Гурвица были бы неотрицательны [22].

Если степень полинома Р(Х) сравнительно большая, то применение

критерия Гурвица становится затруднительным ввиду необходимости подсчета определителей высоких порядков. В этом случае для определения расположения корней полинома Р(\) на комплексной плоскости иногда оказываются

более удобными геометрические признаки [22], эквивалентные критерию Гурвица.

С анализом устойчивости линейных систем связан первый метод Ляпунова, с построением функций Ляпунова связан второй его метод. Следуя [17], охарактеризуем вначале второй метод.

Второй (прямой) метод Ляпунова. Система дифференциальных уравнений возмущенною движения переписывается в виде

-**>

(7)

где анализируется устойчивость невоімущенного решения л'-0 (/(/,0)^0)-

18 Здесь х, как и прежде, рассматривается как точка вещественною п-мерного

пространства ЯГ с нормой ' х \= х{ + .vJ + ...+ хл , или евклидова простран-

ства Е" с нормой ^1 = -/^/+^+.,, + ^. Вещественная вектор-функция

/(/,х) определена и непрерывна в области R-{(t,x): х 0}

(6 = const >0 или = со) и имеет в ней непрерывные частные производные по лРм.,д*я, которые ограничены в каждой замкнутой обласіи

R = {(t9x)s R:\x\Qtt>Q}i 00Это обеспечивает существование,

единственность и нелокальную продолжимость решений системы, непрерывную зависимость их от начальных данных (и t9 трактуемого как время) в области R, 11ри b - со предполагается продолжимость решений при всех г е[05со) так, что возмущенные решения описываются векгор-функцией

А'(/;/0,А"0), Определенной И Непрерывно Дифференцируемой при (/00)є/?,

/є(/сїГ)^0<Г<«.

Для исследования устойчивости Ляпунов [28] ввел вещественные скалярные функции V(t,x), определенные и непрерывно дифференцируемые в

Д,1'(/,0)=0.

Первой производной функции V по времени /, взятой в силу уравнений возмущенного движения (7), называется функция

at ot ax

которая тождественно равна нулю при л=0,т.е, V(t,Q)=0,

Вводятся следующие определения основных свойств функции V. Постоянно положительная (соответственно постоянно отрицательная) функция V{ttx) удовлетворяет неравенству К(/,д*)>0 (соответственно

V{t9x)<0) при (/,д*)єЛ- Функция У{і,х) называется определенно положительной, если существует функция W(x) такая, что V(t9x)>tV(x)>Q при (t,x)eR, х\ >0, Функция V[t9x) допускает бесконечно малый высший предел, если V непрерывна по .v в нач&пе координат равномерно относительно /

и к(/,0) = 0,т.е, Vc>0 36-5(c)>0:!jc[]<5=>|r(^x)|0,оэ).

Классическими являются теорема Ляпунова об устойчивости движения с обращением в [50] и теорема о равномерной устойчивости [29] с обращением в [51 -53], которые приводятся здесь в единой формулировке.

Для устойчивости (соответственно устойчивости равномерной по /0) невозмущенного решения .v = 0 системы (7) необходимо и достаточно, чтобы в некоторой области R существовала непрерывно дифференцируемая определенно положительная (соответственно и допускающая бесконечно малый высший предел) функция V(f,x), производная которой в силу уравнений возмущенною движения (7) постоянно отрицательна (К(/,л-)< 0),

Классическими теоремами о неустойчивости являются две теоремы Ляпунова, первая из которых необратима для произвольной системы (7), как показано в [32], а обращение второй доказано в [52, 54, 55], Они приведены здесь в единой формулировке.

Для неустойчивости невозмущенного решения a- = 0 системы (7) достаточно (соответственно необходимо и достаточно), чтобы в некоторой области R существовала непрерывно дифференцируемая функция K(/Vv), допускающая бесконечно малый высший предел (соответственно ограниченная) и при любом t0 принимающая положительные значения в некоторых точках сколь

угодно малой окрестности начала координат х-0, производная которой в силу уравнений возмущенного движения (7) определенно положительна (соответственно удовлетворяет неравенству (^(/,.v)>?.F(r,x), где X - некоюрое положительное число), (t,x)e R.

Георема об асимптотической устойчшюсти дана Ляпуновым [28], в [56, 57] показано, что при ее условиях имеет место равномерная асимптотическая устойчивость, и предложено ее обращение.

Для равномерной асимптотической устойчивости решения ,v = 0 системы (7) необходимо и достаточно, чтобы в некоторой области R существовала непрерывно дифференцируемая положительно определенная функция К(/,х),

20 допускающая бесконечно малый высший предел, производная которой в силу

системы (7) определенно отрицательна.

В приведенных теоремах предполагается существование функции Ляпунова. Возникает вопрос: существует ли в действительности такая функция, т.е. если решение системы (7) обладает каким-либо свойством устойчивости, го можно ли построить вспомогательную функцию, удовлетворяющую соответствующей теореме?

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

Во-вторых, в ряде случаев устойчивость можно исследовать, рассмотрев сначала упрощенную систему. В случае если устойчивость последней может быть легко установлена, то на основе обратной теоремы делается вывод о су-шествовании соответствующей вспомогательной функции. При определенных предположениях ее можно также использовать в качестве подходящей вегю-моїательной функции для исходной системы, Такой подход типичен для исследования устойчивости при постоянно действующих возмущениях и доказательства устойчивости по первому приближению [25].

Подробно классические способы построения функций Ляпунова приведены, например, в [58,59],

Второй метод Ляпунова широко используется в теории управления и реагирования. Для анализа абсолютной устойчивости нелинейных систем автоматического регулирования используются функции Ляпунова-Лурье, в задачах синтеза оптимальных управлений - функции Ляпунона-Ііеллмана, для оптимизации динамических систем -функции Ляпунопа-Лагранжа [17]. Вообще в настоящее время разработано большое число методов синтеза нелинейных регуляторов и анализа нелинейных систем автоматического управления (методы абсолютной устойчивости, гармонической линеаризации, оптимального управления и др.), которые, как правило, развиты на основе метода функций

21 Ляпунова [60]. Это лишь некоторые приложения прямого метода Ляпунова,

все приложения и направления развития которого в различных областях современной науки не представляется возможным охватить в одном обзоре.

Ряд технических аспектов применения прямого метода Ляпунова освещен в [60-62].

Однако применение рассматриваемого метода Ляпунова, основанное на использовании функций V{t,x), имеет существенную трудность, связанную с тем, что в настоящее время не известен общий способ построения функций Ляпунова.

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

Пусть рассматриваются уравнения возмущенного движения автономной системы

fix JL*.

at где Xі11 (х)^ X а Я' х"1 ..-.v^ , Y,Xm(x) -голоморфныефункции.

Ляпуновым [28] получены следующие теоремы об устойчивости по первому приближению системы (8)-

Пусть все собственные значения матрицы А линейною приближения для системы (8) имеют отрицательные вещественные части, т.е. Re?, (А)<0

( /=1,2,.,,,/7). Тогда невозмущеннос решение х = 0 системы (8) асимптотически устойчиво. Если среди собственных значений матрицы А найдется хотя бы одно с положительно вещественной частью (Ке1ДЛ)>0), то решение

л- = 0 системы (8) неустойчиво. В случае, когда для всех / = 1,2,,.. ,л Re?.;(/1)^0 и существуют ?,;, для которых Re?V;(/l) = Q, то невозмущенное

решение системы (8) может быть как устойчивым, так и неустойчивым. Это зависит от вида нелинейных функций A"*(.v) [17].

Последний случай, когда имеются нулевые вещественные части собст-

22 венных значений матрицы линейного приближения, Ляпуновым был назван

критическим (сомнительным). Им исследованы критические случаи [28];

l)?.I=0,ReX/U)<0(y = 2,...>ii);

2) ?V] =/0), ?v: =-/, Re?v;(/4)<0 ( у = 3,..,,я);

3)?.l=0,?^=0sRe?./(^)<0(y = 31„.fii).

Для неавтономных линейных систем Ляпунов развил теорию характеристических показателей решений, а для нелинейных систем с голоморфными

правыми частями

^A(t)x+ix"(ttx), (9)

at т-2

Ляпуновым даны критерии устойчивости и неустойчивости в случае систем с периодическими правыми частями и в случае правильных систем [17].

Приводимое ниже утверждение (признак Ляпунова) в существенных

чертах заимствовано из [22], Если система первого приближения — = Л(г )v

правильная и все ее характеристические показатели отрицательны, а нелиней-ные члены удовлетворяют условию ' ^X^(t,x)-J\\"\ т>\, то невозму-

іценное решение ,v = 0 полной нелинейной системы (9) асимптотически устойчиво.

При наличии положительного характеристического показателя решение .v = 0 системы (9) неустойчиво.

Для периодических систем исследованы и критические случаи [17]. Развитию теории критических случаев посвящены работы [63 -65],

Тенденции современных исследований устойчивости. Отметим, что во всех актуальных приложениях качественной теории дифференциальных уравнений по-прежнему анализ устойчивости опирается на оба метода Ляпунова. При этом сама качественная теория активно развивается в абстрактном математическом направлении. В частности, согласно [17], состояние вопроса характеризуется следующими результатами: совмещение идеи нескольких функций Ляпунова и метода сравнения с использованием векторных диффе-

23 ренциальньгх неравенств типа Чаплыгина-Важевского позволило получить первые теоремы об устойчивости с вектор-функцией Ляпунова (ВФЛ). Полнее были предложены и матричные функции Ляпунова,

В методе ВФЛ наряду с исходной системой (7) вводится вспомогательная система, называемая системой сравнения (СС) и описываемая конечномерным дифференциальным уравнением размерности т<п, возникающая из мажорирования производной ВФЛ в силу (7). С использованием аналогов условий теорем Ляпунова было доказано, что устойчивость или асимптотическая устойчивость СС в некотором смысле влечет за собой устойчивость, соответственно, асимптотическую устойчивость исходной системы (такие теоремы называются теоремами сравнения с ВФЛ) [17].

Использование теорем сравнения сводит задачу анализа устойчивости различною типа к построению ВФЛ, СС и существенно более простому анализу соответствующих свойств СС [17J.

Как отмечается в [17], прикладная значимость метода ВФЛ вскрыта в работе [66], в которой была предложена идея исследования устойчивости сложных нелинейных систем на основе их декомпозиции и последующего оценочного агрегирования с применением ВФЛ. Позднее это направление выросло в теорию устойчивости сложных систем. Ее изложение дано, например, в [67], где также указаны многочисленные приложения к системам разной природы.

В настоящее время доказаны сотни теорем сравнения с ВФЛ для различных динамических свойств нелинейных ОДУ. Результаты этих теорем объединены в виде основной идеи принципа сравнения: если существуют ВФЛ, удовлетворяющие подходящим условиям, то различные динамические свойства исходного уравнения вытекают из соответствующих динамических свойств системы сравнения [17]. Впервые в такой форме принцип сравнения дан в [68].

Систематическое изложение метода векторных функций Ляпунова в теории з'стойчиности с рассмотрением способов построения ВФЛ и СС, с различными количественными оценками для нелинейных систем и приложения-

24 ми ВФЛ дано в [67].

С друюй стороны, предприняты попытки компьютеризации анализа устойчивости [18-21].

Общая теория Ляпунова, помимо своего основною направления развития, нашла применение в новых отраслях современной математики - теориях особенностей, бифуркаций и катастроф. Бифуркация - качественная перестройка или метаморфоза объекта произвольной природы при изменении параметров, 01 которых он зависит. Катастрофа - скачкообразное изменение, возникающее r виде внезапного ответа системы на плавное изменение внешних условий [69]. Если в качестве системы рассмотреть систему ОДУ, то термин катастрофа будет описывать качественно свойство неустойчивости в смысле Ляпунова решения системы ОДУ.

Необходимо принять во внимание глобальное развитие и фундаментальное значение общей теории Ляпунова, которая нашла свое новое начало в современном интегральном направлении науки - синергетике. Синергетика исследует процессы самоорганизации и охватывает все отрасли знаний о косной и живой природе, технике и экономике. Смысл и содержание этой науки состоит в том, что в открытых системах любой природы возникают процессы самоорганизации, т.е. процессы рождения из биологического, экономического или социальною хаоса некоторых устойчивых упорядоченных структур с новыми свойствами систем. Поведение нелинейных систем вдали от положения равновесия при изменении некоторых управляющих параметров, пространства состояний систем, аттракторы (притягивающие множества в пространстве состояний, т.е. асимптотически устойчивые множества), «странные аттракторы» (аттракторы, отличные от состояний равновесия и периодических колебаний) - лишь часть понятий и вопросов этой науки, к которой наиболее близка по своей идеологии прикладная теория управления [70]. Последняя существенно базируется на общей теории устойчивости Ляпунова.

В 30-х годах прошлого века в выдвинутом Н.Г. Четаевым «постулате устойчивости», указывается, что «устойчивость - явление принципиально общее - как-то должна, по-видимому, проявляться в основных законах при-

25 роды» [70].

«Постулат устойчивости» устанавливал в законах природы теоретические свойства устойчивости как необходимое требование малых отклонений теории от эксперимента по наблюдаемым функциям [17].

Согласно этому постулату, только устойчивые движения являются физически реализуемыми. Однако открытие «странных» (хаотических) аттракторов потребовало корректировки постулата устойчивости. Оказалось, что в природе могут реализовываться и необычные движения на хаотическом аттракторе. В этой связи постулат устойчивости справедлив, в первую очередь, для определенной группы природных процессов, которым всегда присущи порядок и асимптотическая устойчивость [70].

Концепция синергетики в теории управления и ее связь с общей теорией устойчивости Ляпунова детально освещены в [70].

На основе изложенною тема устойчивости по Ляпунову актуальна для современной математической теории и ряда важных научных и технических приложений.

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

Таким образом, тема диссертационной работы актуальна.

Основной целью диссертации является разработка и исследование метода компьютерного анализа устойчивости по Ляпунову решений систем нелинейных ОДУ в нормальной форме.

Метод должен обладать вычислительной устойчивостью и универсальностью в той мерс, которая позволяет его рассматривать как основу для создания компьютерной технологии анализа устойчивости систем нелинейных ОДУ.

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

26 часть системы, начальные условия, шаг и промежуток решения, по и получать

само приближенное решение, попутно моделируя накопление погрешности.

Поскольку предлагаемый подход к компьютеризации анализа устойчивости основан на разностных приближениях решений ОДУ, он по построению отличается от методов символьной обработки и схем анализа устойчивости путем вычисления функций Ляпунова, Предложенные в дальнейшем критерии устойчивости по конструкции не опираются на методы качественной теории дифференциальных уравнений, однако методы этой теории используются для обоснования и исследования разрабатываемых в диссертации методов.

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

Для достижения поставленной цели в диссертационной работе решаются следующие задачи:

  1. выполнить синтез и обоснование программируемых условий устойчивости решений систем нелинейных ОДУ в нормальной форме на основе разностных схем приближенного решения;

  2. исследовать зависимость достоверности условий от замены их предельных выражений на конечные приближения при ограничениях сравнительно общего вида;

  3. исследовать разновидности условий и их отличительные особенности при использовании для их построения различных разностных схем - методов Эйлера» Эйлера-Коши, Рунге-Кутта и Адамса;

4) показать, что в условиях устойчивости существенно ограничивается рост

погрешности разностных схем, и частности, методы Эйлера и Эйлера-Коши имеют не более чем линейное накопление погрешности на произвольном промежутке с коэффициентами, не зависящими от длины промежутка;

  1. сконструировать и отладить программные модели анализа устойчивосги решений систем нелинейных ОДУ, инвариантные относительно правых частей систем нелинейных ОДУ, относительно размерности системы, возмущений начальных данных, разностных схем и выбора численных параметров этих схем, в числе которых длина промежутка и шаг решения;

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

  3. установить границы числовых параметров достоверной работы условий устойчивости, в частности, исследовать предельно допустимые границы шага и промежутка решения, их «оптимальные» значения, временные затраты при условии достоверности программного анализа устойчивости;

  4. исследовать зависимость результатов работы условий устойчивости от величины возмущения начальных данных и от погрешности используемой разностной схемы;

  5. оценить вычислительную устойчивость и универсальность конструируемых программных моделей и предлагаемых схем как основу компьютерной технологии анализа устойчивости решений систем нелинейных ОДУ.

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

Шучная новизна диссертационной работы заключается в разработке программных схем анализа устойчивости по Ляпунову решений систем нелинейных ОДУ в нормальной форме на основе разностных методов. Предлагаемые схемы отличаются от методов качественной теории дифференциальный

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

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

Конкретно, научная новизна результатов может быть охарактеризована следующим образом:

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

  2. предложены и обоснованы необходимые и достаточные условия устойчивости на основе преобразований разностных методов в форму бесконечных произведений; условия отличаются от известных по построению и инвариантностью относительно вида правой части системі.! ОДУ;

  3. исследована достоверность условий устойчивости при приближении бесконечных произведений конечным числом сомножителей при ограничениях сравнительно общего вида;

  4. в случае систем линейных ОДУ с постоянными коэффициентами предложенный программный анализ устойчивости, в отличие от известных методов, не требует построения характеристического мноючлена матрицы постоянных коэффициентов и информации о его корнях;

5) показано, что условия устойчивости решений существенно ограничивают

рост погрешности разностной схемы; в частности, методы Эйлера и Эйле-ра-Коши в рассматриваемых условиях имеют не более чем линейное накопление погрешности на произвольном промежутке с коэффициентами, не зависящими от длины промежутка, что отличает данные оценки от известных;

  1. сочетание данного качества с единообразной гїрограммируемостью предложенных условий устойчивости дает практический достоверный способ компьютерного анализа устойчивости решений нелинейных ОДУ общею вила;

  2. построены программные модели анализа устойчивости решений систем нелинейных ОДУ, инвариантные относительно правых частей ОДУ, размерности системы, возмущений начальных данных, разностных схем и выбора численных параметров этих схем; апробация моделей экспериментально подтверждает достоверность предложенных условий устойчивости;

  3. проведенные на базе данных схем моделирования эксперименты и предшествующее им формальное обоснование условий устойчивости моїут служить основой для разработки компьютерной технологии анализа устойчивости решений систем нелинейных ОДУ; помимо построения, это отличает предложенные схемы результатом практической реализации.

В приложении к решению технической задачи представлены следующие результаты:

  1. предложенные в диссертационной работе схемы программного анализа устойчивости решений систем нелинейных ОДУ апробированы на задаче управления вращением спутника (модель автопилота); получено совпадение результатов анализа с известными теоретическими оценками, выполненными на основе линеаризации; кроме того, известные оценки дополнены анализом в критическом по Ляпунову случае;

  2. в отличие от известного теоретического метода, моделирование по предложенной схеме выполняется непосредственно по виду исходной системы, без применения вспомогательных преобразований, включая линеаризацию;

  1. идентифицирована область устойчивости системы управления вращением спутника, результаты согласуются с теоретическими оценками области устойчивости, данными на основе линеаризации;

  2. на основе программной идентификации экстремумов функций четырех переменных при помощи сортировки оценивается максимальное и минимальное отклонение нелинейной системы от устойчивого состояния при вариации одновременно четырех параметров, включая коэффициенты системы;

  1. с помощью предложенных критериев дана схема компьютерной проверки на устойчивость решения системы при максимальном и минимальном отклонении от устойчивого состояния;

  2. сочетание оптимизационной схемы со схемой проверки на устойчивость программно промоделировано на примере системы Лоренца для случаев ее стабилизации, асимптотической устойчивости и неустойчивости.

Па основе данных схем возможно выполнение мониторинга системы на предмет устойчивости при вариации параметров системы. Основные положения, выносимые на защкгу:

  1. программируемые необходимые и достаточные условия устойчивости на основе мультипликативных преобразований разностных решений систем нелинейных ОДУ в нормальной форме;

  2. обоснование схем анализа устойчивости решений систем нелинейных ОДУ в виде частичных произведений на основе разностных методов Эйлсра-Коши, Рунге-Кутта н Адамса;

  3. линейное накопление погрешности разностных схем в условиях устойчивости с коэффициентами, не зависящими от длины промежутка численною интегрирования;

  4. программные модели для анализа устойчивости решений систем нелинейных ОДУ5 инвариантные относительно правых частей ОДУ, размерности системы, возмущений начальных данных, разностных схем и выбора численных параметров этих схем.

Практическая ценное і ь диссертационного исследования заключается в прикладном характере и компьютерной реализуемости предложенных условий

устойчивости- Исследование устойчивости по Ляпунову актуально для математической теории и практики в ряде важных научно-технических приложений. Предложенный способ позволяете помощью компьютера получить разностное решение, достоверную информацию о характере устойчивости решения ОДУ и одновременно моделировать накопление погрешности разностной схемы. Рассматриваемый в качестве основы компьютерной технолої и и анализа устойчивости, способ может служить эффективным инструментом проверки построения устойчивою управления для системы автоматическою регулирования.

Внедрение и использование результатов работы- Полученные в работе результаты использованы в НИИ механики и прикладной математики им. Воронина И.И. Ростовскою государственною университета; в госбюджетной НИР «Математические методы устойчивой параллельной обработки, поиска и распознавания», код ГРІГГИ 28.23.15, регистрационным номер 01.2.00106436; в учебном процессе кафедры информатики Таганрогского государственною педагогического института в курсах «Численные методы», «Компьютерное моделирование», курсах по выбору и практикуме решения задач на ЭВМ, что подтверждено соответствующими актами об использовании, приведенными в приложении 7 к диссертационной работе.

Апробация работы. Основные результаты работы докладывались на:

ІХ-ХІ Международных конференциях «Математические модели физических процессов и их свойства» (Таганрог, ТГПИ, 2003 -2005гг.);

1-ой Международной научно-практической конференции «Текст в системе высшего профессионального образования» (Таганрог, ТГПИ, 2003к);

Всероссийской научной конференции молодых ученых и аспирантов «Информационные технологии, системный анализ и управление» (Таганрог, ТРТУ, 2003і\);

V Научно-практической конференции преподавателей, студентов, аспирантов и молодых ученых ТИУиЭ (Таганрог, ТИУиЭ, 2004г.);

III Межрегиональной научно-практической конференции студентов, аспирантов и молодых ученых «Молодежь XXI века - будущее российской нау-

32 ки» (Ростов, РГУ, 2005г.);

научно-технических конференциях профессорско-преподавательскою составам аспирантов ТПІИ (Таганрог, 2000-2005гг.);

семинарах «Теоретическая и прикладная информатика» кафедры информатики ТГПИ (Таганрог, 2000 - 2005гг.)-

Публикации. По материалам работы опубликовано 12 печатных рабо і с общим объемом примерно 12 печатных листов.

Структура и объем работы. Диссертационная работа состоит из введения, 3 глав основного раздела, заюпочения, списка литературы и 7 приложений. Основное содержание работы - примерно 7 печатных листов, включая 3 таблицы и список литературы из 87 наименований.

Подобные работы
Буланов Сергей Георгиевич
Разработка и исследование методов программного моделирования устойчивости систем линейных дифференциальных уравнений на основе матричных мультипликативных преобразований разностных схем
Корнеев Петр Кириллович
Численные методы приближения функций и решения уравнений на основе непрерывных дробей
Русин Евгений Владимирович
Разработка и исследование методов и программных комплексов параллельной обработки изображений на основе вертикального представления данных
Кузнецова Юлия Геннадьевна
Разработка методов исследования функционально-технологических свойств пищевых рецептурных смесей на основе теории нечетких множеств
Саядян Дмитрий Левонович
Математическое моделирование стационарных магнитных полей на основе метода интегральных уравнений
Старовойтов Александр Степанович
Компьютерное моделирование взаимодействия заряженных пучков с плазмой на основе самосогласованной системы интегро-дифференциальных уравнений методом независимых частиц
Калиниченко Борис Олегович
Разработка методов выполнения функционально-параллельных программ на основе сетей Петри
Морозов Андрей Владимирович
Разработка фреймово-продукционной модели синтеза цифровых автоматов на основе метода спецификации состояний и ее программная реализация средствами реляционной СУБД
Головин Игорь Михайлович
Разработка алгоритмов, численных методов и программной среды для управления качеством рецептурных смесей на основе методов математического программирования
Авдеенко Татьяна Владимировна
Разработка методов исследования структурной идентифицируемости моделей в пространстве состояний

© Научная электронная библиотека «Веда», 2003-2013.
info@lib.ua-ru.net