CC BY
14ISSN 1814-1196 Научный вестник НГТУ том 76, № 3, 2019, с. 33-56
http://journals.nstu.ru/vestnik Science Bulletin of the NSTU Vol. 7б, No. 3, 2019, pp. 33-5б
ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ
INFORMATICS, COMPPUTER ENGINEERING AND MANAGEMENT
УДК 681.5.865.8 Б01: 10.17212/1814-1196-2019-3-33-56
Аналитический синтез исполнительных систем микророботов в условиях неопределенности
А Р. ГАЙДУК1'", С.Г. КАПУСТЯН2'4
1 347922, РФ, г. Таганрог, пер. Некрасовский, 44, Институт радиотехнических систем и управления, Южный федеральный университет
344006, РФ, г. Ростов-на-Дону, пр. Чехова, 41, Южный научный центр Российской академии наук
а [email protected] ь [email protected]
В последние годы для решения различных задач всё более широко применяются роботы различных типов, в том числе и микророботы. В связи с этим задача разработки систем управления исполнительного уровня микророботов становится весьма актуальной. Особенно остро эта проблема встает в тех случаях, когда порядок и параметры подвижных элементов микророботов - объектов управления - являются априори неопределенными и могут изменяться в процессе эксплуатации. Для управления такими микророботами на исполнительном уровне наиболее целесообразно применять адаптивные системы с идентификацией. Известные методы решения задачи адаптивного управления чаще всего ориентированы на случай, когда порядок объекта известен и не изменяется.
Ниже предлагается новый метод решения задачи синтеза адаптивного управления микророботами в условиях неопределенности, в основе которого лежит марковский метод идентификации и метод аналитического синтеза систем с управлением по выходу и воздействиям. Отличительной особенностью марковского метода идентификации является применение марковских параметров не непрерывного объекта, а марковских параметров виртуального дискретно-подобного объекта, являющегося совокупностью дискретных объектов различных порядков. При этом установлены условия разрешимости задачи идентификации микророботов этим методом. Метод аналитического синтеза систем с управлением по выходу и воздействиям обеспечивает синтез управлений минимальной сложности за счет согласования полюсов замкнутой системы с нулями и полюсами полного объекта; также обеспечиваются порядки астатизма и прямые показатели качества не хуже заданных. Эти свойства достигаются за счет использования стандартных нормированных передаточных функций. Эффективность предложенного похода подтверждается результатами компьютерного моделирования.
Статья получена 15 февраля 2019 г.
Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов № 18-58-00051, № 17-29-07054, № 16-29-04194 и № 18-07-00637, а также в рамках реализации Программы фундаментальных исследований РАН по приоритетным направлениям,
определяемым Президиумом РАН, № 7 «Новые разработки в перспективных направлениях
энергетики, механики и робототехники», № гр. проекта ЛЛЛЛ-Л18-118020190041-1).
Разработанный подход может применяться для создания в условиях неопределенности как адаптивных исполнительных систем микророботов, так и адаптивных систем управления объектами промышленного, сельскохозяйственного производства, а также для решения специальных задач, в том числе задач управления в социальной сфере.
Ключевые слова: микроробот, объект, идентификация, марковские параметры, система управления, синтез, управление по выходу и воздействиям
ВВЕДЕНИЕ
Роботы, и в частности микророботы, применяются для решения всё более широкого круга задач. При этом происходит непрерывное повышение требований к качеству их систем управления. Решение задачи синтеза этих систем требует наличия математических моделей управляемых элементов роботов. Однако последние часто являются неопределенными из-за изменений свойств как роботов, так и микророботов в процессе эксплуатации, что обусловливает необходимость применения адаптивных систем управления [1-7]. Современные адаптивные системы с моделью или системы прямого адаптивного управления ориентированы на случай, когда изменяются параметры, а порядок неопределенного объекта не изменяется с течением времени [8, 9]. В этом плане более эффективными являются адаптивные системы управления с идентификацией, так как они могут применяться и при изменении порядка микророботов в некоторые моменты времени. Такие, по существу, самоорганизующиеся системы создаются с применением цифровых средств автоматизации и нуждаются в алгоритмическом и программном обеспечении специального типа. Это обеспечение может быть создано с применением только аналитических методов идентификации моделей микророботов и синтеза систем управления ими [7, 10-13].
К настоящему времени разработано достаточно много алгоритмов как идентификации объектов, так и синтеза систем управления [5-16], однако большая часть из них ориентирована, во-первых, на одномерные системы, а во-вторых, на применение «вручную». Кроме того, эти методы часто не учитывают дискретность цифровых методов обработки информации.
Одним из наиболее ранних аналитических методов идентификации, по-видимому, является метод Б. Хо [14]. Этот метод разработан в середине XX века и использует свойства марковских параметров динамических систем. Однако он не получил широкого применения, так как ориентирован на использование импульсных переходных функций, которые являются реакцией объектов на дельта-функции, формирование которых практически невозможно.
Известен также аналитический метод идентификации, который базируется на решении уравнения Винера-Хопфа [15]. Однако в этом методе предполагается применение статистической обработки данных, поэтому процесс идентификации этим методом имеет большую длительность. Эта особенность в значительной степени осложняет его применение в системах микророботов.
Существуют методы идентификации неопределенных объектов, представляющие собой процедуру совместного оценивания переменных состояния и параметров, в основе которых лежит алгоритм обобщенного филь-
тра Калмана-Бьюси [13, 14]. Но и в этом случае, как правило, порядок управляемых объектов предполагается известным. Идентификацию и порядка, и параметров объекта можно осуществить, используя наблюдатель, применяемый в самоорганизующемся оптимальном регуляторе с экстраполяцией (СОРЭ) А.А. Красовского [10]. Однако этот наблюдатель строится на основе цепочек интеграторов различных порядков без обратных связей, что снижает эффективность этого наблюдателя, поскольку реальные объекты, в том числе и микророботы, всегда содержат внутренние обратные связи.
В ряде работ применяется метод активной идентификации на основе специальных номограмм, которые строятся заранее экспериментальным путем [16]. При этом предполагается, что порядки и объекта, и системы не изменяются; заранее известны интервалы изменения параметров системы. Эти факторы существенно сужают область применения данного метода активной идентификации в исполнительных системах микророботов.
При адаптивном управлении с идентификацией математическая модель микроробота, необходимая для создания алгоритма управления, формируется в реальном времени путем подачи на него специальных пробных воздействий и фиксации его реакции [8, 11, 12, 15, 16]. Это приводит к необходимости оперативной обработки большого объема данных и решения значительного числа задач. За короткое время эти задачи могут быть решены только вычислительными средствами с высокой производительностью, что требует разработки специальных методов идентификации и синтеза, ориентированных на цифровую реализацию [10, 17].
В настоящей статье предлагается аналитический метод синтеза адаптивных исполнительных систем управления с идентификацией для микророботов в условиях неопределенности. Идентификацию предлагается проводить с применением марковских параметров, являющихся инвариантами динамических систем [17, 18]. Как известно, для экспериментального определения марковских параметров непрерывных объектов необходимо многократное дифференцирование выходных переменных, что практически невозможно. С целью преодоления этой сложности предлагается определять марковские параметры не непрерывных микророботов, а некоторых виртуальных дискретных объектов. Этим обеспечиваются значительные преимущества предлагаемого подхода, так как марковские параметры указанных дискретных объектов могут быть легко найдены по значениям реакций микророботов на постоянные воздействия без дифференцирования. Эта особенность данного подхода позволяет идентифицировать неопределенные микророботы в автоматическом режиме и в реальном времени. При этом синтез исполнительных систем проводится на основе получаемых моделей с применением принципа управления по выходу и воздействиям [19, 20]. Предлагаемый метод ориентирован на цифровую обработку экспериментальных данных и позволяет распараллелить как идентификацию, так и синтез и значительно уменьшить длительность переходных процессов адаптивных систем управления за счет использовании многопроцессорных вычислительных средств [17].
1. ПОСТАНОВКА ЗАДАЧИ
Предположим, некоторый управляемый элемент микроробота (УЭМР) с несколькими каналами вход-выход описывается уравнениями
x = Ax + Bu , у = cTx + dTu , (1)
где x e Rn - вектор состояния размерности n; u e Rm - вектор управлений; y e R1 - выходная переменная этого элемента; A, B, c, d - числовые матрицы и векторы соответствующих размерностей [3, 17, 19]. Порядок УЭМР n и параметры его модели (1) заранее неизвестны и могут претерпевать скачкообразные изменения, оставаясь затем неизменными в течение достаточно длительного интервала времени. Причем во всех случаях УЭМР (1) является полным, а его порядок n < nmax -1, где nmax - заранее известное значение [7, 9]. Выходная управляемая переменная y = y(t) доступна измерению. Пробные воздействия имеют вид: u (t) = u;01(t), Uj (t) = 0, j Ф i, i, j = [1, m], значения u^ = const малы и допустимы.
Задача состоит в разработке метода цифровой идентификации текущей модели УЭМР (1) и синтеза адаптивного цифрового регулятора на основе этой модели. Идентификация должна осуществляться по результатам измерения реакции УЭМР на пробные воздействия. Адаптивный цифровой регулятор (АдЦР) должен быть физически реализуемым. Он должен обеспечивать заданные порядки астатизма к внешним воздействиям и прямые показатели качества замкнутой системы управления не хуже заданных [20, 21].
Переходя к решению поставленной задачи, введем вспомогательные определения. Пусть Ти - период дискретизации управляемой переменной y = y(t) при идентификации, а yik = yi (кТи) - дискретные значения отклонения переменной y(t) объекта (1) на пробные воздействия ui (t) = u;o1(t) при всех i = [1, m] и x(0) = xo = 0.
2. ДИСКРЕТНО-ПОДОБНЫЙ ОБЪЕКТ
Рассмотрим сначала предлагаемый метод идентификации. Имея в виду цифровую обработку данных, введем виртуальный дискретно-подобный объект, соответствующий УЭМР (1), следующим образом.
Допустим, имеется совокупность систем уравнений различных порядков:
+1 = Axk + B uk, yk = cT4 + dTuk . (2)
Здесь xk = [ x1k x2 k --- xvk]T и uk = [u1k u 2 k --- umk]T - вект°ры отклонений переменных состояния и дискретных управлений, а yk - вектор отклонений yik выходной переменной y(t), i = [1, m], вызванных m-вектором uk,
к = 0,1, 2, •.., от соответствующих значений в установившемся режиме; А, В, с, й - числовые матрицы и векторы соответствующих размерностей виртуальных систем (2), порядок которых V = [1, птах].
Пусть при Х0 = 0 и всех ^ = кТи выполняются условия:
% = ui0, У(к = У(к, 1 =[1, к = 0,1,2, ••• (3)
Определение 1. Объект, описываемый системой уравнений (2), называется «дискретно-подобным объектом», соответствующим непрерывному УЭМР (1), если выполняются условия (3).
Подчеркнем, что дискретно-подобный объект (ДПО) (2), соответствующий приведенному определению, как и СОРЭ А.А. Красовского [9], включает птах систем уравнений (2) различных порядков v = [1, птах] [16, 17, 20].
Виртуальный объект (2) возникает в результате рассмотрения непрерывного пробного воздействия и^ (^) = ию ) (как последовательности прямоугольных импульсов) длительностью Ти и амплитудой и^ . Поэтому передаточные функции (ПФ) Жуи. (2, Ти) и Жуи, (я) каждого канала и ^ у ДПО и и^ ^ у УЭМР связаны в силу условий (3) ZT-преобразованиями:
Жуи (2,Ти,щ) = гт {жущ (я)} и жущ (я) = 2ТХ (2,ТИ,щ)}, щ = щ, (4)
где п, щ - степени знаменателей ПФ Жуи. (2, Ти, щ) и Жуи, (я) соответственно; i = [1, т].
Отметим, что функции «с2^> и <^2с» с расширением «^оЬ'» пакета МЛТЬЛБ осуществляют оба преобразования (4) [22]. Однако при наличии «левых» полюсов у элемента (1) и при некоторых значениях Ти информация о некоторых из этих полюсов теряется при формировании ДПО (2), и поэтому возможность идентификации такого элемента (1) на основе ДПО (2) исключается. В связи с этим вводится следующее определение. Пусть ^п - свободный коэффициент знаменателя ПФ Жуи. (2, Ти, щ), а А0 - погрешность вычислений на используемой ЭЦВМ.
Определение 2. Если период дискретизации Ти = Ти°, и ПФ ®уй. ^, Ти1, ^ такова, что
ф » А0, (5)
то ДПО (2) при Ти = Ти является корректным. В противном случае, т. е. при ^ 1п ^ А0, этот ДПО является некорректным [16, 17, 20].
3. МАРКОВСКИЕ ПАРАМЕТРЫ
В случае многомерных дискретных объектов типа (2) вводятся либо
1 х т -матрицы марковских параметров П 1 =
д1 д 2 • • • д т
формулами:
П 0 = СТ , П ?= сТА 1-1В , 1 = 1,2
Т л 1-1;
(6)
либо марковские параметры Д 1 каждого канала й; ^ у аналогичными выражениями:
Д 0 = С , д1 = сТА1-1Ъ*, 1 = 1,2, •.
(7)
где - элементы вектора-строки СТ ; Ъ; - столбцы п х т -матрицы В; I = [1, т] [17, 18].
Система уравнений (2) является дискретной, поэтому вектор у£ ее выходных переменных определяется выражениями:
Уо = СТХо +П0йо, ук = сТАкХо +Х£;=0Пк-]й]-, к = 1,2,... (8)
Из выражений (8) с учетом Хо = 0 и условий (3) выводим соотношения, которые определяют 1х т -матрицы П1 марковских параметров:
П0 =
У10 У20 ут 0 й10 й20 йт 0
П? =
У11 У21 Ут 1 й11 й21 йт
<Н - ■
-£П; , 1 = 1,2, • . (9) } =0
Марковские параметры Д1 ДПО (2) определяются аналогичными выражениями:
д 0=уг-0 / й0, д1 = уг? й-1 - Е/=0д/, 1=12,
(10)
где / = [1, т].
Соотношения (9), (10) позволяют найти любое число Nд 1х т -матриц
П1 или марковских параметров Д1 каждого канала й; ^ у виртуальной системы (2), используя экспериментальные дискретные значения отклонений у;к УЭМР и значения й^ , ; = [1, т], к = 0,1, 2,....
Подчеркнем, что в соответствии с определением ДПО значения входных и выходных переменных всех виртуальных систем (2) одни и те же, независимо от их порядка V = [1, птах ]. Поэтому все векторы П1 и все марковские
параметры ц С , i = [1, m], с = 0,1, 2, •.. тоже одни и те же для всех nmax виртуальных систем (2), соответствующих непрерывному УЭМР (1).
Передаточная 1х m -матрица системы (2) v-го порядка [18, 19] определяется выражением:
Жрй(г, V) = А-1(г, v){cT аС|(zEv - А)В + йТА(г, V)}, (11)
где ас|(zEV - А) - присоединенная матрица, а полином
А(г, V) = ае1 (гЕу- А ) = Л г , + ^1 +... + Л г + П0, ^ = 1. (12)
Как известно, присоединенная -матрица определяется выражением [23, с. 88]
аС|(гЕ - А) =
= ЛЕг^ + (^А + л^-1Е)zV-2 +... + (^-1 + ^—1Ау-2 +... + ЛЕ) . (13) Подставив выражение (13) в (11) и учтя равенства (6), получим
Жуи(г, V) = А-1(г, 0 гу + +п,-1Й0) г+
+ (л№2 +лV-lП1 +^-2П0) гV-2 +
+...+лV-lПV-1 +^-2^-2+...+ЛГП1 +л vп0)}. (14)
Выражение (14) определяет передаточную 1 х т -матрицу многомерной системы (2) v-го порядка. Из этого выражения следует, что V + 1 матриц Пс (6) или марковских параметров ц с , с = 0, V системы (2) и V коэффициентов л Р ,
* г
р = 0, V-1 характеристического полинома (12) полностью определяют все передаточные 1 х т -матрицы Жуи (г, Ти, V), так как л 1.
Но передаточная 1 х т -матрица Жуи (г, Ти, V) (11) может быть разложена в ряд:
Жуи (^ Ти, v) =
=А-1 (г, V) {а; гу + а- гУ-1 + •+а? г+а0у+г-1 + г-2 + •..}, (15)
где О р = [у^р ] - 1х т -матрицы постоянных коэффициентов. Сравнение выражений (14) и (15) приводит к равенствам:
^у-р _ ¿^ Пу+^-р11 ?=0
р = 0,1, 2,..., V-1, V, у +1,
(16)
Вг^р ~ ^0 ~ ^ 1 ~ V Лу-1 ~ ^гтР
частности, матрица О0 = л 0П + Л 1П + ... + ЛучП + ЛуП .
На основе выражений (14) и (15) заключаем, что элементы у V.- 1 хт -мат-
V —
риц О- -
У
■, ]
(16) при ] = [0, V] и ■ = [1, т] фактически являются коэффи-
циентами числителей скалярных передаточных функций Жуй. (г, Ти, V) каналов й^ ^ у виртуальных систем (2) порядка V. Эти функции имеют вид
Жущ (г Ти, V) =
У г,0 +У 1,1г + . + У г
- V - V - V V-1 V
Л 1,0+Л г,1г + . + Л г>-1г + г
(17)
а сами функции Жуй. (г, Ти, V) являются элементами передаточной 1 хт -матрицы Жуй (г, Ти, V) (11). Отметим, что степени полиномов в выражении (17)
могут быть меньше V в частных случаях.
В тех случаях, когда передаточная 1 х т -матрица Жуй (г, Ти, V) ф 0, то в
равенстве (15), как минимум, одна из 1х т -матриц О^ ф 0, ] е [0, V], и все
1хт -матрицы О^ = 0 при ] = -1, - 2,.... Поэтому из выражений (14), (15) и (17) следует система алгебраических уравнений
ПР+1 ПР+2 . .. ПР+V " "Л V" " П Р+v+1 "
ПР+2 ПР+3 . . ПР+v+1 = - ПР+v+2
ПР+v+1 . . ПР+2v-1 Л V ПР+2v
Р е [0, Ртах], (18)
а также система равенств
о;=П 0, GvV-]■ =П]+] , ] = [1, v-l]
(19)
определяющих 1 х т -матрицы О^- .
Порядок идентифицируемого непрерывного УЭМР (1) может быть определен по значению определителя матрицы из системы (18) при различных значениях V и р. Однако размер этой матрицы может быть довольно большим, поэтому вычислить его практически сложно, особенно учитывая
ограниченные возможности микророботов. В связи с этим более целесообразной представляется оценка непосредственно скалярных передаточных функций (г, Ти, V) каналов й- ^ у . При этом порядок п УЭМР (1) оценивается по оценкам щ степеней знаменателей ПФ Жуи. (г, Ти, V) ДПО (2). Эти оценки вычисляются по формулам:
П = тах
{45 V* 0},
5 ' = Ыг
Р,4
(20)
где
"й Р+1 йР+2 . . йг
мР,^ йР+2 йР+3 . . йР+V+1
_й Р+V йР^1 _ . йР^-1
4 [1, ттах],
Ре[0, Р
maxJ • (21)
Учитывая ограниченные возможности вычислительных средств микророботов оценки коэффициентов знаменателя передаточной функции Жуй. (г, Ти, V) (17) при 4 = щ , ДПО целесообразнее вычислять не по (18), а по
вытекающим из этого сложного матричного выражения формулам:
Л п=-( м Р п )-1 X Р
Рп
(22)
где
Л г' =
~ пг ~ пг
Л ',0 Л Л
Л
г,П'- -1
X Р,пг = г
~ Р+п +1 ~Р+пг + 2 й I й '
й г
Р+п +п
(23)
По той же причине оценки коэффициентов у4,р числителя передаточной функции ДПО (17) при V = п вычисляются по формулам:
~пг ~ пг -1
У г] = й I -
Уп'-1-] ]
р=0, п -1, у пп =й.
(24)
В выражениях (18)-(24) ' = [1, т], а в может принимать любое значение из интервала [0, ртах]. После определения численных значений коэффициентов Лпр и уп] они подставляются в выражение (17) при V = п , а полученная
передаточная функция Жуи . (г, Ти, п?-) подвергается обратному 2г-преобра-
зованию. В результате получается искомая передаточная функция Жуи. (5) =
= В. (5)/ А. (^) канала й. ^ у УЭМР. Затем описанные операции повторяются для остальных его каналов.
Матрицы (21) можно представить с учетом (6) следующим образом:
~Т
мР ,v =
сТА
АЪ1
АУ-1Ъ1
■ = [1, т] .
(25)
Из выражений (20) и (25) следует, что если Р = 0,1, 2,... и detМР^ при
некотором V* не равен нулю, а при V = v* +1 равен нулю, то полная часть канала й. ^ у ДПО имеет порядок п. = v* . Причем определитель матрицы МР,У пропорционален det А, если РФ 0, поэтому равенство (22) применимо при этом условии, если только det А Ф 0. Это неравенство соответствует, очевидно, условию корректности ДПО (5) по отношению к идентифицируемому УЭМР (1), так как оно выполняется только при соответствующем
значении периода дискретизации Ти° его выходных переменных.
Отдельные каналы даже полных многомерных УЭМР типа (1) или ДПО типа (2) могут быть неполными [7, 19, 24]. Поэтому в общем случае оценку п порядка идентифицируемого полного УЭМР целесообразно определять по степени его характеристического полинома, пользуясь следующими соотношениями:
А(5) = НОК {А .(5)}, п = deg Л($):
■=[1, т]
(26)
где А {(5) - знаменатели полученных по выражению (4) оценок Жущ (5) = В. (5)/ А ¿(5). Тогда оценку искомой передаточной 1 х т -матрицы Жуй (5) неопределенного УЭМР (1) можно оценить, пользуясь выражением
Жуй (5) = А-1( 5)[ В1( 5) ... Вт ( 5)(21т (5)] :
уй
(27)
где ((уй■■ (5) = А(5) / А ■ (5), ■ =[1, т].
Порядок применения к неопределенному УЭМР марковского алгоритма идентификации, вытекающего из приведенных выражений, будет показан ниже на численном примере. Этот алгоритм позволяет идентифицировать как одноканальные, так и многоканальные УЭМР и другие объекты, причем идентификацию отдельных каналов многомерных объектов можно осуществлять параллельно. Только выражения (26) и (27) реализуются совместно. Распараллеливание процесса идентификации позволяет значительно сократить его длительность за счет применения вычислительных средств, ориентированных на параллельную обработку информации.
Адаптивную систему управления многомерным неопределенным УЭМР с числом управлений, равным числу управляемых переменных, целесообразно синтезировать с применением декомпозирующего управления. Это управление позволяет создавать адаптивные регуляторы сепаратных каналов й{ ^ у., ■ = [1, т] как для одномерных, так и для многомерных объектов, но с учетом требуемых влияний на у. других управлений йу, у Ф ■ [24]. В связи с этим
метод аналитического синтеза систем с управлением по выходу и воздействиям (АССУВВ) [5,19, 20] для наглядности излагается здесь применительно к двухканальному УЭМР с одной управляемой переменной.
4. ЗАДАЧА СИНТЕЗА АДАПТИВНОГО УПРАВЛЕНИЯ
Предположим, в результате идентификации некоторого УЭМР изложенным выше методом получена передаточная 1х 2-матрица Жу^ (5) типа (27).
Эта матрица была подвергнута прямому Хт -преобразованию (4) при Т = Ту >> Ти , в результате чего получено уравнение
А( г) у (7) = В (7 )й( 7) + Е (г) / (7)
(28)
где у(г), й(г) и /(г) - г-изображения управляемой переменной, управления и внешнего неизмеряемого возмущения микроробота; А(г), В(г), Е(г) - некоторые полиномы с известными значениями коэффициентов, причем полиномы А( г) и В (г) не имеют общих множителей.
Метод АССУВВ ориентирован на применение принципа «управления по выходу и воздействиям», поэтому уравнение адаптивного цифрового регулятора (АдЦР) для УЭМР (28) имеет вид
Я( г )й (г) = ( г) g (г) - Ь( г) у( г), (29)
где g(z) - г-изображение задающего воздействия; Я(г), ((г), Ь(г) - искомые
полиномы. Как видно, выход и задающее воздействие УЭМР подаются на АдЦР (29) в соответствии с принятым принципом управления, а возмущение / не подается, так как не измеряется. Величина
ццр = шт{г - д, г - l}, (30)
где г = deg Я( г), д = deg ((г), I = deg Ь (5), называется относительным порядком АдЦР [20]. Полиномы Я(г), ((г), Ь(г) должны удовлетворять условию физической реализуемости
Ццр , (31)
где ц^ - допустимое значение ццр . Обычно ц^ = 1 или ц^ = тзап / Ту , где т зап - суммарное запаздывание по времени в измерительной системе и в цифровом регуляторе.
Замкнутая система (28), (29) должна иметь заданные порядки астатиз-ма ^ и по задающему воздействию g и по возмущению /; длительность
переходного процесса не более ^п , с; перерегулирование не более с , %, и степень устойчивости не хуже Р* < 1 [20, 21].
5. РЕШЕНИЕ ЗАДАЧИ СИНТЕЗА
Условие на степень устойчивости синтезируемой системы
Рсис - max| j>ß*, (32)
где zj - корни ее характеристического полинома, определяет на комплексной плоскости область Q - допустимого расположения корней этого полинома. Полином A(z) из уравнения УЭМР (28) факторизуется относительно
границы области Q : A(z) = A q(z)A q(z), где полином A q(z) включает все корни полинома A(z), удовлетворяющие условию (32), а A q(z) = = A(z)/ A q(z), причем det A(z) = nj, det Aq (z) = «q, а det A q(z) = «q .
Пусть в уравнении (28) B (z) является числом или полиномом степени &, все корни которого удовлетворяют условию (32), т. е. B(z) = ß& или B(z) =
= ßsn&=i(z -zB j . Другими словами, B(z) = Bq (z). В этом случае характеристический полином синтезируемой системы берется в виде D(z) = A q(z)B0(z)D(z), где B0(z) = ß--1B(z), а D(z) - некоторый полином, причем D( z) = Dq (z) [20].
В общем случае полиномы A(z) и F(z) из (28) могут иметь вид
A(z) = (z - 1)njA1(z), F(z) = (z -1)«FF1(z), где A1(1) * 0, F1(1) * 0. В этом случае дискретная система (28), (29) имеет требуемые порядки астатизма, если
v = max{Cg-nj; С/ -«f; 0} (33)
сумматоров имеются в алгоритме АдЦР. С целью учета условия (33) и согласования полюсов системы с нулями и полюсами УЭМР (28) полиномы R( z) и L(z) выбираются в виде
R( z) = B0 (z)(z -1)v R( z), L (z) = A q (z) L (z), (34)
где R(z), L(z)- неопределенные полиномы степеней r = r v и ~l = r-^цр -nQ .
С учетом равенств (29) и (34) характеристический полином &(г) адаптивной системы управления УЭМР (28) определяется выражением
&(2) = А д(7)Во (г)&(г) = А(г)Б0 (г)(г -1)*Я(г) + РзБ0 (г)А д (г)X(г). (35) Из равенства (35) вытекает расчетное полиномиальное уравнение
&( г) = А д (г)( г -1) * Я (г) + раХ( г), (36)
минимальное решение которого с учетом условия (31) определяется следующей системой алгебраических уравнений:
Рз 0 0 Ло 0
0 р3 0 Л1 Л0
0 0 '•. ! л1
! 0 ' Ля :
0 : ' 0 Ля
0 0 ... 0 0
10
Р0
«0 «1
]"В -1
пи
(37)
Матрица системы (37) имеет I +1 столбцов, составленных из коэффициента Из, и г +1 столбцов, составленных из коэффициентов Лр полинома
Н (г) = (г -1)* А д (г) = Хр=0 Лр г Р .
Согласно уравнению (36) deg&(г) = п& = V + яд + г, поэтому система (37) включает = п& +1 уравнений и = Г +1 +1 + 1 неизвестных коэффициентов, т. е.
Nу = п& +1 = яд + * + г +1, = г +I + 2 = 2г - 3 - * - ц^ - п д + 2 .
цр
Так как полиномы А(г) и Б(г) не имеют общих множителей, то система (37) будет иметь единственное решение, если = Му , т. е. если
г = пА + * + ц,цр -1, г = пА + ц,цо -3-1, п = яд + V,
цр
'д
(38)
I = пд + * -1, п& = п^ + пд + цЦр + * -3-1.
Коэффициенты полинома &(г) степени п& из (36), (37) назначаются
с использованием стандартных нормированных передаточных функций (СНПФ) [20, с. 258-260]. С этой целью коэффициенты некоторых СНПФ Уу^ (5) непрерывных систем хранятся в памяти микропроцессора АдЦР.
Как известно, функции У(г), полученные в результате преобразования Хт {V(5)} (4), обычно имеют относительную степень цу г равную нулю или единице. Поэтому степень знаменателя вспомогательной ПФ У(5) берется равной пу = п& - цдоп , где Цдоп = тах{цоу + ^цр - цу,г; 0}; цоу = п2 - 3 -
относительный порядок УЭМР (28). Коэффициенты Др и величина ipXag
подходящей СНПФ V(s) выбираются из памяти микропроцессора по значени-
* *
ям Ç , итаб = nv и с<с , %. Затем коэффициенты знаменателя этой вспомогательной ПФ V(s) вычисляются по формуле
5j ^j®nv~j, j = [nv, 0], (39)
где коэффициент юг = ?таб/(^ - Ty j . В случае астатизма первого порядка
передаточная функция имеет вид V(s) = boj(snv +^n-isnv 1 +••• + ôjs + 5o j ; она подвергается Zj преобразованию (4) при T = Ty, что дает функцию v(z) = Ço (z)/ Ç(z). Подчеркнем: преобразование v(s) ^ v(z) (4) удобно выполнять с помощью функции «c2d» пакета MATLAB [22].
Далее находятся полиномы D(z) = z^допÇ(z) = znD + -1znD 1... +
+ 8iz + ôo и H ( z ) = (z -1)v A q ( z ) = ^n z n + •-- + no; их коэффициенты 5 j, л j, а также P- подставляются в систему (37); ее решение позволяет записать полиномы
R(z) = Xj=oРjzj и L(z) = £j.=0Xjzj . (40)
Затем полиномы R(z) и L(z) записываются по формулам (34), а полином Q(z) из уравнения АдЦР (29) определяется по формуле
Q (z) = h'-1 A Q (z) Ç o( z). (41)
Найденные полиномы Д(z), Q(z) и Л(z) подставляются в уравнение (29), в результате чего находится уравнение искомого АдЦР для идентифицированного УЭМР. Корректность расчетов проверяется с помощью полиномов
A(z ) = H o( z ) Pq ( z ) D( z ), D2(z ) = P (z )Д (z ) + H ( z )Л^ ). (42)
Результаты являются корректными, если коэффициенты при одинаковых степенях z этих полиномов являются близкими друг к другу.
В некоторых случаях при Çg = Çj > 1 полином R(z) = (z -1)vR^z).
В этом случае уравнение АдЦР (29) записывается в виде
R1(z)Ç(z) = Q(z)g(z) - L(z)y(z), (z -1)vu(z) = Ç(z), (43)
где z ) - z-изображение вспомогательной переменной ^. Алгоритм работы микроконтроллера по вычислению значений управляющего воздействия uk, k = o, 1, 2, •.., формируется путем перехода в уравнениях (43) к оригиналам
uk, uk-v, \k-v, Sk-v, yk-v .
Приведенные выражения (28)-(43) являются соотношениями аналитического алгоритма синтеза АдЦР для каждого канала gj ^ y УЭМР (1). Для автоматической реализации этого алгоритма параметры необходимых СНПФ должны храниться в памяти микропроцессора АдЦР. Необходимые перекрестные связи между различными каналами УЭМР могут быть обеспечены, как показано в работе [24].
Покажем работоспособность предложенных алгоритмов на численном примере.
5. ПРИМЕР
Пусть в результате изучения некоторого неопределенного УЭМР установлено, что порядок его модели не превышает 4, т. е. nmax = 4; он имеет
один выход y(t) и два входа: ui по управлению u(t) и U2 по возмущению f (t) . Если Pmax = 2, то Nд = 2nmax + Pmax = 10; если Ти = 1 с, то время измерения tH = NДГи = 10 с; допустимые значения пробных управлений
Uio = U20 = 3 ; ^0 = 0,01; величина = 10 .
Предлагаемый алгоритм идентификации и синтеза адаптивной системы управления включает следующие пункты.
П. 1. Пробное воздействие U1 (t) = 3 • 1(t) при U2 (t) = 0 подается на объект, который находится в установившемся режиме, несколько раз. Отклонения выходной переменной y (t) измеряется с периодом Ти = 1 с, полученные значения усредняются. Усредненные значения y(k) = y^ при U(t) = U1 (t) и значения марковских параметров Д^, q = k, вычисленные по выражениям (10), приведены в табл. 1. Там же приведены значения y2(k) = y2k и Д2, найденные при U2 (t) = 3 • 1(t) и U1 (t) = 0 .
Таблица 1 Table 1
Отсчеты выходных переменных и марковские параметры, Ти = 1 c Samples of output variables and Markov parameters, Ти = 1 s
k 0 1 2 3 4 5 6 7 8 9 10
y1k 0 0,2017 0,3784 0,6547 1,0525 1,5988 2,3264 3,2755 4,4951 6,0452 7,9988
y2 k 4,2 5,2091 5,8091 6,4091 7,0091 7,6091 8,2091 8,8091 9,4091 10,209 10,609
Д k 0 0,0672 0,0589 0,0921 0,1326 0,1821 0,2425 0,3164 0,4065 0,5167 0,6512
д 2 1,4 0,3364 0,2000 0,2000 0,2000 0,2000 0,2000 0,2000 0,2000 0,2000 0,2000
П. 2. При р^ = ; Р = 0, V = 1, 2, 3, 4 вычисляются определители матрицы по (21). Это дает М0^ = 0,0672, ^ М°0,2 = 0,0027,
det M 01'3 = -2,8868-10 5, detMj'4 = -8.6722-10—9 « 0. Отсюда из формул (20) следует, что «1 = 3. П. 3. При i = 1, «1 = 3, р = 0 по (22) и (23) вычисляется вектор ^3= [-3,4072-10—10 1,2214 -2,2214]т. Так как значе-
3 —10 ""
ние л 10 = 3,4072-10 « Д0 = 0,01, то условие (5) не выполняется, поэтому
ДПО (2) при Т = 1 с является некорректным, в связи с чем выполняется переход на п. 7. П. 7. По формуле Ти = Ти / ^т = 0,1c; переход на п. 1. П. 1. Так как время однократного измерения реакции УЭМР теперь ¿и = 10Ти = 1c, то принимаются Ню = 15, a U20 = 0,5 . Повтор действий, описанных выше, дает новые округленные значения, приведенные в табл. 2.
Таблица 2 Table 2
Отсчеты выходных переменных и марковские параметры Ти = 0,1 c Samples of output variables and Markov parameters, Ти = 0,1 s
k 0 1 2 3 4 5 6 7 8 9 10
У1к 0 0,5045 0,5918 0,6360 0,6785 0,7241 0,7733 0,8262 0,8830 0,9436 1,0083
У 2 k 0,7 0,7706 0,7873 0,7981 0,8082 0,8182 0,8282 0,8382 0,8482 0,8582 0,8682
100* p 0 3,3632 0,5821 0,2946 0,2837 0,1821 0,2425 0,3164 0,4065 0,5167 0,6512
Й 1,4 0,1413 0,0334 0,0215 0,0202 0,0200 0,0200 0,0200 0,0200 0,0200 0,0200
П. 2. Определители матрицы M f'v теперь вычисляются по выражению (21) и данным табл. 2 при ц^ = ц1 ; р = 0, v = 1, 2, 3, 4. Это дает: det M 0'1 = 0,0336, detM°0,2 = 6,5177 • 10-5, detMO^ = -8,5611 -10"10, det M0,4 = -2,7667 • 10"26 . Отсюда из формул (20) также следует, что «1 = 3 .
П. 3. Вектор л 1 = [-0,113 1,244 -2,131] вычисляется по выражениям (23) при
Л ~ 3 ~
i = 1, «1 = 3, р = 0. Значение л 10 = 0,113 >Д0 = 0,01, т. е. условие (5) выполняется, и соответствующий ДПО (2) при Ти = 0,1 с является корректным; переход на п. 4. П. 4. вычисляются по выражению (24) оценки коэффициентов: у30 = 0,0324, у31 = -0,0658, у32 = 0,0336, у33 = 0. Следовательно, по (17) оценка передаточной функции W\(z, Ти, щ) имеет вид
z; 0,1; 3) = -3
0,0336z2 -0,0658z + 0,0324
z3 -2,131z2 + 1,244z -0,113
П. 5. ПФ W (z, 0,1; 3) подвергается преобразованию Zt , используя функцию «d2c» из MATLAB с расширением «'zoh'» при Т = 0,1 с. Для наглядности приведем соответствующие команды:
bz1 = [0.0336 -0.0658 0.0324]; dz1 = [1 -2.131 1.244 -0.113]; sysd11= tf(bz1,dz1,0.1); sysc1 = d2c(sysd1,,zoh') Enter
По данным MATLAB записывается оценка
Wi(s) =
0,8s 2 + 0,3s + 0,4
s3 + 21,8s2 -4,4s + 3,54-10"10 '
(44)
П. 6. Аналогично выполняются шаги 2-5 при Т = 0,1 с, 1 = 2 и данным табл. 2, что дает оценку
W2( s) =
1,4sz - 27,03s - 7,8 s2 + 22s + 3,908 • 10-13
(45)
переход на п. 8.
П. 8. Оценки характеристического полинома А (я), порядка п и передаточной матрицы Ж(s) идентифицированного УЭМР определяются по формуЛ 3 Л 3
лам (26) и (27) при очевидных условиях ц о =112 о = 0 в выражениях (44) и (45): Л(5) = 53 + 21,852 - 4,4,ч, П = 3, а (¿^) = 5 - 0,2;
W (s) =
0,8s2 + 0,3s + 0,4
s 3 + 21,8s 2 - 4,4s
1,4s3 - 27,31s2 -2,394s +1,56 s 3 + 21,8s 2 - 4,4s
(46)
Синтез адаптивного регулятора. Найти разностные уравнения АдЦР (29) с периодом Ту = 0,6 с и ц^ = 1 для УЭМР (46). Система управления должна
иметь ^ = С/ = 1, Р*= 0,32; длительность переходного процесса
не более ?пп = 5,6 с; перерегулирование не более о = 5 %; переменные у£ и gk измеряются.
П. 9. Результат идентификации - текущая ПМ (46) объекта подвергается 2т -преобразованию с помощью функции «с2^> МЛТЬЛБ при Ту = 0,6, что
дает
-Т
W (z) =
0,05155z -0,08458z + 0,03998 z3 -2,128z2 + 1,128z-2,087•Ю-6
1,4z3 - 5,804z2 + 7,376z - 2,946 z3 -2,128z2 + 1,128z -2,087 • 10-6
(47)
В (47) малый коэффициент -2,087 -10 6 заменяется нулем, и записываются полиномы уравнения (28):
А(г) = г3 -2,128г2 + 1,128г, В(г) = 0,05155(г2 -1,6407г + 0,7756),
(48)
F(z) = 1,4z3 - 5,804z2 + 7,376z - 2,946.
= 0,6016, zl4 = 0, zA = 1,128, zA = 1,0,
П. 10. В данном случае
р(1) Ф 0, т. е. пд = 3, пп = 1, ий= 2, Пр = 0; 3 = 2, р9 = 0,05155,
в —
ц0у = 3-2 = 1; корни условию (32) удовлетворяют; полином А(г)
факторизуется, что дает А г) = г , А г) = г2-2,128г +1,128, п^ = 2; у= тах{1 -1;1 - 0; 0} = 1. Вычисления по формулам (38) и последующим
дают: г = 3 +1 +1-1 = 4, Г = 3 +1-2-1 = 1, п = 2 +1 = 3, I = 2 +1-1 = 2, п^ = 3 + 2 + 1 +1-2-1 = 4. Пусть цу г = 1, тогда цдоп = 1 +1-1 = 1,
«у = 4-1 = 3 .
П. 11. Коэффициенты Д j и величина ?рХад выбираются из памяти мик-
*
ропроцессора по значениям £ = 1, п таб = «у = 3 и с % = 0% < 5%: Д 0 = 1, Д1 = 3, Д2 = 3, Д3 = 1, ?ртаб = 6,31с. Коэффициент юг = 6,31/(5,6-0,6) « « 1,262, а коэффициенты 5j вычисляются по (39): 83 = 1, 52 = 3 -1,292 = = 3,786, 51 = 3 -1,2922 = 4,7779; 50 = 1 1,2923 = 2,0099. Составленная по
-»—г 3 2
этим коэффициентам вспомогательная ПФ у (я) = 2,0099 / (я + 3,7865 + + 4,77795 + 2,0099) подвергается преобразованию (4) при Ту = 0,6 с,
что дает полиномы £0(г) = 0,04147г2 + 0,09497г + 0,01331 и £(г) = г3-2
- 1,40695 ^ + 0,6598г - 0,1031. Коэффициенты знаменателя округлены
по условию астатизма первого порядка [20, с. 151], причем
< 0,4861,
j = 1,3, т. е. Ç(z) = ÇQ (z) и D(z) = zÇ(z) = z4 -1,40695z3 + 0,6598z2-0,1031z.
3 2
В данном случае полином H(z) = z -3,128z + 3,256z-1,128 и = 0,05155, l = 2, r = 1. Решение системы (37), составленной по полученным данным, позволяет с учетом (34), (40) и (41) записать полиномы: L (z ) = 54,0688z -- 88,8233z + 37,6594; R(z) = z + 1,7211;
R(z) = (z -1)(z3 + 0,0804z2 - 2,0482z +1,3349),
L(z) = 54,0688z3 - 88,8233z2 + 37,6594z, (49)
Q(z) = 0,8045z3 +1,8423z2 + 0,2582z .
Найденные полиномы R(z), Q(z) и L(z) подставляются в (29), что дает уравнение искомого АдЦР. Корректность расчетов проверяется с помощью полиномов:
D1 (z) = z7 - 3,0476z6 + 3,7438z5 - 2,2769z4 + 0,6809z3 - 0,08z2;
D2(z) = z7 -3,0476z6 + 3,7436z5 -2,2765z4 + 0,6805z3 -0,0796z2 -0,0001z.
Результаты являются корректными, так как коэффициенты при одинаковых степенях z этих полиномов являются близкими друг к другу.
П. 12. Уравнение АдЦР записывается по (29) с учетом v = 1 и (49):
(z - 1)u (z) = Кz), (z3 - 0,0804z2 - 2,0482z + 1,3349)^(z) = (0,8045z3 + + 1,8423z2 + 0,2582z)g(z) -(54,0688z3 -88,8233z2 + 37,6594z)y(z). (50)
П. 13. Алгоритм работы микроконтроллера по вычислению значений управляющего воздействия uk для УЭМР (46) формируется путем перехода в уравнениях (50) к оригиналам:
%к-1 = 0,0804^-2 + 2,0482^¿-3 -1,3349^-4 + 0,8045gk-1 + 1,8423gk-2 + + 0,2582gk-3 - 54,0688yk-1 + 88,8233yk-2 - 37,6594yfc-3 ;
uk = uk-1 k-1; k =1» 2» 3, ...
Для краткости здесь приведены округленные значения коэффициентов. С целью получения корректных уравнений АдЦР необходимо при вычислениях учитывать не менее 13-15 значащих цифр после запятой.
На рисунке приведена переходная функция синтезированной дискретной системы управления (47) и (50). Как видно из выражения (50) и рисунка, относительный порядок АдЦР, порядки астатизма, перерегулирование и длительность переходного процесса синтезированной системы соответствуют заданным значениям.
Slcp Ro »ролю
0« joe
и
О 2 4 в 8 юТ, s
График переходной функции синтезированной системы Synthesis system transition function graph
Приведенные выражения (28)-(43) позволяют по идентифицированным моделям аналитически синтезировать алгоритм АдЦР для каждого канала gi ^ yi УЭМР (1), если необходимые СНПФ хранятся в памяти микропроцессора АдЦР.
ЗАКЛЮЧЕНИЕ
Предложенный в данной работе метод позволяет создавать адаптивные системы управления полными многомерными неопределенными элементами микророботов и другими объектами. Приведенные в статье расчетные соотношения охватывают как процедуру идентификации текущих значений порядка и параметров моделей, так и процедуру синтеза цифровых регуляторов по этим моделям. При этом разрешимость задачи идентификации обеспечивается автоматическим выбором периода дискретизации выходных переменных. Получаемые цифровые регуляторы имеют минимальную сложность за счет согласования полюсов системы с достаточно левыми нулями и полюсами микророботов и обеспечивают заданные порядки аста-тизма к задающим воздействиям и внешним возмущениям. Алгоритм расчета цифровых регуляторов является полностью аналитическим, и после настройки под конкретный микроробот реализуется автоматически. При этом используются минимальная информация о микророботах, а также стандартные нормированные передаточные функции, которые хранятся в памяти микроконтроллера.
Предложенный подход может применяться для создания адаптивных систем управления не только неопределенными исполнительными системами микророботов, но и неопределенными объектами промышленного, сельскохозяйственного и специальных производств, а также для решения различных задач управления в социальной сфере.
СПИСОК ЛИТЕРАТУРЫ
1. ИвановА.В., Юревич Е.И. Мини- и микроробототехника: учебное пособие. - СПб.: Изд-во Политехн. ун-та, 2011. - 96 с.
2. Еремин Е.Л., Теличенко Д.А., Чепак Л.В. Модели современных систем автоматического управления: учебное пособие. - Благовещенск: Изд-во АмГУ, 2014. - 66 с.
3. Acar C., ShkelA. MEMS vibratory gyroscopes: structural approaches to improve robustness. - New York: Springer Science & Business Media, 2008. - 260 p.
4. A neuromuscular interface for robotic devices control / I. Kastalskiy, V. Mironov, S. Lobov, N. Krilova, A. Pimashkin, V. Kazantsev // Computational and Mathematical Methods in Medicine. -2018. - Art. 8948145. - P. 1-8. - doi: 10.1155/2018/8948145.
5. Concurrent learning adaptive control of linear systems with exponentially convergent bounds / G. Chowdhary, T. Yucelen, M. Muhlegg, E. Johnson // International Journal on Adaptive Control and Signal Processing. - 2013. - Vol. 27, N 4. - P. 280-301.
6. Круглов С.П. Адаптивное управление неминимально-фазовым скалярным объектом второго порядка с обеспечением заданных характеристик переходного процесса // Научный вестник НГТУ. - 2016. - № 4. - С. 33-53.
7. Гайдук А.Р., Плаксиенко Е.А. Адаптивные системы управления: учебное пособие / Южный федеральный университет. - Ростов н/Д.; Таганрог: Изд-во ЮФУ, 2018. - 120 с.
8. Zhu Y., Hou Z. Controller Dynamic Linearisation-based model-free adaptive control framework for a class of non-linear system // IET Control Theory and Applications. - 2015. - Vol. 9, iss. 7. - P. 1162-172.
9. Красовский А.А., НаумовА.И. Аналитическая теория самоорганизующихся систем управления с высоким уровнем искусственного интеллекта // Известия РАН. Теория и системы управления. - 2001. - № 1. - С. 69-75.
10. Айда-заде К.Р., Рагимов А.Б. Идентификация кусочно-постоянных параметров процесса фильтрации и границ областей их постоянства // Автоматика и телемеханика. - 2017. -№ 8. - С. 60-75.
11. Identification, estimation and control for linear systems using measurements of higher order derivatives / Z. Shao, C. Zheng, D. Efimov, W. Perruquetti // Journal of Dynamic Systems, Measurement, and Control. - 2017. - Vol. 139 (12). - P. 1-6.
12. Jayawardhana R.N., Ghosh B.K. Kalman filter based iterative learning control for discrete time MIMO systems // Proceedings of 30th Chinese Control and Decision Conference. - Shenyang, China, 2018. - P. 2257-2264.
13. Калман Р., ФалбП., АрбибМ. Очерки по математической теории систем. - М.: Мир, 1971. - 237 с.
14. Воскобойников Ю.Е., КрысовД.А. Двухэтапный устойчивый алгоритм непараметрической идентификации системы с высоким уровнем шума входного сигнала // Научный вестник НГТУ. - 2019. - № 1 (74). - С. 21-40. - DOI: 10.17212/1814-1196-2019-1-21-40.
15. ВоеводаА.А., Бобобеков К.М. Активная идентификация параметров модели перевернутого маятника по углу при подаче на вход синусоидальных сигналов // Сборник научных трудов НГТУ. - 2016. - № 2 (84). - С. 21-37.
16. ГайдукА.Р., Плаксиенко Е.А. Анализ и синтез систем управления с применением системных инвариантов. - Saarbrücken: Palmarium Academic Publishing, 2016. - 128 с.
17. Марковские параметры многомерных динамических систем управления / И.А. Каляев, А.Р. Гайдук, С.Г. Капустян, В.Н. Рябченко // Вестник ИГЭУ. - 2013. - Вып. 1. - С. 45-51.
18. Chen C.T. Linear systems theory and design. - New York: Oxford University Press, 1999. -
334 p.
19. Гайдук А.Р. Синтез систем автоматического управления по передаточным функциям // Автоматика и телемеханика. - 1980. - № 1. - С. 11-16.
20. Гайдук А.Р., Плаксиенко Е.А. Анализ и аналитический синтез цифровых систем управления: монография. - СПб.: Лань, 2018. - 199 с. - ISBN 978-5-8114-2813-7.
21. Гайдук А.Р., Капустян С.Г., Шаповалов И.О. Дискретное адаптивное управление движением группы мобильных роботов в неорганизованной среде // Системный анализ, управление и обработка информации: труды IX Международной научной конференции (Дивномор-ское, 20-26 сентября 2018 г.). - Ростов н/Д., 2018. - C. 56-67.
22. Дьяконов В.П. MATLAB R2006/2007/2008 + Simulink 5/6/7. Основы применения. -2-е изд., перераб. и доп. - М.: Солон-пресс, 2008. - 800 с.
23. Гантмахер Ф.Р. Теория матриц. - М.: Наука, 1989. - 549 с.
24. ГайдукА.Р., Колоколова К.В. Синтез систем автоматического управления неустойчивыми многомерными объектами // Научный вестник НГТУ. - 2017. - № 1 (66). - С. 26-40.
Гайдук Анатолий Романович, доктор технических наук, профессор Южного федерального университета, действительный член Академии инженерных наук имени А.М. Прохорова, член IEEE. Основное направление научных исследований - теория систем автоматического управления, анализ и синтез, робототехника. Имеет более 350 научных публикаций, в том числе 19 монографий. E-mail: [email protected]
Капустян Сергей Григорьевич, доктор технических наук, ведущий научный сотрудник Южного научного центра РАН. Основное направление научных исследований -групповая робототехника, многопроцессорные информационно-управляющие и вычислительные системы. Имеет более 280 научных публикаций, в том числе 5 монографий. E-mail: [email protected]
GaidukAnatoly Romanovich, D.Sc. (Eng.), professor at the Southern Federal University; head of the department in the Kislovodsk Humanities Technical Institute; full member of the Academy of Engineering Sciences named after A.M. Prokhorov. The main field of his scientific research is the theory of automatic control systems and robotic systems, analysis and synthesis. He has over 330 scientific publications including 19 monographs. E-mail: [email protected]
Kapustyan Sergey Grigorievich, PhD (Eng.), a leading research fellow in the UNSC RAS; head of the Research Institute of the MVS named after Academician A.V. Kalyaev in SFU. The main field of his scientific research is group robotics, multiprocessor information management and computer systems. He has over 280 scientific publications including 5 monographs. E-mail: [email protected]
DOI: 10.17212/1814-1196-2019-3-33-56 Analytical design of micro-robot executive systems under uncertainty
A.R. GAIDUK1, S.G. KAPUSTYAN2
1 Southeren Federal University, 44 Nekrasovskiy lane, Taganrog, 347922; Russian Federation
2 Federal Research Centre the Southern Scientific Center of the Russian Academy of Sciences, 41 Chekhov Street, Rostov-on-Don, 344006, Russian Federation
a [email protected] b [email protected]
Abstract
Recently parallel structure mechanisms have been more and more often used instead of console type manipulators. These mechanisms allow achieving higher speeds, accelerations and accuracy of movements of work tools caused by such technological operations as laying- out and laser cutting of sheet materials, assembly, welding and many others. In this connection the problem of developing corresponding design methods of control becomes quite relevant. The known design methods results in rather complex control systems both from the point of view of their design and their realization. One of the challenges is taking into account the working area limitation of parallel structure mechanisms, especially when designing control optimal in terms of the square-law criteria. The design method of multi-agent control for parallel structure mechanisms is suggested on the example of tripod control. Control laws and algorithms are generated on the basis of a decoupled approach and a method of analytical system design with control by output and impacts (ASDCOI). Application of this analytical method allows creating multi-agent control systems for parallel structure mechanisms with the necessary quality parameters in transient and steady state modes.
In the given paper the known mathematical tripod model is used. It is shown that desirable trajectories of work tool movements should be formed based on the developed manipulator tasks. The computer simulation results confirm the efficiency of the suggested approach. The developed approach can be applied to create multi-agent manipulators with parallel structure mechanisms for solving a wide scope of problems such as application of fertilizers, environmental observations, improving quality of the system machine vision, and carrying out many industrial and technological operations.
Keywords: Console manipulator, parallel structure mechanism, tripod, movement trajectory, control, system, design, decoupled approach, control by output and impacts
*
Received 15 February 2019. The study was carried out with the financial support of the Russian Federal Property Fund in the framework of scientific projects No. 18 58 00051, No. 17 29 07054, No. 16-29-04194 and No. 18-07-00637, as well as in the framework of the implementation of the Basic Research Program of the Russian Academy of Sciences in priority areas determined by the Presidium of the Russian Academy of Sciences No. 7 "New developments in promising areas of energy, mechanics and robotics", No. gr. Project AAAA-A18-118020190041-1).
REFERENCES
1. Ivanov A.V., Yurevich E.I. Mini- i mikrorobototekhnika [Mini and microrobotics]. St. Petersburg, SPbPU Publ., 2011. 96 p.
2. Eremin E.L., Telichenko D.A., Chepak L.V. Modeli sovremennykh sistem avtomaticheskogo upravleniya [Models of modern automatic control systems]. Blagoveshchensk, AmSU Publ., 2014. 66 p.
3. Acar C., Shkel A. MEMS vibratory gyroscopes: structural approaches to improve robustness. New York, Springer Science & Business Media, 2008. 260 p.
4. Kastal'skiy I., Mironov V., Lobov S., Krilova N., Pimashkin A., Kazantsev V. A neuromuscular interface for robotic devices control. Computational and Mathematical Methods in Medicine, 2018, art. 8948145, pp. 1-8. doi: 10.1155/2018/8948145.
5. Chowdhary G., Yucelen T., Muhlegg M., Johnson E. Concurrent learning adaptive control of linear systems with exponentially convergent bounds. International Journal on Adaptive Control and Signal Processing, 2013, vol. 27, no. 4, pp. 280-301.
6. Kruglov S.P. Adaptivnoe upravlenie neminimal'no-fazovym skalyarnym ob"ektom vtorogo poryadka s obespecheniem zadannykh kharakteristik perekhodnogo protsessa [Adaptive control of a non-minimal-phase scalar object of the second order with the maintenance of the preset transient characteristics]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta -Science bulletin of the Novosibirsk state technical university, 2016, no. 4, pp. 33-53.
7. Gaiduk A.R., Plaksienko E.A. Adaptivnye sistemy upravleniya [Adaptive management systems: a tutorial]. Rostov-on-Don, Taganrog, YuFU Publ., 2018. 120 p.
8. Zhu Y., Hou Z. Controller Dynamic Linearisation-based model-free adaptive control framework for a class of non-linear system. IET Control Theory and Applications, 2015, vol. 9, iss. 7, pp. 1162-172.
9. Krasovskii A.A., Naumov A.I. Analiticheskaya teoriya samoorganizuyushchikhsya sistem uprav-leniya s vysokim urovnem iskusstvennogo intellekta [Analytical theory of self-organizing control systems with a high level of artificial intelligence]. Izvestiya Rossiiskoi akademii nauk. Teoriya i sistemy upravleniya - Journal of Computer and Systems Sciences International, 2001, no. 1, pp. 69-75. (In Russian).
10. Aida-zade K.R., Rahimov A.B. Identifikatsiya kusochno-postoyannykh parametrov protsessa fil'tratsii i granits oblastei ikh postoyanstva [Identification of piecewise constant filtration parameters and boundaries of their constancy domains]. Automation Remote Control, 2017, vol. 78, iss. 8, pp. 1404-1416. Translated from Avtomatika i telemekhanika, 2017, no. 8, pp. 60-75.
11. Shao Z., Zheng C., Efimov D., Perruquetti W. Identification, estimation and control for linear systems using measurements of higher order derivatives. Journal of Dynamic Systems, Measurement, and Control, 2017, vol. 139 (12), pp. 1-6.
12. Jayawardhana R.N., Ghosh B.K. Kalman filter based iterative learning control for discrete time MIMO systems. Proceedings of 30th Chinese Control and Decision Conference, Shenyang, China, 2018, pp. 2257-2264.
13. Kalman R., Falb P., Arbib M. Topics in mathematical system theory. New York, McGraw-Hill, 1969 (Russ. ed.: Kalman R., Falb P., Arbib M. Ocherki po matematicheskoi teorii sistem. Moscow, Mir Publ., 1971. 237 p.).
14. Voskoboinikov Yu.E., Krysov D.A. Dvukhetapnyi ustoichivyi algoritm neparametricheskoi identifikatsii sistemy s vysokim urovnem shuma vkhodnogo signala [Two-step robust algorithm of nonparametric identification for a system with a high-noise input signal]. Nauchnyi vestnik Novosibir-skogo gosudarstvennogo tekhnicheskogo universiteta - Science bulletin of the Novosibirsk state technical university, 2019, no. 1 (74), pp. 21-40. DOI: 10.17212/1814-1196-2019-1-21-40.
15. Voevoda A.A., Bobobekov K.M. Aktivnaya identifikatsiya parametrov modeli perevernutogo mayannika po uglu pri podache na vkhod sinusoidal'nykh signalov [Active identification of the inverted pendulum model data on angle in applied to the input sinusoidal signal]. Sbornik nauchnykh
trudov Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta - Transaction of scientific papers of the Novosibirsk state technical university, 2016, no. 2 (84), pp. 21-37.
16. Gaiduk A.R., Plaksienko E.A. Analiz i sintez sistem upravleniya s primeneniem sistemnykh invariantov [Analyses and synthesis of control systems with application of system invariant in MATLAB]. Saarbrucken, Palmarium Academic Publishing, 2016. 128 p.
17. Kalyaev I.A., Gaiduk A.R., Kapustyan S.G., Ryabchenko V.N. Markovskie parametry mnogomernykh dinamicheskikh sistem upravleniya [Markov parameters of multivariable dynamic control systems]. Vestnik Ivanovskogo gosudarstvennogo energeticheskogo universiteta - Vestnik of Ivanovo State Power Engineering University, 2013, iss. 1, pp. 45-51.
18. Chen C.T. Linear systems theory and design. New York, Oxford University Press, 1999. 334 p.
19. Gaiduk A.R. Sintez sistem avtomaticheskogo upravleniya po peredatochnym funktsiyam [Synthesis of control systems on transfer functions]. Automation and Remote Control, 1980,vol. 41, iss. 1, pp. 6-11. Translated from Avtomatika i telemekhanika, 1980, no. 1, pp. 11-16.
20. Gaiduk A.R., Plaksienko E.A. Analiz i analiticheskii sintez tsifrovykh sistem upravleniya [Analysis and analytical design of digital control systems]. St. Petersburg, Lan' Publ., 2018. 199 p. ISBN 978-5-8114-2813-7.
21. Gaiduk A.R., Kapustyan S.G., Shapovalov I.O. [Discrete adaptive control of a group of mobile robots motion in an unorganized environment]. Sistemnyi analiz, upravlenie i obrabotka infor-matsii: trudy IX Mezhdunarodnoi nauchnoi konferentsii [System analysis, control and information processing: Proceedings of the IX International Scientific Conference], Divnomorskoe, September 20-26, 2018. Rostov-on-Don, 2018, pp. 56-67. (In Russian).
22. D'yakonov V.P. MATLAB R2006/2007/2008 + Simulink 5/6/7. Osnovy primeneniya [MATLAB R2006/2007/2008 + Simulink 5/6/7. Application basics]. 2nd ed. Moscow, Solon-press Publ., 2008. 800 p.
23. Gantmakher F.R. Teoriya matrits [Matrix theory]. Moscow, Nauka Publ., 1989. 549 p.
24. Gaiduk A.R., Kolokolova K.V. Sintez sistem avtomaticheskogo upravleniya neustoichivymi mnogomernymi ob"ektami [Synthesis of control systems by unstable multivariable plants]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta - Science bulletin of the Novosibirsk state technical university, 2017, no. 1 (66), pp. 26-40.
Для цитирования:
Гайдук А.Р., Капустян С.Г. Аналитический синтез исполнительных систем микророботов в условиях неопределенности // Научный вестник НГТУ. - 2019. - № 3 (76). - С. 33-56. -DOI: 10.17212/1814-1196-2019-3-33-56.
For citation:
Gaiduk A.R., Kapustyan S.G. Analiticheskii sintez ispolnitel'nykh sistem mikrorobotov v uslovi-yakh neopredelennosti [Analytical design of micro-robot executive systems under uncertainty]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta - Science bulletin of the Novosibirsk state technical university, 2019, no. 3 (76), pp. 33-56. DOI: 10.17212/1814-11962019-3-33-56.
ISSN 1814-1196, http://journals.nstu.ru/vestnik Science Bulletin of the NSTU Vol. 76, No 3, 2019, pp. 33-56
