CC BY
41
УДК 681.51
ИНТЕРВАЛЬНЫЙ НАБЛЮДАТЕЛЬ ДЛЯ М ОДЕЛИ БИОЛОГИЧЕ СКОГО РЕАК ТОРА1
Т.А. Харьковская3, А.С. Кремлев3, Д.М. Сабирова3, Д.В. Ефимовь, Т. Раиссис
a Университет ИТМО, Санкт-Петербург, Россия, easymedia@mail.ru
ь Национальный институт исследований по информатике и автоматике, Лилль, Франция, efde@mail.ru с Центр исследований в области компьютерных наук и телекоммуникаций, Париж, Франция, tarek.raissi@cnam.fr Рассматривается метод построения интервального наблюдателя для нелинейных систем с параметрической неопределенностью. Задача синтеза интервального наблюдателя для систем с переменными параметрами сводится к следующему: если задан интервал неопределенности для значений состояния системы, ограничивающий начальные условия системы и множество допустимых значений для вектора неизвестных параметров и входов, то условие существования интервала оценок переменных состояния системы, в котором содержится фактическое значение состояния в данный момент времени, также должно выполняться на всем рассматриваемом временном интервале. Показаны условия построения интервальных наблюдателей для рассматриваемого класса систем: ограниченность состояния и входа, существование мажорирующей функции, задающей вектор неопределенностей системы, Липшицева непрерывность или ограниченность этой функции, существование коэффициента усиления наблюдателя с соответствующей матрицей Ляпунова. Основное условие построения подобного устройства оценки связано с кооперативностью динамики ошибки интервальной оценки. Рассматривается вопрос выбора индивидуальной матрицы усиления наблюдателя. Для обеспечения свойства кооперативности динамики ошибки интервальной оценки предлагается статическое преобразование координат. Результат работы метода продемонстрирован с помощью компьютерного моделирования системы биологического реактора. Возможными областями применения подобных систем интервального наблюдения являются области робастного управления, где предполагается наличие различного рода неопределенностей в динамике системы, биотехнологические и экологические системы и процессы, мехатроника, робототехника и др. Ключевые слова: интервальная оценка, наблюдатель, нелинейные системы, системы с переменными параметрами, параметрическая неопределенность, биореактор.
INTERVAL OBSERVER FOR A BIOLOGICAL REACTOR MODEL1 T.A. Kharkovskaya3, A.S. Kremlev3, D.M. Sabirova3, D.V. Efimovb, T. Raissic
3 ITMO University, Saint Petersburg, Russia, easymedia@mail.ru
b INRIA (Institut national de recherche en informatique et en automatique), Lille, France, efde@mail.ru c Centre for Research in Computer Science and Telecommunications (Cedric), Cnam, Paris, France, tarek.raissi@cnam.fr The method of an interval observer design for nonlinear systems with parametric uncertainties is considered. The interval observer synthesis problem for systems with varying parameters consists in the following: if there is the uncertainty restraint for the state values of the system, limiting the initial conditions of the system and the set of admissible values for the vector of unknown parameters and inputs, the existence condition of the interval estimates for the system state variables, which contains the actual state at a given time, must also be performed on all considered time segment. Conditions of the interval observers design for the considered class of systems are shown. They are: limitedness of the input and state, the existence of a majorizing function defining the uncertainties vector for the system, Lipschitz continuity or finiteness of this function, the existence of an observer gain with the suitable Lyapunov matrix. The main condition for design of such a device is cooperativity of the interval estimation error dynamics. An individual observer gain matrix selection problem is considered. In order to ensure the cooperativity of interval estimation error dynamics property a static transformation of coordinates is proposed. The proposed algorithm is demonstrated by computer modeling of the biological reactor. Possible applications of these interval estimation systems are the spheres of robust control, where the presence of various types of uncertainties in the system dynamics is assumed, biotechnology and environmental systems and processes, mechatronics and robotics, etc. Keywords: interval estimation, observer, nonlinear systems, parametric-varying systems, parametric uncertainty, bioreactor.
Введение
Биологический реактор (ферментер) представляет собой резервуар с мешалкой, сконструированный для культивирования клеток животных [1, 2]. Принцип работы биореактора достаточно прост, а его устройство и методики сочетания необходимых условий, наоборот, сложны. Существует несколько различных моделей биологических реакторов [1-4] в зависимости от действия, факторов среды, методов культивирования и др., многие из которых имеют в своей структуре параметрическую неопределенность и нелинейные зависимости, что усложняет процесс синтеза управления данными объектами. Улучшение в методах управления биореактором может привести к значительной экономии в биохимической промышленности и к повышению производительности. В настоящей работе рассматривается метод построения наблюдателя, оценивающего интервал, в котором находится фактическое значение неизмеряемого состояния системы, что может быть полезно для развития алгоритмов управления в биологических системах.
1 Работа выполнена при государственной финансовой поддержке ведущих университетов Российской Федерации (субсидия 074-U01) и гранта Президента Российской Федерации МК-464.2013.8.
1 This work was financially supported by the Government of the Russian Federation, Grant 074-U01, and by grant of the
President of the Russian Federation МК-464.2013.8._
Постановка задачи
Интервальные наблюдатели были введены, чтобы справиться с неопределенностями, которые, как известно, характерны для некоторых классов систем. Это одна из причин, почему данная техника становится все более актуальной [5, 6] и особенно успешной в области биотехнологических процессов, где можно говорить о достаточно больших неопределенностях [7, 8]. В таких случаях использование классических методов построения наблюдателей, оценки которых сходятся к точному значению состояния при отсутствии шума, невозможно. Однако возможно использование методов интервальной оценки, т.е. методов построения интервального наблюдателя, который вычисляет множества допустимых значений (интервал) для вектора состояний системы.
На данный момент существует несколько подходов к построению интервальных наблюдателей [6, 7, 9, 10]. Эта работа рассматривает и продолжает подход к построению интервальных наблюдателей, основанный на теории монотонных систем [6, 7, 11-13]. Одним из самых сложных допущений для построения интервального наблюдателя является требование кооперативности динамики ошибки интервальной оценки, которое было рассмотрено в работах [12, 14-17]. Тем не менее, кооперативность - достаточно специфическое свойство, и большинство систем не кооперативны. Показано, что при некоторых нестрогих условиях, применяя статическое преобразование координат, гурвицева матрица может быть преобразована в гурвицеву и мецлерову матрицу (матрица называется мецлеровой, если она имеет неотрицательные элементы вне главной диагонали). Матрица преобразования - это решение уравнения Сильвестра, и конструктивный порядок решения этого уравнения был приведен в [12].
Целью настоящей работы является продолжение развития данного подхода к построению интервальных наблюдателей для нелинейных систем с переменными параметрами. Результат продемонстрирован на примере компьютерного моделирования системы биореактора.
Общие сведения
Евклидова норма для вектора х£ R" будет обозначаться как |х|, и для измеримого и локального существенно ограниченного входа u: (R+ = {t£R: т >0}) символ ||u||[to,tl] обозначает его Lo
норму:
||«|[t0,t1] = esssup{|u(t)|,t £ [t0,ti]}, если t1 = +oo, тогда мы можем просто написать ||и||. Будем обозначать £о - множество всех входов u со свойством ||u|| < о. Обозначим последовательность целых чисел 1,..., к как 1, к. Символы In, Enxm и Ер обозначают единичную матрицу с размерностью пхп, матрицы, у которых все элементы равны 1, с размерностями пхти рх 1 соответственно. Для матрицы A£Rnxn вектор его собственных значений обозначается как Я(А), ||А||тах = таХ;=^7-=^|Ау| (поэлементно максимальная норма, не полумультипликативная) и ||А||2 = Jmax1=^^À1 (АГА) (индуцированная L2 норма матрицы), соотношение
||А||тах < ||А||2 <п||А||тах выполняется между этими нормами.
Для двух векторов хх,х2 £Е" или матриц АХ,А2 £Rnxn отношения хх <х2 и А± <А2 понимаются поэлементно. Соотношение Р<0 (Р >0) означает, что матрица Р £ Rnxn является отрицательно (положительно) определенной. Для матрицы A£Rmxn определено A = max{0, А}, А = А — А (аналогично для векторов) и обозначена матрица абсолютных значений всех элементов как |А| =А + А. Запись А £ M означает, что матрица А - мецлерова, т.е. имеет неотрицательные элементы вне главной диагонали.
Лемма 1. Пусть х£ R" будет вектором переменных, х <х<х, для некоторых х,"х £R", и А£ Rmxn будет постоянной матрицей, тогда
Ах - Ах <Ах<Ах - Ах.
Доказательство. Отметим, что Ах = (Â — А)х, что для х <х<х дает необходимые оценки.
Матрица A£Rnxn называется гурвицевой, если все ее собственные значения имеют отрицательную вещественную часть. Любое решение линейной системы
X = Ах + w(t),w :
с x£R" и мецлеровой матрицей А поэлементно неотрицательно для всех t> 0, при условии, что х(0) > 0. Такие динамические системы называются кооперативными [18, 19].
Лемма 2. Даны матрицы A£Rnxn, R£Rnxn и C£Rpxn. Если существует матрица L£Rnxp такая, что матрицы A — LC и R имеют одинаковые собственные значения, тогда R = S_1(A — LC)S, где матрица S£Rnxn при условии, что пары (А — LC, ех) и (R, е2) наблюдаемы для некоторых ех £Rlxn, е2 £Rlxn.
Этот результат был использован в [12] для построения интервальных наблюдателей для линейных стационарных систем с мецлеровой матрицей R.
Построение интервального наблюдателя
Будем рассматривать следующий вид системы, имеющей зависимость от неизвестных нестационарных параметров ^(t) 6 0:
fx(t) = A(t,y,u)x(t) +/(t,y,u,0
ly(t) = C(t,u)x(t) , (1)
где x(t) 6 1" - состояние; y(t) 6lp - выходная переменная; u(t) 6 Rm - известное входное воздействие; ^(t) 6 - неизвестное входное воздействие или неизвестные изменяющиеся параметры ^(t) 6 0Vt > 0, множество 0 известно. Отметим, что
х = А(Ох+ B(<;)u = Ах + fit, и, О, fit, и, О = [AQO- А]х + В(^)и.
Допущение 1. ||х|| <Х, ||uH <U и ||у|| <Y, границы X>0,U>0,Y> 0 заданы.
Допущение 2. Пусть х<х<х для некоторых х,х 6R", тогда f(t,x,x, u) < /(t,x, и, Q <
<f(t,x,x,и) для некоторых заданных /: ш2п+т+1 ^Rn, /: R2n+m+1 ^Rn и всех t> 0, ||u|| <U,Z, 6 0.
Допущение 3. Существует матричная функция L: Р : ^Rnxn, P(0 = Р(0Г >0
такая, что для всех t> 0 и ||u|| <U, ||y|| <Y:
Pit) + D(t,y,u)rP(t) + P(t)D(t,y,u) + P(t)2 +Q = 0, D(t, y, u) =A(t, y, u) — L(t, y, u)C(t, u), Q=Qr >0.
Допущение 2 означает, что если даны границы х, х состояния х, то значения нелинейной функции f заключены в интервале [/,/] для всех ^6 0. В допущении 3 представлен коэффициент усиления наблюдателя L(t, у, и), который обеспечивает устойчивость нестационарной матрицы D(t, у, и) с матрицей функции Ляпунова P(t), это допущение определяет условия устойчивости динамики оценки.
При этих допущениях, если существует матрица-усилитель L 6 Rnxp из допущения 3 такая, что матрица D(t, у, u) =А — LC является гурвицевой и мецлеровой, можно построить интервальный наблюдатель [6, 11] следующего вида:
х = A(t,y,u)x + /(t,x,x,u) + L(t,y,u)[y- C(t,u)x],
- I = _ (2)
x = A(t,y,u)x + /(t,x,x,u) + L(t,y,u)[y- C(t,u)x].
Теорема 1. [20, 21] Пусть выполнены допущения 1, 2 и 3 и матрица D(t, у, и) - мецлерова для всех t> 0 и ||u|| <U, ||у|| <Y. Пусть выполнено одно из следующих условий:
1. |/(t,x,x,u)| < +а>, |/(t,x,x,u)| < +с» для любых t> 0, ||u|| <U и всех x,x6l";
2. для любых t> 0, ||х|| <Х, ||u|| <U, ^ 6 9 и всех х,х 6К"
|/(t,x,u,£) -/(t,x,x,u)| +|/(t,x,x,u) -/(t,x,u,^)|2 <p|e|2 + р|ё|2 + a
для некоторых a6R+,p 6 и pi„ -Q + R = 0, R = RT >0.
Тогда в (1) и (2) переменные x(t), x(t) остаются ограниченными для всех t> 0 и х(0) <х(0) <х(0) обеспечивает соотношение X <х<х.
Доказательство. Рассмотрим ошибки интервального оценивания ё = х — х, е=х — х: e = D(t, у, u)e + /(t,х, u, £) -f(t,х,х, и)
е = D(t,y,u)e + /(t,x,x,u) -/(t,x,u,£). Согласно допущению 2, для мецлеровой матрицы D для всех t> 0 свойства fit, x(t),u(t),£) > >/(t,x(t),x(t),u(t)), /(t,x(t),u(t),^) <7(t,x(t),x(t),u(t)) и x(t) <x(t) <x(t) выполняются при условии, что х(0) <х(0) <х(0). Чтобы доказать, что переменные x(t), x(t) ограничены, рассмотрим производную функции Ляпунова V=erP(t)e + erP(t)e: V = er[P(t) + D(t,y,u)rP(t) + P(t)D(t,y,u)]e +
+er[P(t) + D(t,y,u)rP(t) + P(t)D(t,y,u)]e + 2erP(t) [/(t,x,u,^) -/(t,x,x,u)] +
+2erP(t)[7(t,x,x, u) —/(t,x, u, £)]. Согласно допущению 3, это уравнение может быть переписано следующим образом:
V<-erQe-erQe+|/(t,x,u,^) -/(t,x,x,u)|2 + |7(t,x,x,u) -/(t,x,u,^)|2.
Если первое условие теоремы верно, тогда элементы |/(t,x, u, — /(t,x,x, и)| и |/(t,x,x, и) —
-/(t,x,u,^)| ограничены для любых t> 0, ||х|| <Х, ||u|| <U,^ 6 9 и всех х,х 6 1". Таким образом,
ошибки е,е ограничены стандартными аргументами Ляпунова, и поэтому переменные х,х также ограничены (из допущения 1 состояние х ограничено). Если второе условие теоремы выполняется, то это неравенство принимает вид
V <-ërQë -erQe + а, что подразумевает ограниченность х,х теми же аргументами.
Результат теоремы 1 основан на довольно строгом допущении, что матрица D - мецлерова. Все остальные предположения довольно часто встречаются в теории оценивания. Для постоянной матрицы D это допущение снимается в лемме 2, где показано, что в условиях допущения 3 существует вещественная статическая матрица преобразования подобия S с S_1DS, являющаяся гурвицевой и мецлеровой. В нашем случае D(t, у, и) является нестационарной матрицей, расширение леммы 2 для этого случая показано ниже.
Лемма 3 [12, 21]. Пусть Ze£cRnxn - нестационарная матрица, удовлетворяющая интервальным ограничениям S = {Ze Rnxn:Za - Д <Z<Za +Д} для некоторых ZaT = Za £Rnxn и ДеЕ^хп. Если для некоторой постоянной ^e К и диагональной матрицы YeKnxn мецлерова матрица R= —Y имеет те же самые собственные значения, что и матрица Za, тогда существует ортогональная матрица SeR"xn такая, что матрицы SrDS - мецлеровы для всех Zeï при условии, что п||Д||тах Доказательство этой леммы приведено в [21].
Допущение 4. Пусть D(t,y,u)eï для всех t> 0, ||u| <U, ||у|| <Y, где ï = {De Rnxn:Da -Д < < D<Da +Д} для некоторых DaT = Da eR"xn и ДеЕ+хп. Пусть для некоторой постоянной ц>п||Д||тах и диагональной матрицы YeKnxn мецлерова матрица R= цЕ„ —Y имеет те же самые собственные значения, что и матрица Da.
При этом допущении существует ортогональная матрица Se Rnxn такая, что матрицы SrD(t, у, u)S - мецлеровы для всех D(t, у, u) eS. Введем новую переменную состояния z = Srx, тогда система (1) может быть переписана в новых координатах:
z(t) = SrA(t,y,u)Sz(t) + 9(t,z,u,Ç), где ф(£, z, u, Ç) = Sт f(t, Sz, u, Ç). Используя условия леммы 3, мы имеем следующие соотношения: x(t) <x(t) <x(t),
х = S+z-S"z,x = S+z-S"z, (3)
где z <z<z-это интервальные оценки для переменной z. Точно также из допущения 2 получаем:
ф(с,г,г,и) = S+T/(t,x,x,u) - S"T/(t,x,x,u) < ф(с,г,и,£) <S+T/(t,x,x,u) - S"T/(t,x,x,u) = = ф(£, z,z, u).
В новых координатах интервальный наблюдатель принимает вид, аналогичный (2): z = SrA(t,y,u)Sz + ф(£:,г,г,и) + SrL(t,y,u)[y- C(t,u)Sz],
- I _ I (4)
z = SrA(t,y,u)Sz + ф(с,г,г,и) + SrL(t,y,u)[y- C(t,u)Sz].
Теперь можно показать следующий расширенный вариант теоремы 1.
Теорема 2 [21, 22]. Пусть выполнены допущения 1, 2, 3 и 4. Пусть выполнено одно из следующих условий:
1. |/(t,x,x,u)| < +œ, |/(t,x,x,u)| < +œ для любых t> 0, ||u|| <U и всех x,xel";
2. для любых t> 0,||х|| <Х, ||u|| <U,Çe 0 и всех z,zel"
|ф(£,г,и,£) -ф(£,г,г,и)| +|ф(£,г,г,и) - ф(£, z, и, Ç)|2 <ß|z - z|2 + ß|z-z|2 + a
для некоторых ael+,ßeR+ и ßl„ — SrQS + R < 0, R = Rr >0.
Тогда в (4), (3) и (1) переменные x(t), x(t) ограничены для всех t> 0, и z(0) <z(0) <z(0) обеспечивает соотношение x(t) <x(t) <x(t).
Доказательство. Рассмотрим динамику ошибок интервального оценивания e = z-z, e = z-z:
é = SrD(t, y, u)Se + ф(£, z, u, Ç) -ф(£,х,х, u) ë = SrD(t, y, u)Së+9(t,z,z, u) - ф(£, z, u, Ç).
Далее доказательство проводится на основе рассуждений из доказательства теоремы 1 при помощи функции Ляпунова V=erSrP(t)Se + êrSrP(t)Sê [21].
Эта теорема предлагает интервальный наблюдатель для системы с переменными параметрами, явно опуская требование кооперативности матрицы D замкнутого контура. Предложено преобразование координат, обеспечивающее свойство кооперативности исходной некооперативной системе. Вдобавок,
снято допущение о том, что существует коэффициент усиления наблюдателя, который делает динамику ошибки оценки устойчивой и кооперативной. Коэффициент усиления наблюдателя должен, как обычно, обеспечить устойчивость ошибки наблюдения, а далее предлагается соответствующая замена координат.
Моделирование
Рассмотрим модель изотермического непрерывного биологического реактора, состоящую из следующих уравнений массового баланса клеток и ограниченного субстрата [2, 3]: = СТ(5)Х + _£ _£)
сИ V
где X - концентрация клеток; £ - концентрация субстрата в биореакторе; ц(5) - удельная скорость та; а(5) - удельная скорость потребления субстрата; 50 - концентрация субстрата в поступающем потоке; V - объем реактора; ^ - объемная скорость потока через биореактор.
Предполагается, что функция с(8) имеет вид [1, 2]
^ ; У(5) (Кт +5)7(5)' где цт - максимальная удельная скорость роста, Кт - константа насыщения субстрата Монода, тогда выражение выхода 7(5) имеет вид
количество сформированной биомассы
7(5) =-----=а5 + Ь,
количество потребляемого субстрата
что отражает увеличение выхода в ответ на увеличение концентрации субстрата 5. Рассматривается случай постоянного выхода при а = 0.
В дальнейшем мы предполагаем, что [1] скорость разбавления Б = ^ строго положительна и ограничена, т.е. £>т;п (с) <£>тах, скорость подачи 50 также ограничена, и каждая реакция включает по крайней мере один реагент, который не является ни катализатором, ни автокатализатором. Согласно [1, 4], переменная X и скорость разбавления Б ограничены на всем рассматриваемом интервале времени.
Коэффициент усиления наблюдателя был найден как Ь= [0,8309 —1,5726]т, что позволяет обеспечить устойчивость наблюдателя и кооперативность ошибки наблюдения.
Построенные границы вектора состояния (X, £) представлены на рис. 1, 2. Результаты показывают, что измерения состояний гарантированно находятся внутри оцененного интервала, а интервальная оценка сходится к области с шириной, зависящей от динамики ошибок измерений.
2о
*, с
Рис. 1. Результаты моделирования: изменение состояния X концентрации клеток (1) и его интервальная
оценка - верхняя (2) и нижняя (3) границы
3
Рис. 2. Результаты моделирования: изменение состояния концентрации субстрата S (1) и его интервальная оценка - верхняя (2) и нижняя (3) границы
Заключение
В работе показано, что интервальный наблюдатель позволяет получить область оценок переменных состояния системы, гарантированно содержащих фактическое значение состояния в данный момент времени. Рассмотрен вопрос выбора индивидуальной матрицы усиления наблюдателя, который должен гарантировать только устойчивость ошибки наблюдения. Для обеспечения свойства кооперативности динамики ошибки интервальной оценки предложено статическое преобразование координат, с помощью которого устойчивая система с переменными параметрами может быть превращена в устойчивую и кооперативную. Приводится доказательство теоремы об ограниченности траекторий полученной области на основе свойства кооперативности системы. Апробация метода проведена с помощью компьютерного моделирования системы биологического реактора.
Литература
1. Bastin G., Van Impe J.F. Nonlinear and adaptive control in biotecnology: a tutorial // European Journal of Control. 1995. V. 1. N 1. P. 37-53.
2. Ajbar A., Alhumaizi K. Dynamics of the Chemostat: A Bifurcation Theory Approach. CRC Press, 2011. 368 p.
3. Fossas E., Ros R.M., Fabregat J. Sliding mode control in a bioreactor model // Journal of Mathematical Chemistry. 2001. V. 30. N 2. P. 203-218.
4. Гордеева Ю.Л., Гордеев Л.С. Математическая модель непрерывного процесса в биореакторе с рециклом субстрата и биомассы // Вестник АГТУ. Серия: Управление, вычислительная техника и информатика. 2013. № 2. С. 9-18.
5. Mazenc F., Niculescu S.I., Bernard O. Interval observers for linear systems with delay // Proc. of the 48th IEEE conference on decision and control. Shanghai, China, 2009. P. 1860-1865.
6. Moisan M., Bernard O., Gouze J.-L. Near optimal interval observers bundle for uncertain bioreactors // Automatica. 2009. V. 45. N 1. P. 291-295.
7. Bernard O., Gouze J.-L. Closed loop observers bundle for uncertain biotechnological models // Journal of Process Control. 2004. V. 14. N 7. P. 765-774.
8. Moisan M., Bernard O. Interval observers for non monotone systems. Application to bioprocess models // IFAC Proceedings Volumes (IFAC-PapersOnline). 2005. V. 16. P. 43-48.
9. Jaulin L. Nonlinear bounded-error state estimation of continuous time systems // Automatica. 2002. V. 38. N 6. P. 1079-1082.
10. Kiefer M., Walter E. Guaranteed nonlinear state estimator for cooperative systems // Numerical Algorithms. 2004. V. 37. N 1-4. P. 187-198.
11. Raissi T., Videau G., Zolghadri A. Interval observers design for consistency checks of nonlinear continuous-time systems // Automatica. 2010. V. 46. N 3. P. 518-527.
12. Raissi T., Efimov D., Zolghadri A. Interval state estimation for a class of nonlinear systems // IEEE Transactions on Automatic Control. 2012. V. 57. N 1. Art. N 5983407. P. 260-265.
13. Efimov D., Fridman L.M., Raissi T., Zolghadri A., Seydou R. Interval estimation for LPV systems applying high order sliding mode techniques // Automatica. 2012. V. 48. N 9. P. 2365-2371.
14. Mazenc F., Bernard O. Interval observers for linear time-invariant systems with disturbances // Automatica.
2011. V. 47. N 1. P. 140-147.
15. Чеботарев С.Г., Кремлев А.С. Синтез интервального наблюдателя для линейной системы с переменными параметрами // Изв. вузов. Приборостроение. 2013. Т. 56. № 4. С. 42-46.
16. Chebotarev S., Efimov D., Raissi T., Zolghadri A. On Interval observer design for a class of continuous-time LPV systems // IFAC Proceedings Volumes (IFAC-PapersOnline), 2013. V. 9, part 1. P. 68-73.
17. Chebotarev S., Kremlev A. Analysis conditions on interval observer synthesis for linear systems with variable parameters // 18th International Conference on Methods and Models in Automation and Robotics, MMAR 2013. Mi^dzyzdroje, Poland, 2013. Art. N 6669939. P. 390-392.
18. Чеботарев С.Г., Кремлев А.С. Анализ линейных систем с переменными параметрами для синтеза интервальных наблюдателей // Научно-технический вестник информационных технологий, механики и оптики. 2012. № 6 (82). С. 50-53.
19. Smith H.L. Monotone dynamical systems: an introduction to the theory of competitive and cooperative systems. Providence: AMS, 1995. V. 41. 174 p.
20. Ефимов Д.В., Кремлев А. С., Харьковская Т.А., Чеботарев С.Г. Построение системы интервального оценивания для модели регуляции гормона тестостерона // Научно-технический вестник информационных технологий, механики и оптики. 2013. № 6 (88). C. 56-60.
21. Efimov Denis, Raissi Tarek, Chebotarev Stanislav, Zolghadri Ali. Interval state observer for nonlinear time varying systems // Automatica. 2013. V. 49. N 1. P. 200-205.
22. Efimov Denis V., Raissi Tarek, Chebotarev Stanislav, Zolghadri Ali. On set-membership observer design for a class of periodical time-varying systems // Proceedings of the IEEE Conference on Decision and Control.
2012. Art. N 6426474. P. 6767-6772.
Харьковская Татьяна
Александровна
Кремлев Артем Сергеевич
Сабирова Дина Мизхатовна
Ефимов Денис Валентинович
Раисси Тарек
Tatiana A. Kharkovskaya Artem S. Kremlev
Dina M. Sabirova Denis V. Efimov
Tarek Raissi
магистрант, Университет ИТМО, Санкт-Петербург, Россия, easymedia@mail.ru
кандидат технических наук, зам. декана, доцент, Университет ИТМО, Санкт-Петербург, Россия, kremlev_artem@mail.ru
магистрант, Университет ИТМО, Санкт-Петербург, Россия, dy.alastor@bk.ru
доктор технических наук, ответственный исследователь первого ранга, Национальный институт исследований по информатике и автоматике, Лилль, Франция, efde@mail.ru
кандидат технических наук, доцент, Центр исследований в области компьютерных наук и телекоммуникаций, Париж, Франция, tarek.raissi@cnam.fr
student, ITMO University, Saint Petersburg, Russia, easymedia@mail.ru Associate professor, PhD, ITMO University, Saint Petersburg, Russia, kremlev_artem@mail.ru
student, ITMO University, Saint Petersburg, Russia, dy.alastor@bk.ru INRIA (Institut national de recherche en informatique et en automatique), Principal investigator of the first rank, Lille, France, D.Sc., efde@mail.ru Associate professor, PhD, Paris, France, Centre for Research in Computer Science and Telecommunications (Cedric), Cnam, tarek.raissi@cnam.fr
Принято к печати 02.04.2014 Accepted 02.04.2014
