Next Article in Journal
The Cauchy Conjugate Gradient Algorithm with Random Fourier Features
Next Article in Special Issue
A Multi-Stage Homotopy Perturbation Method for the Fractional Lotka-Volterra Model
Previous Article in Journal
Effect of Voids Behind Lining on the Failure Behavior of Symmetrical Double-Arch Tunnels
Previous Article in Special Issue
New Hermite–Hadamard Type Inequalities Involving Non-Conformable Integral Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Stability Theory for Grünwald–Letnikov Inverse Model Control in the Multivariable LTI Fractional-Order Framework

by
Wojciech Przemysław Hunek
* and
Łukasz Wach
Department of Electrical, Control and Computer Engineering, Opole University of Technology, Prószkowska 76 Street, 45-758 Opole, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(10), 1322; https://doi.org/10.3390/sym11101322
Submission received: 4 October 2019 / Revised: 17 October 2019 / Accepted: 18 October 2019 / Published: 22 October 2019
(This article belongs to the Special Issue Recent Advances in Discrete and Fractional Mathematics)

Abstract

:
The new general theory dedicated to the stability for LTI MIMO, in particular nonsquare, fractional-order systems described by the Grünwald–Letnikov discrete-time state–space domain is presented in this paper. Such systems under inverse model control, principally MV/perfect control, represent a real research challenge due to an infinite number of solutions to the underlying inverse problem for nonsquare matrices. Therefore, the paper presents a new algorithm for fractional-order perfect control with corresponding stability formula involving recently given H- and σ -inverse of nonsquare matrices, up to now applied solely to the integer-order plants. On such foundation a new set of stability-related tools is introduced, among them the key role played by so-called control zeros. Control zeros constitute an extension of transmission zeros for nonsquare fractional-order LTI MIMO systems under inverse model control. Based on the sets of stable control zeros a minimum-phase behavior is specified because of the stability of newly defined perfect control law described in the non-integer-order framework. The whole theory is complemented by pole-free fractional-order perfect control paradigm, a special case of fractional-order perfect control strategy. A significant number of simulation examples confirm the correctness and research potential proposed in the paper methodology.

1. Introduction

Problems concerning stability of multivariable control systems described by both transfer function and state–space frameworks have recently attracted substantial research interest. A plethora of methods providing time and frequency domain approaches meets the requirements, nevertheless there is still considerable room and growing demand for the new techniques to overcome the constantly emerging control challenges. One of the remarkable and innovative methodology seems to be strategy based on the multivariable inverse model control, according to an infinite number of solutions underlying inverse problem for nonsquare matrices. The methodology presented in this paper provides a set of stability-devoted tools, derived from inverse model control, particularly MV/perfect control. It implements and develops the idea of recently introduced the so-called control zeros, until now defined in the discrete- and continuous-time integer-order frameworks. The control zeros play here fundamental role in determining stable/unstable properties for non-integer-order systems and allow predefinition of the stable and robust MV/perfect control strategies especially for nonsquare MIMO plants with different numbers of input and output variables. It should also be noticed that control zeros have first been derived for discrete-time integer-order LTI MIMO nonsquare systems described by the input–output model in form of multivariable transfer function [1]. They were defined as the points of the complex plain which determine the stability of the closed-loop minimum variance/perfect control procedure. It has been proven that for stable control zeros the perfect control algorithm is also stable, and the system should be called a minimum-phase plant [2,3,4,5]. On the other hand, the nonminimum-phase system with unstable control zeros is associated with unstable perfect control law. Next, the control zeros were defined for continuous-time plants with a definition of the minimum-phase systems stated in relation to the parallel stable MV/perfect control law [6]. Control zeros form an extension of transmission zeros for nonsquare LTI MIMO systems being under IMC, in particular MV/perfect control/perfect regulation [7,8,9,10,11]. Like the transmission zeros defined for SISO and square MIMO systems, the control zeros are in relation to the stabilizing potential of the MV/perfect control algorithm [12,13,14]. Placement of control zeros to desirable locations guaranteeing the robustness of such law is a difficult problem that has only been solved recently [15]. Due to existence of an infinite number of solutions associated with an application of nonunique inverses involved in the control algorithm, we can freely compose the MV/perfect control signals of the closed-loop plants. The case locating the control zeros at the origin of the complex plane with no transmission zeros seems to be very interesting, and the instance is obviously called a structurally stable pole-free scenario. The paper extends the entire theory of control zeros to cover the LTI MIMO discrete-time fractional-order plants described in the G-L framework [16,17,18]. The new perfect control strategy for such a class of systems is given together with related issue of minimum/nonminimum-phase property. Since the new tool is useful for synthesis of more general non-integer systems, it can provide a strong motivation to extend it for the case of continuous time.
The paper is organized as follows. The second section presents the fractional-order system representation in the state–space domain. Next section outlines the notion of perfect control for systems defined in Section 2. The definitions of right inverses and their application in the perfect control law are indicated in Section 4. Section 5 deals with the stability of discrete-time perfect control for both integer- and fractional-order cases. The following part of the paper provides the definitions of control zeros and minimum-phase property of fractional-order plants. In this unit the fractional-order pole-free perfect control strategy is also given. Vital issues involving robustness of examined systems are discussed in the subsequent section. Simulation examples included in the eighth section confirm the significant potential offered by the proposed new method. The final section of this paper contains summary and statement of the open problems.

2. Fractional-Order State–Space System Representation

Consider an LTI fractional-order plant S F O with n u -inputs, n y -outputs and n-state variable vector in discrete-time k described by
Δ α x ( k + 1 ) = A d x ( k ) + B u ( k ) , x ( 0 ) = x 0 y ( k ) = C x ( k ) ,
where the discrete-time difference operator Δ α is the G-L difference of an arbitrary fractional-order α , with 0 < α < 2 , as follows
Δ α x ( k ) = j = 0 k ( 1 ) j α j x ( k j ) ,
where
α j = 1 j = 0 α ( α 1 ) ( α j + 1 ) j ! j > 0 .
Assume also the system is controllable and observable. Please note that for α = 1 the plant as in Equations (1)–(3) comes down to the classical integer-order one S I O in the following form
x ( k + 1 ) = A x ( k ) + B u ( k ) , x ( 0 ) = x 0 y ( k ) = C x ( k ) . .
In both non-integer- and integer-order cases, we must assume that A R n × n , B R n × n u , C R n y × n with A d = A I n , whereas I n is the identity n-matrix.

