CC BY
13Научный вестник НГТУ. - 2013. -№ 2(51) СООБЩЕНИЯ
УДК 681.513
Полиномиальный метод расчета многоканальных регуляторов заданной структуры
А.А. ВОЕВОДА, В.В. ВОРОНОЙ
Предлагается полиномиальная методика синтеза многоканальных регуляторов, составленная путем объединения модальных методик [1-6] и использования принципа «обратного» дифференцирования для вычисления неизвестных параметров регулятора. В основу методики заложено решение так называемого диофантова уравнения. По сравнению с [1, 2] предлагается использовать другую структуру матрицы Сильвестра, более удобную для поиска и оценки линейно-зависимых строк. Для поиска параметров регулятора используется оптимизационная процедура сдвига полюсов замкнутой системы в заданную область устойчивости, для чего предлагается рассматривать движение корней на каждом шаге оптимизационной процедуры на корневом портрете системы.
Ключевые слова: модальный метод синтеза, полиномиальное разложение, многоканальная система, диофан-тово уравнение, регулятор заданного порядка, «обратное» дифференцирование, матрица Сильвестра.
ВВЕДЕНИЕ
В литературе представлено множество различных работ, посвященных синтезу многоканальных регуляторов модальным методом с использованием полиномиального разложения. Так, например, в [1, 2] можно проследить методику синтеза многоканальных систем, учитывающую различные возмущения (доступных и не доступных измерению), измеренное значение управляемой величины и др.
В [3,4] также можно проследить методику, которая при определенных допущениях в постановке задачи синтеза, описании объекта управления и системы в целом представляет собой частный случай [1], положительной чертой которой является возможность оценки нижней границы порядка регулятора и соответствующей желаемой характеристической матрицы замкнутой системы. К сожалению, методику [3, 4] не совсем удобно применять для синтеза систем, когда регулятор имеет определенную заданную структуру1. В связи с чем, при синтезе возникают определенные проблемы, связанные, например, с невозможностью задания точного желаемого распределения полюсов характеристической матрицы системы. Для решения этой проблемы в [5, с. 73-75] предлагается достаточно формализованная методика, представляющая собой набор пошаговых действий, при помощи которых можно достаточно легко синтезировать регулятор заданной структуры при определенных начальных условиях. К сожалению, в самой методике не отображены способы расчета порядка необходимого регулятора, а также выбор порядка полиномов желаемой характеристической матрицы. Однако в самой работе рассматривается выбор степени регулятора путем простого перебора различных вариантов, и далее рассматриваются различные варианты заданных элементов регулятора. В работах [5-9] для синтеза регуляторов пониженного порядка предлагается использовать метод «обратного» дифференцирования, суть которого заключается в применении дифференцирования «справа», т. е. дифференцирование
полинома P(s) — snP{\! s) [10, с. 113-114], но с возвратом к старой переменной s после дифференцирования Р( s) 2 s"-lP '(1 / s) .Обозначим д^1ую операцию гак Р ( s) := dP(s) \ ds .
Статья получена 14 января 2013 г. Повторно после исправлений 28 марта 2013 г. Работа выполнена по заданию Министерства образования и науки по проекту «Исследование предельных точностей оптических методов измерения параметров движения и мехатронных методов управления движением и разработка новых робототехнических и электромеханических систем», Тем-план, проект № 7.559.2011, гос. рег. номер НИР№01201255056
1 Встречаются также понятия: заданный порядок, пониженный порядок, низкий порядок, фиксированная структура, означающих практически одно и то же.
Далее предлагается модальная методика синтеза многоканальных регуляторов заданной структуры с использованием полиномиального разложения, составленная путем объединения всех представленных выше методик, со значительными доработками использования «обратного» дифференцирования, предложенного в [5-9].
МЕТОДИКА СИНТЕЗА
— Исходя из требований к показателям качества поведения системы, задается желаемая область расположения полюсов системы ÇlD, или иначе Q.D - это область D-устойчивости.
— Из описания системы в пространстве состояний осуществляется переход к полиномиальному матричному описанию в виде левого взаимно-простого полиномиального разложения.
— Переход к правому взаимно-простому полиномиальному разложению.
— Определяются столбцовые степени матрицы D(s): |дг- := dDci, г = I, р.
— Выбирается максимальное значение |д := тах(|д1, |д2, ••• Vvp )
— Формируются матричные полиномы D{s) и N(s) в виде:
D(s ) = Dц ^ + + ...D0, N (s) 21ц s^ + +
— Выполняется первый шаг формирования матрицы Сильвестра2:
■No
5 =
D 9
N, N
ц-i
ц-1
D0
N«
— Вычисляется - количество линейно-независимых М-строк в матрице Сильвестра от всех предыдущих.
— Вычисляется сумма ^ у,- = у,
г=1
-v2
р '
— Если сумма линейно-независимых М-строк меньше порядка объекта п, тогда формируется следующая матрица:
5 =
DM Do 0
0 DM Do
NM No 0
0 NM No
— Если нет, то к матрице, полученной на предыдущем шаге, приписываются еще строки М и 0 '1 со сдвигом вправо. Данный пункт выполняется до тех пор, пока сумма линейно-независимых М-строк не станет равной порядку объекта п, тогда построение результанта Сильвестра закончено.
— Проверяется обусловленность полученной матрицы Если сопс1(5) > \|/ (\|/ - некоторое допустимое значение числа обусловленности), то необходимо вернуться к этапу формирования правого взаимно-простого разложения и ввести замену переменной 1 = М ■ 1 повторить формирование матрицы Сильвестра и проверить обусловленность полученной матрицы 5. Если < \|/, то во введении масштаба по времени нет необходимости.
' Отличается от [3] структурой построения, более удобной для поиска линейно независимых М-строк.
216
АЛ. ВОЕВОДА, В В. ВОРОНОЙ
— Вычисляется строчный индекс результанта V := тах(у1, ••• •
— Задаются столбцовые степени матрицы «знаменателя» регулятора: т1 > V — 1.
— Формируются матричные полиномы регулятора:
У(.у) 2 Ут8т + Ут_18т 1 + ... + (8) 2 Хт8т + Хт_18т 1 +... + Х0. Здесь / = 0,/77 числовые матрицы размерностью рх р.
— Выбирается необходимая структура регулятора, задавая некоторые матрицы У, X¡, г — 0,т нулевыми. Иначе говоря, вводится ограничение на структуру регулятора. Обозначим через / - количество заданных нулевых матриц.
— Составляется матрица К = \ Ут Ут_1 ... Хт Хт_х ... Х0].
— Из матрицы Я исключаются блочные нулевые столбцы. Полученную матрицу обозначим К.
— Из матрицы 5 исключаются блочные строки, соответствующие блочным нулевым столбцам матрицы М.
— Выписывается левая часть диофантова уравнения + А'(л') Д'(л').
— Проводится операция обратного дифференцирования I раз:
а1 (у (8) 9(8)+х (8) N (8)) \ а/.
— Выписывается новая матрица Сильвестра. Ниже представлен пример для случая 1=1:
9,-1 29,-2 ... (, -1) 9о 0 0 0
9, 29м ,9о 0 0
0 0 0 (, -1)9, ,9,-1 2,90
21,-2 ... (, -1) Nо 0 0 0
2 N,-1 ,N0 0 0
0 0 0 (, -1) I, 2,10
— Вычисляются столбцовые степени /■, г—\,р формируемой желаемой характеристической матрицы [11] системы Ь по выражению /■ =т1 + .
— Выбирается максимальное значение из всех г — 1, р : /тах г— \Р-
— Из каждого /у ,7 = 1, к равного /тах вычитается /: / ■ = - /, у = 1, к ,
где к - количество столбцов матрицы /* со столбцовыми степенями равными /П1ах . Все остальные столбцовые степени остаются без изменений.
— Выбирается структура желаемой характеристической матрицы (диагональная3 или не диагональная).
3 Диагонализация «знаменателя» облегчает решение задачи автономного управления или задачи
И.Н. Вознесенского.
— Формируется характеристическая матрица системы в виде матричного полинома
— Формируются диагональные матрицы:
Нс (1) = diag(1^1, ^2, ..., ), Нг (1) = diag(1т, , $тр ).
— Вьшисляется малица Ръ = 11т Н~1(я)р(1)Н—(1).
— Если матрица /•}, получилась невырожденной числовой матрицей, то переходим к следующему шагу. Иначе нужно снова перейти выбору столбцовых степеней матрицы Р .
— Задается желаемое расположение корней4 (с учетом масштабирования, если применялось на предыдущих шагах) и рассчитывается каждый полином желаемой характеристической матрицы5.
— Составляется матрица Р — | Рт ц 1 ... -Ро].
— Выписывается система уравнений в матричном виде: ЯБ = Р.
— В матрице 5" необходимо определить нулевые столбцы и удалить их. Обозначим полученную матрицу .
— В матрице Р необходимо также удалить столбцы, соответствующие нулевым столбцам 5". Полученную матрицу обозначим Р\.
— Определяются линейно-зависимые строки матрицы 5*1 и удаляются. Полученную матрицу обозначим .
— В матрице параметров регулятора М необходимо задать нулевыми столбцы соответствующие удаленным строкам на предыдущем шаге. Полученную матрицу обозначим .
— Составляется матричное уравнение: ^ = Р\.
— Решается матричное равнение: ^ = р18—
— Восстанавливается матрица Я.
— Выписываются матрицы регулятора К(л') и Х(х).
— Находится характеристическая матрица и вычисляются ее корни, которые соответствуют полюсам замкнутой системы. Для удобства можно корни отобразить на комплексной плоскости, и рассматривать корневой портрет.
— Если хоть один полюс не попал в желаемую область устойчивости О (л'у 2 О), то необходимо вернуться назад, задать новое значение для желаемого полинома Р и повторить расчеты. Если все полюса попали в желаемую область, то процедура синтеза регулятора заканчивается.
— При необходимости провести моделирование замкнутой системы.
4 В соответствии с теорией графов, разработанной А.В. Чехонадских, удобно задавать желаемые корни кратными.
5 Требование задается не к полной системе, а системе полученной после обратного дифференцирования.
218
А.А. ВОЕВОДА, В В. ВОРОНОЙ
ЗАКЛЮЧЕНИЕ
Предложена методика синтеза многоканальных регуляторов заданной структуры, сочетающая в себе положительные моменты методик, изложенных в [1-5], и принципа «обратного» дифференцирования [5-9]. Сложность построения заключалась в том, что подробная пошаговая методика расчета многоканальных регуляторов была изложена только в [6], во всей остальной литературе методика изложена неявно, а многие моменты требуют доработки и самостоятельного доосмысливания. Кроме того, отличительной особенностью от работ [5-9] является то, что принцип «обратного» дифференцирования заложен в итерационную процедуру с последовательным вхождением в заданную область устойчивости И.
СПИСОК ЛИТЕРАТУРЫ
[1] Гайдук А.Р. Теория и методы аналитического синтеза систем автоматического управления (полиномиальный подход) / А.Р. Гайдук. - М.: ФИЗМАТЛИТ, 2012. - 360 с.
[2] Гайдук А.Р. Теория автоматического управления: учебник / А.Р. Гайдук. - М.: Высш. шк, 2010. - 415 с.
[3] Chen С.Т. Linear System Theory and Design, Third Edition / C.T. Chen // New York Oxford, 1999. - 334 p.
[4] Doyle J.C. Multivariable feedback design: concept for a classical/modern synthesis / J.C. Doyle, G. Stein // IEEE Trans. Aut. Control. - 1981. - AC-26. - № 1. - P. 4-16.
[5] Шоба E.B. Модальный метод синтеза многоканальных динамических систем с использованием полиномиального разложения: дис. ... канд. техн. наук/Е.Б. Шоба. -Новосибирск: НГТУ, 2013. - 192 с.
[6] Вороной В.В. О методике синтеза регулятора «пониженного» порядка методом «обратного» дифференцирования / В.В. Вороной // Теория и практика современной науки: материалы VIII Международной научно-практической конференции, г. Москва, 26-27 декабря 2012 г. В 3 т.: т. 1. Науч.-инф. издат. центр «Институт стратегических исследований». - Москва: Изд-во «Спецкнига», 2012. - С. 215-220.
[7] Voevoda А.А. Low order controllers synthesis using the «reverse derivative» / A.A. Voevoda, V.V. Voronoy, E.V. Shoba // Proceedings of RFBR and DST Sponsored «The 2-nd Russian-Indian Joint Workshop on Computational Intelligence and Modern Heuristics in Automation and Robotics», 10-13 September, 2011. - P. 12-22. [Синтез регулятора пониженного порядка с использованием «обратного дифференцирования»].
[8] Воевода А.А. Модальный синтез многоканальных регуляторов пониженного порядка с использованием «обратного» дифференцирования характеристической матрицы / А.А. Воевода, В.В. Вороной // Сб. науч. тр. НГТУ. -2011. - № 3 (65). - С. 3-10.
[9] Воевода А.А. Модальный синтез многоканального регулятора пониженного порядка с использованием «обратной» производной на примере трёхмассовой системы / А.А. Воевода, В.В. Вороной, Е.Б. Шоба // Науч. вестник НГТУ.-2012,-№ 1 (46).-С. 15-22.
[10] Немировский А.С. Необходимые условия устойчивости полиномов и их использование / А.С. Немировс-кий, А.А. Воевода, В.В. Вороной, Е.Б. Шоба // Науч. вестник НГТУ. - 2012. - № 1 (46). - С. 15-22.
[11] Петрикевич Я.И. Рандомизированные методы стабилизации дискретных линейных систем / Я.И. Петрике-вич // АиТ. - 2008. - № 11.-С. 103-113.
Воевода Александр Апександрович, доктор технических наук, профессор кафедры автоматики Новосибирского государственного технического университета. Основное направление научных исследований: управление многоканальными объектами. Имеет более 200 публикаций. E-mail: [email protected].
Вороной Вадим Владимирович, аспирант кафедры автоматики Новосибирского государственного технического университета. Основное направление научных исследований: синтез многоканальных систем автоматического управления. Имеет 25 публикаций. E-mail: [email protected].
Voevoda A.A., Voronoy V.V.
Polynomial methodfor calculating multi-channel controllers ofa given structure
Proposed the polynomial procedure for multi-channel controllers synthesis, it is composed by modal methods symbiosis [1-6] and using the "reverse" derivative principle to calculate the unknown controller parameters. The methodology laid the solution of so-called Diophantine equation. Compared to [1, 2] proposed to use a different structure of the Sylvester matrix, more convenient to search for and evaluation of linearly-dependent rows. The optimization procedure of closed system poles shift in the given stability region is used to find the controller parameters, which is proposed to consider the roots movement in each optimization procedure step at the system root portrait.
Key words: modal synthesis method, polynomial decomposition, multi-channel system, Diophantine equation, given order controller, «reverse» derivative, Sylvester matrix.
