Влияние собственных векторов на траектории дискретных систем с задержкой сигнала управления Текст научной статьи по специальности «Физика»

НАУЧНО-ТЕХНИЧЕСКИМ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ март-апрель 2017 Том 17 № 2 ISSN 2226-1494 http://ntv.i1mo.ru/

SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS March-April 2017 Vol. 17 No 2 ISSN 2226-1494 http://ntv.ifmo.ru/en

УДК 519.714.25, 62-52

ВЛИЯНИЕ СОБСТВЕННЫХ ВЕКТОРОВ НА ТРАЕКТОРИИ ДИСКРЕТНЫХ СИСТЕМ С ЗАДЕРЖКО Й СИГНАЛ А УПРАВЛЕНИЯ Н.А. Вундер", А.С. Павлов", А.В. Ушаков"

a Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация Автор для переписки: 79214215187@ya.ru

Информация о статье

Поступила в редакцию 17.12.16, принята к печати 15.02.17 doi: 10.17586/2226-1494-2017-17-2-269-278 Язык статьи - русский

Ссылка для цитирования: Вундер Н.А., Павлов А.С., Ушаков А.В. Влияние собственных векторов на траектории дискретных систем с задержкой сигнала управления // Научно-технический вестник информационных технологий, механики и оптики. 2017. Т. 17. № 2. С. 269-278. doi: 10.17586/2226-1494-2017-17-2-269-278

Аннотация

Предмет исследования. Представлены результаты исследования свободного движения линейных дискретных систем в случае, когда сигнал управления передается с задержкой, не превышающей длительности интервала дискретности, в непрерывный технический объект. Установлено, что задержка вывода сигнала управления, не превышающая длительности интервала дискретности, увеличивает размерность дискретного представления непрерывного объекта на единицу. Рассмотрен класс дискретных систем с матрицей состояния простой структуры и структурой собственных векторов, матрица которых обладает числом обусловленности, заметно превышающим единицу. Метод. Задача решена с использованием векторно-матричного формализма метода пространства состояния, для чего исследуемая система представлена в векторно-матричной форме. Аналитические оценки процессов в системе сформированы по норме вектора состояния свободного движения с вычислением числа обусловленности матрицы собственных векторов матрицы состояния системы. Основные результаты. Обнаружено, что при плохо обусловленной структуре собственных векторов матрицы состояния дискретной системы возможно возникновение заметных отклонений траекторий свободного движения системы от монотонно убывающей кривой. Практическая значимость. Указанные свойства траекторий свободного движения дискретных систем рекомендовано держать под контролем как на этапе проектирования, так и на этапе эксплуатации систем. Ключевые слова

дискретная система, задержка, технический объект, собственные вектора, отклонения, свободное движение, число

обусловленности

Благодарности

Работа выполнена при поддержке правительства Российской Федерации (грант 074-U01), Министерства образования и науки Российской Федерации (проект 14. Z50.31.0031), гранта президента Российской Федерации №14.У31.16.9281-НШ.

INFLUENCE OF EIGENVECTORS ON TRAJECTORIES OF DISCRETE SYSTEMS

WITH CONTROL SIGNAL DELAY N.A. Vunder", A.S. Pavlov", A.V. Ushakov"

a ITMO University, Saint Petersburg, 197101, Russian Federation Corresponding author: 79214215187@ya.ru

Article info

Received 17.12.16, accepted 15.02.17 doi: 10.17586/2226-1494-2017-17-2-269-278 Article in Russian

For cit"tion: Vunder N.A., Pavlov A.S., Ushakov A.V. Influence of eigenvectors on trajectories of discrete systems with control signal delay. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2017, vol. 17, no. 2, pp. 269-278 (in Russian). doi: 10.17586/2226-1494-2017-17-2-269-278

Abstract

Subject of Study. We present research results of the free motion of discrete linear systems in case of the control signal transmissionto continuous plantwith a delay not exceeding the discreteness interval duration. It was foundthat if the output control signal delay does not exceed the discreteness interval durationthen the dimension of the discrete representation of the continuous plant is increased by one. We consider a class of discrete systems with simple structure state matrix and with the eigenvectors structure such that condition number of their matrix is much greater than one.Method. The problem is solved with the use of vector-matrix formalism of state space method;for that reason the system is presented in vector-matrix form. Analytical assessments of processes in the system arefigured on the norm of state vectorfree motion with the calculation of

