Например, Бобцов

КОНСТРУИРОВАНИЕ ВЕЩЕСТВЕННОЗНАЧНОЙ КРИТЕРИАЛЬНОЙ МАТРИЦЫ ДЛЯ ОДНОКАНАЛЬНОЙ СИСТЕМЫ

62 Н. А. Дударенко, М. В. Полякова, А. В. Ушаков
УДК 62.50

Н. А. ДУДАРЕНКО, М. В. ПОЛЯКОВА, А. В. УШАКОВ
КОНСТРУИРОВАНИЕ ВЕЩЕСТВЕННОЗНАЧНОЙ КРИТЕРИАЛЬНОЙ МАТРИЦЫ
ДЛЯ ОДНОКАНАЛЬНОЙ СИСТЕМЫ

Рассматривается проблема контроля вырождения динамической системы типа „многомерный вход—многомерный выход“. Решение задачи базируется на конструировании инструментария контроля вырождения критериальной матрицы.
Ключевые слова: динамическая система, функционал, критериальная матрица, матричная геометрическая прогрессия.

Введение. Сложная динамическая система типа „многомерный вход—многомерный

выход“ (МВМВ) в математической постановке ее функционирования реализует линейный

(локально линейный) оператор (ЛО), отображающий пространство целевых намерений в про-

странство осуществляемых реализаций, при этом размерности этих пространств полагаются

равными.

В математической постановке линейный оператор оказывается вырожденным [1], если

его ранг становится меньше размерности указанных пространств. Множество причин может

привести к вырождению ЛО сложной динамической системы МВМВ-типа. Система может

вырождаться, когда из ее состава выпадает некоторый функциональный элемент, в результате

система становится функционально неполной. Как следствие, сокращается размерность про-

странства осуществляемых реализаций. Причины вырождения могут носить организацион-

ный характер, когда формируемые целевые намерения неудачно распределяются по входам

сложной динамической системы. Вырождаться могут системы вследствие параметрической

неопределенности, когда не должным образом организованы связи между каналами системы

МВМВ-типа, когда неправильно назначены по знаку и величине коэффициенты передачи

этих связей, когда неудачно сформированы полосы пропускания каналов.

Технология контроля вырождения. Для построения технологии контроля вырожде-

ния в соответствии с математическим определением вырождения необходимо представить

cложную динамическую систему МВМВ-типа в виде линейной алгебраической задачи (ЛАЗ),

матрица которой задает оператор, отображающий пространство целевых намерений в про-

странство осуществляемых реализаций

η(w) = N (w, θ)χ(w) ,

(1)

где N (w, θ) — (т×т)-матрица для любых w , θ ; η(w) , χ(w) — р-мерные векторы, θ — р-мерный параметр, изменяющий алгебраические свойства матрицы N .
Инструмент контроля вырождения построен на спектре сингулярных чисел α j , функ-
ционалов Jν (ν = m, 1) с использованием предложенной ЛАЗ (1). Для контроля вырождения воспользуемся функционалами Jν , которые строятся на спектре

{ ( ) }σα {N} = α j = µ1j/ 2 :µi :det µI − NT N = 0: j =1, m

(2)

( )сингулярных чисел α j j =1, m критериальной матрицы N и в силу соотношений

Jν {N} = αν {N} / α1 {N}; ν = m,1.

(3)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

Конструирование вещественнозначной критериальной матрицы для одноканальной системы 63 Воспользуемся сингулярным разложением матрицы N (SVD-процедурой) [1, 2]:

N =U N ΣNVNT ,

(4)

N VN j = α jU N j , j =1, p .

(5)

Векторно-матричное соотношение (5) придает исходной линейной алгебраической за-

даче (1) геометрический смысл такой, что вектор χ =VN j ( j =1, p) отражается в подпростран-

ство, натянутое на j-й элемент U N j левого сингулярного базиса U N так, что соответствую-

щий ему вектор имеет норму, равную α j . Тогда задача контроля вырождения формализуется

как оценка перехода критериальной матрицы N из сферы, расположенной в пространстве,

натянутом на векторы χ , в эллипсоид, натянутый на левый сингулярный базис U N с размерами полуосей, совпадающими с сингулярными числами матрицы N .
При J =1 (см. (3)), матрица N вырождается, т.е. происходит „сплющивание“ этого эл-
липсоида вдоль его p-й полуоси, т.е. вдоль р-го левого сингулярного вектора UN p . Нетрудно

видеть, что если параметр θ изменяет матрицу N (θ) таким образом, что последовательно,

