№ 2 (83)
A, UNI
те)
UNIVERSUM:
технические науки
февраль, 2021 г.
DOI: 10.32743/UniTech.2021.83.2-1.26-30
УСТОЙЧИВОСТЬ МНОГОМЕРНЫХ ДИСКРЕТНЫХ ДИНАМИЧЕСКИХ
СИСТЕМ УПРАВЛЕНИЯ
Сиддиков Исамиддин Хакимович
д-р техн. наук, профессор Ташкентского Государственного Технического Университета,
Республика Узбекистан, г. Ташкент E-mail: isamiddin54@gmail. com,
Измайлова Рената Николаевна
докторант
Ташкентского Государственного Технического Университета,
Республика Узбекистан, г. Ташкент E-mail: verona [email protected]
STABILITY OF MULTIDIMENSIONAL DISCRETE DYNAMIC CONTROL SYSTEMS
Isamiddin Siddikov
Doctor of Technical Sciences, Professor of Tashkent State Technical University,
Uzbekistan, Tashkent
Renata Izmaylova
Doctoral student of Tashkent State Technical University, Uzbekistan, Tashkent
АННОТАЦИЯ
В статье рассмотрены вопросы исследования и получения обобщенных показателей оценки качества процессов систем, вывода прогнозирующих уравнений и определения устойчивости многомерных дискретных систем непосредственно по математическим моделям. Основываясь на свойствах матриц смежности, получены выражения, позволяющие непосредственно осуществлять количественную оценку поведения систем. Показаны способы получения обобщенных характеристик структурных состояний в виде норм взвешенных матриц смежности, и на этой основе даны уравнения, прогнозирующие поведение многомерных систем управления динамическим объектом по однократно вычисляемым нормам матрицы смежности на любое наперед заданное число шагов дискретности.
ABSTRACT
The article considers the issues of research and obtaining generalized indicators for assessing the quality of systems processes, deriving predictive equations and determining the stability of multidimensional discrete systems directly from mathematical models. Based on the properties of the contiguity matrices, expressions are obtained that allows to quantify the behavior of systems. Methods of obtaining generalized characteristics of structural states in the form of norms of weighted matrices of contiguity are shown, and on this basis equations, allows to predict the behavior of multidimensional control systems of dynamic object by the once-calculated norms of the incompatibility matrix for any forward specified number of discreteness steps, are given.
Ключевые слова: норма матрицы, система с динамической структурой, устойчивость, дискретность, оценка поведения.
Keywords: matrix norm, system with dynamic structure, stability, discreteness, behavior assessment.
Введение
Для решения задач исследования структурно -сложных динамических систем управления целесообразно перейти от «координатного» способа их описания к описанию на основе некоторых обобщенных характеристик. За одну из таких обобщенных
характеристик систем предлагается использовать норму матрицы смежности [2, 5]. Матрица смежности совпадает с матрицей перехода соответствующей подсистемы, то есть имеется возможность суждения о динамических свойствах подсистем. Поскольку окончательное суждение о качественных показателях системы можно вынести по оценкам, непосредственно
Библиографическое описание: Сиддиков И.Х., Измайлова Р.Н. Устойчивость многомерных дискретных динамических систем управления // Universum: технические науки : электрон. научн. журн. 2021. 2(83). URL: https://7universum.com/ru/tech/archive/item/11218 (дата обращения: 25.02.2021).
№ 2 (83)
A, UNI
те)
UNIVERSUM:
технические науки
февраль, 2021 г.
связанным с координатами, необходимо дополнительно установить связи между обобщенными показателями в виде норм матриц смежности и конкретными координатами.
Основная часть
Пусть имеется линейная система управления с динамической структурой и заданы последовательные соответствия:
g
, =(Q,XjXj), j = 1,2,...,к.
где
X0 = (x (0),x2 (0),...,xB (0));...; Хк = (X (кТ),x2 (кТ),...,X (кТ))
(1)
(2)
- множества, соответствующие переменным состояния линейной системы с динамической структурой (СДС) п-го порядка; Q - график соответствия.
Поставим каждой паре (х (у - 1Т), хг (]Т)) в соответствие коэффициент ан (Т) - значение весовой функции [2, 3].
Рассмотрим динамику функционирования системы на конечном множестве интервалов ,
где Т - переменный период повторения состояний системы.
Обозначив через ||Д.|| нормы матрицы смежности [1,7], х (0) вектор начальных условий системы, и (]Т) - управление на j-м шаге, получим нормы вектора состояния х^ЛТ) для сепаратного и перекрестного каналов передач в виде:
(nt )<|| än\\ an _J|..| Д|| J X (0 )
A
+ B
ÄN—1 M A2 ■
Bji • i« (0) + | anii • an—i
B
• \u(t) +... +
An\ • BN—1
u (n — 2t ) +
(3)
Для случая, когда период повторения - постоянная величина, выражение (3) принимает более компактный вид:
Ii N
(NTAlf X(NT)+Y||All(N—-1)
i=0
bII • U (iT )+| Bll • U (N - 1T)
(4)
Аналогичным образом можно получить выражения для норм векторов состояний многомерной системы. Например, нормы вектора состояний для
двумерной системы определяется следующим образом:
II*1 (nt ) < j|an II AN-i||...||a;| • ||x1 (0)|+1|an|| • || AN—1..|| a2| jb;|| •«(0)+... ...+1\AN | • \BN j • u1 (n—2~) + bn к • u1 (n—iT)} ;
(5)
(nt )<| AN1| -|| AN—,|..| a-|| • Ix21 (0)1+1 A
,||л21|| У 7)21
+ AN J BN—1
U1 (n—2~ ) + | IBN1I • I«1 (N—1T )
AN— ..J A221 •! B-II • U (0) +...
(6)
(nt )<| a2|| AN—11... a12 [\\x2 (ö|+| A
A2
an—1
b12 • U2(0) +...
... + A
BN—1
• u'
(n — 2t )+| BN\\ • u2 (n — 1t )}
(7)
(n~)<|а£|ш-к2!f(0)k||an2 -|ahl...ä№1• u2(0)+...
. II Л2 512 7 ... + AJ •B. _J • U
(n — 2~ ) + || bn2 ii • U2 (n —1t)}
+
+
3
U
91
9
9
9
9
19
№ 2 (83)
UNIVERSUM:
технические науки
февраль, 2021 г.
В этом случае оценки для векторов выходов многомерной системы могут быть получены из соотношения:
y
(NT)
<
NT
i=0 i^r
(9)
Выражения (3) - (9) получены, исходя из известных свойств норм векторов и матриц [1, 2, 7]. Смысл этих соотношений заключается в том, что они дают картину изменения во времени величины нормы вектора состояния или вектора выхода многомерных систем по известным нормам отдельных структурных состояний. Однако, наряду с этим, необходимо иметь возможность оценки самих координат, то есть отдельных компонент векторов выходов и состояний.
Для решения этой задачи в работе предлагается следующий способ. Чтобы оценить вектор состояния (выхода) системы для #-го момента времени, найдя вектор для (#-1)-го момента времени, а затем, пользуясь скалярным произведением векторов, получим выражение:
r (nt )=/ xr (n - lt l h), (10)
xr
с помощью которого определяются значения нужных координат.
В выражении (10) хГ (ит) - компонента вектора состояния г-го канала; хг (и —1~) - вектор состояния г-го канала для (#-1)-го момента времени; И{ - вектор каналов, инцидентных вершинам
хг (N-1~) и вершине X (иТ) е Xг (иТ ). Применяя неравенство Коши [1], получим:
' n
V i=i
\2 / <
У
nn
V i=i yV i=i у
(11)
Извлечение из обеих частей квадратных корней приведет к другой форме неравенства Коши:
a, b
<1 \a\
(12)
(nt )
Сравнением (12) и (10) для компоненты х1 можно получить оценку, связанную с нормами векторов Xг (И—1Т ) Иг:
xi (nt )< xr (n - it )
h
(13)
Подставляя вместо х "(и — 1Т) его значения в соответствии с выражениями (3) - (9), получим искомый результат, то есть оценку координат через обобщенные характеристики подсистем (структурных
состояний ) - нормы матриц смежности А. Например, для случая, когда матрицы смежности на всех интервалах одинаковы, для произвольной х координаты многомерной системы получим:
x
(nt )<|
(N-i-2)
(N-1)
xr (0)1 hJ + 2 (iT )x
i=0
Br (0)1 h 1 + |Mr (N - 2T )|
(14)
Bih
Подобным образом можно получить выражения координат состояния и для других случаев. Выражения, аналогичные (13), позволяют прогнозировать значения переменных состояния на любое наперед заданное число интервалов рассмотрения по известным однократно вычисляемым нормам матриц и векторов, не вычисляя на каждом шаге самих векторов состояния [1, 4, 6]. Этот результат дает возможность осуществлять анализ и синтез структурно-сложных систем, рассматриваемых с позиций СДС, на уровне обобщенных характеристик, что важно в связи с известной проблемой размерности систем.
x(kT )= A kx (о).
(15)
Рассмотрим условие, при котором X(кТ) ^ 0 ,
при к ^ го и X(0) Ф 0. Из выражения (15) видно,
что это имеет место, если Ак ^ 0 при к ^ го . В свою очередь это возможно, если нормы матриц
< 1 для любого из значений I = 1,2,..., к. Если
учесть, что максимальная из сумм модулей передач дуг, инцидентных вершинам из множеств
X (7 — 1Т)} и {хг (]Т)} совпадает с нормами матрицы а вида:
= тах
n
ук, (t )
i =1
A
= тах
n
XkT)
(16)
(17)
7=1
X (кТ) - векторы состояния СДС, то достаточность доказана.
Необходимость условия устойчивости. Предположим, что X(0) Ф 0 и система устойчива (асимптоти-
чески), но
A
> 1. Вектор состояний x(kT)
в этом случае будет неограниченно расти, то есть
>
r
X
l
b
№ 2 (83)
UNIVERSUM:
технические науки
февраль, 2021 г.
х (кТ) ^ го при к ^ го , что противоречит предположению.
Определение устойчивости, согласно сформулированного утверждения, осуществляется по следующему алгоритму:
1. На каждом шаге вычисляются максимальные значения сумм модулей коэффициентов каналов передач:
{Хг (уТ)}, то есть тах X\ан(Т) ;
г 1=\
п
тах X Иг/ (Т ) •
1 1=1
2. Если на некотором шаге один из максимумов сумм модулей коэффициентов каналов передач становится меньше единицы, то система устойчива.
Пример решения.
Требуется определить устойчивость двумерной несинфазной системы с конечной длительностью
замыкания импульсных элементов и непрерывной части:
/(г) = /2(г) = 1(г); Т = Т2 = 1с ; ^ =г2 = 0,2с ; ^(- )=^; *(* )=^;
с'2 (- )=^; с 21 (^ )=цЬ;
¿3 , { 2? ¿4,
Структурная схема приведена на рисунке 1. Смещение импульсов на выходе импульсного элемента (ИЭ2) по отношению к входным импульсам ИЭ1 - 0,3 с.
Имеем СДС с тремя структурными состояниями
51,со ответствуют замкнутым
состояниям ИЭ1 и ИЭ2, Б2 соответствует разомкнутому состоянию ИЭ1 и ИЭ2.
Рисунок 1. Структурная схема двумерной дискретной динамической системы управления.
Цикл структурных состояний состоит из четырех участков 5 = {5, Б2,5,5 }. Получим для интервала [(]Т), (] ± 1Т)] его взвешенную матрицу смежности:
А =
' 0,2466 0
-0,0736 0,1353
0,0300 —0,0552
0 0225 —0 0414
—0,0986 0
—0,0593 0
0,0740 —0,0299
0 0182 0 0273
n
Откуда имеем max V = 0,3452
что пока-
i=1
зывает устойчивость системы.
Заключение
Предложен алгоритм определения устойчивости систем управления структурно-сложных динамических объектов, рассматриваемых с позиций систем с динамической структурой. Показаны способы получения обобщенных характеристик в виде норм взвешенных матриц смежности, получаемых на основе координатного способа их описания, на основе которой получено управление, позволяющее на основе вычисления норма матрицы прогнозировать пове -дение многорежимной системы управления с динамической структурой. Достоверность полученного алгоритма подтверждена решением конкретного примера. Разработанный алгоритм отличается простотой реализации на вычислительной машине. Полученные результаты могут быть использованы в задачах управления динамическими объектами непрерывно-дискретного характера.
№ 2 (83)
UNIVERSUM:
технические науки
февраль, 2021 г.
Список литературы:
1. Афанасьев В.Н., Колмановский В.Б. и др. Математическая теория конструирования систем управления: Учебное пособие для втузов. - М.: Высшая школа, 1989. - 447 с.
2. Кадыров А.А. Динамические графовые модели в системах автоматического и автоматизированного управления. - Ташкент: Фан, 1984. - 240 с.
3. Солодовников В.В., Коньков В.Г. и др. Микропроцессорные автоматические системы регулирования. Основы теории и элементы: Учебное пособие / Под ред. В.В.Солодовникова. - М.: Высшая школа, 1991. -
4. Солодовников В.В., Семенов В.В. и др. Расчет систем управления на ЦВМ: Спектральный и интерполяционный методы / Под ред. В.В.Солодовникова, М.Пешеля - М.: Машиностроение. Berlin, Verlag Technik, 1979. -664 с.
5. Цыпкин Я.З. Теория линейных импульсных систем. - М.: Гос. изд-во физ-мат литературы, 1963. - 968 с.
6. Siddikov I.X., Umurzakova D.M. Fuzzy-logical Control Models of Nonlinear Dynamic Objects. Advances in Science, Technology and Engineering Systems Journal Vol. 5, No. 4, 419-423 (2020). ASTESJ ISSN: 2415-6698.
7. Siddikov I.X., Umurzakova D.M. The research on the dynamics of the three-impulse system of automatic control of water supply to the steam generator when the load changes. First International Conference on Advances in Physical Sciences and Materials Journal of Physics: Conference Series 1706 (2020) 012196.
255 с.