к. т.н. Юничева Н.Р.

 

Институт проблем информатики и управления МОН РК,

 г. Алматы, Республика Казахстан

 

ОЦЕНКА УСТОЙЧИВОСТИ  СИСТЕМ УПРАВЛЕНИЯ

ИНТЕРВАЛЬНЫМИ ОБЪЕКТАМИ

 

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

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

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

,                                                                                              (1)

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

,                         (2)

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

          Определение 1. Интервальная матрица замкнутой системы управления вида (1) асимптотически устойчива, если асимптотически устойчивы все точечные матрицы .

Как известно [1] , характеристическое уравнение матрицы  представляется следующим образом:

                                                                                                (3)

Тогда, для интервальной матрицы  имеет место следующее определение:

          Определение 2. Семейство характеристических уравнений для всех точечных матриц  называется характеристическим уравнением интервальной матрицы   и формально записывается в следующем виде:

                                                                                             (4)

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

          Для получения матричного критерия асимптотической устойчивости воспользуемся известным дробно- линейным преобразованием [2] следующего вида:

                                                                                                  (5)

Данное преобразование позволяет перевести круг заданного радиуса  в левой части плоскости комплексного переменного , в единичный круг с центром в начале координат плоскости комплексного переменного . Значение  из (5) подставим в (3), опуская промежуточные вычисления, получим:

                                                                                                (6)

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

Как известно [6], если являются собственными значениями матрицы , то собственными значениями матрицы будут числа . Следовательно, если матрица   замкнутой системы управления устойчива, то последовательное возведение матрицы  в k–ую степень позволяет уменьшить абсолютную величину собственных значений , так как все     локализованы внутри круга единичного радиуса и по модулю меньше единицы 

                                                                                            (7)

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

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

          Интервальным аналогом точечной матрицы  является матрица следующего вида:

                                                                                                   (8)

где  .    

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

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

          Введем в рассмотрение вещественную матрицу  с элементами :

                                                                                 (9)

Возможность оценки устойчивости по следу матрицы

                                                                                                               (10)

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

.                                                                  (11)

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

,                                                                                                       (12)

Если при исследовании системы (1) выполняется неравенство (12), тогда среди  найдется хотя бы одно, для которого удовлетворяется условие  . Доказательство данного утверждения получено в [3].

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

                                                                            (13)

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

          Процедура вычисления радиуса, охватывающего все собственные значения матрицы замкнутой системы управления приведена в [3].

 

Литература

1.     Чернецкий В.И., Дидук Г.А., Потапенко А.А. Математические методы и алгоритмы исследования автоматических систем. Л.: Энергия, 1970. -375 с.

2.        Карпелевич Ф.И. О характеристических корнях матрицы с неотрицательными элементами. ИАН СССР, серия  15, 1951, - 132 с.

3.        Исследование динамических свойств дискретных интервальных замкнутых систем управления  // Вестник Национальной инженерной академии РК, №2(20), 2008.C.32-37.