Известия Института математики и информатики Удмуртского государственного университета
2020. Том 55. С. 102-121
УДК 517.977 © В. Е. Хартовский
ОБ ОДНОМ ЛИНЕЙНОМ АВТОНОМНОМ ДЕСКРИПТОРНОМ УРАВНЕНИИ С ДИСКРЕТНЫМ ВРЕМЕНЕМ. II. КАНОНИЧЕСКОЕ ПРЕДСТАВЛЕНИЕ И СТРУКТУРНЫЕ СВОЙСТВА
В статье изучается линейное однородное автономное дескрипторное уравнение с дискретным временем
т
Бод(к + 1) + ^ Бгд(к + 1 - г) = 0, к = т,т + 1,...,
г=1
с прямоугольными (в общем случае) матрицами Бг. Такое уравнение возникает при исследовании задач управления системами с многими соизмеримыми запаздываниями в управлении: задачи 0-управляемости, задачи синтеза регулятора типа обратной связи, обеспечивающего успокоение решения исходной системы, задачи модальной управляемости (управляемости коэффициентов характеристического квазиполинома), задачи спектральной приводимости и задачи синтеза наблюдателей для двойственной системы наблюдения. Основной метод представленного исследования базируется на замене исходного уравнения эквивалентным уравнением в «расширенном» пространстве состояний, которому сопоставили некоторый пучок матриц. Это позволило исследовать ряд структурных свойств исходного уравнения посредством использования канонической формы пучка матриц, а полученные результаты выразить в терминах минимальных индексов и элементарных делителей. В статье получен критерий существования нетривиального допустимого начального условия исходного уравнения, проверка которого основана на вычислении минимальных индексов и элементарных делителей пучка матриц. Изучена следующая задача: требуется построить решение исходного уравнения в виде д(к + 1) = Тф(к + 1), к = 1, 2 ..., где Т — некоторая матрица, последовательность векторов ф(к +1), к = 1, 2,..., удовлетворяет уравнению ф(к +1) = Бф(к), к = 1, 2,..., а квадратная матрица 5 имеет наперед заданный спектр (или часть спектра). Полученные результаты позволяют строить решения исходного дескрипторного уравнения с наперед заданными асимптотическими свойствами, например, равномерно асимптотически устойчивые.
Ключевые слова: линейное автономное дескрипторное уравнение с дискретным временем, подпространство начальных условий, представление решения, пучок матриц.
Э01: 10.35634/2226-3594-2020-56-08 Введение
Настоящая статья является второй частью работы [1], посвященной изучению линейного однородного автономного дескрипторного уравнения с дискретным временем
т
В0д(к + 1) + ^ Вгд(к + 1 - г) = 0, к = т,т +1,..., (0.1)
г=1
где Вг Е Егахг, г = 0,т, д(к) Е Ег, к = 1, 2,.... Интерес к такому уравнению возник в [2] в связи с исследованием задачи 0-управляемости, хотя в дальнейшем уравнение (0.1) было использовано и в других задачах управления линейными системами с последействием [1,3-6].
Одна из основных сложностей при изучении уравнения (0.1) заключается в описании множества начальных условий, для которых уравнение (0.1) имеет решение, а также нахождение одного из решений, порождаемого некоторым начальным условием. В работе [1]
предложен метод описания всех допустимых начальных условий, который основан на решении конечной цепочки линейных однородных алгебраических систем, а также дано описание всего множества решений, порождаемого заданным начальным условием. Основное достоинство разработанного в [1] подхода заключается в отсутствии необходимости приведения исходного уравнения (0.1) к какой-либо канонической форме.
В настоящей статье, в отличие от работы [1], на базе канонической формы матричного пучка изучаются структурные свойства этого уравнения. В частности, получено описание структуры матриц уравнения (0.1), имеющего нетривиальные решения, и предложен способ построения решений, удовлетворяющих некоторым заданным асимптотическим свойствам.
Другие подходы к исследованию дескрипторных уравнений можно найти, например, в [7-9].
§ 1. Предварительные сведения
Объект исследования — линейное автономное дескрипторное уравнение с дискретным временем
т
Вод(к + 1) + ^В*д(к + 1 - г) = 0, к = т,т +1,..., (1.1)
г=1
и начальным условием
д(г) = <7г, г = 1,т, (1.2)
где Вг е Егахг, г = 0,т, д(к) € Ег, к = 1, 2,...; дг € Ег, г = 1,т, — некоторые заданные векторы. Считаем, что среди матриц Вг, г = 0, т, по крайней мере две отличны от нулевой матрицы, причем одна из них — это матрица В0.
Определение 1.1. Под решением начальной задачи (1.1), (1.2) будем понимать бесконечную последовательность векторов д(к) € , к = 1, 2,..., удовлетворяющую соотношениям (1.1), (1.2).
Для некоторого заданного начального условия (1.2) решение начальной задачи (1.1), (1.2) может не существовать или оказаться не единственным. Набор векторов г = 1,т, для которого существует решение задачи (1.1), (1.2), будем называть допустимым набором начальных значений для уравнения (1.1). Если т = 1, то будем просто говорить о допустимом начальном векторе.
Следуя [1], рассмотрим уравнение
В0д(к + 1) + ад(к) = 0, к = т, т +1,..., (1.3)
с начальным условием
д(т) = ?с
(1.4)
где
Во
Во 0 0 L
0 0
0 0 ... Ir
Вг
В1 В2 -Ir 0
00
Вт— 1 Вт
-Ir 0
(1.5)
1Г € Кгхг — единичная матрица; д0 € Штг — заданный вектор. Решение начальной задачи (1.3), (1.4) определяется так же, как и решение начальной задачи (1.1), (1.2).
Лемма 1.1 (см. [1, лемма 1]). Если выполняется равенство
до = col [gm,..., gi],
(1.6)
0
0
то существование решения одной из начальных задач (1.1), (1.2) или (1.3), (1.4) влечет существование решения другой начальной задачи и справедлива формула
g(k) = ET(k), k = m, m +1,...,
где E = [Ir, 0,..., 0], E G Rrxmr.
Пусть начальное условие для уравнения (1.3) задается соотношением (1.4). Тогда, в силу леммы 1.1, вектор То, определяемый формулой (1.6), является допустимым начальным вектором для уравнения (1.3) в том и только том случае, когда его компоненты являются (с сохранением порядка следования) допустимым набором начальных значений для уравнения (1.1).
В силу линейности уравнения (1.3), множество его допустимых начальных векторов, которое обозначим G, образует линейное подпространство в Rmr. Если g(m) G G, то решение T(k), k = m,m + 1,..., уравнения (1.3) существует, и в силу автономности уравнения (1.3) векторы T(k), k = m, m + 1,..., принадлежат подпространству G. На основании леммы 1. 1 подпространство G дает полное описание всех допустимых наборов начальных значений для уравнения (1.1).
Пусть в подпространстве G задан некоторый базис. Введем матрицу T G Rrmxr0, где Го = dim G, в качестве столбцов которой возьмем векторы этого базиса. Тогда
G = T Rr0.
В статье [1] предложен способ построения матрицы Т, основанный на решении конечной цепочки линейных однородных алгебраических уравнений. При этом, если столбцы матрицы T образуют базис подпространства G, то, очевидно, столбцы матрицы TC, где C — произвольная невырожденная матрица, также образуют базис подпространства G. Поэтому матрица T определяется с точностью до постоянного невырожденного матричного множителя, на который матрица T умножается справа.
Замечание 1.Из определения подпространства G вытекает критерий того, что столбцы произвольной матрицы образуют базис в подпространства G. А именно, если dim G = Го, то столбцы матрицы T G Rrmxr0, rankT = Го, образуют базис в подпространстве G тогда и только тогда, когда множество допустимых начальных векторов уравнения
Во TT(k + 1)+ Bi TT(k) = 0, k = m, m + 1,..., (1.7)
совпадает с пространством Rr0.
Пусть T — любая, но фиксированная матрица, построенная указанным выше способом. Рассмотрим алгебраическое уравнение относительно неизвестной матрицы S G Rr0xr0
EoTS + BiT = 0. (1.8)
Лемма 1.2 (см. [1, лемма 2]). Матрица S, удовлетворяющая уравнению (1.8), существует.
Разобьем матрицу T на блоки T = col[Ti,.. . Tm], где матрицы Ti G R 0, i — 1, m, состоят из г строк матрицы T c номерами (i — 1)r + 1 по ir включительно, взятых с сохранением порядка их следования. Зафиксируем матрицу S — решение уравнения (1.8). Положим T = Tm, тогда равенство (1.8) можно переписать в виде системы равенств, соответствующих блочной структуре матриц £>0,
m
Ti = TSm-i, i = 1~m, BoTSm + BiTSm-i = 0.
i=1
Теорема 1.1 (см. [1, теорема 1]). Векторы i = 1,m, образуют допустимый набор начальных значений для уравнения (1.1) тогда и только тогда, когда найдется вектор c £ Rr0 такой, что
Tji = TSi-1c, i =1,m.
Теорема 1.2 (см. [1]). Пусть в начальном условии (1.2) векторы gi имеют вид 7ji = TSi-1c, i = 1,m, с некоторым c £ Rr0. Тогда любое решение начальной задачи (1.1), (1.2) представимо в виде
g(k) = T^(k) + ^(k), k = 1, 2,..., (1.9)
где последовательность ^(k), k = 1, 2,..., определяется уравнением
^(k + 1) = S^(k), k = 1,2,..., ^(1) = c, (1.10)
а последовательность ^(k), k = 1, 2,..., удовлетворяет уравнению (1.1) с нулевым начальным условием.
§ 2. Критерий существования нетривиального решения начальной задачи
Изучим вопрос существования нетривиального решения задачи (1.1), (1.2), то есть последовательности векторов g(k), k = 1, 2,..., удовлетворяющей (1.1), (1.2), у которой по крайней мере один член отличен от нулевого вектора. После этого дадим описание подпространства G, указав его размерность и формулу для вычисления одной из возможных матриц Т, столбцы которой представляют собой базисные векторы подпространства G.
Для выяснения указанных вопросов перейдем от уравнения (1.1) к уравнению (1.3). После этого невырожденным преобразованием приведем матрицы Т0, Т1 к блочно диагональному виду, диагональные блоки которых имеют «упрощенную» структуру. В силу диагональной структуры полученных после преобразования матриц, исходное уравнение (1.3) будет равносильно конечной цепочке «достаточно простых» уравнений, исследование которых не представляет трудностей.
Для того чтобы одновременно привести пару матриц {В0, к блочно диагональному виду с максимально простой структурой диагональных блоков, заменим эту пару матриц так называемым пучком матриц 5(A) = A50 + Т1 и воспользуемся возможностью приведения этого пучка матриц к каноническому виду [10, с. 331]. В связи с этим приведем некоторые понятия и факты, связанные с вопросами эквивалентности двух пучков матриц [10, с. 331].
Рассмотрим уравнение
Б(АМА) = 0 (2.1)
относительно вектора q(A). Ограничимся только теми решениями уравнения (2.1), которые, во-первых, являются полиномами, а во-вторых, среди этих решений рассматриваются только решения наименьшей степени. Под степенью решения q(A) уравнения (2.1) понимается наибольшая степень полиномов, являющихся компонентами вектора q(A).
Предположим, что имеется k нетривиальных решений этого уравнения qi(A), i = 1,k. Будем говорить, что решения qi(A), i = 1, k, линейно независимы, если
rank [qi (A),... ,qfc(A)] = k.
Напомним, что под рангом полиномиальной матрицы понимается [10, с. 335] наибольший порядок ее минора, отличного от тождественного нуля. Если же rank [q1(A),..., qk(A)] < k, то решения qi(A), i = 1,k, уравнения (2.1) будем называть линейно зависимыми. Обратим
внимание, что линейная зависимость решений ^(А), i = 1,k, предполагает существование k одновременно не равных нулю скалярных многочленов ^¿(А), i = 1,k, таких, что
EL&(АМА) = 0.
Среди всех решений уравнения (2.1) выберем ненулевое решение qi(A) минимальной степени а1. Далее выберем ненулевое решение q2 (А) уравнения (2.1) минимальной степени а2, линейно независимое с решением ^1(А). После этого выберем ненулевое решение ^з(А) уравнения (2.1) минимальной степени а3, линейно независимое с решениями ^¿(А), i = 1, 2. Будем осуществлять этот процесс до тех пор, пока это возможно. Так как число линейно независимых решений уравнения (2.1) конечно (ранг матрицы [q1 (А),..., (А)] при любом k не может превышать количества ее строк), то этот процесс должен закончиться. На последнем шаге выберем ненулевое решение дДА) уравнения (2.1) минимальной степени линейно независимое с решениями д^(А), i = 1,s — 1. Понятно, что а1 ^ ... ^ Заметим, что выбранные решения д^(А), i = 1, s, уравнения (2.1) могут определяться неоднозначно, но ряд степеней a1;...,as всегда определяется однозначно. Числа a1;...,as называются [10, с. 343] минимальными индексами для столбцов пучка матриц -В(А).
Пусть числа ß1,..., ß являются минимальными индексами для столбцов транспонированного пучка матриц (-В(А))T. В этом случае числа ß1;..., ß называются [10, с. 343] минимальными индексами для строк пучка матриц -В(А).
Заметим, что уравнение (2.1) может не иметь нетривиальных решений. В этом случае говорят, что у пучка матриц В (А) нет минимальных индексов для столбцов. Если же у пучка матриц (Б(А))T отсутствуют минимальные индексы для столбцов, то у пучка матриц В (А) не будет минимальных индексов для строк.
Напомним, что бесконечными элементарными делителями пучка матриц В (А) называются [10, с. 333] элементарные делители вида матрицы АВ0 + ßB1. Под конечными элементарными делителями пучка матриц В (А) понимаются элементарные делители матрицы В (А).
Два матричных пучка В (А) и В (А) называются строго эквивалентными, если существуют не зависящие от переменной А невырожденные матрицы Г1; Г2 такие, что В (А) = Г1 -0(А)Г2. Известно [10, с. 344], что два пучка матриц строго эквивалентны тогда и только тогда, когда эти пучки матриц имеют одни и те же минимальные индексы для строк и для столбцов и одни и те же элементарные конечные и бесконечные элементарные делители.
Пусть пучок матриц В (А) имеет минимальные индексы для столбцов
а1 = ••• = аТ1 = 0, аТ1+1,..., и минимальные индексы для строк
ß = ••• = ßT2 =0, ßT2 + 1,...,ßl.
Если равные нулю минимальные индексы для столбцов или для строк отсутствуют, то считаем, что т1 = 0 или т2 = 0 соответственно. Обозначим i = 1,p, — степени всех бесконечных элементарных делителей вида пучка матриц -В(А). Рассмотрим матрицу K(А) размера (r(m — 1) + n) х mr (то есть такого же размера, что и матрица -В(А)), которая имеет вид
K(А) = diag [0,L«ti+i(А),.. .,Las(А), (L^+i(А))т,
..., (Lei (А))т ,N(^(А),..., N)(А), J(А)] , (2.2)
где нулевой блок, стоящий на первом месте матрицы в правой части (2.2), имеет размер т2 х ть матрица Ьг(Л) размера г х (г + 1) имеет вид
ад
Л 10 ... 0 0 0 Л 1 ... 0 0
0 0 0 . . . Л 1
1,2,...;
N «(Л) = ЛЯ(i) + H(i) — матрица размера i х i, у которой элементы первой наддиа-гонали равны единице, остальные элементы равны нулю; пучок матриц J(Л) имеет вид J(Л) = Л/0 + J, где J0, J — произвольные матрицы размера $ х но такие, что | J0| = 0, а наборы конечных элементарных делителей пучков матриц В (Л) и J(Л) совпадают. Матрицы ¿¿(Л) и N(i) (Л) не имеют конечных элементарных делителей, бесконечные элементарные делители имеют только блоки вида N(i)(Л). Матрица К(Л) называется канонической формой пучка матриц В (Л).
Замечание 2. В большинстве случае запись вида diag [Ai,..., A¿] используется для обозначения блочно диагональной матрицы с квадратными диагональными блоками A¿ (или просто диагональной матрицы, если A¿ — некоторые скалярные величины). Однако в настоящей работе обозначение diag [Ai,..., A¿] указывает на матрицу с блочно диагональной структурой, диагональные блоки A¿ которой в общем случае могут быть прямоугольными.
Замечание 3. Если тг = 0, а т2 = 0, то в матрице (2.2) первые т2 строк равны нулю, оставшиеся строки образуют блочно диагональную матрицу виду (2.2). Если же Ti = 0, а Т2 = 0, то в матрице (2.2) первые Ti столбцов равны нулю, а оставшиеся столбцы также образуют блочно диагональную матрицу вида (2.2).
Известно [10, с. 342], что для любого пучка матриц В (Л) существуют невырожденные матрицы Гь Г2, приводящие этот пучок матриц к каноническому виду
ri Б(Л)Г2 = K (Л). (2.3)
Пусть имеет место (2.3). Тогда
Г i Б0Г2 = Ко, Г i Bi Г2 = Ki, (2.4) где Ki = К (0), К0 = К (1) — Ki (то есть К (Л) = К0Л + Ki). Введем матрицы
L¿i = L¿(0), ¿¿о = L¿(1) — L¿i, (2.5) (то есть ¿¿(Л) = + Lii) и рассмотрим уравнения
(Ьго)Т g(k + 1) + (L¿i)T g(k) = 0, k = m,m +1,..., (2.6)
H (i)g(k + 1)+ g(k) = 0, k = m,m +1,.... (2.7)
Л е м м а 2.1. Пространства решений уравнений (2.6) и (2.7) состоят только из нулевых векторов.
Доказательство. Рассмотрим уравнение (2.6). Представим векторы g(k), k = 1, 2, ..., в виде g(k) = [gi(k),... ,g¿(k)]T. Тогда уравнение (2.6) в координатной форме будет иметь вид
gi(k + 1) = 0, g2(k + 1)+ gi(k) = 0, ..., (28) gi(k + 1) + gi-i(k) = 0, gi(k) = 0, k = 1, 2,.... ( . )
Из предпоследнего уравнения в (2.8) c учетом равенства g¿(k + 1) = 0, k = 1, 2,..., получаем, что gi-i(k) = 0, k = 1, 2,.... Продолжая эти рассуждения от (i — 1)-го уравнения до второго уравнения системы (2.8), получим, что gj(k) = 0, j = 1, i, k =1, 2,....
Если теперь переписать в координатной форме уравнение (2.7), то получим соотношения (2.8), только без первого уравнения. Проведя такие же рассуждения, приходим к завершению доказательства леммы. □
Замечание 4. Лемму 2.1 можно также доказать, если строить базис пространства решений из уравнений (2.6) и (2.7) по методу статьи [1].
Теорема 2.1. Для того чтобы существовало нетривиальное решение уравнения (1.1) необходимо и достаточно, чтобы выполнялось по крайней мере одно из следующих условий:
(a) пучок В(\) имеет минимальные индексы для столбцов;
(b) пучок В(\) имеет конечные элементарные делители.
Доказательство. Пусть матрицы Г1; Г2 реализуют равенство (2.3). В уравнении (1.3) сделаем замену g(k) = r2jj(k), k = m,m +1,..., и умножим полученное соотношение на матрицу Г1 слева. Тогда последовательность векторов jj(k), k = m, m + 1,..., в силу (2.4) удовлетворяет уравнению
K0jj(k + 1) + K1g(k) = 0, k = m,m +1,.... (2.9)
Очевидно, что существование нетривиального решения уравнения (1.3) равносильно существованию нетривиального решения уравнения (2.9), поэтому далее ограничимся исследованием уравнения (2.9). В векторе jj(k) выделим векторные компоненты, соответствующие блокам матрицы (2.2):
j(k) = col [gri (k), (k).., (k), S^+i(k),
(k),SYi (k),...,SYp (k),jJ (k)], (2.10)
k = m, m + 1 , . . . ,
где jy» G RT1, За, (k) G RTi+1, i = Ti + 1,s, je (k) G Re, i = T2 + 1,1, jj7i (k) G RYi, i = 1,p, jj(k) G R^. Для нахождения векторных компонент (2.10) заменим уравнение (2.9) в соответствии с блочной структурой матриц K0, K1 следующей системой уравнений
Laiogai (к + 1) + 1д« (к) = 0, г = Т1 + 1,5, к = т,т +1,..., (2.11)
(Ьдо)тдд(к + 1) + (^а)Тдв,(к) = 0, г = т2 + 1,1, к = т,т +1,..., (2.12)
Ны д7, (к + 1)+ дъ (к) = 0, г = Т7р, к = т,т +1,..., (2.13)
^ (к + 1) + Лд./ (к) = 0, к = т,т +1,.... (2.14)
Отметим, что в качестве векторной компоненты дТ1 (к), к = т,т + 1,..., (если Т1 = 0) можно взять любые векторы подходящего размера.
Проанализируем уравнения (2.11 )—(2.14). В силу леммы 2.1 пространства решений уравнений (2.12) и (2.13) состоят только из нулевых векторов. Поэтому уравнения (2.12) и (2.13) имеют только тривиальное решение при нулевых начальных условиях д^ (т) = 0, г = т2 + 1,/, д7,(т) = 0, г = 1,р. Уравнения (2.11) и (2.14) разрешимы при любых начальных условиях д«(т), г = т1 + 1,5, д^/(т), поскольку ранги каждой из матриц ¿а,0, г = Т1 + 1,5, и равны количеству ее строк. Действительно, в данном случае у каждой из матриц ¿а,0, г = т1 + 1,5, и всегда можно выбрать базисный минор, порядок которого равен количеству строк соответствующей матрицы. Поэтому каковы бы ни были
векторы да. (к), г = Т1 + 1,5, д.;(к), к = т, т + 1,..., всегда найдутся векторы (к + 1), г = Т1 + 1,5, 'gJ(к + 1), являющиеся решениями алгебраических уравнений (2.11) и (2.14) относительно неизвестных векторов (к + 1), г = Т1 + 1,5, д/(к + 1). Далее для уравнений (2.11)—(2.14) будем рассматривать только допустимые начальные условия, причем считаем, что для уравнений (2.11) и (2.14) начальные условия являются ненулевыми.
Пусть выполняется условие (а) теоремы 2.1. Тогда в матрице (2.2) присутствуют нулевые столбцы (если Т1 = 0) или блоки вида Ьа. (Л) (или и то, и другое). Если выполняется условие (Ь), то в матрице (2.2) имеется блок 3(Л).
Если Т1 = 0, то одно из решений д(к), к = т, т + 1,..., уравнения (2.9) можно получить, если в представлении (2.10) взять в качестве компонент дТ1 (к), к = т, т + 1,..., произвольные ненулевые векторы, а остальные векторные компоненты в (2.10) положить равными нулю. Если в матрице К (Л) присутствуют блоки (Л), г = Т1 + 1,5, то каждому из них будет соответствовать нетривиальное решение д(к), к = т, т + 1,..., уравнения (2.9), у которого в представлении (2.10) векторные компоненты (к), к = т, т + 1,..., есть решение уравнения (2.11), а остальные векторные компоненты равны нулю. Если в матрице К (Л) присутствует блок 3 (Л), то в качестве нетривиального решения д(к), к = т, т +1, ..., уравнения (2.9) можно взять векторы, у которых в представлении (2.10) векторные компоненты (к), к = т,т + 1,..., есть решение уравнения (2.14), а остальные равны нулю.
Теперь предположим, что ни одно из условий (а) и (Ь) теоремы 2.1 не выполняются. Тогда в матрице (2.2) нет нулевых столбцов и блоков вида (Л), 3(Л). В этом случае уравнение (2.9) равносильно уравнениям (2.12), (2.13), которые могут иметь только тривиальные решения при нулевых начальных условиях. Теорема доказана. □
Из представления уравнения (2.9) в виде (2.11)—(2.14) следует, что уравнение (1.3) разрешимо при любом начальном условии (1.4) в том и только том случае, когда в канонической форме К (Л) отсутствуют блоки вида (¿¿(Л))т и N (г)(Л), что равносильно отсутствию ненулевых минимальных индексов для строк и бесконечных элементарных делителей пучка матриц. Действительно, в силу блочной структуры матрицы К (Л), уравнение (1.3) равносильно в этом случае системе уравнений (2.11), (2.14). Сформулируем это утверждение в терминах, не требующих приведения пучка матриц к каноническому виду.
Следствие 2.1. Множество допустимых векторов уравнения (1.3) совпадает со всем пространством Егт (то есть О = Егт) тогда и только тогда, когда пучок матриц В(Л) не имеет ненулевых минимальных индексов для строк и бесконечных элементарных делителей.
Замечание 5. Следствие 2.1 дает идею для построения одной из возможных матриц Т в случае, когда пучок матриц В(Л) имеет каноническую форму (2.2), то есть В(Л) = К (Л). В соответствии с замечанием 1 уравнение
КоТс(к + 1)+ К1Тс(к + 1) = 0, к = т,т + 1,..., (2.15)
должно быть разрешимо при любом векторе с(т). То есть уравнение (2.15) должно быть равносильно системе уравнений вида (2.11) и (2.14). Поэтому матрицу Т можно выбрать так, чтобы в результате умножения матриц К (Л) и Т образовалась матрица, которая получается из матрицы К (Л) путем удаления у нее столбцов, содержащих элементы блоков вида (¿¿(А))т и N(г)(А). Перенесем эту идею на общий случай и дадим точную формулировку.
Рассмотрим соотношение (2.3), приводящее пучок матриц В(Л) к каноническому виду (2.2). Составим матрицу Тг вида
°$х6>1 ^
т =
г
где 0iXj обозначает нулевой блок размера i х j, 01 = т1 + s + aT1+1 + ... + — сумма количеств столбцов всех матриц вида L¿(Л), 02 = eT2+i + ... + А + Y1 + ... + YP — сумма количеств столбцов всех матриц вида (¿¿(Л))Т и N(г)(Л).
Следствие 2.2. В качестве матрицы Т можно взять матрицу Т = Г2Тг.
Доказательство. Матрица К (Л)Тг получается из матрицы К (Л), если в последней удалить столбцы, содержащие элементы блоков (L^(Л))Т, i = т2 + 1,/, и N(7г)(Л), i = 1,p. Поэтому, согласно следствию 2.1, уравнение
КоТг c(k + 1) + KiTr c(k) = 0, k = m,m +1,..., (2.16)
разрешимо при любом начальном векторе c(m) G Rrm. На основании формул (2.4) уравнение (2.16) запишем в виде
Г i ВоГ Тг c(k + 1) + Г1В1 Г2Т c(k) = 0, k = m, m + 1,...
Умножив слева полученное соотношение на матрицу Г-1, приходим к уравнению (1.7), которое будет разрешимо при любом начальном векторе c(m) G Rrm. Обращение к замечанию 1 завершает доказательство. □
Размерность пространства G совпадает с количеством столбцов матрицы Т (или, что то же самое, с количеством столбцов матрицы Тг). Поэтому можно сформулировать следующее утверждение.
Следствие2.3. Размерность го подпространства G есть разность между количеством столбцов пучка матриц В (Л) и суммой всех его ненулевых индексов для строк и всех его степеней бесконечных элементарных делителей: го = mr — 02 (или, что то же самое, Го = 0! +
Введем матрицу В(Л) = Y1 ™о Для того чтобы при вычислении минимальных
индексов пучка матриц В(Л) уменьшить размерности решаемых уравнений, представляется полезным следующее утверждение.
Теорема 2.2. 1) Любое решение д(Л) уравнения (2.1) представимо в виде д(Л) = = col [Лт-1/г, Лт-21г,..., Ir]д(Л), где д(Л) = col[ífi(Л),..., (Л)] — векторный полином, удовлетворяющий уравнению
В (Л)д(Л) = 0.
2) Любое решение д(Л) уравнения д(Л)В(Л) = 0 представимо в виде д(Л) = £(Л) [/„, ВоЛ + Bi, Во Л2 + В1Л + В2,..., Во л™-1 + В^-2 + ■ ■ ■ + Bm_i], где д(Л) = [д1 (Л),..., #П(Л)] — векторный полином, удовлетворяющий уравнению
£(Л)В(Л) = 0.
Доказательство. Докажем первую часть теоремы 2.2. Вектор решения д(Л) уравнения (2.1) представим в виде д(Л) = со1[д1(Л),..., дт(Л)], где %(Л) — r-векторные многочлены. Тогда, в силу блочной структуры матриц (1.5), уравнение (2.1) равносильно m уравнениям ( )
(ЛВо + В1) ?1(Л) + В2^2(Л) + ■ ■ ■ + Втдт(Л) = 0,
Л?2(Л) — 51(Л) = °, (2 17)
Лдт(Л) — дт_1(Л) = 0.
Положим дт(А) = д(А). Тогда последнее уравнение в (2.17) равносильно тому, что дт-1(А) = Ад(А). Предпоследнее уравнение в (2.17) равносильно тому, что дт-2(А) = = Адт-1(А) = А2д(А), и так далее. Продолжая эти рассуждения получим, что последние т — 1 уравнений системы (2.17) равносильны соотношениям дт-ДА) = Агд(А), г = 0,т — 1. Подставляя дт-ДА), г = 0, т, в первое уравнение системы (2.17), приходим к утверждению теоремы.
Перейдем к доказательству второй части теоремы 2.2. Уравнение д(А)5(А) = 0 запишем в виде системы уравнений, советующих блочной структуре матрицы -В(А):
#1 (Д + АД) = #2, ^1^2 + А#2 = дз, ..., д1^т-1 + Адт-1 = дт, д1Вт + Адт = 0,
где д^ — векторные компоненты вектора д, то есть д = [д1,..., дт]. Далее, применив к выписанной системе уравнений рассуждения такого же типа, что и при доказательстве первой части теоремы, придем к справедливости доказываемого утверждения. □
§ 3. Собственные значения матрицы Б
На протяжении настоящего параграфа считаем, что зафиксирован произвольный базис подпространства О, на базе которого построена некоторая матрица Т. Рассмотрим формулы (1.9), (1.10). Асимптотическое поведение последовательности векторов ^(к), к = 1, 2,..., зависит от собственных значений матрицы Б, найденной из уравнения (1.8). В общем случае матрица Б может определяться неоднозначно. Соответственно и решение начальной задачи (1.1), (1.2) также будет не единственным. Однако при решении некоторых задач управления и наблюдения динамическими системами с запаздыванием матрица Б может определять структурные свойства соответствующих регуляторов и наблюдателей [3-6]. Поэтому возникает задача выбора матрицы Б, позволяющей получить одно из решений начальной задачи (1.1), (1.2), которое удовлетворяет некоторому заданному асимптотическому свойству, зависящему от собственных значений матрицы Б. Ниже приводятся результаты исследования этого вопроса.
Если матрица Т равна нулю, то Б — любая матрица размера г0 х г0 ив дальнейшем исследовании этого случая необходимости нет. Поэтому далее считаем, что Т = 0.
Рассмотрим пару матриц [В0Т, ДТ), в соответствие которой поставим пучок матриц
% (А) = АДТ + ДТ.
Лемма 3.1. Пусть матрица Т не равна нулю. Тогда пучок матриц -Тт(А) строго эквивалентен пучку матриц
К (А) = ^ [0, Др (А),..., (А), 7(А)], (3.1)
где размеры нулевого блока, стоящего на первом месте матрицы в правой части (3.1), соответствуют количеству нулевых минимальных индексов пучка Др(А); вид матриц Д(А) описан ранее; 7(А) = Ад + Д, матрицы Д г = 0,1, являются квадратными, причем | д = 0; числа г = + — минимальные не равные нулю индексы для столбцов пучка матриц Д^(А) (некоторые из блоков в представлении (3.1) могут отсутствовать).
Доказательство. Согласно [10, с. 342], существуют невырожденные матрицы Г3, Г4, приводящие пучок матриц -Тт(А) к каноническому виду: Г3Др(А)Г4 = К(А). Матрица К(А) в общем случае имеет вид
К (А) = ^[0,ДР (А),..., ^
Ч+„(А), (Д»(А))т,..., (Д»+г(А))т, N(^(А),..., N(^(А), 7(А)], .
где г = т, и> + г, — минимальные не равные нулю индексы для строк пучка (Л), матрицы N(г)(Л) описаны ранее, г =1,5, — степени всех бесконечных элементарных делителей вида ^ пучка Д^(Л), остальные обозначения указаны при описании формулы (3.1). Предположим, что множество индексов ,..., , £ъ • • •, не является пустым.
Рассмотрим уравнение (1.8). Согласно лемме 1.2 это уравнение всегда имеет решение. Пусть Б — любое фиксированное решение уравнения (1.8). Обозначим
Б = Г-1БГ4 (3.3)
и рассмотрим упорядоченный набор квадратных диагональных блоков матрицы Б
(Бо, БйР, • • • , Бйр+^, БС», • • • , БС»+*, Б?1, • • • , %, Б1) , (3.4)
соответствующих блокам матрицы (3.2):
0, ч(Л), • • •,Ь6р+Л(Л), (Л))т, • • •, (Дго+г(Л))т,N^(л), • • •, N^(Л), 7(Л)
(то есть размеры блоков (3.4) совпадают с количествами столбцов диагональных блоков матрицы (3.2)). Заметим, что матрица Б в общем случае может не иметь блочно диагональной структуры. Однако блоки матрицы Б, отличные от диагональных, для дальнейшего доказательства не нужны. Обозначим
К = К(0), К = К(1) - К (3.5)
(то есть К (Л) = ЛК (0) + К (1)). Из существования матрицы Б следует, что имеет место равенство
КБо Б + К = 0. (3.6)
Матрицы К0 и К 1 имеют блочно диагональную структуру. Поэтому из соотношения (3.6) следует система равенств (заметим, что обратная импликация в общем случае неверна)
Д; oSSi + Да = 0, i = + (LCio) T/Ci + (Lca)T = 0, i = + Z (37)
H+ Д = 0, i = 1, s, JoSj + Ji = 0,
где матрицы L¿i0, Дд в формулах (3.7) определены соотношениями (2.5), остальные обозначения описаны после формулы (2.2). Однако ни при каких матрицах равенства вида
(Lio)TSi + (Lii)T = 0, i G {Zw,... , Zw+z}, H+ = 0, i G {&,..., £J,
в (3.7) выполняться не могут, поскольку у матриц (Д0)T и H(i) последние строки нулевые, а у матриц (L¿1)T и I последние строки нулевыми не являются. Получили противоречие с утверждением о том, что матрица S/ а вместе с ней и матрица S, существуют. Поэтому из существования матрицы S следует, что множество индексов (Zw,..., (w+z, ..., является пустым. А это значит, что в канонической форме пучка КС (Л) отсутствуют блоки вида (¿¿(Л))1,, N(г)(Л). То есть каноническая форма пучка Др(Л) имеет вид (3.1). Лемма 3.1 доказана. □
Замечание 6. Обратим внимание, что доказательство леммы 3.1 было основано на существовании решения матричного уравнения (1.8). Однако лемму 3.1 можно доказать иным способом, опираясь на теорему 2.1 и следствие 2.2. Для этого необходимо привести пучок матриц B(Л) к каноническому виду K(Л), определяемому формулой (2.2), затем умножить полученную матрицу K(Л) на матрицу Тг. После этого, на основании доказательства следствия 2.2, остается только переставить строки полученной матрицы K(Л)Тг так, чтобы результатом этого преобразования
стала матрица такого же вида, что и матрица К (Л). В этом случае будут выполняться равенства aTl+i = ¿р,..., а = ¿p+d- Кроме этого, пучки матриц J (Л) и ./(Л) будут строго эквивалентными.
Замечание 7. Из леммы 3.1 следует, что если T = 0, а матричный пучок Bf (Л) имеет минимальные индексы для строк, то эти индексы равны нулю. Бесконечных элементарных делителей матричный пучок Bf (Л) иметь не может.
Из доказательства следствия 2.2 понятно, что умножение матрицы K(Л) на матрицу Тг «удаляет» у матрицы K(Л) столбцы, содержащие элементы блоков (¿¿(Л)) и Мг)(Л). Поэтому можно сделать следующий вывод.
Следствие3.1. Матрицы В(Л) и B,f (Л) имеют одни и те же минимальные индексы для строк и одни и те же конечные элементарные делители.
Далее по-прежнему считаем, что матрица T = 0. Пусть столбцы матрицы S0 представляют собой фундаментальную систему решений однородного уравнения B0Tx = 0 (см. уравнение (1.8)). То есть B0TS0 = 0, матрица S0 имеет максимально возможный ранг, и столбцы этой матрицы линейно независимы. Тогда множество всех матриц S, удовлетворяющих уравнению (1.8), определяется равенством
S = S0C + S *, (3.8)
где квадратная матрица S* есть любое решение уравнения (1.8), C G Rs°xr° (s0 = rank S0) — произвольная матрица.
Обозначим r1 = rank В^(Л). Заметим, что r1 ^ r0 (напомним, что r0 — количество столбцов матрицы В^(Л)). Введем в рассмотрение множество Л1, состоящее из таких чисел Л. G C, для которых ранг числовой матрицы B,f (Л.) меньше числа ее столбцов:
Л1 = {Л. G C: rankВ^(Л,) < r0}.
Если r1 < r0, то множество Л1 совпадает с множеством C, то есть состоит из бесконечного числа элементов. Такая ситуация может быть, если в матрице B,f (Л) имеются нулевые столбцы или блоки вида ¿(Л) (заметим, что матрица (Л) имеет размер i х (i + 1), а rank ¿¿(Л) = i). А последнее возможно, если пучок матриц B,f (Л) имеет минимальные индексы для столбцов (напомним, что в силу леммы 3.1 матрица В^(Л) не имеет ненулевых минимальных индексов для строк и бесконечных элементарных делителей).
Пусть ^¿(Л), i = 1,7, — система всех конечных элементарных делителей матрицы B,f (Л) и
¥>(Л) = П (Л),
¿=0
а ^ — степень полинома <^(А), ^ = deg<^(А). Если матрица £>т(А) не имеет конечных элементарных делителей, то полагаем <^(А) = 1. Предположим, что пучок матриц £>т(А) имеет минимальные индексы для столбцов
¿1 = ... = ¿р-1 = 0, ¿р, ..., ¿р+^
(если у пучка матриц Др(А) нет нулевых минимальных индексов для столбцов, то считаем, что р = 1).
Обозначим ф3 (А) — характеристический полином матрицы найденной из уравнения (1.8).
Теорема3.1. Пусть Т = 0. Тогда имеют место следующие утверждения.
(a) Множество различных собственных значений всевозможных матриц S, найденных из уравнения (1.8), совпадает с множеством Аь
(b) Для характеристического полинома любой матрицы S, найденной из уравнения (1.8), справедливо равенство
фд (A) = ^(А)п(А), (3.9)
где п(А) — некоторый приведенный полином, deg п(А) = го — d^.
(c) Для любых приведенных полиномов с действительными коэффициентами v0(A), vp+j(A), i = 0,d, deg v0(A) = p — 1, degv¿(A) = #p+i + 1, i = 0,d, существует такое решение S уравнения (1.8), что
ф3 (Л) = p(A)v (Л), (3.10)
где v(A) = vo(A)vp(A)... vp+d(A) и deg v(A) = p + d + óp + ... + íp+d = Го — d^. 5 частности, для любого полинома v(A), deg v(A) = r0 — d^ с действительными коэффициентами, имеющего только действительные корни, существует такое решение S уравнения (1.8), что имеет место (3.10).
Доказательство. 1) Вначале докажем первое утверждение (a) теоремы 3.1. Пусть A0 — собственное значение произвольной фиксированной матрицы S, найденной из уравнения (1.8), а q £ Rro — собственный вектор матрицы S, отвечающий собственному значению A0. Тогда имеет место соотношение
(A0/ro — S) q = 0. (3.11)
Умножим обе части равенства (3.11) слева на матрицу В0ГГ. В силу уравнения (1.8) справедливо равенство B0TS = — B^, поэтому можем записать
(AB0? + BiT)q = 0. (3.12)
Равенство (3.12) по другому можно записать в виде B^(A)q = 0. Следовательно A0 £ Л1. Значит все собственные значения любой матрицы S, найденной из уравнения (1.8), являются элементами множества Л1.
Покажем, теперь, что для любого числа A £ Л1 найдется матрица S, удовлетворяющая уравнению (1.8) такая, что число A есть собственное значение матрицы S. Поскольку A £ Ль то найдется ненулевой вектор q £ Rro, удовлетворяющий алгебраической системе B^(A)q = 0, то есть системе (3.12). Пусть S* — любое решение уравнения (1.8), тогда для него справедливо равенство -B0TS* = —BB1T. На основании этой формулы соотношение (3.12) можно переписать в виде B0T(A/ro — S*)q = 0. Значит вектор (A/ro — S*) q есть линейная комбинация столбцов матрицы S0,
(Airo — S*)q = S0C0, (3.13)
где c0 £ Rso — некоторый постоянный вектор. Для любого вектора c0 £ Rso найдется матрица C £ Ms°xr° такая, что c0 = Cq. Поэтому равенство (3.13) можно переписать в виде
(iroA — S0C — S*) q = 0. (3.14)
Соотношение (3.14) в силу (3.8) говорит о том, что число A есть собственное значение некоторой матрицы S, найденной из уравнения (1.8).
2) Перейдем к доказательству части (Ь) теоремы. Пусть Б — любая матрица, удовлетворяющая уравнению (1.8). Поскольку преобразование подобия не меняет характеристического полинома матрицы, рассуждения будем вести применительно к матрице Б?, определяемой формулой (3.3). В этом случае вместо уравнения (1.8) будем использовать уравнение (3.6), матрицы которого определяются формулами (3.5). В силу того что матрица Т = 0, пучок матриц (Л) строго эквивалентен пучку матриц К(Л), который имеет вид (3.1). При этом упорядоченный набор (3.4) диагональных блоков матрицы Б?, описанный в доказательстве леммы 3.1, будет следующим:
(Бо ,Бйр,... ,Бйр+^, Б!). Представим матрицы Б?, К0, К? 1 следующим образом
ББ
Б?11 Б-Б21 Б? 1
12
Ко
КБ01 0 0 Б
К1
КБ11 0 0 /
(3.15)
(3.16)
где Бу, К01, К:1 — некоторые блоки подходящих размеров, блок Б?/ тот же, что и в (3.15), блоки Б0, определяют в формуле (3.1) матрицу /(Л) = ЛБ0 + / Напомним, что в силу леммы 3.1 |Б0| = 0 (а значит, и | Б(Л)| ф 0). Используя (3.16), выполним перемножение
матриц в соотношении (3.6), откуда получим, что Б0Б;
21
0, ЛоБт + Б
0. Значит все
элементы матрицы Б?, расположенные в тех же строках, что и блок Б?у, но левее его, равны нулю, то есть Б21 = 0, а блок Б/ определяется однозначно формулой Б/ = — Б—1 Бъ Поэтому матрица Б? имеет вид
Б11 Б12 0 — б—Б
б =
Все конечные элементарные делители матрицы К?(Л) определяются конечными элементарными делителями матрицы Б(Л) = Б0(ЛД + Б-^) . Значит для характеристического полинома матрицы Б выполняется равенство
^¿?(Л) = |Л1# + 1 б?:1 |Л1го- # — Б11|
^(Л)|Л/Г0_^ — 'и1
(д — размер блока Б (Л)). Отсюда имеем равенство (3.9) при п(Л) = |Л/Го_^ — Б?11|. Кроме того видно, что ^ = д. А так как из канонической формы (3.1) следует, что г0 = д + р +
+ й + ¿р + ... + , то deg V(Л) = Го — д = Го — ^ = р + й + ¿р + ... + ¿2
'р+й-
3) Докажем часть (с) теоремы. Зададим произвольный приведенный полином V(Л), deg V(Л) = г0 — й^. Покажем теперь, как выбрать блоки матрицы Б так, чтобы выполнялось равенство (3.10).
Введем матрицы Ф^(с^0,..., сц) € Е(г+1)х(г+1), г = 1, 2,..., формулой
Фг(Со, . . . ,Сг) =
0 —1
00
0
Со
0
Сг1
0 —1
, г = 1, 2,
где числа су € К определяются ниже. Для того чтобы выбрать матрицу Б?, для которой имеет место равенство (3.10), полагаем равными нулю все блоки этой матрицы, за исключением диагональных. Легко увидеть, что матрица Б?, диагональные блоки (3.15) которой имеют вид
Бо = Со, = (с^о,... ), г = р,р + й, Б/ = — Б0 Б:,
(3.17)
с
с
где C0 — произвольная матрица соответствующего размера, , j = i = p,p + d, — любые действительные числа, удовлетворяет уравнению (3.6). При указанном выборе диагональных блоков (3.15) характеристический полином матрицы S будет иметь вид
^(А) = |A/p-i - Co| + -Т—^П |A/¿i+i - Фг.(c¿ii,..., c¿i)|,
¿=p
где (p — 1) — размер блока S0. Но | A/# + = <^(A). Поэтому остается только выбрать
подходящие матрицу C0 и числа cj так, чтобы выполнялось равенство
p+d / ¿i+1 \ v(А) = | A/p-i — C0| П Aái+1 + 1)jc*íi-jAái+1-j . (3.18)
¿=p V j=i /
Для этого достаточно обеспечить выполнение следующих соотношений
íp+i + 1
vo(A) = |A/p-i — C0|, Vp+i(A) = Aíp+i+1 + £ (—1)jcíp+iAíp+i+1-j, i = 0Д
j=i
что не представляет труда. Если же равенство (3.10) следует обеспечить для любого приведенного полинома v(A) с действительными коэффициентами, то это можно сделать без труда, воспользовавшись равенством (3.18). Теорема 2.1 доказана. □
Напомним следующие понятия [10, с. 332]. Пучок матриц B,f (A) называется регулярным, если матрицы B¿T, i = 0,1, являются квадратными матрицами и |B,f(A)| ф 0. Во всех остальных случаях (матрицы B0T, BlT, i = 0,1, не являются квадратными или |B,f (A) | ф 0) пучок матриц B,f (A) называется сингулярным.
Из доказательства теоремы 3.1 можно сделать следующие выводы.
Следствие 3.2. Пусть матрица T = 0. Если у пучка матриц B,f (A) отсутствуют минимальные индексы для столбцов, то матрица S определяется уравнением (1.8) однозначно.
Доказательство. Действительно, при выполнении условий следствия 3.2 каноническая форма пучка матриц B,f(A) имеет вид К (A) = col [0, J(A)], где J(A) — регулярный пучок матриц, J(A) = A/ + / Поскольку в силу леммы 3.1 у пучка B,f(A) отсутствуют бесконечные элементарные делители, то | J01 = 0. Поэтому матрицу S/ определяемую уравнением (3.3), можно однозначно найти из уравнения J0S + ./^l = 0. □
Следствие 3.3. Если пучок матриц Bf(A) является регулярным, то матрица S определяется уравнением (1.8) однозначно.
Доказательство. Из разрешимости уравнения (1.7) имеет место следующая импликация: если для некоторого вектора g выполняется равенство gTB0T = 0, то gTBlT = 0. Но в силу регулярности пучка матриц B,f(A) равенство gT B0T = 0 возможно только для нулевого вектора g. Действительно, если существует ненулевой вектор g такой, что gTB0T = 0, то в силу равенства gTBlT = 0 строки определителя |B,f (A)| линейно зависимы, что противоречит условию регулярности.
Поскольку равенство gTB0T = 0 возможно только для нулевого вектора g, то |B0T| = 0, и матрицу S можно определить по формуле S = — (B0T) 1B?1TT. Следствие доказано. □
Рассмотрим формулы (1.9), (1.10). Теорема 3.1 говорит о том, как в ряде случаев можно построить одно из решений начальной задачи (1.1), (1.2), удовлетворяющее некоторым заданным асимптотическим свойствам. Остановимся на этом подробнее.
Предположим имеется необходимость построить равномерно асимптотически устойчивое решение #(&) начальной задачи (1.1), (1.2). Существование такого решения в силу формул (1.9), (1.10) равносильно существованию матрицы Б, собственные значения которой меньше единицы. Проанализируем, опираясь на теорему 3.1, возможные случаи исследования этого вопроса.
1) Собственные значения любой матрицы Б принадлежат множеству Л:. Предположим, что число элементов множества Л: конечно. Это значит, что пучок матриц -Е>^(Л) не имеет минимальных индексов для столбцов и в канонической форме (3.1) имеется только блок Б(Л) (и, возможно, нулевые строки). Тогда, если все элементы множества Л: (или корни полинома <^(Л), что в случае конечной мощности множества Л: тоже самое) по модулю меньше единицы, то нужная матрица Б существует. В данном случае ее можно найти как единственное решение уравнения (1.8). Если же среди элементов матрицы Л: (или корней полинома <^(Л)) есть больше или равные по модулю единицы, то соответствующая матрица Б не существует.
2) Пусть, теперь, число элементов множества Л: бесконечно. Множество собственных значений любой матрицы Б содержит корни полинома <^(Л). Поэтому если среди корней полинома <^(Л) есть такие, модуль которых больше или равен единицы, то подходящая матрица Б также не существует. А если все корни полинома <^(Л) по модулю меньше единицы, то нужная матрица Б существует. Найти такую матрицу Б можно, например, опираясь на доказательство теоремы 3.1.
Проиллюстрируем сказанное примером.
Пример 3.1. Пусть в уравнении (1.1) т =1, а его матрицы имеют вид
Во
0000 0000 0100 0 0 0 1
-1
0 0 0 0
0 0 0 0
0 0 1 0
0 0 0 — 2
Требуется построить матрицу Б, собственные значения которой по модулю меньше единицы. В данном случае матрица Т = 14, матричные пучки В (Л) и —т (Л) равны: -В(Л) = ¿^(Л),
вБ(Л)
0 0 0 0 0 0 0 0 0 Л 1 0
000 л — 2
г0 = 4, г: = 2, множество Л: состоит из бесконечного числа элементов, система элементарных делителей имеет вид ^(Л) = Л — ^ и <^(Л) = Л — ^. Значит матрица Б с подходящими собственными значениями существует. При этом среди собственных значений матрицы Б обязательно одно будет равным ^, а остальные собственные значения за счет выбора подходящей матрицы Б можно сделать произвольными. Выберем матрицу Б такой, чтобы ее
1
характеристический полином имел вид ф3 (Л) = ^Л — ^ , положив V(Л) = ^Л — ^ I . Легко проверить, что в данном случае = 0, ¿2 = 1, р = 2, й = 0, д =1, Б (Л) = Л — 2, а матрица -Б^Л) имеет вид (3.1). Поэтому Б = Б?. Следуя доказательству теоремы 3.1,
4
4
1
в (3.17) полагаем Б0 = С0 = -, Бу = — / ^ = -,
1
2'
2'
¿>¿2 = $¿2 (СЙ2 0,СЙ2 1) = Ф1(С10,СИ)
0 —1 С10 С11
Числа с10, с11 выбираем так, чтобы IЛ12 — Ф1(с10,с11 )|
1
Л — ^ ) • Для этого полагаем
с10
с11 = 1. Далее формируем матрицу Б = diag[Б0,Бй2, Б^ и, учитывая равенство
4
Б = Б?, окончательно получаем
Б=
2 0 0 0 0 0 —10
0 4 00
10 0 2
§4. Заключение
Изучены структурные свойства линейного автономного дескрипторного уравнения с дискретным временем, возникающего в различных задачах управления системами с многими запаздываниями. На базе канонической формы пучка матриц дано описание матриц исходного уравнения, для которого существует нетривиальное решение. Предложен способ выбора решений этого уравнения, удовлетворяющего заданным асимптотическим свойствам. Представляет интерес дальнейшее развитие методов настоящей работы для исследования аналогичного объекта, но зависящего от случайного параметра [11], а также исследования стабилизации [12] неоднородного уравнения (0.1).
СПИСОК ЛИТЕРАТУРЫ
1. Хартовский В. Е. Об одном линейном автономном дескрипторном уравнении с дискретным временем. I. Приложение к задаче 0-управляемости // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2020. Т. 30. Вып. 2. С. 290-311. https://doi.org/10.35634/vm200211
2. Хартовский В. Е. Обобщение задачи полной управляемости дифференциальных систем с соизмеримыми запаздываниями // Известия РАН. Теория и системы управления. 2009. № 6. С. 3-11. https://www.elibrary.ru/item.asp?id= 12988545
3. Хартовский В.Е. Критерий модальной управляемости вполне регулярных дифференциально-алгебраических систем с последействием // Дифференциальные уравнения. 2018. Т. 54. № 4. С. 514-529. https://doi.org/10.1134/S0374064118040088
4. Хартовский В. Е. Приведение к конечному спектру вполне регулярных дифференциально-алгебраических систем с последействием // Дифференциальные уравнения. 2018. Т. 54. № 6. С. 827-841. https://doi.org/10.1134/S0374064118060110
5. Хартовский В. Е. Синтез наблюдателей для линейных систем нейтрального типа // Дифференциальные уравнения. 2019. Т. 55. № 3. С. 409-422. https://doi.org/10.1134/S0374064119030142
6. Хартовский В. Е. Управление спектром линейных вполне регулярных дифференциально-алгебраических систем с запаздыванием // Известия РАН. Теория и системы управления. 2020. № 1. С. 23-43. https://doi.org/10.31857/S0002338820010084
7. Белов А. А., Курдюков А. П. Дескрипторные системы и задачи задачи управления. М.: Физмат-лит, 2015.
8. Бояринцев Ю. Е. Линейные и нелинейные алгебро-дифференциальные системы. Новосибирск: Наука, 2000.
9. Riaza R. Differential-algebraic systems: Analytical aspects and circuit applications. Hackensack, NY: World Scientific, 2008. https://doi.org/10.1142/6746
10. Гантмахер Р. Ф. Теория матриц. М.: Наука, 1988.
11. Родина Л. И., Тютеев И. И. Об асимптотических свойствах решений разностных уравнений со случайными параметрами // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2016. Т. 26. Вып. 1. С. 79-86. https://doi.org/10.20537/vm160107
12. Зайцев В. А., Ким И. Г., Хартовский В. Е. Задача назначения конечного спектра для билинейных систем с несколькими запаздываниями // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2019. Т. 29. Вып. 3. С. 319-331. https://doi.org/10.20537/vm190303
Поступила в редакцию 10.10.2020
Хартовский Вадим Евгеньевич, к. ф.-м. н., доцент, заведующий кафедрой логистики и методов управления, Гродненский государственный университет им. Я. Купалы, 230023, Беларусь, г. Гродно, ул. Ожешко, 22. E-mail: [email protected]
Цитирование: В. Е. Хартовский. Об одном линейном автономном дескрипторном уравнении с дискретным временем. II. Каноническое представление и структурные свойства // Известия Института математики и информатики Удмуртского государственного университета. 2020. Т. 55. С. 102-121.
V.E. Khartovskii
On a linear autonomous descriptor equation with discrete time. II. Canonical representation and structural properties
Keywords: linear descriptor autonomous equation with discrete time, the subspace of initial conditions, representation of the solution, beam of matrixes.
MSC2010: 93B99, 93C55
DOI: 10.35634/2226-3594-2020-56-08
We consider a linear homogeneous autonomous descriptor equation with discrete time
m
B0g(k + 1) + ^ B^k + 1 - i) = 0, k = m, m + 1,..., i=1
with rectangular (in general case) matrices Bj. Such an equation arises in the study of the most important control problems for systems with many commensurate delays in control: the 0-controllability problem, the synthesis problem of the feedback-type regulator, which provides calming to the solution of the original system, the modal controllability problem (controllability to the coefficients of characteristic quasipolynomial), the spectral reduction problem and the synthesis problem observers for dual surveillance system. The main method of the presented study is based on replacing the original equation with an equivalent equation in the "expanded" state space, which allows one to match the new equation of the beam of matrices. This made it possible to study a number of structural properties of the original equation by using the canonical form of the beam of matrices, and express the results in terms of minimal indices and elementary divisors. In the article, a criterion is obtained for the existence of a nontrivial admissible initial condition for the original equation, the verification of which is based on the calculation of the minimum indices and elementary divisors of the beam of matrices. The following problem was studied: it is required to construct a solution to the original equation in the form g(k + 1) = T0(k + 1), k = 1, 2 ..., where T is some matrix, the sequence of vectors ^(k + 1), k = 1, 2,..., satisfies the equation ^(k + 1) = #0(k), k = 1,2,..., and the square matrix S has a predetermined spectrum (or part of the spectrum). The results obtained make it possible to construct solutions of the initial descriptor equation with predetermined asymptotic properties, for example, uniformly asymptotically stable.
REFERENCES
1. Khartovskii V.E. On a linear autonomous descriptor equation with discrete time. I. Appendix to the 0-controllability problem, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2020, vol. 30, issue 2, pp. 290-311 (in Russian). https://doi.org/10.35634/vm200211
2. Khartovskii V. E. A generalization of the problem of complete controllability for differential systems with commensurable delays, Journal of Computer and Systems Sciences International, 2009, vol. 48, no. 6, pp. 847-855. https://doi.org/10.1134/S106423070906001X
3. Khartovskii V. E. Criteria for modal controllability of completely regular differential-algebraic systems with aftereffect, Differential Equations, 2018, vol. 54, no. 4, pp. 509-524. https://doi.org/10.1134/S0012266118040080
4. Khartovskii V. E. Finite spectrum assignment for completely regular differential-algebraic systems with aftereffect, Differential Equations, 2018, vol. 54, no. 6, pp. 823-838. https://doi.org/10.1134/S0012266118060113
5. Khartovskii V. E. Synthesis of observers for linear systems of neutral type, Differential Equations, 2019, vol. 55, no. 3, pp. 404-417. https://doi.org/10.1134/S0012266119030133
6. Khartovskii V. E. Controlling the spectrum of linear completely regular differential-algebraic systems with delay, Journal of Computer and Systems Sciences International, 2020, vol. 59, no. 1, pp. 19-38. https://doi.org/10.1134/S1064230720010086
7. Belov A. A., Kurdyukov A. P. Deskriptornye sistemy i zadachi upravleniya (Descriptor systems and control problems), Moscow: Fizmatlit, 2015.
8. Boyarintsev Yu. E. Lineinye i nelineinye algebro-differentsial'nye sistemy (Linear and nonlinear al-gebro-differential systems), Novosibirsk: Nauka, 2000.
9. Riaza R. Differential-algebraic systems: Analytical aspects and circuit applications, Hackensack, NY: World Scientific, 2008. https://doi.org/10.1142/6746
10. Gantmakher R. F. Teoriya matrits (Matrix theory), Moscow: Nauka, 1988.
11. Rodina L.I., Tyuteev I.I. About asymptotical properties of solutions of difference equations with random parameters, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2016, vol. 26, issue 1, pp. 79-86. https://doi.org/10.20537/vm160107
12. Zaitsev V. A., Kim I. G., Khartovskii V. E. Finite spectrum assignment problem for bilinear systems with several delays, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2019, vol. 29, issue 3, pp. 319-331. https://doi.org/10.20537/vm190303
Received 10.10.2020
Khartovskii Vadim Evgen'evich, Candidate of Physics and Mathematics, Associate Professor, Head of
Department of Logistics and Methods of Control, Y. Kupala State University of Grodno, ul. Ozheshko, 22,
Grodno, 230023, Belarus.
E-mail: [email protected]
Citation: V. E. Khartovskii. On a linear autonomous descriptor equation with discrete time. II. Canonical representation and structural properties, Izvestiya Instituta Matematiki i Informatiki Udmurtskogo
Gosudarstvennogo Universiteta, 2020, vol. 55, pp. 102-121.