Abstract
This paper considers the conversion problem from unstructured linear time-invariant (LTI) controllers to observer-structured LTI controllers, whose structure is similar to but not exactly the same as the so-called “Luenberger observer–based controllers,” for linear parameter-varying (LPV) plant systems. In contrast to Luenberger observer–based controllers, observer-structured LTI controllers can be defined and constructed even if the plant systems are given as LPV systems. In the conversion problem, the state-space matrices of the observer-structured LTI controller are parameterized with those of the given unstructured LTI controller, one free matrix, and a state transformation matrix. We also show a method to obtain the optimal state transformation matrix with respect to the convergence of the discrepancy between the plant state and the observer-structured controller state for a stochastically defined non-zero initial plant state. Several toy examples are included to illustrate the effectiveness and the usefulness of observer-structured LTI controllers, and the utility of the proposed conversion parametrization.
1 Introduction
It is well-known that H∞ control is a powerful tool for controlling plant systems with uncertainties (Zhou et al., 1996). At first, the Riccati equation–based synthesis approach is proposed (Doyle et al., 1989); however, the plant systems should satisfy several assumptions (the so-called “standard assumption”), and this restriction diminishes their applicability and usability. On this issue, after the paper by Sampei et al. (1990) which tackles H∞ controller synthesis with output feedback controllers in terms of linear matrix inequality (LMI), many research studies on LMI-based H∞ controller synthesis have been conducted (Gahinet and Apkarian, 1994; Iwasaki and Skelton, 1994; Masubuchi et al., 1998), and the applicability and usability of H∞ control have been extended. However, it is also well-known that H∞ controllers have a complicated structure; that is, H∞ controllers have no special structures as H2 controllers whose structure is composed of LQR controllers and observers. Due to this property, on-site engineers have a difficulty to understand the structure of H∞ controllers as well as the meaning of the figures of the state-space matrices of the designed H∞ controllers.
On this issue, Alazard has explored a new horizon in the studies of Alazard and Apkarian (1999) and Alazard (2012). In the study by Alazard and Apkarian (1999), for LTI plant systems with no direct feedthrough, the authors have proposed a conversion method from a priori designed unstructured LTI controllers to Luenberger observer–based LTI controllers by finding appropriate state transformation matrices for the controllers. Furthermore, in the study of Alazard (2012), the method in the study of Alazard and Apkarian (1999) is extended to the case in which a direct feedthrough exists. In those papers, the dimensions of the controllers are not restricted to be the same as those of the plant systems; that is, even if the dimensions of the controllers are different from those of the plant systems, appropriate state transformation matrices which render the unstructured LTI controllers to Luenberger observer–based controllers can be found by solving the generalized non-symmetric and rectangular Riccati equation. In summary, LTI controllers (including H∞ controllers) for LTI plant systems can be equivalently represented as Luenberger observer–based controllers using the methods in those papers. This achievement is very helpful for on-site engineers because a priori designed unstructured LTI controllers can be converted to well-known Luenberger observer–based controllers without deteriorating control performance as long as the plant systems are LTI systems. As an application example, the converted controllers can be used as “virtual sensors” (Goupil et al., 2014) for plant health monitoring, fault detection, etc., without any additional controllers or observers because the state of the converted controllers, i.e., Luenberger observer–based controllers, estimates the plant state faithfully. Some application examples can be found in the study of Alazard (2012). By considering the above, it is concluded that the conversion from unstructured LTI controllers to observer-based controllers is useful.
Although the methods in the studies of Alazard and Apkarian (1999) and Alazard (2012) are effective and attractive, there is a drawback that the methods cannot be applied to plant systems which are given as linear parameter-varying (LPV) systems. This is because Luenberger observer–based controllers need the nominal state-space matrices of the plant systems; however, as demonstrated by Peaucelle et al. (2017), Sato (2018), etc., the use of the nominal state-space matrices of plant systems does not always lead to the optimal H∞ control performance. This fact poses a simple question: “Under the supposition that LPV plant systems can be interpreted to be composed of their nominal LTI plant systems and norm bounded variations, do the methods in the studies of Alazard and Apkarian (1999) and Alazard (2012) give the state transformation matrices which minimize the discrepancies between plant systems’ state and converted controllers’ state even for LPV plant systems?” This question motivates us to try to extend the methods in the studies of Alazard and Apkarian (1999) and Alazard (2012) to the case in which the plant systems are given as LPV systems, and also to try to find the counterpart conversion method for LPV plant systems. That is, our addressed problem in this paper is the counterpart problem in the studies of Alazard and Apkarian (1999) and Alazard (2012) for LPV plant systems. To this end, we first propose observer-structured LTI controllers whose structure is similar to but not exactly the same as the so-called Luenberger observer–based controllers. By using the observer-structured controllers, we then propose a method producing appropriate state transformation matrices which convert the a priori designed unstructured LTI controllers to observer-structured LTI controllers. As a consequence, even if plant systems are given as LPV systems, we can use the converted controllers obtained by our proposed method as “virtual sensors” (Goupil et al., 2014) for plant health monitoring, fault detection, etc., without any additional controllers or observers.
This paper is structured as follows: In Section 2, we first give definitions of an LPV plant system with a direct feedthrough and strictly proper LTI controllers (an unstructured controller designed a priori and an observer-structured controller), then show our parameterization of the observer-structured controller with the state-space matrices of the unstructured controller, one free matrix and a state transformation matrix, and finally clarify the relation between our method and the method in the study of Alazard (2012) for the case that plant systems are given as LTI systems with direct feedthrough, and the dimensions of the controllers are the same as those of the LTI plant systems. In Section 3, we propose a method to obtain the optimal state transformation matrix, which gives the minimum convergence of the discrepancy between the plant state and the observer-structured controller state for a stochastically defined non-zero initial plant state, in terms of parameter-dependent linear matrix inequality (LMI). In Section 4, several toy examples are introduced to illustrate our contributions (i.e., proposition and parameterization of observer-structured controllers, and proposition of design method for state transformation matrices minimizing the estimation errors), and finally, we give concluding remarks in Section 5.
We summarize the notation used in this paper. and In, respectively, denote a compatibly dimensional zero matrix and an n × n identity matrix; , , and , respectively, denote the sets of n-dimensional real vectors, n × m real matrices, and n × n positive-definite real matrices; is a shorthand notation of X + XT for a square matrix X; “” in a matrix denotes an abbreviated element by its symmetry; for a square matrix X, denotes its trace. In this paper, we address the continuous-time (CT) case as well as the discrete-time (DT) case simultaneously; therefore, the operator δ [⋅] is used to denote the time-derivative in the CT case and one-step shift operator in the DT case. That is, δ [x] in the CT case denotes for the current time t, and δ [x] in the DT case denotes x(k + 1) for the current step number k. Similarly, Φ in inequalities denotes for the CT case and for the DT case.
2 Proposed Parameterization
2.1 Linear Parameter-Varying Plant System
Let us suppose that the stabilization problem of the following LPV plant system is addressed:where denotes the state, denotes the measurement output, denotes the control input, and matrices A(θ), B(θ), C(θ), and D(θ) are supposed to be compatibly dimensional real matrices which are dependent on parameters . The parameters can be scheduling parameters as well as uncertainty parameters. Here, it is supposed that all parameters are independent from each other. It is also supposed that the existing regions of parameters and their derivatives (CT case) or deviations per single sampling period (DT case) are bounded. Thus, the following are supposed with a priori given polytopes Λi:
Then, the following is also supposed:
2.2 Unstructured Linear Time-Invariant Controller
We next define an unstructured LTI controller which has already been designed:where denotes the state, and matrices , and are constant matrices. Note that the dimension of the controller, i.e., nc, might be different from that of the plant system; that is, it can be possible that nc ≠ n holds in our proposed parameterization.
Then, the closed-loop system comprising G(θ) and CG is straightforwardly derived as follows:
2.3 Observer-Structured Linear Time-Invariant Controller
We next define an observer-structured LTI controller inspired by Sato (2020a). In the DT case, only predictor form (Alazard, 2012) is considered hereafter:where denotes the state, and not only matrices but also matrices , are to be determined.
Then, the closed-loop system comprising G(θ) and CO is straightforwardly derived as follows:
Remark 1As is obvious, the controller does not have the same structure as the so-called “Luenberger observer–based controller” does. However, if there are no parameters in the LPV plant system, that is, the LTI plant system is supposed instead of the LPV plant system , then setting , , , and to be , , , and , respectively, makes the controller to be identical to the conventional “Luenberger observer–based controller.” Therefore, the controller is referred to as the observer-structured controller due to the structural similarity between the observer-structured controller and the “Luenberger observer–based controller.”
2.4 Parameterization of Observer-Structured Linear Time-Invariant Controller
By comparing and under the consideration of the freedom of the state transformation xo = T−1xc with a non-singular matrix T, it is easily confirmed that the two systems are equivalent if the following equations hold:
The last equation in Eq. 8 is equivalently represented as follows:
Then, we give one of our main results.
The matrix satisfying Eq. 9 is parameterized as in Eq. 10 with one free matrix :where matrices with superscript “†” denote the corresponding Moore–Penrose inverse matrices.
Proof 1The corresponding assertion for constant matrices is easily proved by using Theorem 2.3.1 in the study of Skelton et al. (1998), and thus, our claim is also easily proved. However, for completeness, we give the detailed proof.
We now prove that the following two statements are equivalent:
Remark 2As the left-hand sides and the right-hand sides of all equations in Eq. 11 are, respectively, constant and parameter-dependent, in general, complicated constraints consequently arise for Z(θ). A simple solution to escape from the complicated constraints is to set Z(θ) to be constant, i.e., Z(θ) = Z. This setting reduces the generality; however, it also reduces numerical complexity when obtaining the state-space matrices Ao, Bo, Co, and Do in the next section.Note that matrices Ao, Bo, Co, and Do have freedom as in the last equation of Eq. 11; however, they must satisfy Eq. 9, i.e., the last equation of Eq. 8. That is, we would like to emphasize that even if they are chosen as desired matrices within the parameterization of Eq. 11, the input–output property of CO is the same as that of CG. Therefore, the essential freedom of the conversion from CG to CO is only the state transformation represented by T.In the next section, we propose a design method to obtain the optimal state transformation matrix with respect to the convergence of the discrepancy between the plant state and the observer-structured controller state for a stochastically defined non-zero initial plant state, and then we also propose a method to obtain matrices Ao, Bo, Co, and Do as close to as the designated matrices within the parameterization of the last equation of Eq. 11.
2.5 Relation Between Our Method and Existing Method for Linear Time-Invariant Plant Systems
In this subsection, we clarify the relation between our method and the method in the study of Alazard (2012) for LTI plant systems. For simplicity, nc is supposed to be equal to n; that is, only full-order controllers are considered. In the DT case, only predictor form (Alazard, 2012) is considered.
2.5.1 Brief Review of Existing Method in the Study of Alazard (2012)
Now, it is supposed that the plant system is given as an LTI system:where denotes the state, denotes the measurement output, denotes the control input, and matrices A, B, C, and D are supposed to be compatibly dimensional constant matrices.
We define the observer-based controller aswhere denotes the observer state and matrices and are, respectively, the observer gain and the state-feedback gain.
The closed-loop system comprising G and CL is expressed as follows:
In order to make the closed-loop system identical to in Eq. 5 for an LTI plant system G in Eq. 12, the state transformation with a non-singular matrix is considered. Then, the following conditions are straightforwardly derived:If an appropriate is found, then and are, respectively, given as and . Then, the remaining task is to find the state transformation matrix satisfyingi.e.,
The solution of
Eq. 16can be obtained by the following method (
Alazard, 2012):
• Choose n closed-loop eigenvalues which are the eigenvalues of .
• Find n-dimensional invariant subspace, i.e., matrices such that
where the eigenvalues of
are the selected eigenvalues in the previous step.
• If the matrix U1 is confirmed to be non-singular, then the state transformation matrix is obtained as .
By following the procedure for obtaining shown above, the unstructured controller CG can be converted to the Luenberger observer–based controller CL with and .
2.5.2 Comparison of Our Method With Existing Method
Comparing CO in Eq. 6 and CL in Eq. 13 concludes that if matrices Ao, Bo, Co, and Do are set as A, B, C, and D, respectively, and gain matrices K and L are also set as and , respectively, then CO is identical to CL. That is, the observer-structured controller CO encompasses the Luenberger observer–based controller CL.
Furthermore, if the state transformation matrix T is set as in the method of Alazard (2012), then the two closed-loop systems, i.e., and , are also identical to each other. Thus, if the observer-based controller CL can be obtained from the unstructured controller CG by appropriately defined , then it is always possible to have the same controller within the expression of observer-structured controller CO by setting T to be identical to .
In summary, the observer-structured controller encompasses the expression of Luenberger observer–based controller, and the transformation from CG to CO also encompasses the corresponding one from CG to CL. Thus, we conclude that the observer-structured controller is a kind of generalization of the Luenberger observer–based controller, and the conversion from CG to CO is also a kind of generalization of the one from CG to CL.
3 Observer-Structured Controller With Optimal Estimation Error for Non-Zero Initial Plant State
In this section, we give a formulation to obtain the optimal state transformation matrix with respect to the convergence of the discrepancy between the plant state and the observer-structured controller state for a stochastically defined non-zero initial plant state.
We first define performance output to be evaluated. Now, the dimension of the unstructured controller is not always the same as that of the LPV plant system. We thus select nsl state variables of the LPV plant system to be estimated by the observer-structured controller and set the vector containing them as which is expressed as with an appropriately defined constant matrix . Similarly, the corresponding vector in the observer-structured controller is defined as which is expressed as with an appropriately defined constant matrix .
Then, using , we define the following system with performance output z defined as :
Here, .
By considering that if xsl and are close to each other, then it can be concluded that of the observer-structured controller CO plays as a good estimation of xsl of the LPV system G(θ), and we now define the following problem in which the estimation error between the LPV plant state and the converted controller state is minimized.
Problem 1For in Eq. 17, obtain the minimal positive scalar γ and the corresponding non-singular matrix T satisfyingHere, it is supposed that the initial plant state x(0) is set as a stochastic variable which satisfies E[x(0)] = 0 and E[x(0) x(0)T] = In and that the initial controller states xc(0) and xo(0) are both set as 0.To address this problem, we define the following system which is composed of with the performance output z defined in :
Note that the following hold, since xo = T−1xc is applied to to derive Eq. 8:Then, the following theorem is proposed for Problem 1.
If Eq. 21 is solved affirmatively, then calculated with satisfies Eq. 18 with the optimized γ:
Proof 2As the upper-left block of the left-hand side of Eq. 24 is negative-definite, the following are confirmed for all combinations of (θ, δ[θ]) ∈ Λ:Since the positivity of is ensured in Eq. 23, the stability of is confirmed.Note that Eq. 24 is equivalent to the following for all combinations of (θ, δ[θ]) ∈ Λ:After is set, multiplying and its transpose to Eq. 26 from the left and the right, respectively, leads to the following inequality for all combinations of (θ, δ[θ]) ∈ Λ:Then, the satisfaction of Eq. 18 is confirmed after consideration of the stability of and Eq. 22.
Remark 3If is singular, then a small perturbation to is to be conducted to obtain non-singular , similarly to Masubuchi et al. (1998). In particular, if is set as a row full-rank rectangular matrix, then has freedom, which has no effect on the performance index γ, and thus, there is a possibility to have singular . In this case, a small perturbation to is necessary to obtain non-singular .In general, inequality Eq. 24 is not parametrically affine; thus, some relaxation methods (e.g., sum-of-squares (Chesi et al., 2003; Parrilo, 2003), Polya’s theorem–based relaxation (de Oliveira and Peres, 2007), and slack-variable approach (Peaucelle and Sato, 2009)) should be applied to be solved numerically. However, such relaxations usually increase the numerical complexity. On this issue, if the state-space matrices of the LPV plant G(θ) are supposed to be multi-affine with respect to parameters (Amato et al., 2005; Peacelle and Ebihara, 2018; Sato, 2020b) and the parameter-dependent decision matrix in Theorem 2, i.e., , is also set as multi-affine with respect to parameters, then the multi-affine property [i.e., Lemma 3 in the study of Sato (2020b)] can be used for Theorem 2 after extended/dilated LMI technique (de Oliveira et al., 1999; Peaucelle et al., 2000; Pipeleers et al., 2009) is applied to maintain the numerical complexity as small as possible while conservatism in solving LMIs is also kept small.To this end, the following assumption is now made.
Assumption 1All the parameter-dependent state-space matrices of G(θ) in Eq. 1 are supposed to be multi-affine with respect to θi; that is, they are represented as follows:where is the coefficient matrix.Then, the following lemma is straightforwardly derived from Theorem 2 with the use of extended/dilated LMI technique (de Oliveira et al., 1999; Peaucelle et al., 2000; Pipeleers et al., 2009) and multi-affine property [i.e., Lemma 3 in the study of Sato (2020b)].
Lemma 1Under Assumption 1, if Eq. 29 is solved affirmatively with multi-affine , then calculated with (if necessary, a small perturbation in Remark 3 is conducted) satisfies Eq. 18 with the optimized γ:In inequalities (Eqs 30–32), denotes the vertex set.The proof is omitted as it is straightforward in consideration of the following relation:If the initial plant state is given as a deterministic (but unknown) variable instead of a stochastic variable, and its region is supposed to be in a polytope , then the following method is applicable instead of Lemma 1.
Lemma 2Under Assumption 1, if Eq. 33 is solved affirmatively with multi-affine , then calculated with (if necessary, a small perturbation in Remark 3 is conducted) satisfies Eq. 18 with the optimized γ for all :By using Lemma 1 or Lemma 2, we can obtain an optimal constant state transformation matrix T for Problem 1. Although the essential freedom for the conversion from CG to CO is the state transformation matrix T, we still have freedom for the state-space matrices (i.e., Ao, Bo, Co, and Do) in an observer-structured controller CO. To address this issue, we propose the following problem to obtain the state-space matrices Ao, Bo, Co, and Do as close to as the designated ones after obtaining T. In our proposed optimization problem shown below, we set Z(θ) in Eq. 10 to be constant Z by considering Remark 2.
Problem 2If Eq. 35 is solved affirmatively, then all the elements of the state-space matrices in the observer-structured controller CO satisfy .where ϒi,j represents the (i, j) element of the right-hand side of the last equation of Eq. 11 with constant , i.e., represents its designated figure assigned by the designer, and Wi,j represents the weighting coefficient for .In the problem above, is overbounded by η/Wi,j. Similarly, it can be possible to overbound the sum of , i.e., .We would like to emphasize that even if ϒi,j can be set as the same figure as , the state-space matrices of CO satisfy Eq. 9. Thus, the input–output property of CO is the same as that of CG. This will be confirmed by toy examples in the next section.
4 Numerical Examples
Several toy examples are introduced to clearly illustrate our contributions. To this end, we first confirm that our conversion method from CG to CO, i.e., Eq. 11, encompasses the method from CG to CL proposed in the study of Alazard (2012) for LTI plant systems and full-order unstructured LTI controllers. We next show the effectiveness of our method for state transformation matrices, i.e., Lemma 2, and our method for obtaining a priori designated matrices as the state-space matrices in CO, i.e., solving Eq. 35 in Problem 2. We finally show the effectiveness of Lemma 1 in CT and DT cases.
4.1 Confirmation of Relation Between the Method in the Study of Alazard (2012) and Our Method
We first confirm the relation of our method and the method in the study of Alazard (2012) using the following CT LTI plant and CT unstructured LTI controller. The controller is designed to assign the closed-loop poles at −7, −5, −3 ± i; that is, a stabilizing controller is designed:
.
By following the procedure in the study of Alazard (2012), if the poles −3 ± i are set to be the eigenvalues of , then the following state transformation matrix, and observer and state-feedback gains are correspondingly obtained:The state-space representation of CL is consequently obtained as follows:
We now obtain K and L in Eq. 11 by setting . Then, the observer gain and the state-feedback gain in CO are obtained as follows:
Obviously, they are the same as the result obtained by the method in the study of Alazard (2012) because T is set as . That is, it is confirmed that our method is identical to the method in the study of Alazard (2012).
Next, we obtain matrices Ao, Bo, Co, and Do in Eq. 11 with several Z, i.e., constant Z due to an LTI plant system and LTI controller. To this end, and are, respectively, set as and , and matrix Z is set as Z = 100 × I3 (case a), Z = 0 (case b), and (case c). Then, the matrix XO is obtained as follows:
Although the matrix XO above has corresponding figures in accordance with different Z, Ao + BoK + LCo + LDoK is calculated as in all three cases. It is also confirmed that this matrix is the same as T−1AcT.
That is, we have confirmed that our parameterization encompasses the formulation in the study of Alazard (2012) regardless of the choice of Z.
4.2 Effectiveness Demonstration of Our Method for Linear Time-Invariant Plant System
We next consider the same example in the study of Alazard (2012), i.e., the following CT LTI plant:
.We now suppose that the following CT unstructured LTI stabilizing controller, which is also borrowed from the study of Alazard (2012), has already been designed:
.(36)
This controller assigns the closed-loop poles at −3, −4, −10 ± 10i. We now suppose that poles −3 and −4 are assigned by . Then, the state transformation matrix and observer and state-feedback gains are correspondingly obtained as follows:
The consequently calculated CL is given as
.(37)
We now solve Eq. 33 in Lemma 2 for and then the state transformation matrix T is obtained as with γ = 2.154. The numbers of LMI rows and decision variables in LMIs are 20 and 47, respectively. The correspondingly calculated CO with Z being 0 is shown as
(38)
We show the simulation results using controllers CG, CL, and CO in Figure 1 and show the finite-time state estimation performance, i.e., , in Table 1. Although the estimation of the first state of the plant system by CL is faithful, the second state of CL initially moves in the opposite direction to the second state of the plant system; in this sense, the estimation of the second state by CL is poor.
FIGURE 1
In contrast, the observer-structured controller CO has a much better state estimation performance as indicated in Table 1. This clearly illustrates the effectiveness of our method compared to the method in the study of Alazard (2012).
We next address Problem 2 with the following ϒdes and W using the above T:where ∗ denotes the element of no interest. Solving Eq. 35 gives the following Z and the following corresponding state-space matrices of CO:
(40)
The transfer functions of CG in Eq. 36, CL in Eq. 37, and CO in Eq. 38 are all calculated as . The transfer function of CO with Ao, Bo, Co, and Do in Eq. 40 is also calculated as ; however, the state-space matrices of CO can have an artificially assigned special structure as in Eq. 40.
4.3 Effectiveness Demonstration of Our Method for Linear Parameter-Varying Plant System
We finally consider the conversion problem for LPV plant systems. Note that the methods in the studies of Alazard and Apkarian (1999) and Alazard (2012) cannot be applied to the example shown below because the plant system is given as an LPV system.
4.3.1 CT Case
We first consider the following CT LPV system:
(41)where . We set that the two parameters are supposed to be frozen in the interval [0.9, 1.1]; that is, the vertex sets of Λ1 and Λ2 are set as follows:
It is supposed that the following stabilizing controller is given as an unstructured controller:
(42)
For this problem setup, we solve Eq. 29 in Lemma 1 with Csl = [1 0] and and obtain T = 0.333 with γ = 0.052. The numbers of LMI rows and decision variables in LMIs are 44 and 44, respectively. The state-space representation of CO with Z = 0 is consequently given as follows:
(43)
The simulation results with CG in Eq. 42 and CO in Eq. 43 are shown in Figure 2. It is confirmed that the state of CO faithfully represents the plant state x1 in both cases with x(0) = [1 0]T and x(0) = [0 1]T.
FIGURE 2
We next address Problem 2 with the following ϒdes and W using the above T:where ∗ denotes the element of no interest. Solving Eq. 35 gives the following Z and the following corresponding state-space matrices of CO:
.(45)
The transfer functions of CG in Eq. 42 and CO in Eq. 43 are both calculated as . The transfer function of CO with Ao, Bo, Co, and Do in Eq. 45 is also calculated as ; however, the state-space matrices of CO can have an artificially assigned special structure as in Eq. 45.
4.3.2 DT Case
We next consider the discretized system of Eq. 41 by using Euler approximation with the sampling period ΔT = 0.1 (s):
.(46)
We now set that θ1 can vary arbitrarily fast within the interval; however, θ2 is set as frozen in the interval; that is, the following vertex sets are considered:
The controller in Eq. 42 is also discretized by using Euler approximation with the sampling period ΔT = 0.1 (s):
(47)
For this problem setup, we solve Eq. 29 in Lemma 1 with Csl = [1 0] and and obtain T = 0.353 with γ = 1.179. The numbers of LMI rows and decision variables in LMIs are 72 and 44, respectively. The state-space representation of CO is consequently given as follows:
.(48)
The simulation results with CG in Eq. 47 and CO in Eq. 48 are shown in Figure 3. It is indeed better to have simulations with varying θ1 within the interval [0.9, 1.1]; however, there are uncountable combinations for the variation of θ1. Thus, we conduct numerical simulations only with the fixed extreme points of θ1. It is confirmed that the state of CO faithfully represents the plant state x1 in both cases with x(0) = [1 0]T and x(0) = [0 1]T.
FIGURE 3
We next address Problem 2 with the following ϒdes and W using the above T:
In this case, we set ϒdes composed of the upper-left element of A(θ), the first element of B(θ), and the first element of C(θ) and D(θ) in the state-space matrices of Eq. 46. Our aim is to obtain matrices close to the nominal state-space matrices corresponding to the first-input-first-output system. We minimize the Frobenius norm for ϒ−ϒdes, and obtain the following Z and the corresponding state-space matrices:
(50)
In this case, in contrast to the previously shown examples, we cannot obtain the designated figures for the state-space matrices Ao, Bo, Co, and Do due to lack of enough freedom in Eq. 11.
The transfer functions of CG in Eq. 47 and CO in Eq. 48 are both calculated as . The transfer function of CO with Ao, Bo, Co, and Do in Eq. 50 is also calculated as ; however, the state-space matrices of CO can have very close figures as assigned in Eq. 49.
5 Conclusion
We address the conversion problem from unstructured LTI controllers to observer-structured LTI controllers, whose structure is similar to but not exactly the same as Luenberger observer–based controllers, for LPV systems with direct feedthrough. To this end, we first define observer-structured LTI controllers, then parameterize the state-space matrices with a priori designed unstructured LTI controller, one free matrix and a state transformation matrix, and finally propose a method which produces the optimal state transformation matrix with respect to the convergence of the discrepancy between the plant state and the observer-structured controller state for a stochastically defined non-zero initial plant state. Several toy examples are introduced to clearly illustrate the effectiveness and usefulness of observer-structured LTI controllers and the proposed method for obtaining optimal state transformation matrices with respect to the minimization between the discrepancies between the LPV plant state and the converted LTI controller state.
In this paper, we address the controller conversion problem in the case that only the stabilization problem of LPV plant systems is considered. As our next step, we are now tackling the same conversion problem in the case that some control performance criteria are also considered. Then, we will demonstrate the practicality of the conversion by using practical systems including H∞ performance.
In this paper, it is also supposed that LTI controllers are given for LPV plant systems; however, the usefulness and effectiveness of using LPV controllers for LPV plant systems are also well recognized. Thus, the extension of our results to LPV controllers is another future research topic.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
MS was in charge of theoretical development, numerical examples, and writing. NS supervised Sato’s work.
Funding
This work was conducted under the support of KAKENHI Grant Number 22K04176.
Acknowledgments
The authors give a sincere acknowledgement to Ms. Shikada for checking the numerical example results.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AlazardD. (2012). Advances in Missile Guidance, Control, and Estmation. Boca Raton, FL: CRC Press. Chap. Control Designand Gain-Scheduling Using Observer-Based Structures. 77–128.
2
AlazardD.ApkarianP. (1999). Exact Observer-Based Structures for Arbitrary Compensators. Int. J. Robust Nonlinear Control.9, 101–118. 10.1002/(sici)1099-1239(199902)9:2<101::aid-rnc400>3.0.co;2-u
3
AmatoF.MatteiM.PirontiA. (2005). Gain Scheduled Control for Discrete-Time Systems Depending on Bounded Rate Parameters. Int. J. Robust Nonlinear Control.15, 473–494. 10.1002/rnc.1001
4
ChesiG.GarulliA.TesiA.VicinoA. (2003). Solving Quadratic Distance Problems: An LMI-Based Approach. IEEE Trans. Automat. Contr.48, 200–212. 10.1109/TAC.2002.808465
5
de OliveiraR. C. L. F.PeresP. L. D. (2007). Parameter-Dependent LMIs in Robust Analysis: Characterization of Homogeneous Polynomially Parameter-Dependent Solutions via LMI Relaxations. IEEE Trans. Automat. Contr.52, 1334–1340. 10.1109/TAC.2007.900848
6
de OliveiraM. C.BernussouJ.GeromelJ. C. (1999). A New Discrete-Time Robust Stability Condition. Syst. Control. Lett.37, 261–265. 10.1016/S0167-6911(99)00035-3
7
DoyleJ. C.GloverK.KhargonekarP. P.FrancisB. A. (1989). State-Space Solutions to Standard H/Sub 2/and H/Sub Infinity/Control Problems. IEEE Trans. Automat. Contr.34, 831–847. 10.1109/9.29425
8
GahinetP.ApkarianP. (1994). A Linear Matrix Inequality Approach to H∞ Control. Int. J. Robust Nonlinear Control.4, 421–448. 10.1002/rnc.4590040403
9
GoupilP.DayreR.BrotP. (2014). “From Theory to Flight Tests: Airbus Flight Control System TRL5 Achievements,” in Proceedings of the 19th IFAC World CongressCape Town, South Africa, 10562–10567. IFAC Proceedings Volumes. 10.3182/20140824-6-za-1003.0254747
10
IwasakiT.SkeltonR. E. (1994). All Controllers for the General ∞ Control Problem: LMI Existence Conditions and State Space Formulas. Automatica30, 1307–1317. 10.1016/0005-1098(94)90110-4
11
MasubuchiI.OharaA.SudaN. (1998). LMI-Based Controller Synthesis: A Unified Formulation and Solution. Int. J. Robust Nonlinear Control.8, 669–686. 10.1002/(sici)1099-1239(19980715)8:8<669::aid-rnc337>3.0.co;2-w
12
ParriloP. A. (2003). Semidefinite Programming Relaxations for Semialgebraic Problems. Math. Programming96, 293–320. 10.1007/s10107-003-0387-5
13
PeacelleD.EbiharaY. (2018). “Affine versus Multi-Affine Models for S-Variable LMI Conditions,” in Preprints of the 9th IFAC Symposium on Robust Control DesignFlorianopolis, Brazil, 643–648.
14
PeaucelleD.ArzelierD.BachelierO.BernussouJ. (2000). A New Robust - stability Condition for Real Convex Polytopic Uncertainty. Syst. Control. Lett.40, 21–30. 10.1016/S0167-6911(99)00119-X
15
PeaucelleD.EbiharaY.HosoeY. (2017). Robust Observed-State Feedback Design for Discrete-Time Systems Rational in the Uncertainties. Automatica76, 96–102. 10.1016/j.automatica.2016.10.003
16
PeaucelleD.SatoM. (2009). LMI Tests for Positive Definite Polynomials: Slack Variable Approach. IEEE Trans. Automat. Contr.54, 886–891. 10.1109/TAC.2008.2010971
17
PipeleersG.DemeulenaereB.SweversJ.VandenbergheL. (2009). Extended LMI Characterizations for Stability and Performance of Linear Systems. Syst. Control. Lett.58, 510–518. 10.1016/j.sysconle.2009.03.001
18
SampeiM.MitaT.NakamichiM. (1990). An Algebraic Approach to H∞ Output Feedback Control Problems. Syst. Control. Lett.14, 13–24. 10.1016/0167-6911(90)90075-6
19
SatoM. (2020a). Observer-Based Robust Controller Design with Simultaneous Optimization of Scaling Matrices. IEEE Trans. Automat. Contr.65, 861–866. 10.1109/TAC.2019.2919811
20
SatoM. (2020b). One-Shot Design of Performance Scaling Matrices and Observer-Based Gain-Scheduled Controllers Depending on Inexact Scheduling Parameters. Syst. Control. Lett.137, 104632. 10.1016/j.sysconle.2020.104632
21
SatoM. (2018). Simultaneous Design of Discrete-Time Observer-Based Robust Scaled-H8 Controllers and Scaling Matrices. SICE J. Control Meas. Syst. Integration11, 65–71. 10.9746/jcmsi.11.65
22
SkeltonR. E.IwasakiT.GrigoriadisK. (1998). A Unified Algebraic Approach to Linear Control Design. London, UK: Taylor & Francis.
23
ZhouK.DoyleJ. C.GloverK. (1996). Robust and Optimal Control. Upper Saddle River, New Jersey: Prentice-Hall.
Summary
Keywords
observer-based controller, observer-structured controller, LTI controller, conversion, LPV system
Citation
Sato M and Sebe N (2022) Conversion From Unstructured LTI Controllers to Observer-Structured Ones for LPV Systems. Front. Control. Eng. 3:804549. doi: 10.3389/fcteg.2022.804549
Received
29 October 2021
Accepted
26 January 2022
Published
29 August 2022
Volume
3 - 2022
Edited by
Olivier Sename, Grenoble Institute of Technology, France
Reviewed by
Natalia Alexandrovna Dudarenko, ITMO University, Russia
Updates
Copyright
© 2022 Sato and Sebe.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Masayuki Sato, sato.masayuki@jaxa.jp
This article was submitted to Control and Automation Systems, a section of the journal Frontiers in Control Engineering
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.