Next Article in Journal
Research on Non-Pooling YOLOv5 Based Algorithm for the Recognition of Randomly Distributed Multiple Types of Parts
Next Article in Special Issue
Research on Rapid and Low-Cost Spectral Device for the Estimation of the Quality Attributes of Tea Tree Leaves
Previous Article in Journal
SF-Partition-Based Clustering and Relaying Scheme for Resolving Near–Far Unfairness in IoT Multihop LoRa Networks
Previous Article in Special Issue
Research on a Non-Contact Multi-Electrode Voltage Sensor and Signal Processing Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Faculty of Mechanical and Electrical Engineering, Kunming University of Science and Technology, Kunming 650093, China
2
Faculty of Chemical Engineering, Kunming University of Science and Technology, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(23), 9334; https://doi.org/10.3390/s22239334
Submission received: 12 November 2022 / Revised: 25 November 2022 / Accepted: 27 November 2022 / Published: 30 November 2022
(This article belongs to the Collection Multi-Sensor Information Fusion)

Abstract

:
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.

1. 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.
However, bifurcation control sometimes exhibits contradictory outcomes: improper control actions or process–model mismatch might cause the closed-loop system to move across the stable manifold, leading to more complicated dynamic behaviors [13], and one cause is that many control laws are developed on the basis of linear system theory (gain-scheduling, feedback linearization, etc.), but most dynamic systems with bifurcation are un-linearizable, i.e., truncating out the nonlinear terms alters the structural stability of the original system, and the higher-ordered-terms might not shrink to zero as the time approaches infinity. On the other hand, normal form degeneracy provides an effective way to evaluate those nonlinear terms, where types of bifurcations are expressed with quadratic or cubic degeneracy, and state-feedback based on normal forms and invariant sets are developed [14,15].
The main goal of this paper was to address the problem of an un-linearizable system in bifurcation control, where a comparative study is conducted on both the original and linearized models using washout-filter-aided state feedback. Through washout filters, a feedback signal is expected to be well geared with the open-loop dynamics and generate asymptotic outcomes. Abed and co-workers [16,17,18,19] implemented washout filters for bifurcation control, and feedback stabilization theory is proposed through auxiliary state variables that are inter-coupled. For systems with varying bifurcation parameters, the adaptive control scheme is also developed [12]. Since washout-aided feedback is introduced to migrate out of unstable eigenvalues, Marquardt [20] suggested manipulating directly on the unstable eigenspace, and introducing a pair of complex conjugate eigenvalues instead. Zhang [21] further employed eigenvalue assignment to achieve bifurcation control with guaranteed transient dynamics. Chen [22] proposed a systematic way of choosing constants and control parameters for the control of a high-dimensional bifurcating system. However, for a complicated system, i.e., both single- and dual-unstable eigenvalues exist, and controller design could be rather sophisticated.
The aforementioned washout-filter-aided control was developed on the linearized system, and the linearization procedure might change the structural stability of the original system and cause redundant control actions. In this work, the elementary properties of a linearizable system are offered, which provides a comparison with the un-linearizable system when these properties fail. Then, the system under Andronov–Hopf bifurcation is taken as a typical case of an un-linearizable one and the related analysis and bifurcation control strategy are developed. Since washout-filter-aided control can migrate the unstable eigenvalues, numerical bifurcation analysis on the closed-loop system reveals that, by inserting only one washout filter to any of the dual-unstable eigenspaces, the system might be stabilized. Compared with previous developed control strategies [23], our method could downsize the numbers of washout filters, facilitating the application for more practical and general cases, especially for high-dimensional systems.
The inconsistency of feedback stabilization and numerical bifurcation analysis stems from the fact that many nonlinear systems with bifurcation are un-linearizable [24]. This work is organized as follows: Through the analysis of the normal form of Andronov–Hopf bifurcation, the properties of an un-linearizable system are detailed in this paper; then, washout-filter-aided feedback on the un-linearizable system is analyzed in contrast with previously developed washout-filter-aided control strategies; and the continuous ethanol fermentation model is taken as an example of the un-linearizable system and is discussed in detail, and finally, concluding remarks are provided last.

