getPassiveIndex

Compute passivity index of linear system

collapse all in page

Syntax

R = getPassiveIndex(G)
nu = getPassiveIndex(G,'input')
rho = getPassiveIndex(G,'output')
tau = getPassiveIndex(G,'io')
DX = getPassiveIndex(G,dQ)
index = getPassiveIndex(___,tol)
index = getPassiveIndex(___,tol,fband)
[index,FI] = getPassiveIndex(___)
[index,FI,Qout,dQout] = getPassiveIndex(___)

Description

getPassiveIndex computes various measures of the excess or shortage of passivity for a given system.

A linear system G(s) is passive if all its I/O trajectories (u(t),y(t)) satisfy:

∫0Ty(t)Tu(t)dt>0,

for all T > 0. Equivalently, a system is passive if its frequency response is positive real, such that for all ω > 0,

G(jω)+G(jω)H>0

(or the discrete-time equivalent).

example

R = getPassiveIndex(G) computes the relative passivity index. G is passive when R is less than one. R measures the relative excess (R < 1) or shortage (R > 1) of passivity.

For more information about the notion of passivity indices, see About Passivity and Passivity Indices.

example

nu = getPassiveIndex(G,'input') computes the input passivity index. The system is input strictly passive when nu > 0. nu is also called the input feedforward passivity (IFP) index. The value of nu is the minimum feedforward action such that the resulting system is passive.

For more information about the notion of passivity indices, see About Passivity and Passivity Indices.

example

rho = getPassiveIndex(G,'output') computes the output passivity index. The system is output strictly passive when rho > 0. rho is also called the output feedback passivity (OFP) index. The value of rho is the minimum feedback action such that the resulting system is passive.

For more information about the notion of passivity indices, see About Passivity and Passivity Indices.

example

tau = getPassiveIndex(G,'io') computes the combined I/O passivity index. The system is very strictly passive when tau > 0.

For more information about the notion of passivity indices, see About Passivity and Passivity Indices.

DX = getPassiveIndex(G,dQ) computes the directional passivity index in the direction specified by the matrix dQ.

index = getPassiveIndex(___,tol) computes the passivity index with relative accuracy specified by tol. Use this syntax with any of the previous combinations of input arguments. index is the corresponding passivity index R, nu, rho, tau, or DX.

index = getPassiveIndex(___,tol,fband) computes passivity indices restricted to a specific frequency interval.

[index,FI] = getPassiveIndex(___) also returns the frequency at which the returned index value is achieved.

[index,FI,Qout,dQout] = getPassiveIndex(___) also returns the sector matrix Qout for passivity and the directional index matrix dQout.

Examples

collapse all

Open Live Script

Compute passivity indices for the following dynamic system:

G(s)=s2+s+5s+0.1s3+2s2+3s+4

G = tf([1,1,5,.1],[1,2,3,4]);

Compute the relative passivity index.

R = getPassiveIndex(G)
R = 0.9512

The system is passive, but with a relatively small excess of passivity.

Compute the input and output passivity indices.

nu = getPassiveIndex(G,'input')
nu = 0.0250

rho = getPassiveIndex(G,'output')
rho = 0.2583

This system is both input strictly passive and output strictly passive.

Compute the combined I/O passivity index.

tau = getPassiveIndex(G,'io')
tau = 0.0250

The system is very strictly passive as well. A system that is very strictly passive is also strictly positive real. Examining the Nyquist plot confirms this, showing that the frequency response lies entirely within the right half-plane. 

nyquistplot(G)

The relatively small tau value is reflected in how close the frequency response comes to the imaginary axis.

Input Arguments

collapse all

Model to analyze for passivity, specified as a dynamic system model such as a tf, ss, or genss model. G can be MIMO, if the number of inputs equals the number of outputs. G can be continuous or discrete. If G is a generalized model with tunable or uncertain blocks, getPassiveIndex evaluates passivity of the current, nominal value of G.

If G is a model array, then getPassiveIndex returns the passivity index as an array of the same size, where:

index(k) = getPassivityIndex(G(:,:,k),___)

Here, index is any of R, nu, rho, tau, or DX, depending on which input arguments you use.

Custom direction in which to compute passivity, specified as a symmetric square matrix that is 2*ny on a side, where ny is the number of outputs of G.

The rho, nu, and tau indices each correspond to a particular direction in the y/u space of the system, with a corresponding dQ value. (See dQout for these values.) Use this argument to specify your own value for this direction.

