Approximate Prediction-Based Control Method for Nonlinear Oscillatory Systems with Applications to Chaotic Systems

The approximate Prediction-Based Control method (aPBC) is the continuous-time version of the well-known Prediction-Based Chaos Control method applied to stabilize periodic orbits of nonlinear dynamical systems. The method is based on estimating future states of the free system response of continuous-time systems using the solution from the Runge-Kutta implicit method in real-time. Some aspects of aPBC are evaluated in the present work, particularly, its robustness to low future states estimation precision is exempliﬁed

Chaotic sets are composed of an infinite number of unstable periodic orbits (UPOs) [10,11] and the stabilization of one of these orbits leads to periodic oscillation, possibly, reducing amplitude.Different chaos control methods have been developed aiming at stabilizing these orbits with low control effort using the main characteristics of chaos [9].
One classical chaos control method is the Delayed Feedback Control (DFC) proposed by Pyragas (1992) [12].This is a state-feedback method whose control signal is computed through the difference between delayed and current time measured system states.DFC was initially proposed for continuous-time and has many applications on discrete-time systems [13].Nevertheless, DFC has a well-known limitation proofed for discrete-time systems, the odd-number limitation [14][15][16][17], which means that DFC cannot stabilize orbits with an odd number of Floquet multipliers larger than +1.This limitation is questioned in literature for continuous-time systems with counter example [18][19][20].
Different modifications on DFC have been proposed to overcome the odd-number limitation [13] and one of the most interesting modifications is the Prediction-Based Control (PBC) [21] developed for discrete-time systems.Instead of the delayed state, it uses the future state one period of the target UPO ahead computed along the trajectories of the free system response as reference for the control signal.Developments of PBC are proposed always for discrete-time systems; for example, in [22] a method is proposed for tuning control gain without previous knowledge about the target period- UPO that leads to fast convergence rate of trajectories to the stabilized periodic orbit; in [23] it is proved that period-2 UPOs can be stabilized by a pulsed control signal reducing the required control effort for stabilization; and in [24] equilibrium points (or period-1 UPOs) are stabilized in the presence of multiplicative or additive noise.
The literature provides some results on a PBC-like strategy applied to continuous-time systems.In [25] it is proposed to use the system Jacobian matrix at each point of the trajectory, instead of system future states, for stabilization of equilibrium points.The results showed that this strategy does not lead, necessarily, to stabilization of free system equilibrium points or UPOs and control effort does not vanish in steady state, two main advantages of chaos control proposed by Ott et al. [26].In [27], the strategy proposed in [25] is incremented using Neural Networks also for stabilization of equilibrium points.
In fact, PBC has a practical limitation for continuous-time systems if it is applied as proposed for discrete-time systems [21], the need of future state values in real time.In [28], we proposed the approximate Prediction-Based Control (aPBC) with a methodology based on the implicit Runge-Kutta method and state estimation applied to predict future states for the free system in real time based on system model.The authors claim that aPBC is the continuous-time applicable version of PBC because it uses prediction and stabilizes free system UPOs, ideally, vanishing steady state control effort.On the contrary, predicting future states has some drawbacks because the proposed future state prediction scheme requires an increase in the closed-loop system order and consequently computational power for numerical integration.Moreover, it is also subjected to prediction model and real system mismatch and integration method precision for application.
Continuing the aPBC development and study, one of the main questions that remained in [28] is about the robustness of the method for low future state estimation precision.The forced van der Pol (vdP) oscillator (nonautonomous) and Rössler system (autonomous) are used in the present work for numerical examples on the aPBC robustness and tradeoff between estimation accuracy and computational cost.We present a methodology and performance indexes that aid finding this trade-off and show results evidencing that it is possible to find a lower bound for precision.Besides that, this work presents an UPO of the forced van der Pol oscillator with one of the Floquet multipliers larger than +1 that is stabilized by the aPBC and could not be stabilized by the DFC with constant control gain following the procedures proposed in [18][19][20].The optimization procedure proposed in [29] was also applied to the DFC without success.This is evidence that the aPBC maintain the advantage of discrete-time PBC of not being sensitive to the odd-number limitation.
The paper is divided as follows.The aPBC is reviewed in Section 3, generalizing its proposition.In Section 4, aPBC application is presented using orthogonal collocation method [30,31] as implicit Runge-Kutta method.In Section 5, the forced van der Pol (vdP) oscillator and Rössler system and their chaotic behavior are presented.Both systems are used for numerical examples on the aPBC robustness and tradeoff between estimation accuracy and computational cost in Section 6.The example of an UPO of the forced vdP oscillator that could not be stabilized by DFC and is stabilized by aPBC is shown in Section 6.1.4.
The function  : R + × [0, ] → R  is such that (, 0) = () and (, ) = ( + , , (), 0).Condition (6) ensures zero control effort when the trajectory is on the unstable periodic solution  * of the free system.However, applying (, ()) in real systems is based on the exact prediction and real-time computation of future states ( + , , (), 0).Due to these practical constrains, alternatively, we are also interested in stabilizing an orbit close to the UPO  * of (1) resulting in a low control effort.This task is solved using the aPBC reviewed in the sequel.

Principles of the Approximate Prediction-Based Control Method (aPBC)
Computing the future state of the prediction term ( + , , (), 0) requires solving the free system ODE at each time , from time  to  + .This cannot be done exactly in real time and its solution is the basis of the aPBC.This task is divided into two steps, the first one consists in approximating the solution of ( + , , (), 0) by an implicit Runge-Kutta (R-K) ODE integration method, and the second is expressing the solution of the R-K method as a state observer that can be integrated by any explicit method in real time.

Approximation of the Prediction Term: 1st
Step.The first step consists in approximating the solution of (7b) by an implicit R-K ODE integration method [32], in order to estimate the prediction term, that is, the terminal value,  (, ) =  +  ∫ In order to estimate (, ) given by ( 9), the state transition map of the free system ( 2 ,  1 , , 0) is first approximated by the operator z defined by () =  ( +   ,  + where   : R → R  and   ,   ∈ R are weights chosen according to the implicit method used [32].The approximation z( + , , ()) of (, ),  ∈ [0,1], is calculated at the discretization points  =   ,  = 1, . . ., .  shall be chosen accordingly such that the sums in (10a) and (10b) lead to an adequate approximation of the integral in (9) for the given (, , ).
For simplicity, (10b) is written in the vector form where and   : R × R  × R  → R  is defined by In order to compute z( + , , ) through (10a), it is necessary to solve the algebraic system of (11) with unknown () ∈ R  . Writing The real-time solution of the DAE (16a) and (16b) requires the computation of (), its algebraic term, at each time .We therefore introduce an observer equation in the sequel, to transform the controlled system into a system of ODEs.

Approximation of the Prediction Term: 2nd
Step.We now approximate (16b) by solving the -dimensional ODE (17) whose solution L() is an estimation of ().The initial value L(0) is intended to be (precisely) computed off-line to provide good tracking quality for ().
The scalar gain   is chosen positive in order that L() of (17) tends asymptotically towards the solution () of (16b) when  → +∞ and typically in such a way that the estimator dynamics are faster than the controlled system dynamics.If the evolution of L() may indeed be chosen in order to fulfill (17), convergence does occur.
From (17), we obtain +  2   (,  () , L ()) (1  ⊗  (,  () , 0)) where  = (  ) and   is the partial derivative with respect to the th variable.  is the  ×  identity matrix and 1  is the column vector of dimension  with all elements equal to 1. ⊗ is the Kronecker product.
By determining L () from ( 18) we get Clearly, solving (19) in order to obtain (17) requires invertibility of the first factor.
We then define   : From (16a), (16b), (17), and ( 20) and denoting l () as the components of L(), the control law yields the following closed-loop system of ODEs: The solution of ( 21) is an approximation of the solution of the PDE given in (7a) and (7b) and therefore only an approximate stabilization of the initial orbit is expected or rather the stabilization to an orbit close to the initial one.
The ODE (21) has two types of state components, corresponding to the controlled system dynamics and to the dynamical state controller.Once L() stands for a set of unmeasured state variable components, (17) can be interpreted as a state observer.Notice that this estimator introduces dynamical feedback whose state has a dimension equal to the number of points of the R-K method adopted multiplied by the dimension of the initial system to be controlled ().

aPBC Application Issues
Two aspects of the aPBC shall be considered for application.One is the choice of an implicit R-K method to solve (10a) and (10b) and the other is the control gain tuning for stabilization.Here, we use the orthogonal collocation method with Lagrange polynomials and a constant control gain  computed through the closed-loop monodromy matrix described in the sequel.

Orthogonal Collocation Method with Lagrange Polynomials.
The implicit R-K method given in (10a) and (10b) is a general formulation used for the integration of differential equations whose application depends on the choice of a specific implementation.Herein, using the results in [28], we apply the orthogonal collocation [30][31][32] as implicit R-K method.
Collocation methods amount to approximate the prediction term by z( + , , ()), where z( + , , ()) is defined on the whole interval  ∈ [0, 1] by where   : R + → R  ,  = 1, . . ., , is a column vector.The functions   () are the Lagrange polynomials attached to the choice of the points The link with the implicit R-K method is as follows.
Equation ( 21) is obtained applying the substitutions provided by (22) and Theorem 1 with interpolating times   obtained by solving (25).

Closed-Loop Monodromy Matrix.
The aPBC was developed for the general control gain (, ()) of system (21).However, methods to choose the time and state dependent control gain have not yet been studied.Moreover, the control gain is tuned for the PBC method and applied to the aPBC, justified by the precision of the estimation of future states [28].For simplicity, a constant control gain  is applied that depends upon the ability of computing the closed-loop monodromy matrix of  * ().The computation of this matrix requires the integration of the closed-loop system and its variational equation ( 26 ∇  ( + , , , 0) is the free system monodromy matrix.Using ( 26) and ( 27), we compute the closed-loop monodromy matrix of  * () given gain .The Floquet multipliers are computed to measure the local stability of the controlled orbit for the chosen .
In practice, we fix , compute the monodromy matrix by integrating (26) with an explicit R-K method, and ( + , ,  * (), 0) is computed by integrating the free system over a period  at each step of the integration of (26).After that, we obtain the corresponding Floquet multipliers of the closedloop system.

Case Studies
Two continuous-time systems are used as case study for the aPBC-the forced van der Pol oscillator and Rössler system.In this section we provide a brief analysis about them evidencing their chaotic behavior through bifurcation diagrams and chaotic attractors.

Forced van der Pol
Oscillator.The forced van der Pol (vdP) oscillator is described by the nonlinear state space model (28): where  : R + → R 2 and  : In the present work we set  = 0.45 and  = 1 and provide bifurcation diagrams by attractor Poincaré map (Figure 1(a)) and largest nonzero Lyapunov exponent ( max , Figure 1(b)) for  = [0.9850.989] evidencing the route to chaos.Poincaré section is defined as usual for nonautonomous systems using time, and in this case where ( mod ) is the modulo operation or the remainder of the integer division of  by .Form Figure 1 it is possible to note that increasing  we have a series of period-doubling bifurcations leading to two chaotic attractors that are merged by a merging crises bifurcation.Lyapunov exponents are computed for each  evidencing the bifurcations and chaotic attractors.This work follows [29] and adopts  = 0.988 with chaotic attractor in continuous-time and Poincaré map shown in Figure 2.

Rössler
System.The Rössler system is described by the nonlinear state space model where  : R + → R 3 and  : R

Numerical Results
Numerical results are presented providing illustrative examples about the aPBC characteristics.We use the forced van der Pol oscillator as case study to show an UPO that cannot be stabilized by the DFC method due to the odd-number limitation and is stabilized using the aPBC.Using the same dynamical system, the robustness of the aPBC is evaluated for low-precision future states estimation.The aPBC is applied to the autonomous Rössler system and its robustness is also evaluated.For all the cases, the values ,   , and  that characterize the aPBC are tuned prior to application.

Applying aPBC to the Forced van der Pol
Oscillator.In this section, we provide numerical results on the approximated Prediction-Based Control (aPBC) applied to the forced van der Pol (vdP) oscillator according to (32), where (, ()) is defined in (28).
Herein, we apply the aPBC using the orthogonal collocation method (see Section 4.1), obtaining z( +   , , ()) =   (), and estimator (17) that leads to m () ≈   ().The control signal of the closed-loop control method (21) applied to system (32) is with the constant control gain  : R + × R 2 → R 2 is defined with the same dimension of the free system, but only the first component of the control signal vector is different from zero.This choice of matrix  makes only the first state variable influence on feedback and simplifies its tuning; however, better results can be obtained for a full matrix .
The control signal is designed by approximating the feedback term, leading to the closed-loop system (21).The solution of the latter depends on the choice of three parameters: the control gain , the observer gain   , and the number of discretization points for the orthogonal collocation method .The values   and  are directly related to the estimation quality and these two parameters will be tuned first.
The superscript on m(,  ) , () indicates the value of  and   used in the estimation and the subscript  refers to the state variable of the -dimensional vector m ().independent of (, ()).The values of  < 72 lead to integration instability and are not computed.
Repeating the same process, we evaluate (152,   ) for different values of   ; see Figure 5(b).Increasing  and   results in larger computation effort and trade-off between estimation quality and computation effort should be considered.

(𝑡).
Once the computational cost of the aPBC comes from processing the set of  ODEs used for prediction with reasonable error tolerance, m(102,10)  () is also a reference for a good trade-off between future state estimation and computation cost.
6.1.2.Target Unstable Periodic Orbits.Two UPOs are chosen as target UPO (TUPO) to be stabilized by the aPBC.TUPO1 was used in [29] for the application of an optimal DFC and in [28] for comparison between aPBC and optimal DFC.Herein, it is used to exemplify the effect of low-precision future state estimation through a small number of collocation points .TUPO2 is used to exemplify aPBC stabilizing an orbit subjected to the odd-number limitation when applying DFC.
Figure 7 shows the state space of the target periodic orbits, in continuous-time (line) and Poincaré map ( * ), embedded in the chaotic attractor.One point of the orbits with adequate precision for open-loop simulations and their Floquet multipliers are presented in Table 1.Period-1 (where period-1 means one cycle of the forced term or a one-point Poincaré map for the nonautonomous system) orbits of the forced van der Pol oscillator with  = 0.45 leads to  = 2/0.45. = 102 and   = 10 that result in a good estimation of future states according to value (,   ) have been obtained.Results for these  and   and a comparison with DFC are shown in [28].In the present section, the aPBC is applied to stabilize TUPO1 using  < 102 exemplifying the robustness of the method to the use of low-precision future state estimates, aiming at possible experimental applications.

Stabilizing TUPO1 through aPBC for Low-Precision
Before applying aPBC to TUPO1, a control gain  that leads to closed-loop stabilization shall be found.
Figure 8 shows the magnitude of the two Floquet multipliers || for different values of .Stability is achieved when |  | ≤ 1,  = 1, 2. We choose  = 0.25, which results in  || max ≈ 0.01213.All the Floquet multipliers in this case are real numbers.
Applying the aPBC for a Low . Figure 5(a) shows an increase of (, 50) when decreasing  for  = 0.This value refers to the free system and now we are then interested in characteristics of the closed-loop system, particularly the value of steady state error and control effort as a function of .
Steady state error index  is the distance (Euclidean norm) between predicted state m () and current state () for an entire cycle in the close vicinity of target UPO.Steady state control effort index  is the measure of total external signal (norm-1) necessary to keep the trajectory in close vicinity around one cycle of the target periodic orbit.It can be observed in Figure 9 that even for a reduction in the number of collocation points, stability of (an orbit close to) TUPO1 is achieved for  ≥ 67.TUPO1 stabilized with the aPBC for  = 82 and  = 67 and (0) = [−1.5 1.0]  are shown in Figures 10 and 11, respectively.Figures 10(e) and 11(e) show the initial value for the predicted states   (0),  = 1, 2, . . ., , and (, 0, (0), 0).It is observed that   (0) is relatively close to (, 0, (0), 0).

𝑑 = lim
The high values observed in Figure 9 for  = 3 are explained in Figure 12.It exemplifies a situation wherein a stable orbit different from TUPO1 is achieved due to the reduced .High amplitude periodic control signal in Figure 12(b) indicates that this stabilized orbit is not an orbit of the original free system.Figure 12(e) shows that  = 3 leads to   (0) almost unrelated to (, 0, (0), 0) and consequently  2 () = ( m,2 () −  2 ()) is a high amplitude control signal.According to Figure 12(a), a trajectory with initial condition (0) on TUPO1 (Poincaré point) diverges from it showing that it is not stable in this case.
It is noticed from the time-series of trajectories in Figures 10(c), 10(d), 11(c), and 11(d) that in certain torsions of the trajectory of the system controlled by the aPBC there is an oscillation on m ().This oscillation is generated because small  leads to instability on the predicted state and it reflects in time-series of control signals in Figures 10(b) and 11(b) and on the values of  and  in Figure 9.
Another consequence on reducing  is the loss of stability of m (), exemplified in Figure 13 and occurring even after changing the error tolerances of the explicit R-K integrator (the simulations present in this work were performed using the MATLAB5 integration routine ode45 and, in this specific case, the routines ode113 and ode15s were also tried, all of them resulting in interruption of the numerical integration).One of the reasons for this loss of stability is that  is the parameter that rules the implicit R-K integration error and Figure 13(d) shows that   (0) slightly diverges from (, 0, (0), 0) at time  and, contrasting to example of Figure 12, the trajectories of the new system do not converge to a stable periodic orbit.

Stabilizing TUPO2 through aPBC: An Odd-Number Limitation Example.
One of the main purposes when developing PBC for discrete-time systems is the odd-number limitation subjected by the DFC [14,21].Herein, we show a case where the continuous-time DFC with a constant control gain is subjected to the odd-number limitation and aPBC is not.TUPO2 has an odd number of real Floquet multipliers larger than +1 (see Table 1) and is used for a numerical example.m j,1 (0) m j,2 (0)  1 (T, 0, x(0), 0)  2 (T, 0, x(0), 0) time-series of the actual and predicted state variables; (e) time-series of (, 0, (0), 0) and the initial value of the predicted states   (0).x 2 (t) time-series of the actual and predicted state variables; (e) time-series of (, 0, (0), 0) and the initial value of the predicted states   (0).m j,1 (0) m j,2 (0)  1 (T, 0, x(0), 0)  2 (T, 0, x(0), 0) time-series of the actual and predicted state variables; (e) time-series of (, 0, (0), 0) and the initial value of the predicted states   (0).Applying the DFC.The DFC with control signal ∈ R 2×2 , is applied to system (32) with the goal of stabilizing TUPO2.The procedure of finding  for DFC follows [29].
A scalar control gain  is used with  defined in (34) resulting in the root locus chart of Figure 14.This chart contains the free system orbit Floquet multipliers  marked with ×: in black solid line, the fourth largest magnitude (the closed-loop system is infinite-dimensional when applying DFC; see [29] for details) closed-loop Floquet multipliers of the orbit for positive  and, in dark gray, the largest closedloop Floquet Multiplier of the orbit for negative .Arrows indicate increasing ||.It is possible to conclude that for any  there is a Floquet multiplier with magnitude larger than 1 and the orbit is not stabilizable.In addition, all the Floquet multipliers of the orbit tend to +1 when increasing DFC positive scalar gain .The same chart was done for the other three possible scalars  with equivalent (omitted) results.
Following [29], a matrix gain was tried with values obtained by an optimization process using the MATLAB routines fminsearch and fmincon [35][36][37][38].This optimization aims at finding a minimum for the largest magnitude Floquet multiplier through changing matrix .A scan process varying all the values of the matrix gain was also tried.The final trial used a matrix gain composed of a rotational matrix (formed by sines and cosines) multiplied by a real scalar gain equivalent to the gain used in [18][19][20], varying both the rotation angle and the gain.After performing the three methodologies, a matrix control gain that stabilizes TUPO2 was not found.
Applying aPBC and Stabilizing TUPO2. Figure 15 shows the magnitude of the Floquet multipliers as a function of the control gain  for the TUPO2 by applying the aPBC.It was noticed that the aPBC stabilizes orbits with an odd number of real Floquet multipliers larger than +1.One detail is relevant here; the stabilization is achieved only for negative values of the gain  and positive values of  make the largest magnitude Floquet multiplier increase.

Applying aPBC to the Rössler System.
In this section, we provide numerical results on approximated Prediction-Based Control (aPBC) for the Rössler system using control scheme (21) applied to system (30) according to (32).The goals are (i) to show that the method is applicable to autonomous systems without any modification and (ii) to compare the effect of different numbers of collocation points  on stabilization.aPBC is applied following the same procedure shown in Section 6.1 with control signal (, ()) defined in (33).Herein, we use the constant control gain This choice of matrix  makes only the third state variable influence on feedback and simplifies its tuning.In this case  < 32 values lead to numerical integration instability, except when  = 3, and are not computed.Herein we notice that less collocation points are necessary to obtain estimation quality equivalent to the one obtained for the vdP oscillator (see Figure 5).It can be seen that for the vdP there is (102, 50) ≈ 4.4 × 10 −4 , while for the Rössler system there is (52, 50) ≈ 4.4 × 10 −4 .Values (152,   ) for large   of the chosen trajectory of the Rössler system also stabilize to a value much smaller than for the vdP oscillator.
Figure 18 shows  2 () and m,2 ( − ) (beware of the time shift here) for the Rössler system, which can be compared with Figure 6.Notice that  = 5.8439698764 (the period of the TUPO4) for the Rössler system while  = 2/0.45for vdP oscillator, but the reduction on (152,   ) is much more related to the behavior of the trajectories than the period itself, because collocation points   are computed for a normalized time between 0 and 1.

Target Unstable Periodic
Orbit.Asenjo et al. [34] present one UPO of the chaotic attractor for  1 =  2 = 0.2 and  3 = 4.5.We apply the Newton-Raphson method with Poincaré section (31) to refine the position and period of the orbit and call it TUPO4 (Figure 19).One point of TUPO4 and its period  with adequate precision for open-loop simulations and Floquet multipliers  are presented in Table 2.

Stabilizing TUPO4 through aPBC.
This section is dedicated to the stabilization of TUPO4 using aPBC for different values .From Figure 17, values  = 102 and   = 10 are indicated for an accurate estimation of future states and aPBC will be initially applied for these values.Results for low  are also presented, evidencing the robustness of the controller for low-precision future states estimation in autonomous systems.
The control gain  of (38) shall be tuned for TUPO4 prior to the application.Herein, we observe that  ≥ 102 leads to a precise estimation once there is a convergence of  to the value computed using  = 152 and   = 50 (dashed line).Comparing Figure 21 with Figure 9, a great reduction on  and  is observed for the control method applied to the Rössler system, which is in accordance with what was observed for (, 50).In fact,  between 47 and 52 leads to a value of  equivalent to the value obtained for  = 102 for the vdP oscillator.Considering , this reduction is more significant, with  = 37 leading to a value of  close to the value    1 (T, 0, x(0), 0)  2 (T, 0, x(0), 0)  3 (T, 0, x(0), 0)   1 (T, 0, x(0), 0)  2 (T, 0, x(0), 0)  3 (T, 0, x(0), 0) obtained for  = 102 for the vdP.This result is important for practical aspects because it shows that, for certain systems, the augmented state of the closed-loop system controlled with the aPBC can be reduced significantly (in comparison with the initial results shown for the vdP oscillator).

Journal of Control Science and Engineering
Figure 21 shows that stability of TUPO4 is achieved for  ≥ 32.The TUPO4 stabilized by the aPBC for  = 102,  = 52, and  = 32 is shown in Figures 22,23  .Deviation of m,2 () from  2 () is observed in Figure 24(d), which reflects the numerical integration instabilities for  < 32.We verify that stabilization is not achieved for  = 3, which leads to the high values of  and  observed in Figure 21.

Conclusions
Approximate Prediction-Based Control (aPBC), which is the continuous-time version of Prediction-Based Control (PBC), was reviewed in this work and some of its characteristics were analyzed through numerical simulations.It was shown that continuous-time application of PBC requires future state values and aPBC is based on a methodology for estimating them in real-time.
The main drawback in aPBC is that estimation of future states requires extending the order of the closed-loop system by the number of implicit Runge-Kutta solution points, herein the orthogonal collocation points, times the number of state variables.Applications require real-time explicit integration of these extended equations and reducing them implies reducing the computational cost while reducing estimation precision.Both exemplified cases have shown that there is a lower bound on the number of collocation points to achieve stability, to both target UPO and numerical integration.While comparing examples, it was noticed that the number of collocation points depends upon the trajectory torsions in state space.Future studies should be carried out in order to reduce the number of implicit Runge-Kutta solution points through different orthogonal polynomials or different implementation methods.
Another interesting example is an UPO of a continuoustime system that could not be stabilized by DFC with constant control gain because of the odd-number limitation.The same orbit is stabilized by aPBC with a negative scalar constant gain without any modification on the method.Then, the main justification for the application of PBC in comparison to DFC in discrete-time is also valid for continuous-time.
The method proposed for estimation of open-loop future states of nonlinear systems applied to aPBC may be generalized for other applications.As example, model-based predictive control [39] applied to nonlinear systems requires open-loop future state values to optimize control gain at each time step and this routine may be fastened by using the proposed scheme.

Figure 1 :
Figure 1: Bifurcation diagram of forced van der Pol oscillator for  = 0.45 and  = 1 using (a) Poincaré map and (b) largest nonzero Lyapunov exponent.

Figure 5 :
Figure 5: Relative error of the free system response estimated future value using different (a)  and (b)   for the forced vdP oscillator.

6. 1 . 1 .
Tuning the aPBC Parameters   and .In order to evaluate the future state estimation, whose characteristics are related to the parameters   and , we first set  = 0.The system is integrated for  ∈ [0, ] and the points m (+/),  = 0, . . ., , are collected for different values of   and .One then assesses the convergence by comparison with the solution obtained with  = 152 and   = 50, by computing  (,   ) =      m(152,50)

Figure 7 :
Figure 7: Target UPOs (TUPO1 and TUPO2) embedded in the chaotic attractor (CA) of the forced vdP oscillator.The points for the target UPOs and the chaotic attractor in the Poincaré map (PM) are represented by asterisk and gray bullet, respectively.
Future State Estimation.In Section 6.1.1,the values of

Figure 8 :
Figure 8: Magnitude of the two Floquet multipliers for different values of  for the TUPO1 controlled by the aPBC.

Figure 9 :
Figure 9: Indexes (a)  and (b)  when applying aPBC to TUPO1 for different  with  = 0.25 and   = 10.Bottom dashed line are the values of (a)  and (b)  computed with  = 152 and   = 50 used as reference.

𝑡→+∞ 36 )Figure 9
Figure 9(a) shows  when stabilizing the TUPO1 by applying the aPBC for different values of  with  = 0.25 and   = 10. Figure 9(b) shows the same for .It can be observed in Figure9that even for a reduction in the number of collocation points, stability of (an orbit close to) TUPO1 is achieved for  ≥ 67.TUPO1 stabilized with the aPBC for  = 82 and  = 67 and (0) = [−1.5 1.0]

Figure 14 :Figure 15 :
Figure14: Root locus chart for the Floquet multipliers of the TUPO2 using the DFC.In black, the four largest magnitude Floquet multipliers for positive control gain .In dark gray, the largest magnitude Floquet multiplier for negative control gain .The Floquet multipliers of the free system orbit are represented by (×) and the unit circle by dashed line in black.
Figure 20 shows the magnitude of the three Floquet multipliers || for different values of  according to the procedure presented in Section 6.1.3.The value  1 = 1 is characteristic of autonomous system trajectories and stability is achieved when, for the other Floquet multipliers, |  | ≤ 1,  = 2, 3. Choosing  = 1.85 results in | 2 | ≈ 2.994 × 10 −3 (and | 3 | ≈ 2.167 × 10 −14 ).All Floquet multipliers in this case are real numbers.Applying the aPBC for Different Values of .In this section, we evaluate the characteristics of aPBC for different values  applied to the stabilization of TUPO4.Herein, we use  ≤ 102,   = 10, and  = 1.85.

Figure 21 :
Figure 21: Indexes (a)  and (b)  when applying aPBC to TUPO4 for different  with  = 1.85 and   = 10.The bottom dashed lines are the values of (a)  and (b)  computed with  = 152 and   = 50 used as reference.

Figure 17 (
Figure17(a)  shows an increase of (, 50) observed when decreasing  for  = 0. Figure21shows  and  for different values of  with  = 1.85 and   = 10.Values of 3 <  < 32 were not computed because of instability, the same as observed in Section 6.1.3.Herein, we observe that  ≥ 102 leads to a precise estimation once there is a convergence of  to the value computed using  = 152 and   = 50 (dashed line).Comparing Figure21with Figure9, a great reduction on  and  is observed for the control method applied to the Rössler system, which is in accordance with what was observed for (, 50).In fact,  between 47 and 52 leads to a value of  equivalent to the value obtained for  = 102 for the vdP oscillator.Considering , this reduction is more significant, with  = 37 leading to a value of  close to the value

Table 1 :
Floquet multipliers and one Poincaré map point for  = 0.988 of the two target periodic orbits of the vdP oscillator.