начиная с α p , принимают нулевые значения остальные ( p −1 ) сингулярных чисел, кроме α1 ,

то в пространстве, натянутом на левый сингулярный базис, будет наблюдаться последова-

тельное „сплющивание“ эллипсоида вдоль векторов UN p , UN p−1 , …, U N 2 . В итоге сфера

отобразится в отрезок прямой.

Функционалы Jν используются для количественной оценки вырождения, они конструируются на спектре сингулярных чисел αJ критериальной матрицы системы.
Свойство 1. Функционалы Jν критериальной матрицы N обладают свойством

0 ≤ Jν {N}≤1; ν =1, m .

(6)

Свойство 2. Функционалы Jν не зависят от умножения критериальной матрицы N на
скаляр α , в общем случае параметризованный неким параметром γ , так что α = α (γ) , что

представимо в форме

J Gν {α (γ)[N ]}= J Gν {[N ]}= J Gν ; ν =1, m .

(7)

Свойство 3. Функционалы Jν критериальной матрицы N не зависят от умножения

критериальной матрицы слева или справа на ортогональную матрицу.

Свойство 4. Глобальный функционал JG критериальной матрицы N совпадает с вели-

чиной, обратной числу С{N} обусловленности этой матрицы

{ }C{N}= {N} N − 1 ,

(8)

если в (8) использованы спектральные нормы матриц, так что выполняется равенство

J [N ]= C− 1 [N ]= αm α1 .

(9)

Свойство 5. Для глобального функционала JG прямой матрицы N и обратной ей N −1

выполняется соотношение

( )JG ( N ) = JG N − 1 .

(10)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

64 Н. А. Дударенко, М. В. Полякова, А. В. Ушаков
Основной результат. Формирование вещественнозначной критериальной матрицы на основе разложения резолвенты. Для сведения системы МВМВ-типа к ЛАЗ вида (1) воспользуемся аппаратом передаточных матриц, выделив из него передаточную матрицу
Φ(s) отношения „вход—выход“, задаваемую в форме

Φ (s) = arg{y (s) = Φ (s) g (s)} .

(11)

Соотношение (11) задает ЛАЗ в виде (1) и связывает реализации комплекснозначных
векторов входа g(s) и выхода y(s). Для конструирования вещественнозначной версии ЛАЗ, а
следовательно и вещественнозначной критериальной матрицы, ограничимся установившимся значением векторного экзогенного воздействия, установившимся значением вектора скорости этого воздействия и вектора ускорения этого воздействия. Тогда для установившегося значения вектора выхода системы МВМВ-типа оказывается справедливым представление

yуст (t ) = K0 g (t )+ K1g (t )+ K2 g (t )+...K l g(l) (t ) ,

(12)

K0 , K1, K2 , ..., Kl — матричные коэффициенты, связывающие вектор установившегося зна-
чения выхода с векторами производных соответствующих порядков экзогенного воздействия. Нетрудно видеть, что выражение (12) содержит l +1 критериальную матрицу (для физически реализуемых систем l +1= 3 ). Задача исследования вырождения сложной динамической сис-
темы МВМВ-типа, порождаемая установившимся значением g(t ), установившимся значе-
нием скорости его изменения и установившимся значением ускорения, будет решена, если исследователь получит аналитическое представление этих матричных коэффициентов как функции структурных моделей компонентов систем. При этом исследователь получит отдельные критериальные матрицы, имеющие вещественнозначные представления.
Для вычисления введенных матричных коэффициентов Kl (l = 0, 1, 2, ...) воспользуемся
разложением резолвенты, входящей в состав аналитического представления передаточной
функции Φ(s) системы МВМВ-типа, задаваемой с помощью векторно-матричного описания
аппарата метода пространства состояния [4, 5]

x (t ) = Fx (t )+ Eg (t ): x (0) = x (t ) t=0 y (t ) = Cx (t )+ Dg (t ) .

(13)

Произведем над выражениями (12) и (13) аналитические преобразования. Для этого к (12) применим преобразование Лапласа, в результате чего получим

( )ууст (s) = K0 + K1s+ K2s2 +... g (s) .

(14)

В свою очередь, применение преобразования Лапласа к выражению (13) позволяет сформи-
ровать выражение для передаточной функции Φ(s) системы МВМВ-типа в форме

Φ (s) =C (sI − F )−1 E + D .

(15)

Очевидно, подстановка выражения (15) в соотношение (11), примененная для случая ус-

тановившегося режима, позволяет записать