Relative accuracy for the calculated passivity index. By default, the tolerance is 1%, meaning that the returned passivity index is within 1% of the actual passivity index.

Frequency interval for determining passivity index, specified as an array of the form [fmin,fmax]. When you provide fband, then getPassiveIndex restricts the frequency-domain computation of the passivity index to that frequency range. For example, the relative passivity index R is the peak gain of the bilinear-transformed system (I-G)(I+G)^-1 (for minimum-phase (I + G)). When you provide fband, then R is the peak gain within the frequency band.

Specify frequencies in units of rad/TimeUnit, where TimeUnit is the TimeUnit property of the dynamic system model G.

Output Arguments

collapse all

Relative passivity index, returned as a scalar, or an array if G is an array.

The system G is passive when R is less than one.

  • R < 1 indicates a relative excess of passivity.

  • R > 1 indicates a relative shortage of passivity.

When I + G is minimum phase, R is the peak gain of the bilinear-transformed system (I - G)(I + G)^-1.

For more information about the notion of passivity indices, see About Passivity and Passivity Indices.

Input passivity index, returned as a scalar, or an array if G is an array. nu is defined as the largest value of ν for which:

∫0Ty(t)Tu(t)dt>ν∫0Tu(t)Tu(t)dt,

for all T > 0. Equivalently, nu is the largest ν for which:

G(jω)+G(jω)H>2νI

(or the discrete-time equivalent). The system is input strictly passive when nu > 0. nu is also called the input feedforward passivity (IFP) index. The value of nu is the minimum feedforward action such that the resulting system is passive.

Output passivity index, returned as a scalar, or an array if G is an array. rho is defined as the largest value of ρ for which:

∫0Ty(t)Tu(t)dt>ρ∫0Ty(t)Ty(t)dt,

for all T > 0. The system is output strictly passive when rho > 0. rho is also called the output feedback passivity (OFP) index. The value of rho is the minimum feedback action such that the resulting system is passive.

Combined I/O passivity index, returned as a scalar, or an array if G is an array. tau is defined as the largest value of τ for which:

∫0Ty(t)Tu(t)dt>τ∫0T(u(t)Tu(t)+y(t)Ty(t))dt,

for all T > 0. The system is very strictly passive when tau > 0.

Directional passivity index in the direction specified by dQ, returned as a scalar, or an array if G is an array. The directional passivity index is the largest value of D for which:

∫0Ty(t)Tu(t)dt>D∫0T((y(t)u(t))TdQ(y(t)u(t)))dt,

for all T > 0. The rho, nu, and tau indices correspond to particular choices of dQ (see the output argument dQout). To compute DX, the software uses the custom dQ value you supply, dQ.

Frequency at which the returned passivity index is achieved, returned as a nonnegative scalar, or an array if G is an array. In general, the passivity indices vary with frequency (see passiveplot). For each index type, the returned value is the largest value over all frequencies. FI is the frequency at which this value occurs, returned in units of rad/TimeUnit, where TimeUnit is the TimeUnit property of G.

Sector geometry used for computing the passivity index, returned as a matrix. For passivity indices, Qout is given by:

Qout = [zeros(ny),-1/2*eye(ny);-1/2*eye(ny),zeros(ny)];

where ny is the number of outputs of G. For example, for a SISO G,

Qout = [ 0,  -0.5;
     -0.5, 0   ];

For more information about sector geometry, see getSectorIndex.

Direction in which passivity is computed, returned as a square matrix that is 2*ny on a side, where ny is the number of outputs of G. The value returned for dQout depends on what kind of passivity index you calculate:

  • nu — For the input passivity index, dQout is given by:

    dQout = [zeros(ny),zeros(ny);zeros(ny),eye(ny)];

    For instance, for a SISO system, dQout = [0,0;0,1].

  • rho — For the output passivity index, dQout is given by:

    dQout = [eye(ny),zeros(ny);zeros(ny),zeros(ny)];

    For instance, for a SISO system, dQout = [1,0;0,0].

  • tau — For the combined I/O passivity index, dQout is given by:

    dQout = eye(2*ny);

    For instance, for a SISO system, dQout = [1,0;0,1].

  • DX — dQout is the custom value you provide in the dQ input argument.

  • R — The relative passivity index does not involve a direction, so in this case the function returns dQout = [].

For more information about directional indices, see getSectorIndex.

See Also

| | | | |

Topics

Introduced in R2016a

Control System Toolbox Documentation

Support