Например, Бобцов

ЭКСПОНЕНЦИАЛЬНАЯ УСТОЙЧИВОСТЬ НЕПРЕРЫВНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

ЭКСПОНЕНЦИАЛЬНАЯ УСТОЙЧИВОСТЬ НЕПРЕРЫВНЫХ ДИНАМИЧЕСКИХ СИСТЕМ
3 АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ

УДК 681.5
ЭКСПОНЕНЦИАЛЬНАЯ УСТОЙЧИВОСТЬ НЕПРЕРЫВНЫХ ДИНАМИЧЕСКИХ СИСТЕМ
С.В. Быстров, В.В. Григорьев, Е.Ю. Рабыш, Н.А. Черевко

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

Введение

Для оценки качества переходных процессов в системах с одним входом и одним выходом широкое распространение в инженерной практике получили показатели качества, вводимые по переходной функции, т.е. по реакции системы на единичное ступенчатое воздействие при нулевых начальных условиях. Такие показатели качества процессов, как время переходного процесса и перерегулирование, достаточно полно отражают физические процессы, происходящие в системе, и позволяют сравнивать динамические свойства различных систем. Однако вводимые по переходной функции показатели качества не применимы непосредственно к многомерным системам. Множество траекторий в пространстве выходов многомерной системы, вызванных воздействиями, приложенными к различным входам, может состоять из существенно разнящихся визуально кривых. На множестве таких кривых утрачивается очевидность их сравнения, а следовательно, и сопоставимость динамических свойств различных систем. Аналог переходной функции для линейных непрерывных многомерных систем введен в [1]. Дальнейшее развитие подобного подхода нашло отражение во введении понятия совокупной переходной функции. Для этого вводится понятие нормированной совокупной переходной функции многомерной системы [2, 3], по которой, как и для систем с одним входом и одним выходом, определяются типовые показатели качества в виде времени переходного процесса и перерегулирования. Причем эти показатели качества определяются по множеству траекторий, исходящих из области допустимых начальных значений. Чтобы избежать при нахождении значений показателей качества большого объема вычислений, вызванного необходимостью выявления наихудших случаев по нормированию совокупных переходных функций, соответствующих исследуемому множеству траекторий, проще производить оценку показателей качества. В связи с этим в дальнейшем на основе условий экспоненциальной устойчивости были получены показатели качества в виде времени переходного процесса для непрерывных и дискретных динамических систем. Но следует отметить, что лишь для дискретных динамических систем по полученным показателям качества были разработаны алгоритмы синтеза оптимальных регуляторов, обеспечивающие экспоненциальную устойчивость [2].

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

Современная теория линейных систем автоматического управления основана на использовании

метода пространства состояний. Среди различных направлений основанный на нем метод аналитическо-

го конструирования оптимальных регуляторов является одним из самых распространенных [4–7]. Белл-

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

«менее важной задачей» часто заменяют исходную, «более реалистичную задачу» оптимизации, и под-

черкивал: «Это напоминает историю об одном человеке, который, потеряв кольцо посреди улицы, искал

его под фонарем, потому что там светлее», хотя оно «оставалось в той непроглядной темноте, которая

показалась слишком затруднительной для поисков».

Предположим, что поведение динамической системы описывается дифференциальным уравнением

xt  Atxt Btut,

(1)

где x  Rn – вектор состояния динамической системы; x0  x0  Rn – вектор начальных состояний;
u  Rm – управление; A – n n матрица описания; B – n m матрица входа; t  0 – время. Полагая, что объект управления (1) является полностью управляемым, рассмотрим задачу нахождения такого закона управления, который обеспечивает для замкнутой системы экспоненциальную устойчивость с заданными значениями параметров, определяемыми по требуемым значениям времени переходных процессов.

Стабилизация нестационарных непрерывных динамических систем

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

44 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 3 (73)

С.В. Быстров, В.В. Григорьев, Е.Ю. Рабыш, Н.А. Черевко

При изучении условия экспоненциальной устойчивости будем полагать, что при любых x0  Rn

решение x  Rn уравнения (1) существует и единственно, а система является полной в Rn , т.е. под этим

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

определение экспоненциальной устойчивости [3].

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

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

x0  Rn , существуют такие параметры  :  1 и  :   0 , при которых в любой момент времени t  0

выполняется условие

xt    et x0 ,

(1)

где норма вектора задается соотношением

 n
x 

xi

2

1/ 

2

,

i1 

где xi – i-ая компонента вектора состояния x .

Неравенство (2) устанавливает экспоненциальную мажоранту процессов, которая может служить

для оценки качества непрерывной системы в некоторой окрестности ее положения равновесия.

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

лучить оценку времени переходного процесса [3]:

tп



1 

lnп



,

(3)

где п  0,05 .

Для линейных систем (1) подлежит минимизации квадратичный критерий качества:

 
J x,u   xT tQtxt uT tRtut dt ,

(4)

0

где Q  QT – положительно полуопределенная n n матрица штрафов по отклонениям траекторий дви-

жения системы от положения равновесия; R  RT – положительно определенная m m матрица штра-

фов по управлению. Управление u в задаче (1), (4) ищется в виде

ut  Ktxt ,

(5)

удовлетворяющее условию (2). Здесь K t – m  n матрица линейных обратных связей по состояниям

объекта управления. Для решения задач (1), (4) необходимо построить функцию V xt  , которую в рас-

сматриваемом случае естественно искать в виде квадратичной формы:
Vxt  xT tPtxt ,

(6)

где P  PT – положительно определенная n n матрица.
Теперь дадим локальные достаточные условия экспоненциальной устойчивости [3]. Теорема 1. Непрерывная система (1) экспоненциально устойчива в точке x  0 , если существуют
такая квадратичная функция Ляпунова V xt  и такой параметр  :   0 , при которых на всех траекто-

риях движения системы в любой момент времени t  0 выполняется условие
Vxt  2V xt.

(7)

Условие (7) имеет место, если справедливо уравнение
Vxt 2V xt  xT tQtxt uT tRtut,

(8)

подставив квадратичную функцию Ляпунова (6) и ее производную в уравнение (8), получим:
Atxt Btut xtT Ptxt xT tPtAtxt Btut xt xT tPtxt

 xT tQtxt uT tRtut  0 ,

(9)

на основе которого, воспользовавшись методом локальной оптимизации [3], находим оптимальное

управление:

ut  R1tBT tPtxt .

(10)

Таким образом, используя (4), (8), (9), получим систему матричных уравнений типа Риккати:

Pt At BtKt IT Pt PtAt BtKt I KT tRtKt Qt  0 ,

(11)

Kt  R1tBT tPt ,

(12)

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 3 (73)

45

ЭКСПОНЕНЦИАЛЬНАЯ УСТОЙЧИВОСТЬ НЕПРЕРЫВНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

где I – единичная матрица. Если все матрицы A, B, Q, R постоянны, то процесс может достичь устано-

вившегося состояния в том смысле, что становится постоянной матрица P ( P  0 ):
xt  Axt But ,

A  BK  I T P  PA  BK  I  KT RK  Q  0 ,

(13)

K  R1BT P .

(14)

На основании теоремы 1 существование положительно определенного решения системы уравне-

ний (13), (14) является достаточным условием устойчивости синтезируемой матрицы замкнутой систе-

мы:

Re iF  I  0, i  1, 2,..., n ,

и, следовательно, получения степени устойчивости системы:

Re iF  , i  1, 2,...,n .

В заключение рассмотрим последовательность синтеза регулятора для обеспечения экспоненци-

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

данными – в общем случае управляемым нестационарным объектом управления с матрицами описания

At и входа Bt :

1. по заданным оценкам времени переходного процесса определение значения параметра  по (3);
2. решение системы матричных уравнений типа Риккати (11), (12) относительно матрицы Pt с по-

следующим вычислением матрицы линейных обратных связей K t по (12).

В качестве примера рассмотрим непрерывный стационарный объект управления
xt  Axt But, yt  Cxt

A



0 1

15; B  04;C  1

1

со следующими заданными показателями качества:

tп  3,0 с; п  0,01 .

В общем случае, предполагая, что все переменные вектора состояния объекта управления доступ-

ны для измерения, управление ищется как функция состояний объекта управления в виде

u  Kx .

Взяв параметры

Q



0 0

0 0

;

R



1

и решив матричные алгебраические уравнения (13), (14), получим матрицы обратных связей K регуля-

торов при   1,535 (вычисляемая по (2) и обеспечивающая экспоненциальную устойчивость) и при

  0 (обеспечивается асимптотическая устойчивость) соответственно:

K1  6,194 4,035 .

K2  2,447 1,856 .

y

0,5
0 1
–0,5

2

–1
–1,5 0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 t, с
 Рисунок. Переходные процессы при x0T  1 0 : 1 – в экспоненциально устойчивой системе;
2 – в асимптотически устойчивой системе
46 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 3 (73)

А.Д. Громышева, И.Е. Овчинников, А.В. Егоров

На рисунке показаны графики переходных процессов при x0T  1 0 , из которых явно видно, что
на переходном процессе 1 (в отличие от переходного процесса 2), который получается при использовании условий экспоненциальной устойчивости, удовлетворяется желаемое время переходного процесса, что и подтверждает полученные результаты.

Заключение

На основе достаточных условий экспоненциальной устойчивости, метода локальной оптимизации и оценки времени переходного процесса получен алгоритм синтеза оптимальных регуляторов, обеспечивающий синтезируемой непрерывной нестационарной динамической системе экспоненциальную устойчивость.
Работа выполнена при финансовой поддержке РФФИ (грант 09-08-00857-а «Методология применения теории качественной устойчивости при проектировании систем управления адаптивной оптикой»).

Литература

1. Соболев О.С. Однотипные связанные системы регулирования. – М.: Энергия, 1973.– 136 с. 2. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А.В. Синтез дискретных регуляторов при
помощи ЭВМ. – Л.: Машиностроение, Ленингр. отделение, 1983. – 245 с. 3. Григорьев В.В. Аналитический синтез регуляторов на основе качественной устойчивости [Текст]:
дис. ... д-ра техн. наук: 05.13.01/ В.В. Григорьев. ЛИТМО. – Л.: [б.н.], 1988.– 483 с. 4. Афанасьев В.Н., Колмановский В.Б., Носов В.Р. Математическая теория конструирования систем
управления. – М.: Высшая школа, 2003.– 447 с. 5. Мирошник И.В. Теория автоматического управления. Нелинейные и оптимальные системы. – СПб:
Питер, 2006.– 272 с. 6. Фурасов В.Д. Устойчивость движения, оценки и стабилизация. – М.: Наука, 1978.– 248 с. 7. Красовский А.А., Буков В.Н. и Шендрик В.С. Универсальные алгоритмы оптимального управления
непрерывными системами. – М.: Наука, 1977.– 272 с.

Быстров Сергей Владимирович Григорьев Валерий Владимирович Рабыш Евгений Юрьевич Черевко Николай Александрович

– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат технических наук, доцент, sbystrov@mail.ru
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, доктор технических наук, профессор, grigvv@yandex.ru Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант, Rabysh@yandex.ru Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант, epostbox1@mail.ru

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 3 (73)

47