CC BY
133
УДК 517.9
DOI: 10.14529/mmph170306
СТАБИЛИЗАЦИЯ ОБРАТНОГО МАЯТНИКА НА ДВУХКОЛЕСНОМ ТРАНСПОРТНОМ СРЕДСТВЕ
В.И. Ряжских1, М.Е. Семенов2, А.Г. Рукавицын3, О.И. Канищева4, А.А. Демчук4, П.А. Мелешенко3
1 Воронежский государственный технический университет, г. Воронеж, Российская Федерация
2 Воронежский государственный архитектурно-строительный университет, г. Воронеж, Российская Федерация
3 Воронежский государственный университет, г. Воронеж, Российская Федерация
4 Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина», г. Воронеж, Российская Федерация
E-mail: [email protected]
Рассматривается механическая система, состоящая из двухколесной тележки, на оси которой располагается обратный маятник. Задача заключается в формировании такого управляющего воздействия, формируемого по принципу обратной связи, которое, с одной стороны, обеспечивало бы заданный закон движения механического средства, а с другой, стабилизировало бы неустойчивое положение маятника.
Ключевые слова: механическая система; двухколесное транспортное средство; обратный маятник; люфт; стабилизация; управление.
Введение
Возможность управления неустойчивыми техническими системами теоретически рассматривалась уже давно, однако практическая значимость такого управления отчетливо проявилась лишь в последнее время [1-5]. Оказалось, что неустойчивые объекты управления при подходящем управлении обладают рядом «полезных» качеств. Примерами таких объектов могут служить космический корабль на этапе взлета, термоядерный реактор и многие другие. В тоже время при выходе из строя автоматической системы управления неустойчивый объект может представлять собой существенную угрозу, опасность и для человека, и для окружающей среды. В качестве катастрофического примера результатов отключения автоматического управления можно привести аварию на Чернобыльской АЭС. По мере того, как системы управления становятся все более надежными, все более широкий круг технических неустойчивых в отсутствие управления объектов применяется на практике. Одним из самых простых примеров неустойчивых объектов является классический обратный маятник. С одной стороны, задача о его стабилизации сравнительно простая и наглядная, с другой, она может найти практическое применение при создании моделей двуногих существ, а также антропоморфных устройств (роботов, киберов и др.), перемещающихся на двух опорах. В последние годы появились работы, посвященные проблемам стабилизации обратного маятника, связанного с движущимся двухколесным транспортным средством [6-8]. Эти исследования имеют потенциальные перспективы применения во многих областях, таких как транспорт и разведка, в связи с компактной конструкцией, удобством эксплуатации, высокой маневренностью и низким расходом топлива таких устройств. Тем не менее, рассматриваемая задача еще далека от окончательного решения. Известно, что многие традиционные технические устройства имеют как устойчивые, так и не устойчивые состояния и режимы работы. Характерный пример - сегвей, изобретённый Дином Кейменом электрический самобалансирующийся самокат с двумя колёсами, расположенными по обе стороны от водителя. Два колеса скутера расположены соосно. Сегвей автоматически балансируется при изменении положения корпуса водителя; для этой цели используется система индикаторной стабилизации: сигналы с гироскопических и жидкостных датчиков наклона поступают на микропроцессоры, которые вырабатывают электрические сигналы, воздействующие на двигатели и управляющие их движениями. Каждое колесо сегвея приводится во вращение своим электродвигателем, реагирующим на изменения равновесия машины. При наклоне тела ездока вперёд сегвей начинает катиться вперёд, при увеличении же угла наклона тела ездока скорость сегвея увеличивается. При отклонении корпуса назад само-
кат замедляет движение, останавливается или катится задним ходом. Руление в первой модели происходит с помощью поворотной рукоятки, в новых моделях - качанием колонки влево-вправо. Задачи управления колебательными механическими системами имеют значительный теоретический интерес и большое практическое значение.
Известно, что в процессе функционирования механических систем вследствие старения и износа деталей неизбежно возникают люфты, упоры, поэтому для описания динамики таких систем необходимо принимать во внимание влияние гистерезисных эффектов. Математические модели таких нелинейностей в соответствии с классическими представлениями сводятся к операторам, которые рассматриваются как преобразователи на соответствующих функциональных пространствах. Динамика таких преобразователей описывается отношениями «вход-состояние» и «состояние-выход» [1].
Постановка задачи
В настоящей работе рассматривается механическая система, состоящая из двухколесной тележки, на оси которой располагается обратный маятник. Задача заключается в формировании такого управляющего воздействия, которое, с одной стороны, обеспечивало бы заданный закон движения механического средства, а с другой, стабилизировало бы неустойчивое положение маятника. При этом учитываются гистерезисные свойства в управляющем контуре изучаемой системы. Ниже графически представлены элементы, изучаемой механической системы - двухколесного транспортного средства с закрепленным на нем обратным маятником.
Рис. 1. Основные структурные элементы рассматриваемого механического устройства
тут / 1 / I feili / Fr I
" 1 " \ 1 \ 1 i R J
Hr ! / / / / /1 / / /
Рис. 2. Левое и правое колеса механического устройства с управляюшим моментом
Параметры и переменные, которые описывают рассматриваемую систему: j - угол поворота транспортного средства; D - расстояние между двумя колесами вдоль центра оси; R - радиус колес; Jj - момент инерции; Tw - разность крутящих моментов левого и правого колес; v -
продольная скорость транспортного средства; в - угол отклонения маятника от вертикального положения; m - масса перевернутого маятника; l - расстояние между центром тяжести тела и
осью колеса; Ти - сумма крутящих моментов левого и правого колес; х - перемещение транспортного средства по направлению продольной скорости; М - масса шасси; М* - масса колес; И - раствор люфта.
Динамика системы
Динамику системы описывают следующие уравнения:
j
D
2Rj
-T
-L Sil
n = —[J ml q2 sin в - m2l2g sin qcos в] + — Tn, W в á WR n
в = —[Mxmgl sin в - m2l2 в2 sin qcos в] - ml C0S в Tn,
W
WR
]2
(1)
(2) (3)
где Т* = Ть - ТЯ; Тп = Ть + ТЯ; Мх =М + т + 2(М* + ^*); 1в = т/2 + 1С ; 0.=Мх1в-т2/2 соъ2 в;
Я
П2 1
( = 1п+ — (М* + -*).
2 Я
Будем считать заданным желаемый закон изменения параметров системы:
<Р* = Рл С )Л = ^ С № = ^ О. (4)
Модель, описывающую динамику изменения параметров системы, можно представить в виде двух независимых подсистем. Первая подсистема состоит из одного уравнения - р -подсистемы,
определяющего угловые движения транспортного средства:
Р =
(5)
Уравнение (5) можно переписать в виде системы из двух уравнений:
¿1 = ¿2
D
2RJ
-Tw-jd'
(6)
где е1 =Р-Рй , е2 =(Р-(Ра .
Вторая подсистема, описывающая радиальные движения транспортного средства, а также колебания установленного на ней маятника, состоит из двух уравнений - {у,в} -подсистемы:
U =—[ Jqml в2 sin в- m2l2 g sin в cos в] + Jq Tu W в S J WR u
в =—[Mxmgl sin в - m2l2 в2 sin qcos в] - ml С°*в Tv W WR
Систему (7) удобно представить в виде системы уравнений первого порядка:
(7)
e3 = ¿4,
¿4 = ТГ" [ Jqml(qd + e6)2 sin(e5 +qd ) - m¿l2g sin(e5 + qd ) cos(e5 +qd )] + ТЩT v- Xd,
22
Jq
W0 R
TU XH
W
e5 = ч^
¿6 =~^[Mxmglsin(^5 +qd)-m2l2q + ¿б)2 sin(e + вй)cos(e5 + вй)]- ^^^ +в)
(8)
W
W0 R
Tu-qd,
где W0 = MxJq- П121 2cos2(qd + e5), e3 = X - Xd , ¿4 = v - vd , ¿5 =q-qd, ¿6 =q-qd •
Рассмотрим подсистему (6), управлять которой будем по принципу обратной связи. Для этого введем новую переменную и определим поверхность переключения в фазовом пространстве системы как ^ = 0 .
5 = в! + с1е1, (9)
где с - положительный параметр. Непосредственно из определения вытекает:
■Я = е+с1 е1 -срй + с1 е1. (10)
2—Ур
Для стабилизации вращательного движения определим управляющий момент следующим образом:
2—Зр .. .
Т№ Р - ^ в1 - -М§П(51) - к2 (11)
где , - положительно заданные параметры.
Аналогично будем строить управление второй подсистемой (8), управлять которой, будем также по принципу обратной связи. Для этого введем новую переменную и определим поверхность переключения в фазовом пространстве системы, как ■2 = 0 .
■2 = вз + С2вз, (12)
где с2 - положительный параметр, тогда
. .. . 1 . 2 2 2
■2 = е3 + с2 е3 = (в + в6) ^5 + вё ) - т 1 § ^5 + вс1 )С08(е5 + ва )] +
3 . 0. . (13)
+7^Т - + с2 ез
о —
Для стабилизации радиального движения определим управляющий момент:
тт'2/2 ^ к Т =-Кт/ (вй+еб)г^т(еь + вй)+яп^ + вй )е08(е5 + вй)--0- \сг ез - +^п^)+кА ^],(14)
где к3, к4 - положительно заданные параметры.
Для того, чтобы одновременно управлять обеими подсистемами системы, введем дополнительное управляющее воздействие:
= § Хапв--[ва + с3(в-вй) - к588п(^3) - кб 53], (15)
т/ 008 в
где § - ускорение свободного
падения; с3, к5, кб - положительные параметры; 53 - поверхность переключения, определяемая соотношением:
53 = е6 + с3е5 .
Сформулируем основные результаты работы, заключающиеся в принципиальной возможности стабилизации обеих подсистем, в сделанных предположениях относительно управляющих воздействий, в окрестности нулевого положения равновесия.
Теорема 1. Система (6) с управляющим воздействием (11) абсолютно асимптотически устойчива:
Нш || е11|® 0,
® (16)
Нш || е2 ||® 0. t®¥u 2
Доказательство: определим функцию Ляпунова как
1 2
V = ^as?, (17)
где a = Dj 2 RJр .
Очевидно, что функция V > 0, тогда
1 9
V =1М2, (21)
V =Ш1 Si = Si[rw-a(jd - Ci ei)]. (18)
Подставив (14) в V, получим
V = -(£ Sgn(S1) + k2(S1))S1. (19)
Очевидно, что V1 <pVx, где p - положительный параметр, что в силу теоремы Ляпунова означает асимптотическую устойчивость.
Теорема 2. Рассмотрим подсистему (8) с управляющим воздействием (14). В сделанных предположениях эта система абсолютно асимптотически устойчива, т. е. при любых начальных условиях выполняются соотношения:
lim ||e3 ||® 0,
t®¥ (20) lim 11 е41|® о.
Доказательство: определим функцию Ляпунова для системы (8) посредством соотношения
1 2
где b =Wo R!Je .
Очевидно, что функция V2 > 0, и
V2 = М S2 = S2[b(&d + ев)2 sin(e5 +0d) - b Sin(e5 + Od)cos(e5 + O) + Tv, (22)
где b1 = Rml, B2 = .
Je
Подставив (14) в (20), получим
V2 = -(£3 Sgn(S2 ) + k4 (S2 ))S2 . (23)
Очевидно, что V2 <pVx, где p - положительный параметр, что в силу теоремы Ляпунова означает асимптотическую устойчивость.
Динамика системы с гистерезисным управлением
Люфты и упоры возникают благодаря старению и износу механических параметров, поэтому представляет интерес изучение динамики рассматриваемой системы в ситуации, когда управляющее воздействие (Tv и Tw) воздействует на систему опосредованно, т. е. они являются входом на гистерезисный преобразователь - люфт, а его выход является управляющим воздействием на систему. Величина люфта определяет степень поворота элемента управления, которая не приводит к изменениям в управляемой системе. Наличие люфта в управляющем воздействии в общем случае осложняет задачу управления [9-13], так как возникают зоны нечувствительности по отношению к управляющему воздействию. Приведем краткое описание используемого в дальнейшем гистерезисного преобразователя - люфта, основанное на операторной трактовке. Выход преобразователя - люфта на монотонных входах описывается соотношением:
x(t0) при тех t, при которых x(t0) - h < u(t) < x(t0), x(t) = \u(t) при тех t, при которых u(t) > x(t0), (24)
u(t) + h при тех t, при которых u(t) < x(t0) - h,
которое иллюстрирует рис. 3.
С помощью полугруппового тождества действие оператора распространяется на все кусочно-монотонные входы:
Г [u(t1, h]x(t) = Г [ Г [u0, h]x(t1), h]x(t) (25)
и с помощью специальной предельной конструкции на все непрерывные. Так как выход этого оператора не является дифференцируемым, то в дальнейшем используется аппроксимация люфта моделью Боука-Вена [7]. Эта известная полуфизическая модель широко используется для феноменологического описания гистерезисных эффектов. Популярность модели Боука-Вена обу-
славливается ее способностью охватывать в аналитическом виде различные формы гистерезис-ных циклов. Формальное описание модели сводится к системе следующих уравнений:
Fbw (х, ^ = акх() + (1 -a)Dkz(t), = D"1( AX -р\х \\z \п-1 z -ухе | z |п ). (26)
Fbw(x,t) трактуется как выход гистерезисного преобразователя, а x(t) - как вход. Здесь п > 1,
D > 0 k > 0 и 0 <а< 1.
J лу Ги
Si/ U* //
•fyNf * / у
/
Su -l*h /
Ж /a "ic
г /
Рис. 3. Динамика входно-выходных соответствий люфта
Рассмотрим обобщение систем (6) и (8), в которых управляющее воздействие поступает на вход гистерезисного преобразователя, а выход является управляющим воздействием на систему:
¿1 = бо
D
2RJ
~Fbw ~Фс1
(27)
Fbw (х, t) = akx(t) + (1 - a)Dkz(t), z = D_1( Ax-b\x || z \n-1 z - gx | z\n).
e, = e&
1
¿4 = W-J mlQd + еб)2 sin(e5 + q) - m2l2g sin(e5 + ed) cos(e5 + 0d)] +
Wn
JQ
W0 R
-Ft
bw
e* = e.
1
A 9 9/
¿б = W [Mxmgl sin(e5 + Qd) - m2l 2Qd + ¿б)2 sin(e5 + 0d) cos(e5 + Qd)] -Fbw (x, t) = akx(t) + (1 - a)Dkz(t),
ml cos(e5 + Qd)
(28)
W0 R
F
bw
^ z = D_1( A x- b\x\\z\n-1 z-gx \ z\n).
Как и ранее в рассматриваемой системе, основным являлся вопрос о стабилизации, т. е. асимптотическом поведении ее фазовых переменных. Ниже приводятся графики при одних и тех же физических параметрах системы с люфтом и без люфта. Эта система исследовалась посредством численных экспериментов. Данная задача была решена в среде программирования Wolfram Mathematica.
Значения констант и начальные условия приведены ниже:
m = 3; M = 5; Mw = 1; D = 1,5; R = 0,25; l = 0,2; Jw = 1,5; Jc = 5;
Jv = 1,5; j(0) = 0;x(0) = 0; Q(0) = 0,2; y(0) = [ j(0) x(0) Q(0)f = [0,0,0,002f ; k1 = k3 = 0,5; k2 = k4 = 0,2; k5 = 3; k6 = 5; c1 = c2 = c3 = 0,5.
На рис. 8-11 приведены результаты численного моделирования исследуемой системы.
e2 =
d
Рис. 4. Зависимость е1, е2 от времени г = 50
Рис. 5. Зависимость е1 от е2, при г = 50
Рис. 6. Зависимость е1, е2 от времени г с наличием Рис. 7. Зависимость е1 от е2 с наличием люфта в обрат-люфта в обратной связи, при а = 0,5 ной связи, при а = 0,5, г = 50
Рис. 8. Зависимость е3, е4 от времени г
Рис. 10. Зависимость е3, е4, при г = 50
Рис. 11. Зависимость е5 от е6 при г = 50
Численные результаты показывают, что релейное управление, формируемое по принципу обратной связи, позволяет удерживать фазовые координаты в окрестности положения равновесия, при этом асимптотическая устойчивость не имеет места, что обусловлено наличием люфтов в контурах обратной связи исследуемых систем.
Заключение
Исследована стабилизации перевернутого маятника на двухколесной тележке в процессе ее движения по заданной траектории. Предложен алгоритм управления, обеспечивающий устойчивость системы и диссипативность при наличии гистерезисных звеньев в контуре обратной связи.
Работа выполнена при поддержке РФФИ грант № 16-08-003.
Литература
1. Красносельский, М.А. Системы с гистерезисом // М.А. Красносельский, А.В. Покровский.
- М.: Наука, 1983. - 272 с.
2. Стабилизация перевернутого маятника вертикальными осцилляциями с помощью гистере-зисного управления / М.Е. Семенов, Д.В. Шевлякова, О.И. Канищева, Д.В. Грачиков // Наукоемкие технологии. - 2012. - Т. 13, № 3. - С. 027-034.
3. Семенов, М.Е. Адаптивное управление неустойчивым объектом с гистерезисными свойствами / М.Е. Семенов, А.Г. Рукавицын, О.И. Канищева, А.Е. Пигарев // Вестник ВГУ. Серия: Системный анализ и информационные технологии. - 2014. - № 1. - C. 40-44.
4. Капица, П. Л. Маятник с вибрирующим подвесом / П.Л. Капица // УФН. - 1951. - Т. 44, № 3. - С. 7-20.
5. Мирошник, И.В. Теория автоматического управления. Нелинейные и оптимальные системы // И В. Мирошник. - СПб.: Питер, 2006. - 272 с.
6. Черноусько, Ф.Л. Управление колебаниями / Ф.Л. Черноусько, Л.Д. Акуленко, Б.Н. Соколов. - М.: Наука, 1980. - 383 с.
7. Ozturk, i. A novel method for producing pseudo random numbers from differential equation-based chaotic systems / i. Ozturk, R. Kill? // Nonlinear Dynamics, 2015, Vol. 80, Issue 3, pp. 11471157.
8. Решмин, С.А. Оптимальное по быстродействию управление перевернутым маятником в форме синтеза / С. А. Решмин, Ф.Л. Черноусько // Известия РАН. Теория и системы управления. -2006. - № 3. - С. 51-62.
9. Robust stabilization of fractional-order chaotic systems with linear controllers: LMI-based sufficient conditions / M.R. Faieghi, S. Kuntanapreeda, H. Delavari, D. Baleanu // Journal of Vibration and Control. - 2014. - Vol. 20. - Issue 7. - pp. 1042-1051.
10. Semenov, M.E. Inverted pendulum under hysteretic control: stability zones and periodic solutions / M.E. Semenov, D.V. Shevlyakova, P.A. Meleshenko // Nonlinear Dynamics. - 2014. - Vol. 75.
- Issue 1-2. - P. 247-256.
11. On the state feedback control of inverted pendulum with hysteretic nonlinearity / M.E. Semenov, D.V. Grachikov, A.G. Rukavitsyn, P.A. Meleshenko // MATEC Web of Conferences. - 2014. -Vol. 16. - P. 05009.
12. Stephenson, A. On an induced stability / A. Stephenson // Philosophical Magazine. - 1908. -Vol. 15. - Issue 86. - P. 233-236.
13. Неймарк, Ю.И. Динамические модели теории управления / Ю.И. Неймарк, Н.Я. Коган, В П. Савельев. - М.: Наука, 1985. - 400 с.
Поступила в редакцию 31 августа 2016 г
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 3, pp. 41-50
DOI: 10.14529/mmph170306 STABILIZATION OF INVERTED PENDULUM ON A TWO-WHEELED VEHICLE
V.I. Ryazhskih1, M.E. Semenov2, A.G. Rukavicin3, O.I. Kanishcheva4, A.A. Demchuk4, P.A. Meleshenko3
1 Voronezh State Technical University, Voronezh, Russian Federation
2 Voronezh State Architecture and Engineering University, Voronezh, Russian Federation
3 Voronezh State University, Voronezh, Russian Federation
4 Military Educational and Scientific Center of the Air Force "N.E. Zhukovsky and Y.A. Gagarin Air Force Academy", Voronezh, Russian Federation
E-mail: [email protected]
A mechanical system consisting of a two-wheeled vehicle with a fixed inverted pendulum is considered. The objective lies in forming such control action which, on the one hand, would provide movement of the system on the specified trajectory, and on the other hand would provide stabilization of the inverted pendulum in a neighborhood of unstable equilibrium position. Characteristic feature of the set objective is the fact that the control action is performed with the use of a hysteretic connection formalized by equations of the Bouc-Wen model. Equations of the researched system's model are divided into two independent subsystems of the second and the fourth order, which describe rotational and longitudinal modes of motion correspondingly. Control action on each of the subsystems, formed on the basis of the feedback principle, is constructed according to the relay law; at that in the phase space of each subsystem linear manifolds which determine switching surfaces are constructed. Theorems that ascertain the presence of errors in the assumptions made, which determine the discrepancy between the desired and simulated laws of the vehicle motion, as well as the deviation of the pendulum from the unstable equilibrium position asymptotically tend to zero, are formulated and proved in the article. In proving the theorems, the apparatus of Lyapunov functions was used. As it is shown in the article, the presence of the hysteresis component in the feedback loop makes manageability of the system more difficult, and as a result, in this case we can only talk about limitation of errors of the discrepancy and the dissipative motion of the pendulum. Results of computational simulation of dynamics of the system under study both under conditions of hysteresis connections and without them, which illustrate the proved theorems, are given in the article.
Keywords: mechanical system; two-wheeled vehicle; inverted pendulum; stabilization; control.
References
1. Krasnosel'skiy M.A., Pokrovskiy A.V. Sistemy s gisterezisom (Systems with hysteresis), Moscow, Nauka Publ., 1983, 272 p. (in Russ.).
2. Semenov M.E., Shevlyakova D.V., Kanishcheva O.I., Grachikov D.V. Inverted pendulum stabilization by vertical oscillations with hysteresis control. Naukoemkie tekhnologii (Science Intensive Technologies). 2012, Vol. 13, no. 3, pp. 027-034. (in Russ.)
3. Semenov M.E., Rukavitsin A.G., Kanishcheva O.I., Pigarev A.E. Adaptive control of unstable objects with hysteresis properties. Vestnik VGU. Seriya: Sistemnyy analiz i informatsionnye tekhnologii, 2014, no. 1, pp. 40-44. (in Russ.).
4. Kapitsa P.L. Mayatnik s vibriruyushchim podvesom (The pendulum with vibrating suspension). UFN(Phys. Usp.), 1951, Vol. 44, no. 3, pp. 7-20. (in Russ.).
5. Miroshnik I.V. Teoriya avtomaticheskogo upravleniya. Nelineynye i optimal'nye sistemy (The theory of automatic control. Non-linear and optimal systems). St. Petersburg, Piter Publ., 2006, 272 p. (in Russ.).
6. Chernous'ko F.L., Akulenko L.D., Sokolov B N. Upravlenie kolebaniyami (Vibration control). Moscow, Nauka Publ., 1980, 383 p. (in Russ.).
7. Ozturk i., Kill? R. A novel method for producing pseudo random numbers from differential equation-based chaotic systems. Nonlinear Dynamics, 2015, Vol. 80, Issue 3, pp. 1147-1157. DOI: 10.1007/s11071-015-1932-5
8. Reshmin S.A., Chernous'ko F.L. Time-optimal control of an inverted pendulum in the feedback form. Journal of Computer and Systems Sciences International, 2006, Vol. 45, no. 3, pp. 383-394. DOI: 10.1134/S1064230706030063
9. Faieghi M.R., Kuntanapreeda S., Delavari H., Baleanu D. Robust stabilization of fractional-order chaotic systems with linear controllers: LMI-based sufficient conditions. Journal of Vibration and Control, 2014, Vol. 20, Issue 7, pp. 1042-1051. DOI: 10.1177/1077546312475151
10. Semenov M.E., Shevlyakova D.V., Meleshenko P.A. Inverted pendulum under hysteretic control: stability zones and periodic solutions. Nonlinear Dynamics, 2014, Vol. 75, Issue 1-2, pp. 247-256. DOI: 10.1007/s11071-013-1062-x
11. Semenov M.E., Grachikov D.V., Rukavitsyn A.G., Meleshenko P.A. On the state feedback control of inverted pendulum with hysteretic nonlinearity. MATEC Web of Conferences, 2014, Vol. 16, p. 05009. DOI: 10.1051/matecconf/20141605009
12. Stephenson A. On an induced stability. Philosophical Magazine, 1908, Vol. 15, Issue 86, pp. 233-236.
13. Neumark Y., Kogan N., Saveliev V.P. Dinamicheskie modeli teorii upravleniya (Dynamic models of control theory). Moscow, Nauka Publ., 1985, 400 p. (in Russ.).
Received August 31, 2016