3. Fractional-Order Perfect Control

Observe that the Equation (2) can be rewritten into the form
x ( k + 1 ) = Δ α x ( k + 1 ) j = 1 k + 1 ( 1 ) j α j x ( k j + 1 ) .
Consequently, after minimizing the (noise-free) control performance index for the plants with time delay d = 1
J ( u ) = k = 0 y ( k + 1 ) y ref ( k + 1 ) T y ( k + 1 ) y ref ( k + 1 ) ,
where
y ( k + 1 ) = C Δ α x ( k + 1 ) j = 1 k + 1 1 j α j x ( k j + 1 ) ,
and y ref ( k + 1 ) are the one-step deterministic output predictor and reference/setpoint, respectively, we obtain the perfect control law [19]
u ( k ) = ( CB ) R [ y ref ( k + 1 ) C A d x ( k ) + C j = 1 k + 1 ( 1 ) j α j x ( k j + 1 ) ,
where symbol ’ R ’ denotes any right inverse of product of CB [8].
Remark 1.
The Formula (8) is only valid for systems with regard to several input and output variables. We must take into account the so-called right-invertible plants fulfilling the condition n u n y . This problem is explained in detail further on in the paper.
Remark 2.
It is worth noticing that the non-integer-order prefect control law (8) can also be applied to systems blurred by zero-mean uncorrelated disturbance in form of white noise. However, minimum variance control problem goes beyond the scope of this paper.
Remark 3.
For α = 1 the Equation (8) reduces to the perfect control of integer-order systems as follows
u ( k ) = ( CB ) R [ y ref ( k + 1 ) C A x ( k ) ] .
With the notion of issues of both fractional-order system representation and relevant perfect control law, we can proceed to the problem of application of the appropriate inverse of the product CB .

4. An Application of Parameter Matrix Right Inverses into Fractional-Order Perfect Control Law

4.1. Selection of Nonunique Right Inverse

The issues regarding the stability and robustness of MV/perfect control for LTI MIMO integer-order discrete-time systems in state–space domain are subject to increasingly intensive scientific exploration [20,21]. The search for various methods is being conducted in particular to obtain new classes of right inverses. Until now, the mostly used and widely presented in the literature is minimum-norm right T-inverse [22,23]. T-inverse considered, however, as the ‘minimum-energy’ inverse, cannot be applied regarding a wide class of cases where the control inputs remain unstable under perfect control law. The recently introduced σ -inverse and H-inverse with various types of so-called degrees of freedom offer powerful tools to be applied regarding the above class of problems [24,25]. It is interesting that the first inverse constitutes an extension of minimum-norm right T-inverse in which degrees of freedom are organized as separate matrices being both in parameter as well as polynomial form. On the other hand, the second of the inverses is based on well-known SVD factorization. That approach gives the opportunity to impact on the perfect control strategy by influencing matrices obtained from SVD factorization in form of parameter or polynomial degrees of freedom as well. The tools mentioned above, previously employed to integer-order systems, are transferred to fractional-order perfect control of LTI MIMO discrete-time state–space systems [26]. It is remarkable that the concept of engaging the inverses of parameter matrices to increasing the robustness of multivariable discrete-time fractional-order perfect control structures in a state–space forms a new and original idea developed by these authors, not presented so far. Therefore, it is an attempt to synthesis of discrete-time fractional-order perfect control with regard of its stability and robustness. To obtain that objective, the recently devised a fractional-order multivariable discrete-time perfect control algorithm, in particular dedicated to the so-called nonsquare state–space systems, i.e., systems with different number of input and output variables, is used [19]. The simulations performed into perfect control involve parameter σ -inverse and H-inverse. The results of research conducted in MATLAB/Simulink environment show that use of appropriate inverse machinery may determine stability and, consequently, minimum-phase property of the non-integer plants. The definitions of the right inverses presented in the next section constitute the essential element for calculation of the control zeros, as they affect the properties of the fractional-order perfect control systems discussed in this paper.

4.2. Inverses of Nonsquare (Parameter) Matrices

Let us start our considerations with the definitions of the classical minimum-norm inverses of polynomial matrix W ̲ ( q 1 ) in the backward shift operator q 1 , afterwards called T-inverse [8,27,28,29,30,31]. The following two definitions are dedicated to full and non-full normal rank polynomial matrices W ̲ ( q 1 ) , respectively.
Definition 1.
Let the polynomial matrix W ̲ ( q 1 ) = w ̲ 0 + w ̲ 1 q 1 + + w ̲ m q m be of full normal rank n y (or n u ). The (unique) minimum-norm right (or the least-squares left) T-inverse of W ̲ ( q 1 ) is defined as W ̲ 0 R ( q 1 ) = W ̲ T ( q 1 ) [ W ̲ ( q 1 ) W ̲ T ( q 1 ) ] 1 (or W ̲ 0 L ( q 1 ) = [ W ̲ T ( q 1 ) W ̲ ( q 1 ) ] 1 W ̲ T ( q 1 ) ).
Definition 2.
Let the polynomial matrix W ̲ ( q 1 ) = w ̲ 0 + w ̲ 1 q 1 + + w ̲ m q m of non-full normal rank r be skeleton-factorized as W ̲ ( q 1 ) = C ̲ ( q 1 ) D ̲ ( q 1 ) , where d i m [ W ̲ ( q 1 ) ] = n y × n u , d i m [ C ̲ ( q 1 ) ] = n y × r , d i m [ D ̲ ( q 1 ) ] = r × n u . The (unique) Moore-Penrose T-inverse of W ̲ ( q 1 ) is defined as W ̲ 0 # ( q 1 ) = D ̲ 0 R ( q 1 ) C ̲ 0 L ( q 1 ) , where D ̲ 0 R ( q 1 ) = D ̲ T ( q 1 ) [ D ̲ ( q 1 ) D ̲ T ( q 1 ) ] 1 and C ̲ 0 L ( q 1 ) = [ C ̲ T ( q 1 ) C ̲ ( q 1 ) ] 1 C ̲ T ( q 1 ) .
Remark 4.
The systems taken into account are parameter plants which constitute a special case of polynomial ones.
An extension of the unique right T-inverse is nonunique right σ -inverse including degrees of freedom. The newest definition of σ -inverse can be found in Ref. [32] in the form of the following corollary
Corollary 1.
Let the polynomial matrix W ̲ ( q 1 ) = w ̲ 0 + w ̲ 1 q 1 + + w ̲ m q m be of full normal rank n y (or n u ) and let z s β ̲ ( z 1 ) = β ̲ ( z ) R n y × n u [ z ] be arbitrary, including an arbitrary order s. Assume additionally that the product W ̲ ( q 1 ) β ̲ T ( q 1 ) is of full normal rank n y (or n u ). Then a right σ-inverse can be defined as
W ̲ R ( q 1 ) = β ̲ T ( q 1 ) W ̲ ( q 1 ) β ̲ T ( q 1 ) 1 .
Remark 5.
Naturally, the presented polynomial matrix right σ-inverse reduces to the parameter one for W ̲ ( q 1 ) = w ̲ 0 and takes the form
W ̲ R ( q 1 ) = β T W ̲ ( q 1 ) β T 1 ,
with parameter degrees of freedom β. Some issues concerning β with parameter matrices can be found in [33].
Remark 6.
It is clear now that in Equation (11) the W R should be used instead of W ̲ R ( q 1 ) , as well as W instead of W ̲ ( q 1 ) .
Remark 7.
For β ̲ ( q 1 ) = W ̲ ( q 1 ) the σ-inverse specializes to the T-inverse and in the parameter case is as follows
( W ) 0 R = W T [ W W T ] 1 ,
where symbol ( . ) 0 R denotes the unique minimum-norm right T-inverse.
Finally, the recently introduced nonunique parameter right H-inverse based on SVD factorization is shown below.
Theorem 1
([24]). Consider an m × n matrix A being of full rank m under the SVD decomposition A = U Σ V T , where U and V are unitary matrices and Σ = d i a g ( σ 1 , , σ m 0 m × ( n m ) ) includes eigenvalues of A . Then the right matrix H-inverse of A can be given as
A R = ( V T ) 1 Σ R U 1 ,
where Σ R = Σ m 1 L and Σ m 1 = d i a g ( 1 σ 1 , , 1 σ m ) while L R ( n m ) × m ( q 1 ) is an arbitrary matrix polynomial in the backward shift operator domain.
Remark 8.
As in preceding case, the presented polynomial matrix L ( q 1 ) , including degrees of freedom, can be reduced to the parameter form L .
Remark 9.
For L = 0 we obtain the classical minimum-norm right or least-square left T-inverse [34].
With the needed notion of polynomial inverses and their parameter representations, we are able to proceed with the implementation of right T-, σ - and H-inverse into the fractional-order perfect control law as in Equation (8), which is the subject of the subsequent section.

4.3. Employment of Right Inverses into Fractional-Order Perfect Control Law

Starting to consider the implementation outlined in former section right inverses into Equation (8) we must assume that the product of CB is of full rank. First, we take into account the classical minimum-norm T-inverse. It is important to note that the minimum-norm T-inverse represents a special case of a more general one related to σ -inverse (see Remark 7). Therefore, regarding Equation (8), the fractional-order perfect control formula for LTI MIMO discrete-time systems in state–space can be presented in the following form
u ( k ) = ( CB ) T [ CB ( CB ) T ] 1 y ref ( k + 1 ) CA d x ( k ) + C j = 1 k + 1 1 j α j x ( k j + 1 ) .
Secondly, to implement σ -inverse, we rewrite Equation (8) according to the Corollary 1 with corresponding remarks as follows
u ( k ) = β T [ CB β T ] 1 y ref ( k + 1 ) CA d x ( k ) + C j = 1 k + 1 1 j α j x ( k j + 1 ) ,
with parameter degrees of freedom matrix β of n y × n u dimension.
Eventually, the SVD-based H-inverse formula defined in the Theorem 1 is applied to Equation (8). Hence, we obtain fractional-order perfect control law related to H-inverse in the following form
u ( k ) = ( V CB T ) 1 Σ CB R U CB 1 y ref ( k + 1 ) CA d x ( k ) + C j = 1 k + 1 1 j α j x ( k j + 1 ) .
With the notion of the three different forms of fractional-order perfect control law we need to compile necessary conditions, as presented in the next section.

4.4. Necessary Conditions for Application of Right Inverses to the Product of CB

Taking into account the nonsquare MIMO systems with different number of input and output variables, the problem of the inverse of nonsquare matrices arises since the perfect control bases on the IMC structure. The selection of (non)unique right inverses for nonsquare MIMO systems is essential for the problem of the stability of perfect control discussed in the next sections. Because of the fact that we consider plants with n u n y , the product of CB must be of full rank. For non-full rank CB product the prefect control does not exist [8]. Therefore, first it must be assumed that the product of CB is only of full rank n y . In fact, the right-invertible systems in which input variables outnumber the output ones ( n u > n y ) should be taken into account. For left-invertible systems, i.e., plants with n u < n y , the perfect control cannot be established. These two requirements constitute the basic conditions for the adequate application of (non)unique inverses.
Moreover, the issues regarding the regular inverse of the square systems, i.e. systems with n u = n y , are not touched in this paper as the unique T-inverse and nonunique σ and H inverses reduce to the regular one. For square systems, including the SISO ones, the regular unique matrix inverse of CB must be applied. Such an inverse does not have degrees of freedom, which are crucial for the subject of the paper.
We can also note that the minimum-norm T-inverse applied in (non-)integer-order perfect control problems in many cases may cause an instability of the whole control system. Therefore, in perfect control robustness research there is a need to employ other nonunique inverses, for example said σ -inverse or H-inverse implemented in Formulas (15) and (16), respectively. Actually, the application of the other right inverses could be related to the problem of control zeros placement and finally determined by the appropriate selection of the degrees of freedom originated from the above-mentioned inverses. The problem of the instability of multivariable fractional-order perfect control systems defined in discrete-time state–space framework is discussed in the following sections and illustrated by simulation examples.

5. Stability of Fractional-Order Discrete-Time Perfect Control

The problem of stability of multivariable integer-order perfect control systems, described by both transfer function and state–space frameworks, is related to the so-called control zeros [1,6]. Therefore, it is important to consider the control zeros as the vital element connected with the perfect control stability. Thus, the location of different sets of control zeros, being an extension of transmission zeros for the nonsquare systems, in the complex plane can be selected based on the application of various nonunique right inverses. It is crucial that in this study we use a similar paradigm as in the case of integer-order plants to determine the stability of perfect control for LTI MIMO fractional-order systems in state–space framework. As the right inverses are essential during the design of robust perfect control strategy, they can be treated as an effective tool for robustification of MV/perfect control strategies.

5.1. Integer-Order Instance

The stability of classical integer-order discrete-time deterministic perfect control law for state–space right-invertible systems described as follows
x ( k + 1 ) = A x ( k ) + Bu ( k ) y ( k ) = C x ( k ) x ( 0 ) = x 0 ,
where
u ( k ) = K I O x ( k ) ,
while
K I O = ( CB ) R CA ,
(see Equation (9) for y ref = 0 ) is determined by the well-known formula [27]
det ( z I n A + BK I O ) .
Remark 10.
For the simplicity of the conducted research study we take into account the consideration that y ref = 0 , as it does not distort the perfect control stability theorem in the general case. This statement is confirmed by the following theorem.
Theorem 2.
Let us consider systems described by Equation (17). The Formula (20) can be applied to the perfect control laws as in Equation (9) and (18).
Proof. 
Immediately, after substitution Equation (9) into the state equation of Formula (17) we obtain
x ( k + 1 ) = ( A BK I O ) x ( k ) + B ( CB ) R y ref ( k + 1 ) ,
where the second part of the right side of equation does not influence the stability. □
Remark 11.
Observe that the same results can be achieved for plants with a time delay d > 1 .
Following the introduction of issues related to integer-order systems, the more general fractional-order case is presented in the next section.

5.2. Fractional-Order Instance

Consider the discrete-time fractional-order plants defined by Equations (1)–(3). As mentioned above, the stability of such systems, as the superclass of integer-order case, should be analyzed in a corresponding manner to integer-order ones. Thus, by analogy, we assume that the whole machinery of Section 5.1 can be successfully transferred to the non-integer-order instance. Therefore, Equation (8) should be taken into account to identify the matrix K F O . Consequently, after a few manipulations, the Equation (8) with y ref = 0 can be transformed into the following formula
u ( k ) = CB R C A d I n j = 1 k + 1 ( 1 ) j α j q j + 1 x ( k ) ,
where q is the backward shift operator. Finally, the matrix K F O of fractional-order closed-loop perfect control systems is represented as follows
K F O = ( C B ) R C A d I n j = 1 k + 1 ( 1 ) j α j q j + 1 .
With the required tools, let us try to introduce the theorem related to the stability of non-integer-order plants.
Theorem 3.
Let the LTI MIMO discrete-time fractional-order systems S F O as in Equations (1)–(3) is described by G-L framework. Then a stability region of the perfect control defined by Equation (8) is restricted by the following formula
det z I n I n B ( CB ) R C A d I n j = 1 k + 1 1 j α j z j + 1 ,
where z is the Laplace operator.
Proof. 
After substitution of the prefect control law (22) to the plants S F O defined by Equation (1) and altered by Formula (5) we obtain
x ( k + 1 ) + j = 1 k + 1 ( 1 ) j α j x ( k j + 1 ) = A d x ( k ) + B ( CB ) R C A d + I n j = 1 k + 1 ( 1 ) j α j q j + 1 x ( k ) ,
and, finally, a simplifed form is derived, as follows
x ( k + 1 ) = I n B ( CB ) R C A d I n j = 1 k + 1 ( 1 ) j α j q j + 1 x ( k ) .
Thus, transferring Equation (26) to the z-domain we receive the subsequent statement
X ( z ) z = I n B ( CB ) R C A d I n j = 1 k + 1 ( 1 ) j α j z j + 1 X ( z ) ,
which holds the Equation (24). □
Remark 12.
Observe that Equation (24) corresponds to the succeeding expression
det ( z I n A F O + BK F O ) ,
where K F O is defined by the statement (23) and
A F O = A d I n j = 1 k + 1 ( 1 ) j α j x ( j k + 1 ) .
Of course, Equation (28) for α = 1 reduces to the Formula (20) dedicated to the integer-order perfect control law.
Remark 13.
Naturally, the number of roots of Equation (24) or Formula (28) must be understood in terms of simulation horizon k selected by the designer.
Concluding, let us formulate the practical theorem which allows us in a simple way to transfer formulas between integer- and non-integer-order domains.
Theorem 4.
The equations defined by integer- and non-integer-order frameworks can be considered in a similar manner under formula
A = A F O .
Proof. 
Immediately, after using Equation (30). □
Remark 14.
The results, which were obtained, are still valid for stochastic case of perfect control, i.e., the minimum variance control (MVC).
Once we know the major issues involving the stability area of fractional-order perfect control system we can proceed to define the control zeros being the main goal of this paper.

6. Control Zeros and Minimum-Phase Property

The idea of control zeros was first introduced for discrete-time LTI MIMO integer-order plants described by the input–output model [1] and extended to cover the continuous-time systems [6]. Control zeros were defined as the points of the respective complex plane which indicated the stable/unstable property of the MV/perfect control strategy. For square systems, the set of control zeros was calculated by the application of the regular inverse to some matrix involving the dynamic description of the analyzed plant. In the cases of nonsquare plants with different number of input–output variables the above-mentioned control zeros were determined by employing unique/nonunique inverses, e.g., given in Section 4.2. Therefore, a proper set of inverses should be selected to make the IMC structure stable/unstable.
The proposition to apply the same paradigm to perfect control fractional-order plants is presented in the following sections.

6.1. Control Zeros

Definition 3.
Consider a discrete-time LTI MIMO fractional-order system S F O described by G–L state–space framework with product of CB being right-invertible. The complex number ζ C is called a control zero of the plant S F O iff it is a root of Formulas (24) or (28).
Observe that similarly to the case of the LTI MIMO integer-order plants defined by the input–output models, the control zeros determine the stability of MV/perfect control formula. If the control zeros are located outside the unit circle the above-mentioned algorithm is unstable and, contrary, if they are placed inside the unit disk, the said control remains asymptotically stable. The statement fulfills the classical stability paradigm of discrete-time systems. Moreover, according to the Remark 13, several control zeros must be respected regarding time horizon k. This feature is presented by the following theorem and confirmed by simulation examples of Section 8.
Theorem 5.
Consider an LTI fractional-order plant S F O with n u inputs, n y outputs and n-state vector in discrete-time k described by Equations (1)–(3). The total number of control zeros for such system can be calculated according to the formula
N k = n ( k + 1 ) ,
where n denotes the amount of state variables.
Proof. 
After using the Formula (24) or (28) we obtain the following expression det ( F ( z ) ) where
F n × n ( z ) = P 11 a 12 a 1 n a 21 P 22 a 2 n a n 1 a n 2 P n n ,
P i i = a i i ( 1 ) z + a i i ( 2 ) + a i i ( 3 ) z 1 + + a i i ( k + 2 ) z k , i = 1 , , n ,
is a matrix with polynomials on the main diagonal and a i j , i = 1 , , n , j = 1 , , m , are some parameters while a i j ( 1 ) = 1 for i = j . Then, the total number of control zeros depends on dimension n and the time horizon k. □
Remark 15.
We can note that the control zeros constitute a generalization of so-called transmission zeros on the nonsquare systems defined by the input–output model [8]. The problem of the transmission zeros for the state–space approach and their impact on the fractional-order perfect control goes beyond the scope of this paper and should be understood as a mathematical peculiarity.
With the tools for placement of control zeros, let us define the special pole-free fractional-order perfect control case.

6.2. Pole-Free Fractional-Order Perfect Control

Observe that under some assumptions resulted in det ( F ( z ) ) = z n we obtain a special case of fractional-order perfect control algorithm called pole-free perfect control strategy. In such a case, the poles of the closed-loop plants or rather control zeros are located at zero, i.e., are placed at the origin of the unit circle. In this scenario we achieve a maximum-speed/maximum accuracy perfect control law in opposite to the stable-pole instance, in which the control zeros lie inside the unit disk.
Following the Formulas (24) or (28) we can state that pole-free perfect control is established when I n B ( CB ) R C = 0 . Therefore, for B ( CB ) R C = I n we can distinguish subsequent instances: (1) B and C are regular-invertible, (2) B is right-invertible and C is regular-invertible, (3) B is regular-invertible and C is right-invertible, (4) B and C is right-invertible and (5) B is right-invertible and C is left-invertible. The above-mentioned expression holds under the first and second instances. In such cases we consider I n B B 1 C 1 C = 0 and I n B B R C 1 C = 0 . On the other hand, the third and fourth instances cannot be established due to the fact that ( CB ) R B 1 C R and ( CB ) R B R C R , respectively. The fifth instance does not meet our requirements providing product of CB with full normal rank. We should rather use the well-known Moore-Penrose inverse in such a scenario. Of course, we should also respect the case in which [ A d I ] = 0 . Such case occurs if α I n = A d | j = 1 . Summing up, we introduce below the definition of pole-free non-integer-order perfect control.
Definition 4.
Consider LTI MIMO discrete-time fractional-order plants S F O described by the G–L state–space framework as in Equations (1)–(3). The pole-free fractional-order perfect control law has to be understood in context of Formula (8) being under one of the following assumptions: 1) I n B B 1 C 1 C = 0 , 2) I n B B R C 1 C = 0 , where symbol ’ R ’ denotes any right inverse of CB product.

6.3. Minimum-Phase Property of the Discrete-Time Fractional-Order Systems

Taking into consideration the guidelines outlined in the previous sections there is a possibility to introduce in this unit the definition of minimum-phase behavior of the fractional-order LTI MIMO plants described in the discrete-time G-L domain. For the square/nonsquare systems represented by the input–output approach the minimum/nonminimum-phase property has been defined based on stable/unstable IMC or rather perfect control plant. It is connected with the stable/unstable control zeros sets obtained using different types of nonsquare inverses. The idea of the new definition of the minimum-phase S F O systems is proposed hereafter.
Definition 5.
An LTI MIMO fractional-order system defined by the G–L discrete-time state–space framework is called minimum-phase iff it is stably invertible in context of any generalized inverse of the product CB , otherwise the system is called nonminimum-phase.
Remark 16.
According to Definition 5 it can be concluded that the minimum-phase system in parallel is defined as that being stable under MV/perfect control algorithm. This scenario is related to the set of stable control zeros obtained using some generalized inverse; in case of unstable MV/ perfect control procedure connected with the unstable control zeros the nonminimum-phase system is regarded.
Having established the entire control zeros theory dedicated to the fractional-order perfect control, in the next section the robustness issue is touched as the crucial problem that is faced in such a strategy.

7. Robustness of Fractional-Order Perfect Control

It is well-known that the main disadvantage of the perfect control systems described by the integer-order discrete-time state–space framework is associated with the lack of robustness [8]. The high value of control variables beyond technological boundaries can occur in several practical cases. This feature makes the MV/perfect control algorithm limited in engineering applications [7,10,11]. A set of tools has been introduced as the remedy for this obstacle, among them, an important role is taken on by inverses of nonsquare matrices. As mentioned in Section 4 the implementation of e.g., different nonsquare inverses in integer-order perfect control law influences its stability and robustness [20]. The same paradigm can be adopted to the perfect control for LTI MIMO nonsquare discrete-time fractional-order plants defined by the G-L statement. This fact is confirmed in a few reports in the area; however, the formal proof is still to be conducted. It turns out that the minimum-energy-based approach can improve the fractional-order closed-loop perfect control robustness [26]. The chosen degrees of freedom { β ̲ ( q 1 ) } for, e.g., σ -inverse substituted to the fractional-order perfect control law as in (15) can provide its stable/robust instance. Therefore, the selection process of proper β ̲ ( q 1 ) should be realized under procedure subject to the minimum-energy constraint for the fractional-order perfect control signal input, involving performance index as the one below
E 1 = min { β ̲ ( q 1 ) } k = 0 N 1 u ( k ) u ss ( k ) T u ( k ) u ss ( k ) ,
where symbol ( . ) ss stands for steady state of the runs and N 1 denotes a time horizon over all possible β ̲ ( q 1 ) . An alternative way seems to be use of performance index in form of
E 2 = min { β ̲ ( q 1 ) } k = 0 N 1 u T ( k ) u ( k ) .
Generally, the optimization procedure using index (33) or (34) leads the control zeros to the stable region. This observation is given based on simulation studies.
Remark 17.
All the possible β ̲ ( q 1 ) should be considered in terms of their dimensions derived from number of input and output fractional-order perfect control variables. In our research we only analyze parameter β ̲ ( q 1 ) . Therefore, the zero-order of its components does not generate additional possible degrees of freedom.
Remark 18.
Similar techniques providing the robustification of different non-integer-order control laws can be find in Refs. [35,36,37,38,39,40,41].
Familiarity with the whole theory concerning the control zeros for LTI MIMO fractional-order G–L plants defined in the discrete-time state–space domain allows us to proceed to the next simulation step related to confirmation the new ideas introduced in the paper.

