StabilityEquivalence
Authors: T. Cluzeau
I - Content
In the literature of linear 2D discrete systems, two particular classes of systems have been particularly studied:
- The Roesser models of the form
(R) x^h(i+1,j) = A11*x^h(i,j) + A12*x^v(i,j) + B1*u(i,j),
x^v(i,j+1) = A21*x^h(i,j) + A22*x^v(i,j) + B2*u(i,j),
where the Aij's and Bi's (i,j=1,2) are matrices with constant entries;
- The (generalized) Fornasini-Marchesini models of the form
(FM) x(i+1,j+1) = F1*x(i+1,j)+F2*x(i,j+1)+F3*x(i,j)+G1*u(i+1,j)+G2*u(i,j+1)+G3*u(i,j),
where the Fi's and Gi's (i=1,2,3) are matrices with constant entries.
Within the ANR MSDOS (Multidimensional Systems: Digression On
Stability) project, we have studied the following problems related to
these systems:
- The equivalence problem
in the sense of algebraic analysis: namely, can we transform any (R)
model into an equivalent (FM) model and vice versa?
- How can we use these equivalence transformations to study
structural stability and structural stabilization issues for these
classes of systems?
StabilityEquivalence is a related Maple package that can be used to:
- Check the structural stability of linear 2D discrete Roesser and generalized Fornasini-Marchesini models using the IsStable package,
- Transform (if possible) a Roesser model into an equivalent generalized Fornasini-Marchesini model and vice versa.
For more details, see the related publications:
[C2015] T. Cluzeau, Constructive algebraic analysis approach to the equivalence of multidimensional linear systems,
Proceedings of the IEEE 9th International Workshop on Multidimensional
(nD) Systems (IEEE nDS 2015), Vila Real (Portugal) (07-09/09/15).
[BCDY2016a] O. Bachelier, T. Cluzeau, R. David, N. Yeganefar, Structural stability and equivalence of linear 2D discrete systems,
Proceedings of the the IFAC Joint Symposium SSSC 2016 & TDS 2016, Istanbul (Turkey), 22-24/06/2016.
[BCDY2016b] O. Bachelier, T. Cluzeau, R. David, N. Yeganefar, Structural stabilization of linear 2D discrete systems using
equivalence transformations,
Accepted for publication in Multidimensional Systems and Signal Processing.
II - List of the procedures
Here are the procedures of our package. See the example file illustrating how to use them.
- FornasiniToRoesser2D
Input: the constant matrices Fi's and Gi's defining the (FM) and the dimension dx of the state variable x of the (FM) model
Output: the matrices A = (Aij)i,j and B=(B1^T B2^T)^T defining an equivalent (R) model and the respective dimensions dh(=dx) and
dv(=dx+du) of the horizontal and vertical state variables of the (R) model
- RoesserToFornasini2D
Input: the matrices A = (Aij)i,j and B defining the (R) model and the respective dimensions dh and dv of the horizontal and vertical
state variables of the (R) model
Output: the matrices Fi's and Gi's defining an equivalent (FM) model or "Case not handled by this procedure" (see [C2015])
- TestStructuralStabilityRoesser2D
Input: the matrix A of the (R) model and the dimensions dh and dv of the horizontal and vertical state vector
Output: true or false
- TestStructuralStabilityFornasini2D
Input: the matrices F1, F2 and F3 defining the (FM) model and the dimension dx of the state variable
Output: true or false
- FornasiniToOperatorMatrix2D
Input: the constant matrices Fi's and Gi's defining the (FM) and the dimension dx of the state variable of the (FM) model
Output: the matrix R of operators (sigma1 and sigma2) definig the associated linear system Ry=0
- RoesserToOperatorMatrix2D
Input: the matrices A = (Aij)i,j and B defining the (R) model and the respective dimensions dh and dv of the horizontal and vertical
state variables of the (R) model
Output: the matrix R of operators (sigma1 and sigma2) definig the associated linear system Ry=0
- TestStructuralStabilityAutonomousLinearSystem2D
Input: the matrix of operators R defining a square autonomous linear system
Output: true or false
III - Package and example file
You can get the package by downloading the following .mpl file: StabilityEquivalence.mpl
Here is a Maple file illustrating the use of the different procedures on the example considered in
[BCDY2016b]:
StabilityEquivalencePackageExample.mw (StabilityEquivalencePackageExample.pdf)
IV - Bug reports
You can now run your own examples.
If you have any problem with the package, find a bug or want to ask
questions, then contact me.