2. The Un-Linearizable System

2.1. Properties of the Linear System

Previously to the discussion on the nonlinear dynamic system, basic properties concerning the linear system are provided. Consider a system exhibiting outcomes (characters) M, and which is detected with finite numbers of sensible causes x 1 , , x n , related by function f; then, a linear/linearizable system starts as a weighted contribution of individual variables
M = f ( x 1 , x 2 , , x n ) f x 1 x 1 + f x 2 x 2 + f x n x n .
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,
M = M 1 M 2 = f 1 ( x 1 , x 2 ) f 1 ( x 1 , x 2 ) f 1 x 1 f 1 x 2 f 2 x 1 f 2 x 2 x 1 x 2 .
Assume that the Jacobian matrix is non-singular, and a coordinate change y = A x exists that could diagonalize Equation (2)
M 1 M 2 f 1 x 1 f 1 x 2 f 2 x 1 f 2 x 2 x 1 x 2 = A 1 J 11 0 0 J 22 A x 1 x 2 = A 1 J 11 y 1 J 22 y 2 , w h e r e : f 1 x 1 f 1 x 2 f 2 x 1 f 2 x 2 = A 1 J 11 0 0 J 22 A , y = y 1 y 2 = A x 1 x 2 .
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 l i n e a r s y s t e m 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
x = A x + B u + ξ ( x , u ) y = C x + D u L a p l a c e y u = G ( s ) C B s I A + D ,
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 gain-scheduling 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.

