Control for a Class of Unstable High-Order Systems with Time Delay Based on Observer–Predictor Approach

: This work considers the stabilization of a high-order system with time delay; an observer– predictor scheme is designed to estimate an internal signal of the system that is not available for measurement: this internal signal is the output before being delayed. By using the estimated signal, it is possible to design a controller for the delay-free system. The key point to carrying out this estimation strategy is to obtain conditions assuring that the estimated signal converges to the internal variable of the system. A necessary and sufﬁcient condition to achieve an appropriate convergence in the proposed observer–predictor scheme is given. In addition, an analysis of the disturbance rejection and robustness with respect to the delay term is provided. The correct functioning of this scheme is veriﬁed through an example.


Introduction
Time delay systems are often found in dynamic processes.It is well known that when the delay is small enough in comparison with the dynamics of the process, it can be neglected.However, if this is not the case, a careful design of the control stage becomes necessary.Systems with a delay are common in different fields of application, for example, in chemical processes [1], satellite communications [2], hydraulic systems [3], economic systems [4] and technological processes, such as remote operation [5,6].From the control point of view, linear delayed systems pose a challenge because when closing the control loop, the delay term is located in the denominator of the transfer function.As a consequence, the characteristic equation has an infinite number of roots.
On the other hand, when a digital control is implemented, the process to be controlled should be discretized.In the particular case of a discretization by using a zero-order hold (ZOH) of the delayed system, the delay term yields n poles at the origin, where n is an integer and must satisfy T = τ/n, where T is the sampling period and τ is the time delay.This approach allows designing a digital controller using any classical technique in the discrete domain (for instance, P, PD, PI, or PID controllers).However, some robustness problems can occur when the delay of the system differs from the model.Another way to use a discrete model of the system, as mentioned above, is to perform a discrete analysis of the transfer function together with the controller by leading the sampling time T → 0 in the analysis (or equivalently n → ∞) to get stability conditions in continuous time [14,15].From a different perspective, the control of the delayed system has been treated by using the well-known Smith Predictor (PS) [16].This strategy estimates a signal before being delayed, and later this estimated signal is used to apply the control to the delay-free system.However, the main limitation of this strategy is that it only allows the treatment of stable systems.Therefore, numerous modifications have been developed in the structure of the PS.Some modifications of the PS for unstable systems with delay are presented by [17][18][19][20][21][22][23].Recently, in [24], an observer-predictor is presented for the case of one, two, and even three unstable poles with time delay.In particular, the stabilization for a system with one unstable pole, several stable poles, and time delay, the same kind of plant as that analyzed in the present work, is discussed.Stability conditions with respect to the time delay are given based on the use of an observer-predictor scheme and a PID-type controller.
In this work, an observer-predictor control scheme for high-order delayed systems with an unstable pole and several stable poles is proposed.For the convergence of the observer-predictor, a special injection scheme is considered, which is discretized, and later the stability condition of the system in continuous time is obtained under the consideration of a sampling period T = τ n , with T → 0 (or equivalently when n → ∞).Special attention is given to the case of a delayed unstable second order (one stable and one unstable pole) because it is a system commonly treated in the literature [21,[25][26][27].Later, the proposed control strategy is generalized to high-order unstable systems with delay.The proposal presented here improves the size of the delay that can be controlled compared to some related works found in the literature, such as [14,24,27,28].It is important to highlight that the observer-predictor structure proposed in the present work allows stabilizing the kind of systems previously mentioned with a larger delay than the observer-predictor control proposed in [11,24].This improvement is achieved due to the fact that the new proposed observer-predictor in the present work has a different structure, as well as the fact that it has two free parameters for the design of the observer.Moreover, necessary and sufficient conditions are presented for the proposal in contrast with some strategies found in the literature, where necessary and/or sufficient conditions are not derived [19,29].
This paper is organized as follows.Section 2 describes the kind of system to be studied.In Section 3, the necessary preliminary results are shown; these results will be used later in the main results of this paper.Section 4 presents the main results, that is, the estimation strategy and the design of the control strategy.In Section 5, the step tracking and step disturbance rejection of the proposed scheme is analyzed.In Section 6, a robustness analysis with respect to uncertainty in the time delay is presented.In Section 7, three simulation examples implemented in Matlab-Simulink software illustrate the performance of the proposed strategy and finally, conclusions are presented in Section 8.

The Class of Systems
Consider the class of linear Single-Input Single-Output (SISO) system with input/output delay given with where Let us define the product as , and ) and Y(s) are the input and output signals respectively, τ > 0 is the constant known delay, a > 0 is the position of the unstable pole, b i > 0, for i = 1, 2, • • • , p; b > 0 is the gain of the system and G(s) is the delay-free transfer function.Note that in system (1) only real poles are considered.
A traditional control strategy based on a unitary feedback output of the form, produces a closed-loop system given by, where the term e −τs located in the denominator of the transfer function in Equation ( 4) leads to a system with an infinite number of poles and, therefore, the stability analysis is complicated to carry out, and the design of the controller should be carefully analyzed.

