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

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

АНАЛИЗ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ …
УДК 681.51
АНАЛИЗ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ ДЛЯ СИНТЕЗА ИНТЕРВАЛЬНЫХ НАБЛЮДАТЕЛЕЙ1
С.Г. Чеботарев, А.С. Кремлев
Анализируются условия синтеза интервальных наблюдателей для линейных систем с неизвестными нестационарными параметрами. В англоязычной литературе этот класс систем получил наименование Linear Parameter-Varying Systems (LPV systems). Рассматриваются возможные варианты соблюдения свойств кооперативности системы, проводится анализ условий, при которых для данной комбинации свойств системы возможно построение интервального наблюдателя, позволяющего получить область оценок переменных состояния системы, гарантированно содержащую фактическое значение состояния в данный момент времени. Ключевые слова: интервальные наблюдатели, кооперативность, мецлеровость, параметрическая неопределенность.
Введение Проблема оценки неизмеряемого вектора состояния очень актуальна, и ее решение требуется во многих случаях [1–3]. В некоторых ситуациях, в связи с наличием неопределенности (параметрической и/или сигнальной), построение обычного наблюдателя, сходящегося в случае отсутствия шума к идеальному значению состояния, невозможно. Тем не менее, интервальная оценка остается возможной. Под интервальной оценкой здесь понимается работа интервального наблюдателя, который гарантирует оценивание множества допустимых значений для вектора состояний системы. Размер рассчитанного множества должен быть пропорционален неопределенности модели объекта. Неопределенности рассматриваются как детерминированные, но неизвестные функции времени. С этими ограничениями на параметры можно оценить границы ненаблюдаемых переменных: это так называемый «робастный» метод. Обычно интервальные наблюдатели генерируют два вектора оценок – минимальных и максимальных значений для каждого элемента вектора состояний объекта. Есть несколько подходов к построению интервальных наблюдателей [4, 5–7]. В настоящей работе рассматриваются ситуация построения интервального наблюдателя, основанного на теории монотонных систем [4, 7]. Этот подход был расширен в работе [8] на случай нелинейных систем с использованием LPV-представления с известными минорной и мажорной матрицами, а в работе [9] – для наблюдаемых нелинейных систем. Одно из самых сложных допущений для построения интервального наблюдателя, связанное с кооперативностью динамики ошибки интервальной оценки, было решено в [9, 10]. Известно, что построение интервального наблюдателя возможно, если динамика ошибки оценки обладает свойствами кооперативности и устойчивости. Было показано, что при некоторых мягких условиях, применяя подобное преобразование координат, гурвицеву матрицу можно преобразовать в гурвицеву и мецлерову матрицу (кооперативную матрицу). Матрица преобразования – это решение уравнения Сильвестра, и конструктивный порядок решения этого уравнения был приведен в [9].
1 Работа выполнена при финансовой поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009–2013 годы (соглашение № 14.B37.21.1928).
50 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 6 (82)

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

Работа организована следующим образом. Некоторые основные факты из теории интервального оценивания приведены в разделе «Общие сведения». Анализ систем и условия построения интервального наблюдателя для различных комбинаций их свойств описаны в разделе «Основной результат».

Общие сведения

Для матрицы A  Rmn определим A  max0, A, A  A  A . Запись A  M означает, что мат-
рица A – мецлерова, т.е. имеет неотрицательные элементы вне главной диагонали. Лемма 1. Пусть x  Rn будет вектором переменных, x  x  x для некоторых x, x  Rn , и
A  Rmn будет постоянной матрицей, тогда Ax  Ax  Ax  Ax  Ax .
 Доказательство. Отметим, что Ax  A  A x , что для x  x  x дает необходимые оценки.
Напомним, что матрица A  Rmn называется гурвицевой, если все ее собственные значения имеют отрицательную вещественную часть. Любое решение линейной системы
x  Ax  (t),  : R  Rn , с x  Rn и мецлеровой матрицей A поэлементно неотрицательно для всех t0 при условии, что x(0)0 [11]. Такие динамические системы называются кооперативными (монотонными) [11].
Лемма 2. Даны матрицы A  Rmn , R  Rmn и C  R pn . Если существует матрица L  Rn p такая, что матрицы A–LC и R имеют одинаковые собственные значения, тогда R=S-1(A–LC)S, где матрица S  Rnn при условии, что пары (A–LC, e1) и (R, e2) наблюдаемы для некоторых e1R1n , e2R1n. Этот результат был использован в [9] для построения интервальных наблюдателей для линейных стационарных систем с мецлеровой матрицей R (основная трудность состоит в доказательстве существования вещественной матрицы S и обеспечении конструктивного метода ее расчета).
Основной результат

