exponenta event banner

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

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

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

Линейная система 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. , также называется индексом «входной пассивности передачи» (IFP) и соответствует минимальному действию передачи, необходимому для того, чтобы сделать систему пассивной.

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

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

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

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

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

Система является «очень строго пассивной» (VSP), если start> 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 (ы) является пассивным. Использовать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

Поскольку start> 0, система G (ы) очень строго пассивна.

Характеристика частотной области

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

G () + GH () >0∀ω∈R.

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

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

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

(G-1 () + 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), так что индекс пассивности на входе равен Аналогично, крайняя левая точка на кривой Найквиста для G-1 (ей) дает выходное значение индекса пассивности

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

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

G () + GH () >0∀ω∈R

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

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

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

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

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

∫0T||y-u||2dt<r2∫0T||y+u||2dt

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

R = getPassiveIndex(G)
R = 0.5556

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

См. также

|

Связанные темы