Этот пример иллюстрирует свойства последовательного соединения пассивных систем.
Рассмотрим взаимосвязь двух подсистем и последовательно. Взаимосвязанная система задается отображением от входа производить .
В отличие от параллельных соединений и соединений с обратной связью, пассивность подсистем и не гарантирует пассивность для взаимосвязанной системы . Примите за пример
Обе системы пассивны, как подтверждено
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
Однако последовательное соединение и не пассивен:
H = G2*G1; isPassive(H)
ans = logical
0
Это подтверждается проверкой, что Годограф Найквиста не является положительным реальным.
nyquist(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
Точно так же проверьте нижнюю границу для выхода индекса пассивности .
% 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