Preliminary Results
Next, the stability condition for the system (1) using a PID controller is given.
Lemma 1 ([11]).Consider the system given in (1), the PID type controller given by and a control law U(s) = [R(s) − Y(s)]C(s) with R(s) as a new input reference.Then, there exist k p , k i , and k d such that the closed-loop system Y(s)/R(s) is stable if and only if Lemma 2. Consider the subsystem G(s) = Ȳ(s)/ Ū(s) in (1) (delay-free transfer function), a PID controller (5) and a control law Then, there exist k p , k i , and k d such that the closed-loop system Y(s)/R(s) is stable if and only if Proof.This result is a particular case of Lemma 1, by considering a delay term equal to zero.
In this paper, it is assumed that system (1) can be stabilized using a PID controller for τ = 0, i.e., the condition given in Lemma 2 is satisfied.

Main Results
Before providing the estimation strategy, a result that establishes the stability conditions for the injection scheme shown in Figure 1 is presented.This result will be used later in the proof of the convergence of the proposed observer.Lemma 3. Consider the system given by Equation ( 1) with an output injection described by Figure 1.Then, there exist gains g 1 and g 2 such that the closed-loop system, Proof.The proof of this lemma is presented in Appendix A.
For future reference, consider a state space representation of the system (1) with G s(s) = 1, given as Therefore, a closed-loop state space representation (related to the closed-loop transfer function ( 6)) considering G s(s) = 1 of the output injection shown in Figure 1 can be written as On the other hand, consider a state space representation of the system (1) with G s(s) considering p ≥ 2 given with Thus, by using (9), a closed-loop state space representation of the closed-loop transfer function (6) or equivalently the output injection shown in Figure 1 can be written as Remark 1.Notice that the closed-loop state space representations (8) and (10) have the same stability condition established in Lemma (3) due to the fact that (8) is a particular case of (10) and by the assumption b i > 0, for i = 1, 2, • • • , p previously made.
The parameters of the injection scheme depicted in Lemma 3 are presented in the following result.

Corollary 1.
Let us consider the system given by Equation (1) with an output injection as in Figure 1.The parameters g 1 and g 2 such that the closed-loop system Y(s)/U 1 (s) (Equation ( 6)) is stable can be obtained as follows: for some 1 , 2 > 0.
Proof.The proof of this corollary is presented in Appendix B.

