О разрешимости на оси автономных дифференциальных уравнений с ограниченным запаздыванием Текст научной статьи по специальности «Математика»

зань, Российская Федерация, аспирант кафедры вычислительной математики, e-mail: ildar.badriev@kpfu.ru

Singatullin Marsel Talgatovich, Kazan Federal University, Kazan, the Russian Federation, Postgraduate Student of the Computing Mathematics Department, e-mail: ildar.badriev@kpfu.ru

Чебаков Юрий Владимирович, Усольский сельскохозяйственный техникум, с. Усолье, Самарская область, Российская Федерация, преподаватель, e-mail: yshk_chebakov@mail.ru

Chebakov Yurii Vladimirovich, Usolsky Agricultural College, Usolye, Samara region, the Russian Federation, Lecturer, e-mail: yshk_chebakov@mail.ru

УДК 517.929

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

© А.С. Баландин

Ключевые слова: функционально-дифференциальное уравнение; запаздывание; разрешимость на оси.

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

Пусть N = {1, 2, 3,...} , М = (-ж, +ж) , М+ = [0, +ж) , М- = (-то, 0] , С — множество комплексных чисел, — пространство локально абсолютно непрерывных функций.

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

Рассмотрим линейное автономное однородное дифференциальное уравнение с ограниченным запаздыванием

г ш

х(г)+/ х(г - в) ¿ф) = 0, г е м, (1)

Jo

где ш > 0, функция г: [0,ш] ^ М(С) ограниченной вариации, г(0) = 0. Интеграл понимается в смысле Римана-Стилтьеса. Вариацию функции г обозначим р = |^г(в)| . Соответствующей заменой переменных к уравнениям вида (1) можно сводить различные классы сингулярных функционально-дифференциальных уравнений, рассматриваемых на конечном отрезке. Для таких уравнений вопросы однозначной разрешимости оказываются нетривиальными.

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

Известно (см. [2, с. 73], [3, с. 98]), что на М+ решения уравнения (1) имеют оценку:

|х(г)| < х(0)ер4.

При изучении уравнений вида (1) на положительной полуоси эта оценка оказывается очень полезной (например, даёт возможность применять интегральные преобразования). Остаётся ли она справедливой на отрицательной полуоси? Ответ на этот вопрос даёт следующий пример, идея которого принадлежит Е.И. Бравому.

1 t

Пример. Рассмотрим функцию p(t) = jlp et2-1 , t € [—1,1] . Используя [4, c. 1051-1052], запишем

те

^(t) = E L°n(i)t2n, (2)

где Ln(t) = Х]П=о(—1)kCntkr — многочлен Лагерра. Функция р бесконечно дифференцируема на отрезке [-1,1] , р(1 - 0) = р(-1 + 0) = 0, p(n)(1 - 0) = p(n)(-1 + 0) = 0 .Из представления (2) следует, что p(2n+1)(0) = 0, p(2n)(0) = Ln(1)(2n)!. Пользуясь асимптотикой многочленов Лагерра [4, с. 1053], получаем

Ln(1) = -1=e1 n-1 cos (2-n - 4) + O (n-.

Функцию р продолжим на всю отрицательную полуось по правилу: p(t) = dttp(t+2) . Легко видеть, что при любом а ^ 0 выполняется limt^—^(p(t)eat) = +то. С другой стороны, функция р удовлетворяет при t € R- частному случаю уравнения (1):

x(t) = x(t - 2).

Итак, существует пример автономного функционально-дифференциального уравнения, имеющего решения, растущие на -то быстрее любой экспоненты. Следовательно, если мы хотим рассматривать решения уравнения (1) с определённой скоростью роста на отрицательной полуоси, нам придётся задавать скорость роста выбором пространства решений. Пусть а — произвольное вещественное число. Рассмотрим следующее функциональное пространство:

Da = \ x € Dloc: sup (|x(t)|e-at) < то I , teR-

и его подпространство вида

D00 = (x € Dioc: lim (|x(t)|e-at) =Л .

i t^—те i

Поставим задачу: найти условия разрешимости и структуру пространства решений уравнения (1) в пространствах Оа и В°а .

2. Некоторые свойства решений автономных функционально-дифференциальных уравнений на М+

Наряду с уравнением (1) на оси рассмотрим уравнение на положителной полуоси:

ГШ

х(г)+/ х(г - в) ^ф) = 0, г е м+. (3)

Jo

При исследовании уравнения (3) важную роль играет характеристическая функция:

ГШ

д(р) = р + е-р йт(£). Jo

Отметим ряд важных свойств решения уравнения (3) и его характеристической функции, установленных в работах [2, 5, 6]. Предложение 1.

1. Функция д аналитическая на всей комплексной плоскости.

2. Функция g имеет счётное множество нулей.

3. Все корни отделены друг от друга, в любой полуплоскости Re p ^ а, а € R, количество нулей функции g конечно.

Обозначим {pn} , n € N, множество всех корней функции g , причём Re pi ^ Re p2 ^ ^ ... ^ Re pn ^ .... Обозначим в1 = Re p1, в2 = max {Re pk}, •••, fin =

keN,Re pk <^1

= max {Repk} , •... Эта процедура не может прерваться, поскольку количество

keN,Re pk <£n-l

корней с одинаковой вещественной частью конечно. Также в работах [2, 5, 6] показано

Предложение2. Фундаментальное решение уравнения (3) представимо в виде

x(t) = J] qn(t)epn t + z(t), t € R+,

Re Pn^A

где Л € R, сумма берется по всем нулям pn функции g, для которых Re pn ^ Л, qn — многочлен, степень которого на единицу меньше кратности pn , для функции z найдутся такие N > 0, y < Л, что справедлива оценка |z(t)| ^ NeYt.

3. Некоторые вспомогательные результаты

Л е м м а 1. Пусть x € Da (или Da), а > 0 . Если функция x на R+ имеет оценку

sup (|x(t)|e7t) < то, Y > 0, то к правой и левой части уравнения (1) можно применить te R+

преобразование Фурье.

Л е м м а 2. Пусть x € Da (или D'à), а > 0 . Если все нули функции g лежат слева от мнимой оси, то в пространстве Da , а > 0 , уравнение (1) имеет только тривиальное решение.

Если n нулей функции g лежат справа от мнимой оси, а остальные нули находятся слева от мнимой оси, то с помощью замены u(t) = x(t) — ^П=1 qk(t)ePkt получаем уравнение

г ш

U(t) + / u(t — s) dr(s) = 0, t € R, (4)

'

причём u € Da' , а' > 0, и выполняется оценка sup (|u(t)|e7t) < то , y > 0. Значит, для

te R+

уравнения (4) справедливы леммы 1 и 2. Откуда вытекает

Л е м м а 3. Пусть x € Da (или Da), а > 0 . Если функция g не имеет нулей на мнимой оси, n нулей лежат справа от прямой Re p = а , а остальные нули находятся слева от прямой Re p = а, то в пространстве Da , а > 0, 'решение уравнения (1) представимо в виде

n

x(t) = £ qk(t)ePkt, k=i

где qk — многочлен, степень которого на единицу меньше кратности pk .

Пусть x € Da , а может быть любого знака. Если ni нулей функции g лежат справа от прямой Re p = а , n2 нулей лежат на прямой Re p = а , а остальные нули находятся слева от прямой Re p = а , то с помощью замены x(t) = v(t)e-et, получаем уравнение

f ш

V(t) — £v(t) + / v(t — s)ees dr(s) = 0, t € R. (5)

'

Для характеристической функции gv , соответствующей уравнению (5), выполняется gv (p) = g(p + в) . Можно подобрать такое в, чтобы v € Da' , а' > 0, и все n1 + n2

нулей функции лежали справа от мнимой оси, на мнимой оси нулей не было, остальные нули лежали слева от мнимой оси. Следовательно, для уравнения (5) справедлива лемма 3.

4. Основной результат

Поставим в соответствие последовательности {/Зп} последовательность классов Брп (0вп ). Отметим несколько простых свойств этих множеств.

а) С С ... С С ...

б) Неравенство а > в имеет место тогда и только тогда, когда Ба С .

в) Включение а е (вП+1,вп) выполняется тогда и только тогда, когда Б1зп+1 С Ба С С 0вп

Множества О00 обладают такими же свойствами.

Рп

В силу предложения 1, последовательность {вп} разбивает вещественную ось на счётный набор интервалов, одному из которых принадлежит а .

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

Теорема 1. Любое 'решение уравнения (1) в пространстве представимо в виде:

х(г) = £ ^(г)е^ + £ СгеР*, (6)

Яе рк >вп Ке Р1=вп

где дк — многочлен, степень которого на единицу меньше кратности рк , ц е М, вторая сумма берётся только по простым корням.

Заметим, что вследствие предложения 1 в представлении решения (6) количество слагаемых конечно.

Частный случай теоремы 1 был получен в работе [7].

Теорема 2. Любое решение уравнения (1) в пространстве Б0п представимо в виде:

х(г) = £ &(г)еРк*,

Яе Рк >в

п

где дк — многочлен, степень которого на единицу меньше кратности рк .

Перейдём к случаю, когда а лежит внутри одного из интервалов (вП+1,вп) .

Следствие 1. Если а е (вп+1, вп), то сужение Ба и Брп+1 (соответственно, О и О0п ) на пространство решения уравнения (1) совпадают.

Следствие 2. Уравнение (1) имеет только тривиальное решение в пространстве Ба тогда и только тогда, когда а > ¡31 .

Следствие 3. Уравнение (1) имеет только тривиальное решение в пространстве О тогда и только тогда, когда а ^ Д .

5. Приложения

Теорема 1 сводит задачи разрешимости уравнения (1) к исследованию расположения нулей характеристической функции д на комплексной плоскости. Эта задача давно и хорошо известна, более того — для многих классов уравнений вида (1) она неплохо изучена [8-12].

Особенно информативен в этом смысле случай а = 0, когда теоремы 1-2 и их следствия сводят рассматриваемую в данной статье задачу к классической задаче устойчивости.

В частности, следствие 2 утверждает эквивалентность однозначной разрешимости уравнения (1) в пространстве ^о и асимптотической (экспоненциальной) устойчивости уравнения (3).

Обозначим: М^ — область, в которой характеристическая функция д имеет ровно к корней справа от мнимой оси. Для уравнений вида (1) с небольшим количеством параметров область Мо допускает аналитическое и графическое описание.

Im £4

Рис. 1: Множества Mk для уравнения (7). Для примера рассмотрим два уравнения:

x(t) + bix(t — 1) = 0, (7)

где bi = |bi|e^ € C ,

x(t) + box(t) + bix(t - 1) = 0, (8)

где b0, b1 € R .

Область Mo для уравнения (7) описывается следующими неравенствами: — П < ф < П , 0 < |bi| < f — |ф| .

Область Mo для уравнения (8) описывается следующими неравенствами: —bo < bi < snö , где в — наименьший положительный корень уравнения bo = —в ctg в. На рис. 1 и 2 приведены графики областей Mk для уравнений (7) и (8).

ЛИТЕРАТУРА

1. Азбелев Н.В., Максимов В.П., Рахматуллина Л.Ф. Введение в теорию функционально-дифференциальных уравнений. М.: Наука, 1991.

2. Беллман Р., Кук К. Дифференциально-разностные уравнения. М.: Мир, 1967.

3. Азбелев Н.В., Симонов П.М. Устойчивость решений уравнений с обыкновенными производными. Пермь: Изд-во Пермск. ун-та, 2001.

4. Г'радштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Физматгиз, 1963.

h

Рис. 2: Множества Mk для уравнения (8).

5. Зубов В.И. К теории линейных стационарных систем с запаздывающим аргументом // Изв. вузов. Математика. 1958. № 6. С. 86-95.

6. Sabatulina T., Malygina V. On positiveness of the fundamental solution for a linear autonomous differential equation with distributed delay // Electron. J. Qual. Theory Differ. Equ. 2014. № 61. P. 1-16.

7. Баландин А. С. О разрешимости на оси некоторых классов дифференциально-разностных уравнений // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2013. Т. 18. № 5-2. С. 2449-2451.

8. Андронов А.А., Майер А.Т. Простейшие линейные системы с запаздыванием // Автоматика и телемеханика. 1946. Т. 7. № 2, 3. С. 95-106.

9. Рехлицкий З.И. Об устойчивости решений некоторых линейных дифференциальных уравнений с запаздывающим аргументом в банаховом пространстве // ДАН СССР. 1956. Т. 111. № 1. С. 29-32.

10. Эльсгольц Л.Э., Норкин С.Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. М.: Наука, 1971.

11. Кирьянен А.И. Устойчивость систем с последействием и их приложения. СПб.: Издательство Санкт-Петербургск. ун-та, 1994.

БЛАГОДАРНОСТИ: Работа поддержана грантом Российского Фонда Фундаментальных Исследований (грант 13-01-96050 р урал а).

Поступила в редакцию 20 июня 2015 г.

Balandin A.S. ON SOLVABILITY OF AUTONOMOUS DIFFERENTIAL EQUATIONS WITH BOUNDED DELAY ON THE AXIS

The representation of solutions for autonomous differential equations with bounded delay on the axis is obtained. The problem of solvability on the axis is reduced to the study of location of zeros of the characteristic function on the complex plane.

Key words: functional differential equation; delay; solvability on the axis.

Баландин Антон Сергеевич, Пермский национальный исследовательский политехнический университет, г. Пермь, Российская Федерация, младший научный сотрудник научно-исследовательского центра «Функционально-дифференциальные уравнения», e-mail: balandin-anton@yandex.ru

Balandin Anton Sergeevich, Perm State National Research University, Perm, the Russian Federation, Junior Researcher of the Research Center «Functional-Differential Equations», e-mail: balandin-anton@yandex.ru

УДК 517.977, 519.711

ОБ УСЛОВИЯХ НАБЛЮДАЕМОСТИ ПОЭТАПНО МЕНЯЮЩИХСЯ УПРАВЛЯЕМЫХ ЛИНЕЙНЫХ СИСТЕМ

© В.Р. Барсегян

Ключевые слова: поэтапно меняющаяся линейная система, измерение, вполне наблюдаемость, условия вполне наблюдаемости.

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

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

Управляемость и наблюдаемость - два фундаментальных понятия теории управления и являются принципиальными как в задачах для обычных систем [1-5], так и в задачах управления и наблюдения для поэтапно меняющихся систем. Вопросы управляемости и наблюдаемости указанных систем переменной структуры исследованы в работах [6-12]. В работах [7-9] получены необходимое и достаточное условие вполне управляемости поэтапно меняющейся линейной системы. В работах [11, 12] получены необходимое и достаточное условие для управляемости и наблюдаемости кусочно линейных импульсных управляемых систем.

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