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;