exponenta event banner

Сведения о границах секторов и индексах секторов

Конические секторы

В простейшем виде коническим сектором является 2-D область, ограниченная двумя линиями, y = au и y = bu.

Затененная область характеризуется неравенством (y-au) (y-bu) < 0. В более общем случае любой такой сектор может быть параметризован как:

(yu) TQ (yu) < 0,

где Q является 2x2 симметричной неопределенной матрицей (Q имеет одно положительное и одно отрицательное собственное значение). Мы называем Q секторной матрицей. Эта концепция обобщает более высокие размеры. В N-мерном пространстве конический сектор представляет собой множество:

S={z∈RN:zTQz<0},

где Q снова является симметричной неопределенной матрицей.

Границы секторов

Секторные границы - это ограничения на поведение системы. Ограничения усиления и ограничения пассивности являются особыми случаями границ секторов. Если для всех ненулевых входных траекторий u (t), выходная траектория z (t) = (Hu) (t) линейной системы H (s) удовлетворяет:

∫0TzT (t) Qz (t) dt<0,∀T>0,

тогда выходные траектории Н лежат в коническом секторе с матрицей Q. Выбор различных Q-матриц накладывает различные условия на отклик системы. Например, рассмотрим траектории y (t) = (Gu) (t) и следующие значения:

H (s) = (G (s) I), Q = (0-I-I0).

Эти значения соответствуют границам сектора:

∫0T (y (t) u (t)) T (0-I-I0) (y (t) u (t)) dt<0,∀T>0.

Эта граница сектора эквивалентна условию пассивности для G (s):

∫0TyT (t) u (t) dt>0,∀T>0.

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

H = (GI).

Условие частотной области

Поскольку условие временной области должно поддерживаться для всех T > 0, получение эквивалентной границы частотной области требует небольшой осторожности и не всегда возможно. Позвольте следующее :

Q = W1TW1-W2TW2

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

∫0T (Hu) (t) TQ (Hu) (t) dt<0,∀T>0

эквивалентно состоянию частотной области:

H () HQH () <0∀ω∈R.

Поэтому достаточно проверить неравенство секторов на реальные частоты. Используя разложение Q, это также эквивалентно:

(W1TH)(W2TH)-1 ‖∞ <1.

Следует отметить, что W2TH является квадратным, когда Q имеет столько же отрицательных собственных значений, сколько входных каналов в H (s). Если это условие не выполняется, уже недостаточно (вообще) просто посмотреть на реальные частоты. Следует также отметить, что если W2TH (ы) является квадратным, то он должен быть минимальной фазой для сектора, который должен удерживаться.

Эта характеристика частотной области является основой для sectorplot. В частности, sectorplot строит график сингулярных значений (W1TH (джом)) (W2TH (джом)) -1 как функции частоты. Граница сектора удовлетворяется тогда и только тогда, когда наибольшее сингулярное значение остается ниже 1. Кроме того, график содержит полезную информацию о полосах частот, где соблюдается или нарушается граница сектора, и степени, в которой она удовлетворяется или нарушается.

Например, рассмотрим график сектора системы с 2 выходами и 2 входами для конкретного сектора.

rng(4,'twister');
H = rss(3,4,2); 
Q = [-5.12   2.16  -2.04   2.17
      2.16  -1.22  -0.28  -1.11
     -2.04  -0.28  -3.35   0.00
      2.17  -1.11   0.00   0.18];
sectorplot(H,Q)

Figure contains an axes. The axes contains 2 objects of type line. This object represents H.

На графике показано, что наибольшее сингулярное значение (W1TH ()) (W2TH ()) -1 превышает 1 ниже примерно 0,5 рад/с и в узкой полосе вокруг 3 рад/с. ПоэтомуH не удовлетворяет ограничению сектора, представленному Q.

Индекс относительного сектора

Мы можем распространить понятие индекса относительной пассивности на произвольные секторы. Пусть H (s) - система LTI, и пусть:

Q = W1TW1-W2TW2, W1TW2 = 0

быть ортогональным разложением Q на его положительную и отрицательную части, как легко получается из разложения Шура Q. Относительный индекс сектора R, или R-индекс, определяется как наименьший r > 0, так что для всех выходных траекторий z (t) = (Hu) (t):

∫0TzT (t) (W1TW1-r2W2TW2) z (t) dt<0,∀T>0.

Поскольку увеличение r делает W1TW1-r2W2TW2 более отрицательным, неравенство обычно удовлетворяется для r достаточно большим. Однако существуют случаи, когда он никогда не может быть удовлетворен, и в этом случае R-индекс является R=+∞. Очевидно, что первоначальный сектор будет удовлетворен, если и только R≤1.

Чтобы понять геометрическую интерпретацию R-индекса, рассмотрим семейство конусов с матрицей Q (r) = W1TW1-r2W2TW2. В 2D, угол наклона конуса λ соотносится с r на

tan (start) =r‖W2 ‖ W1 ‖

