Портал | Содержание | О нас | Авторам | Новости | Первая десятка | Дискуссионный клуб | Чат Научный форум --> Первая десятка "Русского переплета" Темы дня:

Президенту Путину о создании Института Истории Русского Народа. |Нас посетило 40 млн. человек | Чем занимались русские 4000 лет назад?

| Кому давать гранты или сколько в России молодых ученых?

Дано понятие линейной динамической системы и ее описание в векторно-матричной форме. Определены ее свойства устойчивости и стабилизируемости. Приведен алгоритм получения стабилизирующей статической линейной обратной связи. В качестве примера рассмотрена система, состоящая из двух упругосвязанных масс.

СТАБИЛИЗАЦИЯ ЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

В. А. БРУСИН

Нижегородский государственный архитектурно-строительный университет

ОСНОВНЫЕ ОПРЕДЕЛЕНИЯ

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

Линейной динамической системой (ЛДС) в теории управления называют математическую модель управляемого объекта, описание которой может быть представлено в виде системы дифференциальных уравнений относительно функций времени t

где ai j(t), bi k(t) - кусочно-непрерывные функции, заданные при всех t. Число n называют размерностью ЛДС.

Набор x(t) = (x1(t), x2(t), _, xn(t))называется вектором состояния ЛДС (точнее, его текущим в момент t значением), набор функций u(t) = (u1(t), u2(t). _. um(t)) - управлением или управляющей функцией [1]. В дальнейшем мы будем предполагать, что коэффициенты ai j(t), bi k(t) не зависят от t. Такие ЛДС называются стационарными (СЛДС). Кроме того, мы ограничим наше рассмотрение случаем, когда m = 1 (скалярное управление).

Для исследования СЛДС и их использования для конкретных задач целесообразно перейти от уравнений (1) к их векторно-матричному аналогу. Для этого мы должны ввести понятие матрицы и дать необходимые свойства. Более подробно о матричной алгебре можно узнать в учебниках по высшей алгебре [2].

ВЕКТОРНО-МАТРИЧНАЯ ФОРМА УРАВНЕНИЙ ЛДС

Матрицей размера n i m, где n и m - натуральные числа, называют массив из n i m чисел, расположенных в виде таблицы, состоящей из n строк и m столбцов:

Числа ai j называются элементами матрицы, при этом первый индекс i указывает на номер строки, а второй индекс j - на номер столбца, в которых находится этот элемент. (Обычно матрицы обозначаются прописными латинскими буквами.) Матрица называется квадратной n-го порядка, если n = m, в противном случае - прямоугольной. Прямоугольная матрица при m = 1, n > 1 называется столбцом, а при n = 1, m > 1 - строкой. (Строки и столбцы обозначаются, как правило, строчными латинскими буквами.) Строки и столбцы называют векторами соответствующей размерности, а их элементы - координатами. Наконец, в случае n = m = 1 мы получаем алгебраическое число - скаляр.

Произвольную матрицу A вида (2) можно рассматривать как совокупность ее столбцов: A = (a1 , a2 , _, aj , _ _, am), где aj - j-й столбец либо - как совокупность ее строк: - i-я строка.

Матрицы можно умножать на число: при умножении на какое-либо число все ее элементы умножаются на это число. Матрицы одних и тех же размерностей можно складывать. При сложении двух матриц A и B одной размерности с элементами ai j и bi j получается матрица C той же размерности с элементами ci j = ai j + bi j .

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

1. Произведением строки a = (a1 , a2 , _, ak) на столбец той же размерности называют число, равное Обозначим это число бa, bс. (Заметим, что изменять порядок сомножителей нельзя, ибо первый множитель - это строка, а второй - столбец.)

2. Пусть теперь A - матрица размера n i m вида (2), а B - матрица размера m i l с элементами bjk , j = 1, 2, _ _, m, k = 1, 2, _, l,

Обозначим через bk = col (b1k , b2k , _, bmk) k-й столбец этой матрицы. Произведением C = AB называют матрицу размера n i l, элементы ci k которой (i = 1, 2, _, n, k = 1, 2, _, l ) равны

