Parallel Interconnection of Passive Systems

Consider an interconnection of two subsystems $$G_1$$and $$G_2$$in parallel. The interconnected system $$H$$maps the input $$r$$到输出 $$y$$.

If both systems $$G_1$$and $$G_2$$are passive, then the interconnected system $$H$$is guaranteed to be passive. Take for example

Both systems are passive.

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

ans =logical1

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

ans =logical1

We can therefore expect their parallel interconnection $$H$$to be passive, as confirmed by

H = parallel(G1,G2); isPassive(H)

ans =logical1

Passivity Indices for Parallel Interconnection

There is a relationship between the passivity indices of $$G_1$$and $$G_2$$and the passivity indices of the interconnected system $$H$$. Let $${\nu }_1$$and $${\nu }_2$$denote the input passivity indices for $$G_1$$and $$G_2$$, and let $${\rho }_1$$and $${\rho }_2$$denote the output passivity indices. If all these indices are nonnegative, then the input passivity index $$\nu$$and the output passivity index $$\rho$$for the parallel interconnection $$H$$满足

In other words, we can infer some minimum level of input and output passivity for the parallel connection $$H$$from the input and output passivity indices of $$G_1$$and $$G_2$$. For details, see the paper by Yu, H., "Passivity and dissipativity as design and analysis tools for networked control systems,"Chapter 2, PhD Thesis, University of Notre Dame, 2012. Verify the lower bound for the input passivity index $$\nu$$.

% Input passivity index for G1nu1 = getPassiveIndex(G1,'input');% Input passivity index for G2nu2 = getPassiveIndex(G2,'input');% Input passivity index for Hnu = getPassiveIndex(H,'input')

nu = 0.3777

% Lower boundnu1+nu2

ans = 0.1474

Similarly, verify the lower bound for the output passivity index of $$H$$.

% Output passivity index for G1rho1 = getPassiveIndex(G1,'output');% Output passivity index for G2rho2 = getPassiveIndex(G2,'output');% Output passivity index for Hrho = getPassiveIndex(H,'output')

rho = 0.6443

% Lower boundrho1*rho2/(rho1+rho2)

ans = 0.2098