САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
КАФЕДРА ТЕОРИИ УПРАВЛЕНИЯ
Гайфулин Роман Вячеславович
Выпускная квалификационная работа бакалавра
Построение области асимптотической
устойчивости для одного
дифференциально-разностного уравнения
Направление 010300
Фундаментальная информатика и информационные технологии
Научный руководитель,
доктор физ.-мат. наук
профессор
Жабко А. П.
Санкт-Петербург
2016
Оглавление
Введение ................................................................................................................... 3
Постановка задачи................................................................................................... 7
Обзор литературы.................................................................................................... 8
Глава 1. Основные теоретические сведения......................................................... 9
1.1 Понятие дифференциального уравнения..................................................... 9
1.2 Дифференциально-разностные уравнения ................................................ 10
1.3. Основные понятия теории устойчивости ................................................. 12
1.4. Второй метод Ляпунова ............................................................................. 12
1.5 Применение второго метода Ляпунова к исследованию на устойчивость
дифференциально-разностных уравнений ...................................................... 16
Глава 2. Построение области асимптотической устойчивости ........................ 19
2.1 Постановка начальной задачи и определение устойчивости .................. 19
2.2 Исследование на устойчивость по первому приближению ..................... 22
2.3 Методы нахождения областей асимптотической устойчивости
дифференциально-разностных уравнений, стационарных в первом
приближении ...................................................................................................... 24
Глава 3. Программная реализация ....................................................................... 27
3.1 Метод Адамса численного решения дифференциально-разностного
уравнения ............................................................................................................ 27
3.2 Пример построения области асимптотической устойчивости ................ 29
Выводы ................................................................................................................... 35
Заключение ............................................................................................................ 36
Список литературы ............................................................................................... 37
2
Введение
Решения многих дифференциальных уравнений и их систем не
выражаются в элементарных функциях, поэтому для отыскания их решения
прибегают к приближенным методам интегрирования. Однако чаще всего
существует необходимость не узнать конкретное решение, а исследовать его
свойства. В качестве таких свойств могут выступать: изменение поведения
решения при изменении некоторых параметров, изменение некоторых
решений при незначительных изменениях начальных данных, характер
решения, изменение поведения целой системы при изменении отдельных
параметров. Изучением этих вопросов занимается качественная теория
дифференциальных уравнений.
Одним из важнейших вопросов этой теории является исследование на
устойчивость дифференциальных уравнений и их систем. Важность данного
вопроса обусловлена физическим смыслом дифференциального уравнения:
это
модель
движения
физического
объекта.
При
исследовании
дифференциального уравнения на устойчивость главным вопросом является
то, как поведут себя отдельные решения при незначительном изменении
начальных условий: будут ли малые изменения начальных условий означать
малое изменение всего решения.
Для постановки задачи об устойчивости, необходимы следующие
данные:
1) объект, об устойчивости которого будет идти речь,
2) определение устойчивости.
Большинство реальных динамических объектов, чаще всего это
механические
системы,
могут
быть
описаны
дифференциальными
уравнениями или их системами. Устойчивость таких систем чаще всего
рассматривается в смысле устойчивости по Ляпунову. Это понятие,
3
реализующее идею «малых» отклонений решения дифференциального
уравнения на промежутке времени [0, +∞) при «небольших» вариациях
начальных данных этого решения.
Теория устойчивости занимается разработкой методов, которые бы
позволяли судить об устойчивости некоторого заданного решения, не зная при
этом общего решения дифференциального уравнения. Основоположником
этой теории является российский математик А.М. Ляпунов.
Он разделил все многообразие методов теории устойчивости на два
класса. В первый класс были включены те методы исследования
дифференциальных
уравнений
на
устойчивость,
которым
требуется
определенная информация о решениях данной системы. Этот подход получил
название первого метода Ляпунова. Второй же метод Ляпунова составляет
совокупность приемов и средств, которые позволяют исследовать решения
дифференциальных уравнений и их систем при помощи специальных функций
Ляпунова.
Частное решение, которое исследуется на устойчивость, часто называют
невозмущенным, а любое другое – возмущенным решением.
Устойчивость бывает двух типов:
1) Устойчивость относительно возмущения начальных данных.
Основная
идея
заключается
в
том,
что
решения
𝑥 =
𝑥(𝑡; 𝑡0 , 𝑥0 ) и 𝑥 = 𝑥(𝑡; 𝑡0 , 𝑥1 ) с начальными данными (𝑡0 , 𝑥0 ) и
(𝑡0 , 𝑥1 ) должны быть близкими для любого 𝑡 ≥ 𝑡0 , если 𝑥0 и 𝑥1
достаточно близки.
Важнейшим требованием в этом случае является требование
близости решений при всех значениях 𝑡 ≥ 𝑡0 , так как если
ограничить рассмотрение конечным интервалом изменения времени
𝑡, то указанное свойство сохранения близости будет прямо следовать
4
из классической теоремы о непрерывной зависимости решений
системы дифференциальных уравнений от начальных данных.
2) Устойчивость относительно постоянно действующих внешних
возмущений.
Для конечного временного интервала решения исходной 𝑥 = 𝑓(𝑡, 𝑥)
и возмущенной 𝑥 = 𝑓(𝑡, 𝑥) + 𝑔(𝑡, 𝑥) систем, выходящие из одной и
той же начальной точки, будут близкими на данном конечном
интервале, при условии, что внешние возмущения 𝑔(𝑡, 𝑥) во второй
системе достаточно малы.
Это – прямое следствие классической теоремы о непрерывной
зависимости решений системы дифференциальных уравнений от
правых частей уравнений. Поэтому, когда говорят об устойчивости
при постоянно действующих возмущениях, подразумевают, что
𝑡 ≥ 𝑡0 .
Понятие устойчивости по Ляпунову можно расширить. Решение
называется асимптотически устойчивым, если оно устойчиво по Ляпунову и
при малых отклонениях начальных данных отклонения решения на
промежутке
времени
[0, +∞)
стремятся
к
нулю.
Асимптотическая
устойчивость гарантирует, что при небольшом изменении начальных данных
с течением времени отклонения решения системы станут незначительны.
Большинство процессов, заключающихся в передаче массы, энергии или
информации характеризуются тем, что в них присутствует запаздывание.
Существует множество обуславливающих его причин — например,
ограничение
скорости
распространения
сигнала
или
инерционность
некоторых элементов. Если описывать их без использования параметра,
учитывающего запаздывание, результаты могут оказаться неточными или
вовсе абсурдными. Поэтому целесообразно рассматривать дифференциальные
уравнения, принимающие различные значения аргумента.
5
Одним
из
возможных
вариантов
является
рассмотрение
дифференциально-разностных уравнений вида 𝑥 ′ (𝑡) = 𝑓[𝑡, 𝑥(𝑡), 𝑥(𝑡 – ℎ)], в
котором 𝑥 — неизвестная функция независимого аргумента 𝑡, 𝑓: 𝑅3 → 𝑅, а
ℎ — положительное число – запаздывание. Такие уравнения называются
дифференциально-разностными уравнениями с запаздыванием.
В данной работе решается задача аналитического построения области
асимптотической устойчивости для одного такого дифференциальноразностного уравнения и анализа устойчивости его решения при изменении
начальных параметров с помощью среды MATLAB.
6
Постановка задачи
Целью данной работы является анализ аналитически построенной
области
асимптотической
устойчивости
одного
дифференциально-
разностного уравнения в среде MATLAB. Теоретическим обоснованием для
проверки на асимптотическую устойчивость является модифицированный
метод
Адамса.
В
качестве
решения
задачи
построения
области
асимптотической устойчивости будет использован амплитудно-фазовый
метод.
Таким
образом,
для
достижения
поставленной
цели
можно
сформулировать ряд задач:
Найти решение дифференциально-разностного уравнения при наперед
заданном начальном условии с помощью модифицированного метода
Адамса;
Исследовать найденное решение на асимптотическую устойчивость;
В случае, если оно устойчиво – аналитически построить область
асимптотической устойчивости и провести анализ при разных
отклонениях параметров.
Результатом работы должна стать программа, которая строит численное
решение
дифференциально-разностного
уравнения
с
запаздывающим
аргументом и позволяет пользователю изменять входные изменять входные
параметры уравнения (коэффициенты и запаздывание) для последующего
анализа области асимптотической устойчивости.
7
Обзор литературы
В работе «Дифференциально-разностные уравнения» Беллмана Р. и
Кука К. Л. подробно излагается теория линейных дифференциальноразностных уравнений, теория устойчивости и асимптотическое поведение
решений.
Учебное пособие «Теория устойчивости движения», написанное
Ногиным В. Д. излагает базовые понятия и методы теории устойчивости, в том
числе первый и второй методы Ляпунова, последний из которых будет
рассмотрен в данной работе.
Книга
«Устойчивость
решений
дифференциально-разностных
уравнений» за авторством Л. Э. Эльсгольца
предоставляет обширный
материал по теории дифференциально-разностных уравнений, а так же
способы их решения и исследования на асимптотическую устойчивость.
В
статье
«Построение
области
асимптотической
устойчивости
дифференциально – разностных уравнений в среде MATLAB», написанной
Жабко А. П. и Зараник У. П., был выведен алгоритм метода Адамса для
дифференциально-разностных
модифицирован
в
данной
уравнений,
работе
для
который
впоследствии
дифференциально-разностного
уравнения с запаздыванием и реализован в среде MATLAB.
8
был
Глава 1. Основные теоретические сведения
1.1 Понятие дифференциального уравнения
В работе [2] приводятся общие определения, процитированные ниже.
Дифференциальное уравнение – это уравнение, в которое входит
неизвестная функция под знаком производной или дифференциала.
Если неизвестная функция является функцией одной переменной, то
дифференциальное уравнение называют обыкновенным (ОДУ).
Максимальный порядок производной неизвестной функции, входящей в
дифференциальное уравнение, называется порядком дифференциального
уравнения.
Приведем еще несколько определений из работ [1] и [2]. Обыкновенным
дифференциальным уравнением n-го порядка называется соотношение вида
𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , … , 𝑦 (𝑛) ) = 0
между независимым переменным x, его функцией y и производными
y’,y’’,…,y(n).
Функция 𝑦 = 𝜑(𝑥) называется решением этого дифференциального
уравнения, если после замены y на 𝜑(𝑥), y’ на 𝜑 ’(x),…,y(n) на 𝜑 (n)(x) оно
обращается в тождество.
Общее решение дифференциального уравнения – это множество
решений, содержащее все без исключения решения этого дифференциального
уравнения.
Если же решение дифференциального уравнения удовлетворяет
изначально заданным дополнительным условиям, то его называют частным
решением дифференциального уравнения.
9
Основными задачами теории дифференциальных уравнений являются
задачи Коши, краевые задачи и задачи нахождения общего решения
дифференциального уравнения на каком-либо заданном интервале 𝑋.
Задача
Коши
–
это
задача
нахождения
частного
решения
дифференциального уравнения, удовлетворяющего заданным начальным
условиям
𝑓(𝑥0 ) = 𝑓0 ; 𝑓 ′ (𝑥0 ) = 𝑓1 ; 𝑓 ′′ (𝑥0 ) = 𝑓2 ; … ; 𝑓 (𝑛−1) (𝑥0 ) = 𝑓𝑛−1 ,
где 𝑓0 ; 𝑓1 ; 𝑓2 ; … ; 𝑓𝑛−1 – числа.
Характеристическое
уравнение
линейного
однородного
дифференциального уравнения 𝑛-ой степени с постоянными коэффициентами
– это уравнение 𝑛-ой степени вида
𝑓𝑛 ∗ 𝑘 𝑛 + 𝑓𝑛−1 ∗ 𝑘 𝑛−1 + ⋯ + 𝑓1 ∗ 𝑘 + 𝑓0 = 0.
Стоит упомянуть, что на данный момент не существует единого подхода
для решения дифференциальных уравнений, в силу многообразия их видов.
Отдельного рассмотрения требует такой класс дифференциальных
уравнений как дифференциально-разностные уравнения.
1.2 Дифференциально-разностные уравнения
Теория изложена в работе [3].
Под дифференциально-разностным уравнением мы будем в дальнейшем
понимать уравнение относительно неизвестной функции и её производных,
вычисленных при некоторых значениях аргумента, отличающихся на
постоянные. Все производные будем считать обыкновенными.
Так как наши уравнения будут включать в себя производные и разности
различных порядков, стоит дать пару определений.
10
Дифференциальный порядок – порядок наивысшей входящей в него
производной.
Разностный порядок – число, на единицу меньше числа различных
значений аргумента, встречающихся в уравнении.
Общее дифференциально-разностное уравнение дифференциального
порядка n и разностного порядка m имеет вид:
𝐹[𝑡, 𝑢(𝑡), 𝑢( 𝑡 − 𝜔1 ), … , 𝑢 (𝑡 − 𝜔𝑚 ), 𝑢′ (𝑡), 𝑢′ (𝑡 − 𝜔1 ) , … , 𝑢(𝑛) (𝑡), 𝑢𝑛 (𝑡
− 𝜔1 ), … , 𝑢𝑛 (𝑡 − 𝜔𝑚 )] = 0,
где F – заданная функция от 1 + (𝑚 + 1)(𝑛 + 1) переменных; числа
𝜔1 , 𝜔2 , … , 𝜔𝑚 , называемые запаздываниями, также заданы.
Даже простейшее уравнение вида
𝑎0 𝑢′ (𝑡) + 𝑎1 𝑢′ (𝑡 − 𝜔) + 𝑏0 𝑢(𝑡) + 𝑏1 𝑢(𝑡 − 𝜔) = 𝑓(𝑡)
может приводить к различным заключениям при различных значениях
коэффициентов. Поэтому выделим несколько видов дифференциальноразностных уравнений:
Уравнения запаздывающего типа
𝑎0 𝑢′ (𝑡) + 𝑎1 𝑢′ (𝑡 − 𝜔) + 𝑏0 𝑢(𝑡) + 𝑏1 𝑢(𝑡 − 𝜔) = 𝑓(𝑡), если a0≠0 и a1=0,
Уравнения нейтрального типа
𝑎0 𝑢′ (𝑡) + 𝑎1 𝑢′ (𝑡 − 𝜔) + 𝑏0 𝑢(𝑡) + 𝑏1 𝑢(𝑡 − 𝜔) = 𝑓(𝑡), если a0≠0 и a1≠0,
Уравнения опережающего типа
𝑎0 𝑢′ (𝑡) + 𝑎1 𝑢′ (𝑡 − 𝜔) + 𝑏0 𝑢(𝑡) + 𝑏1 𝑢(𝑡 − 𝜔) = 𝑓(𝑡), если a0=0, a1≠0 и
b0≠0.
Если a0=a1=0, то уравнения является чисто разностным, и в этом случае
оно принадлежит к весьма детально исследованному типу функциональных
уравнений. Если 𝑏0 = 𝑏1 = 0, то оно сводится к чисто разностному уравнения.
Если 𝑎0 = 𝑏0 =0 или 𝑎1 = 𝑏1 =0, то оно является ОДУ.
11
1.3. Основные понятия теории устойчивости
Предположим,
что
некоторое
явление
описывается
системой n дифференциальных уравнений
𝑑𝑥𝑖
= 𝑓𝑖 (𝑡, 𝑥1 , 𝑥2 , … , 𝑥𝑛 ),
𝑑𝑡
с начальными условиями
𝑥𝑖 (𝑡0 ) = 𝑥𝑖0 ,
Будем
считать,
что
𝑖 = 1,2, … , 𝑛
𝑖 = 1,2, … , 𝑛.
функции
𝑓𝑖 (𝑡, 𝑥1 , 𝑥2 , … , 𝑥𝑛 )
определены
и
непрерывны вместе со своими частными производными на множестве
{𝑡 ∈ [𝑡0 , +∞), 𝑥𝑖 ∈ 𝑅𝑛 }. Далее без ограничения общности полагаем, что
начальный момент равен нулю:
𝑡0 = 0.
Систему дифференциальных уравнений удобнее записать в векторной форме:
𝑋′ = 𝑓(𝑡, 𝑋), где𝑋 = (𝑥1 , 𝑥2 , … , 𝑥𝑛 ), 𝑓 = (𝑓1 , 𝑓2 , … , 𝑓𝑛 ).
Если траектория движения системы мало изменяется при малых возмущениях
начального положения, то говорят, что движение системы является
устойчивым. [4]
1.4. Второй метод Ляпунова
Второй метод Ляпунова [5] включает в себя методы исследования
решений систем дифференциальных уравнений на устойчивость.
1.4.1 Некоторые предварительные определения.
Рассмотрим вещественную непрерывную числовую функцию
𝑉 = 𝑉(𝑡, 𝑥 ) ∈ 𝐶𝑡𝑥 (𝑍),
где 𝑍 = {(𝑡, 𝑥) | 𝑡 ≥ 0, || 𝑥 || ≤ 𝐻 } при некотором 𝐻 > 0.
12
Определение
1.
Функция
𝑉(𝑡, 𝑥)
называется
положительно
(отрицательно) определенной в Z, если существует непрерывная числовая
функция 𝑊(𝑥) ∈ 𝐶(||𝑥|| ≤ 𝐻), такая, что
𝑉(𝑡, 𝑥 ) ≥ 𝑊(𝑥) > 0 ∀(𝑡, 𝑥) ∈ 𝑍, 𝑥 ≠ 0,
(𝑉(𝑡, 𝑥) ≤ – 𝑊(𝑥) < 0 ∀(𝑡, 𝑥) ∈ 𝑍, 𝑥 ≠ 0)
𝑉(𝑡, 0) = 𝑊(0) = 0 ∀𝑡 ≥ 0.
Следует заметить, что в качестве 𝑊(𝑥) из приведённого определения
иногда удается использовать функцию вида 𝑊(𝑥) = inf|𝑉(𝑡, 𝑥)|.
𝑡≥0
Определение 2. Положительно или отрицательно определенную
функцию именуют знакоопределенной.
Определение 3. Функция 𝑉(𝑡, 𝑥) имеет (допускает) бесконечно малый
высший предел (б.м.в.п.), если 𝑉(𝑡, 0) = 0 при всех 𝑡 ≥ 0, причём эта функция непрерывна по х в точке х = 0 равномерно относительно 𝑡 ≥ 0, т.е.
(∀𝜀 > 0)(∃𝛿 = 𝛿(𝜀) > 0: |𝑉(𝑡, 𝒙)| < 𝜀
∀𝑡 ≥ 0, ∀𝒙: ‖𝒙‖ < 𝛿.
Рассмотрим систему дифференциальных уравнений
𝑑𝒙
𝑑𝑡
= 𝒇(𝑡, 𝒙) (𝑡 ≥ 0),
(1.1)
где
𝑥1 (𝑡)
𝒙 = 𝒙(𝑡) = ( ⋮ ) ,
𝑥𝑛 (𝑡)
𝑓1 (𝑡, 𝒙)
𝒇(𝑡, 𝒙) = ( ⋮ )
𝑓𝑛 (𝑡, 𝒙)
причем 𝒇(𝑡, 𝒙) непрерывна на 𝒁 = {(𝑡, 𝒙) | 𝑡 ≥ 0, || 𝒙 || ≤ 𝐻} при
некотором 𝐻 > 0. Будем считать, что система (3.1) имеет нулевое решение
𝜉 = 0, т.е. 𝑓(𝑡, 0) = 0 ∀𝑡 ≥ 0.
(1,1)
Определение 4. Пусть имеется функция 𝑉 = 𝑉(𝑡, 𝒙) ∈ 𝐶𝑡𝒙 (𝒁) .
Рассмотрим произвольную пару (𝑡, 𝒙) ∈ 𝒁 и соответствующее этой паре
13
начальных данных решение 𝒙 = 𝒙(𝜏; 𝑡, 𝒙) системы (3.1), так что 𝒙( 𝑡; 𝑡, 𝒙 ) =
𝒙. Производной по времени t функции 𝑉(𝑡, 𝒙) в силу системы (3.1) называют
функцию
𝑛
𝑑
𝜕𝑉
𝜕𝑉
𝜕𝑉
𝑉̇ (𝑡, 𝒙) = [ 𝑉(𝜏, 𝒙(𝜏; 𝑡, 𝒙))]|
=
+∑
∙ 𝑓𝑘 (𝑡, 𝒙) =
+
𝑑𝜏
𝜕𝑡
𝜕𝒙𝑘
𝜕𝑡
𝜏=𝑡
𝑘=1
+ 〈𝛻𝑉, 𝒇〉,
(1.2)
есть градиент функции 𝑉.
1.4.2 Устойчивость по Ляпунову
Решение 𝜑(𝑡) системы дифференциальных уравнений
𝑋′ = 𝑓(𝑡, 𝑋)
с начальными условиями
𝑋(0) = 𝑋0
Устойчиво по Ляпунову, если для любого 𝜀 > 0 найдется число
𝛿 = 𝛿(𝜀) > 0, такое, что если
∥ 𝑋(0) − 𝜑(0) ∥< 𝛿, то ∥ 𝑋(𝑡) − 𝜑(𝑡) ∥< 𝜀
для всех значений 𝑡 ≥ 0 . В противном случае решение 𝜑(𝑡) называется
неустойчивым.
Теорема (теорема Ляпунова об устойчивости)
Пусть дана система (1.1), имеющая нулевое решение. Если существует
(1,1)
положительно определённая функция 𝑉 = 𝑉(𝑡, 𝒙) ∈ 𝐶𝑡𝒙 (𝒁) , называемая
функцией Ляпунова, которая обладает неположительной производной по
времени 𝑉̇ (𝑡, 𝒙) в силу системы (1.1), то нулевое решение этой системы
устойчиво по Ляпунову.
14
1.4.3 Асимптотическая устойчивость
Если решение φ(𝑡) системы дифференциальных уравнений не только
устойчиво в смысле Ляпунова, но и удовлетворяет соотношению
𝑙𝑖𝑚𝑡 → ∞ ∥ 𝑋(𝑡) − 𝜑(𝑡) ∥= 0
при условии
∥ 𝑋(0) − 𝜑(0) ∥< 𝛿,
то говорят, что решение 𝜑(𝑡) является асимптотически. В этом случае
все решения, достаточно близкие к 𝜑(0 ) в начальный момент времени,
постепенно сходятся к 𝜑(𝑡) при увеличении 𝑡.
Графическое представление асимптотической устойчивости изображено
на рисунке 1.1.
Рис. 1.1 – Асимптотическая устойчивость
Теорема (теорема Ляпунова об асимптотической устойчивости)
Если для системы (1.1) существует положительно определенная
(1,1)
функция Ляпунова 𝑉 = 𝑉(𝑡, 𝒙) ∈ 𝐶𝑡𝒙 (𝒁), допускающая бесконечно малый
высший предел и имеющая отрицательно определенную производную
𝑉̇ (𝑡, 𝒙) по времени в силу системы (3.1), то нулевое решение 𝜉 = 0
рассматриваемой системы асимптотически устойчиво по Ляпунову.
15
1.4.4 Теорема о неустойчивости
Теорема (первая теорема Ляпунова о неустойчивости)
Пусть
для
системы
𝑑𝒙
𝑑𝑡
= 𝒇(𝑡, 𝒙) (𝑡 ≥ 0),
удовлетворяющей
предположениям из пункта 1.4.1, сделанным к системе 1.1 (в частности,
(1,1)
𝒇(𝑡, 𝟎) ≡ 𝟎 ), существует функция Ляпунова 𝑉 = 𝑉(𝑡, 𝒙) ∈ 𝐶𝑡𝒙 (𝒁) , такая,
что
1) 𝑉(𝑡, 𝒙) допускает б.м.в.п.;
2) существует число 𝑡̅ ≥ 0 такое, что для любого 𝛿 > 0 найдется точка ̅
𝒙,
̅‖ < 𝛿, 𝑉(𝑡̅, 𝒙
̅) > 0;
удовлетворяющая неравенствам ‖𝒙
3) производная 𝑉̇ (𝑡, 𝒙) по 𝑡 в силу системы 1.1 положительно определена.
Тогда нулевое решение системы 1.1 неустойчиво.
1.5 Применение второго метода Ляпунова к исследованию
на устойчивость дифференциально-разностных уравнений
Второй метод Ляпунова лишь с незначительными изменениями может
применяться
к
переформулируем
дифференциально-разностным
вышеописанные
теоремы
уравнениям.
для
Поэтому
дифференциально-
разностных уравнений.
Теорема (первая теорема второго метода Ляпунова)
Если для системы дифференциальных уравнений
𝑥𝑖̇ (𝑡) = 𝑓𝑖 (𝑡, 𝑥𝑗 (𝑡 − 𝜏𝑗𝑠 (𝑡))),
(1.3)
𝑖, 𝑗 = 1, 2, … , 𝑛; 𝑠 = 1, 2, … , 𝑚,
возмущенного
движения
можно
найти
знакоопределенную
функцию
𝑉(𝑡, 𝑥1 , 𝑥2 , … , 𝑥𝑛 ), для которой производная по времени вдоль интегральных
кривых
16
𝑛
𝑑𝑉
𝜕𝑉
𝜕𝑉
=
+ ∑
𝑓 = 𝛷 (𝑡, 𝑥𝑖 (𝑡), 𝑥𝑗 (𝑡 − 𝜏𝑗𝑠 (𝑡)))
𝑑𝑡
𝜕𝑡
𝜕𝑥𝑠 𝑠
𝑠=1
знакопостоянна, противоположного знака с функцией 𝑉 или тождественно
равна нулю, то невозмущенное движение устойчиво.
Теорема (вторая теорема Ляпунова)
Если в условиях предыдущей теоремы производная
функцией знакоопределенной (то есть
𝑑𝑉
𝑑𝑡
𝜕𝑉
𝜕𝑡
является
= 𝛷 (𝑡, 𝑥𝑖 (𝑡), 𝑥𝑗 (𝑡 − 𝜏𝑗𝑠 (𝑡))) ≤
−𝑊1 (𝑥𝑖 , 𝑥𝑗𝑠 ) , где 𝑊1 = 0 лишь при 𝑥𝑖 = 𝑥𝑗𝑠 = 0 (𝑖, 𝑗 = 1, 2, … , 𝑛; 𝑠 =
1, 2 , … , 𝑚 , причем если некоторые 𝜏𝑖𝑠 ≡ 0, то 𝑥𝑖 = 𝑥𝑖𝑠 ), а функция 𝑉
допускает бесконечно малый высший предел 𝑢 sup max 𝜏𝑗𝑠 (𝑡) < ∞ , то
𝑡≥𝑡0 𝑗, 𝑠
невозмущенное движение асимптотически устойчиво.
Теорема (теорема о неустойчивости)
Если для системы дифференциальных уравнений возмущенного
движения (1.3) существует дифференцируемая функция 𝑉(𝑡, 𝑥1 , 𝑥2 , … , 𝑥𝑛 ) ,
удовлетворяющая в некоторой ℎ-окрестности начала координат условиям:
1) 𝑉(𝑡, 0,0, … ,0) = 0;
2) при 𝑡 ≥ 𝑡0 в сколь угодно малой окрестности начала координат
существует независящая от 𝑡 область 𝑉 > 0;
3) в области 𝑉 > 0 функция 𝑉 ограничена;
4) в области 𝑉 > 0 производная
𝑛
𝑑𝑉
𝜕𝑉
𝜕𝑉
=
+ ∑
𝑓 (𝑡, 𝑥𝑖 (𝑡), 𝑥𝑗 (𝑡 − 𝜏𝑠 (𝑡))) = 𝛷 (𝑡, 𝑥𝑖 (𝑡), 𝑥𝑗 (𝑡 − 𝜏𝑠 (𝑡))) > 0
𝑑𝑡
𝜕𝑡
𝜕𝑥𝑖 𝑖
𝑖=1
(это значит, что 𝛷(𝑡, 𝑥𝑖 , 𝑥𝑗 ) > 0 , если 𝑡 ≥ 𝑡0 , а точки с координатами
(𝑥1 , 𝑥2 , … , 𝑥𝑛 ) и (𝑥1𝑠 , 𝑥2𝑠 , … , 𝑥𝑛𝑠 ) , 𝑠 = 1, 2, … , 𝑚 , лежат в области 𝑉 > 0 ),
17
причем при 𝑉 ≥ 𝑎 > 0 производная будет
𝑑𝑉
𝑑𝑡
≥ 𝑙 > 0, где 𝑙 – постоянная, - то
невозмущенное движение неустойчиво.
Можно доказать и аналоги других теорем о неустойчивости, но значение
второго метода Ляпунова для уравнений с запаздывающим аргументом,
очевидно, не очень велико, поэтому в данной работе они не приводятся.
Подытожив все сказанное выше, можно сказать, что большинство
основных теорем теории устойчивости могут быть перенесены на
дифференциально-разностные уравнения с незначительными изменениями.
Однако стоит заметить, что эффективное исследование устойчивости точек
покоя дифференциально-разностных уравнений сейчас осуществимо лишь для
тех уравнений, которые стационарны в первом приближении в некритических
случаях. [6]
Так
как
второй
метод
для
исследования
на
устойчивость
дифференциально-разностных уравнений недостаточно эффективен, для
определения
и
построения
области
асимптотической
устойчивости
дифференциально разностного уравнения нужно использовать новые методы.
Один из них мы разберем в этой работе далее.
18
Глава 2. Построение области асимптотической
устойчивости
2.1 Постановка начальной задачи и определение
устойчивости
Пусть имеется дифференциально-разностное уравнение с запаздывающим аргументом
𝑥̇ (𝑙) = 𝑓 (𝑓, 𝑥(𝑡), 𝑥(𝑡 − 𝑥(𝜏))) , 𝜏(𝑡) > 0.
Основной
нашей
начальной
задачей
является
(2.1)
определение
непрерывного решения 𝑥(𝑡 ) при 𝑡 ≥ 𝑡0 при условии, что 𝑥(𝑡) = 𝜑(𝑡) на
начальном множестве 𝐸𝑡0 состоящем из точки 𝑡 = 𝑡0 и значений 𝑡 − 𝜏(𝑡) ≤ 𝑡0
при 𝑡 ≥ 𝑡0 . Функция 𝜑(𝑡) называется начальной функцией на множестве 𝐸𝑡0 .
При постоянном запаздывании 𝜏(𝑡) = 𝜏 начальное множество 𝐸𝑡0 является
отрезком 𝑡0 − 𝜏 ≤ 𝑡 ≤ 𝑡0
Как известно, существует единственное решение этой начальной задачи,
если 𝑓, 𝜑 и 𝜏(𝑡) непрерывны и 𝑓 удовлетворяет условию Липшица по второму
и третьему аргументам. [6]
Условие Липшица
Если для любых точек х и х', принадлежащих отрезку [а, Ь], приращение
функции 𝑓 удовлетворяет неравенству
|𝑓(𝑥) − 𝑓(𝑥 ′ )| ≤ 𝑀|𝑥 − 𝑥|′𝛼 ,
Где 0 < 𝛼 ≤ 1 и M – некоторая постоянная, то, говорят, что функция
𝑓(𝑥) на отрезке [a,b] удовлетворяет условия Липшица порядка а, и пишут:
𝑓(𝑥) ∈ 𝐿𝑖𝑝 𝛼 ,
или 𝑓 ∈ 𝐿𝑖𝑝𝑀𝛼 , или 𝑓 ∈ 𝐻𝛼 (𝑀).
Каждая
функция,
удовлетворяющая при каком-либо 𝑎 > 0 условию Липшица. на отрезке [ а, 𝑏],
19
равномерно непрерывна на [ а, 𝑏] , a функции, удовлетворяющие условию
Липшица степени 𝑎 = 1 ,- абсолютно непрерывны. Функция, имеющая
на [а, 𝑏] ограниченную производную, удовлетворяет на [а, 𝑏] условию
Липшица с любым α ≤ 1. [7]
В дальнейшем в данной работе эти условия будем считать
выполненными.
Для уравнения 𝑛-го порядка с запаздывающим аргументом :
𝑥 𝑛 (𝑡) = 𝑓 (𝑡, 𝑥(𝑡), … , 𝑥 (𝑛−1) (𝑡), 𝑥(𝑡 − 𝜏(𝑡)), … , 𝑥 𝑛−1 (𝑡 − 𝜏(𝑡))),
где на 𝜏(𝑡) наложены прежние ограничения.
Основная начальная задача ставится почти так же:
𝑥(𝑡) = 𝜑(𝑡) и 𝑥 (𝑘) (𝑡) = 𝜑 𝑘 (𝑡), 𝑘 = 1,2, … , (𝑛 − 1),
на начальном множестве 𝐸𝑡0 , если определяется решение, имеющее
непрерывные производные до (𝑛 − 1)-го порядка включительно при 𝑡 ≥ 𝑡0 .
В тех процессах, которые описываются ОДУ без отклоняющихся
аргументов, или дифференциально-разностными уравнениями, возмущения
могут происходить в любой момент времени, а не только в какой-то
фиксированный момент времени 𝑡 = 𝑡0 . Но в определении устойчивости по
Ляпунову для дифференциальных уравнений без отклоняющихся аргументов,
которые удовлетворяют обычным условия существования единственности
непрерывной зависимости решения от начальных условий, достаточно
потребовать выполнения устойчивости только для возмущений в некоторый
фиксированный момент времени 𝑡 = 𝑡0 , поскльку теорема о непрерывной
зависимости решения от начальных значений и возможность продолжить
решение в сторону убывающих значений 𝑡 дают нам возможность утверждать,
что устойчивое по отношению к возмущению начальных значений решение
будет устойчиво и к возмущениям в любой другой момент времени 𝑡0̅ > 𝑡0 .
20
Вообще говоря, решения многих дифференциальных уравнений не
могут быть непрерывно продолжены в сторону убывающих значений 𝑡.
Поэтому устойчивые по отношению к возмущению начальной функции на
множестве 𝐸𝑡0 решения могут оказаться неустойчивыми по отношению к
возмущениям на начальном множестве 𝐸𝑡̅0 .
Поэтому,
целесообразно
изменить
определение
устойчивости
дифференциальных уравнений с запаздывающим аргументом вида (2.1).
Как уже говорилось ранее, функция 𝜑(𝑡) называется начальной
функцией на множестве 𝐸𝑡0 . Начальную функцию на множестве 𝐸𝑡̅0 мы будем
обозначать 𝜓(𝑡). Решение 𝑥(𝑡) на множестве 𝐸𝑡0 будем обозначать
решение 𝑥(𝑡) на множестве 𝐸𝑡̅0 будем обозначать
𝑥(𝑡)
𝜓(𝑡)
𝑥(𝑡)
𝜑(𝑡)
, а
.
Решение 𝑥(𝑡) , определяемое начальной функцией 𝜑(𝑡) называется
устойчивым, если для каждого 𝜀 > 0 существует такое 𝛿(𝑡̅0 , 𝜀) , что из
неравенства |𝜓(𝑡) −
неравенство |
𝑥(𝑡)
𝑥(𝑡)
−
𝜓(𝑡)
𝜑(𝑡)
| < 𝛿(𝑡̅0 , 𝜀) на начальном множестве 𝐸𝑡̅0 следует
𝑥(𝑡)
𝜑(𝑡)
| < 𝜀 при 𝑡 ≥ 𝑡̅0 для любого 𝑡0̅ > 𝑡0 .
Если, помимо этого, lim |
𝑡→∞
𝑥(𝑡)
−
𝜓(𝑡)
𝑥(𝑡)
| = 0 , то решение
𝜑(𝑡)
𝑥(𝑡)
𝜑(𝑡)
называется асимптотически устойчивым. Если же мы можем выбрать 𝛿 таким
образом, чтобы она зависела только от 𝜀, то устойчивость будет называться
равномерной.
Для дифференциально-разностных уравнений 𝑛 -го порядка данное
определение будет сохраняться, но близость начальных функций решений
будет пониматься в смысле близости (𝑛 − 1)-го порядка. Впрочем, уравнение
𝑛-го порядка можно заметить системой уравнений 1-го порядка и применять
данное выше определение устойчивости для систем уравнений. [6]
21
2.2 Исследование на устойчивость по первому
приближению
Часто для исследований на устойчивость по первому приближению
дифференциально-разностных
уравнений,
стационарных
в
первом
приближении применяется метод, обоснованный в работах Белмана и Райта.
Для тех уравнений первого приближения, которые имеют переменные
коэффициенты или переменные запаздывания в настоящее время не дано
обоснование метода исследования на устойчивость по первому приближению.
Результаты Белмана в применении к одному уравнению 1-го порядка с одним
запаздывающим аргументом можно сформулировать так:
𝑥̇ (𝑡) = 𝑎𝑥(𝑡) + 𝑏𝑠(𝑡 − 𝜏) + 𝑐(𝑡)𝑥(𝑡) + 𝑑(𝑡)𝑥(𝑡 − 𝜏) +
𝑖
𝑗
+ ∑∞
𝑖+𝑗=2 𝑎𝑖𝑗 (𝑡)𝑥 (𝑡)𝑥 (𝑡 − 𝜏),
(2.2)
где 𝑎 и b это постоянные, c(t) и d(t) – непрерывные функции t, такие, что
lim 𝑐(𝑡) = lim 𝑑(𝑡) = 0 ; все непрерывные функции 𝑎𝑖𝑗 (𝑡) ограничены по
𝑡→∞
𝑡→∞
𝑖 𝑗
модулю, а ряд ∑∞
𝑖+𝑗=2 𝑎𝑖𝑗 (𝑡)𝑦 𝑧 сходится при достаточно малых по модулю y
и 𝑧, 𝑖 и j- целые положительные числа.
Уравнение (2.2) имеет асимптотически устойчивое тривиальное
решение, если действительные части всех корней характеристического
уравнения 𝑠 = 𝑎 + 𝑏𝑒 −𝑠𝑡 удовлетворяют неравенству 𝑅(𝑠𝑖 ) ≤ −𝛼 < 0.
Результат Райта обобщает теорему Белмана.
Теорема Райта
Тривиальное решение уравнения
𝑛
𝜈
𝜈
∑𝑚
𝜇=0 ∑𝜈=0 𝑎𝜇𝜈 𝑥 (𝑡 + 𝜏𝜇 ) + 𝜓(𝑡, 𝑥 (𝑡 + 𝜏𝜇 )) = 0,
(2.3)
асимптотически устойчиво, если все корни 𝑠𝑖 характеристического уравнения
22
𝑚
𝑛
∑ ∑ 𝑎𝜇𝜈 𝑠 𝜈 𝑒 𝜏𝜇𝑠 = 0
𝜇=0 𝜈=0
удовлетворяют неравенству 𝑅(𝑠𝑖 ) ≤ −𝛼 < 0.
В уравнении (2.3) все 𝑎𝜇𝜈 и 𝜏𝜇 – постоянные,
𝑎𝑚𝑛 ≠ 0, 0 ≤ 𝜏0 < 𝜏1 < ⋯ < 𝜏𝑚 ; функция
𝜓 (𝑡, 𝑥 𝜈 (𝑡 + 𝜏𝜇 )) =
= 𝜓(𝑡, 𝑥(𝑡 + 𝜏0 ), 𝑥(𝑡̇ + 𝜏0 ), … , 𝑥 (𝑛−1) (𝑡 + 𝜏0 ), … , 𝑥(𝑡 + 𝜏𝑚 ), … , 𝑥 (𝑛−1) (𝑡 + 𝜏𝑚 ),
𝜈 = 0, 1, … , (𝑛 − 1); 𝜇 = 0, 1, … , 𝑚,
непрерывна по всем аргументам в области R, в которой 𝑡 ≥ 0, а остальные
аргументы по модулю меньше некоторой постоянной 𝜎 > 0. Помимо этого,
𝑛−1
функция ψ должна удовлетворять условию: если ∑𝑚
𝜇=0 ∑𝜈=0 |𝜔𝜇𝜈 | < 𝜔 , то
|𝜓(𝑡, 𝜔𝜇𝜈 )| ≤ 𝜔𝜒(𝜔) , где 𝜒(𝜔) – ограниченная функция 𝜔 , такая, что
𝑎 𝜒2 𝜔
∫0
𝜔
𝑑𝜔 сходится при 𝑎 > 0.
Схема доказательства теоремы Райта.
Для доказательства теоремы достаточно обнаружить, что
|𝑥
(𝜈)
𝜑
(𝑡)| < 𝜀𝑒 −𝑐𝑡 , 𝜈 = 0,1, … , (𝑛 − 1)
(2.4)
для всех 𝑡 ≥ 0, где 0 < 𝑐 < 𝛼, |𝜑 𝜈 (𝑡) < 𝛿| на начально множестве 𝐸0 , а 𝜀
– положительная постоянная. Допустим, что неравенства (2.4) справедливы не
для всех 𝑡 ≥ 0; тогда найдутся числа 𝑡 = 𝜉 и 𝜈 = 𝜆, при которых
|𝑥 (𝜆) (𝜉) = 𝜀𝑒 −𝑐𝜉 |и |𝑥 (𝜈) (𝑡)| ≤ 𝜀𝑒 −𝑐𝑡 ,
0 ≤ 𝑡 ≤ 𝜉,
0 ≤ 𝜈 ≤ (𝑛 − 1)
Заменим теперь уравнение (2.3), которое мы кратко запишем в виде
𝐿(𝑥) + 𝜓 = 0, линейным неоднородным уравнением
23
𝐿(𝑧) = 𝜈(𝑡),
(2.5)
где
− 𝜓̇ (𝑡, 𝑥
𝜈(𝑡) = {
(𝜈)
𝜑
(𝑡 + 𝜏𝜇 )) при 0 ≤ 𝑡 ≤ 𝜉 − 𝜏𝑚 ,
0 при 𝑡 > 𝜉 − 𝜏𝑚 ,
решение которого 𝑧(𝑡), как нетрудно проверить, будет совпадать с решением
𝑥
𝜑 (𝑡) уравнения (2.3) при , 0 ≤ 𝑡 ≤ 𝜉 , если для 𝑧(𝑡) взять следующие
условия:
𝑧(𝑡) =
𝑥
𝜑 (𝑡),
0 ≤ 𝑡 ≤ 𝜏𝑚 .
Далее доказывается, что функция 𝑧(𝜆) (𝑡) удовлетворяет всем условиям,
которым должен удовлетворять оригинал при преобразовании Лапласа; а это
значит, что при применении преобразования Лапласа к (2.5) можно явно
выразить функцию 𝑧(𝜆) (𝑡) и, построив ее оценку, показать, что
|𝑧(𝜆) (𝑡)| <
1
2
𝜀𝑒 −𝑐𝜉 , что противоречит предположению (2.4). [6]
2.3 Методы нахождения областей асимптотической
устойчивости дифференциально-разностных уравнений,
стационарных в первом приближении
Если к дифференциально-разностным уравнениям применим метод
исследования на устойчивость по первому приближению, то исследование
таких уравнений, стационарных в первом приближении, сводится к
определению
знаков
действительных
частей
корней
квазиполинома.
Вычисление корней квазиполинома – трудная задача даже в простейшем
случае, поэтому имеет смысл использовать признаки отрицательности
вещественных частей всех корне квазиполинома.
24
Амплитудно-фазовый метод в применении к линейному однородному
дифференциальному уравнению с запаздывающим аргументом с одним
постоянным запаздыванием 𝜏 и с постоянными коэффициентами сводится к
следующему: характеристическое уравнение имеет вид
𝛷(𝑧) + 𝛹(𝑧)𝑒 −𝜏𝑧 = 0,
или кратко 𝑓(𝑧) = 0, где 𝛷(𝑧) и 𝛹(𝑧) – полиномы, причем 𝛷(𝑧) имеет
более высокую степень, чем 𝛹(𝑧) . Рассмотрим образ 𝛤𝜔 контура 𝛾𝑧 ,
состоящего из отрезка мнимой оси −𝑅 ≤ 𝑦 ≤ 𝑅
и полуокружности
|𝑧| = 𝑅 ,𝑅𝑒(𝑧) > 0, радиуса 𝑅 в плоскости 𝜔: 𝜔 = 𝑓(𝑧); согласно принципу
аргумента, если при однократном обходе контура 𝛾𝑧 в положительном
направлении его образ 𝛤𝜔 обходит точку 𝜔 = 0 в положительном направлении
𝑁 − 𝑃 раз. Если точнее, то коэффициент зацепления контура 𝛤𝜔 с точкой
𝜔 = 0 равен 𝑁 − 𝑃 . 𝑃 – сумма кратностью полюсов функции 𝑓(𝑧) внутри
контура. Тогда сумма кратности нулей функции 𝑓(𝑧), находящихся внутри
контура 𝛾𝑧 , равна 𝑁; при этом предполагается, что на контура 𝛾𝑧 нет нулей и
полюсов функции 𝑓(𝑧).
Если на мнимой оси нет нулей функции 𝑓(𝑧) и при 𝑅 → ∞ число 𝑁
остается равным нулю, то все нули функции 𝑓(𝑧) расположены в левой
полуплоскости 𝑅𝑒(𝑧) < 0 и точка покоя 𝑥 ≡ 0 исходного дифференциального
уравнения асимптотически устойчива. Но вместо отображения 𝑓(𝑧) будет
лучше рассмотреть отображение, определяемое функцией
𝑓(𝑧)
𝛷(𝑧)
=1+
𝛹(𝑧)
𝛷(𝑧)
𝑒 −𝜏𝑧 ,
имеющей те же нули, что и функция 𝑓(𝑧), и полюсы, совпадающие с
нулями функции 𝛷(𝑧) (в данной работе мы предполагаем, что функции 𝑓(𝑧) и
𝛷(𝑧) не имеют общих нулей), или даже отображение
𝜔𝜏 (𝑧) = −
25
𝛹(𝑧)
𝛷(𝑧)
𝑒 𝜏𝑧 .
В последнем случае возникает необходимость считать число обходов
образом контура 𝛾𝑧 при отображении 𝜔𝜏 (𝑧) точки 𝜔 = 1.
Так как при 𝑅 → ∞ образ полуокружности |𝑧| = 𝑅 , 𝑅𝑒(𝑧) > 0 ,
стягивается в точку 𝜔 = 0, то остается лишь построить образ при отображении
𝜔𝜏 (𝑧) мнимой оси, проходимой в направлении убывания 𝑦 . Этот образ
называется амплитудно-фазовой характеристикой. Для ее нахождения строят
так называемую предельную характеристику – образ мнимой оси при
отображении 𝜔0 (𝑧) = −
𝛹(𝑧)
𝛷(𝑧)
, а после учитывают влияние множителя
𝑒 −𝑖𝜏𝑦 : 𝜔𝜏 (𝑖𝑦) = 𝜔0 (𝑖𝑦)𝑒 −𝑖𝜏𝑦 . [6]
26
Глава 3. Программная реализация
В настоящее время вычислительная техника не позволяет находить
аналитическое решение многих задач. К числу таких задач относится и
нахождение решения дифференциально-разностных уравнений. Однако
численные методы решения позволяют находить приближенные решения с
минимальными погрешностями. В данной работе используется метод Адамса
для численного решения дифференциального уравнения. Стандартный
алгоритм
может
дифференциальных
использоваться
уравнений,
для
решения
обыкновенных
поэтому
возникает
необходимость
модифицировать его на случай дифференциально-разностного уравнения с
запаздывающим аргументом. Схема такого модифицированного метода
Адамса приводится далее в работе.
3.1 Метод Адамса численного решения дифференциальноразностного уравнения
Предположим найдено несколько приближений 𝑥𝑗 в моменты времени
𝑥𝑗 = 𝑥(𝑡𝑗 ),
Требуется
построить
𝑗 = ̅̅̅̅
0, 𝑖, 𝑡𝑗 = 𝑡0 + 𝑗∆.
правило
вычисления
𝑥𝑖+1 ≈ 𝑥(𝑡𝑖 + 1)
на
следующем шаге. Для этого проинтегрируем правую и левую часть системы
𝑥̇ (𝑡) = 𝑓(𝑥(𝑡), 𝑥(𝑡 − 𝑇)) по промежутку [𝑡𝑖 ; 𝑡𝑖+1 ]:
𝑡
𝑥(𝑡𝑖+1 ) = 𝑥(𝑡𝑖 ) + ∫𝑡 𝑖+1 𝑓(𝑥(𝑡), 𝑥(𝑡 − 𝑇))𝑑𝑡.
𝑖
Заменив
подынтегральную
функцию
(3.1)
𝑓(𝑥(𝑡), 𝑥(𝑡 − 𝑇))
интерполирующим её многочленом 𝑃𝑘 (𝑡) , будем считать, что дискретные
приближенные значения
𝑓𝑗 = 𝑓(𝑥𝑗 , 𝑥𝑇𝑗 ) ≈ 𝑓 (𝑥(𝑡𝑗 ), 𝑥(𝑡𝑗 − 𝑇))
27
известны.
При интерполировании назад из угла 𝑡𝑖 имеем
𝑞(𝑞 + 1) 2
∆ 𝑓𝑖−2 +
2!
𝑞(𝑞 + 1)(𝑞 + 2) 3
𝑞(𝑞 + 1)(𝑞 + 2) × … × (𝑞 + 𝑘 − 1) 𝑘
+
∆ 𝑓𝑖−3 + ⋯ +
∆ 𝑓𝑖−𝑘 ,
3!
𝑘!
𝑃𝑘 = 𝑃𝑘 (𝑡𝑖 + 𝑞ℎ) = 𝑓𝑖 + 𝑞∆𝑓𝑖−1 +
из узла 𝑡𝑖+1
𝑞(𝑞 + 1) 2
∆ 𝑓𝑖−1 +
2!
𝑞(𝑞 + 1)(𝑞 + 2) 3
𝑞(𝑞 + 1)(𝑞 + 2) × … × (𝑞 + 𝑘 − 1) 𝑘
+
∆ 𝑓𝑖−2 + ⋯ +
∆ 𝑓𝑖−𝑘+1 .
3!
𝑘!
Подставим многочлены 𝑃𝑘 (𝑡), 𝑃̃𝑘 (𝑡) в равенство (3.1), используя
𝑃̃𝑘 = 𝑃̃𝑘 (𝑡𝑖+1 + 𝑞ℎ) = 𝑓𝑖+1 + 𝑞∆𝑓𝑖 +
конечные разности, получаем формулы для вычисления очередного значения
𝑥𝑖+1 ≈ 𝑥(𝑡𝑖+1 )
𝑡
𝑥𝑖+1 = ∫𝑡 𝑖+1 𝑃𝑘 (𝑡)𝑑𝑡,
(3.2)
𝑡
𝑥𝑖+1 = ∫𝑡 𝑖+1 𝑃̃𝑘 (𝑡)𝑑𝑡.
(3.3)
𝑖
𝑖
На основании формул (3.2), (3.3) получаем два семейства методов
Адамса.
Рассмотрим первое семейство методов (3.2). Сделав замену переменных
𝑡 = 𝑡𝑖 + 𝑞ℎ в пределах интегрирования, формулу (3.2) можно переписать в
виде
𝑥𝑖+1 = 𝑥𝑖 + ℎ𝐼𝑘 ,
где
1
𝐼𝑘 = ∫0 𝑃𝑘 (𝑡𝑖 + 𝑞ℎ)𝑑𝑞 = [𝑓𝑖 𝑞 +
1 𝑞4
1
𝑞5
+ ( + 𝑞 3 + 𝑞 2 ) ∆3 𝑓𝑖−1 + ( +
6 4
24 5
+
∆𝑓𝑖−1
2
+
5∆2 𝑓𝑖−2
12
1
251
6
720
+ ∆3 𝑓𝑖−3 +
3𝑞 4
2
+
𝑞2
2
∆𝑓𝑖−1 + (
11𝑞 3
3
𝑞3
6
+
𝑞2
4
) ∆2 𝑓𝑖−2 +
+ 3𝑞 2 ) ∆4 𝑓𝑖−4 + ⋯ ] = 𝑓𝑖 +
∆4 𝑓𝑖−4 + ⋯
Таким образом, ограничиваясь методом Адамса четвертого порядка,
получаем
1
5
2
12
𝑥𝑖+1 = 𝑥𝑖 + ℎ (𝑓𝑖 + ∆𝑓𝑖−1 +
1
251
6
720
∆2 𝑓𝑖−2 + ∆3 𝑓𝑖−3 +
28
∆4 𝑓𝑖−4 ).
(3.4)
Сделаем аналогичную замену в пределах интегрирования в (3.3)
𝑡 = 𝑡𝑖+1 + 𝑞ℎ и, подставляя выражение многочлена 𝑃̃𝑘 (𝑡) под знак
интеграла, получаем равенство
𝑥𝑖+1 = 𝑥𝑖 + ℎ𝐼̃𝑘 ,
где
2
3
2
0
𝑞
𝑞
𝑞
𝐼̃𝑘 = ∫−1 𝑃̃𝑘 (𝑡𝑖 + 𝑞ℎ)𝑑𝑞 == [𝑓𝑖+1 𝑞 + ∆𝑓𝑖 + ( + ) ∆2 𝑓𝑖−1 +
2
6
4
1 𝑞4
1
𝑞5
+ ( + 𝑞 3 + 𝑞 2 ) ∆3 𝑓𝑖−1 + ( +
6 4
24 5
−
∆𝑓𝑖
2
−
∆2 𝑓𝑖−1
12
−
1
24
∆3 𝑓𝑖−2 −
Следовательно,
19
720
схема
3𝑞 4
2
+
11𝑞 3
3
+ 3𝑞 2 ) ∆4 𝑓𝑖−4 + ⋯ ] = 𝑓𝑖+1 −
∆4 𝑓𝑖−3 + ⋯
вычисления
следующего
значения
при
∆4 𝑓𝑖−3 ).
(3.5)
интерполировании из угла 𝑡𝑖+1 выглядит следующим образом
1
1
2
12
𝑥𝑖+1 = 𝑥𝑖 + ℎ (𝑓𝑖+1 − ∆𝑓𝑖 −
∆2 𝑓𝑖−1 −
1
24
∆3 𝑓𝑖−2 −
19
720
Можно заметить, что первое семейство методов (3.4) является явным, а
второе (3.5) неявным.
Предположим, что правые части уравнений (𝑘 + 1) раз непрерывно
дифференцируемые функции. Тогда ошибка на шаге составляет
𝑓 𝑘+1 (𝜉, 𝑥(𝜉))
𝑅𝑘 (𝑥) =
𝑞(𝑞 + 1) × ⋯ × (𝑞 + 𝑘)ℎ𝑘+1 .
(𝑘 + 1)!
Таким образом, локальная погрешность метода (3.5), вычисляемая по
1
формуле ℎ ∫0 𝑅𝑘 (𝑥𝑖 + 𝑞ℎ)𝑑𝑞 , составляет величину порядка 𝑂(ℎ𝑘+2 ) , а
глобальная – 𝑂(ℎ𝑘+1 ). Величина погрешности для второго метода (3)
𝑅𝑘 (𝑥) =
𝑓(𝜉,𝑥(𝜉))
(𝑘−1)!
𝑞(𝑞 + 1) × ⋯ × (𝑞 + 𝑘)ℎ𝑘−1 . [8]
3.2 Пример построения области асимптотической
устойчивости
Рассмотрим применение амплитудно-фазового метода к простейшему
уравнению
29
𝑥̇ (𝑡) + 𝑎𝑥(𝑡) + 𝑏𝑥(𝑡 − 𝜏) = 0.
Характеристическая функция этого уравнения имеет вид:
𝑓(𝑧) = 𝑧 + 𝑎 + 𝑏𝑒 −𝜏𝑧 ,
а построенные на его основе новые отображения:
𝑓(𝑧)
𝛷(𝑧)
=1+
𝑏
𝑧+𝑎
𝑒 −𝜏𝑧 , 𝜔𝜏 (𝑧) = −
𝑏
𝑧−𝑎
𝑒 −𝜏𝑧 , 𝜔0 (𝑧) = −
𝑏
𝑧−𝑎
.
Предельной характеристикой является окружность с центром в точке
𝜔= −
𝑏
2𝑎
𝑏
радиуса | | . Очевидно, что при 𝑎 > 0, |𝑏| < 𝑎 при любом 𝜏
2𝑎
амплитудно-фазовая характеристика не включает в себя точку 𝜔 = 1, а значит,
при этих условиях решение 𝑥 ≡ 0 асимптотически устойчиво. При 𝑎 < 0 у
функции 𝜔𝜏 (𝑧) есть полюс в правой полуплоскости 𝑅𝑒(𝑧) > 0 , поэтому
приходится учитывать влияние множителя 𝑒 −𝑖𝜏𝑦 . В результате этого
исследования получаем область устойчивости, показанную на рисунке 3.1.
𝑏
𝜋
2𝜏
−
1
𝜏
Область устойчивости
0
𝑎
Рис.3.1 - Область устойчивости решения 𝑥 ≡ 0 уравнения
𝑥̇ (𝑡) + 𝑎𝑥(𝑡) + 𝑏𝑥(𝑡 − 𝜏) = 0.
30
Изменяя
параметры
𝑎и𝑏
в
программной
реализации
модифицированного алгоритма Адмса, программа строит графики решения,
которые будут представлены далее на рисунках 3.2 – 3.3.
Рис. 3.2 – b = -0.01, a = 0.1, решение асимптотически устойчиво.
31
Рис. 3.3 – b = 100, a = 0.001, решение неустойчиво
Так как верхняя граница области асимптотической устойчивости
зависит от параметра 𝜏, имеет смысл рассмотреть влияние этого параметра при
одинаковых значениях 𝑎 и 𝑏. Примеры представлены на рисунках 3.4 – 3.7.
Значения 𝑎 = 0.01, 𝑏 = 10.
32
Рис. 3.4 – 𝜏 = 0.1, решение устойчиво
Рис. 3.5 - 𝜏 = 1, решение неустойчиво
33
Рис. 3.6 - 𝜏 = 0.1571, решение устойчиво
Рис. 3.7 - 𝜏 = 0.1572, решение неустойчиво
34
Выводы
В ходе данной работы была подробно рассмотрена проблема построения
области асимптотической устойчивости для одного дифференциальноразностного уравнения с запаздывающим аргументом. Для этого в первой и
второй части работы были рассмотрены основные сведения по теории
обыкновенных
дифференциальных
и
дифференциально-разностных
уравнений.
Существующие алгоритмы численного решения не подходят для
дифференциально-разностных уравнений с запаздывающим аргументом,
поэтому возникает необходимость модифицировать один из них. В третьей
части работы приведен модифицированный алгоритм Адамса [8] для
численного решения.
С его помощью была проанализирована построенная аналитически
область асимптотической устойчивости. Программное построение области
асимптотической устойчивости в настоящее время – достаточно трудоемкий
процесс, так как возникает необходимость вычисления корней квазиполинома.
35
Заключение
В ходе данной работы, были рассмотрены основные методы решения и
исследования на устойчивость дифференциально-разностных уравнений с
запаздывающим аргументом.
В силу того, что стандартный метод Адамса используется для решения
обычных дифференциальных уравнений, был приведен и программно
реализован
метод
Адамса
для
построения
численного
решения
дифференциально-разностных уравнений.
С помощью метода амплитудно-фазовых колебаний была аналитически
построена
область
асимптотической
устойчивости
дифференциально-
разностного уравнения и проведено ее исследование с помощью программной
реализации метода Адамса.
Программная реализация, построенная в ходе данной работы может
быть использована для облегчения исследования области асимптотической
устойчивости дифференциально-разностного уравнения с запаздыванием.
Оперируя изменениями входных параметров уравнения, можно
ускорить процесс анализа и проверить правильность теоретически выкладок,
не затрачивая много времени на точное решение.
Цели и задачи, поставленные в начале работы, были выполнены.
36
Список литературы
1) Петровский И. Г. Лекции по теории обыкновенных дифференциальных
уравнений / Под ред. А. Д. Мышкнса, О. А. Олейник. – М.: Изд-во МГУ,
1984. – 296 с.
2) Эльсгольц Л.Э. Дифференциальные уравнения и вариационное
исчисления. / Эльсгольц Л.Э. – М.: Книга по Требованию, 2012. – 424с.
3) Беллман Р., Кук К. Л. Дифференциально-разностные уравнения. М.: Мир,
1967. – 548 с.
4) Основные понятия теории устойчивости. http://www.math24.ru/основныепонятия-теории-устойчивости.html
5) Ногин В.Д. Теория устойчивости движения. СПбГУ: ф-т ПМ-ПУ, 2008.
6) Л. Э. Эльсгольц, Устойчивость решений дифференциально-разностных
уравнений, УМН, 1954, том 9, выпуск 4(62), 95–112.
7) Условие Липшица.
http://dic.academic.ru/dic.nsf/enc_mathematics/2821/ЛИПШИЦА
8) Жабко А. П., Зараник У. П. Построение области асимптотической
устойчивости дифференциально – разностных уравнений в среде
MATLAB // Сборник избранных трудов V научно-практической
конференции. Москва. 8-10 ноября 2010 г. / Под ред. В. А. Сухомлина.
М.: ИНТУИТ. РУ, 2010. - 640 с.
37
Отзывы:
Авторизуйтесь, чтобы оставить отзыв