Поясним смысл этого действия на примере квадратных матриц 3-го порядка. Предположим, что мы проводим линейную замену старых координат (x, y, z) на "новые" (x', y', z'), а затем эти новые на новейшие (x", y", z") по формулам

Матрицы A с элементами ai j и B с элементами bjk называются матрицами преобразований координат. Тогда, как легко видеть, матрицей преобразования от старых координат к новейшим является матрица AB.

Для произведения матриц справедлив распределительный закон: (AB )C = A(BC ) := ABC.

Матрица называется единичной, если aii = 1, ai j = 0 при i ? j. Единичную матрицу часто обозначают E. Очевидно, EA = AE = A. Для квадратных матриц можно также определить понятие обратной матрицы. Квадратная матрица A, для которой имеется обратная A -1, называется неособой. (В курсе высшей алгебры [2] дается понятие определителя квадратной матрицы det A. Матрица A -1 существует в том и только том случае, если det A ? 0.) Так, в примере с заменой координат обратная матрица B -1 к матрице B будет существовать, если преобразование от старых переменных к новым будет взаимно однозначным, то есть система (3а) дает возможность выразить старые координаты через новые. Матрицей обратного преобразования от новых координат к старым и будет матрица B -1. (Формула прямого вычисления обратных матриц дается в любом курсе высшей алгебры [2].) Очевидно, что произведение прямой и обратной матриц даст единичную матрицу.

Таким образом, если обозначить через B и A матрицы коэффициентов преобразований (3а) и (3б), то эти преобразования с учетом упомянутых выше правил можно записать в виде векторно-матричных равенств:

Введем теперь матрицы A(t) и B(t), составленные из коэффициентов ai j(t) и bik(t) системы (1), а также столбцы x(t) = col (x1(t), x2(t), _, xn(t)), u(t) = col (u1(t), u2(t), _, um(t)). Их элементы будут зависеть от t как от параметра. Тогда, учитывая правила действия с матрицами, систему (1) можно записать в виде одного, но векторно-матричного уравнения

Матрица A(t) называется матрицей коэффициентов ЛДС.

УСТОЙЧИВЫЕ СИСТЕМЫ И МАТРИЦЫ

В дальнейшем мы будем рассматривать стационарные ЛДС, когда элементы матриц A и B не зависят от t. Кроме того, ограничимся случаем скалярного управления m = 1. В этом случае матрица B становится столбцом b = col (b1 , b2 , _, bn). И уравнение (5) принимает вид

Пример 1. Рассмотрим систему из двух движущихся вдоль одной оси масс, связанных упругими элементами, подчиняющимися закону Гука, движение которой управляется перемещением u(t) свободного от масс конца упругого элемента (рис. 1). Пусть положение масс и левого конца упругого элемента определяется координатами s1 и s2 , причем начала отсчета этих величин выбраны так, что при s1 = 0, s2 = 0, u = 0 упругие элементы находятся в нейтральном положении - ни растянуты, ни сжаты. В этом случае уравнения движения масс на основе законов Ньютона и закона Гука будут иметь вид

где mi - величина i-й массы, ki - коэффициент жесткости i-го упругого элемента (i = 1, 2).

Если ввести обозначения:

то эти уравнения можно записать в виде

Теперь если положить

то мы приходим к уравнению вида (6), эквивалентному исходной системе (7). Вектором состояния этой системы служит столбец

Определение. СЛДС вида (6) называется устойчивой, если все решения уравнения

получающегося из (6) при u(t) = 0, стремятся к нулю при t ?: (i = 1, 2, _, n). В этом случае пишут также, что где Q - нулевой вектор-столбец col(0, 0, _, 0). Матрица A устойчивой СЛДС называется устойчивой матрицей.

Из теории дифференциальных уравнений следует, что устойчивость матрицы A связана с корнями ее характеристического многочлена DA(p) = pn + d1 pn - 1 + _ + dn .

Характеристическим многочленом матрицы A является многочлен n-й степени pn + d1 pn - 1 + _ + dn , такой, что все неизвестные функции x1(t), x2(t), _, xn(t), удовлетворяющие системе (1) или (10) при u(t) = 0, удовлетворяют дифференциальному уравнению

Из курса высшей математики [3, 4] известно, что матрица A с постоянными коэффициентами устойчива в том и только том случае, если все корни ее характеристического многочлена имеют отрицательную действительную часть.

Рассмотрим способ получения характеристического многочлена на приведенном выше примере (более распространенный способ связан с вычислением характеристического определителя [3, 4]). Будем исходить из первоначальной системы уравнений (7), положив в ней u = 0. Исключая из этих уравнений переменную s2 (или s1), мы приходим к дифференциальному уравнению 4-го порядка

откуда получаем, что характеристический многочлен матрицы A

К такому же уравнению мы бы пришли, если бы процесс исключения проводился в системе (8)).

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

