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

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

ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ И СИСТЕМЫ
УДК 62-50
Д. С. БИРЮКОВ, А. В. УШАКОВ
КОНТРОЛЬ ЭНЕРГЕТИЧЕСКИХ ЗАТРАТ НА УПРАВЛЕНИЕ ПРИ ВОСПРОИЗВЕДЕНИИ
ПОЛИНОМИАЛЬНЫХ ЭКЗОГЕННЫХ ВОЗДЕЙСТВИЙ: ГРАМИАННЫЙ ПОДХОД
Рассматривается задача контроля энергетических затрат на управление при воспроизведении синтезируемой системой полиномиальных экзогенных воздействий с использованием грамианного подхода. Задача имеет прозрачное решение для начальных состояний источника таких воздействий в виде мажоранты и миноранты затрат на управление как функций распределения мод.
Ключевые слова: источник полиномиального экзогенного воздействия, вектор начального состояния, объект управления, установившаяся составляющая, грамиан затрат на управление, минорантная и мажорантная оценки.
В работах [1, 2] рассматривалась задача выбора желаемой структуры мод, осуществляемого при контроле энергетических затрат на управление начальным состоянием технического объекта (ТО), из множества возможных желаемых структур мод при условии выполнения отношения порядка полиномов [3] применительно к длительности переходного процесса и возможному перерегулированию. Также в работах [1, 2] были выделены некоторые типовые источники экзогенных воздействий, для которых проблему оценивания затрат на управление следовало рассмотреть отдельно. В настоящей статье эта проблема формулируется как задача контроля затрат на управление объектом при воспроизведении синтезируемой системой полиномиальных экзогенных воздействий с использованием грамианного подхода.
Грамианный подход сформировался в рамках современной теории управления, опирающейся на векторно-матричный формализм метода пространства состояний. Подход устанавливает аналитическую связь компонентов движения ТО, входящего в состав синтезируемой системы, по выходу и вектору начального состояния источника полиномиальных экзогенных воздействий. При этом наибольший практический интерес в рассматриваемой задаче представляет установившаяся составляющая движения объекта. Алгоритмически данная связь опирается на концепцию векторно-матричного подобия, матрица которого является решением уравнения Сильвестра. Оценка затрат на управление при воспроизведении системой установившегося движения по выходу формируется на основе закона управления, записанного в аддитивной форме, одно слагаемое которого должно обеспечивать требуемое качество воспроизведения полиномиальных экзогенных воздействий, а второе — требуемое качество выходного сигнала системы при установившейся составляющей движения. Очевидно, в установившемся режиме сохраняется только первая составляющая закона управления, что позволяет связать его с начальным состоянием источника полиномиальных экзогенных воздействий.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9

6 Д. С. Бирюков, А. В. Ушаков

Для получения аналитического соотношения между компонентами движения динамиче-

ской системы и вектором начального состояния источника полиномиального экзогенного воздействия (ПЭВ) рассмотрим технический объект, математическая модель которого описы-

вается выражением

x (t ) = Ax (t ) + Bu (t ), y (t ) = Cx (t ), x (t ) t−0 = x (0), ε (t ) = g (t ) − y (t ) ,

(1)

где x(t) — вектор состояния объекта, и(t) — вектор входного воздействия (управление), у(t) —

вектор выходного воздействия, g(t) — вектор экзогенного воздействия, ε(t) — вектор ошибки:

x ∈ Rn; u ∈ Rr ; y, g, ε ∈ Rm; A, B, C — соответственно матрицы состояния, управления и

выхода, размерности которых согласованы с перемножаемыми переменными.

Математическая модель источника полиномиального экзогенного воздействия g (t )

задается в форме

z (t ) = Ez (t ), z (t ) t=0 = z (0), g (t ) = Pz (t ) ,

(2)

где z(t) — вектор состояния источника, z ∈ Rl ; E, P — матрицы состояния и выхода,

согласованные по размерности с размерностями векторов z и g .

Для источника ПЭВ, в отличие от источника гармонического экзогенного воздействия [2], невозможно привести одну универсальную реализацию матриц E и P , поскольку

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

⎡0 1 0

0⎤

⎢⎢0 0 1

0⎥⎥

E = ⎢0 0

0 0⎥ , P = [1 0 … 0] .

(3)

⎢⎢0 0 0

1⎥⎥

⎢⎣0 0 0

0⎥⎦

Например, для источника ПЭВ третьего порядка математическая модель будет

выглядеть следующим образом:

⎡0 1 0⎤
E = ⎢⎢0 0 1⎥⎥ , P = [1 0 0] .
⎣⎢0 0 0⎥⎦

(4)

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

u (t ) = uλ (t ) + ug (t ) = Kg g (t ) − Kx (t ) ,

(5)

где uλ (t ) — составляющая управления, обеспечивающая желаемое распределение мод;

ug (t ) — составляющая управления, обеспечивающая соответствие выходного сигнала за-

дающему воздействию; K , Kg — матрица обратных связей.

Дальнейший ход вычислений приведен в работе [2].

В результате объединения выражений (1) и (5) с получением модели замкнутой системы

и конструирования агрегированной системы с вектором состояния x (t ) = ⎣⎡xT (t ) zT (t )⎤⎦T и

матрицами состояния F и G , которые заданы соотношениями F = A − BK, G = BKg , полу-

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

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9

Контроль энергетических затрат на управление

7

x (t ) = Ix (t ) + 0z (t ) = [I 0] x (t ), x (t ) = Cxx (t ) = CxeFt x (0) ;

(6)

y (t ) = Cx (t ) + 0z (t ) = [C 0] x (t ), x (t ) = Cy x (t ) = CyeFt x (0) ;

(7)

ε(t ) = −Cx (t ) + Pz (t ) = [−C P] x (t ), x (t ) = Cεx (t ) = CεeFt x (0) ;

(8)

z (t ) = 0x (t ) + Iz (t ) = [0 I ] x (t ), x (t ) = Cz x (t ) = CzeFt x (0) .

(9)

В целях дальнейших исследований сформулируем следующее утверждение.

Утверждение 1. Если матрицы F , E, G, P связаны матричным уравнением Сильвестра

TE − FT = GP ,

(10)

то матричная экспонента eFt может быть представлена в виде

eFt

=

⎡eFt ⎢

⎢⎣ 0

TeEt − eFtT eEt

⎤ ⎥ ⎦⎥

.

□ (11)

Доказательство утверждения приведено в работе [4].



Используя выражение (11), уравнения движения агрегированной системы можно

записать следующим образом:

x

(t

)

=

⎡ ⎢ ⎣

x z

(t (t

)⎤ )⎦⎥

=

⎡eFt ⎢ ⎢⎣ 0

TeEt − eFtT eEt

⎤ ⎥ ⎥⎦



⎡ ⎢ ⎣

x z

(0)⎤ (0)⎦⎥

.

Тогда уравнения движения замкнутой системы приобретают следующий вид:

( )x (t ) = eFt x (0) + TeEt − eFtT z (0), y (t ) = Cx (t ) ,

(12) (13)

что позволяет выделить в движении системы по вектору состояния x (t ) установившуюся со-

ставляющую:

x

y

(

t

)

=

lim
t→∞



(t

)

=

Te

Et

z

(

0

)

,

(14)

где xв (t ) — вынужденная составляющая движения системы.

Утверждение 2. Матрица T — решение уравнения Сильвестра (10) — представляет

собой матрицу подобия, связывающую установившуюся составляющую xy (t ) движения

системы и вектор состояния z (t ) источника полиномиального экзогенного воздействия (2). □

Доказательство утверждения приведено в работе [2].



Утверждение 3. Решение уравнения Сильвестра (10) для случая полиномиального экзо-

генного воздействия порядка п представимо в форме