8. Simulation Examples

We start the scenarios with stability of fractional-order perfect control being the crucial goal of our considerations. The sets of control zeros are presented for the assumed fractional-order plants with various operators α and different degrees of freedom of parameter σ -inverse.

8.1. Control Zeros and Fractional-Order Perfect Control Stability

Consider the fractional-order system described by the G-L framework S F O ( A d , B , C , α ) with three inputs and two outputs represented by the following initial matrices
A d = 1.0074 0.0904 0.2564 0.0390 0.0112 1.0188 0.0722 0.2614 0.2558 0.1174 0.9684 0.0172 0.0995 0.2274 0.0924 1.0378 , B = 1.2711 0.4173 0 0.2963 0.1041 0.9926 0 0.5119 0 1.4916 0.5796 0.2039 , C = [ 0 0.5355 0.5290 0.7724 0.1737 0 0 1.2248 ] , X 0 T = [ 1.4851 0.1159 1.5747 2.4256 ] .
We assume that during the simulations, the matrices B and C are the subjects to change, to finally derive several various plants. For every pair ( B , C ) the fractional-order perfect control schemes are established through an application of σ -inverse with arbitrary degrees of freedom β . Figure 1 presents the stability areas restricted by the sets of stable control zeros calculated according to the different operators α . Every set of stable control zeros is obtained in the MATLAB/Simulink environment based on the stable fractional-order perfect control law. Naturally, the fractional-order perfect control systems prepared in Simulink have been supported by randomly chosen matrices B , C and β resulted from MATLAB script.
Figure 2 and Figure 3 show the stable fractional-order perfect control runs for initial plants S ( A d , B , C , α = 0 . 75 ) and S ( A d , B , C , α = 1 . 25 ) together with the respective stability areas.

