new benchmark |
mNo edit summary |
||
| Line 1: | Line 1: | ||
[[Category:benchmark]] | [[Category:benchmark]] | ||
[[Category:linear]] | [[Category:linear]] | ||
| Line 30: | Line 28: | ||
<references /> | <references /> | ||
==Contact== | ==Contact== | ||
[[User:Himpe|Christian Himpe]] | [[User:Himpe|Christian Himpe]] | ||
Revision as of 14:55, 22 May 2013
Description
The Inverse Lyapunov Procedure (ilp) is a synthetic random linear system generator. It is based on reversing the Balanced_Truncation procedure and was developed in [1], where a description of the algorithm is given.
A central point is the solution of the Lyapunov equations for the system matrix instead of the gramian matrix. This is feasable due to the symmetric (semi-)positive definiteness of the gramians and the requirement of a stable system. The solution will not be unique and include a symmetric system matrix, yet can be solved efficiently using empirical gramians.
Usage
To generate a random system using the Inverse Lyapunov Procedure download the M-file ilp.m. The function call requires three parameters; the number of inputs , of states and outputs . Optionally, a symmetric system can be enforced with the parameter . The return value consists of three matrices; the system matrix , the input matrix and the output matrix .
[A,B,C] = ilp(J,N,O,s);
The required Empirical Gramian Framework can be obtained from http://gramian.de. The ilp generator is compatible with wikipedia:MATLAB MATLAB and wikipedia:GNU_Octave OCTAVE.
References
- ↑ S.C. Smith, J. Fisher, "On generating random systems: a gramian approach", Proceedings of the American Control Conference, 2003.