T = ⎡⎣−F −1G −F −2G … −F −nG⎤⎦ .

□ (15)

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

статическую, кинетическую и динамическую составляющие, матрицы E и P модели (2) ис-

точника ПЭВ имеют вид (3). Если матрицу T представить в столбцовой форме, записав ее

как

T = [T1 T2 … Tn ] ,
то уравнение Сильвестра (10) преобразуется к виду

[T1 T2 … Tn ] E − F [T1 T2 … Tn ] = GP .

(16)

Подстановка в уравнение (16) матриц E и P в форме (3) позволяет для столбцов матрицы T

записать T1 = −F −1G, T2 = −F −2G, … Tn = −F −nG , что приводит к выражению (15).



ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9

8 Д. С. Бирюков, А. В. Ушаков

Отметим, что наличие известной из метода модального управления связи M Γ = FM

матрицы Γ состояния модальной модели и матрицы F состояния системы, а также

полученное соотношение (15) для матрицы T в уравнении Сильвестра (10) позволяют

существенно модифицировать сформированный в работе [2] алгоритм оценивания затрат на

управление установившимся движением проектируемой системы в функции от заданного

распределения мод матрицы состояния F . Матрица подобия M определяется из уравнения

Сильвестра: M Γ − AM = −BH , где Н — матрица выхода модальной модели.

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

управление установившимся движением системы, определенный в работе [2]:

( ) ( )∫WU (t ) = t eET τ Kg P − KT T Kg P − KT eEτdτ .

(17)

0

При t → ∞ справедливо соотношение tl→im∞WU (t ) = WU . Грамиан WU затрат на

управление удовлетворяет матричному уравнению Ляпунова:

( ) ( )ETWU +WU E = − Kg P − KT T Kg P − KT .

(18)

В соответствии с выражением (17) можно вычислить оценку затрат U ∞ на управление

( )для начальных состояний источника ПЭВ

z (0) = const

в форме мажоранты

U[0,∞)

и
max

миноранты U[0,∞) min с использованием сингулярного разложения грамиана WU :
( )α1m2in {WU } z (0) = U[0,∞) min ≤ U[0,∞) = zT (0)WU z (0) 1 2 ≤ U[0,∞) max =
= α1m2ax {WU } z (0) ,

(19)

где αmin {WU } , αmax {WU } — минимальное и максимальное сингулярные числа грамиана WU
соответственно. Таким образом, контроль затрат на управление при полиномиальном экзогенном
воздействии может быть осуществлен с использованием следующего алгоритма. 1. Сформировать векторно-матричное описание технического объекта в форме (1). 2. Задать источник полиномиального экзогенного воздействия с помощью векторно-
матричного описания в форме (2). 3. Реализовать требования к качеству процессов в проектируемой системе в переходном и
установившемся режимах для источника ПЭВ (2) в виде структуры мод матрицы F, назначив ее
( )носителем матрицу Γ , задаваемую наблюдаемой парой матриц Γ, H : dim H = dim BT .

4. Сформировать матрицу K обратной связи по вектору состояния x(t) с использовани-
ем методов модального управления: K = HM −1 .
( )5. Сформировать матрицу прямых связей Kg = − CM Γ−1M −1B −1 .

6. Сформировать в форме (15) матрицу T в уравнении Сильвестра (10).
7. Решить уравнение Ляпунова (18) относительно грамиана WU .
8. Вычислить оценки затрат на управление системой в установившемся режиме при заданном экзогенном воздействии в форме мажорант и минорант этих затрат с использованием соотношения (19).

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9

Контроль энергетических затрат на управление

9

9. Провести анализ полученных оценок. Вернуться к п. 3 алгоритма в целях модифи-

кации структуры мод. По достижении минимальных значений мажоранты U[0,∞) max

завершить реализацию алгоритма.

Пример. Рассмотрим исходную систему с характеристическим полиномом Баттерфорта

второго порядка:

Φ1 ( s) = s2 +

2

(nω0 )2
⋅ nω0s +