8.2. Energy-Based Approach to Robustness of Fractional-Order Perfect Control

It is clear now that the proper selection of control zeros set affects the stability of the fractional-order perfect control algorithm. The second feature, strictly related to the control zeros, is robustness of such control strategy. In the next example this property is examined by a relatively simple simulation scenario. To show the fractional-order perfect control robustification the H-inverse, defined in Equation (13), and σ -inverse, given in Equation (10), are effectively employed as in Equations (16) and (15), respectively. This operation is essential in the cases where T-inverse (see Definition 1 and Equation (14)) is unstable.
Let us take into account the fractional-order system (35). After application of a right T-inverse we obtain unstable perfect control runs for chosen α 1 = 0 . 01 and α 2 = 1 . 8 . However, for specific σ - and H-inverses, selected under optimization tasks in the MATLAB/Simulink environment, the control signal achieves the minimum-energy value according to the performance indices (33) or (34). Again, the fractional-order perfect control plants performed in Simulink have been assisted by degrees of freedom β and L according to the MATLAB’s genetic algorithm mechanism. The results are piloted in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.

9. Conclusions and Open Problems

The entire theory concerning the perfect control dedicated to multivariable fractional-order plants in the Grünwald–Letnikov state–space framework is presented in the paper. Due to the application of the newly introduced control zeros, associated with the discrete-time perfect control formula, the minimum/nonminimum-phase behavior of analyzed non-integer-order systems is defined. The results reported in this paper were confirmed by several simulation examples performed in the MATLAB/Simulink environment. The further research study can be focused on the implementation of the postulated ideas to cover the discrete-time minimum variance control algorithm. An extension of the introduced theory to the continuous-time fractional-order case is another open problem to be taken into account in the future research.

Author Contributions

The authors contributed equally to this work.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
G-LGrünwald–Letnikov
IMCInverse Model Control
LTILinear Time-Invariant
MIMOMultiple-Input/Multiple-Output
MVCMinimum Variance Control
S F O Fractional-Order Plant S ( A d , B , C )
S I O Integer-Order Plant S ( A , B , C )