2.2. 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,
x = f ( x , a ) ,
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
z = i ω z + g ( z , z * ) z C 1 ,
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
ω z * S z ( z , z * ) z S z * ( z , z * ) = f 3 , , n ( S ( z , z * ) , p ( z , z * ) ) ,
and the n-2 solutions are exponentially convergent to S ( t )
x 3 , , n ( t ) = S ( t ) + O ( e β t ) , a 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 , x 2 ) f ( x , α ) R 2 , α R , linearization about the equilibrium gives
x = A ( α ) x + h ( x , α ) , w h e r e : A ( α ) = a b c d ,
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
λ 2 t r ( A ) λ + d e t ( A ) = 0 λ = 1 2 ( t r ( A ) + t r ( A ) 2 4 d e t ( A ) ) λ * = 1 2 ( t r ( A ) t r ( A ) 2 4 d e t ( A ) ) ,
where t r ( A ) and d e t ( A ) are the trace and determinant, respectively. Andronov–Hopf bifurcation takes place at α = 0 satisfying t r ( A ) = 0 and d e t ( A ) = ω 0 2 . Assume in a small neighborhood α ( 0 , δ ) that the eigenvalues are written as λ 1 ( α ) = μ + i ω and λ 2 ( α ) = μ i ω , where μ = t r ( A ) and ω 2 = t r ( A ) 2 d e t ( A ) ; then, T ( α ) exists to transform A ( α ) into canonical real form as follows
J = T ( α ) A ( α ) T 1 ( α ) = μ ( α ) ω ( α ) ω ( α ) μ ( α ) ,
where J 1 = J T / ( μ 2 + ω 2 ) . Let y = T ( α ) x and Equation (9) turns to
y 1 = μ ( α ) y 1 ω ( α ) y 2 + o ( | y | 2 ) y 2 = ω ( α ) y 1 + μ ( α ) y 2 + o ( | y | 2 ) .
Introducing z = y 1 + i y 2 and z * = y 1 i y 2 and Equation (12) presents one-dimensional equivalence as follows
z = λ ( α ) z + g ( z , z * ) ,
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 λ 2 A q * = 1 λ 2 A T p , q * = λ 1 λ 2 p , q * p , q * = 0 .
Equation (11) after transposition has T T ( α ) J T = A T ( α ) T T ( α ) , and if p ( α ) has the following form
p ( α ) = p 1 p 2 = T 11 i T 21 T 12 i T 22 ,
where T i j is i-th line, j-th row element, then z has the explicit form as follows
z = p ( α ) , x .
Considering q , p = 1 and Equation(13), the original variable x has the only equivalent transform
x = q q * ( z z * ) .
Substituting Equation (17) into (16), one obtains
z = λ ( α ) z + p ( α ) , F ( z q + z * q * , α ) .
Here, the nonlinear term F ( x , α ) in Equation (9) after Taylor expansion gives
F ( x , 0 ) = 1 2 B ( x , x ) + 1 6 C ( x , x , x ) + O ( | | x | | 4 ) , w i t h B 1 , 2 ( x , y ) = j , k = 1 2 2 F 1 , 2 ( ξ , 0 ) ξ j ξ k | ξ = 0 x j y k C 1 , 2 ( x , y , m ) = j , k = 1 , l = 1 2 3 F 1 , 2 ( ξ , 0 ) ξ j ξ k ξ m | ξ = 0 x j y k m l .
Substituting Equation (19) into (18) and the second-order term B gives
z = λ z + p , B ( q , q ) 2 z 2 + p , B ( q , q * ) z z * + p , B ( q * , q * ) 2 z * 2 + O ( | z | 3 ) .
Through a coordinate change in Equation (21) the second-order term vanishes and gives Equation (22)
z = ω + p , B ( q , q ) 2 λ ω 2 + p , B ( q , q * ) λ * ω ω * + p , B ( q * , q * ) 2 ( 2 λ * λ ) ω * 2
ω = λ ω + O ( | ω | 3 )
Similarly, the third-order term after coordinate change obtains
ω = λ ω + σ ω 2 ω * + O ( | ω | 4 ) , w h e r e : σ = p , C ( q , q , q * ) 2 .
Next, let ω = ρ e i ϕ and Equation (23) turns to
ρ = ρ ( α ρ 2 ) + Φ ( ρ , φ ) φ = 1 + Ψ ( ρ , φ ) , w h e r e : Φ = O ( | ρ | 4 ) ; Ψ = O ( | ρ | 3 ) .
Suppose that φ performs constant rate rotation and Equation (24) mapped to the rotating plane gives
d ρ d φ = ρ ( α ρ 2 ) + Φ 1 + Ψ = ρ ( α ρ 2 ) + ( ρ , φ ) , w h e r e : R = O ( | ρ | 4 ) .
Assume that ρ 0 = 0 is the initial condition, solving Equation (24) that provides
ρ = e α φ ρ 0 + e α φ 1 e 2 α φ 2 α ρ 0 3 + O 1 ( | ρ 0 | 4 ) .
For a fixed α > 0 , there exists a limit cycle, and when φ = 2 π , the state ρ 1 would turn back to the initial point
ρ 1 = e 2 π α ρ 0 e 2 π α ( 2 π + O ( α ) ) ρ 0 3 + O 2 ( | ρ 0 | 4 ) .
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
x ^ 1 x ^ 2 = α 1 1 α x ^ 1 x ^ 2 ± ( x ^ 1 2 + x ^ 2 2 ) x ^ 1 x ^ 2 .
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 variables x ^ 1 and x ^ 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 either x ^ 1 or x ^ 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.

3. Stabilizing Normal Form Andronov–Hopf Bifurcation

3.1. 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.

3.2. 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
x = A x + B u + h ( x , u ) ,
where the control u = [ u 1 , u 2 ] T takes u = v + k y , 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 d z / d t = y , and the closed-loop system is provided as follows
x = A x + B u + h ( x , u ) z = P ( x z ) u = K ( x z ) .
Stabilizing Equation (30) requires the following augmented matrix A C being Hurwitz
x z = A + B K B K P P x z = A C x z .
 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
T 1 A C T 1 1 = I 0 0 P 1 A + B K B K P P I 0 0 P 1 1 = A + B K K B P I P , T 2 T 1 A C T 1 1 T 2 1 = I M 0 I A + B K B K P P I M 0 I 1 = A + B K + M A M K B M M 2 K B P M P I M P .
Through the following small gain analysis, the control ( P , K ) is obtained
A M + K B M + M 2 + K B P + M P = 0 , s . t . P = ε P 1 + O ( ε 2 ) M = M 0 + ε M 1 + O ( ε 2 ) .
Canceling the first- and second-order perturbational terms provides
( A + B K + M 0 ) M 0 = 0 ( A + B K ) M 1 + A P 1 = 0 , w h e r e : M 0 = A B K ,
which dictates,
M 1 = A P 1 ( A + B K ) 1 .
Then, around M 0 in a small region, the only M satisfying Equation (33) is given
M = M 0 + ε M 1 + O ( ε 2 ) = A B K ε A P 1 ( A + B K ) 1 + O ( ε 2 )
with A C satisfying
T 2 T 1 A C T 1 1 T 2 1 = ε A P 1 ( A + B K ) 1 + O ( ε 2 ) 0 I A C 2 ( 2 , 2 ) , w h e r e : A C 2 ( 2 , 2 ) = A + B K + ε ( A P 1 ( A + B K ) 1 P 1 ) + O ( ε 2 ) .
When ( A , B ) is stabilizable, A + B K can be designed Hurwitz. Since A is non-singular, setting P 1 = A 1 ( A + B K ) , then A P 1 ( A + B K ) 1 is Hurwitz, so as A C 2 . 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.
Proof. 
From Equation (31), one has det( A C ) = det(A)det( P ); hence, sign(det( A C )) = sign(det(A)det( P )). If all eigenvaluses of A C have a negative real part, then sign(det( A C )) = ( 1 ) 2 n , which leads to sign(det(A)det( P )) = ( 1 ) n l ( 1 ) n r , where l represents the number of unstable eigenvalues of A, and r represents the corresponding numbers of P , i.e., l + r is an even number.
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 ) x ( s ) = K s I s I + P .
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
s I s I + ε p * = s Δ s + ε p 22 * ε p 12 * ε p 21 * s + ε p 11 * , w h e r e : Δ = s 2 + ε ( p 11 * + p 22 * ) s + ε 2 ( p 11 * p 22 * p 12 * p 21 * ) .
In accordance with the form of washout filters, the higher-ordered terms in Equation (39) need to be simplified. Considering only the transient dynamics are concerned with washout filters, and based on the L a p l a c e e x t r e m e t h e o r e m , s in Equation (39) takes maximal value, which leads to the idea of the form of washout filters
s I s I + ε P * = s s + ε ( P 11 * + P 22 * ) 0 0 s s + ε ( P 11 * + P 22 * ) .
Meaning two washout filters with d = t r ( P ) aid the controller design, and Equation (40) is decoupled. □

3.3. 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
x = A x + B u + h ( x , u ) z = x i d z = y x = a 1 1 a x 1 x 2 ( x 1 2 + x 2 2 ) x 1 x 2 + K i y y = f i ( x i , K i y ) d y ,
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
A ˜ = a 1 1 a ( B K ) i f i x 1 f i x 2 ( B K ) i d .
Expending Equation (42) along the last column gives the characteristic polynomial P c
P c = d e t ( A ˜ λ I 2 + 1 ) = ( ( B K ) i d λ ) ( λ i λ ) + f i x i n · ( B K ) i ( λ i λ ) = ( ( B K ) i d λ ) ( λ λ n ) ( f i x 1 ) n · ( B K ) i ( λ i λ ) .
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.

4. Case Study

4.1. 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 Z y m o m o n a s m o b i l i s 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
d C s d t = μ · C S · C e Y s x · ( K S + C S ) m S · C x + D · ( C s f C s ) d C X d t = μ m a x · C s · C e K s + C s D C x d C e d t = ( k 1 k 2 C p + k 3 C P 2 ) C s · C e ( K s + C s ) D C e d C p d t = μ m a x · C s · C e Y p x · ( K s + C s ) + m p · C x D C p ,
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 two-parameter 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 S 1 and S 2 . Any designing point at S 1 generates a self-sustained oscillator, while the point at S 2 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.

4.2. 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
y C p = s s + d y = μ m a x · C s · C e Y p x · ( K s + C s ) + m p · C x D C p d y ,
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
D = D 0 + k y ,
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 upper-left 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 washout-filter-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.

5. 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 dual-unstable 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.

Funding

The authors gratefully acknowledge the National Natural Science Foundation of China-Yunnan basic research project (Grant No. 202001AU070048) for their support.

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

Ddilution 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
ttime, h
xbiomass concentration, g/L
Y s x , Y p x yield factor of biomass on substrate/product, g/g
m s maintenance item, g/g h
Greek letters
δ Greek letters
μ m a x maximum specific growth rate, h 1
Subscripts
ssubstrate
xbiomass
eK-compartment
pproduct

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 2 A Θ I n satisfies ( μ i + μ j ) n > i > j > 1 , and the function
φ H ( x , α ) = d e t ( 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 ) ( 2 A Θ I n ) ( P Θ P ) 1 = 2 B Θ I n
where B is the transform of A , B = P A P 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,
( 2 A Θ I n ) ( i , j ) ( k , l ) = α i k α i l i f k = j α i k i f k i a n d l = j α i i + α j j i f k = i a n d l = j α j l i f k = i a n d l j α j k i f l = i 0 e l s e
2 A Θ 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. □

References

  1. Wang, X.; Zhai, C.; Zhao, Z.; Palazoglu, A.; Sun, W.; Ai, J. Analysis of the Onset of Chaos for the Belousov-Zhabotinsky Reaction. Comput. Aided Chem. Eng. 2018, 44, 1375–1380. [Google Scholar]
  2. Pashin, D.; Bastrakova, M.; Satanin, A.; Klenov, N. Bifurcation Oscillator as an Advanced Sensor for Quantum State Control. Sensors 2022, 22, 6580. [Google Scholar] [CrossRef] [PubMed]
  3. Ambrożkiewicz, B.; Syta, A.; Georgiadis, A.; Gassner, A.; Meier, N. Experimental Verification of the Impact of Radial Internal Clearance on a Bearing’s Dynamics. Sensors 2022, 22, 6366. [Google Scholar] [CrossRef] [PubMed]
  4. Zhai, H.; Sands, T. Comparison of Deep Learning and Deterministic Algorithms for Control Modeling. Sensors 2022, 22, 6362. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, C.; Liu, Z.; Liao, M.; Li, P.; Xiao, Q.; Yuan, S. Fractional-order bidirectional associate memory (BAM) neural networks with multiple delays: The case of Hopf bifurcation. Math. Comput. Simul. 2021, 182, 471–494. [Google Scholar] [CrossRef]
  6. Yang, R.; Nie, C.; Jin, D. Spatiotemporal dynamics induced by nonlocal competition in a diffusive predator-prey system with habitat complexity. Nonlinear Dyn. 2022, 110, 879–900. [Google Scholar] [CrossRef]
  7. Krishan, R.; Verma, A. Assessment and Enhancement of Hopf Bifurcation Stability Margin in Uncertain Power Systems. Electr. Power Syst. Res. 2022, 206, 107783. [Google Scholar] [CrossRef]
  8. Alfifi, H.Y. Stability and Hopf bifurcation analysis for the diffusive delay logistic population model with spatially heterogeneous environment. Appl. Math. Comput. 2021, 408, 126362. [Google Scholar] [CrossRef]
  9. Shi, J.; He, K.; Fang, H. Chaos, Hopf bifurcation and control of a fractional-order delay financial system. Math. Comput. Simul. 2022, 194, 348–364. [Google Scholar] [CrossRef]
  10. Fiedler, B.; Flunkert, V.; Georgi, M.; Hövel, P.; Schöll, E. Refuting the odd-number limitation of time-delayed feedback control. Phys. Rev. Lett. 2007, 98, 114101. [Google Scholar] [CrossRef] [Green Version]
  11. Zhai, C.; Sun, W. Analytical approximation of a self-oscillatory reaction system using the Laplace-Borel transform. Chaos Solitons Fractals 2021, 142, 110508. [Google Scholar] [CrossRef]
  12. Cheng, Z. Anti-control of Hopf bifurcation for Chen’s system through washout filters. Neurocomputing 2010, 73, 3139–3146. [Google Scholar] [CrossRef]
  13. Hassouneh, M.A.; Lee, H.C.; Abed, E.H. Washout filters in feedback control: Benefits, limitations and extensions. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004; Volume 5, pp. 3950–3955. [Google Scholar]
  14. Kang, W. Bifurcation control via state feedback for systems with a single uncontrollable mode. SIAM J. Control Optim. 2000, 38, 1428–1452. [Google Scholar] [CrossRef]
  15. Krener, A.J.; Kang, W.; Chang, D.E. Control bifurcations. IEEE Trans. Autom. Control 2004, 49, 1231–1246. [Google Scholar] [CrossRef]
  16. Abed, E.H. A simple proof of stability on the center manifold for Hopf bifurcation. SIAM Rev. 1988, 30, 487–491. [Google Scholar] [CrossRef] [Green Version]
  17. Abed, E.H.; Fu, J.H. Local feedback stabilization and bifurcation control, I. Hopf bifurcation. Syst. Control Lett. 1986, 7, 11–17. [Google Scholar] [CrossRef] [Green Version]
  18. Abed, E.H.; Fu, J.H. Local feedback stabilization and bifurcation control, II. Stationary bifurcation. Syst. Control Lett. 1987, 8, 467–473. [Google Scholar] [CrossRef] [Green Version]
  19. Abed, E.H.; Wang, H.O.; Chen, R.C. Stabilization of period doubling bifurcations and implications for control of chaos. Phys. D Nonlinear Phenom. 1994, 8, 154–164. [Google Scholar] [CrossRef]
  20. Laiou, M.C.; Mönnigmann, M.; Marquardt, W. Stabilization of nonlinear systems by bifurcation placement. IFAC Proc. Vol. 2004, 37, 675–680. [Google Scholar] [CrossRef]
  21. Zhang, N.; Qiu, T.; Chen, B. Bifurcation control and eigenstructure assignment in continuous solution polymerization of vinyl acetate. Chin. J. Chem. Eng. 2015, 23, 1523–1529. [Google Scholar] [CrossRef]
  22. Wang, H.; Zhang, N.; Qiu, T.; Zhao, J.; He, X.; Chen, B. Optimization of a continuous fermentation process producing 1,3-propane diol with Hopf singularity and unstable operating points as constraints. Chem. Eng. Sci. 2014, 116, 116. [Google Scholar] [CrossRef]
  23. Ding, L.; Hou, C. Stabilizing control of Hopf bifurcation in the Hodgkin-Huxley model via washout filter with linear control term. Nonlinear Dyn. 2010, 60, 131–139. [Google Scholar] [CrossRef]
  24. Suo, L.; Ren, J.; Zhao, Z.; Zhai, C. Study on the nonlinear dynamics of the continuous stirred tank reactors. Processes 2020, 8, 1436. [Google Scholar] [CrossRef]
  25. Hurwitz, A. On the conditions under which an equation hasonly roots with negative real parts. Math. Ann 1895, 46, 273–284. [Google Scholar] [CrossRef]
  26. Jobses, I.M.L.; Egberts, G.T.C.; Luyben, K.C.A.M.; Roels, J.A. Fermentation kinetics of Zymomonasmobilis at high ethanol concentrations: Oscillations in continuous cultures. Biotechnol. Bioeng. 1986, 28, 868–877. [Google Scholar] [CrossRef]
  27. Douglas, J.M.; Rippin, D.W.T. Unsteady state processoperation. Chem. Eng. Sci 1966, 21, 305–315. [Google Scholar] [CrossRef]
  28. Bailey, J.E. Periodic operation of chemical reactors: A review. Chem. Eng. Commun. 1974, 1, 111–124. [Google Scholar] [CrossRef]
Figure 1. Phase portrait of (a) linearized system, (b) supercritical; and (c) subcritical Andronov–Hopf bifurcations, where α = 0.1 , and the arrows point out the time evolution direction.
Figure 1. Phase portrait of (a) linearized system, (b) supercritical; and (c) subcritical Andronov–Hopf bifurcations, where α = 0.1 , and the arrows point out the time evolution direction.
Sensors 22 09334 g001
Figure 2. Block diagram for bifurcation control with washout-filter-aided stabilization.
Figure 2. Block diagram for bifurcation control with washout-filter-aided stabilization.
Sensors 22 09334 g002
Figure 3. Analysis on the closed system of normal form Andronov–Hopf bifurcation.
Figure 3. Analysis on the closed system of normal form Andronov–Hopf bifurcation.
Sensors 22 09334 g003
Figure 4. Time series of the closed-loop system, where the washout-filter-aided feedback parameters are d = 0.5 , k = −5.
Figure 4. Time series of the closed-loop system, where the washout-filter-aided feedback parameters are d = 0.5 , k = −5.
Sensors 22 09334 g004
Figure 5. Bifurcation analyses on the ethanol fermentation process: two-parameter continuation diagram (left); and one-parameter continuation diagram (right).
Figure 5. Bifurcation analyses on the ethanol fermentation process: two-parameter continuation diagram (left); and one-parameter continuation diagram (right).
Sensors 22 09334 g005
Figure 6. Two–parameter continuation diagram for the control variables, where GH1 and GH1 are the general Hopf bifurcation, ZH is the zero Hopf bifurcation and HH is the Hopf–Hopf bifurcation.
Figure 6. Two–parameter continuation diagram for the control variables, where GH1 and GH1 are the general Hopf bifurcation, ZH is the zero Hopf bifurcation and HH is the Hopf–Hopf bifurcation.
Sensors 22 09334 g006
Figure 7. Two-parameter continuation diagram for the designing parameters.
Figure 7. Two-parameter continuation diagram for the designing parameters.
Sensors 22 09334 g007
Figure 8. Stabilized outcomes for the closed-loop system, where d = 2 and k = 1 .
Figure 8. Stabilized outcomes for the closed-loop system, where d = 2 and k = 1 .
Sensors 22 09334 g008
Table 1. Parameters used for the continuous bio-ethanol reactor model.
Table 1. Parameters used for the continuous bio-ethanol reactor model.
ParameterDefinitionValue
μ m a x Maximum specific biomass growth rate (h 1 )1.0
Y s x Yield factor for substrate to biomass (g/g)0.00244498
Y p x Yield factor for substrate to product (g/g)0.00526315
m s Maintenance substrate consumption rate (g/g h 1 )2.16
K s Monod constant (g/L)0.5
k 1 Empirical constant (h 1 )16
k 2 Empirical constant (L/g h 1 )0.497
k 3 Empirical constant (L 2 /g 2 h 1 )0.000383
m p Maintenance product consumption rate (g/g h 1 )1.1
C s f Substrate feed concentration (g/L)140
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhai, C.; Yang, C.; Na, J. Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters. Sensors 2022, 22, 9334. https://doi.org/10.3390/s22239334

AMA Style

Zhai C, Yang C, Na J. Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters. Sensors. 2022; 22(23):9334. https://doi.org/10.3390/s22239334

Chicago/Turabian Style

Zhai, Chi, Chunxi Yang, and Jing Na. 2022. "Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters" Sensors 22, no. 23: 9334. https://doi.org/10.3390/s22239334

APA Style

Zhai, C., Yang, C., & Na, J. (2022). Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters. Sensors, 22(23), 9334. https://doi.org/10.3390/s22239334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop