LinearAlgebra
DARE
solves the discrete algebraic Riccati equation
Calling Sequence
Parameters
Options
Description
Examples
DARE(A, B, Q, R, S, options, outopts)
A
-
Matrix(square)
B
Matrix
Q
R
S
(optional) Matrix
options
(optional) constructor options for the result objects and/or equation(s) of the form output = value
outopts
(optional) equation(s) of the form outputoptions[o] = list, where o is one of X, L or G
output=value
Specifies which results are returned, where value can be a single name or a list of names.. The choices of names for value are
X - the Matrix which solves the equation
L - the Vector of closed-loop eigenvalues of the symplectic Matrix H
rcond - the reciprocal of the condition number of the system solved by X
G - the gain Matrix where G=B+ ·X·B+R-1·S+A+ ·X·B+ .
The default for value is X.
The constructor options provide additional information (readonly, shape, storage, order, datatype, and attributes) to the Matrix or Vector constructor that builds the result(s). These options may also be provided in the form outputoptions[o]=[...], where [...] represents a Maple list. If a constructor option is provided in both the calling sequence directly and in an outputoptions[o] option, the latter takes precedence (regardless of the order).
The DARE command solves the discrete algebraic Riccati equation,
A+ ·X·A−X−S+A+ ·X·B·B+ ·X·B+R-1·S+A+ ·X·B+ +Q=0
The optional Matrix argument S defaults to the zero Matrix.
This routine operates in the real floating-point domain. Hence, the entries in the Matrix arguments must necessarily be of type numeric.
with⁡LinearAlgebra:
a≔Matrix⁡2,2,0,1,0,0
a≔0100
b≔Matrix⁡2,2,0.5,0.7,1,0
b≔0.50.710
q≔Matrix⁡3,2,1,0,0,1,0,0
q≔100100
r≔Matrix⁡1,1,1
r≔1
DARE⁡a,b,q,r
1.000000000000001.74309927745499×10−161.74309927745499×10−161.92116460960662
DARE⁡a,b,q,r,output=X
DARE⁡a,b,q,r,1,0|0,1
0.626691964175562−0.0851772290376543−0.08517722903765431.60725718203777
DARE⁡a,b,q,r,output=X,L,G,rcond
1.000000000000001.74309927745499×10−161.74309927745499×10−161.92116460960662,1.64673469979885×10−16+0.⁢I−0.157670780786755+0.⁢I,0.0.157670780786755,0.705446591923427513
See Also
CARE
copyright
LyapunovSolve
SylvesterSolve
Vector
Download Help Document