References

  1. Latawiec, K.; Bańka, S.; Tokarzewski, J. Control zeros and nonminimum phase LTI MIMO systems. Annu. Rev. Control. 2000, 24, 105–112. [Google Scholar] [CrossRef]
  2. Borisson, U. Self-tuning regulators for a class of multivariable systems. Automatica 1979, 15, 209–215. [Google Scholar] [CrossRef]
  3. Åström, K.J. Introduction to Stochastic Control Theory; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  4. Åström, K.J.; Wittenmark, B. On self-tuning regulators. Automatica 1973, 9, 185–199. [Google Scholar] [CrossRef]
  5. Ryan, E.P.; Sangwin, C.J.; Townsend, P. Controlled functional differential equations: Approximate and exact asymptotic tracking with prescribed transient performance. ESAIM COCV 2009, 15, 745–762. [Google Scholar] [CrossRef]
  6. Hunek, W.P. Control Zeros for Continuous-Time LTI MIMO Systems and Their Application in the Theory of Circuits and Systems. Ph.D. Thesis, Opole University of Technology Press, Opole, Poland, 2003. [Google Scholar]
  7. Dadhich, S.; Birk, W. Analysis and control of an extended Quadruple tank process. In Proceedings of the 13th IEEE European Control Conference (ECC’2014), Strasbourg, France, 24–27 June 2014; pp. 838–843. [Google Scholar] [CrossRef]
  8. Hunek, W.P.; Latawiec, K. A study on new right/left inverses of nonsquare polynomial matrices. Int. J. Appl. Math. Comput. Sci. 2011, 21, 331–348. [Google Scholar] [CrossRef] [Green Version]
  9. Lin, Z.; Saberi, A.; Sannuti, P.; Shamash, Y.A. Perfect regulation of linear discrete-time systems: A low-gain-based design approach. Automatica 1996, 32, 1085–1091. [Google Scholar] [CrossRef]
  10. Zhang, T.; Li, H.G.; Zhong, Z.Y.; Cai, G.P. Hysteresis model and adaptive vibration suppression for a smart beam with time delay. J. Sound Vib. 2015, 358, 35–47. [Google Scholar] [CrossRef]
  11. Zhang, T.; Yang, B.T.; Li, H.G.; Meng, G. Dynamic modeling and adaptive vibration control study for giant magnetostrictive actuators. Sens. Actuators A Phys. 2013, 190, 96–105. [Google Scholar] [CrossRef]
  12. Tokarzewski, J. Zeros in Linear Systems: A Geometric Approach; Warsaw University of Technology Press: Warsaw, Poland, 2002. [Google Scholar]
  13. Tokarzewski, J. A note on some characterization of invariant zeros in singular systems and algebraic criteria of nondegeneracy. Int. J. Appl. Math. Comput. Sci. 2004, 14, 149–159. [Google Scholar]
  14. Tokarzewski, J. Finite Zeros in Discrete Time Control Systems; Lecture Notes in Control and Information Sciences; Springer: Berlin, Germany, 2006; Volume 338. [Google Scholar] [CrossRef]
  15. Hunek, W.P.; Krok, M. Pole-free perfect control for nonsquare LTI discrete-time systems with two state variables. In Proceedings of the 2017 13th International Conference on Control & Automation (ICCA), Xi’an, China, 20–23 August 2017; pp. 329–334. [Google Scholar]
  16. Nastac, S.; Debeleac, C.; Vlase, S. Hysteretically Symmetrical Evolution of Elastomers-Based Vibration Isolators within α-Fractional Nonlinear Computational Dynamics. Symmetry 2019, 11, 924. [Google Scholar] [CrossRef]
  17. Jajarmi, A.; Hajipour, M.; Mohammadzadeh, E.; Baleanu, D. A new approach for the nonlinear fractional optimal control problems with external persistent disturbances. J. Frankl. Inst. 2018, 355, 3938–3967. [Google Scholar] [CrossRef]
  18. Wei, Y.Q.; Liu, D.Y.; Boutat, D. Innovative fractional derivative estimation of the pseudo-state for a class of fractional order linear systems. Automatica 2019, 99, 157–166. [Google Scholar] [CrossRef]
  19. Wach, Ł.; Hunek, W.P. Perfect control for fractional-order multivariable discrete-time systems. In Theoretical Developments and Applications of Non-Integer Order Systems; Springer: Cham, Switzerland, 2016; pp. 233–237. [Google Scholar]
  20. Hunek, W.P. An application of new polynomial matrix σ-inverse in minimum-energy design of robust minimum variance control for nonsquare LTI MIMO systems. IFAC-PapersOnLine 2015, 48, 150–154. [Google Scholar] [CrossRef]
  21. Hunek, W.P.; Krok, M. A study on a new criterion for minimum-energy perfect control in the state-space framework. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2019. [Google Scholar] [CrossRef]
  22. Wang, G.; Wei, Y.; Qia, S. Generalized Inverses: Theory and Computations; Science Press: Beijing, China, 2018. [Google Scholar]
  23. Wei, Y.; Stanimirović, P.S.; Petković, M.D. Numerical and Symbolic Computations of Generalized Inverses; World Scientific Publishing Co. Pte Ltd.: Singapore, 2018. [Google Scholar]
  24. Hunek, W.P. New SVD-based matrix H-inverse vs. minimum-energy perfect control design for state-space LTI MIMO systems. In Proceedings of the 2016 20th International Conference on System Theory, Control and Computing (ICSTCC), Sinaia, Romania, 13–15 October 2016; pp. 14–19. [Google Scholar]
  25. Noueili, L.; Chagra, W.; Ksouri, M. New Iterative Learning Control Algorithm Using Learning Gain Based on σ Inversion for Nonsquare Multi-Input Multi-Output Systems. Model. Simul. Eng. 2018, 2018, 4195938. [Google Scholar] [CrossRef]
  26. Hunek, W.P.; Wach, Ł. Towards a New Stability Criterion for Fractional-Order Perfect Control of LTI MIMO Discrete-Time Systems in State Space. In Proceedings of the 2017 3rd IEEE International Conference on Cybernetics (CYBCONF), Exeter, UK, 21–23 June 2017; pp. 1–6. [Google Scholar]
  27. De la Sen, M. On Cauchy’s Interlacing Theorem and the Stability of a Class of Linear Discrete Aggregation Models Under Eventual Linear Output Feedback Controls. Symmetry 2019, 11, 712. [Google Scholar] [CrossRef]
  28. Ben-Israel, A.; Greville, T. Generalized Inverses, Theory and Applications, 2nd ed.; Springer: New York, NY, USA, 2003. [Google Scholar]
  29. Bronnikov, A.; Borovkov, A. Inverse control of discrete-time multivariable systems. J. Frankl. Inst. 2002, 339, 335–345. [Google Scholar] [CrossRef]
  30. Karampetakis, N.P.; Tzekis, P. On the computation of the generalized inverse of a polynomial matrix. IMA J. Math. Control Inf. 2001, 18, 83–97. [Google Scholar] [CrossRef]
  31. Penrose, R. A generalized inverse for matrices. Math. Proc. Camb. Philos. Soc. 1955, 51, 406–413. [Google Scholar] [CrossRef] [Green Version]
  32. Hunek, W.P.; Latawiec, K.; Stanisławski, R.; Łukaniszyn, M.; Dzierwa, P. A new form of a σ-inverse for nonsquare polynomial matrices. In Proceedings of the 2013 18th International Conference on Methods and Models in Automation and Robotics (MMAR), Miedzyzdroje, Poland, 26–29 August 2013; pp. 282–286. [Google Scholar]
  33. Hunek, W.P.; Latawiec, K.J.; Majewski, P.; Dzierwa, P. An application of a new matrix inverse in stabilizing state-space perfect control of nonsquare LTI MIMO systems. In Proceedings of the 2014 19th International Conference on Methods and Models in Automation and Robotics (MMAR), Miedzyzdroje, Poland, 2–5 September 2013; pp. 451–455. [Google Scholar]
  34. Majewski, P. Research towards Increasing the Capacity of Wireless Data Communication Using Inverses of Nonsquare Polynomial Matrices. Ph.D. Thesis, Opole University of Technology Press, Opole, Poland, 2017. [Google Scholar]
  35. Dabiri, A.; Butcher, E.A.; Poursina, M.; Nazari, M. Optimal Periodic-Gain Fractional Delayed State Feedback Control for Linear Fractional Periodic Time-Delayed Systems. IEEE Trans. Autom. Control 2018, 63, 989–1002. [Google Scholar] [CrossRef]
  36. Kaczorek, T. Minimum energy control of fractional descriptor positive discrete-time linear systems. Int. J. Appl. Math. Comput. Sci. 2014, 24, 735–743. [Google Scholar] [CrossRef] [Green Version]
  37. Kaczorek, T. Minimum energy control of fractional positive continuous-time linear systems using Caputo-Fabrizio definition. Bull. Pol. Acad. Sci.-Tech. Sci. 2017, 65, 45–51. [Google Scholar] [CrossRef] [Green Version]
  38. Klamka, J. Controllability and Minimum Energy Control Problem of Fractional Discrete-Time Systems. In New Trends in Nanotechology and Fractional Calculus Applications; Baleanu, D., Guvenc, Z.B., Machado, J.A.T., Eds.; Springer: Cham, The Netherlands, 2010; pp. 503–509. [Google Scholar] [CrossRef]
  39. Merrikh-Bayat, F.; Karimi-Ghartemani, M. Method for designing PIλDμstabilisers for minimum-phase fractional-order systems. IET Control Theory Appl. 2010, 4, 61–70. [Google Scholar] [CrossRef]
  40. Yige, Z.; Meirong, X. Stability and stabilization for a class of fractional-order linear time-delay systems. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017; pp. 11340–11344. [Google Scholar] [CrossRef]
  41. Yuan, T.; Zheng, M.; Zhang, K.; Huang, T. Fractional-order PID controllers for stabilization of fractional-order time delay systems based on region stability. In Proceedings of the 2018 Chinese Control and Decision Conference (CCDC), Shenyang, China, 9–11 June 2018; pp. 6633–6638. [Google Scholar] [CrossRef]
