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

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

42 С. Г. Чеботарев, А. С. Кремлев
УДК 681.51
С. Г. ЧЕБОТАРЕВ, А. С. КРЕМЛЕВ
СИНТЕЗ ИНТЕРВАЛЬНОГО НАБЛЮДАТЕЛЯ ДЛЯ ЛИНЕЙНОЙ СИСТЕМЫ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ*
Проанализированы условия, при которых возможно построение интервального наблюдателя, позволяющего определить область оценок переменных состояния системы.
Ключевые слова: интервальный наблюдатель, мецлеровость, параметрическая неопределенность.
Введение. В теории автоматического управления важную роль играет идентификация параметров системы. Поэтому задача оценки неизмеряемого вектора состояния очень актуальна, и ее решение требуется во многих случаях [1—3].
В некоторых ситуациях, например при наличии в описании системы параметрических неопределенностей, невозможно использовать классические методы построения устройства оценки неизвестных неизмеряемых параметров исследуемой системы (наблюдателей), оценки которых сходятся к точному значению состояния при отсутствии шума. В связи с этим необходимо синтезировать интервальный наблюдатель, обеспечивающий оценивание множества допустимых значений вектора состояния системы.
Задача оценки неизмеряемого вектора состояний может быть решена различными способами. Говоря об интервальных методах оценки, можно отметить различные подходы к построению подобных наблюдателей [4—7]. В настоящей статье рассматривается подход к построению интервального наблюдателя, основанный на теории монотонных систем [4, 7]. Для нелинейных систем данный подход был расширен в работе [8] на случай использования LPVпредставления (Linear Parameter Varying) с известными минорной и мажорной матрицами, а в [9] — на случай наблюдаемых нелинейных систем.
Рассмотрим случай, когда неопределенности в описании системы являются детерминированными неизвестным функциями, зависящими от времени. Тогда, при определенных условиях, возможно оценить границы ненаблюдаемых переменных, используя „робастный“ метод оценки параметров. На выходе интервального наблюдателя получаются два вектора оценок: минимальных и максимальных значений для каждого элемента вектора состояний объекта.
Решение задачи возможно при соблюдении некоторых допущений. Одним из самых сложных при построении интервального наблюдателя является требование мецлеровости матрицы [9, 10].
Однако для построения и корректной работы интервального устройства оценки неизвестных неизмеряемых параметров исследуемой системы необходимо выполнять условие ограниченности траекторий, получаемых на выходе верхнего и нижнего наблюдателей. Очевидно, что требуется доказать ограниченность этих траекторий с помощью методов теории автоматического управления.
В связи с этим развитие методов интервальной оценки и алгоритмов построения интервальных наблюдателей для систем, обладающих параметрической неопределенностью, попрежнему остается актуальной задачей.
* Исследование выполнено при поддержке министерства образования и науки Российской Федерации, соглашение 14.В37.21.0421 „Разработка автономной бортовой системы навигации и управления многофункциональными мультиротационными летательными аппаратами“.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 4

Синтез интервального наблюдателя для линейной системы с переменными параметрами 43

Постановка задачи. При выполнении условия мецлеровости [11—13] для матриц, описывающих рассматриваемую систему, полученный интервал оценок переменных состояния системы содержит в себе реальные значения переменных состояния. Однако это не гарантирует ограниченности решений уравнений, описывающих наблюдатели.
Применительно к устойчивым линейным системам кроме условия мецлеровости матрицы, описывающие динамику наблюдателей, также должны удовлетворять условиям гурвицевости. Обеспечение и проверка выполнения указанных условий является основной задачей при построении интервального наблюдателя.
Основной результат. Рассмотрим устойчивую систему общего вида:

x = A( p (t )) x + d (t, y,u),⎫⎪

y = Cx,

⎬ ⎭⎪

p (t) ∈ R,

где x ∈ R+ , y ∈ R p , u ∈ Rm — это векторы переменных состояния, выход системы и управ-

ляющее воздействие соответственно; A( p (t )) — матрица динамики системы, зависящая от

неизвестного параметра p (t ) :

d (t, y,u) ≤ d (t, y,u) ≤ d (t, y,u).
Граничные значения параметра p известны:

(1)

A ≤ A( p(t)) ≤ A.

(2)

В рассматриваемом случае существует вектор коэффициентов усиления наблюдателя L

такой, что:

L : A( p (t)) − LC ∈ M ,

(3)

эта система обладает свойством гурвицевости. Отсюда следует, что системы

A − LC, A − LC ∈ M

(4)

также обладают свойством гурвицевости. Будем рассматривать случай, когда состояния системы неотрицательны, т.е. x(t) ≥ 0.

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

x (t ) ∈ Rn , требовались бы дополнительные действия для сохранения относительной ориен-

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

x = Ax + d (t, y,u) + L (Cx − y),⎫⎪

x

=

Ax

+

d

(t,

y,u)

+

L(

y



Cx),

⎬ ⎪⎭

где x, x — переменные состояния верхнего и нижнего наблюдателей соответственно.

(5)

Проверка условий мецлеровости. Как уже отмечалось, для построения интервального наблюдателя необходима мецлерова матрица. Проверим это условие, проанализировав уравнения ошибки для наблюдателей (5):

e = x − x,⎫ e = x − x,⎬⎭ где e , e — векторы ошибки для верхнего и нижнего наблюдателя соответственно.

(6)

Используя системы (5) и (6), найдем уравнения, описывающие динамку ошибки для рассматриваемых наблюдателей:

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 4

44 С. Г. Чеботарев, А. С. Кремлев

e = x − x,⎫ e = x − x.⎭⎬ Производные ошибок для системы (7) будут иметь следующий вид:
e = Ax + d (t, y,u) + L ( y − Cx ) − A( p (t )) x − d (t, y,u) =

(7)

= ⎣⎡ A( p (t )) − LC⎦⎤ e + [A − A( p (t ))]x + ⎡⎣d (t, y,u) − d (t, y,u)⎤⎦ , e = A( p(t)) x + d (t, y,u) − Ax − d (t, y,u) − L( y − Cx) =

(8)

= ⎣⎡ A( p (t )) − LC⎦⎤ e + [A( p (t )) − A]x + ⎡⎣d (t, y,u) − d (t, y,u)⎦⎤ .

(9)

Проанализируем полученные выражения. Слагаемое A( p (t )) − LC положительно, с

учетом исходных данных (3); слагаемые A − A( p (t )) и A( p (t )) − A в (8) и (9) также облада-

ют данным свойством (см. (2)); слагаемые d (t, y,u) − d (t, y,u) и d (t, y,u) − d (t, y,u) положи-

тельны исходя из соотношения (1). Учитывая, что над слагаемыми производится операция сложения, можно сделать вывод о том, что производные ошибки как для верхнего, так и для нижнего наблюдателя положительны. Тогда, если дополнительно указать, что начальное значение ошибок неотрицательно
e (0) ≥ 0; e (0) ≥ 0,
можно заключить, что для любого t
e (t) ≥ 0 и e (t) ≥ 0.

Проверка устойчивости границ интервала оценки. Представим систему наблюдателей (5) в матричном виде:

X

=

⎡ ⎢ ⎣

x x

⎤ ⎥ ⎦

;

X

=

⎡A ⎣⎢

− LC 0

Введем обозначения:

A

0⎤ − LC ⎥⎦

X

+

⎡L⎤ ⎣⎢ L ⎥⎦

y

+

⎡d ⎢⎣d

(t, (t,

y,u)⎤ y, u )⎥⎦

.

(10)

A*

=

⎡A ⎢⎣

− LC 0

A

0 − LC

⎤ ⎦⎥

;

Λ

=

⎡ ⎣⎢

L⎤ L⎥⎦

;

D

(t,

y,

u

)

=

⎡d ⎢⎣d

(t, (t,

y,u)⎤ y, u )⎥⎦

,

тогда выражение (10) с учетом (11) будет иметь вид

(11)

X = A* X + Λy + D (t, y,u).

(12)

Задача построения сводится к анализу (12) устойчивости и ограниченности решений. Выходной сигнал y неявно зависит от значения x :

y ≤ C ( x + x ).

Проанализируем следующую систему:
X = A* X + Λy, где y ≤ C ( x + x ) ,

для доказательства ее устойчивости рассмотрим функцию Ляпунова вида V = X T PX ,
где матрица P = PT > 0 , тогда ее производная будет вида:
( )V = X T A*T P + PA* X + 2X T PΛy ≤

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 4

Синтез интервального наблюдателя для линейной системы с переменными параметрами 45

( )≤ X T A*T P + PA* + P2 X + yT ΛT Λy ≤

( )≤ X T A*T P + PA* + P2 + 4 c 2 Λ2 I X .

Используя решение уравнения Риккати, можно найти такие значения Λ,C, A* , что для

рассматриваемой системы будут выполняться условия устойчивости.

Пример. Построим интервальный наблюдатель на примере линейной системы второго

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

x = A( p,t ) x + Bu(t),⎫

y = Cx,

⎬ ⎭

A(

p,

t

)

=

⎡ ⎢ ⎣

−4 p (t )

1⎤ −4⎥⎦

,

B

=

⎡1⎤ ⎣⎢0⎥⎦

,

C

= [0

1] , pmin = 3, pmax = 5 .

Значение параметра p изменяется следующим образом:

p (t ) = pmin + 0,5( pmax − pmin )(0, 5sin 2t + 0,5sin 0,38t + 1).

Выберем

L

=

⎡−0, 05⎤

⎢ ⎣

0,1

⎥ ⎦

.

Решением уравнения Риккати для заданных значений является

матрица P , собственные значения которой положительны:

σ( P) = [0, 0037; 0, 0042; 0, 0122; 0,0191].

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

x 2

x1 x2

x1 1 x2

x2 x1 0 5 10 t, с
Полученные результаты позволяют сделать вывод о возможности эффективного использования предложенного подхода при построении интервального наблюдателя для линейных систем с переменными параметрами. Целесообразно дальнейшее его развитие для применения в системах с разными комбинациями свойств и входных данных.
СПИСОК ЛИТЕРАТУРЫ
1. Meurer T., Graichen K., Gilles E.-D. Control and Observer Design for Nonlinear Finite and Infinite Dimensional Systems // Lecture Notes in Control and Information Sciences. Springer, 2005. Vol. 322. P. 422.
2. Fossen T. I., Nijmeijer H. New Directions in Nonlinear Observer Design. Springer, 1999. 525 p.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 4

46 С. Г. Чеботарев, А. С. Кремлев

3. Nonlinear Observers and Applications. Lecture Notes in Control and Information Sciences / Ed. by G. Besançon Springer, 2007. Vol. 363. 224 p.

4. Bernard O., Gouzé J. L. Closed loop observers bundle for uncertain biotechnological models // J. Process Control. 2004. Vol. 14. P. 765—774.

5. Jaulin L. Nonlinear bounded-error state estimation of continuous time systems // Automatica. 2002. Vol. 38, N 2. P. 1079—1082.

6. Kieffer M., Walter E. Guaranteed nonlinear state estimator for cooperative systems // Numerical Algorithms. 2004. Vol. 37. P. 187—198.

7. Moisan M., Bernard O., Gouzé J. L. Near optimal interval observers bundle for uncertain bio-reactors // Automatica. 2009. Vol. 45, N 1. P. 291—295.

8. Raïssi T., Videau G., Zolghadri A. Interval observers design for consistency checks of nonlinear continuous-time systems // Automatica. 2010. Vol. 46, N 3. P. 518—527.

9. Raïssi T., Efimov D., Zolghadri A. Interval state estimation for a class of nonlinear systems // IEEE Trans. on Automatic Control. 2012. Vol. 57, N 1. P. 260—265.

10. Mazenc F., Bernard O. Interval observers for linear time-invariant systems with disturbances // Automatica. 2011. Vol. 47, N 1. P. 140—147.

11. Smith H. L. Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems of Surveys and Monographs // AMS. Providence, 1995. Vol. 41. P. 174.

12. Efimov D., Raïssi T., Chebotarev S., Zolghadri A. Interval State Observer for Nonlinear Time Varying Systems // Automatica. 2013. Vol. 49(1). P. 200—205.

13. Чеботарев С. Г., Кремлев А. С. Анализ линейных систем с переменными параметрами для синтеза интервальных наблюдателей // Науч.-техн. вестн. информационных технологий, механики и оптики. 2012. № 06(82). С. 50—53.

Станислав Геннадьевич Чеботарев Артем Сергеевич Кремлев

Сведения об авторах — аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: freest5@gmail.com — канд. техн. наук; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра систем управления и информатики; доцент; E-mail: kremlev_artem@mail.ru

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

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

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 4