УНИВЕРСИТЕТ итмо
НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ март-апрель 2022 Том 22 № 2 http://ntv.ifmo.ru/
SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS March-April 2022 Vol. 22 No 2 http://ntv.ifmo.ru/en/
ISSN 2226-1494 (print) ISSN 2500-0373 (online)
ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ. МЕХАНИКИ И йПТИКИ
doi: 10.17586/2226-1494-2022-22-2-239-245 УДК 62-50
Эллипсоидные оценки траекторной чувствительности
многомерных процессов на основе обобщенной проблемы сингулярных чисел
Роман Оморович Оморов1, Таалайбек Абакирович Акунов2Н, Акылбек Орузбаевич Айдралиев3
1>2'3 Институт машиноведения и автоматики Национальной академии наук Кыргызской Республики, Бишкек, 720071, Кыргызская Республика
1 romano_ip@list.ru, https://orcid.org/0000-0003-3555-1323
2 takunov@mail.rus, https://orcid.org/0000-0002-0923-9777
3 aidraliev@outlook.com, https://orcid.org/0000-0003-2889-3202
Аннотация
Рассмотрена задача исследования чувствительности процессов управления к вариации параметров. Для решения задачи использован аппарат траекторной чувствительности, применение которого совместно с методом пространства состояний позволило построить модели чувствительности. На основе моделей определены эллипсоидные оценки функций траекторной чувствительности по состоянию, выходу и ошибке линейных многомерных непрерывных систем в форме мажорант и минорант. Выполнены вычисления с использованием обобщенного сингулярного разложения матриц, составленных из функций траекторной параметрической чувствительности. Полученные эллипсоидные оценки в силу содержательных возможностей обобщенного сингулярного разложения матриц обладают свойством минимальной достаточности. Метод оценки дает возможность с помощью левого сингулярного базиса, соответствующего экстремальным обобщенным сингулярным числам, выделить в пространствах состояния, выхода и ошибки подпространства, которые характеризуются в каждый момент времени наибольшим и наименьшим по норме дополнительным движением. Правый сингулярный базис позволяет в пространстве параметров выделить подпространства, порождающие наибольшее и наименьшее по норме дополнительное движение. Предложенный подход решает проблему «оптимального номинала», то есть проблему выбора номинального значения вектора параметров агрегатов объекта управления, доставляющих многомерному управляемому процессу наименьшее значение эллипсоидных оценок функций траекторной чувствительности, а также осуществить сравнение протекания многомерных управляемых процессов по эллипсоидным оценкам траекторной параметрической чувствительности. Ключевые слова
линейная многомерная система, эллипсоидная оценка, траекторная параметрическая чувствительность, модель чувствительности, обобщенное сингулярное разложение
Ссылка для цитирования: Оморов Р.О., Акунов Т.А., Айдралиев А.О. Эллипсоидные оценки траекторной чувствительности многомерных процессов на основе обобщенной проблемы сингулярных чисел // Научно-технический вестник информационных технологий, механики и оптики. 2022. Т. 22, № 2. С. 239-245. doi: 10.17586/2226-1494-2022-22-2-239-245
Elliposoidal estimates of trajectory sensitivity of multi-dimensional processes based
on generalized singular values problems
Roman O. Omorov1, Taalaibek A. Akunov2«, Akylbek O. Aidraliev3
i'2'3 Institute of Machine Science and Automation of the National Academy of Sciences of the Kyrgyz Republic, Bishkek, 720071, Kyrgyz Republic
1 romano_ip@list.ru, https://orcid.org/0000-0003-3555-1323
2 takunov@mail.ru«, https://orcid.org/0000-0002-0923-9777
3 aidraliev@outlook.com, https://orcid.org/0000-0003-2889-3202
© OMOPOB P.O., AKyHOB T.A., AngpameB A.O., 2022
Abstract
The problem of studying the sensitivity of control processes to parameter variations is considered. To solve the problem, the trajectory sensitivity apparatus was used, the use of which, together with the state-space method, made it possible to construct sensitivity models. Based on the models, ellipsoidal estimates of the trajectory sensitivity functions in terms of the state, output, and error of linear multidimensional continuous systems in the form of majorants and minorants are determined. Calculations are performed using the generalized singular value decomposition of matrices composed of trajectory parametric sensitivity functions. The resulting ellipsoidal estimates, due to the meaningful possibilities of the generalized singular value decomposition of matrices, have the property of minimal sufficiency. The estimation method makes it possible, using the left singular basis corresponding to the extremal generalized singular values, to determine subspaces in the state, output, and error spaces that are characterized at each moment of time by the largest and smallest additional motion in terms of the norm. The right singular basis allows us to determine subspaces in the parameter space that generate the largest and smallest additional motion in the norm. The proposed approach solves the problem of "optimal nominal", that is, the problem of choosing the nominal value of the parameter vector of the plant units that provide the multidimensional controlled process with the smallest value of the ellipsoidal estimates of the trajectory sensitivity functions, as well as to compare the flow of multidimensional controlled processes according to the ellipsoidal estimates of the trajectory parametric sensitivity. Keywords
linear multivariable system, ellipsoidal estimate, trajectory parametric sensitivity, sensitivity model, generalized singular value decomposition
For citation: Omorov R.O., Akunov T.A., Aidraliev A.O. Elliposoidal estimates of trajectory sensitivity of multi-dimensional processes based on generalized singular values problems. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2022, vol. 22, no. 2, pp. 239-245 (in Russian). doi: 10.17586/2226-1494-2022-22-2-239-245
Введение
Аппарат траекторной чувствительности [1-3] широко используется в фазе становления теории чувствительности динамических систем как самостоятельной ветви общей теории систем управления. В сочетании с техникой метода пространства состояний аппарат строит модели чувствительности, выполняющие функции наблюдателя нормированного дополнительного движения динамической системы, порожденного вариациями параметров последней относительно их номинальных значений. Аппарат стал идейной основой метода сравнительной чувствительности Круца-Перкинса [4]. Однако в последнее время работы по теории параметрической чувствительности управляемых процессов, использующие методологию траекторной чувствительности, появляются, за некоторым исключением [5], достаточно редко в отечественных и зарубежных научных работах. Тем не менее процесс алгебраиза-ции современной теории управления обнаружил новые возможности аппарата траекторной чувствительности, позволившие сконструировать скалярные оценки дополнительного движения многомерных управляемых процессов, названных эллипсоидными [6].
Постановка задачи
Основная задача исследования чувствительности [1, 2] — изучение свойств дополнительных движений и определение их связи со свойствами исходной системы. Под дополнительным движением Дф(/, Дq) при вариации параметров q относительно номинального значения q0 при условии, что q = qo + Дq, понимается вариация вида:
Дф(/, qo, Дя) = Дф(/, Дц) = ф(/, qo + Дя) - ф(0,
где ф соответствует векторам: состояния х, выхода у и ошибки £. Согласно классической теории чувствительности [1, 2] первым приближением для дополнитель-
ного движения является выражение Дф(/, Д я) = 0(ОДц, где А — приращение, t — имеет смысл времени, я — вектор параметров, 0(0 — матрица, составленная из соответствующих функций чувствительности.
В силу многомерности задача исследования параметрической чувствительности векторных процессов достаточно сложна и связана с рассмотрением множества вариантов дополнительных движений.
Отметим, что основная задача — оценка параметрической чувствительности континуума траекторий многомерных систем управления в виде соотношения для векторных норм по состоянию, выходу и ошибке. С помощью соотношения выполнена оценка подпространств максимальной и минимальной чувствительности. С этой целью введены эллипсоидные оценки функций траекторной чувствительности для отношения евклидовых векторных норм 5 в виде скалярных мажорант и минорант Ъ^У) этих функций в форме
^(0 < 5(0 < ^м(0.
Конструирование оценок осуществлено с использованием обобщенного сингулярного разложения [7] соответствующих матриц траекторной чувствительности [8].
Подход позволяет решить проблему выбора номинального значения вектора параметров агрегатов объекта управления, доставляющих многомерному управляемому процессу наименьшее значение эллипсоидных оценок функций траекторной чувствительности, выполнить сравнение протекания многомерных управляемых процессов по эллипсоидным оценкам траекторной параметрической чувствительности.
Схема использования обобщенного сингулярного разложения матриц
Пусть задача исследования многомерных систем управления сводится к векторно-матричному представлению вида
к(т) = П(т)х(т = 0), Vi,
(1)
где П(т) — матрица, связывающая вектор к(т), зависящий от переменной т, с вектором х(т = 0) = х(0), стационарным по т. Тогда могут быть записаны оценочные неравенства вида
ст(т) < ||к(т)||/||х(0)|| < ам(т), Ут,
где ат(т), ам(т) — экстремальные элементы алгебраического спектра са{П(т)} сингулярных чисел матрицы П(т) (1), однозначно определяющие на матрице правых сингулярных векторов У(т) те из них, которые на сфере ||х(0)|| = у = const отображаются в наибольшую и наименьшую полуоси эллипсоида с полуосями уам(т) и уат(т) соответственно, где у — численное значение векторной нормы ||х(0)||.
Пусть задача исследования многомерных управляемых процессов, сведенная к (1), дополнена вектор-но-матричным соотношением
С(т) = ¥(т)Х(т),
где матрица ¥(т) связывает векторы £(т) и х(т) в форме линейной алгебраической задачи.
Пусть при этом в задаче исследования многомерных управляемых процессов представляет интерес изучение поведения отношения
5к?(т) = ||к(т)||Щт)||
евклидовых векторных норм. Тогда в силу обобщенного отношения Релея [9-12], которое является скалярной формой обобщенной задачи собственных значений, оказываются справедливыми оценочные неравенства
WW < §к?(т) = ||к(т)||/||«т)|| < Хптм(т), (2)
где ^п^т(т) и ^пш(т) — скалярные оценки при условии, что ¥г(тЩт), Ут невырожденная матрица, вычисляются для любых т с помощью решения обобщенного характеристического уравнения
det(ц(т)¥г(тЩт) - Пг(т)П(т)) = 0
так, что
Хпт(т) = 1д1/2(т)|.
(3)
Представим матрицы задач П(т) и ¥(т) в форме
П(т) = ф(т), ¥(т) = ^(т),
где п(т) — произвольная квадратная матрица, тогда оценочные неравенства (2) становятся стационарными по т и принимают вид
V™ < §к?(т) = 1|к(т)||/Ш|| < ^QRM, ^т. (4)
Экстремальные оценки ^йт(т), ^щ(т) в формуле (4) вычисляются при условии, что матрица ЯЯ обратима, с помощью решения обобщенного характеристического уравнения, стационарного по т
- =0 при условии, что = |ц1/2|.
Таким образом, задача сведена к обобщенной проблеме сингулярных чисел [7], а оценки вида (2) в силу
свойства симметричности матриц и QTQ могут быть получены с помощью экстремальных элементов алгебраического спектра обобщенных сингулярных чисел матриц Я б Ятхп с т > п и Q б ЯРхп, а^Шт, ^ы 6 оа(К, Q) = {ай1^„, ..., айп^}, удовлетворяющих условиям и^ЯХ = Ёк = diag{аRi, г = 1, п}, ак, > 0, У^Х = = diag{аQi, г = 1, q}, OQj > 0, где q = т1п(р, п).
Конструирование эллипсоидных оценок траекторной чувствительности многомерных процессов при внешних конечномерных воздействиях
Утверждение 1. Пусть при произвольном значении параметра q модель для непрерывного многомерного процесса имеет вид
x(t, q) = F(q)x(t, q) + G(q)g(t); x(0, q) = x(0); y(t, q) = C(q)x(t, q),
(5)
где {Р^), G(q), C(q)} — зависимая от параметра q тройка матриц состояния, входа и выхода, векторы состояния х 6 Яп, выхода у 6 Ят и ошибки е 6 Ят зависят от вектора параметров q 6 Яр с номинальным значением q0 так, что х = х(/, q), у = у(/, q), £ = е(/, q).
Для номинального значения вектора параметров q0 введем обозначение
x(t, q),
q=qo
А А
= x(tl y(X q)q=q0 = y(t\ q\
q=qo
А
= e(t),
где х(/), у(/), £(/) — номинальные траектории по состоянию, выходу и ошибке многомерного процесса.
Пусть вектор параметров q в модели (5) претерпевает вариацию в форме q = q0 + Д q. Тогда такая вариация Дq параметров порождает в (5) дополнительные движения по состоянию х, выходу у и ошибке £ в формах
Дх(/, q0, Дq) = Дх(/, Дq) = х(/, q0 + Дq) - х(/), Ду(/, qo, Дq) = Ду(/, Дq) = у(/, qo + Дq) - у(/), Д£(/, q0, Д q) = Д£(/, Дq) = £(/, q0 + Дq) - £(/).
Тогда для дополнительных движений Дх(/, Дq) по состоянию Ду(/, Дq) по выходу и ошибке Д£(/, Дq) получим
Дх(/, Дq) = £(/)Дц, Ду(/, Дц) = Н(/)Дц, Д£(/, Дц) = Е(/)Дц.
В нормах можно записать оценочные неравенства ОЕтш(0||Да|| < ||Дх(/, Дя)|| < а£тах(^)||Дц||, V/; (6) аЯт1п(^)||Дц|| < ||Ду(/, ДЯ)|| < аЯтах(/)||Дц||, V/; (7)
аетт(/)||ДЦ|| < ||Д£(/, Дя)|| < ЦетахОНД^, V/,
где а£т1д(/), аЯт1п(/), а£т1п(/) — минимальные элементы алгебраических спектров сингулярных чисел матриц £(/), Н(/), Е(/) для любых /; а^ах®, антах(0, аЯтах(/) — максимальные элементы этих же спектров.
Запишем матрицы траекторной чувствительности соответственно по состоянию х, выходу у и ошибке £
£(t) = row^/t), j = 1p}, H(t) = row{vj(t), j = TTp}, E(t) = rowj), j = TTp},
(8) (9) (10)
где g., Tj, ej — функции траекторной чувствительности многомерного процесса по состоянию x, выходу y и ошибке £ к параметру qj. □
Доказательство утверждения приведено в работах [8, 13-15].
Очевидно, задача получает полное решение при наличии матриц траекторной чувствительности £(t) по состоянию x, H(t) по выходу y и E(t) по ошибке £. Конструирование матриц осуществим с помощью модели траекторной чувствительности [2, 3, 5, 8, 13], в основу которых положены следующие утверждения [5, 8, 14].
Утверждение 2. Пусть многомерный управляемый процесс с вектором x состояния и вектором y выхода при произвольном значении вектора q параметров имеет векторно-матричное описание (1). Тогда, используя модели чувствительности процесса (1), можно построить реализации матриц (8)-(10), столбцы
j) =
д ôx(t, q)
dqj
; j) =
Д Sy(t, q)
; e,(t) =
д дф, q)
q=qo
Ôqj
q=qo
q=qo ' Sqj j = 1, p которых определяются с помощью векторно-ма-тричных соотношений
X(t) = Fx(t) + Gg(t); x(0) = [0T, x(0)T]T;
x(t) = Cxx(t); y(t) = CyX(t); G(t) = Cax(t); (11) n(t) = C^X(t); e(t) = CeX(t),
где x = [gt xT]T
T =
[Gt gT ... gt]t; y = [nT yT]T
T =
F =
[ПТ nT . npT]T, £ = [eT £T]T, eT = [e\ eT2 ... eT]T;
; Fq = col {F qj; j = 1, p};
I ® F px p F q ; G = G
0 F Gq
Gq = col{Gqj; j = 1, p}; Cq = col{Cq,.; j = 1, p}; Cx =
[O„
Cy [Omxpn C]; Cg [In
nxpn xnxnJ; Cy [Omxpn C]; cg [Inpxp °pnxn];
Cn = = [Ipxp ®C CqJ; Ce = [-Ipxp ®C 0], ® — символ
_ д d(o)(g)
кронекеровского произведения матриц; (°)q = ;
j OQi\
д Ъ'Ч=Чо
(°) = {(°)(q)}| q=q0; (°) — имеет смысл матриц F, G и C, а также векторов x, y и £. □
Доказательство утверждения приведено в [8, 14]. Замечание 1. Осуществленная скаляризация задачи анализа чувствительности многомерных управляемых процессов с помощью эллипсоидных оценок траек-торной параметрической чувствительности, задающих мажоранту и миноранту дополнительного движения многомерного процесса, порожденного конечной вариацией Дq вектора параметров, с одновременным использованием содержательных возможностей сингулярного разложения матриц траекторной чувствительности позволяет:
1) рассматривать сингулярные числа aSmin(t), aHmin(t),
«fining и ^max^-b ^max^X «ЯшЬ^Х aEmin(t) матриц
£(t), H(t) и E(t) в качестве элементов функциональ-Г^Т с у-ичной нормой элемента
ного пространства LT
a(t), при этом наиболее употребительными являются у-ичные нормы при у = 2 и у ^ œ||a(t)||y^œ = sup a(t)
t<=T
так, что последняя из них задает максимальный размер трубки, внутри которой располагается дополнительное движение по состоянию и выходу;
2) осуществить сравнение реализаций регуляторов, образующих с исходным объектом управления многомерный управляемый процесс, по предложенным эллипсоидным оценкам;
3) осуществить решение проблемы «оптимального но-
минала» — сделать выбор номинального значения q0 вектора параметров агрегатов объекта управления, доставляющих многомерному управляемому процессу наименьшее значение эллипсоидных оценок функций траекторной чувствительности;
4) с использованием правого сингулярного базиса V в пространстве Rp параметров выделить подпространство Z{Vmax}, порождающее наибольшее по норме дополнительное движение и подпространство Z{Vm^}, порождающее наименьшее по норме дополнительное движение Ax(t) по состоянию, Ay(t) по выходу и As(t) по ошибке;
5) с использованием левого сингулярного базиса U в пространствах Rn и Rm состояния и выхода выделить подпространство Z{Umax} и Z{Umin}, характеризующиеся в каждый момент времени t наибольшим и наименьшим по норме дополнительным движением;
6) осуществить сравнение протекания многомерных управляемых процессов по эллипсоидным оценкам траекторной параметрической чувствительности.
Конструирование эллипсоидных оценок
траекторной чувствительности многомерных процессов на основе обобщенной проблемы сингулярных чисел
Рассмотрим особый случай, когда предметом исследования являются относительные дополнительные движения 5Х(0 и 5/0 по состоянию х(0 и выходу у(0, порожденные вариациями параметров и определенные соотношениями
5Х(0 = ||Дх(^ Дя)|И|х(0||, Vt, (12)
5/0 = ||Ду& Дя)|И1У«||, Vt. (13)
При решении поставленной задачи обратим особое внимание на контроль за рангом матриц. Выделим две ситуации применительно к рассматриваемым процессам в системе (1). В первой ситуации задача решается в автономной версии, во второй — в предположении, что внешнее воздействие g(t) генерируется конечномерной автономной непрерывной системой.
Применительно к первой ситуации справедливы положения следующего утверждения.
Утверждение 3. Из оценочных неравенств
а2-™« < (1|х(0||2 + ||О(0|р)/||х(0|р < а2~ххМ(Ъ Vt,
где 0-^(0, о—МО — экстремальные для любых t элементы алгебраического спектра {о—.¿(0; I = 1, п} обобщенных сингулярных чисел пары матриц
а~хх(/) 6 Са{ехр^)[1 О]Т, ехр(ГЧ)}, V/,
здесь I - (п х п)-единичная матрица; О - (п х рп) — нуль-матрица, из которой следует, что для относительных дополнительных движений 5Х(/) (12) и 5у(/) (13) оказываются справедливыми оценки
5Х(/) < (а2ххы(/) - 1)1/2||Дц||, V/;
5у(/) < аы{Cn}аXxм(t)«m1{c}||Дq||, ^
аы{(°)}, ат{(°)} — соответственно максимальные и минимальные сингулярные числа матрицы (°)л
Доказательство утверждения приведено в [8, 13, 14]. Применительно ко второй ситуации справедливы положения следующего утверждения.
Утверждение 4. Рассмотрим автономные многомерные системы с векторами состояния
х = [хТ zT]T,
образованные агрегированием соответственно номинальной версии системы (1), а также составной системы (11) с источником внешнего воздействия
¿(0 = Ez(/); z(0) = Zo; %(/) = Pz(t) ,
где х 6 Яп; g, у 6 Ят; z 6 Я1; Р 6 Япхп, Е 6 Яы, G, СТ 6 Япхт, Р 6 Ятх1, £(/) = g(t) - у(/) — вектор ошибки, £ 6 Ят, так, что их векторно-матричные описания
х(/) = FX(/); х(0) = [О„Тх1 Ilхl]Tz(0);
X(t) = FX(t); x(0) = [OipXl
характеризуются матрицами F11 F12
O„x/ I/x/]rz(0);
F =
F21 F22
■ [O GP ]; F ■
F11 F12 F21 F22 J
rF GP1 Lo E j
Тогда для относительных дополнительных движений 5Х(/) (12) и 5у(/) (13) оказываются справедливыми оценки
5Х(/) < (а2йы(/) - 1)1/2а2х , Л)||Дц||, V/; (14)
5V(0 < aM{Cn}(a\u(i) - 1)1/2a7xu(t)0xrrU(t)
xzM
"zxMy
xxM
a-1m{Ac}||Aq||, Vt,
(15)
где аХхМ(/), аххМ(/), ОкМ^ ^хМ(/) — максимальные для любых / элементы алгебраических спектров обобщенных сингулярных чисел соответствующих пар матриц:
01
axx(t) б Oajexp(Ft) aXx(t) e Oajexp(Ft) aXz(t) e Ca exp(Ft)
0
LIJ 0
I
0 0 I
exp(Ft)
, exp(Ft)i2 , exp(Ft)22 h
a zX(t) e Ca (exp(Ft))22, exp(Ft)
Доказательство утверждения приведено в [8, 13].
Замечание 2. В связи с приведенными Утверждениями 3, 4 отметим, что мажорирующие оценки 5X(t) и 5y(t) относительных дополнительных движений 5X(t) (12) и 5y(t) (13), определяемые правыми частями уравнений (6), (14), (7) и (15), позволяют для всех по норме ||Aq|| вариаций Aq вектора параметров и нормам ||x(t)|| и ||y(t)|| номинальных траекторий оценить нормы ||Ax(t, Aq)|| и ||Ay(t, Aq)|| дополнительных движений. Однако в общем случае эти оценки могут обладать большей достаточностью, чем оценки, полученные на спектрах обобщенных сингулярных чисел матриц тра-екторной чувствительности £(t) и H(t), вычисляемых с помощью GSVD-процедуры в программной среде MATLAB [16].
Пример. Рассмотрим двухканальную фотоэлектрическую следящую систему (ФЭСС), имеющую матричное описание
F(q) =
0
1
0
0
-149,0(1 + q2) -9,46 -149,0(1 + q2)qx 0
0
. 2,68(1 + q2)q! 0 0
149 - 0(1 + q2) 0
L -2,68(1 + q2)qi
G(q) =
0 1 -2,68(1 + q2) -1,55 J 0
149,0(1 + q2q 0
2,68(1 + q2)
C(q)=с=[0 0? 01
где вектор параметров q = q q2]T имеет номинальную реализацию q0 = [0 0]T и характеризуется вариациями Aq;, принимающими значения |AqJ < 0,3 (i = 1, 2). При условии, что параметр q! характеризует прямые перекрестные антисимметричные связи между каналами ФЭСС, а q2 — крутизну пеленгационной характеристики в общем тракте системы. Сепаратные каналы номинальной реализации ФЭСС существенно неидентичны и характеризуются показателями колебательности Mmax1 = 1,4 и Mmax2 = 1,2, а также эффективными полосами пропускания на уровне 0,707 Ао^ = 16,84 с-1 и Ао2 = 2,11 с-1 соответственно. Построим эллипсоидные оценки траекторной чувствительности по состоянию и выходу при ненулевом начальном состоянии и ступенчатом конечномерном внешнем воздействии.
На рис. 1 приведены графики процессов по норме ||x(t)||, ||y(t)||, их минорант aSinin(t), a^minW и мажорант a^max(t), aHSmax(t) векторных функций траекторной чувствительности по состоянию x(t) и по выходу y(t) при ненулевом начальном состоянии и ступенчатом внешнем воздействия g(t) = const.
Если оценки функций траекторной чувствительности в виде минорант и мажорант умножить на ||Aq||, тогда получим оценки дополнительных движений ||Ах(к, Aq)|| и ||Ay(K, Aq)|| в виде трубок, центрированных относительно номинальных траекторий по состоянию и выходу. Для примера на рис. 2 приведены центрированные относительно номинальных состояний
X
X
НОН
tt£max(i)
azmin (t)
lly(OII
Olmax(0
Рис. 1. Графики процессов: нормы ||x(t)|| (кривая 1),
мажоранты a2max(t) (кривая 2) и миноранты a2mm(t) (кривая 3) векторных функций траекторной чувствительности по состоянию x(t) (а); нормы ||y(t)|| по выходу y(t) (кривая 1) и мажоранте аЖтах(Г) (кривая 2) (b) при ненулевом начальном состоянии и ступенчатом внешнем воздействия g(t) = const Fig. 1. Graphs of processes: norm ||x(t)|| (curve 1), majorants a2max(t) (curve 2) and minorants a2mm(t) (curve 3) of vector functions of trajectory sensitivity with respect to state x(t) (a); norm of output y(t) (curve 1) and majorant aH2max(t) (curve 2) (b) for a non-zero initial state and a stepped exogenous signal g(t) = const
Рис. 3. Кривые относительных дополнительных движений
Sx(t) по состоянию x(t) (a) и Sy(t) по выходу y(t) (b) Fig. 3. Curves of relative additional motions Sx(t) by state x(t) (a) and Sy(t) by output y(t) (b)
и выходов трубки, содержащие все дополнительные движения, порожденные вариациями вектора Aq параметров, характеризующимися нормой ||Aq|| < 0,3. На рис. 3 показаны кривые относительных дополнительных движений 5x(t) и 5/t) по состоянию x(t) и выходу y(t), порожденные вариациями параметров.
Рис. 2. Графики центрированных относительно номинальных состояний: x(t) (кривая 1) (а) и выходов y(t) (кривая 1) (b) трубок, ограниченных кривыми 2 и 3, содержащие все дополнительные движения, порожденные вариациями вектора Aq параметров, характеризующимися
нормой ||Aq|| < 0,3 Fig. 2. Graphs of tubes centered relative to the nominal states
x(t) (curve 1 (a)) and outputs y(t) (curve 1 (b)) bounded by curves 2 and 3, containing all additional movements generated by variations of the parameter vector Aq characterized by the norm ||Aq|| < 0,3
Заключение
Решена задача исследования чувствительности процессов управления к вариации параметров с помощью построения эллипсоидных оценок функций траектор-ной чувствительности по состоянию, выходу и ошибке линейных многомерных непрерывных систем. При расчете использовано обобщенное сингулярное разложение матриц, составленных из функций траекторной параметрической чувствительности. В зависимости от постановки задачи исследования все множество изучаемых траекторий исследуемого процесса генерируется с помощью задания вектора начальных состояний x(0) с нулевой составляющей g(0) последнего.
Литература
1. Розенвассер Е.Н., Юсупов Р.М. Чувствительность систем управления. М.: Наука, 1981. 464 с.
2. Томович Р., Вукобратович М. Общая теория чувствительности: пер. с англ. М.: Советское радио, 1972. 240 с.
3. Оморов Р.О. Модальная чувствительность, робастность и грубость динамических систем (обзорная статья) // Научно-технический вестник информационных технологий, механики и оптики. 2021. Т. 21. № 2. С. 179-190. https://doi.org/10.17586/2226-1494-2021-21-2-179-190
4. Perkins W.R., Cruz I.B., Gonsales R.L. Design of minimum sensitivity systems // IEEE Transactions on Automatic Control. 1968. V. AC-13. N 6. P. 159-167. https://doi.org/10.1109/TAC.1968.1098853
5. Акунов Т. А., Ушаков А.В. Оценки функций траекторной параметрической чувствительности систем управления // Известия вузов. Электромеханика. 1992. № 1. С. 87-92.
6. Акунов Т.А. Анализ и синтез многомерных систем управления с использованием аппарата эллипсоидных оценок качества процессов: диссертация на соискание ученой степени кандидата технических наук / Санкт-Петербургский институт точной механики и оптики. СПб., 1993. 185 с.
7. Хорн Р., Джонсон Ч. Матричный анализ. М.: Мир, 1989. 655 с.
8. Акунов Т.А., Алишеров С., Оморов Р.О., Ушаков А.В. Модальные оценки качества процессов в линейных многомерных системах. Бишкек: Илим, 1991. 59 с.
9. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.: Наука, 1984. 160 с.
10. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1988. 548 с.
11. Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М.: Мир, 1980. 280 с.
12. Голуб Дж., Ван Лоун Ч. Матричные вычисления: пер. с англ. М.: Мир, 1999. 548 с.
13. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А.В. Синтез дискретных регуляторов при помощи ЭВМ. Л.: Машиностроение, 1983. 245 с.
14. Акунов Т.А., Сударчиков С.А., Ушаков А.В. Анализ алгоритмических проблем при исследовании чувствительности дискретных систем // Известия вузов. Приборостроение. 2008. Т. 51. № 7. С. 17-21.
15. Оморов Р.О., Акунов Т.А. Робастность интервальных динамических систем: устойчивость и эллипсоидные оценки качества // Проблемы автоматики и управления. 2021. № 3(42). С. 47-57.
16. Андриевский Б.Р., Фрадков А.Л. Элементы математического моделирования в программных средах MATLAB 5 и Scilab. СПб.: Наука, 2001. 286 с.
References
1. Rozenvasser E.N., Iusupov R.M. Sensitivity of Automatic Control Systems. Moscow, Jenergija Publ., 1981, 464 p. (in Russian)
2. Tomovic R., Vucobratovic M. General Sensitivity Theory. North-Holland, 1972, 258 p.
3. Omorov R.O. The modal sensitivity, robustness and roughness of dynamic systems (review article). Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2021, vol. 21, no. 2, pp. 179-190. (in Russian). https://doi.org/10.17586/2226-1494-2021-21-2-179-190
4. Perkins W.R., Cruz I.B., Gonsales R.L. Design of minimum sensitivity systems. IEEE Transactions on Automatic Control, 1968, vol. AC-13, no. 6, pp. 159-167. https://doi.org/10.1109/TAC.1968.1098853
5. Akunov T.A., Ushakov A.V. Estimates of the trajectory parametric sensitivity functions for control systems. Russian Electromechanics,
1992, no. 1, pp. 87-92. (in Russian)
6. Akunov T.A. Analisys and Synthesis of Multidimensional Systems Using the Tools of Ellipsoidal Assessments of the Processes Quality. Dissertation for the degree of candidate of technical sciences. St. Petersburg, Saint Petersburg Institute of Fine Mechanics and Optics,
1993, 185 c. (in Russian)
7. Horn R.A., Johnson Ch.R. Matrix Analysis. Cambridge University Press, 1985, 561 p.
8. Akunov T.A., Alisherov S., Omorov R.O., Ushakov A.V. Modal Estimations of Quality ofProcesses in Linear Multivarible Systems. Bishkek, Ilim Publ., 1991, 59 p. (in Russian)
9. Voevodin V.V., Kuznetcov Iu.A. Matrices and Calculations. Moscow, Nauka Publ., 1984, 160 p. (in Russian)
10. Gantmacher F. R. The Theory of Matrices. AMS Chelsea Publishing: Reprinted by American Mathematical Society, 2000, 660 p.
11. Forsythe G.E., Malcolm M.A., Moler C.B. Computer Methods for Mathematical Computations. Prentice-Hall, 1977, 259 p.
12. Golub G.H., Van Loan Ch.F. Matrix Computations. JHU Press, 1996, 756 p.
13. Grigorev V.V., Drozdov V.N., Lavrentev V.V., Ushakov A.V. Synthesis of Discrete Controllers Using Computers. Leningrad, Mashinostroenie Publ., 1983, 245 p. (in Russian)
14. Akunov T.A., Sudarchikov S.A., Ushakov A.V. Analysis of algorithmic problems at research of sensitivity of lumped parameter systems. Journal of Instrument Engineering, 2008, vol. 51, no. 7, pp. 17-21. (in Russian)
15. Omorov R.O., Akunov T.A. Robustness of interval dynamic systems: stability and ellipsoid quality estimates. Automation and Control Problems, 2021, no. 3(42), pp. 47-57. (in Russian)
16. Andrievskii B.R., Fradkov A.L. Elements of Mathematical Modeling in MATLAB 5 and Scilab Computing Environments. St. Petersburg, Nauka Publ., 2001, 286 p. (in Russian)
Авторы
Оморов Роман Оморович — доктор технических наук, профессор, член-корреспондент Национальной академии наук Кыргызской Республики, заведующий лабораторией «Синергетики и хаоса динамических систем», Институт машиноведения и автоматики Национальной академии наук Кыргызской Республики, Бишкек, 720071, Кыргызская Республика, 6602708366, https://orc1d.org/0000-0003-3555-1323, romano_1p@l1st.ru
Акунов Таалайбек Абакирович — кандидат технических наук, главный научный сотрудник, Институт машиноведения и автоматики Национальной академии наук Кыргызской Республики, Бишкек, 720071, Кыргызская Республика,£3 6508211498, https://orc1d.org/0000-0002-0923-9777, takunov@ma1l.ru
Айдралиев Акылбек Орузбаевич — инженер, Институт машиноведения и автоматики Национальной академии наук Кыргызской Республики, Бишкек, 720071, Кыргызская Республика, https://orc1d. org/0000-0003-2889-3202, a1dral1ev@outlook.com
Authors
Roman O. Omorov — D. Sc. (Engineering), Professor, Associate Member of the National Academy of Sciences of the Kyrgyz Republic, Head of Laboratory of Synergetics and Chaos of Dynamic Systems, Institute of Machine Science and Automation of the National Academy of Sciences of the Kyrgyz Republic, Bishkek, 720071, Kyrgyz Republic, 6602708366, https://orcid.org/0000-0003-3555-1323, romano_ip@list.ru
Taalaibek A. Akunov — PhD, Chief Researcher, Institute of Machine Science and Automation of the National Academy of Sciences of the Kyrgyz Republic, Bishkek, 720071, Kyrgyz Republic, gg 6508211498, https://orcid.org/0000-0002-0923-9777, takunov@mail.ru
Akylbek O. Aidraliev — Engineer, Institute of Machine Science and Automation of the National Academy of Sciences of the Kyrgyz Republic, Bishkek, 720071, Kyrgyz Republic, https://orcid.org/0000-0003-2889-3202, aidraliev@outlook.com
Статья поступила в редакцию 25.01.2022 Одобрена после рецензирования 03.02.2022 Принята к печати 30.03.2022
Received 25.01.2022
Approved after reviewing 03.02.2022
Accepted 30.03.2022
Работа доступна по лицензии Creative Commons «Attribution-NonCommercial»