Рассматривается следующий вид линейной системы, имеющей зависимость от неизвестных неста-

ционарных параметров pƤ

x  A(p,t)x  d (t, y,u),

 

y  Cx,

(1)

где xRn, uRm, yRp – это состояние, вход и выход системы (1) соответственно, x Ƥ Rq – это вектор неизвестных сигналов или параметров. В настоящей работе рассматривается случай без шума в кана-

ле измерения, предлагаемый результат может быть расширен на случай с шумом, это расширение опу-

щено для краткости изложения.

Допущение 1. Состояние x и вход u ограничены – стандартные допущения в теории оценивания.

Допущение 2. Если даны границы x, x состояния x, то значения функции d(t, y, u) заключены в

интервале d (t, y,u), d (t, y,u) .
Допущение 3. Существует коэффициент усиления наблюдателя L, который обеспечивает устойчивость нестационарной матрицы динамики наблюдателя с матрицей P(t) функции Ляпунова V; это допущение определяет условия устойчивости динамики оценки:
V  eT P(t)e  eT P(t)e.
Множество работ по построению интервального наблюдателя [4, 7–9] посвящено случаю с постоянной матрицей A (или в результате некоторых преобразований ошибка оценки может быть представлена в форме с постоянной матрицей A), и далее может быть найден такой коэффициент усиления наблюдателя L, что A–LC является гурвицевой и мецлеровой. В дальнейших работах [12] снимается предположение о том, что существует коэффициент усиления наблюдателя, который делает динамику ошибки оценки устойчивой и кооперативной. Коэффициент усиления наблюдателя должен, как обычно, обеспечить устойчивость ошибки наблюдения, а далее предлагается статическое преобразование координат, которое обеспечивает требуемое свойство кооперативности.
Рассмотриваемый метод построения интервального наблюдателя для линейных систем с переменными параметрами тесно связан с теорией дифференциальных неравенств и кооперативными системами. Основное условие построения подобного устройства оценки, работающего в условиях интервальности и неизвестного характера изменения параметров, от которых зависит система, связано с кооперативностью динамики ошибки интервальной оценки. Далее рассматриваются возможные варианты выполнения данного ограничения, на основе чего строится дальнейший анализ.

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 6 (82)

51

АНАЛИЗ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ …
1. (t)0; A (p, t)M. 2. x(t)0; A (p, t)M; L: (A (p, t) –LC )M. 3. x(t)0; A (p, t) M; L: (A (p, t) –LC )M; S: R= S-1(A (p, t) –LC) S M. 4. x(t)  n ; A (p, t)M; X: x*=x–X0.
5. x(t)  n ; A(p,t)  A  M ;A : Ax  d (t, x)  Ax .
6. x(t)  n ; A(p,t)  M ;L :  A(p,t)  LC  M ;X : x*  x  X  0.
7. x(t)  n ; A(p,t)  A  M ;L :  A(p,t)  LC  M ; A : Ax  d (t, x)  Ax .
8. x(t)  n ; A(p,t)  M ; L : A(p,t)  LC  M ;S : R  S1 A(p,t)  LCS  M .
  S1x; : *      0.
9. x(t)  n ; A(p,t)  M ;L :  A(p,t)  LC  M ;S : R  S1 A(p,t)  LCS  M ;
   R  const;   S1x;  S1  S1 A S  S  A* : A *   d (t, )  A * .
