ON PRESERVATION OF THE QUADRATIC LYAPUNOV FUNCTION OF A LINEAR DISCRETE SYSTEM UNDER STATIONARY PERTURBATIONS OF THE SYSTEM COEFFICIENTS

Research article
DOI:
https://doi.org/10.23670/IRJ.2024.140.47
Issue: № 2 (140), 2024
Suggested:
10.01.2024
Accepted:
08.02.2024
Published:
16.02.2024
410
7
XML
PDF

Abstract

For a discrete system described by a linear pointwise mapping, sufficient conditions on the smallness of perturbations of the system coefficients are obtained, under which the quadratic Lyapunov function constructed for the original system will be a Lyapunov function for the perturbed system as well. The fact of negativity of the first difference in the power of the system is determined by the fact of negativity of the roots of the auxiliary polynomial, whose coefficients depend on the coefficients of the system. Sufficient conditions for preserving the fact of negativity of the roots of the constructed polynomial at finite perturbations of the system coefficients are obtained. One of the methods for determining the coefficients of the quadratic Lyapunov function possessing the given properties is examined: a Lyapunov function satisfying the restriction on the value of its first difference and convenient for use in evaluating the quantitative characteristics of the system.

1. Введение

Развитие теории дискретных динамических систем в направлении исследования макроструктуры пространства состояний является перспективным как с теоретической, так и с прикладной точек зрения, т.к. такая постановка вопроса оказывается естественной при решении задач различного технического содержания

. Для решения задачи исследования макроструктуры пространства состояний дискретных динамических систем имеет смысл использовать прямой метод Ляпунова, распространенный на дискретные системы
,
и позволяющий оценивать макроструктуру пространства состояний на основе принципа достаточности
.

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

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

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

).

2. Квадратичная функция Ляпунова, удовлетворяющая заданному ограничению для линейного точечного отображения

Для фиксированного img рассмотрим дискретную динамическую систему, описываемую линейным точечным отображением с постоянными коэффициентами

img
(1)

и квадратичную форму с постоянными коэффициентами

img
(2)

Элементы постоянных img-матриц img и img вещественные. Первая разность квадратичной формы img (2) в силу дискретной системы (1) также является квадратичной формой:

img
(3)

Далее считаем, что собственные значения img матрицы коэффициентов img дискретной системы (1), т.е. корни характеристического уравнения img, лежат внутри единичного круга и что квадратичная форма (2) положительно определена (здесь и ниже img – единичная img-матрица). Это предположение означает, в частности, что у дискретной системы (1) имеется

квадратичная функция Ляпунова. В работе
доказана.

Теорема 1. Пусть квадратичная форма (2) является функцией Ляпунова дискретной системы (1) и на заданной поверхности уровня img функции img максимальное значение её первой разности (3) в силу (1) равно img Тогда для собственных значений img матрицы коэффициентов img справедливо неравенство img, а коэффициенты квадратичной формы (2) удовлетворяют равенству

img
(4)

Замечание 1. В дальнейшем матрицу img квадратичной формы (3) будем обозначатьimg Очевидно, что матрица img симметрична: img. Пусть img. Если img и img, то для элементов матрицы img очевидны равенства

img
(5)

в которых учтена симметричность матрицы img. В этих обозначениях (5) запишется в виде img

Доказательство теоремы, приведённое в работе

, основывается на решении соответствующей экстремальной задачи. Дадим другое её доказательство.

Доказательство. Так как (2) – положительно определённая квадратичнаяформа, то корректно задан регулярный пучок img квадратичных форм

, где img – числовой параметр. Пусть img характеристические числа этого пучка
, т.е. корни уравнения

img
(6)

(все корни этого уравнения вещественные

). Тогда, согласно
, имеет место равенство 

img
(7)

Так как отношение указанных в (7) квадратичных форм на любой прямой, проходящей через начало координат, принимает, за исключением точки img одно и то же значение (вообще говоря, своё для каждой прямой), то

img
(8)

Из равенств (7) и (8) заключаем, что img, а значит, вследствие того, что img – корень уравнения (6), справедливо равенство (4). По условию теоремы квадратичная форма (2) является функцией Ляпунова дискретной системы (1), поэтому квадратичная форма (3) отрицательно определена, а значит, imgСледовательно, для завершения доказательства теоремы остаётся установить справедливость неравенства img. Пусть, без нарушения общности, модуль собственного значения img – матрицы img не меньше, чем модули остальных её собственных значений. Могут представиться только два случая: либо 1) img, либо img. Рассмотрим каждый из них отдельно.

