Индексы пассивности

Этот пример показывает, как вычислить различные меры пассивности для линейных инвариантных по времени систем.

Пассивные системы

Линейная система G (s) пассивна, когда все траектории ввода-вывода(u(t),y(t)) удовлетворить

0TyT(t)u(t)dt>0,T>0

где yT(t) обозначает транспонирование y(t).

Чтобы измерить «насколько пассивна» система, мы используем индексы пассивности.

  • Входной индекс пассивности определяется как самый большой ν таким, что

0TyT(t)u(t)dt>ν0TuT(t)u(t))dt

Система G является «входом строго пассивным» (ISP), когда ν>0. ν называется также индексом «входной пассивности с feedforward» (IFP) и соответствует минимальному действию с feedforward, необходимому для того, чтобы сделать систему пассивной.

  • Индекс пассивности выхода определяется как самый большой ρ таким, что

0T(yT(t)u(t)dt>ρ0TyT(t)y(t))dt

Система G «выводит строго пассивно» (OSP), когда ρ>0. ρ называется также индексом «пассивности обратной связи выхода» (OFP) и соответствует минимальному действию обратной связи, необходимому для того, чтобы сделать систему пассивной.

  • Индекс пассивности ввода-вывода определяется как самый большой τ таким, что

0TyT(t)u(t)dt>τ0T(uT(t)u(t)+yT(t)y(t))dt

Система является «очень строго пассивной» (VSP), если τ>0.

Пример схемы

Рассмотрим следующий пример. Берем ток I как вход и напряжение V в качестве выхода. Основываясь на законе тока и напряжения Кирхгофа, мы получаем передаточную функцию для G(s),

G(s)=V(s)I(s)=(Ls+R)(Rs+1C)Ls2+2Rs+1C.

Давайте R=2, L=1 и C=0.1.

R = 2; L = 1; C = 0.1; 
s = tf('s');
G = (L*s+R)*(R*s+1/C)/(L*s^2 + 2*R*s+1/C);

Использование isPassive чтобы проверить, является ли G(s) является пассивным.

PF = isPassive(G)
PF = logical
   1

Поскольку PF = true, G(s) является пассивным. Использование getPassiveIndex для вычисления индексов пассивности G(s).

% Input passivity index
nu = getPassiveIndex(G,'in')
nu = 2
% Output passivity index
rho = getPassiveIndex(G,'out')
rho = 0.2857
% I/O passivity index
tau = getPassiveIndex(G,'io')
tau = 0.2642

С тех пор τ>0, систему G(s) очень строго пассивно.

Характеристика частотного диапазона

Линейная система пассивна тогда и только тогда, когда она «положительно реально»:

G(jω)+GH(jω)>0ωR.

Наименьшее собственное значение левой стороны связано с входом индексом пассивности ν:

ν=12minωλmin(G(jω)+GH(jω))

где λmin обозначает наименьшее собственное значение. Точно так же, когда G(s) является минимальной фазой, выход индекс пассивности задается:

ρ=12minωλmin(G-1(jω)+G-H(jω)).

Проверьте это для примера схемы. Постройте график Годографа Найквиста передаточной функции схемы.

nyquist(G)

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

Весь годограф Найквиста лежит в правой-половинной плоскости так G(s) положительно реально. Крайняя левая точка кривой Найквиста (x,y)=(2,0) поэтому входной индекс пассивности ν=2, то же значение, которое мы получили ранее. Точно так же самая левая точка на кривой Найквиста для G-1(s) задает выход значение индекса пассивности ρ=0.286.

Относительный индекс пассивности

Можно показать, что «положительное реальное» условие

G(jω)+GH(jω)>0ωR

эквивалентно малому условию усиления

||(I-G(jω))(I+G(jω))-1||<1ωR.

Относительный индекс пассивности (R-индекс) является пиковым коэффициентом усиления по частоте (I-G)(I+G)-1 когда I+G является минимальной фазой, и + в противном случае:

R=(I-G)(I+G)-1.

Во временном интервале R-индекс является наименьшим r>0 таким, что

0T||y-u||2dt<r20T||y+u||2dt

Система G(s) пассивен тогда и только тогда, когда R<1, и меньшее R тем более пассивна система. Использование getPassiveIndex для вычисления R-индекса для примера схемы.

R = getPassiveIndex(G)
R = 0.5556

Получающееся R значение указывает, что схема является очень пассивной системой.

См. также

|

Похожие темы