Q = StateSpace.Analysis.observabilityMatrix(ss, method)
This function calculates the observability matrix
Q = [C; C*A; ...; C*A^(n-1)]
of the system
der(x) = A*x + B*u; y = C*x + D*u;
Modelica_LinearSystems2.StateSpace ss=Modelica_LinearSystems2.StateSpace( A=[-1, 1, 1;0, 1, 1;0, 0, 1], B=[0; 0; 1], C=[0, 1, 0], D=[0]); Real Q[3,3]; algorithm Q := Modelica_LinearSystems2.StateSpace.Analysis.observabilityMatrix(ss); // Q = [0, 1, 0; 0, 1, 1; 1, 1, 2]
encapsulated function observabilityMatrix import Modelica; import Modelica_LinearSystems2; import Modelica_LinearSystems2.StateSpace; input StateSpace ss "State space system"; output Real om[size(ss.A, 1)*size(ss.C, 1), size(ss.C, 2)] "Observability matrix"; end observabilityMatrix;