getSectorIndex

Вычисление индекса конического сектора линейной системы

Описание

пример

RX = getSectorIndex(H,Q) вычисляет относительный индекс RX для линейной системы H и конический сектор, заданный Q. Когда RX < 1, все выходные траектории y (t) = Hu (t) лежат в секторе, заданном:

0Ty(t)TQy(t)dt<0,

для всех T ≥ 0.

getSectorIndex можно также проверить, G ли все траектории ввода-вывода {u (t), y (t)} линейной системы лежать в секторе, заданном:

0T(y(t)u(t))TQ(y(t)u(t))dt<0,

для всех T ≥ 0. Для этого используйте getSectorIndex с H = [G;I], где I = eyes(nu), и nu - количество входов G.

Для получения дополнительной информации об границах сектора и относительном индексе смотрите О границах сектора и индексах сектора.

RX = getSectorIndex(H,Q,tol) вычисляет индекс с относительной точностью, заданной tol.

RX = getSectorIndex(H,Q,tol,fband) вычисляет индекс пассивности путем ограничения неравенств, которые определяют индекс, заданным частотным интервалом. Этот синтаксис доступен только при Q имеет столько отрицательных собственных значений, сколько входные входы в H.

[RX,FX] = getSectorIndex(___) также возвращает частоту, при которой значение индекса RX достигается. FX установлено в NaN когда количество отрицательных собственных значений в Q отличается от количества входов в H. Можно использовать этот синтаксис с любой из предыдущих комбинаций входных параметров.

[RX,FX,W1,W2,Z] = getSectorIndex(___) также возвращает разложение Q в его положительную и отрицательную части, а также спектральный фактор Z когда Q является динамическим. Когда Q является матрицей (постоянными границами сектора), Z = 1. Можно использовать этот синтаксис с любой из предыдущих комбинаций входных параметров.

DX = getSectorIndex(H,Q,dQ) вычисляет индекс в направлении, заданном матрицей dQ. Если DX > 0, затем выходные траектории H подходит в коническом секторе, заданном Q. Для получения дополнительной информации об указателе направления смотрите О границах сектора и индексах сектора.

Направленный индекс недоступен, если H является данными частотной характеристики (frd) модель.

DX = getSectorIndex(H,Q,dQ,tol) вычисляет индекс с относительной точностью, заданной tol.

Примеры

свернуть все

Проверяйте, в среднем ли траектории ввода-вывода G(s)=(s+2)/(s+1) относятся к сектору, определяемому:

S={(y,u):0.1u2<uy<10u2}.

В пространстве U/Y этот сектор является затененной областью следующей схемы.

The Q матрица, соответствующая этому сектору, задается:

Q=[1-(a+b)/2-(a+b)/2ab];a=0.1,b=10.

Траектория y(t)=Gu(t) лежит в пределах сектора S, когда для всех T > 0,

0.10Tu(t)2<0Tu(t)y(t)dt<100Tu(t)2dt.

Чтобы проверить, являются ли траектории G удовлетворить секторную границу, представленную Q, вычислите R-индекс для H = [G;1].

G = tf([1 2],[1 1]); 

a = 0.1;  b = 10;
Q = [1 -(a+b)/2 ; -(a+b)/2 a*b];

R = getSectorIndex([G;1],Q)
R = 0.4074

Это приводит к R меньше 1, что указывает на то, что траектории соответствуют сектору. Значение R подсказывает, насколько плотно траектории подгонки в секторе. Это значение, R = 0,41, означает, что траектории будут подгонки в более узком секторе с основой 1/0,41 = в 2,4 раза меньше.

Входные параметры

свернуть все

Модель для анализа против секторных границ, заданная как динамическая системная модель, такая как tf, ss, или genss модель. H может быть непрерывным или дискретным. Если H - обобщенная модель с настраиваемыми или неопределенными блоками, getSectorIndex анализирует текущее, номинальное значение H.

Чтобы проанализировать, все ли траектории ввода-вывода {u (t), y (t)} линейной системы G лежать в конкретном секторе, используйте H = [G;I].

Если H является массивом моделей, тогда getSectorIndex возвращает индекс пассивности как массив того же размера, где:

index(k) = getSectorIndex(H(:,:,k),___)

Здесь, index является либо RX, или DX, в зависимости от того, какие входные параметры вы используете.

Геометрия сектора, заданная как:

  • Матрица, для постоянной секторной геометрии. Q - симметричная квадратная матрица, которая ny со стороны, где ny - количество выходов H.

  • Модель LTI для частотно-зависимой геометрии сектора. Q удовлетворяет Q (s) "= Q (- s). Другими словами, Q (s) оценивает в эрмитову матрицу на каждой частоте.

Матрица Q должен быть неопределенным, чтобы описать четко определенный конический сектор. Неопределенная матрица имеет как положительные, так и отрицательные собственные значения.

Для получения дополнительной информации смотрите О секторных границах и секторных индексах.

Относительная точность для вычисленного индекса сектора. По умолчанию допуск равен 1%, что означает, что возвращенный индекс находится в пределах 1% от фактического индекса.

Частотный интервал для вычисления индекса сектора, заданный как массив вида [fmin,fmax]. Когда вы предоставляете fband, getSectorIndex ограничивает заданный частотный интервал неравенствами, которые определяют индекс. Задайте частоты в единицах rad/TimeUnit, где TimeUnit является TimeUnit свойство модели динамической системы H.

Направление, в котором вычислить индекс направленного сектора, заданный как неотрицательная определенная матрица. Матрица dQ - симметричная квадратная матрица, которая ny со стороны, где ny - количество выходов H.

Выходные аргументы

свернуть все

Относительный индекс системного H для сектора, заданного Q, возвращенный в виде скалярного значения или массива, если H - массив. Если RX < 1, затем выходные траектории H помещается внутри конуса Q.

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

Q=W1W1TW2W2T,W1TW2=0.

(Такое разложение легко получить из разложения Шура Q.) Затем, RX - наименьший R, который удовлетворяет:

0Ty(t)T(W1W1TR2W2W2T)y(t)dt<0,

для всех T ≥ 0. Изменение R эквивалентно регулированию угла наклона конуса, заданного Q, пока конус не будет плотно помещаться вокруг выходных траекторий H. Отношение основания конуса к высоте пропорционально R.

Для получения дополнительной информации о интерпретации относительного индекса, смотрите О границах сектора и индексах сектора.

Частота, при которой RX индекс достигается, возвращается как неотрицательный скаляр или массив, если H - массив. В целом индекс изменяется с частотой (см. sectorplot). Возвращенное значение является самым большим значением на всех частотах. FX - частота, на которой происходит это значение, возвращаемая в единицах rad/TimeUnit, где TimeUnit является TimeUnit свойство H.

Положительные и отрицательные факторы Q, возвращается как матрицы. Для постоянного Q, W1 и W2 удовлетворить:

Q=W1W1TW2W2T,W1TW2=0.

Бистабельная модель в факторизации Q, возвращается как:

  • Если Q является постоянной матрицей, Z = 1.

  • Если Q зависит от частоты, тогда Z является пространством состояний (ss) моделировать таким образом, чтобы:

Q(jω)=Z(jω)H(W1W1TW2W2T)Z(jω).

Индекс направленного сектора системы H для сектора, заданного Q в направлении dQ, возвращенный в виде скалярного значения или массива, если H - массив. Направленный индекс является самой большой τ, которая удовлетворяет:

0Ty(t)T(Q+τdQ)y(t)dt<0,

для всех T ≥ 0.

Введенный в R2016a