dlyap (Control System Toolbox)

Solve discrete-time Lyapunov equations

Syntax

X = dlyap(A,Q)
X = dlyap(A,B,C)
X = dlyap(A,Q,[],E)

Description

X = dlyap(A,Q) solves the discrete-time Lyapunov equation

where and are -by- matrices.

The solution is symmetric when is symmetric, and positive definite when is positive definite and has all its eigenvalues inside the unit disk.

X = dlyap(A,B,C) solves the Sylvester equation

where A, B, and C must have compatible dimensions but need not be square.

X = dlyap(A,Q,[],E) solves the generalized discrete-time Lyapunov equation

where Q is a symmetric matrix. The empty square brackets, [], are mandatory. If you place any values inside them, the function will error out.

Algorithm

dlyap uses SLICOT routines SB03MD and SG03AD for Lyapunov equations and SB04QD (SLICOT) for Sylvester equations.

Diagnostics

The discrete-time Lyapunov equation has a (unique) solution if the eigenvalues of satisfy for all .

If this condition is violated, dlyap produces the error message

Solution does not exist or is not unique.

See Also
covar Covariance of system response to white noise

lyap Solve continuous Lyapunov equations

dlqr dlyapchol