Субоптимальное управление линейным нестационарным объектом с запаздыванием по состоянию Текст научной статьи по специальности «Математика»

DOI: 10.24143/2072-9502-2019-3-123-132 УДК 519.7

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

Е. Р. Галяув, Е. В. Франгулова

Астраханский государственный технический университет, Астрахань, Российская Федерация

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

Ключевые слова: робастное управление, субоптимальное управление, внешнее возмущающее воздействие, неопределенный нестационарный линейный объект, интегральный критерий качества, наблюдатель, функция Ляпунова.

Для цитирования: Галяув Е. Р., Франгулова Е. В. Субоптимальное управление линейным нестационарным объектом с запаздыванием по состоянию // Вестник Астраханского государственного технического университета. Серия: Управление, вычислительная техника и информатика. 2019. № 3. С. 123-132. DOI: 10.24143/2072-9502-2019-3-123-132.

Введение

Синтез систем высокой точности в условиях неопределенности является классической проблемой в теории управления. Актуальной является проблема разработки алгоритмов управления, обеспечивающих не только устойчивость системы в этих условиях, но и оптимальное ее функционирование по некоторым критериям качества. Решению этих задач посвящено множество работ. На данный момент достаточно полно и систематизированно робастная теория изложена в [1]. Так, в условиях полной определенности в теории оптимального управления предложен -подход. Однако при наличии произвольных внешних возмущений, невозможности точно определить параметры модели объекта, изменении динамических свойств системы в процессе функционирования оптимальные системы, синтезированные по квадратичному критерию качества, часто теряют работоспособность. Большими возможностями в этих условиях обладает робастный подход к построению систем управления. К настоящему времени получено достаточно много алгоритмов построения субоптимальных систем для линейных и нелинейных объектов, подверженных параметрическим и внешним возмущениям: -оптимизация, -оптимизация, ц-синтез, ХМ1-подход и др. Следует отметить, что большинство имеющихся методов синтеза предназначены для стационарных систем. Однако практика в изобилии доставляет объекты управления, которые описываются параметрически неопределенными нестационарными с запаздыванием по состоянию дифференциальными уравнениями.

Особое внимание на сегодняшний день уделяется робастному и робастно-субоптимальному управлению системами, когда полный вектор состояния недоступен измерению, а измеряется только вектор выходных переменных [2-7]. Компенсация внешних возмущений путем объединения квазиинвариантной стратегии с классической предложена в [2]. Одним из распространенных подходов к задаче компенсации внешних детерминированных возмущений является метод внутренней модели [3, 4]. Робастные системы с компенсацией возмущений, построенные на основе их оценок, рассмотрены в [6, 7].

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

В предложенной работе решается задача построения робастного субоптимального регулятора для параметрически неопределенного линейного нестационарного объекта, подверженного воздействию внешних ограниченных возмущений. Цель управления - минимизация интегрального критерия качества с малой погрешностью отклонения от номинального значения. Синтез системы управления основан на результатах [8]. Компенсация неопределенностей базируется на подходе, предложенном в [6]. Результаты репрезентативного моделирования иллюстрируют работоспособность предложенной системы управления.

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

Рассматривается линейный нестационарный объект управления, динамические процессы в котором описываются в пространстве состояний уравнением

x(t) = A(t)x(t) + C(t)x(t-t(t)) + B(t)u(t) + D(t) f (t) , y (t) = Lx(t) , x(0) = x0, (1)

где функции x (t)e Rn - вектор состояния, y (t)e R - управляемый выход объекта, u (t)e R -управляющее воздействие, f (t) e R - ограниченное внешнее возмущающее воздействие,

x(t) e R+ - ограниченное неизвестное время запаздывания; —< 1; A (t) e RKXK, B (t) e Rnx1 ,

0 ... 0

C (t)e RKXK и D (t)e Rnx1 - неизвестные функциональные матрицы, C (t) =

< (t) - < (t)

L = [1; 0; —; 0], x0 - начальные условия.

Предположения

1. Коэффициенты матриц А(0, В{(), С(0, и D(t) зависят от некоторого вектора неизвестных параметров ^ е 3, где Е - известное ограниченное множество.

2. Пара (А(0, В(ф - управляема, а пара (А^), L) - наблюдаема.

3. Измерению доступны выходной сигнал у(0 и управляющее воздействие и(1). Использование производных этих величин в системе управления не допускается.

4. Выполнены условия структурного согласования А(t) = Ан + ВнуТ (t) ,

В (t) = Вн + Вн3(t), С (г) = ВнсТ (t) и D (t) = Внk (t) , где Ан е Михи, Вн е Мих1 - известные номинальные матрицы; Ан - гурвицева матрица в форме Фробениуса; )е М, c(t)е М -ограниченные неизвестные векторы; S(t) е М+, k (t) е М - ограниченные неизвестные функции.

5. Для любого фиксированного параметра t объект управления (1) минимально-фазовый. Требуется получить алгоритм функционирования системы управления, который обеспечит минимум функционалу

J = f[( У2 (t) q+u2 (t) r )] —t, (2)

0

где q > 0, r > 0 - весовые коэффициенты, выбираемые разработчиком.

со

Метод решения

Чтобы получить субоптимальную систему управления для нестационарного объекта (1), управление н(0 формируется в виде суммы двух составляющих [8]:

и ^) = и0 ^) + ик ^) .

Сигнал ик(0 необходим для компенсации неопределенностей в объекте (1). Слагаемое и0(0 определяет оптимальный закон управления, обеспечивающий минимум функционалу качества (2).

В силу условия 4 Предположений перепишем систему (1):

х (t ) = Ан х ^ ) + Вни (t) + Вн ф( х, и, t); у (t ) = Ьх (t), (3)

где функция ф(х, и, ^ = ({)х(^ + ат (t)х(^-т(+ и(^ + k(t) f (t) содержит все неопределенности объекта (1) и слагаемое с запаздыванием; Ь = [1; 0; ...; 0].

Рассмотрим систему (3) в предположении, что в ней отсутствуют неопределенности, т. е. при ф(х, и, t) = 0 . Для полученного номинального объекта

хн (t) = Ат хн (t) + ВНи0 (t) ; Ун (t) = Ьхн (t) (4)

рассмотрим критерий качества

да

J = |[хт (t)Мх(t) + и02 (t)г]dt (5)

0

с матрицей М = qЬ' Ь, преобразованный из (2) подстановкой второго уравнения (1).

Согласно [8] оптимальный закон управления и0 ^) в форме обратной связи по состоянию, гарантирующий минимум функционалу (5), определяется следующим образом:

и ^ ) = - К0 х (t), (6)

где матрица К0 = г 1В1тН, где Н = Нт > 0 является решением матричного уравнения Лурье - Риккати А1аН + НАн - г 1НВтВттН = -М; весовой коэффициент г > 0 .

Однако в силу того, что вектор состояния х(0 недоступен измерению, построим наблюдатели

х ^) = Ат х (t) + Вни (^ + S (у (t)-у (t)) + В^ (t) ; у (t) = Ьх (t). (7)

Здесь х (t) е Мя - оценка вектора состояния х(); v(t) - сигнал управления, компенсирующий влияние возмущений на процесс наблюдения; S - матрица коэффициентов усиления, такая, что А1 = Ат - БЬ - гурвицева; Ь = [1; 0; ...; 0].

Тогда оптимальный сигнал сформируем в виде

и (t) = -К0х (t). (8)

Но для качественной оценки вектора х(0 требуется скомпенсировать влияние параметрических и внешних возмущений на эти оценки. Поэтому необходимо выполнить следующие преобразования. Перепишем систему (3), добавляя и вычитая в правой ее части оптимальное управление (8):

х(^ = Атх+ Втик (t) + Втф(х, и, t)-ВтК0х (t); у (^ = Ьх(t). (9)

Прибавим и вычтем теперь значения (6) в (9), в результате получаем

х[{) = А,х+ Втик (t) + Втф(х, и, ^ + ВтК0ф); у (t) = Ьх(10)

Матрица А0 = Ат - ВтК0, с ^) = х (t) - х ^) - вектор ошибки, полученный вычитанием уравнения (7) из (3):

С (t) = Агф) + Втф(х, и, t)-Bтv (t); е (t) = Lс(t). (11)

Выберем вспомогательный контур [6]

С V (t) = Ас. (0-Вт V (0; ?у (0) = с(0); е (t) = (t), (12)

где qv (t) е Rn - вектор состояния, который позволяет выделить в отдельный сигнал все нежелательные воздействия. Принимая во внимание уравнения (11) и (12), составим уравнение рассогласования i (t) = q (t) — qv (t) , которое в переменных «вход-выход» запишется как

Q0 (p)S(t)= R0 (p)ф(X u, t) .

Здесь s(t) = e(t) — ev (t) ; Q0 (p) и R0 (p) - линейные дифференциальные операторы, коэффициентами которых являются коэффициенты полиномов Q0 (Я) = det (XI — A0) и R0 (Я) ;

degQ0 (X) = n ; degR0 (X) = m ; y = n — m > 1, p = d ~ оператор дифференцирования. Тогда

управляющие воздействия uk (t) и v (t) будем формировать следующим образом:

v(t) = |Мs(t) = ф(X, u, t); Uk (t) = —v (t).

Однако в силу Предположения 3 производные сигнала s(t) недоступны измерению. Поэтому выполним «операторное деление»

ом=T (p), aqq ( p ) R0 (p Г 0 () R0 (p)'

где T0 (p) = g0pY + g1 p7-1 + g2p7—2 +... + gY; degAQ0 (p) = m — 1. И, чтобы скомпенсировать все неопределенности системы (1), сигнал управления v (t) зададим в виде

v «=gTi(t)+ws(t )=*X- * 0; (i3)

Uk (t ) = —v (t) ,

где g - вектор, компонентами которого являются коэффициенты оператора T0(p), записанные в обратном порядке; ф(х, u, t) - оценка сигнала ф(х, u, t), содержащего параметрические и внешние

неопределенности исходного объекта; i(t) - оценка вектора 0T(t) = |s(t),s'(t),s''(t), ..., s(Y)(t)J,

полученная с помощью наблюдателя [9]

i(t) = G0I(t) + F0 (s(t) —s(t)); s(t) = Д(t). (14)

Здесь i (t) е ; G0 =

0 1Y 0 0

; = [-/ОмЛ-/1М-2, •••, -/М-7-1 ]; числа /0,/х, ...,/ выбираются так, чтобы матрица G = G0 - FL была гурвицевой; ¥Т =[ /о, /1, ..., /]; ц - малое положительное число; в (t) - оценка переменной в ^).

Для оценки точности наблюдения введем вектор отклонений

) = Г-1 (| (t )-0( t))

с матрицей Г-1 = diag{мУ-1, ..., м, 1}. Продифференцировав п(t) по времени с учетом (14), составим уравнение для нормированных отклонений:

1

П (t) = eW(t) .

П (г) = - г) + Ь 8(у) (г); д( г) = 8 (г) -8(г) = цу-1Ьп(г), ц

где Ь = [0, 0, ..., 1]. Преобразуем последнее уравнение в эквивалентное относительно выхода Л(г):

П (г) = - Gn(t) + Ь8'( г); Л(г) = м^-1^! (г) = цу-1Ьп(г), (15)

где п(г)е ®7; ЬТ =[1, 0, ..., 0]. Последние два уравнения эквивалентны относительно переменных п1 (г) и п1 (г), т. к. они являются различными формами записи одного уравнения

р1 +—р1,-1 + ...+ — * - М

Тогда, принимая во внимание уравнения закона управления (8), наблюдателя (14), уравнение для нормированных отклонений (15) и сигнал управления (13), уравнения (10) и (11) преобразуем к виду

х (г) = 4 х (г) + Ьц^У Л (г) + £н ^ (г); у (г) = Ьх (г);

— (16) & (г ) = Л<( г)-Ьц^Т Л (г); е (г ) = Ь<(г),

где Л (г) = [п1 (г), п1 (г), •••, П((у)(г.

Утверждение. Пусть выполнены условия предположений для объекта управления (1). Тогда существует число ц0 > 0, такое, что при ц < ц0 система управления (8), (12), (13)-(16) диссипативна и минимизирует функционал качества (2) с достаточно малой погрешностью отклонения от значения для номинального объекта.

Доказательство утверждения. Для доказательства утверждения введем новую переменную ю (г) е М, равную разности фазовых векторов (16) и (4), и рассмотрим систему

ю (г) = д)Ю(г) + Ь-Г1 gTл(0; ю(0) = 0; -п (г ) = Gn( г ) + ц2Ь 8'(г),

где ц1 = ц2 = ц.

Воспользуемся леммой [10].

Лемма [10]. Если система описывается уравнением х = f (х, ц1, ц2) , х е М™ , где f (х, ц1, ц2) - непрерывная функция, липшицева по х, и при ц2 = 0 имеет ограниченную замкнутую область диссипативности = {х | F (х) < С}, где ^х) - положительно определенная, непрерывная, кусочно-гладкая функция, то существует ц0 > 0 , такое, что при ц1 < ц0 и ц2 < ц0 исходная система имеет ту же область диссипативности П1, если для некоторых чисел С и ц1 при ц2 = 0 выполнено условие

sup

к

))dF ( х )

дх

f (x, ki, 0)

<-C1 при F(x) = C.

При k2 = 0 последняя система асимптотически устойчива по переменным ) и n(t), т. к. матрицы A0 и G - гурвицевы. Следовательно, переменные ), n(t) ограничены. Тогда из уравнений (17) и условий леммы [10] для некоторого к0 > 0 сигналы г'(t) и A(t) - ограниченные функции, т. е. х — sup|s'(t)| , ф — sup A(t), где х и ф - некоторые положительные числа.

Из вышесказанного и уравнений (4) и (14) следует ограниченность векторов х(0 и (t). Рассмотрим уравнение

ф(х, и, ^ = уТ (^х(^ + аТ (^х^--ф)) + 3(^и+ к(t)/(t) =

= уТ (t)х(t) + аТ (t)х(t - т(t)) + )(ио (t) + ик (t)) + к(t)/^) =

= уТ (t)х(t) + сТ (t)х(t - т(t)) + )(ик ^) - К0х (t)) + к(t)/(t).

Принимая во внимание (13), разрешим его относительно ик ^) :

Щ (t ) = -"^У+Г (УТ (t) х (t ) + °Т (t) х ^))-^) Кох (t) + к ^ ) / ^ )).

Ввиду ограниченности всех сигналов в правой части последнего равенства, что следует из доказанного выше и условия 4 Предположений, управление ик ^) также ограничено. Таким образом, все сигналы в замкнутой системе (8), (12), (13)-(16) ограничены. Однако из этого не следует асимптотическая устойчивость этой системы при ц2 ф 0 . Выберем функцию Ляпунова в виде

V(t) = юТ + (0Рп(0, (18)

положительно определенные симметрические матрицы N и Р определяются из уравнений

ЫА0 + A¡N = -21п; PG + GTP = -21п. Производная от функции V(t) вдоль траектории (17) при м1 = ц2 = м примет вид

'Л _

V(г) = -212--|п(t)|2 + 2юТ (t)N6М0-1 gT A(t) + 2пТ (t)Р6е'^). (19)

М0

Воспользуемся оценками

2ют (t)Nbtf;lgT А(t) < Ц |Nb|2 |ю(t)|2 + Ц0-1 |A(t)|2 < Ц |Nb|2 |ю(t)|2 + Ц0-1ф2;

__1 . .2 \ _|2 I |2 Ii i2 l _|2

2nT (t) Pb e'(t)<—|n(t)| \Pb\ +ц2 |e'(t)| <—|n(t)| \Pb\ + ц2х2. Ц0 Ц0

Подставив полученные оценки в (19), получаем

2 1 |_/ ч|2 I / \|2Л Y-ll,T, |2\ 1 |_/ \|2

V (t ) = -|o(t )|2--|n(t )|2-b(t ^(1 - ц0-11 Nb |2 )-^|n(t )|2+1 --| Pb |2 *+ц2 (ц0-3ф2+X2).

Ц0 ' Ц0

- Ц0 .

Очевидно, всегда найдется такое ц0, обеспечивающее положительность чисел в скобках последнего равенства. Тогда, в силу структуры функции Ляпунова (18), ее производную (19) можно оценить как

V(t) < aV(t) + ц2 (ц0-3Ф2 +х2), где a = min --7777; --, Xmin (•) - минимальное собственное число соответствующей

^n (N )' X min (P ) матрицы.

Решением последнего неравенства является

2

V (t )< e-^V (0) + (1 - )Ц- (ц0-3ф2 +x2)

из чего следует, что

НшV('(иГФ2 +Х2). а

Значит,

ш('(^У +х2)

•М/ (^Ф2 +Х2) при ' ^<х>. а

Пример. Рассмотрим объект управления вида (1), в котором А (') =

ах (') а2 (') а3 (')

" 0 0 0 " " 0 " " 0 "

С (' ) = 0 0 0 , В (' ) = 0 , я (' )= 0

_С (') с 2 (') С3 (')_ _ь (')_ d (')

L =[1; 0; ...; 0], х(0) = [1 1 1]' .

Класс неопределенности Е задан множеством Е = {а;, с, Ь, d :а1 е [-30; 30]; с е [-30; 30];

Ь е[1; 8]; d е[-10; 10]; / = 1, 2, 3}.

Измерению доступен выходной сигнал у('). Функция возмущений, действующая на объект управления, ограничена: /(')| < 10 .

Решим для номинальной системы (4)

хн ( ') =

"0 1 0 " "0"

0 0 1 хн ( 0 + 0

-1 -2 -2 1

.('); Ун (')=[1 0 0] Хн (')

задачу оптимального управления с критерием качества (2) при q = 3, г = 2.

Оптимальный закон управления (8) сформируем следующим образом:

и0 (') = -[0,5811 0,7144 0,3300]х ('),

где Х(') - оценка вектора х('), сформированная наблюдателем (7) с матрицей S = [10 50 29]' . Решением уравнения Лурье - Риккати А1нН + НАн - г 1НВнВ'Н = -М является матрица

Н =

4,5837 3,3681 1,1623 3,3681 3,4868 1,4288 1,1623 1,4288 0,6600

Зададим вспомогательный контур (12)

¿V (') = (')-ВнУ ('); ^ (0) = ?(0); е (') = (')

"10 1 0 " с матрицей А1 = Ан - SL = 50 0 1

28 -2 -2

Тогда уравнение относительно переменной ) = ') - (') примет вид

4 (' ) =

"-10 1 0 " "0"

-50 0 1 ) + 0

-30 -2 -2 1

ф(х, и, '); е(') = е(')-е('); 0)=[0 0 0]'

2

и

Для оценки производных сигнала s(t) воспользуемся наблюдателем (14):

15

Si (t ) = S, (t)+^^^Y (в( t )-g (t)),

i(t) = S3(t)+^^(6(0-ё(t)); é(t) = [Y o o]|(t), S 3 (t )=ol05t W t и (t ))•

Таким образом, закон управления (13) формируется следующим образом: v(t) = |з (t) +12|2 (t) + 72|i (t) + 150s(t);

uk (t) = -v(t)•

На рис. приведены графики переходных процессов по y(t) и ошибке слежения У (t)- Ун (t) при следующих значениях параметров в рассматриваемой системе: a1 (t) = 1 + sin0,2t, a2 (t) = 3 + sin0,5t, a3 (t) = 2 + sin0,7t, c1 (t) = 1 + 2sin0,6t, c2 (t) = -2 + sin0,2t, c3(t) = 3 + 2sin0,2t, b(t) = 2 + sint, d(t) = 2sin2t-1, f (t) = 1 + sint + sin(2t + n/6), x(t) = e- .

y(t) - yH(t)

y(t)

xio-4

Переходные процессы по у(7) и ошибке слежения у (7) - ун (t)

Значение интегрального критерия качества (2) для объекта управления (1) равно Jp(t) = 25,500, для номинального объекта (4) критерий принимает значение JН(t) = 25,501. Нетрудно видеть, что абсолютная погрешность отклонения этих значений очень мала и составляет 0,0001.

Заключение

Предложен принцип построения робастного субоптимального регулятора для параметрически и функционально неопределенного линейного нестационарного объекта по выходу. Цель управления состоит в синтезе системы регулирования, минимизирующей интегральный критерий качества (2) с некоторой погрешностью отклонения от номинального значения. Формирование управляющего воздействия базируется на идее, предложенной в [8], где управляющий сигнал предложено разложить на оптимальное управление, минимизирующее заданный функционал качества (2), и компенсирующее неопределенности объекта управления. Компенсация параметрических и функциональных неопределенностей основана на подходе [6].

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

1. Поляк Б. Т., Щербаков П. С. Робастная устойчивость и управление. М.: Наука, 2002. 303 с.

2. Неймарк Ю. И. Синтез и функциональные возможности квазиинвариантного управления // Автоматика и телемеханика. 2008. № 10. С. 48-56.

3. Никифоров В. О. Наблюдатели внешних детерминированных возмущений 1. Объекты с известными параметрами // Автоматика и телемеханика. 2004. № 10. С. 13-24.

4. Никифоров В. О. Наблюдатели внешних детерминированных возмущений 2. Объекты с неизвестными параметрами // Автоматика и телемеханика. 2004. № 11. С. 40-48.

5. Бобцов А. А. Алгоритмы робастного управления линейным объектом по выходу с компенсацией неизвестного детерминированного возмущения // Изв. РАН. Теория и системы управления. 2003. №2 2. С. 93-97.

6. Цыкунов А. М. Алгоритмы робастного управления с компенсацией ограниченных возмущений // Автоматика и телемеханика. 2007. № 7. С. 103-115.

7. Цыкунов А. М. Робастное управление с компенсацией ограниченных возмущений и помех // Изв. РАН. Теория и системы управления. 2014. № 3. С. 19-26.

8. Буков В. Н. Вложение систем. Аналитический подход к анализу и синтезу матричных систем. Калуга: Изд-во науч. лит. Н. Ф. Бочкаревой, 2006. 720 с.

9. Atassi A. N., Khalil H. K. Separation principle for the stabilization of class of nonlinear systems // IEEE Trans. Automat. Control. 1999. V. 44. N. 9. P. 1672-1687.

10. Брусин В. А. Об одном классе сингулярно-возмущенных адаптивных систем. I // Автоматика и телемеханика. 1995. № 4. C. 119-127.

Статья поступила в редакцию 15.05.2019

ИНФОРМАЦИЯ ОБ АВТОРАХ

Галяув Елена Романовна - Россия, 414056, Астрахань; Астраханский государственный технический университет; канд. техн. наук; доцент кафедры высшей и прикладной математики; [email protected].

Франгулова Елена Владимировна - Россия, 414056, Астрахань; Астраханский государственный технический университет; старший преподаватель кафедры высшей и прикладной математики; [email protected].

SUBOPTIMAL CONTROL OVER LINEAR NON-STATIONARY PLANT WITH DELAY

E. R. Galiauv, E. V. Frangulova

Astrakhan State Technical University, Astrakhan, Russian Federation

Abstract. The article focuses on the fact that the problem of suboptimal control over parametri-cally and functionally undefined linear non-stationary plant has been solved. It is supposed that only scalar input and output of plant are determinable. The purpose of control consists in submini-mization of the integral with an infinite top limit from square-law sub-integral function dependent on a scalar input of the plant and a control signal. The algorithm received is simple and does not require difficult analytical accounts of parameters of a control system. The serviceability of the received algorithms is illustrated on numerical examples.

Key words: robust control, suboptimal control, external disturbance, unspecified non-stationary linear plant, integrated criterion of quality, observer, Lyapunov function.

For citation: Galiauv E. R., Frangulova E. V. Suboptimal control over linear non-stationary plant with delay. Vestnik of Astrakhan State Technical University. Series: Management, Computer Science and Informatics. 2019;3:123-132. (In Russ.) DOI: 10.24143/2072-9502-2019-3-123-132.

REFERENCES

1. Polyak B. T., Shcherbakov P. S. Robastnaya ustojchivost' i upravlenie [Robust stability and control]. Moscow, Nauka Publ., 2002. 303 p.

2. Nejmark Yu. I. Sintez i funkcional'nye vozmozhnosti kvaziinvariantnogo upravleniya [Synthesis and functionality of quasi-invariant control]. Avtomatika i telemekhanika, 2008, no. 10, pp. 48-56.

3. Nikiforov V. O. Nablyudateli vneshnih determinirovannyh vozmushchenij 1. Ob"ekty s izvestnymi parametrami [Observers of external deterministic disturbances 1. Objects with known parameters]. Avtomatika i telemekhanika, 2004, no. 10, pp. 13-24.

4. Nikiforov V. O. Nablyudateli vneshnih determinirovannyh vozmushchenij 2. Ob"ekty s neizvestnymi parametrami [Observers of external deterministic disturbances 2. Objects with known parameters]. Avtomatika i telemekhanika, 2004, no. 11, pp. 40-48.

5. Bobcov A. A. Algoritmy robastnogo upravleniya linejnym ob"ektom po vyhodu s kompensaciej neizvest-nogo determinirovannogo vozmushcheniya [Robust control algorithms for output linear plant with compensation of unknown deterministic disturbance]. Izvestiya RAN. Teoriya i sistemy upravleniya, 2003, no. 2, pp. 93-97.

6. Cykunov A. M. Algoritmy robastnogo upravleniya s kompensaciej ogranichennyh vozmushchenij [Robust control algorithms with limited disturbance compensation]. Avtomatika i telemekhanika, 2007, no. 7, pp. 103-115.

7. Cykunov A. M. Robastnoe upravlenie s kompensaciej ogranichennyh vozmushchenij i pomekh [Robust control with limited disturbance and interference compensation]. Izvestiya RAN. Teoriya i sistemy upravleniya, 2014, no. 3, pp. 19-26.

8. Bukov V. N. Vlozhenie sistem. Analiticheskij podhod k analizu i sintezu matrichnyh sistem [Nesting systems. Analytical approach to the analysis and synthesis of matrix systems]. Kaluga, Izd-vo nauch. lit. N. F. Bochkarevoj, 2006. 720 p.

9. Atassi A. N., Khalil H. K. Separation principle for the stabilization of class of nonlinear systems. IEEE Trans. Automat. Control., 1999, vol. 44, no. 9, pp. 1672-1687.

10. Brusin V. A. Ob odnom klasse singulyarno-vozmushchennyh adaptivnyh sistem. I [On a class of singularly perturbed adaptive systems. I]. Avtomatika i telemekhanika, 1995, no. 4, pp. 119-127.

The article submitted to the editors 15.05.2019

INFORMA TION ABOUT THE A UTHORS

Galiauv Elena Romanovna - Russia, 414056, Astrakhan; Astrakhan State Technical University; Candidate of Technical Sciences; Assistant Professor of the Department of Higher and Applied Mathematics; [email protected].

Frangulova Elena Vladimirovna - Russia, 414056, Astrakhan; Astrakhan State Technical University; Senior Lecturer of the Department of Higher and Applied Mathematics; [email protected].