CC BY
20
3
АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ. ЭЛЕМЕНТЫ И УСТРОЙСТВА ТЕХНИЧЕСКИХ
СИСТЕМ
УДК 62.50
АЛГЕБРАИЧЕСКАЯ ОРГАНИЗАЦИЯ УСЛОВИЙ ОБОБЩЕННОЙ СИНХРОНИЗИРУЕМОСТИ МНОГОАГРЕГАТНЫХ ДИНАМИЧЕСКИХ
ОБЪЕКТОВ
Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
Рассматриваются проблемы алгебраической организации условий обобщенной синхронизируемости многоагрегатных динамических объектов, которые формируются на основе геометрической постановки задачи синхронного движения агрегатов как движения в линейной оболочке, натянутой на задающий вектор. Ключевые слова: алгебраическая организация, многоагрегатные объекты, обобщенное синхронное движение, собственный вектор, линейная оболочка.
В современных технологических процессах по обработке материальных и обслуживанию гуманитарных потоков можно выделить задачи, в которых многоагрегатные компоненты технологических ресурсов участвуют в формировании совокупного финального продукта (результата), функционируя как «единое целое». Характерными примерами [1] таких технологических процессов по обработке материальных потоков являются процессы формирования и подачи ленточного материала в листопрокатном производстве, в производстве бумаги и тканей, в организации заготовительных процессов в составе «бесскладовых» технологических производств, процессы согласованного движения автономных больших полноповоротных радиотелескопов в составе радиоинтерферометра типа «Квазар» и т.д. Примерами технологических процессов по обслуживанию гуманитарных потоков являются процессы движения строем подвижных технических средств, управляемых антропокомпонентами-операторами, в виде колонны автомобилей, строя самолетов и вертолетов, и т.п., и автономных антропокомпонентов в виде строя военнослужащих на марше, в забеге и т.п.
Основным требованием к системам управления, встраиваемых в техническую среду отмеченных выше технологических процессов, является обеспечение синхронности функционирования динамических агрегатов в составе многоагрегатного процесса. Причем условие обеспечения синхронного функционирования агрегатов системно не должно быть доминирующим, чтобы не накладывались сильные ограничения на возможность назначения динамики технологического процесса с желаемыми показателями.
Исследования показали, что при геометрической постановке задачи управления многоагрегатными комплексами, доставляющего им синхронное функционирование, она обнаруживает несколько алгебраических формулировок. Так, в первой формулировке в рамках системной парадигмы вырождения сложных динамических систем [2, 3] условие синхронного функционирования агрегатов технологических ресурсов можно трактовать как вырождение многоагрегатной системы, при котором эллипсоид отношения «вход-выход» системы вырождается в отрезок прямой. В рамках второй формулировки задача сводится к алгебраическому требованию обеспечения определенной структуры собственных векторов линейного оператора динамической системы «многомерный вход - многомерный выход» (МВМВ) [2], при которой один из собственных векторов принадлежит линейной оболочке обобщенного синхронного функционирования.
Введение
Сравнивая две возможные алгебраические формулировки задачи синхронного функционирования агрегатов по степени развитости алгоритмического обеспечения процедур синтеза законов управления агрегатами, приходится отдать предпочтение второй формулировке. Она имеет развитое алгоритмическое обеспечение, базирующееся на возможностях обобщенного модального управления [4], позволяющего доставлять матрице состояния многоагрегатной системы требуемые собственные значения и собственные векторы. Результаты, полученные в статье, базируются на этом алгебраическом подходе.
Завершая введение в проблему, следует сказать, что попытки решить задачу обеспечения синхронного функционирования многоагрегатных объектов МВМВ-типа были предприняты в работах [5-8]. Они опирались на возможности алгебраических свойств линейного оператора, представленного передаточной матрицей МВМВ-типа. Использование для этих целей векторно-матричных представлений метода пространства состояний существенно расширило возможности решения поставленной задачи.
Постановка задачи
Постановка задачи синтеза управления многоагрегатным техническим объектом (ТО) МВМВ-типа, доставляющего процессам в ТО по выходу синхронное функционирование, формулируется на основе использования векторно-матричного модельного представления этих процессов в форме
х(г) = Ах(г) + Ви (г), у (г) = Сх(г), (1)
где и, х, у - соответственно векторы управления, состояния и выхода ТО; и, у е Ят,
х е Яп; А, В, С - соответственно матрицы состояния, управления и выхода, согласованные по размерности с переменными объекта.
Сформулируем геометрическую трактовку задачи обеспечения синхронного функционирования агрегатов многоканального технического объекта (1).
Геометрическая постановка задачи синхронного изменения компонентов У1 (I = 1, т) т -мерного вектора выхода у в линейной по времени и покоординатно объекта (1) под действием т -мерного вектора задающего воздействия g( г) имеет вид
у(г) е Ьс {go } для V г; (2)
где Ьс ^о} - линейная оболочка, натянутая на вектор go, порождающий задающий вектор g ( г) , параметризованный временем г и обладающий свойством вида (2)
g(г) е ^ ^о} для Vг. (3)
Требуемая динамика синхронного движения компонентов вектора у( г), заданного в форме (2), будет определяться динамикой изменения задающего воздействия g (г) в линейной оболочке, натянутой на go, задаваемой в форме
g (г) = у(г) go (4)
где у(г) - скалярная функция, задающая динамику желаемого развития процессов в ТО (1), удовлетворяющих условию (2). Введем два определения.
Определение 1. Будем называть задачу синхронизации задачей синхронизации в общепринятом смысле, если все компоненты вектора go равны друг другу. □
В этом случае значения скоростей и ускорений изменения всех компонентов, подчиненных условию (4), оказываются равными друг другу.
Определение 2. Будем называть задачу синхронизации обобщенной задачей синхронизации, если компоненты вектора go произвольны. □
В этом случае значения скоростей и ускорений изменения всех компонентов, подчиненных условию (4), могут отличаться друг от друга.
Для раскрытия в полной мере алгебраических возможностей матричных компонентов векторно-матричного описания объекта (1) и формируемой системы управления, которая строится на агрегировании ТО (1) и формирователя сигнала управления (ФСУ) этим объектом, необходимо направить усилия на анализ алгебраических свойств матричных функций от матриц. Решение задачи синтеза ФСУ будем искать в классе аналитических представлений сигнала управления и (г) в форме
и (г) = Kgg (г) - Кх(г), (5)
где Kg, К - соответственно матрицы прямой связи по задающему воздействию и отрицательной обратной связи по вектору состояния ТО (1). При построении ФСУ в форме (5) использована априорная гипотеза о полной измеримости вектора задающего воздействия g(г) и вектора состояния ТО х(г). Агрегирование модели ТО (1) и алгоритма функционирования ФСУ в форме (5) образует систему управления:
х(г) = ^(г) + Gg (г); х(0), у (г) = Сх(г), (6)
где F = А - ВК, G = ВК§ .
Для анализа алгебраических свойств матричных функций от матриц целесообразно перейти к описаниям полученных выше векторно-матричных представлений с использованием лапласовых образов их переменных. В этой связи применение преобразования Лапласа к (4) дает
g (5) = ) go. (7)
В свою очередь, переход в (6) к лапласовым образам позволяет записать:
у(5) = С(*1 -F)-1 х(0) + С(*1 -^)-^)|х(0)=0=С(*1 -х(0) = х(г)\=о . (8)
Подстановка в (8) условия (7) дает
у (5) = С (51 - F)-1 Gg (5) |, (^,) ^ = С (51 - F)-1 Gу( 5) go, (9)
где у( 5), g (5) - лапласовы образы вектора выхода и вектора задающего входного воздействия; у (г), 5) - скалярная функция соответственно времени г и комплексной переменной 5 , задающая динамику желаемого развития процессов в ТО (1).
Основной результат. Поиск достаточных алгебраических условий обобщенной синхронизируемости многоканальных динамических объектов
Решение задачи конструирования матриц Kg и К в алгоритме (5) формирования
сигнала управления, опирающейся на алгебраические свойства матричных компонентов соотношения (9), использует гипотезу о том, что выбором матриц К^ и К выражение (9) можно свести к представлениям вида
у(5) = С(51 - F)-1GW(5)go = С (51 - ^)-1 BKgW(S)go =С (5)go, (10)
где £( 5) - скалярная функция, имеющая вещественнозначное обратное преобразование Лапласа, определяющая динамику изменения вектора выхода ТО (1) в линейной оболочке Ьс }, тем самым доставляет системе вида (6) обобщенное синхронное изменение компонентов вектора выхода у(г) в форме (3).
Решение этой задачи содержится в следующих утверждениях.
Утверждение 1 (У1). Чтобы y(t) е Lc { g0 } или чтобы выполнялось соотношение
y(s) = Z(s)go, где Z(s) - скалярная функция переменной s, имеющая вещественно-значное обратное преобразование Лапласа, достаточно, чтобы выполнялись условия:
1. вектор BKggo был бы правым собственным вектором матрицы состояния F,
т. е. выполнялось соотношение
F (BKggo) = X,- (BKggo), (11)
где X;- (i = 1, n) - одно из вещественных собственных значений матрицы F ;
2. вектор go был бы правым собственным вектором матрицы CBKg, т.е. выполнялось равенство
(CBKg)go = Xglgo, (12)
где Xgl (l = 1, m) - одно из вещественных собственных значений матрицы CBKg , так что матрица прямой связи Kg ищется из условия
Kg = arg {(CBKg)go = Xglgo}. □ (13)
Доказательство. Рассмотрим векторно-матричное выражение y(s) = C(sI - F)-1Gv(s)go = C(sI - F)-1 BKggo^(s).
По свойству [8] матричной функции f (F) сохранять геометрический спектр исходной матрицы F Fi>i = Xi^, где Xi : det(XI - F) = o, в форме
f (FЪ = f (Хг Ъ , (14)
(11) и (14) позволяют для матричной функции от матрицы f (F) = (sI - F)-1 записать цепочку соотношений
(sI - F )-1 BKggo = (s -Xi )-1 BKggo. (15)
Подстановка (15) во второе представление выражения (1o) дает y(s) = CBKggow(s) (s-Xi )-1. (16)
Если в (16) учесть (12), то получим:
y(s) = XgMs) (s-Xi )-1 go Z(s) = XgMs) (s-Xi --1. 417)
Заметим, что собственные числа Xt выбираются из требований к динамике развития процессов по выходу в подпространстве Lc {go}. Альтернативой полученному результату являются положения следующего утверждения.
Утверждение 2(У2). Чтобы y(t) е Lc {go } или чтобы выполнялось соотношение y(s) = Z(s)go, где Z(s) - скалярная функция переменной s, имеющая вещественно-значное обратное преобразование Лапласа, достаточно, чтобы выполнялись условия:
1. вектор BKggo удовлетворял соотношению (11), а, следовательно, и (15);
2. матрица CBKg была бы единичной, т. е. применительно к матрице K g выполнялось соотношение
Kg = arg {CBKg = I} . ^(18)
Доказательство. На основании подстановки (15) в (Ю) получим выражение (16). Если в (16) учесть (18), то получим
X,, i = 1, n
в силу со-
y(s) = go(s-Xi) V(s) =(s-Xi) 1 V(s)go ^^^
где Z(s) = (s-Xj )"V(s). ■(IQ)
Таким образом, если задачу понимать расширительно, то матрицы C и B при конструирования объекта управления и матрица Kg при синтезе закона управления
должны формироваться из условия
(с, B, Kg ) = arg {CBKg = I}. (20)
Как и в случае утверждения 1, собственные числа X, выбираются из требований динамики развития процессов в подпространстве Lc {go } .
Утверждение 3(У3). Матрица M, приводящая произвольную nxn-квадратную матрицу F простой структуры к диагональной форме Л = diag отношения
Л = M~lFM , (21)
имеет своими столбцами [8] собственные векторы матрицы F . □
Доказательство. Запишем базовое уравнение матричного подобия для рассматриваемого случая,
MЛ = FM, (22)
в столбцовой форме
M [Л1 Л 2 ... Л, ... Лп ] = F [M1 M 2 ... Ыг ... Mn ], (23)
где л, , Mi - i-ые столбцы соответственно матриц Л и M (i = 1, п) . Перейдем теперь от матричного уравнения (22) к п векторно-матричным уравнениям вида
MЛ, = FM,; i = , (24)
где i -ый столбец Л, диагональной матрицы Л имеет представление
Л, =[0 ... 0 X, 0 ... о] . (25)
Легко видеть, что с учетом (25) векторно-матричное уравнение (24) принимает вид XiMi = FMi; i = \п. (26)
Векторно-матричное соотношение (26) является определением собственного вектора матрицы F , откуда следует, что Mi - собственный вектор матрицы F . ■
Положения утверждений У1-У3, дополненные концепцией обобщенного модального управления [9], являются алгоритмической основой для синтеза управления и доставляют системе МВМВ-типа достаточные алгебраические условия обобщенной синхронизируемости многоканальных динамических объектов.
Приведем алгоритмы конструирования матричных компонентов закона управления, присвоив им номера 1 и 2, а также сопроводим их аббревиатурой А1 и А2 соответственно.
Алгоритм 1 (А1) конструирования матричных компонентов закона управления (5), опирающийся на положения утверждения 1.
1. Задать go, определяющий Lc {go}- линейную оболочку, натянутую на вектор go, единичной размерности, которой в процессе функционирования ОУ должен принадлежать вектор выхода y (t) .
2. Сформировать требования к показателям качества в переходном и установившемся режимах процессов, протекающих в подпространстве Lc {go } .
3. Сформировать векторно-матричное представление технического объекта в форме (1).
4. Задать Xi и решить уравнение Kg = arg {CBKg =XglI};
5. Сформировать собственный вектор матрицы F -i = BKggo ;
6. Сформировать модальную модель с матрицами Л, H, образующих наблюдаемую пару, где диагональная матрица Л является носителем желаемого спектра собственных значений матрицы F состояния проектируемой системы, в состав которого входит собственное значение X,, удовлетворяющее требованию
Xi = arg{ (I - A)-i е JmB}, (27)
где JmB - образ матрицы B, т.е. пространство ее столбцов;
7. Записать матричное уравнение Сильвестра MЛ-AM =-BH в декомпозированной форме
[м;м] Л Л- a [ M M ] = - B [я; H ] , (28)
где в силу утверждения 3 M = —, Л = Хг- : F -i = X;--.
8. Выделить из системы (28) матричное уравнение Сильвестра,
(v - A )-i =- BH (29)
и найти его решение относительно матрицы H, которое в случае rangB = m = n = dim x принимает вид H = B_1(A - X;I)— , а в случае rangB = m < n находится из соотношения
H = (вТ B)-1 вТ ( a -x,i )-i.
9. Выделить из системы (28) матричное уравнение Сильвестра
MЛ- AM = -BH , (3o)
в котором (Л, H) - наблюдаемая пара, Л = diag {Xj : j ^ i; j = 1, n}, и решить уравнение
(3o) при известных Л, A, B и H относительно матрицы M .
10. Сформировать матрицы M = [ M :M ] и H = [Й \ H ] , при этом обеспечить
максимальную близость собственных векторов матрицы F - j = Mj (j ^ i) для обеспечения модальной робастности системы к ортогональности с учетом возможной минимизации затрат на управление.
11. Вычислить матрицу K обратных связей с помощью соотношения
K = HM - =[Я:H][М!M]-1. (31)
12. Построить реализационную версию закона управления в форме u(t) = K g(t) - Kx(t) = K g(t) - Kyy(t) - Kxx(t) |kv = Kg =
y g (32)
= Kg (g(t) - y(t)) - Kxx(t) = Kgs(t) - Kxx(t) |Kg =k£ = Kes(t) - Kxx(t)
где Kx = K - KgC.
Алгоритм 2(А2) конструирования матричных компонентов закона управления (5), опирающийся на положения утверждения 2
1. Выполнить пп.1-3 алгоритма 1.
2. Выполнить п. 4 алгоритма 1 при условии выполнения равенства X i = 1;
3. Выполнить пп. 5-11 алгоритма 1.
Заключение
Основным результатом работы являются алгоритмы конструирования матричных компонентов закона управления, позволяющего обеспечить обобщенную синхронность функционирования динамических агрегатов, входящих в состав технологического процесса. Однако следует отметить, что при решении поставленной задачи использован только один собственный вектор матрицы отношения «вход-выход». Свобода назначения остальных (т -1) векторов ожидает своего системного использования. Предстоит
разработать критерии их назначения и адаптировать разработанную алгоритмическую базу под эту задачу.
Литература
1. Мирошник И.В. Согласованное управление многоканальными системами. - Л.: Энергоатомитздат, 1990.
2. Бочков А.Л., Дударенко Н. А., Ушаков А.В. Синтез многомерных функционально вырожденных динамических систем // Известия вузов. Приборостроение. - 2008. -Т. 51. - № 1. - С. 25-29.
3. Дударенко Н.А., Ушаков А.В. Вырождение сложных динамических систем с равно-темповыми структурными компонентами // Научно-технический вестник СПбГУ ИТМО. - 2005. - № 19. - С. 44-51.
4. Ушаков А.В. Обобщенное модальное управление // Известия вузов. Приборостроение. - 2000. - Т. 43. - № 3. - С. 8-16.
5. Лаврентьев В.В., Ушаков А.В. К вопросу о синтезе перекрестных связей, обеспечивающих синхронную работу параллельно включенных агрегатов // Приборы и системы автоматики. Труды ЛИТМО, вып. 85. - Л.: ЛИТМО, 1975.
6. Мирошник И.В., Ушаков А.В. Синтез алгоритма синхронного управления системой квазиоднотипных объектов // Автоматика и телемеханика. - 1977. - № 11.
7. Никифоров В.О., Ушаков А.В. Управление в условиях неопределенности: чувствительность, адаптация, робастность. - СПб: СПб ГИТМО (ТУ), 2002.
8. Гантмахер Ф.Р. Теория матриц. - М.: Наука, 1973.
9. Власенко В.А., Мирошник И.В., Сабинин Ю.А., Ушаков А.В. и др. Система управления лентопротяжным механизмом // Электротехн. пром. Серия «Электропривод». - 1977. - Вып. 5 (58).
Дударенко Наталия Александровна -
Полякова Майя Вячеславовна
Ушаков Анатолий Владимирович
Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат технических наук, доцент, [email protected] Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант, [email protected]
Санкт-Петербургский государственный университет информационных технологий, механики и оптики, доктор технических наук, профессор, [email protected]
