🔥 emgr 4 COMBINED REDUCTION - Combined State and Parameter Reduction via Empirical System Gramians

emgr - Empirical Gramian Framework

Model Reduction Software

Name

- EMpirical GRamian framework for (nonlinear) input-output systems.

Description

(image) #emgr box

System gramians are matrices associated to linear input-output systems and quantify their controllability, observability and minimality. Empirical gramian matrices (empirical gramians) are computable for linear but also nonlinear (state-space) systems and have, among others, application in model order reduction (MOR) or uncertainty quantification (UQ). Model reduction using empirical gramians can be applied to the state-space, to the parameter-space or both through combined reduction. For state reduction, empirical balanced truncation by the empirical controllability gramian and the empirical observability gramian, or alternatively, direct truncation (approximate balancing) by the empirical cross gramian (or the empirical linear cross gramian for large-scale linear systems) is available. For parameter reduction, parameter identification and sensitivity analysis by the empirical sensitivity gramian (controllability of parameters) or the empirical identifiability gramian (observability of parameters) are provided. Combined state and parameter reduction is enabled by the empirical joint gramian, which computes controllability and observability of states (cross gramian) and observability of parameters (empirical cross-identifiability gramian) concurrently. The empirical gramian framework - emgr - is a compact open-source toolbox for (empirical) GRAMIAN-based model reduction and compatible with OCTAVE and MATLAB. This mathematical software provides a common interface for the computation of empirical gramians and empirical covariance matrices.

Scope

Download

Get emgr here: emgr.m (Version: 5.4)    [mirror]    [source]    [meta]    [emgr_oct.m]    [emgr_lgc.m]

(emgr is written in the Matlab programming language and requires Octave or MATLAB, but emgr has no dependencies on toolboxes or packages.)

License

All source code is licensed under the open source BSD-2-clause license:

Copyright (c) 2013--2018, Christian Himpe

All rights reserved.

Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:

1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.

2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.

THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

Disclaimer

This is research software!

Usage

(flowchart) gramian interrelations

General Usage: W = emgr(f,g,s,t,w,pr,nf,ut,us,xs,um,xm,dp);

Minimal Usage: W = emgr(f,g,s,t,w);

About Info Usage: v = emgr('version');

Mandatory Arguments

Optional Arguments

Empirical Gramian Types

(chart) #emgr application matrix The empirical non-symmetric cross gramian WZ can be selected by setting option flag nf(7) = 1 and is compatible with the empirical cross gramian WX, the empirical linear cross gramian WY and the empirical joint gramian WJ.

Option Flags

Return Values

Configuration

SSP32 Runge-Kutta schematic Custom solver (via handle in global variable: ODE)

Utility

The platform framework curios.m (Clearing Up Reducibility of Input-Output Systems) enables testing and evaluating methods based on empirical gramians. curios has the following capabilities:

Mini Code Example

A = -0.5*eye(4) % System matrix B = [0;1;0;1] % Input matrix C = [0,0,1,1] % Output matrix WX = emgr(@(x,u,p,t) A*x + B*u, @(x,u,p,t) C*x,[1,4,1],[0.1,10.0],'x') % ≈ B*C

Tests

Tests are defined in emgrtest.m and are evaluated using curios. The tests are performed on a time-invariant, linear, state-space symmetric MIMO system, with a negative Lehmer matrix A as a system matrix and optional linear parametrization: ẋ = A*x + B*u + F*p y = C*x

Demos

The demo codes can be found in examples.m and are evaluated using curios.

Combined Reduction: Nonlinear System (WJ)

Here a parametrized nonlinear control system is reduced in states and parameters using the empirical joint gramian. The nonlinearity is given by a hyperbolic tangent with parametrized argument. This is a parametrized hyperbolic network model from "Modeling and control based on a new neural network model" by Y. Quan et al (2001). ẋ = A*tanh(p.*x) + B*u y = C*x

(plot) combined_wj reduction #emgr
Relative L2 output error for varying state and parameter dimension.

Benchmark: Inverse Sylvester Procedure (WC+WO)

The inverse lyapunov procedure can be used to generate random linear systems. Here, the required solution of the lyapunov equation is accomplished by the empirical gramians after factoring (-1). The generated random system is then renduced by using balanced truncation. The ilp.m control system generator is available from the MORwiki. 0 = A*WC + WC*A' - B'*B = (-1*a)*WC + WC*(-1*a') - B'*B = a*(-1*WC) + (-1*WC)*a' - B'*B

(plot) benchmark_ilp reduction #emgr
Relative L1, L2, L output error for varying state dimension.

Benchmark: Flexible Space Structures (WX)

A procedural benchmark modelling modal structural dynamics for applications such as truss masts for space environments (i.e.: COFS-1) based on "Model Reduction for Flexible Space Structures" by W. Gawronski and T. Williams (1991). This is a linear second order model, also listed at the MORwiki, which is reduced by approximate balancing via the empirical cross gramian.

(plot) benchmark_fss #emgr
Relative L1, L2, L output error for varying state dimension.

Benchmark: Nonlinear Model Reduction (WX)

An implementation of a model reduction benchmark from the MORwiki. Here the nonlinear RC-ladder benchmark, a resistor-capacitor cascade with nonlinear resistors (diodes), is tested. The empirical cross gramian is utilized to reduce this nonlinear system.

(plot) benchmark_non #emgr
Relative L1, L2, L output error for varying state dimension.

PDE Reduction: Advection Equation (WY)

Discretized with a first order upwind finite difference scheme, the 1D linear transport equation is a hyperbolic partial differential equation which is reduced by the empirical linear cross gramian and direct truncation. The left boundary of the domain is set as input, the right boundary is set as output. du/dt = -a du/dx x(0,t) = v(t) y(t) = x(1,t)

(plot) reduced advection #emgr
Relative L1, L2, L output error for varying state dimension.

Nonlinear Second Order Reduction: 5-body Choreography (WO)

A 5-body choreography, based on the gravitational N-body problem, resembling the figure eight is reduced in states; based on "New Families of Solutions in N-Body Problems" by C. Simo (2001).

(plot) reduced nbody #emgr
Reduced order 5-body orbits.

Parameter Identification: Stable Orbits Inside Black Holes (WS)

This is an implementation of: "Is there life inside black holes?" by V. Dokuchaev (2011), describing stable orbits of photons and particles inside the event horizon of a black hole. The underlying system is nonlinear and the output function is a coordinate transformation from Boyer-Lindquist to Cartesian coordinates. In the computation of particle and photon orbits the associated parameters are identified by the empirical sensitivity gramian.

(plot) blackhole orbit #emgr
Quasi-stable orbits inside the event horizon of a black hole.

About

A gramian matrix W is the result of all inner products of a set of vectors V = (v1 ... vn), in other words: W = VT V. Properties of (linear) control systems can be assessed by the system gramians, which are based on the controllability and observability operators. Classically, the controllability gramian and observability gramian are utilized in balancing model reduction methods. The cross gramian encodes controllability and observability information into a single matrix and thus minimality.

Empirical gramians extend this approach to nonlinear control systems and thus enable nonlinear model reduction. For linear systems, the empirical gramians are equal to the classic gramians. Yet, empirical gramians contain more information about the underlying system; and the empirical cross gramian conveys even additional information. Also, the Koopman controllability gramian and Koopman observability gramian are very similar to the empirical gramian counterparts. This makes empirical gramians a versatile tool for mathematical engineering.

The (discrete) empirical cross gramian encloses information on the input-output behavior of the associated control system as well as approximate Hankel Singular Values and can be computed very efficiently. For extreme-scale systems, the linear empirical cross gramian, related to Balanced POD, can be utilized. And for parametrized systems, the (empirical) joint gramian, derived from the cross gramian, is available for combined reduction. In case of custom input, an empirical covariance matrix can also be computed.

Model reduction or model order reduction is a research field in applied mathematics, which investigates the algorithmic construction of low-order surrogate models for high-dimensional differential equation models. The resulting reduced order models allow a substantially more economical evaluation in terms of computational or memory resources compared to the original large-scale full order model. Model reduction is essential as developments in computing technology are constantly outpaced by the requirements for model simulation in dimensionality or complexity.

References

Contact

Send feedback to: ch@gramian.de

Cite

gramian.de quick response code

Notes

Troubleshooting

See Also

Category

(plot) cross gramian surface #emgr
Λ