DynamicSystems - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Science and Engineering : Dynamic Systems : System Manipulation Tools : DynamicSystems/NormH2

DynamicSystems

  

NormH2

  

Compute the H2 norm of a linear system

 

Calling Sequence

Parameters

Options

Description

Examples

Compatibility

Calling Sequence

NormH2(sys)

Parameters

sys

-

System; system object

opts

-

(optional) equation(s) of the form option = value; specify options for the NormH2 command

Options

• 

checkstability = truefalse

True means check whether the system is stable; if it is not stable, an error occurs. False means skip the check. The default is true.

Description

• 

The NormH2 command computes the H2 norm of a linear system sys. Both continuous-time and discrete-time systems, and both single-input single-output (SISO) and multiple-input multiple-output (MIMO) systems are supported.

Continuous-time

• 

For a stable SISO linear system with transfer function Hs, the H2 norm is defined in the frequency domain as:

  

‖H‖2=12∫Hjω2ⅆωπ 

• 

For a MIMO linear system with transfer function Matrix Hs, the definition of H2 norm in the frequency domain is generalized to:

  

‖H‖2=12∫TraceHjω%H.Hjωⅆωπ

  

where A%H is the Hermitian transpose of Matrix A.

• 

In the time domain, the H2 norm of a transfer function is calculated assuming that the stable transfer function Hs has a state-space representation:

  

x.=Ax+Bw

  

y=Cx

  

so that Hs=YsWs and Hs=C. sIA1. B.

  

where the feedforward matrix D=0 is necessary for the H2 norm to be finite. It follows that, for non-strictly-causal continuous-time linear time-invariant (LTI) systems (D0), the H2 norm is infinite.

  

From the above definitions, it can be demonstrated that the H2 norm of a continuous-time LTI is equivalent to:

  

‖H‖2=TraceC.P.CT

  

where the Matrix P0 is calculated by solving a continuous Lyapunov equation:

  

A.P+P.AT+B.BT=0 

Discrete-time

• 

In the frequency domain, the H2 norm of a discrete-time LTI system is defined by:

  

‖H‖2=12∫TraceHⅇjω%H.Hⅇjωⅆωπ

  

where A%H is the Hermitian transpose of Matrix A.

• 

In the time domain, the H2 norm of a transfer function is calculated assuming that the stable transfer function Hz has a state-space representation:

  

xk+1=Axk+Bwk

  

yk=Cxk+Dwk

  

so that Hz=C. zIA1. B+D.

  

From the above definitions, it can be demonstrated that the H2 norm of a discrete-time LTI is equivalent to:

  

‖H‖2=TraceC.P.CT+D.DT

  

where the Matrix P0 is calculated by solving a discrete Lyapunov equation:

  

A.P.ATP+B.BT=0

• 

For both continuous and discrete-time systems, the H2 norm is finite if the LTI system is asymptotically stable. It follows that for unstable systems, the H2 norm is infinite.

• 

A deterministic interpretation of the H2 norm is that it measures the energy of the impulse response of the LTI system.

• 

A stochastic interpretation of the H2 norm is that it measures the energy of the output response to unit white Gaussian noise inputs. A white noise process wt has an expected or mean value 𝔼wt=0 and covariance matrix 𝔼wt.wt+τT=𝕀.δτ, where 𝕀 is the Identity Matrix and δ is the Dirac delta function. It follows that the H2 norm is equivalent to: ‖H‖2=TraceCovariance'sys',𝕀 from the interpretation above and DynamicSystems[Covariance].

Examples

withDynamicSystems:

Example 1 : Find the H2 norm of a system with discrete-time transfer function shown below.

sys1TransferFunction102z+110z2+2z+5,discrete,sampletime=0.1:

PrintSystemsys1

Transfer Functiondiscrete; sampletime = .11 output(s); 1 input(s)inputvariable=u1zoutputvariable=y1ztf1,1=20z+1010z2+2z+5

(1)

h2norm1NormH2sys1

h2norm1:=2.46238673166698

(2)

Example 2 : Find the H2 norm of a continuous state-space MIMO system.

sys2StateSpace5,3|3,4,2,3|1,1,1,2|12,1,0,0|0,0:

PrintSystemsys2

State Spacecontinuous2 output(s); 2 input(s); 2 state(s)inputvariable=u1t,u2toutputvariable=y1t,y2tstatevariable=x1t,x2ta=−533−4b=2131c=112−21d=0000

(3)

h2norm2NormH2sys2

h2norm2:=2.52637601270590

(4)

Example 3 : Find the H2 norm of the following discrete system.

sys3Coefficients1,2.841,2.875,1.004,1,2.417,2.003,0.5488,discrete,sampletime=0.1:

PrintSystemsys3

Coefficientsdiscrete; sampletime = .11 output(s); 1 input(s)inputvariable=u1zoutputvariable=y1znum1,1=1,−2.841,2.875,−1.004den1,1=1,−2.417,2.003,−0.5488

(5)

h2norm3NormH2sys3

h2norm3:=1.24382062647607

(6)

Example 4: Find the H2 norm of the system given by the following differential equation.

sys4DiffEquationⅆⅆtⅆⅆtxt=10xtⅆⅆtxt+wt,wt,xt:

PrintSystemsys4

Diff. Equationcontinuous1 output(s); 1 input(s)inputvariable=wtoutputvariable=xtde=x..t=10xtx.t+wt

(7)

h2norm4NormH2sys4

h2norm4:=0.223606797749979

(8)

Example 5 : Find the H2 norm of a non-strictly-causal continuous state-space MIMO system.

sys5StateSpace5,3|3,4,2,3|1,1,1,2|12,1,2,1|3,7:

PrintSystemsys5

State Spacecontinuous2 output(s); 2 input(s); 2 state(s)inputvariable=u1t,u2toutputvariable=y1t,y2tstatevariable=x1t,x2ta=−533−4b=2131c=112−21d=2317

(9)

Since the H2 norm is infinite, an error message is displayed.

h2norm5NormH2sys5

Error, (in DynamicSystems:-NormH2) H2 norm is infinite for continuous 'sys' with D<>0 (system is not strictly causal).

Example 6: Find the H2 norm of an unstable system given by the continuous transfer function G(s).

sys6TransferFunction4s&plus;35s4&plus;7s3&plus;4s2&plus;3s&plus;1&colon;

PrintSystemsys6

Transfer Functioncontinuous1 output(s); 1 input(s)inputvariable&equals;u1soutputvariable&equals;y1stf1,1&equals;4s+35s4+7s3+4s2+3s+1

(10)

Since the H2 norm is infinite, an error message is displayed.

h2norm6NormH2sys6

Error, (in DynamicSystems:-NormH2) H2 norm is infinite for unstable systems. Unstable eigenvalues of 'sys': .324596325e-1-.6550790710*I, .324596325e-1+.6550790710*I

Compatibility

• 

The DynamicSystems[NormH2] command was introduced in Maple 18.

• 

For more information on Maple 18 changes, see Updates in Maple 18.

See Also

DynamicSystems

DynamicSystems[Covariance]

DynamicSystems[Grammians]

LinearAlgebra[HermitianTranspose]

LinearAlgebra[LyapunovSolve]

LinearAlgebra[SylvesterSolve]

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam