Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
5.3.1 - English

Change language to:
Français - 日本語 - Português

Please note that the recommended version of Scilab is 2024.0.0. This page might be outdated.
However, this page did not exist in the previous stable version.

Scilab help >> CACSD > ric_desc

ric_desc

Riccati equation

Calling Sequence

X=ric_desc(H [,E))
[X1,X2,zero]=ric_desc(H [,E])

Arguments

H,E

real square matrices

X1,X2

real square matrices

zero

real number

Description

Riccati solver with hamiltonian matrices as inputs.

In the continuous time case calling sequence is ric_descr(H) (one input):

Riccati equation is:

(Ec)   A'*X + X*A + X*R*X -Q = 0.

Defining the hamiltonian matrix H by:

H = [A  R;
     Q -A']

with the calling sequence [X1,X2,zero]=ric_descr(H), the solution X is given by X=X1/X2.

zero = L1 norm of rhs of (Ec)

The solution X is also given by X=riccati(A,Q,R,'c'))

In the discrete-time case calling sequence is ric_descr(H,E) (two inputs):

The Riccati equation is:

(Ed)  A'*X*A-(A'*X*B*(R+B'*X*B)^-1)*(B'*X*A)+C-X = 0.

Defining G=B/R*B' and the hamiltonian pencil (E,H) by:

E=[eye(n,n),G;               H=[A, 0*ones(n,n);
   0*ones(n,n),A']             -C, eye(n,n)];

with the calling sequence [X1,X2,err]=ric_descr(H,E), the solution X is given by X=X1/X2.

zero= L1 norm of rhs of (Ed)

The solution X is also given by X=riccati(A,G,C,'d') with G=B/R*B'

See Also

<< repfreq CACSD ricc >>

Copyright (c) 2022-2023 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Mar 03 10:59:41 CET 2011