Последовательное соединение пассивных систем

Этот пример иллюстрирует свойства последовательного соединения пассивных систем.

Последовательное соединение пассивных систем

Рассмотрим взаимосвязь двух подсистем G1 и G2 последовательно. Взаимосвязанная система H задается отображением от входа u производить y2.

В отличие от параллельных соединений и соединений с обратной связью, пассивность подсистем G1 и G2 не гарантирует пассивность для взаимосвязанной системы H. Примите за пример

G1(s)=5s2+3s+1s2+2s+1,G2(s)=s2+s+5s+0.1s3+2s2+3s+4.

Обе системы пассивны, как подтверждено

G1 = tf([5 3 1],[1,2,1]);
isPassive(G1)
ans = logical
   1

G2 = tf([1,1,5,.1],[1,2,3,4]);
isPassive(G2)
ans = logical
   1

Однако последовательное соединение G1 и G2 не пассивен:

H = G2*G1;
isPassive(H)
ans = logical
   0

Это подтверждается проверкой, что Годограф Найквиста G2G1 не является положительным реальным.

nyquist(H)

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

Индексы пассивности для последовательного соединения

В то время как последовательное соединение пассивных систем не является пассивным в целом, существует связь между индексами пассивности G1 и G2 и индексы пассивности H=G2G1. Давайте ν1 и ν2 обозначить входные индексы пассивности для G1 и G2, и пусть ρ1 и ρ2 обозначить выходные индексы пассивности. Если все эти индексы положительны, то входной индекс пассивности ν и индекс выхода пассивности ρ для последовательного соединения H удовлетворить

ν-0.125ρ1ρ2,ρ-0.125ν1ν2.

Другими словами, нехватка пассивности на входах или выходах H не хуже, чем выражения правой стороны. Для получения дополнительной информации см. статью Arcak, M. and Sontag, E.D., «Диагональная устойчивость класса циклических систем и его связь с секретным критерием», Automatica, Vol 42, № 9, 2006, pp. 1531-1537. Проверьте эти нижние границы для примера выше.

% Output passivity index for G1
rho1 = getPassiveIndex(G1,'output');
% Output passivity index for G2
rho2 = getPassiveIndex(G2,'output');
% Input passivity index for H=G2*G1
nu = getPassiveIndex(H,'input')
nu = -1.2886
% Lower bound
-0.125/(rho1*rho2)
ans = -2.4194

Точно так же проверьте нижнюю границу для выхода индекса пассивности H.

% Input passivity index for G1
nu1 = getPassiveIndex(G1,'input');
% Input passivity index for G2
nu2 = getPassiveIndex(G2,'input');
% Output passivity index for H=G2*G1
rho = getPassiveIndex(H,'output')
rho = -0.6966
% Lower bound
-0.125/(nu1*nu2)
ans = -6.0000

См. также

|

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте