受動システムのフィードバック相互接続
この例では、受動システムのフィードバック相互接続のプロパティについて説明します。
受動システムのフィードバック相互接続
2 つのサブシステム と のフィードバックでの相互接続について考えます。相互接続されたシステム は入力 を出力 にマッピングします。
システム と がともに受動的である場合、相互接続されたシステム は受動的となることが保証されます。たとえば次の場合
次で確認されるとおり、両方のシステムは受動的です。
G1 = tf([1,1,1],[1,1,4]); isPassive(G1)
ans = logical
1
G2 = tf([1,2],[1,5]); isPassive(G2)
ans = logical
1
したがって、相互接続されたシステムは受動的です。
H = feedback(G1,G2); isPassive(H)
ans = logical
1
これは、 のナイキスト線図が正の実数であることを検証して確認されます。
nyquist(H)
フィードバック相互接続の受動性インデックス
と の受動性インデックスと、相互接続されたシステム の受動性インデックスとの間には関係性があります。 と で と の入力受動性インデックスを表し、 と で出力受動性インデックスを表します。すべてのインデックスが正の場合、フィードバック相互接続 の入力受動性インデックス と出力受動性インデックス は次を満たします。
つまり、 と の入出力受動性インデックスから、閉ループ システム の最小レベルの入出力受動性を一部推測できます。詳細については、Zhu, F.、Xia, M および Antsaklis, P.J. による論文 "Passivity analysis and passivation of feedback systems using passivity indices," American Control Conference , 2014, pp. 1833-1838. を参照してください。入力受動性インデックス の下限を確認します。
% Input passivity index for G1 nu1 = getPassiveIndex(G1,'input'); % Output passivity index for G2 rho2 = getPassiveIndex(G2,'output'); % Input passivity index for H nu = getPassiveIndex(H,'input')
nu = 0.1293
% Lower bound
nu1*rho2/(nu1+rho2)
ans = 7.1402e-11
同様に、 の出力受動性インデックスの下限を確認します。
% Output passivity index for G1 rho1 = getPassiveIndex(G1,'output'); % Input passivity index for G2 nu2 = getPassiveIndex(G2,'input'); % Output passivity index for H rho = getPassiveIndex(H,'output')
rho = 0.4441
% Lower bound
rho1+nu2
ans = 0.4000