CC BY
19
УПРАВЛЕНИЕ, МОДЕЛИРОВАНИЕ, АВТОМАТИЗАЦИЯ
DOI: 10.24143/2072-9502-2017-4-7-20 УДК 681.513
А. А. Воевода, К. М. Бобобеков
СИНТЕЗ ЛИНЕЙНЫХ МНОГОКАНАЛЬНЫХ РЕГУЛЯТОРОВ С ИСПОЛЬЗОВАНИЕМ СТРУКТУРНЫ1Х ПРЕОБРАЗОВАНИЙ
Рассматривается задача использования левого и правого матричных полиномиальных разложений при синтезе регуляторов модальным методом. Предполагается, что задано левое матричное полиномиальное разложение. При синтезе регуляторов сложность вычислений существенно возрастает при неравенстве столбцовых/строчных степеней полиномиальной матрицы «знаменателя» объекта. После перехода от левого матричного полиномиального разложения к правому предлагается осуществить ряд унимодальных преобразований матрицы «знаменателя» объекта, которые сводятся к введению унимодальных матриц в различные участки структурной схемы. При этом осуществляются эквивалентные преобразования, которые «выравнивают» столбцовые/строчные степени полиномиальной матрицы «знаменателя» объекта (для сохранения эквивалентности вводятся сразу две унимодальные матрицы, взаимно обратные друг другу). Предлагаемая методика иллюстрируется на примере двухка-нальной системы, соответствующей динамической системе, которая, в свою очередь, соответствует трехмассовой системе с пружинами.
Ключевые слова: линейные многоканальные системы, системы автоматического управления, синтез регуляторов, матричное полиномиальное представление, структурные преобразования, унимодальные матрицы.
Введение
В последнее десятилетие наряду с традиционными методами синтеза многоканальных линейных систем, таких как синтез в пространстве состояний и на основе матричных передаточных функций, получили развитие методы синтеза на основе полиномиального матричного левого/правого разложения объекта и регулятора [1, 2]. Синтез регуляторов в многоканальных системах при использовании полиномиальных представлений существенно усложняется в случаях, когда столбцовые/строчные степени полиномиальной матрицы знаменателя существенно отличаются. Для упрощения синтеза регуляторов предлагается перед этапом синтеза регулятора произвести серию эквивалентных преобразований полиномиального представления объекта с целью выравнивания столбцовых/строчных степеней полиномиальной матрицы знаменателя. В упрощенном варианте это можно представить как введение единичной матрицы на входе/выходе структурной схемы объекта, а также между матрицами «знаменатель - числитель» объекта. Затем единичную матрицу представляем как произведение полиномиальной матрицы на ее обратную и, далее, эти матрицы объединяем с предыдущей или последующей матрицами. Для того чтобы не выходить из множества полиномиальных матриц при разложении единичной матрицы на полиномиальные матрицы, используем унимодальные матрицы, отличительной особенностью которых является равенство детерминанта вещественному числу. Предлагаемая методика иллюстрируется на примере двухканальной системы, соответствующей динамической системе, которая, в свою очередь, соответствует трехмассовой системе с пружинами.
Постановка задачи
Цель исследования - существенно облегчить процедуру синтеза многоканальных регуляторов, используя унимодальные матрицы (unimodal matrices) как средство, позволяющее осуществлять эквивалентные преобразования объекта и регулятора, выравнивать столбцо-
вые/строчные степени полиномиальных матричных описаний. К такой постановке задачи прийти несложно, если проанализировать примеры синтеза линейных многоканальных систем в случаях, когда вышеупомянутые степени полиномиальной матрицы знаменателя объекта значительно отличаются [3-9]. Для того чтобы преобразования были эквивалентными, при введении унимодальной матрицы необходимо ввести обратную унимодальную матрицу. Использование такого вида эквивалентных преобразований возможно для многоканальных линейных систем, использующих полиномиальное матричное представление. Возможность использования таких преобразований - это отличительная особенность полиномиального представление именно для многоканальных систем. Сложность введения таких унимодальных матриц состоит в необходимости «угадать» вид этих матриц. Эта проблема легко решается, если иметь некоторый опыт по использованию унимодальных матриц. Например, достаточно проанализировать несколько примеров вычисления правого полиномиального разложения по левому разложению [3-11].
Использование унимодальных матриц в многоканальных системах как средство осуществления эквивалентных преобразований
Пусть задан объект в виде левого матричного полиномиального разложения объекта [4-7]:
Dl (s) y (s) = N, (s )u( s). Структурная схема, приведенная на рис. 1, соответствует формуле (1).
(1)
Рис. 1. Структурная схема многоканального объекта, представленная в виде левого матричного полиномиального разложения
Далее, для экономии места и наглядности в формулах и на рисунках, переменные, например q(s), и^), е(s), будем записывать без аргумента 5. По левому разложению найдем правое матричное полиномиальное разложение [3-7, 10]. Для этого составляется блочная матрица из ^ (5) и N (5):
D (s)
N (s)
(2)
над строчками и столбцами которой выполняются элементарные операции, соответствующие умножению матрицы (2) слева на унимодальные матрицы. В результате получим полиномиальную матрицу, являющуюся наибольшим общим матричным полиномиальным делителем Dt (s) и Nl (s). Кроме того, совокупность унимодальных матриц позволяет вычислить правое матричное полиномиальное разложение объекта Dr (s) и Nr (s).
Таким образом, нами получена матричная передаточная функция объекта в виде
Wob (s) = Nr (s) D-'( s), откуда можно записать связь «вход - выход»:
y (s) = Nr ( s) D;\s)u( s).
(3)
Следует отметить, что получено взаимно простое правое разложение. Формуле (3) соответствует структура, приведенная на рис. 2.
Рис. 2. Структурная схема многоканального объекта, представленная в виде правого матричного полиномиального разложения
В качестве объекта возьмем левое матричное полиномиальное разложение:
^2 + 4
-2
-2( ^2 + 4) ^4 +10 ^2 + 8
V л Л
V У 2 у
0
Д (*)
0 * +4^
УМ
V и2 у
соответствующее объекту, который состоит из трех масс и соединяющих их пружин с нулевым коэффициентом демпфирования. Для синтеза необходимо привести правое полиномиальное матричное разложение объекта, для вычисления которого мы воспользовались процедурой, описанной в [4-7].
Г У1 ^ Г16 0 1
V У 2 у V0 *2+4 J
N (*)
16( *2 + 4) -2( *2 + 4)
-32 *4+10*2+8
*)
\-1
V и2 у
Для данного объекта имеем два управляющих сигнала - силы и1 и и2, приложенные к массам щ и т2 , и два выхода, соответствующие координатам у1 и у2 грузов щ и т2 .
В матрице (*) столбцовые/строчные степени равны 2 и 4, что, по-видимому, может приводить к значительному усложнению процедуры расчета регуляторов, обеспечивающих, например, заданное расположение полюсов системы. Попытаемся при помощи структурных преобразований многоканальной линейной системы выравнять столбцовые/строчные степени матрицы (*).
Выравнивание столбцовых/строчных степеней А.^)
В структурную схему объекта (рис. 2) добавим две унимодальные матрицы - R( и R-1(*) . В результате получим структурную схему (рис. 3).
и
и
Рис. 3. Структурная схема объекта после введения единичной матрицы I = R( R 1 (
На рис. 3 блок слева, выделенный штриховой линией, соответствует передаточной функции R, которую обозначим как , т. е.
D1-1(s) = R~\s)D;\s) = ( Dr (s) R( s) )-1.
Для второго выделенного блока (по аналогии) передаточная функция равна Nr (s)R(s), обозначим ее как Nr (s), т. е.
Nr (s) = Nr (s) R( s).
Таким образом, структурная схема на рис. 3 преобразована в структуру, приведенную на рис. 4.
Рис. 4. Структурная схема объекта после введения унимодальных матриц Я( 5) и Я 1 (5)
Проведем масштабирование матриц Бг (5) и N (5), а именно первые столбцы матриц Бг (5) и N (5) поделим на 16 и первые столбцы этих матриц умножим на - 5 и прибавим ко вторым столбцам, что в дальнейшем позволит избавиться от четвертой степени 5 в матрице Бг (5). Соответствующая матрица Я(5) приведена в вычислениях ниже:
nr (s)
/ \ У1
V У 2 У
'16 0 Л/
16-1 -(16-1) s
-к V ((лг, 2
0 s + 4 Д 0
R (s)
16( sz + 4) -2( sz + 4)
VV -32 s4 +10 s2 + 8 у
Dr (s)
-1
V u2 j
R(s)
В результате умножения матриц получим
(y1 W1 -s Vs2 + 4 -s3 -2s2 -4s-8V1 (u1 Л VУ2J 1,0 s2 + 4-2 s4+10s2+2s + 8j Vu2j
Nr (s)
Dr (s)
В структурную схему, представленную на рис. 4, добавим единичную матрицу, соответствующую произведению двух унимодальных матриц:
I = L-1(s) L(s).
Соответствующая структурная схема приведена на рис. 5.
d-1 (s)
u
X
1
Рис. 5. Преобразованная структурная схема на рис. 4 после введения двух унимодальных матриц
Передаточная функция части структурной схемы, выделенной штриховой линией, равна
D;l( 5) г\5) = (L (5) Dr (5))-1 = D;1( 5).
Выше мы перемножили матрицы L(s) и (*) и полученное произведение назвали (*). Таким образом, в матрицу (*) вошли матрицы R(*) и L(s) . В соответствии с рис. 5 передаточная функция объекта после введения унимодальных матриц R(*) и L(s) равна
у (*) = N (*) О-Ч *) ь( * )и (5). Тогда формула (4) преобразуется следующим образом:
(4)
о;'( *)
/ \ /1
У11 '
* 1(( 1 01 Г *2 + 4
У2 у V 0 *2 + 4 XV * 1
— \ —
2*2 - 4* - 8 Л (1 0и1
2
*4 +10*2 + 2* + 8
N (*) ь(*)
После перемножения матриц (5) получим
Г У11
V У2.
- \
0 * + 4,
N00
*2 + 4
— \ —
Ог (*)
2*2 - 4* - 8Л
V* 1 у V и2 у
ь (*)
(5)
3 + 4* - 2 -2*3 + 6*2 -62 + 8у
Дт'( *)
10
цТ)
и2 у
(6)
что соответствует структурной схеме, представленной на рис. 6. Это легко обосновать, если вернуться к рис. 5.
Рис. 6. Структурная схема объекта после выполнения унимодальных операций
Таким образом, в результате ряда унимодальных операций, соответствующих эквивалентным структурным преобразованиям многоканальной системы, которая представлена в виде матричного полиномиального разложения, получен требуемый результат: степени столбцов
и строк матрицы (*) оказались равными трем (6).Таким образом, задача выравнивания степеней решена. Однако описание объекта отличается от правого матричного полиномиального разложения дополнительным множителем Ь(*).
Синтез регулятора
В соответствии с рис. 6 передаточная функция объекта равна
Wob (* ) = Мг (* ) * ) Ь( *).
Временно забудем про матрицу Ь(*) (на входе объекта) и будем считать, что передаточная функция объекта, представленная в виде правого матричного полиномиального представления, равна
Wob (*) = N (*) *).
Пусть матричная передаточная функция регулятора в виде левого матричного полиномиального представления имеет следующий вид:
Жг (* ) = Y *) X (*).
3
и
Как известно, характеристический матричный полином системы (на входе объекта не учитываем матрицу L(s) !) может быть записан следующим образом [3-7]:
(Y(s) ! X(s))
(Dr (s) ^
N~ (s)
= C (s),
(7)
где как С(5) обозначен характеристический матричный полином системы. Если учтем, что объект двухканальный (откуда следует, что и регулятор двухканальный), можем получить матрицы - полиномы знаменателя и числителя регулятора:
Y (s) =
'УпСО У12 (s) У21(s) У22(s) j
, X (s) =
X„(s) X^s)'
X21(s) X22(s) у
где элементы матриц У(5) и X(5)- полиномы. Потребуем автономности каналов системы автоматического управления, т. е. характеристическая матрица системы должна быть диагональной:
4(5) 0 ^
V 0 С22(5)у '
с (s) =
Из уравнений (6) несложно получить матрицы Dr (5) и Ыг (5) в виде матричных полиномов
Dr (s) = D3s3 + D2s2 + D1s + D0, Nr (s) = N2s2 + Nxs + N0.
(8)
Воспользуемся уравнением (6) и запишем матрицы Di и N.-, где i = 0, 3 и j = 0, 2 :
D3 =
(0 -Г\ =
v 1 -2 у
D2 =
(1 -2 ^ 0 6
= (0 -4^ = ( 4 -8 ^
, D1 =
4 -6
, D0 =
-2
8
N 2 =
( 0 0 ^ _
01
N1 =
(0 -П -
00
N0 =
(1 0 04
Учтем размеры матриц в (7) и подставим значения матриц Dr (s) и Nr (s) :
(9)
( Уп( s) У12(s) У21 (s) У22 (s)
x„(s) x12(s) X21(s) X22( s) j
( s2 + 4 -s3 - 2s2 - 4s - 8 ^ s3 + 4s - 2 -2s3 + 6s2 - 6s + 8
-s s2 + 4
(cu(s) 0 ^ 0 c22(s) у
Зададим степени матриц Y(s) и X(s) минимальными:
/* *\ /* *
Y (s) = Y,s + Y0 =
* *
s+
* *
v
X (s) = Xjs + X0 =
f* *
f* *
s+
(10)
(11)
Так как степень полиномиальной матрицы Dr (5) равна трем, а степень полиномиальной матрицы У (5) равна единице, выбираем степень характеристической матрицы равную четырем:
С (s) = C4s4 + С/ + С2 s2 + Qs + С0 =
(* 0 0*
s4 +
(* 0 0*
s3 +
(* 0*
s2 +
(* 0 ^ (* 0 ^
0*
s+
0*
(12)
* *
* *
В формуле (12) звездочками помечены элементы, которые определяют характеристические полиномы первого и второго каналов. Эти полиномы будут заданы ниже. Подставим (8), (10)-(12) в (7):
(7^ + 70 ) ф353 + Д52 + Д 5 + Д) + (Х^ + Х0 ) (Й252 + ^ + Ш0 ) =
= С4/ + С353 + С252 + С15 + С0. (13)
Раскроем уравнение (13):
4 + (7Д + 70 А3 + Х1И2 >3 + (7Д + 70 А2 + + Х0И2 >2 + (7150 + 70 Д + ХД, + X Д > + + 70А0 + X0N0 = С454 + С353 + С252 + Сгч + С0.
Приравняв коэффициенты членов полиномов в левой и правой частях, получим систему линейных уравнений, которую можно записать в матричном виде:
ЗЭТ = К. (14)
Здесь
3 = (7 ¡70 | Х,\ X0), К = (С4 ¡С3 ¡С2 | С ¡С0), (15)
эт =
А А,
А
а
А
А 0
А А
N2 N1 N0 .0
0 N2 N N0,
(16)
Таким образом, нами получена система матричных уравнений с четырьмя матричными неизвестными и пятью матричными уравнениями. Подставим значения Di и N ■ из (9) в (16):
2 3
10
(0 -1 ! 1 -2!0 -4! 4 -8|0 0 ^
эт=
-2
8 ! 0 0
.1.
0 -4! 4 -8 4 -6 ! -2 0
0 4 | 0 0
Т ~-П~Т"0
0 0 | 0 0 | 0 1 | 0 0 | 0 4
(17)
В строчке над матрицей указаны номера столбцов, а в столбце справа от матрицы - номера строчек матрицы. Ранг матрицы ЭТ равен восьми, откуда следует, что есть линейно зависимые столбцы. Если, например, выбросим первый столбик из матрицы ЭТ, то произведение ЗЭТ из матрицы размерами 2 х 10 станет матрицей размерами 2 х 9 и, соответственно, матрица К в правой части уравнения (14) размерами 2 х 10 станет матрицей размерами 2 х 9 . Можно проверить, действительно ли, если у матрицы ЭТ выбросить столбцы 1 и 2, то ранг остается равным восьми. Матрица ЗЭТ имеет размеры 2 х 8. При этом матрица К «укоротится» на первый и второй столб-
4
5
6
7
8
9
цы и будет размерами 2 х 8 . После вычеркивания столбцов у матрицы ЭТ и матрицы К назовем их ЭТХ и соответственно. Таким образом, уравнение (14) преобразовано
в невырожденное (¿еКЭТ^ Ф 0) уравнение
=К1. (18)
Обусловленность и значение детерминанта матрицы ЭТХ удовлетворительные:
со^(ЭТх) = 82,3, ¿е^) = 1024.
В уравнение (18) подставим 3, К из (15) и матрицу , полученную
первого и второго столбцов матрицы ЭТ
3 4 5 6 7 8 9 10
(1 -2 j 0 -4 4 -8 ! 0 0 > 1
0 6 ! 4 J_ -6 -2 8 ! 0 i 0 2
0 -1 ! 1 1 -2 0 -4 ! 4 1 -8 3
[Y1 Y0 \X\ ! x 0 ] 1 -2 ! 0 6 4 -6 ! -2 8 4 (С 1 С 1 С VC3 1 С2 1 С1
ö~ "-Г _____ "0" ¡"о" "Ö 5
0 1 i0 0 0 4 i 0 0 6
0 1 0 0 -1 ! 1 i 0 7
v 0 0 0 1 0 0 i 0 4 у 8
Со).
В строчке над матрицей ЭТХ отсутствуют номера 1 и 2, что соответствует вычеркиванию этих столбцов из матрицы ЭТ. Осталось подставить матрицы С0 — С3.
Зададим полюса у первого и второго контуров равными {—1, — 1, —1}, что соответствует
характеристическим полиномам (я +1)3. Кроме того, потребуем автономности контуров, из чего следует, что желаемая характеристическая матрица системы имеет следующий вид:
вд =
(1 0^ , (3 о'к (3 0
v0 1 у
s3 +
V 0 3 у
s2 +
03
s+
(1 0 01
Найдем параметры регулятора из уравнения (18):
3 = К1ЭТ—11.
В результате вычислений получим
3 =
(0,28125 0,34375 -0,375 0,03125
1,625 0,71875 -0,125 0,375
-0,3125 1,5625 0,0625 0,6875
-4,0625 1,8125 ^ 1,25 -0,75
Нашли матрицы регулятора Х1, Х0 и Yl, Y0. Выпишем полиномиальные матрицы Y(s) и Х(я) с точностью два знака после запятой:
Y (s) =
( 0,28s +1,63 0,34s + 0,72Л ч-0,38s - 0,13 0,03s + 0,38у
X (s) =
(-0,31s - 4,06 1,56s + 1,81 0,06s +1,25 0,69s - 0,75
Кратко эти матрицы могут быть записаны следующим образом:
X(я) = Ххя + Х0 , Y(s) = ^я + Y0 . Найдем полюса регулятора, для чего вычислим определитель матрицы Y(я) :
y
Y
)) = 0Д38я2 + 0,469^ + 0,7,
откуда получаем полюса
-1,70 ± 1,48/ .
Сделаем проверку, а именно, зная матрицы 3 и ^ , вычислим их произведение:
=
(1 0 3 0 3 0 1 01
V0 1 0 3 0 3 0 1,
Другими словами,
= (Сз С2 С1 Со) •
Но = , т. е. мы вычислили характеристическую матрицу системы, в которой числа, существенно меньшие единицы, заменены на нули.
Матрицы Y(s) и Х(5) приведены с точностью два знака после запятой. Такое округление допустимо, т. к. при учете большего числа знаков коэффициенты полинома ёе1(У (5)) и значения корней меняются незначительно.
Вспомним, что на входе объекта была матрица L(s) (рис. 6). Для того чтобы учесть эту матрицу, но сохранить эквивалентность этих схем, на выходе регулятора добавим матрицу ЬТ1^). На структурной схеме (рис. 7) приведен регулятор с учетом унимодальной матрицы Ь_1(5) . Сначала выпишем регулятор с учетом формулы (19):
(^ + Yо)q = (Х^ + X о)е.
Здесь выходной сигнал - q , а если учесть матрицу Ь1 (5), то формула реализации регулятора будет иметь следующий вид:
и = Г1^ )7Х-11 Ххе + 1(-70и + X0е) |.
(20)
Данная схема реализуемая, т. к. det(Y1) = 0,14 Ф 0 . Реализация алгоритма управления (20) приведена на структурной схеме (рис. 7).
Рис. 7. Структурная схема регулятора с учетом Ь х(я)
Последний блок на структурной схеме Ь х(5) вставлен из-за того, что в структурной схеме на рис. 6 слева стоит матрица Ь(5), а для компенсации на выходе регулятора нами было вставлено звено Ь_1(5) . Выходной сигнал регулятора и(5) равен
L x(s)q(s) = u(s) .
(21)
Преобразуем уравнение (21)таким образом, чтобы представить и(5) через линейную комбинацию q(s) и sq(s):
( 1 0 V q11 ( и11
-s 1 Мг
у V 2 у
откуда
•а гил
или
sqi + q2 у V u2 у
Г 0 ^ (q ^
+
qi
V -21 у Vq2 у q
qi
2 у
Для окончательного решения введем числовую матрицу С, тогда
( 0 01(1 ( 0 1
sqi
-1 0уVsq2
УЧ-12У V
sqi
Получили требуемое представление:
u(s) = q(s) + Gsq(s) .
(22)
Преобразуем структурную схему на рис. 7 с учетом (22) и с учетом элементарных структурных преобразований. Преобразованная структурная схема показана на рис. 8.
u
s
u
G
Рис. 8. Структурная схема регулятора с учетом Ь х(я) после структурных преобразований
Вернемся к модели объекта (рис. 1-6). Для моделирования объекта воспользуемся левым разложением. Но ввиду неравенства столбцовых/строчных степеней матрицы Dl (я) матричное структурное представление объекта для левого разложения нереализуемо. Для выравнивания степеней Dl (я) можно воспользоваться эквивалентными структурными преобразованиями, аналогичными преобразованиям над Ог (я) . После ряда эквивалентных преобразований можно получить еще один вид структуры, соответствующей данному объекту (рис. 9).
Рис. 9. Структурная схема объекта в виде левого матричного полиномиального представления после ряда эквивалентных преобразований
Матричное структурное моделирование объекта стало возможным после выравнивания столбцовых степеней. Другими словами, det(D3) Ф 0. Значения матриц Q, ) и Nимеют следующие значения:
Q =
(о 1 1
V0
)з =
(0 1 -1 о
А =
N =
(1
V0 6, (0 0 1 -1 0
) =
( 0 4 -4 2
) =
(4 -2 04
N 0 =
(1 0 24
NN 2 =
( 0 01 V 0 1,
В структурной схеме системы (рис. 10) для моделирования объекта используется структура, приведенная на рис. 9, а для моделирования регулятора - на рис. 8.
Рис. 10. Структурная схема системы регулирования Переходные процессы системы на рис. 10 для v(t) = (1 0)' и v(t) = (0 1)' приведены на рис. 11.
Рис. 11. Переходные процессы системы
На рис. 11 графики, отмеченные цифрой 1, соответствуют первому выходу системы, а цифрой 3 - второму выходу системы. Для сравнения приведен выход системы с передаточной
функцией + 1)3 по каждому каналу. Из рис. 11 следует, что время переходного процесса соответствует заданному. Следует обратить внимание на то, что статический режим является неудовлетворительным ввиду статизма системы. Для удовлетворительного поведения системы в статическом режиме надо либо ввести интегратор в регулятор, либо повысить степень регулятора с целью обеспечить достаточно большой коэффициент усиления разомкнутой системы.
Заключение
При решении задачи синтеза регуляторов в многоканальных линейных системах при использовании полиномиальных матричных разложений, как следует из анализа примеров расчета регуляторов, задача существенно усложняется при существенном различии столбцовых степеней матрицы знаменателя объекта. Чтобы упростить решение этой задачи, предложено использовать новый вид структурных преобразований многоканальных систем, сводящийся к введению унимодальных матриц в различные участки структурной схемы. Для сохранения эквивалентности вводятся сразу две унимодальные матрицы, взаимно обратные друг другу. Унимодальные матрицы при работе с матричными полиномами используются давно (например, вычисление правого полиномиального разложения по левому и наоборот). Работ с использованием унимодальных матриц с целью эквивалентных структурных преобразований обнаружено не было, что подтверждает актуальность предпринятого исследования.
СПИСОК ЛИТЕРА ТУРЫ
1. Гайдук А. Р., Колоколова К. В. Синтез систем автоматического управления неустойчивыми многомерными объектами // Науч. вестн. Новосибирск. гос. техн. ун-та. 2017. № 1 (66). С. 26-40.
2. Гайдук А. Р. Теория и методы аналитического синтеза систем автоматического управления (полиномиальный подход). М.: Физматлит, 2012. 360 с.
3. Воевода А. А. Матричные передаточные функции (Основные понятия). Конспект лекций по курсу «Проектирование систем управления» для 4-5 курсов АВТФ; под ред. А. С. Вострикова. Новосибирск: НГТУ, 1994. 94 с.
4. Воевода А. А., Шоба Е. В. О разрешимости задачи автономизации многоканальной системы. Ч. 1 // Сб. науч. тр. НГТУ. 2010. № 2 (60). С. 9-16.
5. Воевода А. А., Шоба Е. В. О разрешимости задачи автономизации многоканальной системы. Ч. 2 // Сб. науч. тр. НГТУ. 2010. № 3 (61). С. 41-50.
6. Воевода А. А., Шоба Е. В. О разрешимости задачи автономизации многоканальной системы. Ч. 3 // Сб. науч. тр. НГТУ. 2010. № 4 (62). С. 3-12.
7. Воевода А. А., Шоба Е. В. О разрешимости задачи автономизации многоканальной системы. Ч. 4 // Сб. науч. тр. НГТУ. 2011. № 3 (65). С. 11-18.
8. Шоба Е. В. Модальный метод синтеза многоканальных динамических систем с использованием полиномиального разложения: дис. ... канд. техн. наук. Новосибирск: НГТУ, 2013.192 с.
9. Вороной В. В. Полиномиальный метод расчета многоканальных регуляторов пониженного порядка: дис. ... канд. техн. наук. Новосибирск: НГТУ, 2013. 173 с.
10. Шоба Е. В. Расчет двухканального регулятора для стабилизации положения кабины лифта с учетом жесткости каната и двумя управляющими воздействиями // Сб. науч. тр. НГТУ. 2012. № 4 (70). С. 45-52.
11. Шоба Е. В. Расчет многоканального регулятора для поддержания заданной температуры в камере полимерной покраски // Сб. науч. тр. НГТУ. 2012. № 4 (70). С. 53-60.
Статья поступила в редакцию 21.06.2017
ИНФОРМАЦИЯ ОБ АВТОРАХ
Воевода Александр Александрович - Россия, 630073, Новосибирск; Новосибирский государственный технический университет; д-р техн. наук, профессор; профессор кафедры автоматики; ucit@ucit.ru.
Бобобеков Курбонмурод Мулломиракович - Россия, 630087, Новосибирск; Новосибирский государственный технический университет; аспирант кафедры автоматики; kurbon_111@mail.ru.
A. A. Voevoda, K. M. Bobobekov
SYNTHESIS OF LINEAR MULTICHANNEL REGULATORS USING STRUCTURAL TRANSFORMATIONS
Abstract. The paper considers the problem of using left and right matrix polynomial decompositions in the synthesis of regulators using modal method. The left matrix polynomial expansion is assumed. In the synthesis of regulators, the complexity of computations essentially increases when the inequality of the column/row degrees of the polynomial matrix of the "denominator" of the object. After transition from the left matrix polynomial expansion tothe right one, the authors propose to accomplish a number of unimodal transformations of the matrix of the object's "denominator", which will result in bringing unimodal matrices into different sections of a structure chart. In this case, equivalent transformations are performed. They "equalize" the column/row degrees of the polynomial matrix of the "denominator" of the object (in order to retain equivalence there are introduced at once two unimodal matrices, opposite to each other). The proposed technique is illustrated by the example of a two-channel system corresponding to a dynamical system, which, in its turn is corresponding to a three-mass system with springs.
Key words: linear multi-channel systems, automatic control systems, regulator synthesis, matrix polynomial representation, structural transformations, unimodal matrices.
REFERENCES
1. Gaiduk A. R., Kolokolova K. V. Sintez sistem avtomaticheskogo upravleniia neustoichivymi mnogo-mernymi ob"ektami [Synthesis of control systems by unstable multivariable plants]. Nauchnyi vestnikNovosibir-skogo gosudarstvennogo tekhnicheskogo universiteta, 2017, no. 1 (66), pp. 26-40.
2. Gaiduk A. R. Teoriia i metody analiticheskogo sinteza sistem avtomaticheskogo upravleniia (polinomi-al'nyi podkhod) [Theory and methods of analytical synthesis of automatic control systems (polynomial approach)]. Moscow, Fizmatlit, 2012. 360 p.
3. Voevoda A. A. Matrichnye peredatochnye funktsii (Osnovnye poniatiia) [Matrix transfer functions (Basic concepts)]. Konspekt lektsii po kursu «Proektirovanie sistem upravleniia» dlia 4-5 kursov AVTF; pod redaktsiei A. S. Vostrikova. Novosibirsk, NGTU, 1994. 94 p.
4. Voevoda A. A., Shoba E. V. O razreshimosti zadachi avtonomizatsii mnogokanal'noi sistemy. Chast' 1 [About diagonally decoupling for multi-output systems. Part 1]. Sbornik nauchnykh trudov Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta, 2010, no. 2 (60), pp. 9-16.
5. Voevoda A. A., Shoba E. V. O razreshimosti zadachi avtonomizatsii mnogokanal'noi sistemy. Chast' 2 [About diagonally decoupling for multi-output systems. Part 2]. Sbornik nauchnykh trudov Novosibirskogo gosu-darstvennogo tekhnicheskogo universiteta, 2010, no. 3 (61), pp. 41-50.
6. Voevoda A. A., Shoba E. V. O razreshimosti zadachi avtonomizatsii mnogokanal'noi sistemy. Chast' 3 [About diagonally decoupling for multi-output systems. Part 3]. Sbornik nauchnykh trudov Novosibirskogo gosu-darstvennogo tekhnicheskogo universiteta, 2010, no. 4 (62), pp. 3-12.
7. Voevoda A. A., Shoba E. V. O razreshimosti zadachi avtonomizatsii mnogokanal'noi sistemy. Part 4 [About diagonally decoupling for multi-output systems. Part 4]. Sbornik nauchnykh trudov Novosibirskogo gosu-darstvennogo tekhnicheskogo universiteta, 2011, no. 3 (65), pp. 11-18.
8. Shoba E. V. Modal'nyi metod sinteza mnogokanal'nykh dinamicheskikh sistem s ispol'zovaniem po-linomial'nogo razlozheniia: dis. ... kand. tekhn. nauk [The modal method for the synthesis of multi-channel dynamic systems using a polynomial expansion: diss. ... cand. tech. sci.]. Novosibirsk, NGTU, 2013. 192 p.
9. Voronoi V. V. Polinomial'nyi metod rascheta mnogokanal'nykh reguliatorov ponizhennogo poriadka: dis.... kand. tekhn. nauk [A polynomial method for calculating the multi-channel controllers low order: diss. ... cand. tech. sci.]. Novosibirsk, NGTU, 2013. 173 p.
10. Shoba E. V. Raschet dvukhkanal'nogo reguliatora dlia stabilizatsii polozheniia kabiny lifta s uchetom zhestkosti kanata i dvumia upravliaiushchimi vozdeistviiami [Calculation of the 2-nd channel regulator for stabilization of provision of the cabin of the lift taking into account rigidity of the rope and two managing directors of influences]. Sbornik nauchnykh trudov Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta, 2012, no. 4 (70), pp. 45-52.
11. Shoba E. V. Raschet mnogokanal'nogo reguliatora dlia podderzhaniia zadannoi temperatury v kamere polimernoi pokraski [Calculation of the multichannel regulator for maintenance of the set temperature in the chamber of polymeric painting]. Sbornik nauchnykh trudov Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta, 2012, no. 4 (70), pp. 53-60.
Voevoda Alexander Aleksandrovich — Russia, 630073, Novosibirsk; Novosibirsk State Technical University; Doctor of Technical Sciences, Professor; Professor of the Department of the Automation; ucit@ucit.ru.
Bobobekov Kurbonmurod Mullomirakovich — Russia, 630087, Novosibirsk; Novosibirsk State Technical University; Postgraduate Student of the Department of Automatics; kurbon_111@mail.ru.
The article submitted to the editors 21.06.2017
INFORMATION ABOUT THE AUTHORS