В случае 1) обозначим через img собственный вектор матрицы img отвечающий собственному значениюimg. Тогда img и img. Поэтому верно равенство img в силу которого получаем оценку img

В случае 2) собственному значению img матрицы img отвечает в пространстве img собственный вектор img, где img Раскрывая в равенстве img скобки и отделяя вещественную и мнимую части, получаем img и, следовательно, img В силу двух последних соотношений несложно убедиться в том, что имеет место равенство

img
(9)

Выберем тот вектор img, для которого второе слагаемое в правой части равенства (9) неотрицательно. Тогда при таком img вследствие (9) справедливо неравенство img, в силу которого и положительной определённости квадратичной формы (2) получаем оценку img.

Теорема доказана.

Заметим, что уравнение (6) может быть записано в виде img, где 

img
(10)
img
(11)

(причем img), а матрица img получается из матрицы img заменой столбцов с номерами img на одноименные столбцы матрицы img. Многочлен (10)-(11) получается из (6), если раскрыть определитель и разделить полученный многочлен на коэффициент при старшей степени img. Как отмечено в доказательстве теоремы 1, все корни многочлена img (см. (10)) вещественные

.

Предположим, что построена квадратичная форма (2), которая является функцией Ляпунова для отображения (1). Пусть коэффициенты (1) изменились, т.е. функции последования отображения приняли вид 

img
(12)

где img – постоянная вещественная матрица. Тогда для матрицы img первой разности img квадратичной формы (2) в силу (12) очевидно равенство

 img 

Поэтому величины img (см. (5)) преобразуются к виду img, где 

img
(13)

уравнение (6) примет вид

img
(14)

а вместо многочлена (10) получим многочлен 

img
(15)

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

img
(16)

здесь

img
(17)
img
(18)

где img – алгебраическое дополнение элемента img в матрице img, а img – алгебраическое дополнение минора матрицы img, построенного на строках с номерами img  и столбцах с номерами img. Многочлен (15)-(16) получается (14) при учете (13), если раскрыть определитель и разделить полученный многочлен на коэффициент при старшей степени img.

Задача, рассматриваемая в работе, состоит в нахождении условий, которым должны удовлетворять возмущения, чтобы квадратичная функция Ляпунова (2), построенная для дискретной (1), оставалась бы функцией Ляпунова и для возмущённой дискретной системы (12). Уточним постановку задачи.

2.1. Об условиях, дающих решение задачи

Функция Ляпунова (2) построенная для системы (1), будет оставаться функцией Ляпунова и для системы (12), если и только если все собственные значения матрицы img будут отрицательными. Так как матрица imgсимметричная, то все её собственные значения вещественны

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

Теорема 2. Квадратичная форма (2), являющаяся функцией Ляпунова для дискретной системы (1), будет оставаться функцией Ляпунова и для возмущенной системы (12), если и только если для матрицы img выполняются неравенства img

Таким образом, чтобы проверить, будет ли квадратичная форма (2), являющаяся функцией Ляпунова дискретной системы (1), также и функцией Ляпунова возмущенной системы (12), нужно в дополнение к найденным положительным коэффициентам img (см. формулы (11)) найти величины img (см. формулы (16)) и проверить выполнимость неравенств img. Для проверки этого необходимого и достаточного условия необходимо вычислить все величины img, что связано с большим количеством довольно объёмных вычислений. Можно предложить достаточное условие, для проверки которого требуется вычислить фактически только элемент img. Чтобы сформулировать это достаточное условие, рассмотрим наряду с точечным отображением (12) также параметрическое семейство 

img
(19)

параметр img в котором изменяется на отрезке [0,1] и в которое (1) и (12) входят при img и imgсоответственно. Обозначим через img величину img, вычисленную при фиксированном img для (19) по формуле (16), т.е. для (19) нужно в формулах (17) и (18) заменить img на img и img на img. Тогда очевидно, что величина img представляет собой многочлен степени img переменной img с нулевым свободным членом, а именно, 

img
(20)

Обозначим через img многочлен img (см. (15)), построенный для (19). Свободный член многочлена img равенimg и является многочленом степени img с положительным свободным членом. Тогда, согласно [8], будут иметь место следующие утверждения.

Лемма 1. Если многочлен img при всех img положителен, то квадратичная форма (2), являющаяся функцией Ляпунова для системы (1), будет функцией Ляпунова и для системы (12).

Заметим, что матрица img  будет невырожденной, если норма матрицы-возмущения img  не превосходит радиуса невырожденности [12] матрицы img. Другими словами, имеет место

Лемма 2. Если матрица img отрицательно определена, то матрица img будет отрицательно определённой при любой матрице img, удовлетворяющей неравенству 

img,

где img – собственные значения матрицы img. (Через img обозначена спектральная норма матрицы.)

Из леммы 2 очевидно вытекает

Лемма 3. Если матрица img отрицательно определена и выполняется неравенство img то отрицательно определенной будет и матрица img.

2.2. Примеры построения в двух частных случаях квадратичной функции Ляпунова с заданным максимумом на её поверхности уровня первой разности

Пусть квадратичная функция Ляпунова для системы (1) строится в соответствии с методикой работы

, основанной на переходе к каноническим координатам.

1. Рассмотрим случай, когда все корни характеристического уравнения дискретной системы (1) вещественны и различны; без нарушения общности считаем, что img Тогда существует линейное невырожденное преобразование координат img приводящее точечное отображение к виду

 

img
(21)

(У матрицы img её img-й столбец является собственным вектором матрицы img, соответствующим собственному значению img). При этом квадратичные формы (2) и (3) перейдут соответственно в квадратичные формы img и img, причём

img.

Таким образом, не ограничивая общности, можно предполагать, что система (1) имеет канонический вид (21). Тогда равенства (5) примут вид

img,

а уравнение (6) запишется следующим образом:

img
(22)

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

Квадратичную функцию Ляпунова с img зададим следующим образом:

img
(23)

где img,

а img.

То есть матрица img построенной квадратичной формы имеет вид 

img.

Несложно убедиться, что при условии img величина img принадлежит полуинтервалу img, а квадратичные формы с матрицами img и img (матрицы img и img  определены в (21) и (24)) и соответственно положительно и отрицательно определены. Корнями уравнения (22) img являются: (то, что первые img из этих значений корни уравнения (22), очевидно, в том же, что два последних из них корни, проще всего убедиться непосредственной подстановкой в уравнение (22)). При этом в рассматриваемом случае при выборе (24) матрицы квадратичной формы наименьшим и наибольшим корнями уравнения (22) будут img и img.

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

Если img, и только img, то для канонической системы дифференциальных уравнений

квадратичную функцию Ляпунова, удовлетворяющую условию img можно искать в виде (23), где 

img
(25)

а imgсчитая равными между собой коэффициенты img, если переменные img соответствуют паре комплексно-сопряжённых корней характеристического уравнения.

В самом деле, пусть корни характеристического уравнения img действительны и различны, и только img. Тогда существует линейное невырожденное преобразование координат, приводящее точечное отображение к каноническому виду

img
(26)

При этом квадратичные формы (2) и (3) перейдут соответственно в квадратичные формы img и img, причём

 

img.

Таким образом, не ограничивая общности, можно предполагать, что система (1) имеет канонический вид (26).

Нетрудно убедиться, что корнями уравнения (6) являются: img Таким образом, в рассматриваемом случае при выборе (25) матрицы квадратичной формы наибольшим корнем уравнения (6) будет img

Рассмотрим теперь случай, когда среди корней img есть пара комплексно-сопряженных корней, скажем img. Тогда линейным невырожденным преобразованием координат система может быть приведена к виду, подобному (26), но с матрицей

img
(27)

Если искать квадратичную функцию Ляпунова для канонической системы в виде можно искать в виде (23), где img связаны в силу (25), считая равными между собой коэффициенты img, то корнями уравнения (6) будут: imgimg. Таким образом, в рассматриваемом случае при выборе (25) матрицы квадратичной формы наибольшим корнем уравнения (6) будет img.

Случай, когда среди корней характеристического уравнения имеется несколько пар комплексно-сопряженных корней, рассматривается аналогично.

Таким образом, если наибольшему значению действительных частей корней характеристического уравнения соответствует пара комплексно-сопряженных корней, корнями уравнения (6) являются: img И в рассматриваемом случае при выборе в силу (24), (25) матрицы квадратичной формы наибольшим корнем уравнения (6) будет img.

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

3. Заключение

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

.

Article metrics

Views:410
Downloads:7
Views
Total:
Views:410