Abstract
Double integral plants under relay feedback are studied. Complete results on the uniqueness of solutions, existence, and stability of the limit cycles are established using the point transformation method. Analytical expressions are also given for determining the amplitude and period of a limit cycle from the plant parameters.
1. Introduction
Relay feedback forms one important class of nonlinear systems and can cause complex nonlinear behaviors. Its early analysis can be traced to 1950–1960s, and afterwards, two basic approaches, the time domain approach [1, 2] and the frequency domain approach [3], emerged. Although they are almost identical, the frequency approach is more popular because of its ease of manipulation. However, as a general method for relay analysis, such a frequency approach also has some limits in itself. Firstly, only necessary conditions on the existence and stability of possible limit cycles are obtained. Secondly, the conditions are usually expressed as a summation of infinite items. It is certainly desirable to find both sufficient and necessary conditions of the existence and stability of limit cycles for relay feedback systems, as well as to give such conditions explicitly in terms of system parameters without any requirement on numerical computation.
In our recently published paper [4], we presented a complete relay analysis for a class of servo plants, , (), including the uniqueness of solutions, existence, and stability of limit cycles, and its amplitude and period, using a similar analysis to [5]. Naturally, it is desirable to investigate whether such an analysis can be extended to double integral plants, that is, the case of . Most people believe that double integral plants have no limit cycle under relay feedback [1, 3], but our analysis shows that this is not true. Actually, a limit cycle for double integral case is observable, which completely depends on the system initial conditions. Moreover, we also find that the analysis for the case of in [4] is not applicable to the case of . This is because for , the derivative of system output is stable and bounded as . While for , and the derivative of system output becomes unstable and not bounded anymore. Hence, the extension is not straightforward and the results for double integral plants are quite different from our previous ones.
This paper is organized as follows. Section 2 gives the main result followed by the proofs and a few remarks, where more details and explanations are given. Conclusions are drawn in Section 3.
2. The Result
Consider a double integral plant whose state-space representation in the controllable canonical form is given by where , are the state, output, and input of the system, respectively. The plant is under the relay feedback control: where and are the relay amplitudes, , and are the relay hysteresis with . We assume since otherwise (4) becomes a constant but no longer a relay control. The relay control is depicted in Figure 1. The initial function for is where is the initial time and . We call (1)–(5) a relay feedback system (RFS) which is depicted in Figure 2, where represents the controller.
If the RFS generates a limit cycle, let and be the half period and the extreme value corresponding to , respectively, and and be the half period and the extreme value corresponding to , respectively, as shown in Figure 3. (resp., ) is the rising (resp., decreasing) period between the switching time, (resp., ), when the limit cycle trajectory traverses (resp., ) and the peak (resp., valley) time.
We are now in a position to state the sufficient and necessary conditions for the existence of solutions, the existence and stability of limit cycles, and the amplitudes and periods of limit cycles.
Theorem 1. Consider the RFS for double integral plant in (1). (i)A unique solution exists for any initial condition.(ii)A limit cycle exists if and only if , and any of the following holds: (a);(b) but . If this is the case, the limit cycle is formed after the first switch, which is unique with two switchings per period.(iii)If a limit cycle exists, it is globally marginally stable. (iv)If a limit cycle exists, its amplitude and period are described by
Proof. For a RFS, the input is a piecewise constant function. Without loss of generality, suppose
as shown in Figure 4, where
Define the switching planes:
If the trajectory of traverses (resp., ), that is, (resp., ) at some instant with (resp., ) and (resp., ), then the instant is called a switching time. In particular, denotes the switching time when the th switching takes place.
The state response of (2) to in (7) is given by
where , . It is easy to show that both and are continuous at switching time , as for . However, , while indicating is not continuous at due to by our assumption. Straightforwardly, this discontinuity persists for the general case of .
If considering the time between two consecutive switchings only, that is, , (10) is simplified as
where , and
Since and , it yields
It can be seen from (13) that , and will affect behavior of and thus determine whether or not will reach the switching level of or . Without loss of generality, suppose the initial condition , thus . To see how evolves with time, the following four cases are considered which are mutually exclusive and cover all possible cases.
Case 1. and . For , , which implies so that never traverses for all . The trajectory of is governed by (11) with and .Case 2. and . For , it follows from (14) that , which implies there always exists some such that and , that is, traverses at . For , it follows from (14) that if , which implies so that never traverses for all . The trajectory of is governed by (13) and (14) with and . If , it follows from (14) that , which implies there always exists some with and so that traverses at . Afterwards, with the same analysis as above, it is straightforward to verify that will traverse and alternatively and consecutively.Case 3. and . For , it follows from (14) that , which implies there always exists such that and , that is, will traverse at . From onwards, following the same analysis in Case 2, if , never traverses for all , the trajectory of is governed by (13) and (14) with and . Otherwise, if , will traverse and alternatively and consecutively.Case 4. and . For , it follows from (14) that . If , there always exists such that for all . If , then so that never traverses for . The trajectory of is governed by (13) and (14) with and . On the contrary, if , there always exists satisfying such that and , that is, traverses at . If , then , which also implies traverses at some . For , it follows from (14) that if , which implies so that never traverses for all . The trajectory of is governed by (13) and (14) with and . If , it follows from (14) that , which implies there always exists some such that and , that is, traverses at . For , it follows from (14) that , which implies never traverses again for all . The trajectory of is governed by (14) with and .
In view of the above analysis, the solution to (2) always exists. This completes (i).
We can also see from the analysis in (i) that consecutive switchings between and take place if and only if . Under this condition, let be the switching instant of on the plane , that is, and , for sufficiently close . This implies . The case can occur only when and then is a fixed point of . But a fixed point is unique. Hence, no limit cycle exists in this case.Assume , then switches from to at . In Case 2, we have shown that decreases and becomes negative until traverses at . Consider the Poincare mapping
The function is defined on the ray and takes values in the ray . Let and . By definition, and . Since both and are continuous at , it follows from (13) and (14), respectively, by taking that
It follows from (17) that
Substituting (18) into (16) yields
which implies
If , “>” holds in (20) so that the Poincare mapping (15) is not a contraction and no limit cycles exist. If , “” holds in (20), which implies a limit cycle is formed after the first switching. As , must be satisfied so that (a) and (b) can be easily derived with the help of (11). This completes (ii).
Once a limit cycle is formed, if is changed due to some disturbance, cannot remain at the original limit cycle, but form a new one with different amplitude and period, which means the limit cycle is marginally stable. Since the initial condition is arbitrarily chosen, the limit cycle is globally marginally stable. This completes (iii).
It has been shown in (ii) that for , a limit cycle always exists after the first switching. For , it follows from (14) and with the help of (20) that
Thus,
Let , then . Substituting it into (13) yields
In the similar way, we have
For , it follows from (13) that
which gives
where negative root is ignored since . Substituting (27) into (14) yields
Substituting (28) into (22)–(24), respectively, yields
which completes (iv).
Remark 1. In our previous relay analysis for a class of servo plant as , (), we proved the existence of a limit cycle by showing the Poincare map is a strict contraction, that is, for all . Here for the case of , differentiating both sides of (16) yields Differentiating both sides of (16) and using the help of (20) yield If , for all , but no limit cycle exists in this case as we show. It seems to be contradiction there and the reason is that for , is not bounded any more. This also shows the previous results for cannot be simply extended to the double integral case here by just letting .
Remark 2. The limit cycle forms after the first switching but is marginally stable. This means the trajectory will not remain at the original limit cycle even if a very small disturbance is introduced. Nevertheless, it will form a new limit cycle with different amplitude and period immediately after the next switching.
Remark 3. The limit cycle is a parabolic curve with the symmetry of , where is the switching time for . Its amplitude and period are only decided by the initial conditions, which is also consistent with the result from [6].
Example 1. Consider a double integral plant under the relay feedback with , and . The initial conditions are assumed to be and . The simulated response of the system is shown in Figure 5. By Theorem 1(iv), the amplitude and period of the limit cycle occurred are also calculated. The comparison of the results from the theory and the simulation as shown in Table 1, where one can see the error, is quite minor.
3. Conclusion
In this paper, a double integral plants under relay feedback is addressed. Complete results have been established on uniqueness of solutions, existence, and stability of limit cycles and its amplitude and period. Comparing with the plant , in [4], analytical expressions for the amplitude and period of limit cycles with terms of system parameters are available, but the conditions on the existence of limit cycles are quite different. Reasons for such differences as well as some explanations are provided.