Estimation Strategy
The proposed observer-predictor scheme for system (1) is presented in Figure 2. The following results present the convergence conditions for two cases: (i) unstable delayed second order (G s(s) = 1).(ii) unstable delayed higher order (G s(s) with p ≥ 2).A state space representation of system (1) considering G s(s) = 1 is given with Theorem 1.Consider system (1) with G s(s) = 1 with the state space representation (13) and the observer-predictor scheme shown in Figure 2.Then, there exist constants g 1 and g Proof.Consider the observer-predictor scheme of Figure 2. Its complete dynamics can be described as Defining the prediction error as e ω (t) = ω(t) − ω(t), where e ω (t) = e ω1 (t) e ω2 (t) and the estimate of the output error as e y (t) = ŷ(t) − y(t), it is possible to describe the behavior of the dynamic error signal error ėω Whether e ω (t) → 0 as t → ∞ in ( 16) depends on the stability of the dynamics depicted by (16).Notice that if the dynamic error ( 16) is stable and the initial conditions of the error e ω (0) are different from zero, the corresponding trajectories of the error e ω (t) → 0 as t → ∞.In this way, consider now a state-space realization of the system given by (8).It is clear that the stability condition for system (8) given in Lemma 3 (see Remark 1) is the same as the prediction error given in ( 16) since both systems have equivalent dynamic behaviors.Therefore, the condition for the convergence of the observer is τ < Now, to extend these results about the design of the observer to a more general class of systems, consider an unstable delayed higher order system given by (1) (G s(s) with p ≥ 2) with its state space representation given by where , C 1 = 1 0 ... 0 0 , and The following result presents how to extend the result of the proposed observer scheme in Theorem 1 to a system given by (1).Theorem 2. Let us consider the high-order unstable system in (1), its state space representation (17), and the observer scheme given in Figure 2.Then, there exist constants g , with and the estimate of ω(t) as, Proof.The proof of this theorem is similar to the proof of Theorem Since the behavior of the dynamic error signal ( 20) is equivalent to the dynamic system (10), the stability condition is also equivalent.Thus, there are gains g 1 and

Controller Design
Previously, the convergence of the prediction stage was ensured.Thus, designing the controller based on the estimated variables is possible.Therefore, in what follows, the stability conditions based on the observer for the closed-loop system are stated.Theorem 3. Let us consider the system (13) or (17) and the observer scheme presented in Figure 2 with an adequate convergence (Theorems 1 and 2, respectively).Then there is a control C(s) of type PID such that the closed loop system Y(s Proof.The adequate convergence of the observer is provided by Theorems 1 and 2: for systems (13) or (17), respectively.Additionally, by using the separation principle for linear observers, a controller can be designed by using the estimated variable ŷ(t + τ) as if the variable were taken from the original system.Therefore, a PID controller C(s) stabilizing the transfer function without delay, G(s), for both cases (13) or ( 17) can be found (see Lemma 2).Therefore the most restrictive condition between the convergence and control stages is the convergence stage, from which the conclusion of the theorem follows.

Improved Stability Conditions
Considering a partition of the time delay, it is possible to improve the stability conditions for some cases, [14].In this way, the result presented in the previous sections can be improved as follows.Consider the system (1) with the time delay partition given by, The observer scheme given in Figure 2 can be modified as shown in Figure 3.

Lemma 4.
Let us consider the system (21), G s(s) = 1 and the observer scheme presented in Figure 3.Then, there exist gains g 1 , g 2 , and a PID type controller C(s) such that the closed-loop Proof.Note that Theorem 1 can be applied directly to a system with delayed input u(t − τ 1 ) since the convergence of the observer is not affected by adding the delay term e −τ 1 s to the observer stage shown in Figure 3.That is, there are gains g 1 and g 2 such that . Invoking the principle of separation for linear observers, and based on the adequate convergence of ω(t), Lemma 1 can be applied considering the estimated signal ŷ(t + τ 2 ) to the design of the control stage.That is, there exists a control C(s) of type PID that stabilizes the system with the delay τ 1 if and only if τ . Therefore, it can be concluded that the original time Figure 3. Proposed observer-predictor scheme for longer delays.
Corollary 2. Let us consider the system (21), G s(s) with p ≥ 2, and the observer scheme presented in Figure 3.There are gains g 1 and g 2 and a control C(s) of type PID such that the closed-loop system Y(s Proof.Theorem 1 can be applied directly to the observer of Figure 3 since the convergence of the observer is not affected by the delayed input u(t − τ 1 ).Thus, there are gains g 1 and . Taking into account the principle of separation for linear observers and based on the convergence of ŷ(t + τ 2 ), Lemma 1 can be applied to the control stage by considering the estimated signal ŷ(t + τ 2 ).Then, there is a PID type control C(s) that stabilizes the plant G(s) with a delay . Therefore, it can be concluded that the original delay time limit is Remark 2. Since, in Lemma 4, the size of the delays τ 1 and τ 2 is bounded in the same proportion by convergence and control conditions, it is suggested to set the partition as τ 1 = τ/2 and τ 2 = τ/2.For the case of Corollary 2, it is suggested to set the delay partition as τ 1 < τ 2 since the proportion of the bounds in convergence and control conditions is different.

Disturbance Rejection and Tracking Reference
The control strategy for disturbance rejection and tracking reference is analyzed in this section.Note that once an adequate estimation is achieved, there can be proposed any control strategy for the delay-free system (1) (Theorem 3) or with a part of the time-delay (Lemma 4 and Corollary 2).To make the controlled system follow references and reject disturbances, both step types, the following PID control is used, where k p , k i , and k d are the tuning parameters of the controller.The control parameters in (22) are obtained considering the system G(s) (Theorem 3), using a root locus diagram for G(s)C(s).It is recommended to locate the controller's zeros in the left half-plane to ensure the existence of a region of stability.The control parameters can be obtained by using the root locus analysis and considering a damping factor of ζ ≥ 1 as a specification performance of the delay-free process transfer function.

Tracking Reference
The following result guarantees that the controller is able to generate a response allowing reference tracking.
Proof.The transfer function in Figure 2 Y(s)/R(s) is given by Applying the final value theorem to (23), one obtains that That is, the scheme shown in Figure 2 is capable of following a step type reference.

Disturbance Rejection
The following result addresses the problem of rejecting disturbances.Lemma 6.Let us cosider system (1 ; the predictor scheme given by Figure 2, C(s) described by (22), the input reference R(s) = 0 and H(s) = β s a step-type disturbance.Then, lim t→∞ y(t) = 0.
Proof.The transfer function in Figure 2 Y By applying the final value theorem to (25), one obtains This shows that the scheme shown in Figure 2 is able to reject step-type disturbances.
Note that the design of the controller C(s) involved in the control strategy developed in Theorem 3 is carried out by using the delay-free system G(s).However, in the proposals of Lemma 4 and Corollary 2, the controller C(s) must be carried out considering the delayed system G(s)e −τ 1 s .Therefore, for the latter case, the parameters of the PID controller corresponding to C(s) can be calculated using results related to Lemma 1. Furthermore, reference tracking and disturbance rejection issues, both of the step type, can be proven for the control strategies proposed in Lemma 4 and Corollary 2 by using a similar development as that presented in the proofs of Lemmas 5 and 6.

Robustness with Respect to Uncertainty in the Time Delay
In practice, a control strategy that provides stability with respect to the uncertainties of the model is desired.In this section, uncertainty with respect to the time delay is considered for the proposed control strategy.
Following [30], the robustness properties of the control strategy with respect to the time delay are analyzed.For this purpose, let us consider a characteristic quasi-polynomial of the form p(s) = p 0 (s) + p 1 (s)e −τs + p 2 (s)e −τ 0 s , (27) where its stability properties are established as a function of the delay times τ and τ 0 .Using the results presented in [30], it is possible to give a general framework for our particular case.Let T be the set of all points (τ, τ 0 ) ∈ R 2 + such that p(s) (Equation ( 27)) has at least one zero on the imaginary axis.Any (τ, τ 0 ) ∈ T is known as a crossing point, and T is the collection of all stability crossing curves.Consider now system (1) and the proposed scheme shown in Figure 2 with G s(s) = 1 and C(s) given by ( 22), τ being the delay time at the observer, and τ 0 being the delay time of the process.The closed-loop characteristic equation is given by with Note that the characteristic Equation ( 28) has form (27), so it is possible to identify the regions of (τ, τ 0 ) ∈ R 2 + such that p A (s) is stable.Figure 4 shows a numerical example of the stable region (τ, τ 0 ) for the characteristic Equation ( 28).This figure illustrates the interval of values τ 0 such that the proposed observer-based controller remains stable for a nominal delay τ at the observer stage.
The GP scheme has performed better than many previous approaches, so the observer-predictor scheme presented in this work is compared with the GP scheme.For this example, a = 0.0096, b = 0.033, and τ = 20 are the parameters of the process.Thus, using the observer-predictor scheme presented in this work, it is possible to stabilize the process given by ( 29) due to the condition of Theorem 1 being satisfied.Using the Corollary 1, the value g 1 = 0.7 is chosen; this value guarantees the estimation of the signal of interest.Additionally, to improve the output response, a two-degree-of-freedom controller is used, which is defined by s + βk d s, with k p = 14.679, k d = 6.99, k i = 1.398, α = 0.6, and β = 0.9.
For the simulation, a step-type input R(s) of amplitude 1, initial conditions for the plant state of 0.01 and a perturbation acting at 250s are considered.Figure 5 shows the output response of both strategies considering white noise effects; it is observed that both strategies show a stable, steady-state output.It is also important to note that the scheme proposed in this work improves the time to reach a steady state and shows better performance with respect to disturbance rejection.In Figure 6, the behavior of the strategy reported in [22] and the proposed observerpredictor scheme under uncertainties in the plant are shown; the uncertainty is set at the unstable pole of the system with +200% of its nominal value.
Finally, Table 1 presents the evaluation results for different performance indicators: the integral quadratic error criterion (ISE), the integral of quadratic error multiplied by time (ITSE), the integral absolute error (IAE), and the integral of absolute error multiplied by time (ITAE).It is observed that the results of the quantitative evaluation of the strategy proposed in this work have better results in the performance indexes ITSE and ITAE for the control action u(t) and in the performance indexes ISE, ITSE, and ITAE for the output signal y(t) (nominal values) with respect to the work presented in [22].For simulations where uncertainties in the plant are considered, the evaluation results show that the proposed observer-predictor scheme performance indexes ITSE, IAE, and ITAE are better in the control action u(t); likewise, for the output signal y(t) it is observed that in all the performance indexes, the proposed strategy presents minor error magnitude with respect to [22].

Example 2
Consider the unstable process with delay given by For this example, b 1 = 1, a = 0.5, b = 1, and τ = 1.4 are the parameters of the process.Thus, by substituting and in condition of Theorem 1, we obtain τ < 3.236.This means that the maximal delay-time size allowed by the proposed control strategy is τ < 3.236, and since the delay time of the process (τ = 1.4) is less than this bound, therefore the proposed control strategy can stabilize this process.Then, Corollary 1 can be used to obtain the values of g 1 and g 2 .For the simulation, g 1 = 0.5 and g 2 = 0.03 are used.Further, for the simulation experiment, the parameters of the control C(s) given by ( 22) are chosen as k p = 27.29,k i = 9.10, and k d = 18.19.An initial condition of 0.1 is used for each of the plant states.A step input reference, R(s), with an amplitude 5, is set.Figure 8 shows the behavior of the estimation error, and we see that e ω (t) = 0 in the steady state.To test the robustness with respect to uncertainty in the delay time, the analysis presented in Section 6 is used.The interval of time delays that guarantees closed-loop stability is given by 1.33 < τ 0 < 1.467 where the nominal value of the delay is τ = 1.4.This is shown in Figure 9. Figure 10 shows the response of the system to uncertainties in the model parameters, and it can be seen that the system is stable with a parameter a = 0.5 greater than +34% of its nominal value of the unstable pole.Figure 11 shows the response in which the tracking of the steady state reference can be observed.Likewise, it shows that the system is robust to some uncertainty with respect to the nominal delay (τ = 1.4).Finally, Figure 12 shows the output response when white noise effects are considered at the output.In this figure, it can be observed that system stability prevails in the face of such effects.
On the other hand, for this example, b 1 = 0.5, b 2 = 2, a = 0.2, b = 0.2, and τ = 0.5 are the parameters of the process.Considering the observer-predictor presented in this paper, the stability condition given in Theorem 3 is τ < 8.3, which is satisfied.Using Corollary 1, we calculate the values of g 1 and g 2 .In this example, the values of g 1 = 10 and g 2 = −5 are used.For the control C(s) given by Equation ( 22), the parameters used are k p = 3.94, k i = 0.32, and k d = 6.57.
Additionally, in order to achieve a better response, a two-degree of freedom controller of PID type [31] defined by .57, k i = 0.32, α = 0.8, and β = 0.8.For the simulation, a unitary step input reference R(s) = 1/s and a perturbation of the actuators at 40 s is used.The output response of the control scheme proposed in [19] and the output response of the observer-predictor scheme presented in this work are shown in Figure 13 to noise effects.It can be seen that in the presence of a perturbation acting at 40 s, both strategies allow perturbation rejection and reference tracking.It is important to highlight that in [19], necessary and sufficient conditions are not provided to stabilize systems with time delay, while in the observer-predictor presented in this work, it is possible to control the system of example 2 with a delay τ < 8.3.Figure 14 shows the output response under an uncertainty in the delay τ.It can be seen that with a delay uncertainty of +100%, both strategies still have a stable output response.However, at a delay uncertainty τ of +120%, it can be seen in Figure 15 that the strategy presented in [19] is no longer stable, while the proposed observer-predictor scheme preserves a stable response.Figure 16 shows the control signal u(t) response under an uncertainty in the delay τ.
Finally, Table 2 shows the evaluation results for different performance indices.It can be observed that the observer-predictor scheme presented in this work has better performance indexes in the ITSE, IAE, and ITAE for the control action u(t) and ISE, ITSE, IAE, and ITAE indexes for the output signal y(t) (nominal case) with respect to the work reported in [19].When a delay-time uncertainty is considered in the system, the proposed observer-predictor scheme has better performance indexes ITSE, IAE, and ITAE in the control action u(t) and the indexes ISE, ITSE, IAE, and ITAE for the output signal y(t) with respect to the work presented in [19].

Practical Contributions
In this section, we highlight some of the major practical contributions of the proposed control strategy developed in this work.

1.
When theoretical results (such as those provided in the present manuscript) are taken in order to implement real-practical experiments, it is very important to know quickly if the theoretical proposed control strategy can be used for solving the stability problem or not.In this way, providing necessary and sufficient conditions for the existence of the stabilizing control strategy allows one to obtain this information immediately by making simple computations of the stated conditions.This facilitates the task of the control design for the engineers, in contrast with control strategies that are developed under heuristic methods or just present necessary or sufficient conditions.2.
In general, the time-delays phenomenon due to measurement of variables, material transportation, or teleoperation complicates the design of control strategies, however in the particular case when the size of time-delay is greater than the dominant dynamic of the system, the control design has an additional degree of complexity.This is because, under this situation, the time delay has greater adverse effects than the behavior of the controlled system when the delay time is small.As an example of this problem, there are many proposed control strategies [11,24] that are well performed for a "small" time delay, but these control strategies cannot be used for controlling systems with large delay terms.This problem is illustrated and partially solved in the new version of the work; see Example 3 (Figure 15).Therefore, in this work, it is possible to control systems with greater delay size with respect to previously reported works ( [11,24]).For example, in [17], the maximum size of the 1 for a second order delayed system, in the present work, this bound is τ < 2 1/a − 1/b 1 + 1/a 2 − 1/b 2 1 , i.e., the double of the previous delay size.This practical advantage allows the controlling of a greater variety of systems, ensuring closed-loop stability.

3.
In practice, a control strategy that provides stability with respect to the uncertainties of the model is desired.These problems arise mainly from model mismatching.In this way, analytical robustness with respect to time delay is developed in Section 6.Additionally, the proposed control strategy has been evaluated via simulations on face-to-model mismatching of the process, obtaining positive results such as the preservation of closed-loop stability under the adverse mentioned condition.4.
Many industrial solutions (such as chemical engineering applications) require the regulation problem, which consists in keeping a variable at a desired value even in the presence of disturbances.In this way, the capability of a control system for tracking desired step-references and disturbance rejection is essential for solving the regulation problem.Thus, the results provided in Section 5 are related to this issue and illustrated in numerical simulations.

5.
The control action and the output measurable variable performances are qualitative and quantitatively evaluated using a numerical simulation; the control action can be seen as the spent energy of the control stage to achieve the stabilization of the process, which is a very important issue from a practical viewpoint.

Conclusions
Dynamic systems with relatively large delays with respect to the dynamics of the system are tackled in this work.A control proposal based on a delay predictor (observer) scheme is presented.A fundamental part of this work focuses on achieving the convergence of the observer-predictor scheme.The approaches used in the study allow the obtaining of necessary and sufficient conditions for the convergence of the scheme presented.Using this scheme, it is possible to control systems with longer delays than those treated in previous works [11,24].
The general idea of the proof is presented below.Later, a formal analysis using a root locus approach is presented.
Let us take n = 5, for simplicity and without loss of generality, in the analysis.Then, in the following steps, consider the root locus diagram generated by the characteristic Equation (A2).
Step 1. Assume g 1 > 0 and g 2 = 0.When g 1 varies from zero towards positive values, n − 1 poles located at the origin move towards infinity: one pole moves in the direction of a breaking point located at the real axis, and the pole located at z = e aT also moves from right to left on the real axis in the direction of to the breaking point, and for some values of g 1 > 0, the trajectory enters the stability region (unit circle).These trajectories are depicted in Figure A3.In Figure A3, the gray points in the trajectories indicate the location of the roots for a specific value of g 1 > 0, denoted by ḡ1 .Notice that the pole on the right real axis (z = e −b 1 T ) does not move since g 1 has no influence on this pole.
Step 2. With g 1 = ḡ1 and g 2 > 0 variable.This step is explained by using Figure A4 when g 2 varies from zero towards positive values, the pole located at z = e −b 1 T is involved in the locus of the roots, which moves from left to the right on the real axis, causing a breaking point with the relocated pole (originally the pole outside the unit circle) in step one.Note that under these circumstances, a breaking point inside the unit circle is obtained.Gray points in the trajectories in Figure A4 indicate the location of the roots for a specific value of g 2 > 0, denoted by ḡ2 .
The facts depicted in Steps 1 and 2 suggest that it is possible to obtain g 1 = ḡ1 and g 2 = ḡ2 such that all the roots of the locus are inside the unit circle and therefore the discrete injection scheme shown in Figure A2 is stable.

Figure A4.
LGR with g 1 > 0, g 2 > 0. Now, considering the previous developments, the breakpoint derived from the aforementioned analysis is calculated.Consider Step 2, where a breaking point is located in the real axis between the pole located at z = e −b 1 T and the relocated pole that was moved into the unit circle in Step 1.Then, to achieve closed-loop stability, the breaking point should be inside the unit circle.Then, the location of the breaking point is calculated.This can be done by solving for g 2 from Equation (A1), and then its derivative is The breakpoint of interest is obtained by calculating the roots of Equation (A3).Indeed, the breaking point of interest is one of the roots of the polynomial (A3).These roots are analyzed by considering that the parameter g 1 is variable.This allows studying the behavior of the roots of (A3) through a root locus diagram.Note that this root locus diagram does not depict the movement of the roots of the closed-loop injection scheme but the movement of the breaking points related to the system.Then, it is possible to associate a fictitious closed-loop system with Equation (A3).Let us consider an open loop system given by Then, in this new fictitious system (A5), it is possible to analyze the behavior of the breaking points.Taking into consideration an analysis of the locus of the roots for the system (A4) with g 1 as a variable parameter, Figure A5 shows that there exists a breaking point between the pole located within the unit circle and the pole located outside.
The location of the breaking point is obtained by solving for g 1 from the characteristic polynomial Equation (A5) and its corresponding derivative, provided by Performing the same analysis as that used to obtain g 1 , the value of g 2 is given by for some 2 > 0.
The values of g 1 and g 2 ((A14) and (A18)) stabilize the injection scheme shown in Figure 1.

Figure 5 .
Figure 5. System output response with perturbation of example 1.

Figure 6 .
Figure 6.System output response with uncertainty in the plant of example 1.

Figure 7 .
Figure 7. System output response with uncertainty in the plant of example 1.

Figure 8 .
Figure 8. Convergence of the estimated signal from example 2.

Figure 10 .
Figure 10.Output response to plant uncertainty in example 2.

Figure 12 .
Figure 12.Output response to noise effects of example 2.

Figure 13 .
Figure 13.Output response to noise effects of example 3.

Figure 14 .Figure 15 .
Figure 14.Stable system output response with time delay uncertainty τ of example 3.

Table 1 .
Comparative table of quantitative evaluation of the two strategies of example 1.

Table 2 .
Comparative table of quantitative evaluation of the two strategies of example 3.