The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold

Zero dynamics have crucial effect on system analysis and controller design. In the control analysis process, system performance is influenced by the unstable zero dynamics, greatly. This study concerns with the properties of limiting zero dynamics when the signal of controlled continuous-time systems was reconstructed by forward triangle sample-and-hold (FTSH). FTSH is a newly sample and hold method in the signal reconstruction area. However, more theoretical details about the limiting zero dynamics of the resulting discrete-time systems need to be revealed. Firstly, the framework for the limiting zero dynamics in the situation of sufficiently small or large sample period is introduced. Furthermore, this study provides the stable conditions of limiting zeros in the two different sampling situations. The results indicate that one can select a suitable variable parameter value of FTSH to replace the sampling zeros of discrete-time system locate inside the stable region. This paper also reveals the truth that FTSH has the outstanding advantage compared with the BTSH using theoretical analysis method. Finally, example simulations verify the effectiveness of the results in this study.


Introduction
In practical control engineering area, the continuous-time system can not be used to design and realize the control strategy, directly. We should deduce a corresponding discrete-time system to represent the original continuous-time system by using a sample and hold device. The zero properties of the continuous-time linear system affect the behavior of the controlled system. Moreover, for discrete-time system, the minimum phase (MP) can be reflected by the location of sampling zero, and it plays an important role in system analysis and control strategy design process [1,2]. In the controller design process of system, the appearance of unstable zero can be seen as a great unfriendly obstacle to construct control strategy, such as model reference adaptive controller and model matching system [3]. Therefore, it is a meaningful work to avoid the unstable zero arise from the discretization process. A considerable amount of researchers have poured many efforts to study this question in recent decades [4][5][6][7][8][9][10][11][12][13][14]. Discrete-time model of the continuous-time system relay on the reconstruction method about the input and output signals. During the process of the discretization, the relationship of the pole has a single mapping: p i $ e p i T , p i is the continuous-time part and e p i T is the discretetime part, T represents sampling period. However, this simple relation cannot be preserved for the zeros, more complicated transformations about zeros in the discretization process need to be concerned. The dynamics of discrete-time system is complex, only for the classification about the zeros have two categories: intrinsic zeros and limiting (or sampling) zeros, respectively [4,15]. The properties of limiting zeros mainly depend on the method of discretization, the scale of sampling rate and the relative degree of the original system. Limiting zeros have no continuous-time counterpart and simple map relation in the z-domain. Many previous literature have revealed that the limiting zeros converge, as the sampling rate goes to infinity, to the roots of special polynomials (such as Euler-Fröbenious polynomials [4,6,16] and its' modified form [17]). Because of the truth that the limiting zeros is the function about sampling period T, only select some special case of sampling period one can deduce the expression of limiting zeros. For example, the sufficiently small or large sampling period attracts a lot of researches in [6,10,[17][18][19][20]. Åström et al. provide the groundbreaking work to research the exiting question in the limit case [4]. Hagiwara et al. focus on the properties of limiting zeros when using zero order hold (ZOH) or first order hold (FOH) cascade into the control framework with sufficiently large or small sampling rate [6]. In [18], authors extended this work into fractional order hold (FROH). The stability of limiting zeros of discrete-time system in the case of backward triangle sample and hold (BTSH) for sufficiently small or large sampling period was discussed in [9]. These previous literature has revealed the truth that to research the properties of zeros with limit sampling rate value is an important work to explore more potential value about one specific sample and hold method, which will provide a theoretical reference in the practical engineering applications.
Two new sample-and-hold functions: forward triangle sample and hold (FTSH) and BTSH were investigated in [21,22], recently. Two triangular sample and hold functions were considered as an alternative to the ZOH and shown that they can replace the sampling zeros of the corresponding sampled-data systems locate in unite circle by using BTSH or FTSH while ZOH can not to do. However, the mainly results in [21] lay emphasis on numerical investigation. Therefore, it is natural that one should analytical investigate the properties of zeros of a discrete-time system with the above two sample and hold functions. Consequently, Ou et al. [9,23] were analytically investigated the properties of zeros and zero dynamics of sampleddata systems with BTSH for linear and nonlinear systems, respectively. They were also researched the sampling zero properties of discrete-time system with FTSH in [22], but only the case of sufficiently small sampling period was studied. Moreover, points to the fact that the properties of limiting zeros should be considered completely in all the case of sufficiently limiting sampling periods [4,6]. It is a significant research work to explore the capabilities of sampling zero and advantages of FTSH with limiting sampling period.
The main purpose of this study is systematically to reveal the properties of limiting zeros of discrete-time system with FTSH in details. The present study clarifies the corresponding results of limiting zeros for sufficiently small and large sampling period when the FTSH is applied to discretize the continuous-time control system. Further, we also present the stable conditions of the limiting zeros.
The contributions of current study are marked as shown in the following three aspects: (I) FTSH was used as the signal reconstruction, authors deduced a novel polynomial. The stable conditions of limiting zeros were provided, which is a vital complement to the results in [22].
(II) This paper analytically reveals the truth that FTSH has the outstanding advantage compared with the BTSH, which is coinciding with numerical research in [21].
(III) The theoretical framework of properties of limiting zeros for some eligible systems with FTSH was provided. Authors verified that by using the results in this study can replace the limiting zero in the stable region.

Preliminaries
As mentioned in the introduction, mainly purpose of this paper is to reveal the sampling zeros, for FTSH discretization, under limiting value of the sampling period (i.e. T ! 0 and T ! 1). For purpose of better understand the corresponding results of this paper, some wellknown preliminaries knowledge are reviewed in this section. The transfer function of a continue-time system is defined as G c (s). G F (z) represents the corresponding discrete-time transfer function of G c (s) under FTSH as the signal reconstruction, where s denotes a complex variable for Laplace transform and z = e Ts .
the transfer function (1) is irreducible. The integer r = n − m is the relative degree (i.e. pole excess). Then, one can obtain the state space expression about the controllable, observable and single-input single-output linear continuous-time system.
where A c 2 R n×n , b 2 R n×1 and c 2 R 1×n represent n-th order state matrix, column and row vectors, respectively. Moreover, u(t) 2 R, x(t) 2 R n×1 and y(t) 2 R represent the input, state vectors and output of the controlled system, respectively. The following formula provide a bridge between transfer function and state space expression.
Based on the different expression of linear system (1) and (2), those system zeros (also can be defined as transmission zeros or invariant zeros, those three names are coincide) can be deduced by computing the part of numerator polynomial in (1).
In this study, the discrete-time system model is composed of a hold circuit, the continuoustime system and a sampler in cascade, where the signal reconstruction method diagram as shown in Fig 1, and the mathematical formula of FTSH is where k 2 N, N denotes the set of natural, f 2 (0, 1] represents the duty cycle of the switched input, and u Z (t) denotes zero-order hold (ZOH) value in each sampling interval [kT, kT + T). Thus, one can obtain the relation.
Then, the state space of corresponding discrete-time system with FTSH is Analogously, we can use the following equation to get the discrete-time domain transfer function G F (z).
the denominator polynomials D F (z) and numerator polynomial N F (z) can be deduced by the following equation.
Obviously, zeros and poles of discrete-time system G F (z) could be obtained from (8) by calculating the denominator and numerator polynomials.

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold In this section, we have reviewed some preliminaries knowledge and previous results of two types systems with FTSH as sample and hold function. The zeros of the discrete-time system, for the case of sufficiently small or large sampling period, are often called limiting zeros [4]. Based on the preliminaries knowledge, in the following section, some more details about the limiting zeros of discrete-time system using FTSH as method will be provided, firstly.

Limiting zeros of the discrete-time system with FTSH
The purpose of this part is mainly to reveal the properties of limiting zeros about discrete-time system in the case of FTSH. Before showing the corresponding results of limiting zeros of discrete-time system, which are generated by using FTSH as the sample and hold function to construct the input signal. An interesting polynomial in (11) need to be introduced.
Definition 1 [22]: Suppose that f 2 (0, 1], r 2 N ⩾ 0, exist a number p such that f + p = 1. A novel polynomial has the following properties : Remark 1: Considering the mainly results of this paper, the polynomial in Definition 1 plays a crucial role in revealing the properties of limiting zero of discrete-time system with FTSH. More relevant results about the novel polynomial will be presented in the following part of this section.

Remark 2:
More details about the method to obtain the result in Definition 1 can refer the process of Theorem 1 and 2. Moreover, Definition 1 reveals that the interrelationship between limiting zeros and relative degree when using the novel sample-and-hold method FTSH.
In order to provide more theoretical reference knowledge about Definition 1 for the potential reader of this paper research topic, several properties of the novel polynomials (11) will be established, firstly. The corresponding results mirror the extensive literature available on the properties of the standard Euler-Fröbenius polynomials-one can see [16,24,25].

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold Proof: (1): Based on the results in Definition 1, we can rewrite the expression of Z F,r using the following form Hence, the Schur determinant Lemma [26] used to compute the above determinant, one can obtain

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold where adjZ F,r−1 is the adjoint-matrix of Z F,r−1 , i.e., from the theory of matrices, Calculating the determinant and subsequent sub-determinants along with the first column j − 1 times, the following results about the cofactor matrix obtained Ground on the above properties, using the recursive relation of them, one can get the following results Using the relation (11) in Definition 1 can complete the proof.
(2): Combined with the Definition 1, it is not hard to find that proving the result of this part is equivalent to proving the determinant det D F,r = T r det Z F,r . Here, we use the induction method to prove this part.
The detail steps of induction are as follows: First, we note that det D F,1 = T 1 det Z F,1 is true. Then, assuming the result of determinant det D F,i = T i det Z F,i holds for i = 1, � � �, r.
Finally, if we want to prove the result, only need to proof det D F,r+1 = T r+1 det Z F,r+1 is also true.
Considering the proof process as shown in the above part (1), we can obtain Using the same compute method in the part (1), it is not difficult to get the following equation Thus, using the above equation replace the corresponding part in (15) and the determinant of D F,r+1 as shown in the following

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold From the conditions of induction hypothesis, we have assumed that detD F,i = T i detZ F,i , i = 1, � � �, r. is true, and then, we have where the last equality can be deduced from the relation of (14). Until now, the second part of Theorem 1 was completely proofed. Thereafter, the properties of limiting zero dynamics of discrete-time system with FTSH in the case of limit sampling period (i.e. T ! 0 or T ! 1, respectively.) will be discussed in the remaining part of this section. Firstly, we reveal the limit expression of sampling zero about sampling period when it is sufficiently small.
Theorem 2 [22]: When a FTSH as the signal reconstruction to construct the input value of a linear system, one can obtain the following two results.
(1) If the controlled system is composed of a r-th order integrator with continuous-time transfer function G c ðsÞ ¼ 1=s r . Obviously, the relative degree is r. We can get the corresponding discrete-time state space expression: where ; and the output of the system is y k = x 1,k . Then, the corresponding transfer function G F (z) of

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold the exact sampled-data system is where B F,r (z, f) is the novel polynomials.
(2) If the continuous-time linear system with transfer function G c (s), which is a strictly proper n-th order. Then, the limit expression when T ! 0 of the corresponding exact sampled-data model can be shown as follows: where r = n − m is the relative degree.

Proof:
Because of that the proof details of case (1) in Theorem 2 has been provided in [22], for the convenience of potential reader to better understand this Theorem, here, a detailed analysis was carried out for the case (2) in Theorem 2.
the proof process of case (2): Following from the method in [4], one can obtain the ideas to proof the result (20). Specifically, suppose that the continuous-time system G c (s) is a rational strictly proper (i.e. strictly physical realizable) transfer function expressed as G c ðsÞ ¼ Kðs À z 1 Þðs À z 2 Þ � � � ðs À z m Þ ðs À p 1 Þðs À p 2 Þ � � � ðs À p n Þ ; K 6 ¼ 0 with r = n − m is the relative degree, where z i , i = 1, 2, � � �, m. and p j , j = 1, 2, � � �n. represent zeros and poles, respectively. Then, the exact discrete-time system model can be obtained by using inverse Laplace transform and Z-transform, Therefore, as T ! 0, one can characterize the asymptotic model above the equation as follows

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold the integration path of (21) has an infinitesimal detour around origin. Based on the results in part (1), the discrete-time system of an r-th order integrator system has the following transfer function From the above deduce process, we also have the following equation Using s = ω/T as the method of variable substitution again, we can get Finally, combining (21), (22) and (23), one can use the following limit expression to characterize the asymptotic model lim T!0 T À r G F;r ðzÞ ¼ K � B F;r ðz; f Þ ðn À mÞ! � ðz À 1Þ r Remark 3: In Theorem 2, we notice that for the general strictly proper linear system (1) when the signal is reconstructed by FTSH. It is easy to know that n poles converge to 1 as e p i T of the corresponding discrete-time system when the sampling period T ! 0, and m zeros close to 1, moreover, other zeros are equivalent to the roots of B F,r (z, f). Obviously, the sampling zeros are influenced by the values of relative degree and the duty cycle value f of FTSH.
Next, a result of limiting zeros for a sufficiently large sampling period in the case of FTSH will be deduced in this part of this section. Åström et al. [4] have shown that all the zeros of discrete-time system with ZOH approach z = 0 as T ! 1 when the original continuous-time system G c (s) is stable and satisfy the initial condition G c (0) 6 ¼ 0. However, the situation of system exists unstable part is not considered. Hagiwara et al. [6] have continued in-depth to study the case when original continuous-time controlled system has unstable poles and have shown the properties of limiting zeros of sampled-data system with ZOH and FOH, analytically. The corresponding research about limiting zeros in the case of FROH and BTSH has been explored in [9,18], respectively. In next part of this section, we will explore the capabilities of limiting zero dynamics of discrete-time system with FTSH under large sampling period which will greatly enrich the theoretical framework of discrete-time system with FTSH.
Theorem 3: Assuming that G c (s) in (1) is not having any pole on the imaginary axis. n s and n a (= n − n s ) are the number of stable and unstable poles, respectively. Thereafter, suppose that the continuous-time system G c (s) can be decomposed into the following expression: Then, based on the decomposed realization form and (10), one can obtain (26) (as shown in the top pf this page) Using the same steps as shown in [6], n a À 1 x n a À 1 ðax þ 1Þ n s G a 0 ð Þ as T ! 1. Thus, if G a (0) 6 ¼ 0, the sampled-data system has n a −1 roots diverge to x = 0, n s roots approach to x ¼ À 1 a . Based on the relation of them in (27), one can obtained the results in case (ii). And the proof is complete.

Criteria for the limiting zero stabilization
In the above section, we have provided the limiting expression of sampling zeros of discretetime system with limit sampling period. But we only gave the expression about the sampling zeros. Therefore, to reveal the stability of sampling zeros is an important work. Thereafter, the stable conditions of limiting zeros of discrete-time system with FTSH need to be researched and provided precisely. Firstly, the stability condition of limiting zeros of sampled-data system in the case of sufficiently small sampling period is shown using the following theorem.
Theorem 4 [22]: Assuming that the original continuous-time system G c (s) has no zeros on the imaginary axis. And its corresponding discrete-time system is generated by using FTSH in the case of sufficiently small sampling period. All the limiting zeros are stable when r = 2 if all the zeros of G c (s) are stable and the selectable variable parameter f satisfies 0 < f ⩽ 1.
Remark 4: Generally speaking, the relative degree of many linear or nonlinear mechanical systems in practical field is two [27,28]. As the results in Theorem 4 shown, one can find the feasible solution of the variable parameter of FTSH to ensure the corresponding discrete-time limiting zeros are stable.
Remark 5: An insightful observation in the Theorem 4 is that the stability properties and conditions are not given when the relative degree of continuous-time system is greater than two. Using Jury's test can get the stability conditions of the limiting zeros by computing the roots of corresponding numerator polynomial B F,r (z, f) = 0. Thus, we notice that it cannot find a proper f to let limiting zeros stable when the relative degree of continuous-time system is three.
Next, the stable condition of limiting zeros of discrete-time system in the case of sufficiently large sampling period is also provided using the following theorem.
Theorem 5: Suppose that original continuous-time system G c (s) has no poles on the imaginary axis, and let the number of its unstable poles is n a . Further assume that G c (s) is decomposed into (24). we have Case(i) n a = 0: The system G c (s) is stable, when the sampling period T is sufficiently large, if want to make the zeros of G F (z) are stable, the controlled system should satisfy Case(ii) n a = 1: When using large sampling period T, all the discrete-time zeros of G F (z) are stable and converge to zero if and only if G a ð0Þ 6 ¼ 0 Case(iii) n a ⩾ 2: The corresponding sampled discrete-time system G F (z) always has an unstable limiting zeros when the sampling period is sufficiently large.
Proof: Based on the previous analysis and proof process in Theorem 3, the corresponding results in this Theorem can be deduced easily. Thus, the proof details of this Theorem are omitted.
Remark 6: Comparing the results of this paper and the corresponding results with BTSH in [9], when the sampling period is sufficiently large, the FTSH has better performance for limiting zeros than BTSH. For example, the resulting discrete-time system has minimum phase (MP) characteristics, using FTSH as signal reconstruction method, when the original continuous-time systems have an unstable pole. However, the BTSH cannot provide this property.
Remark 7: Compared with the results of limiting zeros between FTSH and BTSH for sufficiently small sampling period. The feasible region (or interval) of the FTSH is larger than BTSH.
In the previous research, it is shown that FTSH has the outstanding advantage compared with the BTSH with selecting the suitable sample-and-hold parameters based on the simulation studies [21]. And in the mainly part of this paper, it has analytically revealed that FTSH can provide a better performance of limiting zeros and a larger feasible region than BTSH. Furthermore, the advantages of FTSH will be validated by some numerical simulations in the next section.

Numerical simulations
From the above analyses we can know that the properties of limiting zeros for discrete-time system model with FTSH are connected with the relative degree, sampling period and the number of unstable poles. Thereby, we provide some interesting numerical examples to verify the results in the main part of this study.
Firstly, the continuous-time system with relative degree two and its' continuous-time transfer function is Obviously, the original controlled continuous-time system is stable and minimum phase. Based on the corresponding results in Theorem 2 and Theorem 4, we can select the appropriate parameter f of FTSH to replace the limiting zero and let it is stable. Therefore, we choose the parameter of FTSH and the location of limiting zeros of (31) simulation result is shown in Fig 2 (This numerical can also refer [22]. But in our previous research work, system with unstable part is not mentioned, which will be analysized in the next part.). As shown in simulation diagram in Fig 2, the discrete-time zeros can be placed inside the unite circle no matter that the sampling period is small or large. Moreover, the phenomenon in simulation are coinciding with the results in Theorem 2 and 3.
Next, we need to consider the situation that the original controlled system has unstable real poles. Combining with the results in Theorem 3 and Theorem 5, two examples with different number unstable real poles are given to verify the results. First of all, one can select a Finally, in order to confirm the result that FTSH has the outstanding advantage compared with the BTSH, which need to illustrate by a numerical simulation, we also introduce the Disk drives serves (see Fig 5) [9, 29] as a didactic practical example to verify the advantages.
We select parameters' value of Disk drives are same as in [9,29]. Then, the transfer function of Disk read-write head model is the relative degree two and has the following form Combined with the feasible interval of sample and hold parameter, we select the device parameter of BTSH and FTSH with f = 0.55 and T 2 (0, 0.2]. Then, the limiting zeros of the corresponding sampled-data system with two sample and hold devices (i.e. FTSH and BTSH) are presented in Fig 6. Evidently, the limiting zeros always can be replaced in the unite circle region by selecting the value of BTSH and FTSH. And interestingly, the value of limiting zeros in the case of BTSH is larger than FTSH (i.e. the location of sampling zeros in the case of FTSH are closer to origin than BTSH).

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold Furthermore, we also select the traditional ZOH as the alternative signal reconstruction method to generate the inputs of the Disk drivers system (34). The absolute values of limiting zeros of the system as shown in Fig 7. Obviously, we can not to let the corresponding discrete-time system is MP by selecting a proper sampling period when use the traditional ZOH as the signal reconstruction method (i.e. the discrete-time system always has a zero locate outside the unite circle.).
In Fig 6, we only have shown the location of sampling zeros with respect to sampling period T under a fixed duty cycle f. However, it is not clear to know the location of sampling zeros with different f. Then, one can select a fixed sampling period T = 0.1, and show the comparative simulations of BTSH and FTSH about the location of sampling zeros for system (34) with different duty cycle f in Fig 8. Based on the previous research in [9], the feasible interval of f is (0, 0.75] in the case of BTSH. As shown in Fig 8, the simulation result of BTSH is coincide with [9]. Another side, the sampling zero locates on the boundary of unite circle when f = 0.75. Both in the case of BTSH and FTSH, the simulation results and the theoretical results are coincide.

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold

Conclusion
This paper has dealt with the properties of limiting zero dynamics of sampled-data system with FTSH as the sample and hold device. The expression of limiting zero polynomial of continuous-time linear system in the case of FTSH has been given. And the stable condition of sampling zeros of the corresponding system has also been provided for the case of sufficiently small sampling period, where FTSH can replace the sampling zeros of the sampled-data system into the stable region while ZOH fails to do so. Moreover, this paper has also derived the properties of sampling zeros of sampled-data system with FTSH for sufficiently large sampling period, FTSH can let all zeros of the discrete-time systems are stable when the original continuous-time system has one unstable poles while BTSH can not provide this properties. Compared with the previous results, this paper was analytical extending work from the simulation research in [21]. It has also analytically revealed that FTSH can provide a better performance of limiting zeros and a larger feasible region than BTSH. We have verified the results in this

PLOS ONE
The limiting zero dynamics of discrete-time system based on forward triangle sample-and-hold paper are correct through numerical simulation. Further work will concentrate on the characteristics of zero dynamics for nonlinear system with FTSH.