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

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

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

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

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

G1(s)=0.1s+1s+2;G2(s)=s2+2s+1s2+3s+10

Обе системы пассивны.

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

G2 = tf([1,2,1],[1,3,10]);
isPassive(G2)
ans = logical
   1

Поэтому мы можем ожидать их параллельного соединения H быть пассивным, как подтверждено

H = parallel(G1,G2);
isPassive(H)
ans = logical
   1

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

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

νν1+ν2,ρρ1ρ2ρ1+ρ2.

Другими словами, мы можем вывести некоторый минимальный уровень входа и вывести пассивность для параллельного соединения H из входа и выхода индексов пассивности G1 и G2. Для получения дополнительной информации смотрите статью Ю, Х. «Пассивность и рассеянность как инструменты проекта и анализа для сетевых систем управления», Глава 2, PhD Thesis, Университет Нотр-Дам, 2012. Проверьте нижнюю границу для входного индекса пассивности ν.

% Input passivity index for G1
nu1 = getPassiveIndex(G1,'input');
% Input passivity index for G2
nu2 = getPassiveIndex(G2,'input');
% Input passivity index for H
nu = getPassiveIndex(H,'input')
nu = 0.3777
% Lower bound
nu1+nu2
ans = 0.1474

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

% Output passivity index for G1
rho1 = getPassiveIndex(G1,'output');
% Output passivity index for G2
rho2 = getPassiveIndex(G2,'output');
% Output passivity index for H
rho = getPassiveIndex(H,'output')
rho = 0.6443
% Lower bound
rho1*rho2/(rho1+rho2)
ans = 0.2098

См. также

|

Похожие темы

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