(

nω0

)2

,

где п — радиус области расположения желаемых мод полинома, ω0 — характеристическая

частота системы.

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

отрицательной обратной связью:

Φ1′

(s)

=

s2

(nω0 )2
+ 2 ⋅ nω0s

.

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

E

=

⎡0 ⎣⎢0

1⎤ 0⎥⎦

;

P

=

[1

0];

z

(0)

=

⎡1⎤ ⎣⎢1⎥⎦

.

Таким образом, полиномиальное экзогенное воздействие представляет собой воздействие с

постоянными статической и кинетической составляющими: g (t ) = t +1.

Оценим затраты на воспроизведение ПЭВ g (t ) = 1t +1 при изменении n от 0,1 до 10
с шагом 0,1. Результаты моделирования рассматриваемой задачи в среде MatLab приведены на рис. 1. Анализ представленной зависимости показывает, что затраты на управление значительно снижаются при увеличении п.
U ∞ ⋅107
18

16

14

12 α1ma2x {U} z(0)
10

8

6

4 U[0,∞) 2 2

0 1 2 3 4 5 6 7 8 9 10 п

Рис. 1
Модифицируем систему, придав ей астатизм второго порядка:

Φ2

(s)

=

s2

2 +

⋅ nω0s + (nω0 )2 2 ⋅ nω0s + (nω0

)2

,

Φ′2 (s) =

2

⋅ nω0s + s2

( nω0

)2

.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9

10 Д. С. Бирюков, А. В. Ушаков
Оценим затраты на воспроизведение ПЭВ g (t ) = 1t +1 n от 0,1 до 1 с шагом 0,01.
Результаты моделирования представлены на рис. 2.
U ∞ ⋅10−8 2,5
2 α1ma2x {WU } z(0)

1,5
U[0,∞) 2 1

0,5

0
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 п
Рис. 2
Анализ зависимости показывает, что в данном случае затраты на управление изначально невелики вследствие того, что второй порядок астатизма обеспечивает в установившемся режиме воспроизведение полиномиального экзогенного воздействия первого порядка с нулевой ошибкой.
Таким образом, предложенный алгоритм формирования затрат на управление установившимся движением системы, порождаемым полиномиальным экзогенным воздействием, построен на использовании экстремальных элементов алгебраического спектра сингулярных чисел грамиана затрат на управление. Однако, по мнению авторов, „тонкое“ решение задачи требует, кроме того, исследования структуры элементов геометрического спектра левого сингулярного базиса этого грамиана.

СПИСОК ЛИТЕРАТУРЫ

1. Бирюков Д. C., Слита О. В., Ушаков А. В. Оценка затрат на управление в задаче обеспечения желаемой структуры мод и их робастности// Изв. вузов. Приборостроение. 2009. Т. 52, № 11. С. 32—37.

2. Бирюков Д. С., Ушаков А. В. Контроль затрат на управление при воспроизведении гармонических экзогенных воздействий: грамианный подход // Науч.-техн. вестн. СПбГУ ИТМО. 2011. № 2 (72). С. 117—123.

3. Кострикин А. И. Введение в алгебру. Основы алгебры. М.: Физматлит, 2004.

4. Ушаков А. В. Модальные оценки качества процессов в линейных многомерных системах при внешних конечномерных воздействиях // Автоматика и телемеханика. 1992. № 10. С. 72—82.

Сведения об авторах

Дмитрий Сергеевич Бирюков

— аспирант; Санкт-Петербургский национальный исследовательский

университет информационных технологий, механики и оптики, ка-

федра систем управления и информатики; E-mail: dbiryukov@list.ru

Анатолий Владимирович Ушаков — д-р техн. наук, профессор; Санкт-Петербургский национальный ис-

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

и оптики, кафедра систем управления и информатики;

E-mail: ushakov-AVG@yandex.ru

Рекомендована кафедрой систем управления и информатики

Поступила в редакцию 10.04.12 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 9