Построение области
асимптотической устойчивости
дифференциально – разностных уравнений
в среде MATLAB
, Санкт – Петербургский государственный университет, профессор, *****@; , Санкт – Петербур- ский государственный университет, аспирантка, *****@***ru.
Рассматривается приближение решения дифференциально-разностной системы решением разносной системы. Построено приближение области асимптотической устойчивости дифференциально –разностного уравнения. Описаны способы описания области в функциональном пространстве. Приведен конструктивный пример.
Введение.
Рассмотрим дифференциально-разностную стационарную систему вида:
(1)
где 
– непрерывная функция своих аргументов, 
– время запаздывания, с начальной функцией 
, где 
. Будем предполагать, что система (1) имеет нулевое решение 
и линейное приближение 
экспоненциально устойчиво по Ляпунову.
Для построения области асимптотической устойчивости дифференциально – разностной системы уравнений запаздывающего типа был разработан подход, состоящий из двух этапов. Первый этап заключается в приближении решений дифференциально-разностной системы уравнений решениями разностной системы и оценки близости полученных решений. Второй этап состоит из построения области асимптотической устойчивости разностной системы уравнений и описания приближаемой области асимптотической устойчивости дифференциально-разностной системы, используя оценки приближения соответствующих решений, полученных на первом этапе.
Схема метода Адамса [1].
Предположим найдено несколько приближений 
в моменты времени 
, 
, 
. Требуется построить правило вычисления 
на следующем шаге. Для этого проинтегрируем правую и левую часть системы (1) по промежутку 
:
(2)
Заменив подынтегральную функцию 
интерполирующим её многочленом 
, будем считать, что дискретные приближенные значения 
известными. При интерполировании назад из угла 
имеем
из узла 
Подставим многочлены 
в равенство (2), используя конечные разности, получаем формулы для вычисления очередного значения
(3)
(4)
На основании формул (3), (4) получаем два семейства методов Адамса.
Рассмотрим первое семейство методов (3). Сделав замену переменных 
в пределах интегрирования, формулу (3) можно переписать в виде
,
где
Таким образом, ограничиваясь методом Адамса четвертого порядка, получаем
(5)
Сделаем аналогичную замену в пределах интегрирования в (4) 
и, подставляя выражение многочлена 
под знак интеграла, получаем равенство
,
где
Следовательно, схема вычисления следующего значения при интерполировании из угла 
выглядит следующим образом
(6)
Можно заметить, что первое семейство методов (5) является явным, а второе (6) неявным.
Предположим, что правые части уравнений 
раз непрерывно дифференцируемые функции. Тогда ошибка на шаге составляет
Таким образом, локальная погрешность метода (6), вычисляемая по формуле 
, составляет величину порядка 
, а глобальная – 
. Величина погрешности для второго метода (5)
Следовательно, справедливо следующее утверждение.
Утверждение 1. Локальная погрешность методов Адамса (5), (6) составляет и 
соответственно.
Основной результат
Рассмотрим семейство разностных систем
(7)
Линейное приближение системы (7) экспоненциально устойчиво, то есть справедливо предельное соотношение
при 
где 
решение системы (7) при 
Лемма 1. Для любого 
, существует 
такое, что линейное приближение системы (7) является экспоненциально устойчивым по Ляпунову.
Рассмотрим приближение области асимптотической устойчивости разностного уравнения 
и дифференциально-разностного уравнения
.
В область асимптотической устойчивости разностного уравнения входят кусочно-непрерывные функции с ограниченной первой производной 
, для которых выполнено условие стыковки
где 
– точки разбиения отрезка 
Лемма 2. Пусть 
где 
. Тогда для любого 
такое, что справедливо неравенство
.
Доказательство. 
при 
Тогда 
Теорема. 
такое, что 
справедливо включение
Доказательство. 
линейное приближение системы (7) экспоненциально устойчиво, следовательно 
такое, что 
Согласно Лемме 2 
, следовательно 
.
Рассмотрим свойства области асимптотической устойчивости дифференциально-разностного уравнения 
:
· ПОЛОЖИТЕЛЬНАЯ ИНВАРИАНТНОСТЬ
Доказательство. Для любой начальной функции 
, следует что решение 
определенно при 
и 
при 
.
· 
– ОТКРЫТОЕ МНОЖЕСТВО
· – СВЯЗНОЕ МНОЖЕСТВО
Доказательство. Для каждого 
и 
из 
следует, что существует 
при 
такая, что 
при 
.
Методы описания области асимптотической устойчивости дифференциально-разностных систем.
· Сечение области асимптотической устойчивости при 
· Сечение области асимптотической устойчивости при 
· Огибающая области асимптотической устойчивости
· 
;
· Пространство среднеквадратичных значений
· 
· Эмиттанс
· 
Пример.
Рассмотрим реактор, состоящий из двух активных зон[2]. Пренебрежем запаздывающими нейтронами и температурными эффектами. Будем считать, что реактивность в зонах реактора зависит только от уровня мощности в зонах. Уравнения динамики имеют вид
(8)
где 
– время жизни нейтронов; 
– коэффициент, характеризующий связь между первой и второй активными зонами; 
, где 
– величина отрицательного мощностного коэффициента реактивности, а 
– мощность; – время запаздывания (
).
Приведем результаты исследования данной системы уравнений на устойчивость двумя различными способами.
Первый способ заключается в аналитическом исследовании корней квазиполинома линейной части рассматриваемой системы (8), который имеет вид
Таким образом, разрешая уравнения относительно 
, получаем следующий результат
Здесь 
, где 
и 
.
Очевидно, устойчивость решения системы зависит от параметров системы (
). Данный вывод улучшает результат исследования области асимптотической устойчивость полученный в [2], с помощью теоремы Красовского и функций Ляпунова. На рисунке 1 показана зависимость параметров 
системы при исследовании на устойчивость с помощью квазиполинома системы (8).
Второй способ заключается в построении области асимптотической устойчивости системы с запаздыванием с помощью вышеизложенного метода Адамса в среде MATLAB. Построение графиков и их анализ подтверждают существование зависимости области асимптотической устойчивости от параметров системы.
Описание программы.
Задавая определенные параметры () системы (8), программа обрабатывает данные, строит график решения с помощью метода Адамса и выводит заданные проекции области асимптотической устойчивости.
Рис. 1. Зависимость параметров системы Покажем на рисунках некоторые
графики, получаемые с помощью программы реализованной в MALAB, а также некоторые способы визуализацию области асимптотической. На рис 2. изображено сечение области асимптотической устойчивости системы (8), где область асимптотической устойчивости дифференциально-разностного уравнения без ограничения на производную (множество B) и область асимптотической устойчивости дифференциально-разностного уравнения с ограничением на производную (множество C)
.
На рис. 3 представлена огибающая области асимптотической устойчивости, по осям 
, 
- координаты начальной функции 
в точке 
Множества точек B и C определяются также,
Рис 2. Сечение области асимптотической Рис 3. Огибающая области асимптотической устойчивости при устойчивости.
На рис. 4 изображен эмиттанс, который показывает распределение плотности точек на поверхности области асимптотической устойчивости. Чем светлее поверхность, тем больше плотность точек, которые принадлежат области асимптотической устойчивости.
Рис 4. Эмиттанс области асимптотической устойчивости
дифференциально – разностной системы (8).
Второй способ построения области асимптотической устойчивости был реализован в среде MATLAB. Таким образом, подтвердился результат, полученный при аналитическом исследовании квазиполинома системы (8). Приведенные способы описания области асимптотической устойчивости дифференциально-разностных системы дают более полную картину при исследовании системы на устойчивость. Полученный результат можно рассматривать как основу для написания раздела в MATLAB, посвященной реализации исследования дифференциально – разностных систем.
Литература
1. , , Кобельков методы. М.: Наука, 19с.
2. Горяченко теории устойчивости в динамике ядерных реакторов. М.: Атомиздат, 19с.
