CC BY
15
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
МАТЕМАТИКА
УДК 517.977
В. М. Марченко, доктор физико-математических наук, профессор, заведующий кафедрой высшей математики (БГТУ); И. К. Асмыкович, кандидат физико-математических наук, доцент (БГТУ)
О ПРЕДСТАВЛЕНИИ РЕШЕНИЙ ДЕСКРИПТОРНЫХ СИСТЕМ С ЗАПАЗДЫВАНИЕМ*
Рассматривается проблема аналитического представления решений линейных дескриптор-ных систем с запаздывающим аргументом. При некоторых предположениях дескрипторная система сводится к гибридной системе специального вида, что позволяет получить интегральное представление решений дескрипторной системы. Дается представление решений рассматриваемой дескрипторной системы в виде ряда по решениям определяющих уравнений соответствующей гибридной системы.
The paper considers the problem of an analytical representation for solutions of linear descriptor time-delay systems. Under some assumptions, we reduce such a system to a special hybrid system with delay and, as a result, we obtain an integral representation for solutions of the descriptor system. Using the defining equation techniques developed for the hybrid system, we give a series representation for solutions of the descriptor system under consideration.
Введение. Активное изучение в качественной теории управления динамическими системами математических моделей, в которых наряду с динамическими (дифференциальными) встречаются и алгебраические (функциональные) зависимости, продолжается в течение последних 30 лет. Системы, которые описывают такие процессы, носят различные наименования, из которых выделим сингулярные [2, 3, 8, 12], дескрипторные [5, 7, 10, 15] и дифференциально-алгебраические (DAE) [6, 11, 13]. Теория таких систем в регулярном случае развита достаточно полно, им посвящено более 1000 работ. Для них предложены различные формы явной и неявной записи решений, изучена зависимость решений от начальных условий и разобраны классы совместимых начальных условий для регулярных систем, рассмотрены задачи управляемости, наблюдаемости и двойственности [8]. Большое число работ посвящено исследованию структуры дескрипторных систем [4, 8, 15], в частности задачам стабилизации, модального управления, расщепимости, регуляризации и нормализации, и изучению влияния линейной обратной связи на качественные характеристики систем. Необходимость учета влияния эффекта запаздывания, часто возникающего в прикладных задачах,
приводит к дескрипторным системам с запаздыванием [2-14], которые изучены гораздо слабее. Вопрос об аналитическом представлении их решений был поставлен в монографии [1], где было отмечено существенное отличие таких систем от невырожденного случая. Для дескрипторных систем в конкретных частных случаях описаны представления решения [2, 3, 6, 8, 13], в частности, для гибридных систем [6], рассмотрены задачи стабилизации [5, 7, 14], особенно подробно для двумерных систем [7]. Был разработан специальный разрешающий оператор [13], позволяющий свести систему к обыкновенной, и предложен алгоритм вычисления переходной матрицы-функции [12].
Постановка задачи. Предварительные результаты. Рассмотрим объект управления, описываемый дескрипторной дифференциально-разностной системой с запаздыванием вида
¿¡х(Г) = /(Г, х(0, х(1 -т), и(Г)), (1)
где 5 - постоянная п х т-матрица с det S = 0, если п = т; х^) е Ят - фазовый вектор, t е +<х>), т> 0; /"х, у, г) есть отображение К+х Ят х Ят х К ^ Я"; и^) е К - управляющее воздействие, достаточное число раз дифференцируемое. Начальные условия заданы таким
Работа выполнена в рамках сотрудничества с Белостокским техническим университетом (S/WI/00.2011).
образом, чтобы система имела единственное решение. Поскольку в практических исследованиях обычно систему линеаризуют вдоль изучаемых режимов, то в общем случае получаем линейную дескрипторную функционально-дифференциальную систему с последействием нейтрального типа
d_
dt
f
\
x(t) - J dsG(t, s)x(t + s)
-h
0
= J dsA(t, s)x(t + s) +
-h
0
+ J dsB(t, s)u(t + s), t > t0.
-h
(2)
Здесь О(V, 5), Л(^ я), Вя) - матрицы-функции соответствующих размеров с ограниченным изменением по второму аргументу на промежутке [—к, 0]; к - некоторое положительное число; и(-) - управляющее внешнее воздействие.
В работе [8] изучен вопрос о разрешимости дескрипторной системы с запаздыванием в простейшей форме
) = ЛхЦ) + Л х(Г — к) + / (V) (3)
с начальными условиями
Хо = ф(х) е С ([—к, 0], Я"), т е [—к,0), (4)
где - постоянная п х п особая матрица с гапкБ = т; Л, Л1 - постоянные матрицы, к > 0.
Эта система является гибридной, т. е. обычно она содержит как дифференциально-разностные, так и алгебраические уравнения.
Выполнив правые и левые элементарные преобразования над матрицей 5, можно получить ее каноническую форму Смита в виде
TlSGl =
Ii 0 0 0
где Т1,01 - невырожденные матрицы.
Тогда умножив систему (3) слева на матри-
цу
T
одновременно с заменой переменных
x = G1y, получим следующую формулу:
7j 0" y (t) = " Ajj A12 " yj(t)"
0 0 _ A21 A22 _ _ У2 (t)"
Г C C 1 Mi 12 " yj(t - h)" Г fj(t ) "
C C y^ 21 22 _ _ У2 (t - h) _ _ f2(t) _
(5)
Выполнив аналогичные преобразования с матрицей Л22 и матрицей в правом нижнем углу новой системы при — к), запишем исходную систему в виде [8]
Ij 0 0 0 1 ~zx(t) 1 rBjj B12 B13 B14
0 02 0 0 ¿2(t) B21 12 0 0
0 0 0з 0 ¿3(t) B31 0 0з 0
0 0 0 04 _ _ZT4(t)_ _B41 0 0 04
Г Zj(t ) "
Z2(t)
Z3(t)
_ Z4 (t)"
Cjj C 12 C M3 C 1 14 Г Z1(t - h)
C 21 C 22 C 23 C 24 Z2(t - h)
C C 32 I3 0 Z3(t - h)
C 41 C 42 0 Г fKt ) 1 04" _ Z4 (t - h)
+ f2(t)
f3(t) f4 (t)
(6)
где векторы z1(t), ), z3(t), z4(t) имеют размерности т1, т2, т3, т4, а все подматрицы имеют соответствующие согласованные размеры. Аналогичным преобразованиям подвергаются начальные условия исходной системы.
Выдвигая различные предположения о соотношениях между числами т1, т2, т3, т4, а также о конкретных видах матриц В{-,С{-, получен ряд достаточных условий существования решений системы (3), единственности решений, а также соотношения, обеспечивающие совместимость начальных условий (4) для разрешимости исходной системы. Например, пусть в (6)
т3 = 0, т4 = 0. (7)
Тогда имеем систему
' ¿1 (V) = В11 г1 (V) + В12 г2 (V) + С11 z1(t — к) +
+ С12 Z2(t — к) + ),
<
0 = В21 г1 (t) + г2 (t) + С2г (t — к) +
+ С22 Z2(t — к) + /2 (),
(8)
и при выполнении условий совместимости 0 = В21ф1(0) + Ф2 (0) + С21Ф1 ( к) +
+ С22Ф2 ( к) + /2(0 (9)
данная система имеет единственное решение [8].
Аналогично могут быть записаны достаточные условия отсутствия непрерывно-дифференцируемого решения системы (6), а также требования к функции / ^) и начальные условия, при которых такие решения существуют. Но данный метод обладает тем недостатком, что для получения вышеотмеченных условий необходимо проведение достаточно громоздких преобразований. Достоинством же предложенной методики является отсутствие любых предположений типа регулярности, а также возможность перенесения полученных результатов на системы со многими запаздываниями, а также на системы с отклоняющимся аргументом нейтрального типа.
Математика
5
Основной результат. Получим явные формулы решения управляемой дескрипторной дифференциально-разностной системы с запаздыванием вида
Ш(Г) = Аг^) + А1 г^ - И) + Ви ^) (10) с начальными условиями
Нг(+0) = Нг(0) = Нг0е К", А1г(т) = Дф(т), те[-И, 0),
(11)
где Н - постоянная I х I особая матрица с гапкН = п, п < I; А, А1, В - постоянные матрицы соответствующих размеров, И > 0. Эта система является гибридной, т. е. она содержит как дифференциально-разностные, так и алгебраические уравнения.
Известно [8], что существуют невырожденные матрицы 5 и Р такие, что
SHP =
!п 0 0 0
(12)
Умножим обе части системы (10) слева на матрицу 5 и выполним замену переменной г^) = Pw(t). Получим
5ИРМ^) = SAPw(t) + 5А1РМ(( - И) + 5Ви^). (13) Введем обозначения
w(t ) =
" x(t ) " , SAP = " A11 A12
_ y (t )_ _A21 A22 _
Г A A 11 A 1 A 12 Г в, 1
SA1P = , SB = 1
A1 21 A1 Л22 _ _ B2 _
(15)
(14)
где х(() - п - вектор; у^) -1 - п - вектор. Перепишем систему (13) в виде
Х(() = А11 Х(() + А12у^) + А^ x(t - И) + + А112у^ - И) + Вхи ^), 0 = А21Х(() + А22у^) + А21 x(t - И) + + А22y(t - И) + В2и^).
Допустим, что
det А22 Ф 0. Тогда получим гибридную систему [4, 6]
х(() = А11 x(t) + А12y(t) + А111 Х(( - И) +
+ А112у(t - И) + Вхи ^), у^) = -А- А21 х(() - А- А21 Х(( - И) -. - А-2 А22y(t - И) - A-2B2u(t),
которую перепишем в нормальной форме [4, 6]
(16)
x (t ) = Au x(t ) + A12 y(t ) + A1,, x(t - h) + + A112 y(t - h) + B,u(t ), y(t ) = A21x(t ) + A121 x(t - h) + + A122y(t - h) + B2U(t).
(17)
Известно [4, 6], что система (17) при любых начальных условиях вида
x(+0) = xoe Rn, х(т) = ф(х), y(T) = у(т) т e [-h, 0)
(18)
x ( t ) =
имеет единственное решение.
Утверждение 1. Решение
= х((, 0, х0, ф, у, и), t > 0, у(() = y(t, 0, х0, ф, у, и), t > 0, системы (17), (18) существует единственно и может быть вычислено по следующей формуле (обобщенная формула Коши):
t
х(() = x(t, 0, х0, ф, у, и) = |-т)В1 +
0
+ У*(( -т)В2)и(т)^т + х((, 0, х0, ф, у, и), t >0 с начальными условиями
X*(0) = X *(0) = 1п,
t
У(() = у(, 0, х0, ф, у, и) = | (X* ( -т)В1 +
+
7y*(t - т) B2)u(T)d т + £ Z* ( jh) B2u(t)
j=0
+ У0, х0, ф, у, и> 0
с начальными условиями
X* (0) = X* (0) = А21, 2* (0) = 2* (0) = 1т,
где матрицы-функции X*(.), У*(.), 2*(.) являются решениями сопряженной системы
dX *(t )
dt
+ X *(t) A11 + Y *(t) A12 +
+ X*(t - h) A111 + Y*(t - h)A, = 0, t > 0, t ф kh; Y * (t ) = X * (t ) A21 + X *(t - h) A22 + Y *(t - h) A,, t > 0, X * (kh) - X * (kh - 0) = Z * (kh) A12, Z *(kh) = Z *((k - 1)h) A12, причем X*(t) = 0, Y*(t) = 0, Z*(t) = 0, t < 0.
Такое представление справедливо для нестационарных гибридных систем, а значит, и для нестационарных дескрипторных систем при выполнении аналога условия (16).
В качественной теории управления динамическими системами с отклоняющимся аргументом широко используется определяющее уравнение [4, 6], которое позволяет записывать решения систем в виде рядов по решениям определяющего уравнения.
Утверждение 2. Решение x(t), y(t) системы (17), (18) существует единственно и может быть вычислено по формулам
x(t) = x(t, 0, x0, Ф, v, u) = x(t, x0, ф, v, 0) +
t-iht. iU\k
+ Z Z Kjh) í (t - ) u (T)dт +
k=0 i 0 k • (t-ih> 0)
+ Z xo(ih) u(t - ih), (19)
i
(t-ih>0)
y(t) = y(t, 0, x0, ф, V, u) = y(t, x0, ф, V, 0) +
+ Z Z Xk2+i(ih)t|h(t-Tk-ih)ku(x)dт +
k=0 i 0 k • (t-ih>0)
+ Z X02(ih) u(t - ih),
i
(t-ih >0)
где x(t, x0, ф, V, 0), y(t, x0, ф, V, 0) - решение системы (17), (18) в случае нулевого управляющего воздействия: u(t) = 0, t > 0; матричные функции Xv (t), v = 1, 2, являются решением определяющих уравнений вида
Xk+1 (t) = A! X¡ (t) + A12 X¡ (t) + A¡1 X¡ (t - h) +
+A Xk2(t - h)+BUk (t),
Xk2 (t) = A21 X¡ (t) + A^ X¡ (t - h) + A^2 Xk2 (t - h) +
+ B~Uk(t), k = -1,0,1,...;t>0
с начальными условиями Xk (t) = 0, Yk (t) = 0, если k < 0, или t < 0; U0(0) = Ir, Uk (t) = 0, если
k2 + t2 * 0.
Заключение. Поскольку решения стационарных ГДР-систем вида (17) имеют рост не выше чем экспоненциальный при аналогичном росте управляющих воздействий, то аналогичное свойство имеет место и для решений дескрипторных дифференциально-разностных систем при выполнении условия (16), что позволяет применять к таким системам преобразование Лапласа.
Литература
1. Беллман, Р. Дифференциально-разностные уравнения / Р. Беллман, К. Л. Кук. - М.: Мир, 1967. - 548 с.
2. Kociecki, M. On solutions of singular equations with delays / M. Kociecki // V Polish-English on «Real-time procese control» Radziejowice. -Warsaw, 1986. - P. 227-235.
3. Campbell, S. L. Singular systems of differential equations with delays / S. L. Campbell // Applicable analysis. - 1980. - Vol. 11, № 2. - P. 129-136.
4. Марченко, В. М. Вполне регулярные системы с последействием / В. М. Марченко // Труды Ин-та математики НАН Беларуси. -2001. - Т. 7. - С. 98-105.
5. Асмыкович, И. К. О задачах качественной теории управления для дескрипторных систем с запаздыванием / И. К. Асмыкович // Труды БГТУ. Сер. VI, Физ.-мат. науки и информатика. - 1997. - Вып. V. - С. 3-14.
6. Марченко, В. М. Представление решений и относительная управляемость линейных дифференциально-алгебраических систем со многими запаздываниями / В. М. Марченко, О. Н. Поддуб-ная // Докл. РАН. - 2005. - Т. 404, № 4. - С. 465-469.
7. Марченко, В. М. К вопросу об устойчивости дескрипторных систем с запаздыванием / В. М. Марченко, И. М. Борковская, А. А. Якименко // Математика и математическое образование. Теория и практика: межвуз. сб. науч. тр. / ЯГТУ. - Ярославль, 2010. - Вып. 7. - С. 34-39.
8. Li, Y. Basic theory of singular systems of linear differential difference equations / Y. Li, Y. Liu // Preprints of the 13th World Congress IFAC, San Francisco, USA, 30th June - 5th Jule 1996 / V.D. - San Francisco, 1996. - P. 79-84.
9. Власенко, Л. А. Теоремы существования и единственности для одного неявного дифференциального уравнения с запаздыванием / Л. А. Власенко // Дифференциальные уравнения. - 2000. - Т. 36, № 5. - С. 624-628.
10. Игнатенко, В. В. Управляемость линейных дескрипторных систем с запаздыванием / В. В. Игнатенко, В. В. Крахотко // Вестник БГУ. Сер. I. - 1993. - № 3. - С. 70-73.
11. Крахотко, В. В. H-управляемость каузальных дифференциально-алгебраических динамических систем / В. В. Крахотко, Г. М. Размыслович // Вестник БГУ. Сер. I. - 2006. - № 1. - C. 123-125.
12. Размыслович, Г. П. Алгоритм вычисления передаточной матрицы для сингулярных систем с запаздыванием / Г. П. Размыслович // Вестник БГУ. Сер. I. - 1996. - № 1. - С. 52-54.
13. Щеглова, А. А. Левый регуляризирующий оператор для алгебро-дифференциальной системы с запаздыванием / А. А. Щеглова // Известия вузов. Математика. - 2003. - № 4(491). - С. 73-85.
14. Pandolfi, L. Generalized control systems, boundary control systems and delayed control systems / L. Pandolfi // Math. Control Systems. -1990. - № 3. - P. 165-181.
15. Марченко, В. М. О структуре дескрип-торных систем / В. М. Марченко // Труды БГТУ. Сер. VI, Физ.-мат. науки и информатика. -2004. - Вып. XII. - С. 3-6.
Поступила 28.02.2011
