Calculate the upper bound of the CARE, i.e. Q + A'*X + X*A - X*G*X = 0 condition number using Lyapunov equations
function k_care_u extends Modelica.Icons.Function; import Modelica.Math.Matrices.norm; import Modelica_LinearSystems2.Math.Matrices.lyapunov; input Real A[:, size(A, 1)] "care-matrix A"; input Real Q[:, size(Q, 1)] "care-matrix Q"; input Real G[size(A, 1), size(A, 2)] "care-matrix G"; input Real X[size(A, 1), size(A, 2)] "solution of care"; output Real ku "upper bound of the care condition number"; end k_care_u;