Figure 1. Control zeros for fractional-order systems: time horizon k = 100 , y ref ( k + 1 ) = 0 (see Equation (24)).
Figure 1. Control zeros for fractional-order systems: time horizon k = 100 , y ref ( k + 1 ) = 0 (see Equation (24)).
Symmetry 11 01322 g001
Figure 2. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 25 ) , k = 100 , y ref ( k ) = 0 , T-inverse.
Figure 2. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 25 ) , k = 100 , y ref ( k ) = 0 , T-inverse.
Symmetry 11 01322 g002
Figure 3. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 75 ) , k = 100 , y ref ( k ) = 0 , T-inverse.
Figure 3. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 75 ) , k = 100 , y ref ( k ) = 0 , T-inverse.
Symmetry 11 01322 g003
Figure 4. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 T , E 1 = 8 . 18 , E 2 = 1889 . 80 , stable H-inverse with L = 0 . 7231 1 . 7809 .
Figure 4. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 T , E 1 = 8 . 18 , E 2 = 1889 . 80 , stable H-inverse with L = 0 . 7231 1 . 7809 .
Symmetry 11 01322 g004
Figure 5. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 1 = 5 . 75 , E 2 = 468 . 81 , minimum-energy-based H-inverse with L = [ 10.1378 0.2495 ] .
Figure 5. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 1 = 5 . 75 , E 2 = 468 . 81 , minimum-energy-based H-inverse with L = [ 10.1378 0.2495 ] .
Symmetry 11 01322 g005
Figure 6. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 44 . 92 , E 2 = 473 . 40 , stable σ -inverse with β = 1 . 0724 1 . 4588 0 . 9329 1 . 5464 0 . 6466 1 . 9729 .
Figure 6. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 44 . 92 , E 2 = 473 . 40 , stable σ -inverse with β = 1 . 0724 1 . 4588 0 . 9329 1 . 5464 0 . 6466 1 . 9729 .
Symmetry 11 01322 g006
Figure 7. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 6 . 21 , E 2 = 423 . 95 , minimum-energy-based σ -inverse with σ = 243 . 0453 73 . 0769 56 . 3869 890 . 1077 321 . 4291 1545 . 2586 .
Figure 7. Control zeros and time-domain signals, S ( A d , B , C , α = 0 . 01 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 6 . 21 , E 2 = 423 . 95 , minimum-energy-based σ -inverse with σ = 243 . 0453 73 . 0769 56 . 3869 890 . 1077 321 . 4291 1545 . 2586 .
Symmetry 11 01322 g007
Figure 8. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 1 = 216 . 13 , E 2 = 10117 . 60 , stable H-inverse with L = [ 1.8883 2.4348 ] .
Figure 8. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 1 = 216 . 13 , E 2 = 10117 . 60 , stable H-inverse with L = [ 1.8883 2.4348 ] .
Symmetry 11 01322 g008
Figure 9. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 2 = 680 . 61 , minimum-energy-based H-inverse with L = [ 0.7347 2.4624 ] .
Figure 9. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = [ 1 2 ] , E 2 = 680 . 61 , minimum-energy-based H-inverse with L = [ 0.7347 2.4624 ] .
Symmetry 11 01322 g009
Figure 10. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 40 . 56 , E 2 = 1335 . 14 , stable σ -inverse with β = 1 . 5432 0 . 7173 2 . 5449 1 . 5224 0 . 4069 2 . 6763 .
Figure 10. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 40 . 56 , E 2 = 1335 . 14 , stable σ -inverse with β = 1 . 5432 0 . 7173 2 . 5449 1 . 5224 0 . 4069 2 . 6763 .
Symmetry 11 01322 g010
Figure 11. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 10 . 88 , E 2 = 176 . 35 , minimum-energy-based σ -inverse with β = 2 . 5685 0 . 9337 1 . 3965 1 . 7471 0 . 3845 1 . 0463 .
Figure 11. Control zeros and time-domain signals, S ( A d , B , C , α = 1 . 80 ) , k = 100 , y ref ( k ) = 1 2 , E 1 = 10 . 88 , E 2 = 176 . 35 , minimum-energy-based σ -inverse with β = 2 . 5685 0 . 9337 1 . 3965 1 . 7471 0 . 3845 1 . 0463 .
Symmetry 11 01322 g011

Share and Cite

MDPI and ACS Style

Hunek, W.P.; Wach, Ł. A New Stability Theory for Grünwald–Letnikov Inverse Model Control in the Multivariable LTI Fractional-Order Framework. Symmetry 2019, 11, 1322. https://doi.org/10.3390/sym11101322

AMA Style

Hunek WP, Wach Ł. A New Stability Theory for Grünwald–Letnikov Inverse Model Control in the Multivariable LTI Fractional-Order Framework. Symmetry. 2019; 11(10):1322. https://doi.org/10.3390/sym11101322

Chicago/Turabian Style

Hunek, Wojciech Przemysław, and Łukasz Wach. 2019. "A New Stability Theory for Grünwald–Letnikov Inverse Model Control in the Multivariable LTI Fractional-Order Framework" Symmetry 11, no. 10: 1322. https://doi.org/10.3390/sym11101322

APA Style

Hunek, W. P., & Wach, Ł. (2019). A New Stability Theory for Grünwald–Letnikov Inverse Model Control in the Multivariable LTI Fractional-Order Framework. Symmetry, 11(10), 1322. https://doi.org/10.3390/sym11101322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop