ISSN 1814-1196 Научный вестник НГТУ том 58, № 1, 2015, с. 40-61
http://journals. nstu. ru/vestnik Science Bulletin of the NSTU Vol. 58, No. 1, 2015, pp. 40-61
АВТОМАТИЧЕСКОЕ УПРАВЛЕНИЕ AUTOMATIC CONTROL
И ИДЕНТИФИКАЦИЯ AND IDENTIFICATION
УДК 681.511.2
Локальные экстремумы квадратичного критерия оптимальности системы управления
*
с характеристическим полиномом второй степени
А Н. КОРЮКИН1, А.А. ВОЕВОДА2
1 630090, РФ, г. Новосибирск, пр. академика Коптюга, 4, Институт математики им. Соболева СО РАН, кандидат физико-математических наук, старший научный сотрудник. E-mail: koryukin@ngs.ru
2 630 0 73, РФ, г. Новосибирск, пр. К. Маркса, 20, Новосибирский государственный технический университет, доктор технических наук, профессор. E-mail: ucit@ucit.ru
Пусть система управления одноканальная, причем ее характеристическое уравнение имеет степень два и у него нет вещественных корней, пусть также регулятор имеет полный порядок.
Можно считать, что любой характеристический полином имеет старший коэффициент «1». Каждый допустимый регулятор задает характеристический полином. В данной ситуации эти полиномы задаются парами своих коэффициентов (при степенях 0, 1). Геометрически эти пары можно рассматривать как точки плоскости. Обозначим ее через P. Сопоставляя каждому допустимому регулятору его характеристический полином, получим отображение множества допустимых регуляторов на некоторое множество точек плоскости P (область допустимых точек плоскости P). В данной ситуации фазовое пространство (x, v, где x - регулируемая величина, v -скорость ее изменения) является двумерным и геометрически является плоскостью. В качестве критерия оптимальности рассматривается положительно определенная квадратичная форма
x2 +v2. Выбран также «финальный» момент времени. Значение квадратичной формы в финальный момент времени является функцией коэффициентов характеристического полинома.
С помощью Maple показано, что эта функция не имеет точек, подозрительных на экстремум (т. е. нет точек с нулевым градиентом), причем во всей плоскости P. Это позволяет при поиске регуляторов, имеющих ограничения и минимизирующих выбранный квадратичный критерий оптимальности, ограничиться регуляторами, для которых пара из P лежит на границе области допустимых точек плоскости P (по крайней мере для компактной области).
Ключевые слова: локальный экстремум, автоматическое управление, одноканальные системы, регуляторы полного порядка, задача стабилизации, Maple, матричная экспонента, характеристический полином
DOI: 10.17212/1814-1196-2015-1-40-61
* Статья получена 21 октября 2014 г.
Работа выполнена при финансовой поддержке Минобрнауки РФ по государственному заданию № 2014/138. Тема проекта «Новые структуры, модели и алгоритмы для прорывных методов управления техническими системами на основе наукоемких результатов интеллектуальной деятельности».
ВВЕДЕНИЕ
Рассмотрим произвольную одноканальную систему управления [12-19]. Запишем дифференциальное уравнение движения ее объекта управления:
D(p) ■ x = N(p) ■ u . (1)
Здесь x = x(t) - регулируемая величина; t - время; u = u(t) - управляющая величина; D(p), N(p) - полиномы от p с вещественными коэффициентами. Предполагаем, что полиномы D(p), N(p) взаимно простые.
При модальном управлении u выбирается как решение дифференциального уравнения Y(p) ■ u = X(p) ■ (v - x) . Здесь v = v(t) - задание; X(p), Y(p) -полиномы от p c вещественными коэффициентами.
При v тождественно равном нулю, говорят о задаче стабилизации. Для нее v = 0. Поэтому u удовлетворяет дифференциальному уравнению
Y (p) ■ u = -X (p) ■ x. (2)
Фиксируем степень полинома Y(p): deg(Y(p)) = k для k > 0 ; введем обозначение для степени полинома X(p): deg(X(p)) = m для m > 0 ; на коэффициенты полиномов введем только следующие ограничения: старший коэффициент полинома Y(p) равен единице; полиномы X(p), Y(p) взаимно простые. Фиксированное таким образом множество регуляторов, определяемое числами k, m, назовем полным классом регуляторов.
Исключим из рассмотрения старший коэффициент полинома Y(p) (он не меняется и равен единице). Оставшиеся k + m + 1 коэффициентов будем называть параметрами регулятора. Пронумеруем их и будем рассматривать
геометрически - как точки пространства К'' '" ' . Будем называть его пространством регуляторов.
Не обязательно рассматривать все точки пространства регуляторов. Можно ограничиться каким-то множеством. Назовем его (допустимым) множеством регуляторов и обозначим через Or .
В качестве допустимого множества регуляторов разумно рассматривать
все '" ' , или же какое-то его компактное (т. е. замкнутое и ограниченное) подмножество. Разумно также предполагать, что компактная область регуляторов является замыканием открытого множества. В качестве таких областей регуляторов интересно рассматривать шары, прямоугольники и т. д.
Подействуем на обе части равенства (1) полиномом Y(p); на обе части равенства (2) - полиномом N(p). Получим Y(p) ■ D(p) ■ x = Y(p) ■ N(p) ■ u, N( p) ■ Y (p) ■ u = -N (p) ■X (p) ■ x . Теперь правая часть первого из равенств совпадает с левой частью второго равенства. Значит, Y(p) ■ D(p) ■ x = = -N(p) ■ X(p) ■ x, или (Y(p) ■ D(p) + N(p) ■ X(p)) x = 0 .
Напомним, что полином Y(p) ■ D(p) + N(p) • X(p) называется характеристическим полиномом системы управления. Коэффициенты этого полинома являются полиномами первой степени от параметров регулятора.
Обозначим через n степень характеристического полинома. Обычно его старший коэффициент равен единице. Остальные n коэффициентов пронумерованы по степеням (последовательность n элементов). Геометрически их
можно рассматривать как точки пространства Ми. Это пространство характеристических полиномов.
Выписывая для конкретного регулятора характеристический полином,
тем самым задаем отображение ср: М^ "' ' —> К" пространства параметров регулятора в пространство коэффициентов характеристического полинома. «Пробегая» все допустимые регуляторы, при этом отображении получим область возможных характеристических полиномов О к = (Ог ). Если регулятор полного порядка, то отображение Ф: МА' '" ' М" сюрьективно («накрывает» все пространство Ки).
На пространстве Шп рассмотрим произвольный непрерывный критерий оптимальности /•': Ми —» М. .
Задача 1. Пусть задан объект управления; выбран полный класс регуляторов; задана область и с М" и непрерывная функция Т7: и —> К (критерий оптимальности). Нужно найти минимум (максимум) этой функции в области U (т. е. нужно подобрать характеристический полином так, чтобы функция F достигала минимума (максимума)).
Конечно, для компактной области регуляторов эта задача имеет решение. Можно рассматривать также I / М" и искать локальные минимумы (максимумы). Конечно же, потом (после нахождения «оптимального» характеристического полинома) можно «вернуться» к регуляторам (т. е. найти соответствующий «оптимальный» регулятор).
В линейной теории критерий оптимальности - обычно квадратичная форма. В качестве этого критерия рассмотрим произвольную положительно определенную квадратичную форму Т7: II —> К .
В данной работе поставленная задача решается в следующей наиболее простой ситуации: характеристический полином системы управления имеет степень два и не имеет вещественных корней; допустимое множество регуля-
2
торов «накрывает» всю плоскость К коэффициентов характеристического полинома; регуляторы полного порядка, т. е. за счет подбора регуляторов коэффициенты характеристического полинома можно сделать любыми; в качестве критерия оптимальности рассматривается положительно определенная 2 2
квадратичная форма х + V , где x - регулируемая величина, v - скорость ее изменения; на фазовой плоскости выбирается точка, отличная от начала координат; для этой точки в плоскости коэффициентов характеристического полинома ищутся локальные экстремумы критерия оптимальности.
1. ПОСТАНОВКА ЗАДАЧИ
Рассмотрим задачу стабилизации произвольной одноканальной системы управления, т. е. предполагаем, что задание нулевое. Предполагаем также, что характеристическое уравнение этой системы имеет вторую степень. Уравнение движения системы выглядит так:
х + с1х + с0х = 0. (3)
Здесь х = х(г) - регулируемая величина. Выпишем характеристический полином системы:
р2 + ср + е0 . (4)
Следовательно, нормальная форма уравнения движения системы:
х = у, 1> = -с0х-с1у. (5)
Здесь у = у(г) - скорость изменения регулируемой величины. Запишем сопровождающую матрицу [9, с. 203] характеристического полинома системы и затем нормальную форму уравнения движения системы в матричном виде:
( х Л ( о> 1 1
Х=МХ, где х = , М =
V V 1—С0, —С1)
Хорошо известно, что решение уравнения (6) в матричном виде выглядит так:
X(г) = ехр(Мг) ■ Хо , (7)
т
где Х0 = (Х0, у0) - начальная точка фазовой плоскости (здесь Т - транспонирование); Х0 = х(0), = у(0) - регулируемая величина и скорость ее изменения в начальный момент времени ^ = 0.
Рассмотрим регулятор полного порядка. В этом случае коэффициенты Су, Со за счет параметров регулятора можно сделать любыми.
Рассмотрим положительно определенную квадратичную форму Е(Х) = (X, X). Скобками обозначено скалярное произведение; X=Х{() - решение уравнения (6).
Решение X = Х{() уравнения (6) зависит от коэффициентов Су, Со характеристического полинома. Цель работы - поиск локальных минимумов для формы Е(Х) = (X, X), рассматриваемой как функция коэффициентов характеристического полинома (более общо - поиск точек, подозрительных на экстремум, т. е. точек с нулевым градиентом).
2. РЕГУЛЯТОРЫ, ДЛЯ КОТОРЫХ КРИТЕРИЙ ОПТИМАЛЬНОСТИ ИМЕЕТ НУЛЕВОЙ ГРАДИЕНТ
Из вида матрицы М в равенствах (6) найдем след матрицы М: 8р (М) = —Су . Запишем матрицу М как сумму скалярной матрицы и матрицы с нулевым следом: М = дЕ +А, где д - некоторое вещественное число; Е -единичная матрица; А - матрица с нулевым следом. Тогда 8р (М) = д 8р (Е) + 8р (А), Бр(М) = 2д, —Су = 2д,
д = —Су/2.
(8)
Значит, М =
( 0,1 ^ -с0, 2д
, А = М - дЕ =
( 0,1 ^ 2д,
- дЕ =
( -д , 1
д.
А =
-д, 1 д.
(9)
Теперь
еМг = ехр (дгЕ + Аг) = ехр (дгЕ) ехр (Аг) = едгЕ ехр(Аг) = едг • ехр(Аг),
ехр (Мг) = ехр (дг) ехр(Аг). (10)
Разложим экспоненту от матрицы At в ряд Тейлора:
еАг = Е + Аг +
2!
■ +... =
Е+
А2г2 А4г4 ^ г
2!
+
4!
+
Е+
А2г2 А4г4 ^
3!
+
5!
Аг. (11)
Нетрудно убедиться, что для произвольной матрицы (2 х 2) с нулевым следом ее квадрат является скалярной матрицей. Вычислим квадрат матрицы A:
А2 = (д2 - со)Е.
Чтобы понять, что такое число д - со, вычислим дискриминант D характеристического полинома из формулы (4): О = С1 - 4со = 4(д - со) . Значит,
2 О л2 О Т7 / - сп = —, А = — Е. 0 4 4
(12)
Рассмотрим случай, когда характеристический полином не имеет вещественных корней (корни характеристического полинома образуют комплексную пару), т. е. 0<0. Рассмотрим вещественное число d такое, что d > 0,
= -О / 4. Теперь равенства (12) можно записать так:
с0 = д2 + а2 , А2 =-Л2Е .
(13)
Из равенства (8) следует: q - полусумма корней характеристического полинома. Отсюда и из первого равенства (13) получим: d > 0 - мнимая часть комплексной пары корней характеристического полинома. Используя первое из равенств (13), матрицу A из равенства (9) можно записать через параметры q, d:
А =
-д
1
-(д2 + а2), д
(14)
Из (11) и из второго равенства (13) получим:
р (Аг ) =
ех
\_(й)^+(йг)
2!
4!
* А2г2 А4г4 ^
Е +-+ -
2!
4
Е+
4!
+
А2г2 А4г4 Е +-+ -
3!
5!
Аг =
Е_^+(йг)4 л
3!
5!
Аг = со$(йг) Е + ^п(йг) Аг, йг
ехр(Аг) = й~1(й соъ(йг)Е + ът(йг)А) . (15)
Из равенств (10) и (15) следует
ехр(Мг) = ехр(дг)ехр( Аг) = й_1 ехр(дг) (й со$(йг)Е + $т(йг) А).
Теперь используем равенство (7) и перепишем решение X(t) уравнения (6):
X(г) = ехр(Мг)Х0 = ед ехр(Аг)Х0 = = й — ■ ед (й со$(йг) Е + $т(йг) А)Х0.
Далее, учитывая (14), получим
Х (г) =
й
( й соъ(йг) — $т(йг )д ) Хо + $т(йг )уо
2 2
—$т(йг)(й + д )х0 + (й соъ(йг) + $т(йг)д)у0
(16)
Выпишем теперь скалярное произведение (X, X):
¿Ж
(Х, Х) = -
й
2
((йС — sд) хо + .у )2 + (—^(й2 + д2) хо + (йе + sд)vо)
2
где
С = сов(йг), s = $т(йг).
(17)
Рассмотрим скалярное произведение (X, X) как функцию от параметров управления д, Л. В точках экстремума частные производные этой функции по д, Л равны нулю. Выпишем частную производную от (X, X) по д; умножим ее 2
на ненулевое й и поделим на ненулевое ехр(2д0; синусы и косинусы от Л приведем к двойному углу; соберем слагаемые с хо , У) , с синусами и косинусами от приравняем к нулю:
((((й2 + д2)2 + д2) г + 2д (й2 + д2) + д) (1 — С) + й2 г (1 + С) + —й (2дг +1) 5) х2 +((2гд (й2 + д2 +1) + й2 + 3д2 + 1)(С — 1) + (18) +2й (—(й 2 + д2) г — д + г) 5) уо хо +
+(((й2 — д2)г — д — г)С + й(2дг +1)5 + (й2 + д2)г + д + г)= о,
где
C = cos(2t), S = sin(2dt). (19)
Выпишем теперь частную производную от (X, X) по d; умножим ее на
3
ненулевое d и поделим на ненулевое exp(2qt)); синусы и косинусы от dt приведем к двойному углу; соберем слагаемые с Х0 , Vq , с синусами и косинусами от 2dt; приравняем к нулю:
((-d4 - 2d2qt + q4 + q2 ) C + (((d2 + q2 )2 - d2 + q2 ) t + q) dS + +d4 - q4 - q2 ) Xq +(2((-d2 - q2 +1) d2t - (q2 + 1)q) C +
9 9 9 9 9 (20)
+(-2qt (d2 + q +1) - d + q2 -1) dS + 2q (q2 +1)) vqxq + +((2d2qt + q2 +1) C + ((-d2 + q2 +1) t - q) dS - q2 -1) vq = 0.
В последнем равенстве C, S также означают косинус и синус от 2dt. Замечание 1. cos(2dt) Ф1 для решений уравнения (18) при d ф Q, t ф Q , если дополнительно Xq Ф Q или Vq Ф Q. Поэтому cos(2dt) < 1 (ведь косинус не превосходит единицы).
Действительно, при cos(2dt) = 1 имеем sin(2dt) = 0. Подставим эти значения в равенство (18) и получим 2d2t(vQ + xQ ) = 0 .
3. НУЛЕВОЕ НАЧАЛЬНОЕ ЗНАЧЕНИЕ х0
Рассмотрим сначала случай Х0 = 0 . Выпишем для него формулы (18), (20): 0
(-С +1) гд2 + (2Баг - С +1) д + ((а2 - 1)С + а2 +1) г + Ба = 0,
(Баг+с-1) д2 + (2са2г - Ба) д+(-а3 + а) Бг+с-1 = 0. (21)
Рассмотрим левые части равенств (21) как полиномы от q. Они второй степени. Чтобы оба равенства (21) были выполнены, необходимо, чтобы результант этих двух полиномов от q был равен нулю. Вычислим этот результант.
Вспомним, что С = С08(2аг) и Б = 8т(2аг) . Введем обозначение 2 = ехр(2аг • I) (здесь и далее I - это мнимая единица). Выразим C, S через г:
С = , Б = 2 - 2 _1
2 21
Выразим в результанте С, S через г и разложим получившийся полином на множители. Приравняем к нулю и получим равенство
(2 -1)4 (4а2г22 + 22 - 22 + 1)2 =0 (22)
16 24
Имеем z Ф1 (ведь C ф 1 согласно замечанию 1 и z = C + SI). Значит,
g(z) = 0, где g(z) = z2 + 2(2d2t2 -1) z +1. (23)
Далее, комплексное число, сопряженное к z, совпадает с обратным; числа t, d вещественные. Поэтому
gz) = z-2 + 2 (2d2t2 -1) z-1 +1 (чертой над g(z) обозначено сопряжение). Поделим полином g(z) на сопря-
- 9
женный и разложим на множители. Получим g(z) / g(z) = z , или
g(z) • z_1 = g(z) • z , т. е. g(z) • z-1 = g(z) • z_1 . Последнее равенство означает, что число g (z) • z_1 вещественное.
Сделаем в g(z) • z-1 подстановку z = C + IS, z_1 = C - IS. Получим:
_1 J
g (z) • z = 2(dt) + C -1. Поэтому равенство g(z) = 0 равносильно
2(dt )2 + cos(2dt) -1 = 0. (24)
Так как 2sin(dt )2 = 1 - cos(2dt), то равенство (24) равносильно
2(dt)2 - 2sin(dt)2 = 0,
или ±dt = sin(dt) . (25)
Но уравнение ±y = sin(y) имеет только одно вещественное решение: y = 0. Поэтому равенство (25) выполнено только при dt = 0, но t > 0, d > 0.
Замечание 2. При td ф 0 уравнение (24) не имеет вещественных решений.
Поэтому уравнение g(z) = 0 не имеет решений таких, что | z |= 1. Значит, уравнение (22) не имеет решений таких, что | z |= 1. Поэтому уравнения (21) совместно не решаются (ведь их результант по q не может обратиться в 0). Итак, для совместных решений уравнений (18) и (20) получим Xq ф 0 .
4. РЕЗУЛЬТАНТ
Обозначим через i<j , Fq левые части равенств (18) и (20). Рассмотрим их как функции от Xq , Vq : Fi = Fi (xq,Vq) , где i = 1, 2. Это однородные полиномы второй степени. При этом Xq ф 0 . Значит,
fi(xq, vq) = x02f (1 VqIXQ ), где i = 1, 2. Уравнения (18) и (20) имеют вид Fi(xq, Vq) = 0; поэтому они равносильны Fi(1, V) = 0, где V = Vq/Xq , i = 1, 2. (26)
Левые части в (26) являются полиномами второй степени от V. Обозначим результант этих полиномов через Res. Выпишем его:
Res/2 = (8dЬ4 -8sin(2dt)d713 + (-2sin(4dt)d5 + 4sin(2dt)d5)t --d4 cos(4dt) + 4cos(2dt) d4 - 3d4) q4 +((-\6cos(2dt) d8 + \6d8) t3 +
+(-4d6 cos(4dt) + \6cos(2dt)d6 -\2d6)t)q3 + ((\6dl° + \6d8)t4 -\6sin(2dt)d713 + (-4sin(4dt)d5 + 8sin(2dt)d5)t + +(-2d6 - 2d4 ) cos(4dt) + (8d6 + 8d4 )cos(2dt) - 6d6 - 6d4 ) q2 + +(((-\6d10 -\6d*)cos(2dt) + \6dl° + \6d8)P + +((-4d8 - 4d6)cos(4dt) + (\6d8 + \6d6)cos(2dt) - \2d8 - \2d6) t)q + +(8d12 - \6dl° + 8d8)t4 + (8dn - 8d 1)sin(2dt)t3 + +(-\6cos(2dt)d8 + \6d8) t2 +((2d9 - 2d5)sin(4dt) + (-4d9 + 4d5)sin(2dt)) t + +(-d8 - 2d6 - d4 )cos(4dt) +(4d8 + 8d6 + 4d4 )cos(2dt) - 3d8 - 6d6 - 3d4.
Уравнение Res = 0 является необходимым условием разрешимости равенств (26) (т. е. из уравнений (26) следует уравнение Res = 0). Для лучшей читаемости в правой части равенства (27) произведения синусов и косинусов выражены через синусы и косинусы большего аргумента; слагаемые сгруппированы по степеням q, t, по синусам и косинусам.
Опять воспользуемся обозначением z = exp(2dtI); выразим синусы и косинусы, участвующие в записи Res, через z:
n . -n n -n Z + z z — z cos (2ndt) =---, sin (2ndt) = ———
(эти формулы верны для любого целого n). Выразим в Res синусы и косинусы по этим формулам через z. Полученный элемент обозначим через Res z. Это элемент алгебры Щс/.с/\. Это алгебра полиномов с комплексными коэффициентами от q, d, t, z, расширенная элементом z \ :
R[q,d, t,z] с К[q,d,t,z, z"']cl(^d, t, z),
где R(c[,d,t,z} - поле частных алгебры полиномов F¿ [<г/, б/,/, z].
Ассоциативные коммутативные кольца с единицей и без делителей нуля, в которых выполнена теорема о единственности разложения на простые множители, называются факториальными; любая алгебра полиномов над произвольным полем является факториальной [10, 11]. Поэтому алгебра полиномов R[q,d,t,z] факториальная. Любой отличный от нуля элемент ее поля частных однозначно представляется в виде дроби (с точностью до умножения числителя и знаменателя дроби на отличное от нуля комплексное чис-
ло). Поэтому в алгебре V'\q,d ,l,z,z~^ | выполнена теорема об единственности разложения на простые множители. Разложим элемент Res z алгебры M[g,<i,f,z,z_1] на простые множители:
Res z = -d4 gl g2 g3/z2 , (28)
где
g1 = ((2dt +1) q2 + 4Id2qt - 2td3 + Id2 + 2dt +1) z - Id2 - Iq2 -1, g2 = (-Id2 - Iq2 -1) z + (-2dt +1) q2 + 4Id2qt + 2td3 + Id2 - 2dt +1,
g3 = 1 + z2 + (4d2t2 - 2) z . (29)
Вспомним, что переменные q, d, t принимают вещественные значения; при любых вещественных значениях переменных q, d, t элемент Res также принимает только вещественные значения (это гарантирует равенство (27)).
Рассмотрим сопряжение как автоморфизм алгебры М|с/.с/./.г. z~' | такой,
_ — _ _ -1 _i -
что q = q, d = d , t = t, z = z , z = z. Тогда Res z = Res z =
4_ _ _ / 2 4_ _ _ 2
= -d gj g2 g3/ z =-d gj g2 g3 z . В частности, элементы gi (i = 1, 2, 3) алгебры d,/. z, z"1 ] являются делителями элемента Res z.
Элементы gi можно получить из формул (29) заменой I на -I, z на z 1 . После этого, разлагая на множители, можно убедиться, что
gl/g2 =-z, g3/g3 = z2 . (30)
2_ _i _ _i _i _
^этш; g3 = z g3 , z g3 = zg3 , z g3 = z g3 , gj =-z• g2 , а поэтому
gjg2 =-zg2 g2 . Отсюда и из равенства (28) получим Res z = -d4g1g2g3/z2 = -d4 gl1g2 ^ = -d 4 -zg2I2 Ü . Отсюда справедлива
z z z z
-1
Лемма 1. Res z = d4g2g2g3z 1 ,где g2g2 , g3z 1 вещественные. Введем обозначение: T = 2dt. Тогда
z = cos(r) +1 sin(T), z-1 = cos(T) -1 sin(T) .
Пользуясь этими формулами, легко убедиться, что
g3z-1 = 2(2d 2t2 + cos(T)-1). (31)
Имеем Res = 0. Поэтому Res z = 0. При этом d ф 0 . Поэтому по лемме 1
g2 §2 g3 z-1 = 0 . (32)
Из замечания 2 и равенства (31) получим gзz 1 ф 0 . Отсюда и из (32) следует
Лемма 2. g2 • g2 = 0.
Рассмотрим второе из равенств (29) (формула полинома g2 ). Сделаем в ней замену: г = С + I Б, где С = ео8(Т), Б = 8т(Т). Получим
g2 = (с2 + д2 +1) - 1С) + (-2С+1) д2 +1 (4ё 2дг+с2+1)+2ё?(с2 -1).
Вспомним, что I - это мнимая единица, а числа ё, д, I, Б, С - вещественные; найдем вещественную и мнимую части для g2 :
Ж g2) = (с2 + д2 +1) б+ш (с2 - д2 -1), з(g2) = (с2 + д2+1) (1 - с)+лс V.
Так как g2g2 есть сумма квадратов вещественных и мнимых частей элемента g2 , то теперь из леммы 2 получим: вещественная и мнимая часть элемента g2 равны нулю, т. е.
(ё2 + д2 +1) Б + 2Л(ё2 -д2 -1) = 0, (ё2 + д2 +1)(1 -С) + 4ё2дг = 0 .
В этих двух равенствах сделаем замену: I = Т/(2ё). Получим
(а2 + д2 +1)Б + Т(ё2 -д2 -1) = 0 , (а2 + д2 +1)(1 -С) + 2ёдТ = 0. (33)
5. КВАДРАТИЧНЫЕ УРАВНЕНИЯ ДЛЯ МНИМОЙ И ВЕЩЕСТВЕННОЙ ЧАСТЕЙ КОРНЕЙ ХАРАКТЕРИСТИЧЕСКОГО ПОЛИНОМА
Замечание 3. В уравнениях (33) имеем 1 - С > 0, ё > 0, Т > 0, д < 0 .
Действительно, С = ео8(Т), Т = 2А. При этом 1 - С > 0 согласно замечанию 1. Отсюда и из второго из равенств (33) получим 2ёдТ < 0. При этом ё > 0 (это мнимая часть комплексной пары); Т > 0, ведь Т = 2ё, ё > 0, ^ > 0 (контролируем ситуацию после начального момента времени ^ = 0). Значит, д < 0.
9 9
Далее, Т ф 0, ведь Т > 0. При этом ё + д +1 ф 0. Это позволяет собрать в равенствах (33), с одной стороны, ё, д, а с другой - Т; С = ео8(Т), Б = 8т(Т):
1 + д2 - С2 _ Бт(Т) 2дё _1-ео8(Т)
1 + д2 + ё2 " Т ' _1 + д2 + ё2 " Т "
Фиксируем параметр Т > 0 и рассмотрим равенства (34) как уравнения относительно параметров управления д, ё.
Замечание 4. 8т(Т) < Т при Т > 0.
Действительно, | 8т(Т) | < | Т | при любом вещественном Т; равенство достигается только при Т = 0. Поэтому 8т(Т) < Т при Т > 0.
Относительно параметров д < 0, Л > 0 рассмотрим более общие уравнения:
2 2
1 + д — й 2дй ,
~+д2-2 = ° ,--2-2 = Ь (35)
1 + д + й 1 + д + й
(числа а, Ь произвольные, но фиксированы). Понятно, что уравнения (34) являются частным случаем уравнений (35) при
а = $т(Т) / Т, Ь = (1 — со$(Т)) / Т. (36)
Замечание 5. Если уравнения (35) имеют решения д < 0, Л > 0, то а < 1, Ь > 0.
Действительно, из первого равенства (35) видим, что а < 1; так как д < 0, Л > 0, то из второго из равенств (35) получим Ь > 0. Замечание 6. й/д = —(1 — а)/Ь .
Действительно, из первого равенства (35) видно, что 1 — а = = 2й2/(1 + д2 + й2) . Отсюда и из второго равенства (35) получим (1 — а)/Ь = 2й2/(—2дй) = —й/д .
Лемма 3. а2 + Ъ2 < 1, д = —Ь/лД — а2 — Ь2 .
12, (1 + д2 — й2)2 (1 + д2 + й2)2 — (1 + д2 — й2)2 Действительно, 1 — а2 = 1 — ±—^-= ^—2-\ 4 \-<—.
(1 + д2 + й2)2 (1 + д2 + й2)2
При этом (1 + д2 + й2)2 —(1 + д2 — й2)2 = ((1 + д2') + й2)2 — ((1 + д2') — й2)2 = = 4(1 + д2 ) й2 .
Поэтому 1 — а2 = 4(1 + д2)й2 / (1 + д2 + й2)2 . Значит,
1 — а2 — Ь2 = 4й2/(1 + д2 + й2)2 . (37)
При этом й ф о . Поэтому 1 —
а2 — Ь2 > о , или а + Ь < 1.
Ь2 (2дй )2 2 Из второго равенства (35) и из (37) получим ---— = —-—— = д .
1 — а2 — Ь2 4 й 2
При этом Ь > 0, 1 — а2 — Ь2 > о, д < 0. Поэтому д = —ь/л/1 — а2 — Ь2 . Лемма доказана.
Следствие 1. й = (1 — а)/лД — а2 — Ь2 .
В самом деле, й = —(1 — а) д / Ь согласно замечанию 6. При этом д/Ь = 1Д/1 — а2 — Ь2 согласно лемме 3. Поэтому й = (1 — а)/^1 — а2 — Ъ2 .
6. ДИСКРИМИНАНТ УРАВНЕНИЯ F1 (l, V) = 0
Вернемся к уравнению (26) F (1, V) = 0. Напомним, что F - это левая части равенства (18) (частная производная по q исходной квадратичной формы); F является полиномом от V второй степени. Обозначим его дискриминант через Dis и вычислим его:
2Dis = (-4q4 + (8d2 - 8)q2 - 4d4 - 8d2 - 4)cos (T) + +(q4 + (-6d2 + 2) q2 + d4 - 2d2 + 1)cos(2T) +(8dq3 + (-8d3 + 8d) q)sin(T) + +(-4dq3 + (4d3 - 4d) q) sin (2T) +3q4 + (-2d2 + 6) q2 + 3 + 3d4 + (-8T2 +10) d2
Подставим в Dis вместо q, d их выражение через a, b (лемма 3 и следствие 1); разложим получившееся выражение на множители. Получим
Dis = -^>1_ H,
(1 - a2 - b 2)2
где
H = (4b2 - 4) cos(T ) + (a2 - b2) cos(2T ) - 4 sin(T )ab + +2sin(2T)ab + (2T2 - 1)a2 + (2T2 - 3)b2 - 2T2 + 4.
Подставим в H вместо a, b их выражение через T (равенства (36)). Получим
H = -2(T2 + 2cos(T) - 2)2 / T2 .
Итак, справедлива
Лемма 4. Dis = -2 ('-(T2 + 2cos(T) - 2)2 .
(1 - a2 - b2)2 T2
Вспомним, что T = 2dt; из замечания 2 следует, что уравнение Т2 + 2cos(r ) - 2 = 0 имеет только одно вещественное решение: T = 0, а оно нас не интересует. Отсюда вытекает Лемма 5. Dis < 0.
Следствие 2. Частная производная по q исходной квадратичной формы
(уравнение (18)) при q = -b/V' - a2 - b2 , d = (1 - a)/V' - a2 - b2 , a = sin(T)/T, b = (' - cos(T))/T не может обращаться в ноль.
Следствие 3. В любой момент времени t > 0, для любой исходной точки фазовой плоскости (Хо, Vo )^(0,0), частные производные исходной квадратичной формы по q, d, не могут одновременно обращаться в ноль.
Теорема 1. Рассмотрим произвольную одноканальную систему управления. Предположим, что ее характеристический полином имеет степень два, его корни не являются вещественными, а регулятор имеет полный порядок. В качестве фазового пространства рассмотрим плоскость (x,v), где x - регулируемая величина; v - скорость ее изменения. В качестве критерия опти-
2 2
мальности выбрана положительно определенная квадратичная форма х + V . Тогда значение этой формы в любой отличной от начала координат точке фазовой плоскости, как функция коэффициентов характеристического полинома, не имеет точек, подозрительных на экстремум (т. е. точек с нулевым градиентом).
7. ОДНОМАССОВАЯ СИСТЕМА
Рассмотрим одномассовую систему - материальную точку на пружинке. На пружинке жесткости к закреплена масса т. Материальная точка осуществляет одномерные колебания. При отсутствии управляющего воздействия объект движется согласно следующему закону:
ту+ ку = 0.
Здесь у - отклонение массы от положения равновесия. Введем обозначения: ^ - скорость изменения регулируемой величины; (у,м>)Т - вектор-столбец с
т 2
координатами у, м>. Векторы (у, м?) образуют фазовую плоскость М .
При этом движении в этой фазовой плоскости сохраняется энергия. Обозначим через Е удвоенную энергию:
Е(У) = ку2 + шм2 .
Энергия является положительно определенной квадратичной формой от
т 2
векторов У = (у, м>) из фазовой плоскости М . Рассмотрим энергию как
критерий оптимальности.
По формуле (У\,У2)Е = (Е(7\ + ^2) — Е(У1) — Е(У2))/2 выпишем линеаризацию этой квадратичной формы:
(¥1,¥2)Е = кУ\У2 + . (38)
Здесь У1 = (У1,м>1)Т для / = 1, 2. Линеаризацию будем рассматривать как скалярное произведение. Будем называть его энергетическим.
Фазовая плоскость М вместе с энергетическим скалярным произведением будет евклидовым пространством. Для вычислений в этом пространстве нужно выбрать какой-нибудь ортонормированный базис. Обозначим через \\У\\е длину вектора при энергетическом скалярном произведении.
Выберем сначала произвольный ненулевой вектор, например е = (1, о)Т .
Для него = (е1,е\)Е = к . Значит, ЦеЦ^ = ^к. Нормируем его:
е^ 4к
= 1.
Е
Пользуясь формулой (38), выпишем подпространство е^ векторов, ор-
I Т
тогональных к вектору в\ : е1 =(0,1) • М. Это одномерное подпространство,
порожденное вектором = (0,1) . Вычислим длину этого вектора (при энергетическом скалярном произведении): ^Ц^. = (е2,e2)E = m. Значит,
? = л/т , ||е2/ >/т||Е = 1 •
||е2|| E = Ыm , ||е2/Ыт^д-
т т
Итак, векторы ё?1 = ^Д/^, о) , ^2/л/к = (ОдД/т ) являются ортонорми-
рованным базисом фазовой плоскости с энергетическим скалярным произведением.
Для вектора Y = (у,введем обозначения: х, V - его координаты в ба-
____т"> ___
зисе е1, e2 ; (х, v) E - запись этого вектора в базисе el, e2 . Тогда (у,w)T = у • е1 + w • e2 = у*-[к• + • ^2 = x -ё|+у е . Значит,
х = уу[к , у = • (39)
2 2
Поэтому Е = x2 + У2, ведь (Y,Y)E = ку2 + mv2 = (ул/К) + (ул/т) =
2 2 = X2 + V2 •
Запишем уравнение движения объекта при наличии управляющего воздействия: ту + ку = и . Здесь и - управляющее воздействие. В качестве управляющего воздействия рассмотрим ПД-регулятор: и = (а + Ьр)е , где р - дифференцирование по времени; е - ошибка. Рассмотрим задачу стабилизации: е = -у (при нулевом задании). Тогда и = ~{а + Ьр)у, ту + Ьу + (а + к)у = 0. у + фу + (а + к)у)/т = 0 . Запишем последнее уравнение в нормальной форме:
у = м/,м> = -((а + к)у + Ьм>)1т. (40)
Запишем формулы, обратные к (39): у = х/*,[к , w = ^^л/т . Используя их, перейдем в уравнениях (40) к параметрам х, V:
ах ¡к Жу а + к Ь
— = У,/— , — = —-¡= X--V • (41)
Ж \ т & у кт т
Размерность жесткости к пружины есть масса, деленная на квадрат времени. Поэтому размерность параметра к/т обратна квадрату времени. Значит,
размерность параметра Vк / т обратна времени. Поэтому величина
т = ^к / т (42)
безразмерна. Это безразмерное время (время, отмеряемое в единицах л/т / к ). Исключим из уравнений (41) время t по формуле I = т4т / к . Получим
Жх Жу а + к Ь
— = V , — =--х —-== V .
Ж т Жт к укт
Сравнивая с формулами (6), видим
со = (а + к)/к, с1 = ь/\[кт . (43)
Вычислим теперь дискриминант В характеристического полинома из
2 2
л /„ч^ 2Л Ь ла + к Ь — 4(а + к)ш _ ^
формулы (4): и = с — 4со =--4-=-. Тогда В < 0 равно-
кт к кт
сильно
Ь2 < 4(а + к)т . (44)
Из формулы (8) и из второй формулы (43) получим
д = —Ь(2^кт ) . (45)
Теперь из первых равенств (13), (43) и из равенства (45) следует
,2 2 а + к Ь1 4т(а + к) — Ь1 ^ , Л
а = со — д =---=-. При этом а > 0.
к 4кт 4кт
Поэтому
, 1 й = —. 2\
4т(а + к) — Ь1
кт
(46)
8. ПРИМЕР
Рассмотрим конкретную одномассовую систему. Например, к = 1, т = 1. Чтобы воспользоваться теоремой 1, должно быть выполнено неравенство (44). Сейчас это
Ъ2 < 4(а +1).
(47)
Из последнего равенства следует а > -1. Для наглядности нарисуем параболу
Ь2 = 4(а +1) . Неравенство (47) равносильно тому, что точка (а, Ь) лежит внутри параболы (правее кривой). Выберем какой-нибудь прямоугольник, точки которого удовлетворяют этим условиями. Например,
о < а < 1, —1 < Ь < 1.
(48)
Задача 2. Одномассовая система имеет массу т = 1 и жесткость к = 1; ее начальное отклонение Уо = 1; начальная скорость изменения отклонения = 1. Найти параметры ПД-регулятора а, Ь из прямоугольной области (48), обеспечивающие минимум энергии в момент времени 11 = 1.
Решение. Шаг 1. Подставим т = 1, к = 1 в формулу (39). Получим х = у, V = w. Эти равенства при ^ = 0 даютх0 = Уд , Уд = Wд .Вспомним о начальном состоянии объекта: Уд = 1, = 1. Получим Хо = 1, ^ = 1.
Шаг 2. Подставим т = 1, к = 1 в формулу (42). Получим т = ^ . Отсюда в «финальный» момент времени ^ = 1 получим «финальное» «безразмерное время» т = 1.
Шаг 3. Чтобы в данной задаче применить теорему 1, нужно в формуле (17) для критерия оптимальности (X, X) вместо времени ^ рассматривать безразмерное время т . При этом нужно рассматривать финальное время. В итоге в формулу (17) нужно подставить ^ = 1 (сначала вместо ^ подставить т, вместо т подставить т ; наконец вспомним, что т^ = 1).
Подставим также в формулу (17) Хд = 1, Уд = 1. Получим
Е = ехр(2д)((йс - sq + 5)2 + (-5(й2 + д2) + йс + 5д)2)/й2 , где с = соз(й), 5 = 8т(й). (49)
Здесь Е - критерий оптимальности (удвоенная энергия объекта) в финальный момент времени.
Шаг 4. Подставим т = 1, к = 1 в формулы (45) и (46). Получим
д = -Ь, й = ^^/Що + Ъ-Ь2 . (50)
Шаг 5. Подставим в формулы (49) вместо д, й их выражение через а, Ь:
--л2 / ,-
е"Ь
Е = -
с^^м'а + Т)—^2 + 5Ь + 25 1 +1) + с*£4(а+Г)-Ь2 -5Ь
\
4(а +1) - Ь2
где с = 008 \ 2
^^¡Ца + Ц-Ь2^, 5 = 8т+ . (51)
Шаг 6. Вычислим значение функции Е(а, Ь) на границе а = 0 прямоугольной области (48). Для этого подставим а = 0 в формулы (51). Получим
с = 0081 41 - Ь2/4 |, 5 = 8Ш | V1 - Ь2 /4
Е(0, Ь) = е
-Ь
с44 -Ъ2 + 5(Ь + 2)^ + (-25 + с^4 -Ь2 - 5Ь^ ^4 -Ь2) . (52)
Рассмотрим Е(0, Ь) как функцию параметра Ь. Нарисуем ее на интервале -1 < Ь < 1.
Как видим, на этом интервале функция Е(0,Ь) монотонно убывает и достигает наименьшего значения при Ь = 1 . Вычислим его:
Е(0,1) «1,59.
Шаг 7. Вычислим значение функции Е(а, Ь) на границе а = 1 прямоугольной области (48). Для этого подставим а = 1 в формулы (51). Получим
■78
Е (1, Ь) = е"
с = собI 2~14^Ъ2 1, 5 = бш! 2_1л/8-Ъ2
сл]8-Ь2 + ьЬ + 25^1 + Г-4ь + сл/8-Ъ2 -ьЬ
2
(8 - Ь2).(53)
Е(1, Ь) есть функция параметра Ь. Нарисуем ее на интервале -1 < Ь < 1.
Как видим, на этом интервале функция Е(1, Ь) монотонно убывает и достигает наименьшего значения при Ь = 1 . Вычислим его:
Е(1,1) «1,591.
Шаг 8. Вычислим значение функции Е(а, Ь) на границе Ь = 1 прямоугольной области (48). Для этого подставим Ь = 1 в формулы (51). Получим
с = соб (2 1^4а + 3 ) , 5 = БШ ( 2 1^4а + 3 ), Е(а,1) = е"1 с>/4а + 3 + 35)2 + (-25 (а +1) + о/4а + 3 - 5)2^^4а + 3) .
(54)
Е(а, 1) есть функция параметра а. Нарисуем ее на интервале 0 < а < 1.
Как видим, на этом интервале функция Е(а, 1) достигает минимума не на концах отрезка, а в точке локального минимума. Вычислим эту точку:
Е(а,1) «1,498 при а « 0,45 . (55)
Шаг 9. Вычислим значение функции Е(а, Ь) на границе Ь = -1 прямоугольной области (48). Для этого подставим Ь = -1 в формулы (51). Получим
Ь
с =
cos (2 V4a + 3 j , 5 = sin (2 1л]4а + 3 j,
Е (а, -1) = е [^(с-ч/4а + 3 + 5 )2 + (-2.ф +1) + с^ 4а + 3 + 5 )2 ^ у/(4а + 3) (56)
Е(а, -1) есть функция параметра а. Нарисуем ее на интервале 0 < а < 1.
Как видим, на этом интервале функция Е(а, -1) достигает минимума не на концах отрезка, а в точке локального минимума. Вычислим эту точку: Е(а, — 1) « 2,282 при а « 0,51.
Шаг 10. Сравним наименьшие значения функции Е(а,Ъ) на границах прямоугольной области (48) (предыдущие четыре шага).
Вывод. На границе прямоугольника наименьшее значение достигается при Ъ = 1 (формула (55)). Значит, это наименьшее значение во всей прямоугольной области (по теореме 1).
ЗАКЛЮЧЕНИЕ
Теорема 1 позволяет находить минимумы критерия оптимальности (положительно определенной квадратичной формы) на границе допустимой области коэффициентов характеристического полинома (по крайней мере, при условии компактности этой области).
СПИСОК ЛИТЕРАТУРЫ
1. Математическая теория оптимальных процессов / Л.С. Понтрягин, В.Г. Болтянский, Р.В. Гамкрелидзе, Е.Ф. Мищенко. - М.: Наука, Главная редакция физико-математической литературы, 1976. - 392 с.
2. Корюкин А.Н. Предел устойчивости двухмассовой системы с обобщенным ПИД-регулятором // Научный вестник НГТУ. - 2012. - № 3 (48). - С. 178-184.
3. Корюкин А.Н. Наибольший запас устойчивости для одноканальной двухмассовой системы с обобщенным ПИД-регулятором // Научный вестник НГТУ. - 2012. - № 4 (49). -С. 178-185.
4. Корюкин А.Н. Обобщенный ПИД-регулятор двухмассовой системы с наибольшим запасом устойчивости // Научный вестник НГТУ. - 2013. - № 3 (52). - С. 10-17.
5. Корюкин А.Н., Воевода А.А. Наибольшая степень устойчивости двухмассовой системы для регуляторов пониженного порядка // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2014. - № 1 (31). - С. 229-237.
6. Корюкин А.Н. Степень устойчивости ПИД-регулируемой двухмассовой системы с двукратным вещественным корнем и комплексной парой на одной вертикали // Современные технологии. Системный анализ. Моделирование. -2014. - № 3 (43). - С. 120-126.
7. Корюкин А.Н. Разложение на простые множители невырожденных полиномиальных матриц // Научный вестник НГТУ. - 2014. - № 3 (56). - С. 23-36.
I. Корюкин А.Н., Воевода А.А. Наибольшая степень устойчивости трехмассовой системы с регулятором пониженного порядка // Известия Томского политехнического университета. -2014. - Т. 325, № 5: Информационные технологии. - С. 52-59.
9.МальцевА.И. Основы линейной алгебры. - СПб.: М.: Краснодар: Лань, 2009. - 470 с. 102. Бурбаки Н. Алгебра. Т. 3. Модули, кольца, формы: пер. с фр. - М.: Наука, Главная редакция физико-математической литературы, 1966. - 554 с.
II. Бурбаки Н. Коммутативная алгебра. - М.: Мир, 1971. - 706 с.
12. Anthaklis P.J., Michel A.N. Linear systems. - Boston: Birkhauser, 2006. - 670 p.
13. Wang Q.-G. Decoupling control. - Berlin; Heidelberg: Springer-Verlag, 2003. - 356 p.
14. Vidyasagar M. Control system synthesis: a factorization approach. Pt. 1. - San Rafael, California: Morgan & Claypool, 2011. - 184 p.
15. Vidyasagar M. Control system synthesis: a factorization approach. Pt. 2. - San Rafael, California: Morgan & Claypool, 2011. - 227 p.
16. Kailath T. Linear systems. - Englewood Cliffs, New Jersey: Prentice-Hall, 1980. - 704 p.
17. Wolovich W. Linear multivariable systems. - Berlin; Heidelberg: Springer-Verlag, 1974. -
372 p.
18. Воевода А.А. Матричные передаточные функции (основные понятия): учебное пособие. - Новосибирск: Изд-во НГТУ, 1994. - 94 c.
19. Воевода А.А. Матричные передаточные функции (синтез): конспект лекций. -Новосибирск: Изд-во НГТУ, 1995. - 94 c.
Корюкин Анатолий Николаевич, кандидат физико-математических наук, старший научный сотрудник Института математики им. Соболева СО РАН. Основное направление научных исследований - теория автоматического управления. Имеет около 50 публикаций. E-mail: koryukin@ngs.ru
Воевода Александр Александрович, доктор технических наук, профессор кафедры автоматики Новосибирского государственного технического университета. Основное направление научных исследований - многоканальные системы в теории автоматического управления. Имеет более 300 публикаций. E-mail: ucit@ucit.ru
Local extrema of the quadratic criterion of control system optimality with a characteristic polynomial of degree *
A.N. KORYUKIN1, A.A. VOEVODA2
1 Sobolev Mathematical Institute, 4, Academician Koptyg Prospekt, Novosibirsk, 630090, Russian Federation, PhD (Phys&Math), senior research worker. E-mail: koryukin@ngs.ru
2 Novosibirsk State Technical University, 20, K. Marx Prospekt. Novosibirsk, 630073, Russian Federation, D. Sc. (Eng.), professor. E-mail: ucit@ucit.ru
Let the control system be a single-channel systeml, and its characteristic equation have degree two and no material roots. Let the regulator also have a full order.
It is possible to consider that any characteristic polynomial has a unity senior coefficient. Each admissible regulator sets a characteristic polynomial. In this situation, these polynomials are set by pairs of their coefficients (degrees 0, 1). Geometrically these pairs can be considered as a plane point. Let's designate it by P. Comparing each admissible regulator with its characteristic polynomial, we obtain the mapping of a set of admissible regulators on some set of
* Received 21 October 2014.
Work is executed at financial support of the Minobmauka Russia state job № 2014/138, the theme of the project "New patterns, models and algorithms for breakthrough methods of control of technical systems based on high results intellectuality activity".
points of the plane P (area of admissible points of the plane P). In this situation the phase space (x,v) (where x is an adjustable size, v is a speed of its change) is two-dimensional and geometrically is a plane. A positively defined quadratic form is considered as an optimality criterion. Also, a "final" time point is chosen. The value of a quadratic form at a final time point is a function of characteristic polynomial coefficients.
By means of Maple it is shown that this function has no points, suspicious at the extre-mum point (that is, there are no points with a zero gradient), and on the whole plane P. It makes it possible to limit the search of the regulators having restrictions and minimizing the chosen optimality quadratic criterion to regulators for which a pair from P lies on the border of the admissible points area of the plane P (at least for a compact domain).
Keywords: local extremum, automatic control, single-channel systems, full order, problem of stabilization, Maple, matrix exhibitor, characteristic polynomial
DOI: 10.17212/1814-1196-2015-1-40-61 REFERENCES
1. Pontryagin L.S., Boltyanskii V.G., Gamkrelidze R.V., Mishchenko E.F. Matematicheskaya teoriya optimal'nykh protsessov [The mathematical theory of optimal processes]. Moscow, Nauka, Glavnaya redaktsiya fiziko-matematicheskoi literatury, 1976. 392 p.
2. Koryukin A.N. Predel ustojchivosti dvuhmassovoj sistemy s obobschennym PID-regulyatorom [Stability limit of the two-mass system for the generalized PID-controller]. Nauchnyj vestnik NGTU - Science bulletin of the Novosibirsk state technical university, 2012, no. 3 (48), pp. 178-184.
3. Koryukin A.N. Naibol'shii zapas ustoichivosti dlya odnokanal'noi dvukhmassovoi sistemy s obobshchennym PID-regulyatorom [Maximum stability degree for two-mass SISO system with generalized PID-controller]. Nauchnyj vestnik NGTU - Science bulletin of the Novosibirsk state technical university, 2012, no. 4 (49), pp. 178-185.
4. Koryukin A.N. Obobshchennyi PID-regulyator dvukhmassovoi sistemy s naibol'shim zapa-som ustoichivosti [The generalized PID-regulator of the two-mass system with the greatest stock of stability]. Nauchnyj vestnik NGTU - Science bulletin of the Novosibirsk state technical university, 2013, no. 3 (52), pp. 10-17.
5. Koryukin A.N., Voevoda A.A. Naibol'shaya stepen' ustoichivosti dvukhmassovoi sistemy dlya regulyatorov ponizhennogo poryadka [The greatest degree of the stability of a two-mass system for regulators of the lowered order]. Doklady Tomskogo gosudarstvennogo universiteta sistem uprav-leniya i radioelektroniki - Proceedings of Tomsk state university of control systems and radioelec-tronics, 2014, no. 1 (31), pp. 229-237.
6. Koryukin A.N. Stepen' ustoichivosti PID-reguliruemoi dvukhmassovoi sistemy s dvukratnym veshchestvennym kornem i kompleksnoi paroi na odnoi vertikali [The stability degree of a PID-control two-mass system with the double real root and complex pair on the right vertical both]. Sovremennye tekhnologii. Sistemnyi analiz. Modelirovanie - Modern Technologies. System analysis. Modeling, 2014, no. 3 (43), pp. 120-126.
7. Koryukin A.N. Razlozhenie na prostye mnozhiteli nevyrozhdennykh polinomial'nykh matrits [The decomposition on simple multipliers of nondegenerate polinomialny matrixes]. Nauchnyj vestnik NGTU - Science bulletin of the Novosibirsk state technical university, 2014, no. 3 (56), pp. 23-36.
8. Koryukin A.N., Voevoda A.A. Naibol'shaya stepen' ustoichivosti trekhmassovoi sistemy s regulyatorom ponizhennogo poryadka [Maximum stability degree of a threemass system with a lower-order controller]. Izvestiya Tomskogo politekhnicheskogo universiteta. Informatsionnye tekhnologii -Bulletin of the Tomsk polytechnic university. IT Technologies, 2014, vol. 325, no. 5, pp. 52-59.
9. Mal'tsev A.I. Osnovy lineinoi algebry [Foundations of linear algebra]. St.-Petersburg, Moscow, Krasnodar, Lan' Publ., 2009. 470 p.
10. Bourbaki N. Elements de mathematique. Premiere partie. Les structures fondamentales de L'analyse. Livre II. Algebre [Elements of mathematics. Algebra II]. Paris, Hermann et Cie., 1964 (Russ. ed.: Burbaki N. Algebra. Moduli, kol'tsa, formy. T. 3. Translated from French. Moscow, Nauka Publ., Glavnaya redaktsiya fiziko-matematicheskoi literatury, 1966. 554 p.).
11. Bourbaki N. Elements de mathematique. Livre VII. Algebre commutative [Elements of mathematics. Commutative algebra]. Paris, Hermann, 1961 (Russ. ed.: Burbaki N. Kommutativnaya algebra. Translated from French. Moscow, Mir Publ., 1971. 706 p.).
12. Anthaklis P.J., Michel A.N. Linear systems. Boston, Birkhauser, 2006. 670 p.
13. Wang Q.-G. Decoupling control. Berlin, Heidelberg, Springer-Verlag, 2003. 356 p.
14. Vidyasagar M. Control system synthesis: a factorization approach. Pt. 1. San Rafael, California, Morgan & Claypool, 2011. 184 p.
15. Vidyasagar M. Control system synthesis: a factorization approach. Pt. 2. San Rafael, California, Morgan & Claypool, 2011. 227 p.
16. Kailath T. Linear systems. Englewood Cliffs, New Jork, Prentice-Hall, 1980. 704 p.
17. Wolovich W. Linear multivariable systems. Berlin, Heidelberg, Springer-Verlag, 1974.
372 p.
18. Voevoda A.A. Matrichnye peredatochnye funktsii (osnovnye ponyatiya) [Matrix transfer functions (the main concepts)]. Novosibirsk, NSTU Publ., 1994. 94 p.
19. Voevoda A.A. Matrichnye peredatochnye funktsii (sintez) [Matrix transfer functions (synthesis)]. Novosibirsk, NSTU Publ., 1995. 94 p.
ISSN 1814-1196, http://journals.nstu.ru/vestnik Science Bulletin of the NSTU Vol. 58, No.1, 2015, pp. 40-61