{ }yуст (s) = C (sI − F )−1 E + D g (s) .

(16)

Из сравнения выражений (14) и (16) видно, что задача конструирования матричных коэффициентов Kl (l = 0, 1, 2, ...) , являющихся отдельными критериальными матрицами, будет
решена, если в (16) удастся разложить резолвенту ( sI − F )−1 по положительным степеням

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

Конструирование вещественнозначной критериальной матрицы для одноканальной системы 65

комплексной переменной s . Для этого, осуществляя преобразование, приведем резолвенту
(sI − )F −1 к виду

(

sI



F

)−1

=



(

F



sI

)−1

=



⎡⎣(

F



sI

)−1

⎤ ⎦

=



⎡⎣(

F



sI

)−1

⎤ ⎦

=

( ) ( )=



⎡ ⎣

sI − F −1

F

⎤ ⎦

−1

=

−F

−1

I − sF −1

−1 .

(17)

Представление (17) резолвенты позволяет записать ее в мультипликативной форме мат-

рицы F −1 и суммы членов матричной геометрической прогрессии с показателем sF −1

( ) ( )I −sF −1 −1 = −F −1 I + sF −1 + s2 F −2 +...s p F − p = −F −1 − sF −2 − s2 F −3 −...− s p F − p . (18)

Воспользуемся соотношением (18) для представления выражения (16) по положительным

степеням s , тогда получим

(( ) )yуст (s) = I +CF −1E + D +CF −2 Es+CF −3Es3 +...+CF − p+1Es p g ( s) .

(19)

Сравнение выражений (19) и (14) позволяет получить представления матричных коэф-

фициентов Kl в форме

( )K0 = I +CF −1E + D ; Kl = CF −(l+1) E (l =1, 2, …) .

(20)

Нетрудно видеть, что матричные коэффициенты вида (20) представляют собой искомые вещественнозначные критериальные матрицы семейства линейных алгебраических задач вида

yl уст (t ) = Kl g(l) (t ); l = 0, 1, 2 ...

(21)

сложной динамической системы МВМВ-типа, связывающие l-й компонент ylуст (t ) значения
вектора выхода yуст (t ) с установившимися значениями вектора скорости изменения внешне-
го воздействия g(t ), установившимся значением вектора его ускорения и остальных высших
производных. Заключение. Полученное семейство вещественнозначных критериальных матриц по-
зволяет для случая установившегося режима функционирования сложной системы МВМВтипа, который наиболее характерен для решения задач исследования вырождения, оценивать покомпонентно влияние скорости, ускорения и прочих высших производных на возможность вырождения системы указанного типа.

СПИСОК ЛИТЕРАТУРЫ
1. Дударенко Н. А., Слита О. В., Ушаков А. В. Математические основы современной теории управления: аппарат метода пространства состояний: Учеб. пособие. СПб: СПбГУ ИТМО, 2008.
2. Дударенко Н. А., Полякова М. В., Ушаков А. В. Контроль вырождения сложных динамических систем созидательного типа с антропокомпонентами // Науч.-технич. вестн. СПбГУ ИТМО. Вып. 55. Управление, моделирование, информационная безопасность. 2008. С. 25—30.
3. Бочков А. Л., Дударенко Н. А., Ушаков А. В. Синтез многомерных функционально вырожденных динамических систем // Изв. вузов. Приборостроение. 2008. Т. 51, № 1. С. 25—29.
4. Акунов Т. А., Ушаков А. В. Синтез систем гарантированной модальной стабильности // Изв. РАН. Теория и системы управления. 2003. № 4. С. 9—17.
5. Бесекерский В.А., Попов Е.П. Теория систем автоматического управления. СПб: Профессия, 2007.
6. Дудаpенко Н. А., Полякова М. В., Ушаков А. В. Достаточные алгебраические условия обобщенной синхронизируемости многоканальных динамических объектов // Мехатроника, автоматизация, управление. 2009. № 5. С. 17—20.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

66 Н. А. Дударенко, М. В. Полякова, А. В. Ушаков

Наталия Александровна Дударенко Майя Вячеславовна Полякова Анатолий Владимирович Ушаков

Сведения об авторах — канд. техн. наук, доцент; Санкт-Петербургский государственный
университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: dudarenko@yandex.ru — аспирант; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: 12noch@mail.ru — д-р техн. наук, профессор; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail:Ushakov-AVG@yandex.ru

Рекомендована кафедрой систем управления и информатики

Поступила в редакцию 01.07.09 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11