functionlinearize
Extends from Modelica.Icons.Function (Icon for functions).
Information
This function initializes a Modelica model and then simulates the model with its default experiment options until time instant "t_linearize". If t_linearize=0, no simulation takes place (only initialization). At the simulation stop time, the model is linearized in such a form that
- all top-level signals with prefix "input" are treated as inputs u(t) of the model,
- all top-level signals with prefix "output" are treated as outputs y(t) of the model,
- all variables that appear differentiated and that are selected as states at this time instant are treated as states x of the model.
Formally, the non-linear hybrid differential-algebraic equation system is therefore treated as the following ordinary equation system at time instant t_linearize:
der(x) = f(x,u)
y = g(x,u)
Taylor series expansion (linearization) of this model around the simulation stop time t_linearize:
u0 = u(t_linearize) y0 = y(t_linearize) x0 = x(t_linearize)
and neglecting higher order terms results in the following system:
der(x0 + dx) = f(x0,u0) + der(f,x)*dx + der(f,u)*du
y0 + dy = g(x0,u0) + der(g,x)*dx + der(g,u)*du
where der(f,x) is the partial derivative of f with respect to x, and the partial derivatives are computed at the linearization point t_linearize. Re-ordering of terms gives (note der(x0) = 0):
der(dx) = der(f,x)*dx + der(f,u)*du + f(x0,u0)
dy = der(g,x)*dx + der(g,u)*du + (g(x0,u0) - y0)
or
der(dx) = A*dx + B*du + f0
dy = C*dx + D*du
This function returns the matrices A, B, C, D and assumes that the linearization point is a steady-state point of the simulation (i.e., f(x0,u0) = 0). Additionally, the full Modelica names of all inputs, outputs and states shall be returned if possible (default is to return empty name strings).
Inputs
| Type | Name | Default | Description |
|---|---|---|---|
| String | modelName | Name of the Modelica model | |
| Modelica.Units.SI.Time | t_linearize | 0 | Simulate until T_linearize and then linearize |
Outputs
| Type | Name | Default | Description |
|---|---|---|---|
| Real[nx,nx] | A | ABCD[1:nx, 1:nx] | A-matrix |
| Real[nx,nu] | B | ABCD[1:nx, nx + 1:nx + nu] | B-matrix |
| Real[ny,nx] | C | ABCD[nx + 1:nx + ny, 1:nx] | C-matrix |
| Real[ny,nu] | D | ABCD[nx + 1:nx + ny, nx + 1:nx + nu] | D-matrix |
| String[nu] | inputNames | xuyName[nx + 1:nx + nu] | Modelica names of inputs |
| String[ny] | outputNames | xuyName[nx + nu + 1:nx + nu + ny] | Modelica names of outputs |
| String[nx] | stateNames | xuyName[1:nx] | Modelica names of states |