(см. схему ниже). Таким образом, учитывая конический сектор с матрицей Q, значение R-индекса R < 1 означает, что мы можем уменьшить tan (start) (сузить конус) на коэффициент R, прежде чем некоторая выходная траектория H покинет конический сектор. Точно так же, значение R > 1 означает, что мы должны увеличить tan (start) (расширить конус) на коэффициент R, чтобы включить все выходные траектории Н. Это ясно делает R-индекс относительной мерой того, насколько хорошо отклик Н вписывается в конкретный конический сектор.

На схеме

d1|W1Tz|‖W1 , d2|W2Tz|‖W2 , R = | W1Tz | | W2Tz |,

и

tan (start) =d1d2=R‖W2 W1 ‖.

Когда W2TH (ы) является квадратной и минимальной фазой, R-индекс также может быть охарактеризован в частотной области как наименьший r > 0, так что:

H () H (W1TW1-r2W2TW2) H () <0∀ω∈R.

Используя элементарную алгебру, это приводит к:

R = maxλ (W1TH ()) (W2TH (jλ)) -1 ‖.

Другими словами, R-индекс является пиковым усилением (стабильной) передаточной функции (s): = (W1TH (s)) (W2TH (s)) -1, и сингулярные значения (jw) могут рассматриваться как «основные» R-индексы на каждой частоте. Это также объясняет, почему график зависимости R-индекса от частоты выглядит как график сингулярного значения (см.sectorplot). Существует полная аналогия между относительным индексом сектора и системным усилением. Однако следует отметить, что эта аналогия сохраняется только тогда, когда W2TH (ы) является квадратной и минимальной фазой.

Индекс направленного сектора

Аналогично, мы можем распространить понятие индекса направленной пассивности на произвольные секторы. Учитывая конический сектор с матрицей Q и направлением δQ, индекс направленного сектора является наибольшим, так что для всех выходных траекторий z (t) = (Hu) (t):

∫0TzT (t) (Q + δQ) z (t) dt<0,∀T>0.

Индекс направленной пассивности для системы G (систем) соответствует :

H (s) = (G (s) I), Q = (0-I-I0).

Индекс направленного сектора измеряет, насколько мы должны деформировать сектор в направлении δQ, чтобы он плотно вписывался вокруг выходных траекторий Н. Граница сектора удовлетворяется тогда и только тогда, когда индекс направления является положительным.

Общие секторы

Существует множество способов определения границ секторов. Затем мы рассмотрим часто встречающиеся выражения и дадим соответствующую систему H и матрицу сектора Q для стандартной формы, используемой getSectorIndex и sectorplot:

∫0T (Hu) (t) TQ (Hu) (t) dt<0,∀T>0.

Для простоты эти описания используют нотацию:

x‖T=∫0T‖x (t) 2dt,

и опустить ∀T>0 требование.

Пассивность

Пассивность - это сектор, связанный с:

H (s) = (G (s) I), Q = (0-I-I0).

Ограничение усиления

Ограничение усиления G‖∞<γ является сектором, связанным с :

H (s) = (G (s) I), Q = (I00-γ 2I).

Соотношение расстояний

Рассмотрим «внутреннее» ограничение,

‖y-cu‖T <r‖u‖T

где c, r - скаляры и y (t) = (Gu) (t). Этот сектор связан с:

H (s) = (G (s) I), Q = (I-cI-cI (c2-r2) I).

Нижележащий конический сектор симметричен относительно y = cu. Аналогично, «внешнее» ограничение,

‖y-cu‖T> r‖u‖T

является сектором, связанным с:

H (s) = (G (s) I), Q = (-IcIcI (r2-c2) I).

Двойное неравенство

При работе со статическими нелинейностями принято учитывать конические секторы вида

au2 < yu < bu2,

где y = (u) - выход нелинейности. Хотя эта взаимосвязь сама по себе не связана с сектором, она ясно подразумевает:

a∫0Tu (t) 2dt<∫0Ty (t) u (t) dt<b∫0Tu (t) 2dt

вдоль всех траекторий ввода-вывода и для всех T > 0. Это условие, в свою очередь, эквивалентно сектору, связанному с:

H (s) = (start(.) 1), Q = (1- (a + b )/2- (a + b )/2ab).

Форма продукта

Обобщённые секторные границы вида:

∫0T (y (t) -K1u (t)) T (y (t) -K2u (t)) dt < 0

соответствуют:

H (s) = (G (s) I), Q = (2I- (K2 + K1T) - (K1 + K2T) K1TK2 + K2TK1).

Как и прежде, статический сектор связан:

(y-K1u) Т (y-K2u) < 0

подразумевает интегральный сектор, связанный выше.

Диссипативный QSR

Система y = Gu является QSR-диссипативной, если она удовлетворяет:

∫0T (y (t) u (t)) T (QSSTR) (y (t) u (t)) dt>0,∀T>0.

Этот сектор связан с:

H (s) = (G (s) I), Q = - (QSSTR).

См. также

| |

Связанные темы