Система (8) при u = 0 примет следующий вид:

(m1 p2 + k1 - k2)S1 + k2S2 = 0,

(m2 p2 + k2)S2 - k2S1 = 0.

Исключая из этой системы одно из неизвестных, например S2 , мы получаем уравнение

{m1m2 p4 + [(k1 + k2)m2 + k2m1]p2 + k1k2}S1 = 0.

Теперь, если возвратить символу p его первоначальное значение, мы снова приходим к уравнению (12).

Существует критерий Рауса-Гурвица, который дает ответ, устойчив многочлен или нет, не находя его корней [3]. Необходимым условием устойчивости многочлена является условие, чтобы все его коэффициенты имели одинаковые знаки. Таким образом, многочлен DA(p) (12) не является устойчивым. Это отражает физически понятный факт, что при отклонении масс от их положения равновесия они при отсутствии управления будут совершать незатухающие колебания.

СТАБИЛИЗАЦИЯ СЛДС С ПОМОЩЬЮ ЛИНЕЙНОЙ ОБРАТНОЙ СВЯЗИ

ПО СОСТОЯНИЮ

Неустойчивую СЛДС во многих случаях можно стабилизировать, если ею надлежащим образом управлять. Управление может быть программным или может осуществляться с помощью обратной связи [1, 5].

Закон управления вида

u(t) = k1x1(t) + k2x2(t) + _ + knxn(t),

где x = col (x1 , x2 , _, xn) - вектор-столбец состояния СЛДС, а k = (k1 , k2 , _, kn) - вектор-строка коэффициентов, называют линейной (статической и стационарной) обратной связью по состоянию, а регулятор, осуществляющий такую зависимость, - линейным (стационарным) регулятором по состоянию [1]. (Естественно, для того чтобы иметь возможность управлять объектом по закону (13), необходимы приборы, измеряющие текущие значения координат вектора состояния.)

Возникает вопрос: можно ли сделать СЛДС устойчивой с помощью такой обратной связи? Этот вопрос можно облечь в точную математическую формулировку. Закон управления (13) можно записать в векторной форме как

u = kx.

Подставив соотношение (14) в векторно-матричное управление ЛДС (6), приходим к уравнению замкнутой системы (рис. 2)

Матрица n-го порядка

называется матрицей коэффициентов замкнутой системы [1]. Тогда наш вопрос принимает следующую форму.

Задача. Существует ли для заданной пары (A, b) вектор-строка k, такая,что матрица устойчива? Если да, то как его найти. (Если матрица А гурвицева, то ответ, естественно, положителен и ki = 0 (i = 1, 2, _, n).)

Определение (Р. Калман). Пара (A, b) называется стабилизируемой, если существует такой вектор k, при котором матрица (16) устойчива [1].

Ниже мы дадим ответ на поставленную задачу.

Замечание. Представленная здесь формулировка задачи стабилизации по состоянию является частной. Проблема в общей форме формулируется для случая матричных величин B и K.

АЛГОРИТМ НАХОЖДЕНИЯ

ВЕКТОРА КОЭФФИЦИЕНТОВ СТАБИЛИЗИРУЮЩЕЙ ОБРАТНОЙ СВЯЗИ

Зададимся желаемым устойчивым характеристическим многочленом матрицы

Например, мы могли бы взять этот многочлен в виде (p + l1)(p + l2)_(p + ln), li > 0, (- li) - его желаемые корни.

Составим (n i n)-матрицу K (матрицу Р. Калмана)

K = (b, Ab, _, A n - 1b).

Предположим, что эта матрица K имеет обратную. В этом случае пара (A, b) называется полностью управляемой по Р. Калману [1] и по известной теореме [1] стабилизируемой.

Образуем квадратную матрицу S:

где di - коэффициенты многочлена DA(p). Как произведение неособых матриц матрица S также будет неособой.

Сделаем замену вектора-столбца x = col (x1 , x2 , _, xn) в СЛДС (6) на вектор-столбец y = col (y1 , y2 , _, yn) по закону

y = S -1x.

Такая замена переменных будет, конечно, взаимно однозначной. Уравнение СЛДС в новых переменных можно получить умножая левую часть уравнения (6) на S -1. Тогда с учетом (19) приходим к уравнению

Оказывается, что при выборе матрицы S ' в виде (18) матрицы и будут иметь простой вид:

При этом последнее уравнение системы (21) с учетом первых n - 1 уравнений, как легко заметить, принимает вид дифференциального уравнения

где di - коэффициенты характеристического многочлена исходной матрицы A .

Положим

и определим

Тогда из (22)-(24) будет следовать

поэтому закон изменения y1(t) будет определяться желаемым характеристическим многочленом. Так как он выбран устойчивым, то этот процесс и все его производные будут стремиться к нулю при t ? со скоростью не меньшей, чем скорость убывания экспоненциальной функции вида e- lt, .

Производя в (24) обратный переход к старым переменным, получаем

(Очевидно, что при таком переходе асимптотические свойства при t ? не изменятся.)

Обратная связь с вектором-строкой и решает поставленную задачу. Все решения замкнутой системы (6), (25) будут стремиться к нулю при t ? со скоростью не меньшей, чем скорость убывания экспоненты e- lt.

Пример 2. Применим данный алгоритм к системе из примера 1. Положим m1 = m2 = k1 = k2 = 1. Таким образом,

Характеристический многочлен будет равен p4 + 3p2 + 1, так что d1 = 0, d2 = 3, d3 = 0, d4 = 1.

После вычислений получаем

Предположим, что желательным характеристическим многочленом является (p + 1)4 = p4 + 4p3 + 6p2 + + 4p + 1. В этом случае, согласно (23), и искомая обратная связь, согласно (25), будет иметь вид

ЗАКЛЮЧЕНИЕ

Итак, мы привели решение задачи стабилизации линейных динамических систем с помощью линейной обратной связи по состоянию. Это решение позволяет помимо главной цели получать заданное распределение корней характеристического многочлена замкнутой системы. Положительной чертой рассмотренной обратной связи является ее простота. Однако она имеет и существенные недостатки. Далеко не во всех задачах управления реальными объектами текущие значения всех координат вектора состояний могут быть измерены. Гораздо чаще возникает ситуация, когда только часть этих координат может быть измерена по техническим причинам. Поэтому возникла задача о стабилизации с помощью линейной обратной связи (14), которая зависела бы только от некоторых координат вектора состояния. (Так, в рассмотренном примере может оказаться, что доступны измерению положение и скорость только правой массы или что можно измерять только положения масс, но не их скорости.) При этом корни характеристического многочлена замкнутой системы должны попасть в заданную область в левой части комплексной плоскости. Ясно, что решение такой задачи не всегда возможно. Полного решения этой задачи, несмотря на ее долгую историю, в настоящее время нет.

ЛИТЕРАТУРА

1. Андреев Ю.Н. Управление конечномерными линейными объектами. М.: Наука, 1976. 423 с.

2. Смирнов В.И. Курс высшей математики. М.: ГИТТЛ, 1953. Т. 3, ч. 1. 339 с.

3. Понтрягин Л.С. Обыкновенные дифференциальные уравнения. М.: Наука, 1974. 331 с.

4. Мышкис А.Д. Лекции по высшей математике. М.: Наука, 1969. 640 с.

5. Брусин В.А. Об управлении динамическими системами в условиях неопределенности // Соросовский Образовательный Журнал. 1996. ╧ 6. С. 115-121.

Рецензент статьи В.Б. Колмановский

* * *

Владимир Александрович Брусин, доктор физико-математических наук, профессор, зав. кафедрой высшей математики Нижегородского государственного архитектурно-строительного университета, член-корреспондент РАЕН. Область научных интересов - математические проблемы теории устойчивости и теории управления. Автор более 180 научных статей и двух учебных пособий.