Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters

Information fusion integrates aspects of data and knowledge mostly on the basis that system information is accumulative/distributive, but a subtle case emerges for a system with bifurcations, which is always un-linearizable and exacerbates information acquisition and presents a control problem. In this paper, the problem of an un-linearizable system related to system observation and control is addressed, and Andronov–Hopf bifurcation is taken as a typical example of an un-linearizable system and detailed. Firstly, the properties of a linear/linearized system is upon commented. Then, nonlinear degeneracy for the normal form of Andronov–Hopf bifurcation is analyzed, and it is deduced that the cubic terms are an integral part of the system. Afterwards, the theoretical study on feedback stabilization is conducted between the normal-form Andronov–Hopf bifurcation and its linearized counterpart, where stabilization using washout-filter-aided feedback is compared, and it is found that by synergistic controller design, the dual-conjugate-unstable eigenvalues can be stabilized by single stable washout filter. Finally, the high-dimensional ethanol fermentation model is taken as a case study to verify the proposed bifurcation control method.


Introduction
In nonlinear systems covering ecology, economics, sociology and many other fields, complex dynamic evolutionary behaviors, such as multiplicity, self-oscillation, and even chaos arise in the vicinity of bifurcation points [1], and various methods are adopted to enhance sensory ability for accurate manipulation and control: the bifurcation oscillator is applied as a quantum state amplifier [2], the oscillating clearance of the rolling-element bear is taken as the bifurcating parameter with repetitive errors fused into the control system [3], and deep learning is implemented for the control of the Van del Pol system [4]. Mathematically, a bifurcation takes place at a parameter value where the system loses structural stability with respect to parameter variations, i.e., a phase portrait around the equilibrium is not locally topologically conjugate to the phase portraits around the equilibrium at nearby parameter values [5,6], which leads the solution manifold to change significantly, and exacerbates system information acquisition and control implementation [7,8].
Bifurcation control [9] refers to the task of designing a controller to suppress or reduce existing bifurcating dynamics, thereby achieving desirable dynamical behaviors, i.e., to obtain an asymptotically stable phase portrait. Feedback is often applied on a system with bifurcations, where the methods of linear/nonlinear feedback, washout-filter-aided feedback, impulse control, vibrational control, delayed feedback [10], etc., are frequently adopted. Nonlinear analysis and control based on harmonic balance approximations or nonlinear frequency response is also studied [11]. Recently, the idea of the anti-control of bifurcations [12] is proposed to segregate unwanted dynamics by deliberately introducing simple bifurcations into the system. tion f ; then, a linear/linearizable system starts as a weighted contribution of individual variables Furthermore, the coupling relation between causes and outcomes is decoupled through variable substitution and combination, e.g., a 2 × 2 input-output system with a linearized form is given as follows, Assume that the Jacobian matrix is non-singular, and a coordinate change y = Ax exists that could diagonalize Equation (2) where : Equation (3) indicates that for a linear/linearizable system, outcomes can be departed as a sum of the weighted causes separately, and sets of causes can be found to correlate every outcome exclusively. The above linear system properties are held for the linear-time-invariant (LTI) dynamic system. Given in state-space form, a multi-input-multi-output (MIMO) system undergoes Laplace transform and obtains similar distributive and decoupling properties where ξ is the nonlinear terms that are truncated out in the linearizable dynamic system, whilst y and u represent outcomes and causes, respectively. Mention that many nonlinear controls are developed on the linearized system. In practical applications, the presence of system nonlinearities causes the dynamic behavior to be qualitatively different from one operating regime to another; hence, one might obtain linear (at a specific set-point) models valid within a "small" region about the linearization point, and perform local designs for a set of operating conditions and then construct a gainscheduling scheme that interpolates controller gains as the process traverses the operating region.
Accordingly, un-linearizable dynamic systems are those truncating out ξ as it alters the structural stability of the original system, and many dynamic systems with bifurcation are un-linearizable. Since a system under Andronov-Hopf bifurcation has dual-unstable eigenvalues with planar phase portrait, which is a typical multivariate, un-linearizable case, the following analysis and control is focused on the systems with Andronov-Hopf bifurcation.

Normal Form Andronov-Hopf Bifurcation: A Typical Un-Linearizable System
A standard approach to study the behavior of ordinary differential equations (ODE) around a bifurcation point is through normal form analysis, which starts by center manifold degeneracy. For a parameterized dynamical system, where x → f (x, α) ∈ R n and α is the bifurcation parameter, with equilibrium x 0 satisfying f (x 0 , 0) = 0. Assume eigenvalues of the Jacobian about x 0 are λ 1,2 = ±iω, λ k < 0(k = 3, 4, . . . , n), then, there exists n-2 invariant/center manifold that is tangent at the bifurcation point to the eigenspace of the neutrally stable eigenvalues λ 1,2 . Specifically, for a small neighborhood of α ∈ (0, δ) and within the invariant manifold exists a 2-dimensional attracting manifold and the coordinate change of the first two eigenspace; by introducing complex conjugates z and z * , Equation (5) around x 0 can be given as where g(z, z * ) is the degenerated terms at α = 0 and can be given explicitly through rules of implicit function theorem. Assume S(z, z * ) ∈ R (n−2) represents the stable manifold, the n-2 stable equivalence is provided as follows and the n-2 solutions are exponentially convergent to S(t) where β is the positive convergent rate. Considering that x 3,...,n are locally asymptotic and therefore can be neglected in a local stability analysis around the bifurcation point, in the following discussion, about the Andronov-Hopf bifurcation, a two-dimensional topologically equivalent eigenspace containing dual-unstable (λ 1,2 = ±iω) is discussed instead of Equation (5). For a general two-dimensional parameterized system, i.e.,x = (x 1 , where : where α, b, c and d are elements of the Jacobian A(α), which is simplified as A; h(x, α) represents the higher-ordered term. The characteristic equation corresponding to A gives where tr(A) and det(A) are the trace and determinant, respectively. Andronov-Hopf bifurcation takes place at α = 0 satisfying tr(A) = 0 and det(A) = ω 2 0 . Assume in a small neighborhood α ∈ (0, δ) that the eigenvalues are written as λ 1 (α) = µ + iω and λ 2 (α) = µ − iω, where µ = tr(A) and ω 2 = tr(A) 2 − det(A); then, T(α) exists to transform A(α) into canonical real form as follows where J −1 = J T /(µ 2 + ω 2 ). Let y = T(α)x and Equation (9) turns to Introducing z = y 1 + iy 2 and z * = y 1 − iy 2 and Equation (12) presents one-dimensional equivalence as follows z = λ(α)z + g(z, z * ), (13) with g = o(z 2 ). Let q(α) be the eigenvector of λ 1 about A(α) and one can give an eigenvector p(α) corresponding to λ 2 . Since λ 1 and λ 2 are complex conjugates, p and q can be chosen so that the inner product satisfies q, p = 1, and p, q * = p, 1 Equation (11) after transposition has T T (α)J T = A T (α)T T (α), and if p(α) has the following form where T ij is i-th line, j-th row element , then z has the explicit form as follows Considering q, p = 1 and Equation (13), the original variable x has the only equivalent transform Substituting Equation (17) into (16), one obtains Here, the nonlinear term F(x, α) in Equation (9) after Taylor expansion gives Substituting Equation (19) into (18) and the second-order term B gives Through a coordinate change in Equation (21) the second-order term vanishes and gives Equation (22) Similarly, the third-order term after coordinate change obtains Next, let ω = ρe iφ and Equation (23) turns to , where : Suppose that ϕ performs constant rate rotation and Equation (24) mapped to the rotating plane gives where : Assume that ρ 0 = 0 is the initial condition, solving Equation (24) that provides For a fixed α > 0, there exists a limit cycle, and when ϕ = 2π, the state ρ 1 would turn back to the initial point Comparing Equations (26) and (27), it is found that the 4th-or higher-order terms are unaffected by the formation of limit cycles, and hence can be truncated out. Furthermore, an Andronov-Hopf bifurcation is structurally equivalent to the following standard form As shown in the following schematic diagram, if Equation (28) is linearized by truncating out the nonlinear terms, the phase portrait shows that the outcome is asymptotically stable to the origin ( Figure 1a); if the nonlinear terms at Equation (28) are preserved, the phase portraits of the Andronov-Hopf bifurcation would be super-critical (Figure 1b), where the periodic solution is stable, or sub-critical (Figure 1c), where the periodic solution is unstable, and the criterion for super-or sub-criticality is judged by the sign of the third-order term of Equation (28).
It is an obvious standard form of the Andronov-Hopf bifurcation system that differs from the linearized system as it either attracts or repels an oscillatory trajectory, while the linearized one exhibits a diverging trend in the whole phase portrait. Hence, the nonlinear terms in Equation (28) are indispensable, and a system with Andronov-Hopf bifurcation is un-linearizable. Address that there exists two variablesx 1 andx 2 interacting with one another to form a limit cycle around the Andronov-Hopf bifurcation point (α > 0). Furthermore, Equation (28) can be represented one-dimensionally using complex variables, which means manipulating eitherx 1 orx 2 to affect the dynamic behavior of the system directly (or equally). Hence, introducing one state feedback suffices stabilizing a system with Andronov-Hopf bifurcation. Notion, the sign of the third-order term in Equation (28), decides whether the bifurcation is super-critical (−) or sub-critical (+), and one can refer to [21] for more details.

Washout Filter as a State Manipulator
Since washout filters reject steady-state signals and passing transient ones, it could be used as a state manipulator and by fine-tuned feedback, eliminates the instabilities originating from open-loop bifurcations. Figure 2 is the washout-filter-aided controller design scheme: by inserting washout filters between states and outputs, transient dynamic behaviors are separated out that can be used to cancel out the perturbations of the system through feedback; while outputs preserve the steady state information, this is important for a control system, meaning washout filters are inoperative until vibrational signals are detected. Assume that x k is the k-th element of the state variables, and d is the filter time constant. The control law used in this paper is quite simple, and v is the reference input used to determine the actual operating point and k is a controller parameter.

Feedback Stabilization Based on the Linearized System
To stabilize the normal-form Andronov-Hopf bifurcation, washout-filter-aided control is constructed. Without loss of generality, the linearized system is extended as follows where the control u = [u 1 , u 2 ] T takes u = v + ky, and k is the feedback gain; B is an arbitrary input-state transfer matrix and h(x, u) represents the higher-ordered-terms. Adding the auxiliary state satisfies dz/dt = y, and the closed-loop system is provided as follows Stabilizing Equation (30) requires the following augmented matrix A C being Hurwitz Theorem 1. If A is non-singular and (A, B) is stabilizable, there exists (P, k), which makes A C Hurwitz.
Proof. A C after similarity transforms of T 1 and T 2 gives Through the following small gain analysis, the control (P, K) is obtained Canceling the first-and second-order perturbational terms provides where : which dictates, Then, around M 0 in a small region, the only M satisfying Equation (33) is given with A C satisfying where : When (A, B) is stabilizable, A + BK can be designed Hurwitz. Since A is non-singular, setting P 1 = A −1 (A + BK), then −AP 1 (A + BK) −1 is Hurwitz, so as A C2 . Plus, proper feedback gain k can promise the good transient dynamics of the closed-loop system.
One can design P based on above perturbation theory, but how many washout filters needed to stabilize the system is still a problem, the following lemma provides some insights. Lemma 1. The parity of the unstable eigenvalues of −P needs to be consistent with that of A to make A C in Equation (31) be Hurwitz.
Therefore, to stabilize the normal form Andronov-Hopf bifurcation based on information of Equation (29), −P needs to be a 2 × 2 matrix. The Laplace transform between state and control gives u(s) Obviously, when the state x(s)is one-dimensional, the washout filter has d = P, but for a multi-dimensional system, e.g., P is 2 × 2, and the parameters of the washout filters are difficult to provide where : In accordance with the form of washout filters, the higher-ordered terms in Equation (39)  .
Meaning two washout filters with d = tr(P) aid the controller design, and Equation (40) is decoupled.

Bifurcation Analysis on the Closed-Loop System
Consider only one washout filter is introduced to x i with feedback u = K i y, and (i, j) pairing generates four possible cases for a two-input-two-output system. The normal form Andronov-Hopf bifurcation is augmented and h(x, u) explicitly presented where y is the washout-filter-output corresponding to x i . Based on Routh-Hurwitz stability criterion [25], one can deduce the control gain K i of the washout filter that stabilizes Equation (41), and the Jacobean of the augmented system gives Expending Equation (42) along the last column gives the characteristic polynomial P c One can design pairs of (K i , d) to make eigenvalues of [∼] in Equation (43) stable, but λ −i is kept unchanged. Meaning that, for the linearized system, dual-unstable eigenvalues need two washout filters.
However, Equation (41) is a nonlinear control, and the dynamics of f i could be compensated through feedback. The Jacobian from Equation (42) ignores nonlinear correlations between system and control feedback, which leads to redundant control loops. When numerical bifurcation analysis is implemented, it is found that only one washout filter can stabilize the closed-loop system represented by Equation (41).
One can design controllers through the identification of critical points for the emergence/elimination of unstable eigenvalues. For the case of Equation (41), the limit cycles when α > 0 need to be eliminated (or stabilized); hence, with the aid of a washout filter framed in Figure 2, Andronov-Hopf bifurcation points in the (K i , d) domain can belocated numerically (refer to Appendix A for more details).
Substantial configurations are provided in Figure 3, where the solid lines are Andronov-Hopf bifurcation curves. BT represents Bogdanov-Takens bifurcation, and ZH represents zero-Hopf bifurcation, both of which are codimension-2 bifurcations originating from Andronov-Hopf bifurcation curves. Since x 1 and x 2 in Equation (41) are somewhat symmetric, configurations C1 and C4 are identical, and so are C2 and C3. For each stabilizable region, d > 0 implies the analog signal after washout filters are stable, while d < 0 means that the feedback signal is unstable. As shown in Figure 3, the Andronov-Hopf bifurcation curves separate out parameter space (K i , d) that stabilizes the closed-loop system. Mention d = 0 in cases of C1 and C4 indicates that traditional PID control stabilizes the system. It is observed d < 0 is able to stabilize the system from a nonlinear analysis perspective, but the unstable washout filter generates diverging signals which are dangerous for real-time applications where interference is unavoidable and is not recommended. Exclusively for d > 0, it is found C2 and C3 have broader stable margins than C1 and C4, which implies manipulating variables and control feedback could be chosen such that stabilizing control is easy to achieve.
Based on an analysis in Figure 3, numerical simulations on the Simulink platform were conducted, where α = 0.1 for Equation (41). As is shown in Figure 4, when the parameter pair (−5, 0.5) is adopted for case C1, it is stabilizable. The perturbation output after washout filter damps down and so as the states x 1 and x 2 , and a positive phase shift between y and x 1 implies that the oscillatory dynamics of the control feedback and forward signal are well geared to stabilize the closed-loop system. Since x 1 and x 2 are coupled in the self-oscillatory system, stabilizing x 1 would cause x 2 to damp down. Further study of the information acquisition by the washout filter shows that, the stable one obtains oscillatory signals, while the unstable one is exponentially divergent, where the overall stable character of the closed-loop system is obtained by the subtraction of the forward and backward signals. Hence, to stabilize the self-oscillatory dynamic system, stable washout filters are suggested.

Bifurcation Analysis on the Continuous Ethanol Fermentation Model
One of the most important bio-industrial sectors is the production of ethanol, which can be used as a gasoline additive to improve vehicle performance and reduce carbon emission. However, complex microorganism growth kinetics and ethanol inhibition present challenges for continuous production at the industrial scale, which limits the contribution of bio-fuel to current energy supply, and hence an exquisite controller design is critical to stabilizing the process. The structured model for Zymomonasmobilis fermentation [26] in a continuous stirred reactor is provided as follows, which considers the mass balance of the substrate (C s , g/L), microorganism (C x , g/L), K-compartment biomass (C e , g/L), product (C p , g/L) and the related parameters are given in Table 1 dC where the dilution rate D and the substrate feed concentration C s f are set as the design /operation parameters because of practical concerns. D is the reciprocal of the residence time for the reaction stream, which has a fundamental impact on the bioreactor behavior; C s f always comes from upstream, which may be interfered with due to upstream uncertainties. With this model, nonlinear analysis is conducted on the open-loop system. A twoparameter bifurcation diagram decomposes the parameter space (D, C s f ) into different sections, and each has distinct dynamic properties. Furthermore, as is shown in Figure 5, the Andronov-Hopf bifurcation curve divides the parameter space into S1 and S2. Any designing point at S1 generates a self-sustained oscillator, while the point at S2 is in stable steady equilibrium. Particularly, when C s f is set as 140 g/L, H is detected at D H = 0.05024 h −1 . H is a supercritical Hopf bifurcation point due to its negative first Lyapunov coefficient (l p), and operating the process with D at the left hand side of H leads to self-oscillations, while stable equilibriums are obtained on the right-hand side. It is obvious that the highest C p is achieved when D is infinitesimal, which means that the reactor size needs to be very large. Since a high C p of the bioreactor outlet cuts down the separation cost of later processing, one should design the reactor by setting D as the small reasonable equipment cost. We set the design point for Equation (44) as P d (0.022, 140) in this work, which falls into the self-oscillatory section. The idea of designing a self-oscillatory process was inspired by [27,28]. Taking P d point as example, time average C p is 8.8% higher than the corresponding steady-state output of the same design condition. However, the presence of self-oscillatory dynamics brings a challenge for the control of the continuous system; therefore, bifurcation control with the aid of washout filters is demanded to stabilize this multi-dimensional system.

Washout-Filter-Aided Control of the Fermentation Model
The model in Equation (44) is four-dimensional, and dual-unstable eigenvalues exist when Andronov-Hopf bifurcations emerge. For P d , the equilibrium for C s , C x , C e and C p are 0.8509 g/L, 21.6547 g/L, 0.01915 g/L and 58.71 g/L, respectively; and the corresponding eigenvalues are −2, −0.07136, 0.01232 ± 0.1411. Dual-unstable eigenspace emerges between states C e and C p , indicating that the controlling variable would be chosen among them, C p is set as the washout filter input, and it is assumed that C p can be detected with online measuring devices. Then, the washout filter output is provided where d is the filter parameter, and to promise that y is stable and convergent, d needs to be positive. The manipulating variable is set to D, which is related to the input flow rate. Compared to the nominal D 0 = 0.022 h −1 , the feedback by P control is provided where k is the feedback gain.
Combining Equations (44)-(46), one could implement bifurcation analysis on the closed-loop system. As is shown in Figure 6, when the washout-filter-aided P-control is facilitated on the continuous fermentation model, the controlling space (k, d) is separated into three sections by the Andronov-Hopf and neutral saddle curves, where only the upperleft one is stabilizable. Note that both Andronov-Hopf and neutral saddle bifurcations are codimension-1 bifurcations, and the bifurcations happening on these two curves are codimension-2 bifurcations. For the Andronov-Hopf bifurcation curve, three codimension-2 bifurcations, GH1 (0.1745, 0.0024), GH2 (0.0623, 0.0006) and ZH (0, 0.0005) are detected; and for the neutral saddle bifurcation curve, HH (0.02221, −0.03037) is detected. Mention for a given filter parameter, e.g., d = 2, with the decrease in k, that the dynamics of the closed-loop system evolves from stable to saddle, and then diverges. Further study is conducted by setting the control parameters to constant, e.g., d = 2 and k = 1, and observing the changes in the dynamics. As shown in Figure 7, between the two Andronov-Hopf bifurcation curves is the segment that can be stabilized by washoutfilter-aided P-control. When the control is facilitated, stabilizing the designing point demands that C p reaches the steady state, and then the washout filter output y has a tendency to converge towards zero, and Equation (45) tends towards zero; hence, the closed-loop system has the same equilibrium as the open-loop one, and C p = 58.71 g/L is preserved at operational point P d . As provided in Figure 8, the simulation results show converged dynamics after control d = 2 and k = 1, where C p and D reach P d as time approaches infinity.

Conclusions
In the presented study, we have focused on a subtle case with bifurcations, which is always un-linearizable and exacerbates information acquisition and a control problem. Taking the system with Andronov-Hopf bifurcation as a typical example, normal form analysis reveals that the cubic terms are the integral part of the whole system, and with the consideration of the un-linearizable characteristics, the stabilization procedure adopting washout-filter-aided control could realize the migration of dual-unstable eigenvalues with only one washout filter, which is not possible by previous control strategies. Moreover, the continuous ethanol fermentation process is taken as a case study for stabilizing a high-dimensional Andronov-Hopf bifurcation system: (1) variables mapping to the dualunstable eigenspace are chosen as the manipulated variable; and (2) a stable washout filter parameter is demanded to stabilize the dual-unstable eigenvalues. By washout-filter-aided P-control, the ethanol fermentation process is stabilized and the equilibrium is preserved.
Author Contributions: All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by C.Z., C.Y. and J.N. The first draft of the manuscript was written by C.Z. and both authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: All data generated or analyzed during this study are included in this published article.

Conflicts of Interest:
The authors declare that they have no conflict of interest.

Nomenclature
D dilution rate, h −1 K S saturation constant, g/L k 1 , k 2 , k 3 inhibition constant, L/g T 0 inherent period, h P d design point for unforced system r e growth rate of the K-compartment, g/L h C s f substrate feed concentration, g/L t time, h x biomass concentration, g/L Y sx , Y px yield factor of biomass on substrate/product, g/g m s maintenance item, g/g h Greek letters δ Greek letters µ max maximum specific growth rate, h −1 Subscripts s substrate x biomass e K-compartment p product

Appendix A
The following test function (Stephanos theorem) locates the critical points numerically. Stephanos theorem: If a n × n matrix A has eigenvalues µ 1 , µ 2 . . . µ n ordered in an arbitrary order, then 2AΘI n satisfies (µ i + µ j ) n>i>j > 1, and the function ϕ H (x, α) = det(2 f x (x, α)ΘI n ) vanishes at an Andronov-Hopf bifurcation point, where "det" represents the determinant, and Θ is the bi-alternate product.
Proof. The eigenvalues of a matrix are preserved by similarity transforms as follows (PΘP)(2AΘI n )(PΘP) −1 = 2BΘI n where B is the transform of A, B = PAP −1 . Hence, one can assume that A is in upper triangular form (by reducing to the Jordan form if necessary). Then, the eigenvalues of A become its diagonal elements. By applying the equation below, α ik i f k = i and l = j α ii + α jj i f k = i and l = j α jl i f k = i and l = j −α jk i f l = i 0 else 2AΘI n is an upper matrix if its indices are ordered lexicographically, and its diagonal elements are given as Stephanos theorem claims. By defining the test equation ϕ H (x, α) = 0, the critical point α c can be located and inserted into the continuation sequence.