Отдельно рассматривается ситуация, когда состояния системы неотрицательны (п.п. 1–3), и ситуация, когда состояния знакопеременны (п.п. 4–9). Такое разделение обусловлено различным характером поведения границ полученного интервала оценок состояния. Это видно, если записать уравнение ошибки для рассматриваемой системы.
Пусть x, x обозначают оценки состояния системы для нижней и верхней границ интервала соответственно. Тогда производные оценок состояния системы для каждой из двух границ примут вид
x  Ax  d (t, y,u)  LCx  x, x  Ax  d (t, y,u)  LC x  x.
Запишем уравнение ошибки для нижней границы интервала как разность между состоянием и его оценкой:
e  x  x, e  x  x  A(p,t)x  d (t, y,u)  Ax  d (t, y,u)  LCx  x.
Для верхней границы интервала уравнение ошибки записывается аналогично. После преобразований динамика ошибки принимает следующий вид:
e  A  LC e  A(p,t)  A x.
Исходя из вида уравнения, можно сделать вывод о том, что второе слагаемое оказывает влияние на значение ошибки, следствием чего является возможность смены относительной ориентации границ интервала оценок состояния системы, если состояние меняет знак. Следовательно, системы с неотрицательными и знакопеременными значениями состояния стоит рассматривать отдельно и во втором случае предпринимать действия по исключению влияния данного фактора. Это можно осуществить двумя способами – сдвигом координат в область с положительными значениями состояний (п.п. 4, 6, 8) или записью уравнений, описывающих динамику наблюдателей, в другом виде, с использованием постоянной матрицы А и матрицы А, состоящей из неотрицательных элементов (п.п. 5, 7, 9).
Как уже отмечалось, необходимым условием построения интервального наблюдателя является особый вид матрицы, описывающей динамику объекта и наблюдателей. Они должны обладать свойством мецлеровости. Кроме того, так как на данном этапе рассматриваются устойчивые линейные системы, эти матрицы также должны удовлетворять условиям гурвицевости. Здесь можно отметить три варианта: 1. система сразу обладает необходимыми свойствами мецлеровости и гурвицевости (п.п. 1, 4, 5); подоб-
ный вариант встречается крайне редко; 2. может быть найден такой коэффициент усиления наблюдателя L, что А – LC является гурвицевой и
мецлеровой (п.п. 2, 6, 7); 3. коэффициент усиления наблюдателя L должен обеспечить лишь устойчивость ошибки наблюдения, а
далее предлагается статическое преобразование координат, которое обеспечивает требуемое свойство кооперативности (п.п. 3, 8, 9).
Таким образом, при анализе условий синтеза интервальных наблюдателей для линейных систем с неизвестными нестационарными параметрами необходимо принять во внимание знакопеременность значений состояний рассматриваемой системы и способ получения матрицы, описывающей динамику наблюдателей, в мецлеровом виде. Предполагается также устойчивость системы и обеспечение устойчивости траекторий наблюдателей с помощью выбора коэффициентов усиления L. Все возможные комбинации указанных условий приведены в п.п. 1–9. Следовательно, для любой системы из рассматриваемого класса можно построить интервальный наблюдатель, позволяющий получить область оценок переменных состояния системы, гарантированно содержащую фактическое значение состояния в данный момент времени, при соблюдении определенных условий, которые были проанализированы в данной работе.
52 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 6 (82)

А.Ф. Ломакин, С.И. Школьный

Заключение

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

Литература

1. Meurer T., Graichen K., Gilles E.-D. (Eds). Control and Observer Design for Nonlinear Finite and Infinite Di-

mensional Systems // Lecture Notes in Control and Information Sciences. – Springer, 2005. – V. 322. – P. 422.

2. Fossen T.I., Nijmeijer H. New Directions in Nonlinear Observer Design. – Springer, 1999. – 525 p.

3. Nonlinear Observers and Applications. Lecture Notes in Control and Information Sciences. – Besançon G.

(Ed.), Springer, 2007. – V. 363. – 224 p.

4. Bernard O., Gouzé J.L. Closed loop observers bundle for uncertain biotechnological models // J. Process

Control. – 2004. – V. 14. – P. 765–774.

5. Jaulin L. Nonlinear bounded-error state estimation of continuous time systems // Automatica. – 2002. –

V. 38. – № 2. – P. 1079–1082.

6. Kieffer M., Walter E. Guaranteed nonlinear state estimator for cooperative systems // Numerical Algorithms.

– 2004. – V. 37. – P. 187−198.

7. Moisan M., Bernard O., Gouzé J.L. Near optimal interval observers bundle for uncertain bio-reactors // Au-

tomatica. – 2009. – V. 45. – № 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. – V. 46 – № 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. – V. 57. – № 1. – P. 260–265.

10. Mazenc F., Bernard O. Interval observers for linear time-invariant systems with disturbances // Automatica. –

2011. – V. 47. – № 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. – V. 41. – 174 p.

12. Efimov D., Raïssi T., Chebotarev S., Zolghadri A. Interval State Observer for Nonlinear Time Varying Sys-

tems // Automatica.

– 2012 [Электронный ресурс]. – Режим доступа:

http://www.sciencedirect.com/science/article/pii/S0005109812004761, своб., Яз. англ. (дата

обращения 02.11.2012).

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

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

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 6 (82)

53