the condition number of the state matrix eigenvectors matrix. Main Results. It has been found that the ill-conditioned eigenvectors structure of the state matrix may cause significant deviations of free motion trajectories from a monotone decreasing curve. Practical Relevance. These properties of the free motion of discrete systems trajectories are recommended for being kept under supervision, both in the design and the operational phases of systems. Keywords

discrete system, delay, plant, eigenvectors, deviations, free motion, condition number Acknowledgements

This work was supported by the Government of the Russian Federation, Grant 074-U01; the Ministry of Education and Science of the Russian Federation, Project 14.Z50.31.0031; the Russian Federation President Grant No.14.Y31.16.9281-HO.

Введение

Прежде чем исследовать проблему, вынесенную в заголовок статьи, остановимся на понятии дискретного объекта управления, которое установилось к настоящему моменту [1]. Под дискретным объектом управления понимается функциональное объединение цифрового устройства, формирующего цифровой сигнал управления, цифро-аналогового преобразователя и непрерывного технического объекта. Потребность получения единого математического описания указанной выше функциональной цепи приводит к необходимости исследования процессов в моменты t = k-At, где k - целое положительное число, выполняющее функцию дискретного времени, выраженного в числе интервалов дискретности длительности At. Таким образом, под дискретным объектом управления принято понимать дискретную по времени выборку из непрерывных переменных состояния и выхода при кусочно-постоянном управлении длительности At. Следует отметить, что вывод управления из цифрового устройства может происходить как без задержки, так и с задержкой на величину т . Это обстоятельство порождает две возможные версии дискретных представлений дискретного объекта с одним и тем же непрерывным техническим объектом, что проявляется в увеличении размерности матриц дискретной модели непрерывного объекта и в модификации структуры собственных векторов матрицы состояния дискретной системы. Таким образом, обнаруживается, что не только кратность собственных значений [2-4], но и способ вывода сформированного сигнала управления в технический объект может быть причиной заметных отклонений [5-11] нормы вектора состояния свободного движения дискретной системы.

Влияние собственных векторов на траектории дискретных систем в случае отсутствия задержки вывода сигнала дискретного управления в технический объект

Рассматривается непрерывный технический объект (ТО)

x(t) = Ax (t) + Bu (t); x(0) = x(t)|^, y (t) = Cx (t) , (1)

где x, u, y - соответственно вектора состояния, управления, выхода; A, B, C - соответственно матрицы состояния, управления и выхода ТО; при этом перечисленные векторно-матричные компоненты имеют размерности dim (x) = n , dim (u) = r , dim (y ) = m , dim (A ) = (n x n) , dim (B) = (n x r ) ,

dim (C) = (n x m ) . Если управление непрерывным объектом (1) при t = kAt осуществляется без задержки с помощью цифрового устройства (ЦУ) и объединенного с ним цифро-аналогового преобразователя (ЦАП), то оно представимо в форме

u(t)= u(k), kAt < t <(k + 1)At (2)

так, что объединением (1) и (2), следуя [12], получаем дискретное описание непрерывного ТО (1) или просто дискретный ТО, имеющее вид

x(k +1) = Ax(k) + Bu(k);x(0) = x(k)|k 0, y (k) = Cx(k) , (3)

где k = arg (t = kAt) - дискретное время, выраженное в числе интервалов дискретности длительности At,

dim (A )= dim (A ) , dim (B) = dim (B ) , dim (C ) = dim (C ) , при этом A = exp (AAt), B = (A -1) A-1B,

C = с .

Аналитически управление (3) записывается в виде

u(k)= Kgg(k)-Kx(k), (4)

где g - внешнее воздействие; Kg, K - соответственно матрица прямых связей по внешнему воздействию и матрица обратных связей по вектору состояния ТО, при этом перечисленные векторно-матричные компоненты имеют размерности dim (g) = m , dim (Kg )=(r x m), dim (K) = (r x n).

Объединение формирователя сигнала (4) и дискретного ТО (3) дает замкнутую дискретную систему вида

x (k +1) = Fx (k) + Gg (k); x (0) = x (k)k=0; y (k) = Cx (k ) ; e(k) = g(k)-y (k), (5)

где _

F = А - BK, G=ВКв, (6)

£(к) - ошибка воспроизведения системой внешнего воздействия.

В случае, когда внешнее воздействие недоступно непосредственному измерению, закон управления (4) модифицируется и представляется в форме

u (к) = К g(k)- Кх( к) = К,г (к)- Кг х (к),

где Ke = Kg = arg(c(l-f) BKg = i) = (c(I-f) b)

Kr = k-kyq_ _

r y ik ,, =к g

(7)

Kr - матрица обратной связи

по части компонентов вектора состояния.

Структурное представление физической реализации системы (5) с формирователем сигнала управления в форме (7) в ее составе представлено на рис. 1.

Ф) т к£ u(k) в x(k+1)

ЭЗ

x(k)

А

y(k)

Рис. 1. Схема дискретной системы для случая отсутствия задержки в выводе сигнала управления Заметим, что модификация управления в форме (7) не меняет математического представления системы в форме (5). Спектр собственных чисел и собственные вектора матрицы F задаются в виде

о{F} = {X,. = arg(det((I-F) = o):Im(I,.) = 0,I,. *Ij;i, j = Щi * jj, Ff, = I, Ü; i = ~n.

Выделим в дискретной системе (5) автономную составляющую

x (к +1) = Fx (k); x (0) , (8)

решение которой принимает вид

x(к) = Fkx(0). (9)

Для матрицы состояния F системы (5) выполняется [13] отношение матричного подобия вида МЛ = FM , (10)

где Л = exp(ЛД/) = diag|ä,. = exp(Ä,.Д);i = 1,nj, M = row^Mi = 4:F4,- = (¡4,-;i = 1,nj, представленная конечным рядом из собственных векторов матрицы F .

Использование (9) и (10) позволяет записать матричное выражение

МЛк = Fk М.

Композиция соотношений (9) и (11) порождает цепочку матрично-векторных равенств: X (к) = ¥к X (0) = МЛкМ-1х (0). Перейдем в (12) к нормам, тогда получим

'(k )|| = | |Fк x (0)|<| |F k||| |x

l||x(0)|=l

= Fk

(11) (12)

(13)

Поставим задачу формирования покрытия процесса, заданного в форме (13):

||х (к )|| = | |М Л к М-1х (0 )||< |МЦ |Л к||| |М-1||| |х (0 )| = С {М}|| ( Хк ); / = 1, „|| |х (0 ) = С {М}Х Ц |х (0 ) ,(14) где С {М} = ||М|||М'Ц - число обусловленности [14, 15] матрицы М собственных векторов матрицы

состояния дискретной системы (5); Хшгк - максимальное собственное значение матрицы F, которое

удовлетворяет условиям 1т (тах) = 0, > 0.

Анализ качества процессов свободного движения дискретной системы для случая, когда сигнал управления направляется в непрерывный объект без задержки, будем проводить на основе соотношения

(14) применительно к оценке Бир (||х (к)Ц) сверху для нормы ||х (к)||, которая при ||х (0) = 1 получает представление

Kr

Бир

К«*

1 = C {м}х,

k

шах *

(15)

Выражение (15) позволяет констатировать:

1. информационным компонентом, определяющим качество процессов по норме свободного движения, является число обусловленности матрицы М собственных векторов;

2. если собственные вектора матриц состояния системы (5) ортогональны, то число обусловленности С {М} будет близко к единичному значению, и процессы в системе, стартуя в точке ||х (0)||, будут

иметь монотонно сходящийся характер, покрытие которых кривой Бир (||х (к)Ц) будет стартовать из

точки С {М}|| х (0)|| = | |х (0)||;

3. если собственные вектора матриц состояния системы (5) не будут ортогональными и будут содержать хотя бы пару векторов, близких к коллинеарной структуре, то число обусловленности С {М} может

быть чрезмерно большим. При этом процессы в системе (5), стартуя в точке ||х (0)||, будут иметь чрезмерно большие отклонения от монотонно сходящихся процессов, покрытие которых кривой Бир (||х (к)||) будет стартовать из точки С {М}||х (0)|| >>||х (0)||, имеющей чрезмерно большое удаление от начала координат.

Влияние собственных векторов на траектории дискретных систем в случае наличия задержки вывода сигнала дискретного управления в технический объект

Если управление (2) непрерывным объектом для t = к А/ осуществляется с задержкой на величину т , не превышающую длительности А/ интервала дискретности, с помощью ЦУ и объединенного с ним ЦАП, то управление и () представимо в виде [16]

[и (к -1), к А/ < / < к А/ + т;

[и (к), к А/ + т < / <(к + 1)А/.

При этом объединение (16) с моделью непрерывного объекта (1) дает описание [16] дискретного объекта в виде

х (к + 1) = Ах (к) + Б; (т)и (к-1) + В (т) и (к), у (к) = Сх (к), (17)

где В; (т) = А (I - е-Ат) А-1В , В (т) = (Ае-Ат -1) А-1В .

Структурное представление объекта (17) приведено на рис. 2.

и(к)

и (/) = '

(16)

Рис. 2. Схема технического объекта для случая наличия задержки в выводе сигнала управления Сделаем дополнительную разметку первого элемента задержки (ЭЗ), приписав его выходу переменную х (к), а его входу - х (к +1) . Тогда в соответствии со структурной схемой и с учетом сделанной

разметки первого ЭЗ для объекта управления, при наличии задержки вывода сигнала дискретного управления в технический объект, можно записать:

1(к +1) =

[х (к +1)] А В, (т) х (к) В (т)

X (к +1) 0 0 _х (к)_ 1 I

где А =

"А В. (т)] [В (т)]

1 V ) ; В = V )

_ 0 0 _ I

и (к) = А х (к) + В и (к); у (£) = С х (к),

; С = [С 0], <Иш (х) = <Иш (и) = г .

(18) (19)

Рассмотрим два случая учета влияния на динамику системы появления задержки в выводе сигнала управления из устройства управления в объект при условии, что г = 1.

Первый случай назовем «случай непредсказуемого появления задержки» (иначе случай неучтенной возможности появления задержки), при котором закон управления сохраняется в форме (4), но в силу модификации объекта управления в форме (18), (19) происходит расширение размерности матрицы состояния системы так, что система получает представление

i(k +1) = F(т)x(k) + G(т)g(k);x(0) = x(k)|k=0;y (k) = Cx(k),

где

F (т ) =

[A -B (т)К B1 (т)] • G (т ) = B (т) K g "

_ K 0 _ • V / _ K g _

; C = [C 0] .

(20) (21)

Структурное представление физической реализации системы (20) с формирователем сигнала управления в форме (7) в ее составе представлена на рис. 3.

g(k) s(k) Ke u(k)

X(k+1)

В(т)

x(k+1)

ЭЗ Bi(T)^^—- ЭЗ %(k)

x(k)

А

y(k)

Рис. 3. Схема дискретной системы для случая появления непредсказуемой задержки в выводе сигнала

управления

Для свободного движения дискретной системы (20) по аналогии с (12) можно записать

x (k) = (F (т))k x (0) = МЛkM-1x (0) , (22)

что, в свою очередь, по аналогии с (14) позволяет сформировать представление

||x(k)|| = ||(F (т))k x(0)| = ||МЛkM-1x(0)| < С{MJX^ ||x(0)| , (23)

где С {MJ - число обусловленности (n +1) х (n +1) -матрицы

М = row^MM-i = arg (FF (т) ММf = X..M..); i = 1, n + 1J собственных векторов матрицы состояния системы (20);

Xmax = max {X.. = arg (И - F (т) = 0); i = 1, n + 1J . Таким образом, для свободного движения системы (20) сформировано покрытие процессов по норме вектора состояния системы вида

SUP ( (k»U = С {ММ} Xmkax.

Следует отметить, что изменение числа обусловленности С {МJ произойдет уже при т = 0, хотя спектр собственных чисел матрицы F (т) = F (0) пополнится по сравнению со спектром собственных чисел матрицы F на нулевое (n +1) -ое собственное число Xn+1 = 0 . Для подтверждения этого предположения рассмотрим F (0), записанную в развернутой форме:

[ A-BK 0 ] -K 0

В силу свойства [13] блочно-треугольной матрицы спектр собственных чисел матрицы (24) будет составлен из спектра собственных чисел матрицы о {F = A - BKJ = {X= X. ; i = 1, nJ и Xn+1 = 0 . Для собственных векторов £t матрицы состояния F (0) запишем

F(0)£i = Xi£i ^(F(0)-Xii)£. = 0.

Нетрудно видеть, что собственный вектор £. принадлежит нуль-пространству n{F (0)- XiiJ характеристической матрицы (F (0) - XiI), которое имеет свою реализацию для каждого собственного числа X.. , в том числе для X^ = 0, содержащую ненулевой собственный вектор. В силу наличия в матрице F (0) внедиагонального элемента -K, несмотря на то, что в спектре собственных чисел матрицы F (0) имеется спектр собственных чисел матрицы F , в собственных векторах £., соответствующих собственным числам {X. = X..; i = 1, nJ, не будут сохраняться компоненты собственных векторов £. матрицы F , что

F (0 ) =

(24)

Kr

влечет за собой заметное изменение числа обусловленности С{м| уже при т = 0. При т Ф 0 матрица

Р (т) принимает вид (21), в котором наличествуют два внедиагональных элемента и претерпевает изменение ненулевой диагональный элемент, что влечет за собой изменение спектра собственных чисел, собственных векторов и числа обусловленности С {М|.

Второй случай назовем «случай предсказуемого появления задержки» (иначе случай учитываемой возможности появления задержки, для которого т Ф 0: т - неопределенная фиксированная величина, не

превышающая интервала дискретности). Здесь закон управления и (к) формируется как функция состояния объекта (18) и задающего воздействия § (к) в форме

и (к) = % g (к)- К х (к)- ¿х (к) = ¿г (к)- Кг X (к)- ¿х (к). (25)

Объединение объекта управления (18) и закона управления (25) порождает замкнутую дискретную систему, записываемую в виде

X (к +1) = Р (т) X (к) + 6 (т) g (к); X (0) = X (к)|к=о; у (к) = СX (к), (26)

; С = [С 0],

причем управление (25) формируется из условия сохранения в спектре собственных чисел матрицы Р всех элементов спектра собственных чисел матрицы Р (6), при этом дополнительное (п +1) -ое собственное число 1 п+1 выбирается из условия 1 п<< 1,.,1 = 1, п .

Структурное представление физической реализации системы (26) с формирователем сигнала управления в форме (25) в ее составе представлено на рис. 4.

A - B (т)K, B1 (т) - B (т )K X B (т )K g

где F = V / X ; G = V / g

-к , KX . _ Kg _

**) «(k) к

Рис. 4. Схема дискретной системы для случая появления предсказуемой задержки в выводе сигнала

управления

Для свободного движения дискретной системы (26) по аналогии с (22) можно записать X (к) = (Р (т))* X (0) = МЛкМ-^ (0), что, в свою очередь, по аналогии с (23) позволяет сформировать представление

< с {мл11 т

М Лk М"

где

■x (о )||< с {мм j а ц |x (о )||,

обусловленности

число

(n +1) х (n +1) -матрицы

|x (*)| = |(F (т)) x (0) с {мм j

М = row{Mt = arg(FF(т)Mt = O..M..);i = 1,n + lj собственных векторов матрицы состояния системы (25);

а max = max {а t = arg (OI - FF (т ) = 0); i = 1, n + lj. Таким образом, для свободного движения системы (25) сформировано покрытие процессов по норме вектора состояния системы вида

sup (l|x

x (0)1=1

{MF

= с {М5 O

Компьютерное исследование дискретной системы при различных способах вывода сигнала

управления в технический объект

Компьютерное исследование свободного движения дискретной системы с матрицей состояния простой структуры по норме вектора состояния как функции структуры собственных векторов для различных случаев вывода сигнала управления в ТО, дискретное представление которого имеет вид (3), при

А/ = 0,01с с матрицами А =

1 0,01" [0,0001"

А = ; в =

0 1 _ _ 0,01 _

; С = [1 0] проводилось в полном соответствии с

текстом статьи.

На рис. 5 представлена кривая свободного движения автономной дискретной системы (8) по норме ее вектора состояния в форме (13). Система имеет размерность п = 2, матрицу состояния - Г 0,999 0,0094"

Ё = , которая доставляется системе с помощью матрицы обратной связи

-0,1884 0,886

К = [18,84 11,4], имеющей алгебраический спектр собственных значений

г_, г_ - , Г 0,4472 " [0,0995"

о{Ё}={Х1 = 0,9802;= 0,9048} и собственные вектора £. = ;£2 = . Матрица послед-

1 > 1 1 2 > 1 -0,8944 2 -0,995

них характеризуется числом обусловленности С {М} = С {[£1, £2 ]} = 5,434 . Число обусловленности матрицы собственных векторов заметно больше единицы, и в траектории системы по норме ||х (к)Ц вектора

состояния должно наблюдаться заметное отклонение от монотонного характера ее изменения, что обнаруживает кривая рис. 5.

1|х(/)||

1

0,5

0 0,5 1 1,5 2 2,5 /=А/к Рис. 5. Кривая процесса ||х (к )||

На рис. 6 представлена кривая свободного движения автономной версии дискретной системы (20)

по норме ее вектора состояния в

форме ||х (к)| = (Ё (т))к х (0) < (Ё (т))к ||х (0)|

II* (0)1=1

которая характери-

Ё (т) =

А - В(т)К Б; (т)

-К 0

зуется размерностью п = 3 и матрицей состояния

1-18,84(0,5А/-1,5т)(А/-т) 0,01 -11,4(0,5А/-1,5т)(А/-т) (А/-0,5т)т~ -18,84 (А/-т) 1-11,4(А/-т) т

-18,84 -11,4 0

где А/ = 0,01 с, т = 0,5 А/, имеющей алгебраический спектр собственных значений о {Ё (т)} = {X! = 0,9801; X2 = 0,8987; X3 = 0,0644 } и собственные вектора

матрица которых характеризуется числом обусловленности

С {М} = С {[£;, £2, £3 = 70,674 . Число обусловленности матрицы собственных векторов заметно больше единицы, и в траектории системы по норме || х (к) || вектора состояния должно наблюдаться заметное

отклонение от монотонного характера ее изменения, что обнаруживает кривая рис. 6.

На рис. 7 представлены кривые свободного движения автономной версии дискретной системы (26)

[-0,2179" [ 0,0089 " [ 0 "

1г = 0,4352 ; 12 = -0,0933 ; 13 = 0,0057

-0,8735 0,9956 -1

по норме ее вектора состояния в

форме ||х (к)|= (Ё (т))к х (0) < (Ё (т))к ||х (0)

которая характери-

IIх (°)Н

зуется размерностью п = 3 , матрицей состояния Р (т)

т=0,5Д/ = 0,005

для

[ 1 0,01 0 " -0,0171 0,9888 0,0086 -3,4157 -2,2468 0,715

рис. 7, а, которая доставляется системе с помощью матрицы обратной связи К = [3,4157 2,2468 -0,715], имеющей алгебраический спектр собственных значений

о { (т)} = { = 0,9802; 12 = 0,9048; 13 = 0,81871 и собственные

$1 =

сти

Р (т )

[ 0,2192 " [-0,0102" [-0,0026"

-0,4363 ; $ 2 = 0,0996 ; $ 3 = 0,05

0,8727 -0,995 -0,9987

С

{М) = С {

$1, $ 2, $ 3

т=0,5Д/ = 0,005

вектора

матрица которых характеризуется числом обусловленно-= 131,9. Для кривой (рис. 7, б) матрица состояния вида которая доставляется системе с помощью матрицы обрат-

1 0,01 0 -0,0009 0,9985 0,0094 -0,1875 -0,3014 0,8766

ной связи К = [0,1875 0,3014 -0,8766], имеющей алгебраический спектр собственных значений

5{Р (т )) = {{

= 0,9802; 12 = 0,9048; 13 = 0,99

собственные

вектора

[ 0,2192 " [-0,0102" [-0,5783"

$1 = -0,4363 ; $ 2 = 0,0996 ; $ 3 = 0,5769

0,8727 -0,995 -0,5769

матрица которых характеризуется числом обусловленно-

сти С

{М)= С {

$ 1, $ 2, $$3

= 17,5787 .

нон

15 10 5

0 0,2 0,4 0,6 0,8 1 1,2 t=Дtk Рис. 6. Кривая процесса || X(к) || для случая непредсказуемого появления задержки

1М*)|| 10

8

6

4

2

t=Дtk

l|x(t)|| 2,5

2

1,5 1

0,5 0

2 3 б

4 t=Дtk

Рис. 7. Кривые процесса || X(к) || в случае предсказуемого появления задержки для различных значений третьего собственного числа: 13 = 0,8187 (а); 13 = 0,99 (б)

1

0

а

Полученные кривые обнаруживают, что с увеличением модуля третьего собственного значения 13 уменьшается число обусловленности матрицы собственных векторов, что свидетельствует о демпфировании отклонения траектории дискретной системы по норме вектора состояния, при этом сходимость процессов в системе к неподвижному состоянию x (да) = 0 существенно замедляется.

Заключение

Установлено, что устойчивые дискретные системы с простым спектром собственных чисел матрицы состояния системы и структурой ее собственных векторов, матрица которых характеризуется числом обусловленности, большим единицы, обнаруживают заметные отклонения траекторий от монотонно убывающего характера своего поведения в режиме свободного движения. Это отклонение усиливается в случае, если вывод сигнала управления из цифрового формирователя этого сигнала в технический объект осуществляется с запаздыванием, не превышающем интервала дискретности. Обнаруженные отклонения могут быть столь значительными, что могут привести к разрушительным последствиям. Указанные свойства траекторий свободного движения дискретных систем следует держать под контролем как на этапе проектирования, так и на этапе эксплуатации систем.

Литература References

1. Ту Ю. Современная теория управления. М.: Машиностроение, 1971. 472 с.

2. Акунов Т.А., Дударенко Н.А., Полинова Н.А., Ушаков А.В. Исследование колебательности процессов в апериодических непрерывных системах, порождаемой фактором кратности собственных чисел // Научно-технический вестник информационных технологий, механики и оптики. 2013. № 3 (85). С. 55-61.

3. Вундер Н.А., Ушаков А.В. Исследование отклонения траекторий линейных устойчивых дискретных систем от монотонно убывающей кривой свободного движения в случае кратных вещественных положительных собственных чисел их матриц // Проблемы управления и информатики. 2016. № 5. С. 1-11.

4. Vunder N.A., Ushakov A.V. On the features of the trajectories of autonomous discrete systems generated by a multiplicity factor of eigenvalues of state matrices // Proc. IEEE Int. Symposium on Intelligent Control. Sydney, Australia, 2015. P. 695-700. doi: 10.1109/ISIC.2015.7307291

5. Whidborne J.F., McKernan J. On the minimization of maximum transient energy growth // IEEE Transactions on Automatic Control. 2007. V. 52. P. 1762-1767. doi: 10.1109/TAC.2007.900854

6. Polyak B.T., Smirnov G. Large deviations for non-zero initial conditions in linear systems // Automatica. 2016. V. 74. P. 297-307. doi: 10.1016/j.automatica.2016.07.047

7. Whidborne J.F., Amar N. Computing the maximum transient energy growth // BIT Numerical Mathematics. 2011. V. 51. N 2. P. 447-457. doi: 10.1007/s10543-011-0326-4

8. Halikias G., Papageorgiou A., Karcanias N. Non-overshooting stabilization via state and output feedback // International Journal of Control. 2010. V. 83. N 6. P. 1232-1247. doi: 10.1080/00207171003664844

9. Francis B.A., Glover K. Bounded peaking in the optimal linear regulator with cheap control // IEEE Transactions on Automatic Control. 1978. V. 23. N 4. P. 608-617. doi: 10.1109/TAC.1978.1101792

10. Kimura H. A new approach to the perfect regulation and the bounded peaking in linear multivariable control systems // IEEE Transactions on Automatic Control. 1981. V. 26. N 1. P. 253-270. doi: 10.1109/TAC.1981.1102573

11. Vunder N.A., Ushakov A.V. Investigation of continuous systems oscillatory processes created with the multiplicity factor of the eigenvalues of the state matrices // Journal of Automation and Information Sciences. 2015. V. 47. N 11. P. 18-35. doi: 10.1615/JAutomatInfScien.v47.i11.20

12. Zadeh L.A., Desoer C.A. Linear System Theory: The State Space Approach. 4th ed. NY: Dover Publ., 2008. 656 p.

13. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1973. 575 с.

14. Голуб Дж., Ван Лоун Ч. Матричные вычисления: Пер. с англ. М.: Мир, 1999. 548 с.

1. Tou J.T. Modern Control Theory. NY, McGraw-Hill, 1964, 343 p.

2. Akunov T.A., Dudarenko N.A., Polinova N.A., Ushakov A.V. Process oscillativity study in aperiodic continuous systems, generated by eigenvalues multiplication factor. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2013, no. 3, pp. 55-61. (In Russian)

3. Vunder N.A., Ushakov A.V. Investigation of deviation of trajectories of linear stable discrete systems from monotonously decreasing curve on free motion in the case of multiple real positive eigenvalues of their matrices. Journal of Automation and Information Sciences, 2016, vol. 48, no. 10, pp. 13-24. doi: 10.1615/JAutomatInfScien.v48.i10.20

4. Vunder N.A., Ushakov A.V. On the features of the trajectories of autonomous discrete systems generated by a multiplicity factor of eigenvalues of state matrices. Proc. IEEE Int. Symposium on Intelligent Control. Sydney, Australia, 2015, pp. 695-700. doi: 10.1109/ISIC.2015.7307291

5. Whidborne J.F., McKernan J. On the minimization of maximum transient energy growth. IEEE Transactions on Automatic Control, 2007, vol. 52, pp. 1762-1767. doi: 10.1109/TAC.2007.900854

6. Polyak B.T., Smirnov G. Large deviations for non-zero initial conditions in linear systems. Automatica, 2016, vol. 74, pp. 297-307. doi: 10.1016/j.automatica.2016.07.047

7. Whidborne J.F., Amar N. Computing the maximum transient energy growth. BIT Numerical Mathematics, 2011, vol. 51, no. 2, pp. 447-457. doi: 10.1007/s10543-011-0326-4

8. Halikias G., Papageorgiou A., Karcanias N. Non-overshooting stabilization via state and output feedback. International Journal of Control, 2010, vol. 83, no. 6, pp. 1232-1247. doi: 10.1080/00207171003664844

9. Francis B.A., Glover K. Bounded peaking in the optimal linear regulator with cheap control. IEEE Transactions on Automatic Control, 1978, vol. 23, no. 4, pp. 608-617. doi: 10.1109/TAC. 1978.1101792

10. Kimura H. A new approach to the perfect regulation and the bounded peaking in linear multivariable control systems. IEEE Transactions on Automatic Control, 1981, vol. 26, no. 1, pp. 253-270. doi: 10.1109/TAC.1981.1102573

11. Vunder N.A., Ushakov A.V. Investigation of continuous systems oscillatory processes created with the multiplicity factor of the eigenvalues of the state matrices. Journal of Automation and Information Sciences, 2015, vol. 47, no. 11, pp. 18-35. doi: 10.1615/JAutomatInfScien.v47.i11.20

12. Zadeh L.A., Desoer C.A. Linear System Theory: The State Space Approach. 4th ed. NY, Dover Publ., 2008, 656 p.

13. Gantmakher F.R. Teoriya Matrits [Matrix Theory]. Moscow, Nauka Publ., 1973, 575 p. (In Russian)

14. Golub G.H., Van Loan C.F. Matrix Computations. Baltimore, Johns Hopkins University Press, 1996,. 728 p.

15. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.: Наука, 1984. 320 с.

16. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А. В. Синтез дискретных регуляторов при помощи ЭВМ. Л.: Машиностроение, 1983. 245 с.

Авторы

Вундер Нина Александровна - аспирант, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, polinova_nina@mail.ru

Павлов Андрей Сергеевич - аспирант, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, a.s.pavlov@email.su

Ушаков Анатолий Владимирович - доктор технических наук, профессор, профессор, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, 79214215187@ya.ru

15. Voevodin V.V., Kuznetsov Yu.A. Matritsy i Vychisleniya [Matrices and Calculations]. Moscow, Nauka Publ., 1984, 320 p.

16. Grigor'ev V.V., Drozdov V.N., Lavrent'ev V.V., Ushakov A.V. Sintez Diskretnykh Sistem Regulyatorov pri Pomoshchi EVM [Synthesis of Digital Systems Using Computer Controls]. Leningrad, Mashinostroenie Publ., 1983, 245 p.

Authors

Nina A. Vunder - postgraduate, ITMO University, Saint Petersburg, 197101, Russian Federation, polinova_nina@mail.ru

Andrey S. Pavlov - postgraduate, ITMO University, Saint Petersburg, 197101, Russian Federation, a.s.pavlov@email.su

Anatoliy V. Ushakov - D.Sc., Professor, ITMO University, Saint Petersburg, 197101, Russian Federation, 79214215187@ya.ru