Processing math: 100%
A journal of IEEE and CAA , publishes high-quality papers in English on original theoretical/experimental research and development in all areas of automation
Dianwei Qian, Hui Ding, SukGyu Lee and Hyansu Bae, "Suppression of Chaotic Behaviors in a Complex Biological System by Disturbance Observer-based Derivative-Integral Terminal Sliding Mode," IEEE/CAA J. Autom. Sinica, vol. 7, no. 1, pp. 126-135, Jan. 2020. doi: 10.1109/JAS.2019.1911834
Citation: Dianwei Qian, Hui Ding, SukGyu Lee and Hyansu Bae, "Suppression of Chaotic Behaviors in a Complex Biological System by Disturbance Observer-based Derivative-Integral Terminal Sliding Mode," IEEE/CAA J. Autom. Sinica, vol. 7, no. 1, pp. 126-135, Jan. 2020. doi: 10.1109/JAS.2019.1911834

Suppression of Chaotic Behaviors in a Complex Biological System by Disturbance Observer-based Derivative-Integral Terminal Sliding Mode

doi: 10.1109/JAS.2019.1911834
Funds:  This work was supported by the Fundamental Research Funds for the Central Universities (2018MS29)
More Information
  • Coronary artery systems are a kind of complex biological systems. Their chaotic phenomena can lead to serious health problems and illness development. From the perspective of engineering, this paper investigates the chaos suppression problem. At first, nonlinear dynamics of coronary artery systems are presented. To suppress the chaotic phenomena, the method of derivative-integral terminal sliding mode control is adopted. Since coronary artery systems suffer from uncertainties, the technique of disturbance observer is taken into consideration. The stability of such a control system that integrates the derivative-integral terminal sliding mode controller and the disturbance observer is proven in the sense of Lyapunov. To verify the feasibility and effectiveness of the proposed strategy, simulation results are illustrated in comparison with a benchmark.

     

  • BIOLOGICAL systems are diverse and they are different from each other. This paper focuses on a kind of biological systems, called coronary artery systems. The coronary artery systems are very complex and they are directly related with heart health. Coronary artery disease has dominated mortality for most of the past century worldwide [1]. Its treatments, both pharmaceutical and surgical, have become leading sectors of the healthcare economy. In 2015, coronary artery disease affected 110 million people and resulted in 8.9 million deaths [2]. In order to explore new treatment options, both the medical and engineering communities try to provide some insight into coronary artery disease.

    From the aspect of medicine, coronary arteries act as the blood vessels that supply oxygen-rich blood to the heart muscle. In case plaque builds up inside the coronary arteries, the buildup of plaque gradually obstruct the coronary arteries. Such obstruction may lead to vascular spasm. The vasospasms constitute the basic cause of a variety of coronary artery disease. Some drugs have been developed via the mechanism [3]–[7].

    In the field of engineering, a coronary artery system consists of the bio-mathematical model of blood vessels, where the vasospasms are the chaotic states of the blood vessels [8]. Since chaotic systems are sensitive to perturbations, the chaotic phenomena may result in fatal chaos in the coronary arteries [9]. The chaotic phenomena must be suppressed to avoid serious health problems and illness development. In order that treatments can be effective, the chaos suppression has to achieve state synchronization of the blood vessels with pathological changes and the normal blood vessels.

    The problem of state synchronization of a chaotic system was proposed by Pecora and Carroll [10]. Since then, the solution of this problem has become emerging because of the potential applications. Many synchronization approaches have been developed for various systems. The approaches can be roughly cataloged into two classes, that is, linear or partial linearisation methods [11], [12] and nonlinear approaches [13], [14]. Especially, linear or partial linearisation methods are just appropriate for stable linear systems that are characterised by a few distinct peaks in their power spectrum.

    Concerning the coronary artery system, its complex bio-system is inherently nonlinear and calls for nonlinear methods to achieve its state synchronization. The research topic has been paid more and more attention and some reports can be seen. Some synchronization methods have been reported. Gong et al. [15] developed a backstepping-based controller to synchronize the spastic vessel with a normal vessel. Li [9] designed a nonlinear tracking controller and the controller could robustly and adaptively drive a chaotic coronary artery system into the normal orbit. Wu et al. [16] considered the drug absorption time and presented the time-delay state feedback control synthesis.

    As a nonlinear design tool, the methodology of sliding mode control (SMC) is advocated for its invariance, meaning that a SMC system on the sliding-mode stage is completely insensitive to parametric uncertainties and external disturbances under certain matching conditions [17]. The characteristic makes the SMC exhibit better performances than the ordinary conventional robust control methods. In practice, Lin et al. [18] have investigated that the SMC-based method can be carried out for the chaos suppression problem of coronary artery systems. Zhao et al. also have applied the higher-order sliding mode control [19] and the terminal sliding mode control [20] for the same chaos suppression problem.

    One of the SMC design methods is entitled “terminal sliding mode control (T-SMC)”. Such a design method is characterized by its terminal sliding surface. Usually, the design method has faster convergence speed and higher steady-state accuracy than other SMC design methods [21]. As an extension, the integral T-SMC method introduces sign and fractional integral terminal sliding mode concepts and can track relative-degree-one systems with uncertainty and disturbance [22]. The integral T-SMC design method allows a system to start on the integral terminal sliding surface. Further, the derivative-integral T-SMC design method, a progress of the extension, is derived to track high-order systems because it can guarantee the exact estimation of finite error convergence time [23].

    Motivated by the merits of the progress, this paper investigates the derivative-integral T-SMC design method for the chaos suppression problem in a coronary artery system. The purpose of the paper is to reduce the occurrence of coronary artery disease by suppressing the abnormal chaotic behaviors of the coronary artery system. Even when they have different initial conditions, the abnormal chaotic behaviors can still synchronize with the normally unstable periodic orbit. For this purpose, the derivative-integral T-SMC design method is adopted to achieve the behavior synchronization.

    On the other hand, the coronary artery system suffers from external disturbances such as the perturbations of blood pressure and body temperature in reality. The external disturbances challenge the chaos suppression problem very much because they are neither known nor predictive. As proven, the disturbance-observer technique can effectively deal with disturbances with guaranteed robustness [24], [25]. Here the technique is employed in order to attack the challenging issue.

    In the paper, both the derivative-integral T-SMC design method and the disturbance-observer technique are integrated to suppress the abnormal chaotic behaviors in the coronary artery system with disturbances. In the sense of Lyapunov, the stability of the synchronization control system is analyzed. Compared with some published results, the feasibility and effectiveness of the integration solution is illustrated via some numerical results.

    The remainder of this paper is organized as follows. Section II introduces the bio-mathematical model of the coronary artery system. Section III formulates the control problem. The synchronization control via the derivative-integral T-SMC-based controller and the disturbance observer is derived in Section IV for the coronary artery system. To show the synchronization performance, numerical simulations are carried out in Sections V. Finally, some conclusions are drawn in Section VI.

    Inherently, coronary arteries are a kind of muscular blood vessels. Concerning a coronary artery system, the lumped-parameter model that describes its dynamics [18] has the form of

    ˙x1=bx1cx2˙x2=λ(1+b)x1λ(1+c)x2+λx31+Ecosωt (1)

    where x1 is the inradius change of the vessel, x2 means the pressure change in the vessel, t indicates the time variable, b, c and λ denote lumped parameters of the coronary artery system and Ecosωt represents a periodical disturbance term.

    As one factor that results in myocardial infarction, coronary artery spasm is caused by partial or complete occlusion of coronary arteries [1]. Equation (1) reveals the medical phenomenon from the perspective of mathematics. In (1), the initial condition [x1(0) x2(0)]T=[0.2 0.2]T is taken into consideration, where the parameters E and ω are determined by 0.3 and 1, respectively.

    Fig. 1 reveals the the bifurcation diagram of the coronary artery system (1) with respect to the change of λ, where another parameters are set by b=0.15, c=1.7, E=0.3 and ω=1. Fig. 2 reveals the the bifurcation diagram of the coronary artery system (1) with respect to the change of b, where another parameters are set by λ=0.65, c=1.7, E=0.3 and ω=1. Fig. 3 reveals the the bifurcation diagram of the coronary artery system (1) with respect to the change of c, where another parameters are set by λ=0.65, b=0.15, E=0.3 and ω=1. All the bifurcation diagrams share the same initial condition [x1(0),x2(0)]T=[0.2, 0.2]T. The time series are set by 0, 2π, 4π, 6π, 10π and 12π, where the system period is 2π.

    Figure  1.  Bifurcation diagrams with respect to the change of λ from –1 to 0.
    Figure  2.  Bifurcation diagrams with respect to the change of b from 0 to 1.
    Figure  3.  Bifurcation diagrams with respect to the change of c from –2 to –1.5.

    The bifurcation diagrams in Figs. 13 illustrates that the coronary artery system under some parameters is disordered and will lead to chaos. As far as the coronary artery system is concerned, the following lumped parameters are considered, that is, λ=0.65, b=0.15, c=1.7, E=0.3 and ω=1. With the group of lumped parameters, Fig. 4 illustrates the coronary artery system has very complex dynamics in the phase plane.

    Figure  4.  Phase plane of the coronary artery system.

    From Figs. 14, the coronary artery system is apparently disordered with respect to the changes of these lumped parameters. The coronary artery system may descend into chaos under some initial conditions with respect to the change of the lumped parameters. The chaotic behaviors in clinic medicine indicate that the vascular spasm of the coronary artery system may exhibit a series of coronary artery disease, including but not limited to angina, myocardial infarction, and sudden cardiac death. The chaos extremely hazard health and it must be suppressed immediately. Since the chaos suppression problem of such a coronary artery system is challenging and interesting, this paper explores the problem from the aspect of control design via derivative-integral terminal sliding mode.

    From the bio-mathematical model (1), the coronary artery system with uncertainties (2) can be drawn

    ˙x1=bx1cx2+d1˙x2=λ(1+b)x1λ(1+c)x2+λx31+Ecosωt+u(t)+d2. (2)

    Here x=[x1, x2]T is defined as the state vector and u(t) generated by the designed controller is the control input. The vector d=[d1, d2] describes the uncertainties, meaning unmodelled dynamics, external disturbances and structural variations.

    Assumption 1: d is bounded, that is, ||d||d, where d is positive but unknown.

    Assumption 2: d is slowly time varying, that is, ˙d02×1, where 02×1=[0, 0]T.

    From (2), the nominal system of the coronary artery system with no control input can be written by

    ˙ˉx1=bˉx1cˉx2˙ˉx2=λ(1+b)ˉx1λ(1+c)ˉx2+λˉx31+Ecosωt. (3)

    Concerning the bio-mathematical model (1), (3) can display different dynamics, highly sensitive to its parameters. This fact means that even small differences of the parameters yield widely diverging outcomes of the chaotic behaviors. Consequently, the purpose of chaos suppression in the coronary artery system is to synchronize the uncertain coronary artery system (2) with the nominal system (3) via derivative-integral terminal sliding mode, that is

    limte=02×1. (4)

    Here e=[e1, e2], e1=x1ˉx1 and e2=x2ˉx2.

    Considering (2), (3) and (4), the error system can be described by

    ˙e1=be1ce2+d1˙e2=λ(1+b)e1λ(1+c)e2+λe31+3x1ˉx1e1+u+d2. (5)

    Further, define ˜x1=e1 and ˜x2=˙e1 and differentiate ˜x1 and ˜x2 with respect to the time variable t. Then, substituting (5) into ˙˜x2 has the form of

    ˙˜x1=˜x2˙˜x2=λ(cb)˜x1(b+λ+cλ)˜x2λc˜x313cx1ˉx1˜x1cucd2+˙d1+λ(1+c)d1. (6)

    According to Assumption 2, ˙d10 in (6). Then, (6) can be written by (7) in the form of state space.

    ˙˜x=F(˜x,d)+Buy=H(˜x). (7)

    Here ˜x=[˜x1, ˜x2]T, d=cd2+λ(1+c)d1, H(˜x)=˜x1, B = [0, c]T, F(˜x,d)=[˜x2, f(˜x)+d]T and f(˜x)=λ(cb)˜x1 (b+λ+cλ)˜x2λc˜x313cx1ˉx1˜x1.

    Assumption 3: The functions F(˜x,d) and H(˜x) are sufficiently smooth so that all Lie derivative calculations are well-defined.

    Proven by Chiu [23], the derivative-integral T-SMC design can be available for high relative-degree systems that have higher relative degree than 1. According to Assumption 3, the relative degree of (7) can be calculated as (8), (9) and (10) in order to develop the derivative-integral T-SMC design for (7).

    LBH(˜x,d)=H(˜x,d)˜xB=[10]×[0c]=0 (8)
    LFH(˜x,d)=H(˜x,d)˜xF=[10]×[˜x2f(˜x)+d]=˜x2 (9)
    LBLFH(˜x,d)=LB(LFH(˜x,d))=(LFH(˜x,d))˜xB=[01]×[0c]=c0. (10)

    From (8), (9) and (10), the relative degree of (7) can be determined as 2 so that it is possible to develop a derivative-integral terminal sliding mode controller for the chaos suppression [23]. Furthermore, the input-output dynamics can be expressed as

    ¨y=L2FH(˜x,d)+LBLFH(˜x,d)u. (11)

    Here L2FH(˜x,d)=(LFH(˜x,d))˜xF=f(˜x)+d. Then, (11) can be written by

    ¨y=f(˜x)+dcu. (12)

    This paper investigates the derivative-integral T-SMC method for the chaos suppression problem. The so-called derivative-integral terminal sliding mode means that the derivative and integral terms exist in a sliding surface. Motivated by this purpose, the following sliding surface s is taken into consideration.

    s=eD1+αeI1. (13)

    Here α>0 is a design parameter, eD1=˙ep11/q11D0+βeD0, eD0=e1, eI1=t0eq21/p21D1(τ)dτ and eI1(0)=eD1(0)α, where β>0 is pre-defined and p11>q11>0 chosen by designers are odd integers, as well as p21>q21>0.

    Theorem 1: If the sliding surface is defined as (13), then [˜x1, ˜x2]T in (6) will reach 02×1 in the finite convergence time τDI, where τDI is formulated by

    τDI=|eD1(0)|1q21p21α(1q21p21)+|eD0(t11)|1q11p11β(1q11p11). (14)

    Here t11 is the reaching time of the sliding mode eD1=0.

    Proof: From (13), the sliding mode starts at t=0. From then on, eD1=αeI1 can always hold true via control design. Subsequently, substituting eD1=αeI1 into ˙eI1=eq21/p21D1(t) yields

    ˙eI1=αq21p21eI1q21p21 (15)

    where eI1(0)=eD1(0)/α.

    Solving (15), the convergent time of eI1 can be gotten as

    t11=|eD1(0)|1q21p21α(1q21p21). (16)

    On the sliding surface s=0, eD1=αeI1 indicates that the sliding mode of eD1 takes place at the same time t11 to zero. When eD1=0, eD0 will successively converge to zero. At t=t11, eD0 can be formulated by

    ˙eD0(t11)=βq11p11eq11p11D0(t11). (17)

    From (17), the time spent from eD0(t11) to eD0=0 can be calculated as

    t01=|eD0(t11)|1q11p11β(1q11p11). (18)

    According to (16) and (18), the time spent from s(0)=0 to e1=0 is their summation because the integral and derivative terms in (13) are independent. Consequently, the finite convergence time can have the form of (14), that is, τDI=t11+t01. ■

    In order to obtain the derivative-integral T-SMC law, a Lyapunov candidate is selected as V0=12s2. Differentiate V0 with respect to the time variable t. ˙V0=s˙s can be gotten. Further, the derivative of s can be formulated by

    ˙s=p11q11˙e(p11q111)D0¨eD0+β˙eD0+α(˙ep11q11D0+βeD0)q21p21. (19)

    Substituting (12) into (19) yields

    ˙s=p11q11˙e(p11q111)D0[f(˜x)+dcu]+ψ (20)

    where ψ=β˙eD0+α(˙ep11q11D0+βeD0)q21p21.

    In the sense of Lyapunov, a derivative-integral T-SMC law can be deduced from ˙V0<0, determined by

    u=1c[ke|ψ|(q11p11)A1s|s|f(˜x)[κsgn(s)+ηs]d0] (21)

    where A={˙e(p11q111)D01,|˙e(p11q111)D01|εε,otherwise, κ>0 and η>0 are chosen by designers, d0=[|c|+|λ(1+c)|]d and ke>0 is also a design parameter, determined later.

    Concerning the expression of A, a parameter ϕ>0 exists such that

    |(A˙e(p11q111)D01)A1|ϕ<1. (22)

    In order to have ˙V0<0, ke>11ϕ is picked up.

    Replace u by (21) in (20) and consider the condition (22). ˙V0<0 can be guaranteed if κ>d0 holds true. Unfortunately, d in Assumption 1 is unknown so that d0 is also unknown, meaning that designers can not guarantee ˙V0<0 by selecting an appropriate value of κ.

    From another conservative perspective, the value of κ can be assigned as large as possible in order to have the guaranteed stability. However, such a solution is not prosperous because a large value of κ may induce the chattering phenomenon and hazard the control performance.

    To deal with issue raised by the unknown d0, a disturbance observer is adopted in order that it can estimate d0. Meanwhile, the control system can also possess the guaranteed stability as well as the good control performance. The disturbance observer is formulated by

    ˙z=LB2zL(B2L˜x+A˜x+B1u+M)ˆd=z+L˜x (23)

    Here zR1×1 is the internal state variable of the observer, ˆdR1×1 is the the estimate of d, LR1×2 pre-defined by designers is the observer gain vector, the matrix A is written by [01λ(cb)(b+λ+cλ)], the vectors B1, B2 and M are determined by [0c], [01] and [0λc˜x313cx1ˉx1˜x1].

    Define an estimate error as

    ed=dˆd. (24)

    Assumption 4: ed is bounded, that is, |ed|ed, where ed>0 is unknown.

    The derivative of the estimate error can be written by

    ˙ed=˙d˙ˆd. (25)

    From Assumption 2, the assumption indicates that the uncertainties in (2) are slowly time varying compared with the dynamic characteristics of the observer (23). Further, it is concluded that ˙d0 from Assumption 2. Finally, (26) can be obtained by substituting (23) into (25).

    ˙edLB2P+L(B2L˜x+A˜x+B1u+M)L(A˜x+B1u+B2d+M)=LB2(ˆdL˜x)+L2B2˜xLB2d=LB2(ˆdd)=LB2ed. (26)

    The solution of (26) can be calculated by

    ed=exp(LB2)ed(0). (27)

    Select a vector L such that the eigenvalue of LB2 is a positive constant. This selection means that the estimate error ed is able to exponentially converge as t. Since ed is exponentially convergent from (27), Assumption 4 is mild. Consequently, the unknown d0 in (21) can be replaced by the disturbance observer output ˆd, that is, the estimate of d.

    Theorem 2: Consider the bio-mathematical model of the coronary artery system (2); take Assumptions 1–4 into account; define the error system (5) and have the input-output dynamics (12); formulate the derivative-integral terminal sliding surface (13); design the disturbance observer (23). If the derivative-integral T-SMC law is formulated by (28), the closed-loop chaos suppression system is then asymptotically stable in the presence of uncertainties.

    u=1c[ke|ψ|(q11p11)A1s|s|f(˜x)[κsgn(s)+ηs]ˆd] (28)

    where κ>ed.

    Proof: The aforementioned Lyapunov candidate function V0=12s2 are taken into consideration again. Differentiate V0 with respect to the time variable t and replace ˙s by (20). Then, we can obtain

    ˙V0=sp11q11˙e(p11q111)D0[f(˜x)+dcu]+sψ. (29)

    Substituting (28) into (29) yields

    ˙V0=sψs˙e(p11q111)D0ke|ψ|A1s|s|+sp11q11˙e(p11q111)D0[dˆd(κsgn(s)+ηs)]. (30)

    Owing to p11>0, q11>0 and p11>q11, p11q11>1 exists such that ˙e(p11q111)D01>0 holds true in (30) for all ˙e(p11q111)D010. Consequently, the first and second terms in (30) can have the form of

    sψs˙e(p11q111)D0ke|ψ|A1s|s|=sψke|ψ|s+ke|ψ|AA1ss˙e(p11q111)D0ke|ψ|A1s|s||ψs|ke|ψs|+ke|ψ(A˙e(p11q111)D0)A1s||ψs|ke(1ϕ)|ψs|. (31)

    In (31), ke>(1ϕ) can be picked up in advance such that sψs˙e(p11q111)D0keψA1ss<0.

    Further, the third term in (30) can be written by

    sp11q11˙e(p11q111)D0[dˆd(κsgn(s)+ηs)]=p11q11˙e(p11q111)D0[edsκ|s|ηs2]p11q11˙e(p11q111)D0[(edκ)|s|ηs2]. (32)

    In (32), κ>ed can be selected by designers such that sp11q11˙e(p11q111)D0[dˆd(κsgn(s)+ηs)]<0.

    From (31) and (32), ˙V0<0 holds true by the derivative-integral T-SMC design, that is, the control design can asymptotically stabilize the uncertain coronary artery system by the chaos suppression. ■

    From Assumption 4, ed is unknown as well, indicating that κ in (28) has to be assigned in advance. From the perspective of the system stability, the value of κ should be large enough. In this sense, (28) makes no difference over (21) at a first glance. However, ed is able to exponentially converge via the disturbance observer (23). It is concluded that its maximum exists at t=0, i.e., ed(0). Correspondingly, d is just bounded according to Assumption 1. It is completely unknown when, where and how many its maximum is. Consequently, κ in (21), associated with d0, should be more conservative in order to guarantee the stability, while κ in (28), concerned to ed, can contribute to the alleviation of the chattering phenomenon.

    This section presents some numerical simulation results that illustrate the feasibility and validity of the proposed strategy in comparison with a benchmark. Concerning the uncertain system (2) and the nominal system (3), the lumped parameters are fixed to λ=0.65, b=0.15, c=1.7, E=0.3 and ω=1, which are kept unchanged from the model presented by Zhao et al. [19], [20]. The initial states are assigned to [ˉx1(0), ˉx2(0)]T=[0.2, 0.2]T, which are also the same as the initial states presented in [19], [20]. Considering the derivative-integral T-SMC law, the controller parameters are determined by α=2, β=2, p11=p21=11, q11=q21=13, ke=2, η=0.1 and κ=2. As far as the disturbance observer is concerned, its sole parameter L, the gain vector, is given by [0, 5]T.

    The uncertain system (2) contains two disturbance sources, that is, d1 and d2. From the viewpoint of sliding mode control, the two disturbances have different characteristics that d2 is matched and d1 is mismatched. Further, d1 and d2 are merged into d in (6). From the expression of d, both u and d2 share the same coefficient c. This fact indicates that d2 is matched. Meanwhile, d1 can hardly be matched because its coefficient λ(1+c) is almost impossible to be c.

    Pointed out by Utkin [17], the matched uncertainties enter the coronary artery system through the control channel u, so that the control system is insensitive to the uncertainties when the sliding mode is reached. Fig. 5 displays the proposed strategy can effectively suppress the chaotic phenomena under the effects of matched uncertainties, where d2 is set by 2sin(t) and d1 is assigned to 0.

    Figure  5.  Error states and control input in (7) in comparison with the benchmark under the effects of matched uncertainties.

    The sliding surface variable is shown in Fig. 6. From Fig. 6, the sliding mode of the proposed strategy can be reached at the limited time as proven in Theorem 1. Thereafter, the coronary artery system can suppress the uncertainties in light of the invariance of SMC [17]. Consequently, the chaos suppression by the proposed strategy in Fig. 5 has the best performance, compared with the benchmark and the derivative-integral T-SMC without observer.

    Figure  6.  Sliding surface variable in comparison with the benchmark under the effects of matched uncertainties.

    Fig. 7 depicts the phase plane trajectories of the error system (7), the nominal system (3) and the uncertain system (2). As proven in Theorem 2, the uncertain coronary artery system with matched uncertainties can synchronize the nominal coronary artery system by the proposed strategy in Fig. 7, while the tracking errors can be convergent to zero in the limited time given in (14).

    Figure  7.  Phase plane trajectories of the error system (7), the nominal system (3) and the uncertain system (2) under the effects of matched uncertainties.

    Further, the estimate error of the disturbance observer is given in Fig. 8. In Fig. 8, the estimate error is very large at the outset. But it can be exponentially convergent till the observer output can track the disturbances. The exponential convergence of ed makes the control performance of the proposed strategy superior to that of the sole derivative-integral T-SMC without observer.

    Figure  8.  Estimate error of the disturbance observer under the effects of matched uncertainties.

    Since the mismatched uncertainties can not be suppressed by the invariance of SMC, it is indispensable to illustrate their effects on the control performance. Here the controller parameters are kept unchanged, as well as the observer gain and the initial conditions. The sole change is that d1 is set by 2sin(t) and d2 is assigned to 0.

    Fig. 9 shows that the proposed strategy can also suppress the chaotic phenomena even if the mismatched uncertainties exist in the coronary artery system. Although the invariance of SMC makes no difference to the mismatched uncertainties, the proposed strategy can deal with them via the integration of the derivative-integral terminal sliding mode controller and the disturbance observer.

    Figure  9.  Error states and control input in (7) with mismatched uncertainties under the effects of mismatched uncertainties.

    The sliding surface variable is demonstrated in Fig. 10. Compared with Fig. 6 and Fig. 10, the sliding mode of the proposed strategy can be reached at the shortest time and keep on sliding on the surface even if the mismatched uncertainties have adverse effects.

    Figure  10.  Sliding surface variable under the effects of mismatched uncertainties.

    Fig. 11 shows the phase plane trajectories of the error system (7), the nominal system (3) and the uncertain system (2) under the effects of mismatched uncertainties. Shown in Fig. 11, the coronary artery system with mismatched uncertainties can achieve the synchronization with the nominal coronary artery system as proven in Theorem 2.

    Figure  11.  Phase plane trajectories of the error system (7), the nominal system (3) and the uncertain system (2) under the effects of mismatched uncertainties.

    The estimate error of the disturbance observer is given in Fig. 12. Similar to the curve in Fig. 8, the estimate error is very large at the outset but it exponentially converge as t.

    Figure  12.  Estimate error of the disturbance observer under the effects of mismatched uncertainties.

    No matter what kinds of the uncertainties are, the proposed strategy in Figs 512 can solve the adverse effects of the uncertainties and synchronize the uncertain coronary artery system with the nominal one. In clinical medicine, the solution can contribute to the reduction and elimination of chaotic vasospasms, where the mechanism is to synchronize any chaotic blood vessels with a nominal one. The proposed strategy can also benefit the research and development of clinical pharmacy for therapeutic purposes because medications can be treated as a kind of control input to suppress the chaotic vasospasms.

    In the field of medicine, coronary artery disease in the early stage is often treated with pharmacies, meaning that a suitable pharmacodynamics has to be drawn for patients. Concerning the victims of advanced coronary artery disease, they have to suffer a lot from invasive procedures such as stent and bypass. These invasive procedures may result in extra pains and risks for patients.

    The proposed strategy can benefit the research and development of clinical pharmacy and medical device for therapeutic purposes because therapies can be treated as a kind of control input to suppress the chaotic vasospasms. Considering a certain therapy, it can be treated as the control input of this chaotic system. Once the therapy is effective, the chaotic phenomenon has to achieve state synchronization of the blood vessels with pathological changes and the normal blood vessels. In this sense, the designed control can be utilized by a pharmacodynamics or a medical device, which is the potential application of the presented control design.

    Coronary artery disease can be induced by blood vessel spasms. Since the vasospasms can be considered as chaos in the aspect of engineering, this paper has explored the chaos suppression problem of coronary artery systems. In order to have some insight into the bio-mathematical model of coronary artery systems, their nonlinear behaviors are presented by a series of bifurcation diagrams with respect to the changes of lumped parameters. The degree of complexity is described by the phase plane trajectory of coronary artery systems. The purpose is to synchronize any chaotic coronary artery with a nominal one. Motivated by the purpose, the method of derivative-integral terminal sliding mode control is taken into consideration. Since coronary artery systems suffer from uncertainties, the uncertainties trouble the control design. To attack this issue, the technique of disturbance observer is employed. The proposed strategy integrates the derivative-integral terminal sliding mode controller and the disturbance observer to achieve the chaos suppression. In the sense of Lyapunov, the stability of such an integrated control system has been proven. The simulation results have demonstrated the effectiveness and feasibility of the proposed strategy via some comparisons. The proposed strategy is beneficial to the research and development of clinical pharmacy.

  • [1]
    G. K. Hansson, " Mechanisms of disease-Inflammation, atherosclerosis, and coronary artery disease,” N. Engl. J. Med., vol. 352, no. 16, pp. 1685–1695, 2005. doi: 10.1056/NEJMra043430
    [2]
    K. Ozaki and T. Tanaka, " Molecular genetics of coronary artery disease,” J. Hum. Genet., vol. 61, no. 1, pp. 71–77, 2016. doi: 10.1038/jhg.2015.70
    [3]
    G. W. He and D. P. Taggart, " Spasm in arterial grafts in coronary artery bypass grafting surgery,” Ann. Thorac. Surg., vol. 101, no. 3, pp. 1222–1229, 2016. doi: 10.1016/j.athoracsur.2015.09.071
    [4]
    X. L. Liu, F Hou, H Qin, and A. M. Hao, " Robust optimization-based coronary artery labeling from X-ray angiograms,” IEEE J. Biomed. Health Inform., vol. 20, no. 6, pp. 1608–1620, 2016. doi: 10.1109/JBHI.2015.2485227
    [5]
    A. Hernandez-Vela, C. Gatta, S. Escalera, L. Igual, V. Martin-Yuste, M. Sabate, and P. Radeva, " Accurate coronary centerline extraction, caliber estimation, and catheter detection in angiographies,” IEEE T. Inf. Technol. Biomed., vol. 16, no. 6, pp. 1332–1340, 2012. doi: 10.1109/TITB.2012.2220781
    [6]
    Z. Qian, I. Marvasty, S. Rinehart, and S. Voros, " A Lesion-specific coronary artery calcium quantification framework for the prediction of cardiac events,” IEEE T. Inf. Technol. Biomed., vol. 15, no. 5, pp. 673–680, 2011. doi: 10.1109/TITB.2011.2162074
    [7]
    C. C. Wang and H. T. Yauuu, " Chaos analysis and synchronization control of coronary artery systems,” Abstract Appl. Anal., pp. 209718, 2013. doi: 10.1155/2013/209718,2013
    [8]
    T. Schauer, N. O. Negard, F. Previdi, K. J. Hunt, and M. H. Fraser, " Ferchland E and Raisch J online identification and nonlinear control of the electrically stimulated quadriceps muscle,” Control Eng. Practice, vol. 13, no. 9, pp. 1207–1219, 2005. doi: 10.1016/j.conengprac.2004.10.006
    [9]
    W. L. Li, " Tracking control of chaotic coronary artery system,” Int. J. Syst. Sci., vol. 43, no. 1, pp. 21–30, 2012. doi: 10.1080/00207721003764125
    [10]
    L. M. Pecora and T. L. Carroll, " Synchronization in chaotic systems,” Phys. Rev. Lett., vol. 64, no. 8, pp. 821–824, 1990. doi: 10.1103/PhysRevLett.64.821
    [11]
    M. Rafikov and J. M. Balthazar, " On control and synchronization in chaotic and hyperchaotic systems via linear feedback control,” Commun. Nonlinear Sci. Numer. Simul., vol. 13, no. 7, pp. 1246–1255, 2008. doi: 10.1016/j.cnsns.2006.12.011
    [12]
    Y. G. Yu, H. X. Li, and J. Duan, " Chaos synchronization of a unified chaotic system via partial linearization,” Chaos Solitons Fractals, vol. 41, no. 1, pp. 457–463, 2009. doi: 10.1016/j.chaos.2008.02.010
    [13]
    D. Ghosh and A. R. Chowdhury, " Nonlinear observer-based impulsive synchronization in chaotic systems with multiple attractors,” Nonlinear Dyn., vol. 60, no. 4, pp. 607–613, 2010. doi: 10.1007/s11071-009-9618-5
    [14]
    A. Chithra and I. R. Mohamed, " Synchronization and chaotic communication in nonlinear circuits with nonlinear coupling,” J. Comput. Electron., vol. 16, no. 3, pp. 833–844, 2017. doi: 10.1007/s10825-017-1013-8
    [15]
    C. Gong, Y. Li, and X. Sun, " Backstepping control of synchronization for biomathematical model of muscular blood vessel,” J. Appl. Sci., vol. 24, no. 6, pp. 604–607, 2006.
    [16]
    W. S. Wu, Z. S. Zhao, J. Zhang, and L. K. Sun, " State feedback synchronization control of coronary artery chaos system with interval timevarying delay,” Nonlinear Dyn., vol. 87, no. 3, pp. 1773–1783, 2017. doi: 10.1007/s11071-016-3151-0
    [17]
    D. W. Qian, C. D. Li, S. G. Lee, and C. Ma, " Robust formation maneuvers through sliding mode for multi-agent systems with uncertainties,” IEEE/CAA J. Autom. Sinica, vol. 5, no. 1, pp. 342–351, 2018.
    [18]
    C. J. Lin, S. K. Yang, and H.T. Yau, " Chaos suppression control of a coronary artery system with uncertainties by using variable structure control,” Comput. Math. Appl., vol. 64, no. 5, pp. 988–995, 2012. doi: 10.1016/j.camwa.2012.03.007
    [19]
    Z. S. Zhao, J. Zhang, G. Ding, and D. K. Zhang, " Chaos synchronization of coronary artery system based on higher order sliding mode adaptive control,” Acta Phys. Sin., vol. 64, no. 21, 2015. doi: 10.7498/aps.64.210508,2015
    [20]
    Z. S. Zhao, X. M. Li, J. Zhang, and Y. Z. Pei, " Terminal sliding mode control with self-tuning for coronary artery system synchronization,” Int. J. Biomath., vol. 10, no. 3, 2017. doi: 10.1142/S1793524517500413,2017
    [21]
    A. M. Zou, K. D. Kumar, Z. G. Hou, and X. Liu, " Finite-time attitude tracking control for spacecraft using terminal sliding mode and Chebyshev neural network,” IEEE Trans. Syst. Man Cybern,Part B-Cybern., vol. 41, no. 4, pp. 950–963, 2011. doi: 10.1109/TSMCB.2010.2101592
    [22]
    C. S. Chiu and C. T. Shen, " Finite-time control of DC-DC buck converters via integral terminal sliding modes,” Int. J. Electron., vol. 99, no. 5, pp. 643–655, 2012. doi: 10.1080/00207217.2011.643493
    [23]
    C. S. Chiu, " Derivative and integral terminal sliding mode control for a class of MIMO nonlinear systems,” Automatica, vol. 48, no. 2, pp. 316–326, 2012. doi: 10.1016/j.automatica.2011.08.055
    [24]
    M. Chen and S. Z. Ge, " Direct adaptive neural control for a class of uncertain nonaffine nonlinear systems based on disturbance observer,” IEEE T. Cybern., vol. 43, no. 4, pp. 1213–1225, 2013. doi: 10.1109/TSMCB.2012.2226577
    [25]
    T. Poloni, I. Kolmanovsky, and B. Rohal’-Ilkiv, " Simple input disturbance observer-based control: case atudies,” J. Dyn. Syst. Meas. Control Trans. ASME, vol. 140, no. 1, 2018. doi: 10.1115/1.4037298,2018
  • Related Articles

    [1]Xiang Chen, Yujuan Wang, Yongduan Song. Unifying Fixed Time and Prescribed Time Control for Strict-Feedback Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2025, 12(2): 347-355. doi: 10.1109/JAS.2024.124401
    [2]Kamal Rsetam, Zhenwei Cao, Zhihong Man, Xian-Ming Zhang. GPIO-Based Continuous Sliding Mode Control for Networked Control Systems Under Communication Delays With Experiments on Servo Motors[J]. IEEE/CAA Journal of Automatica Sinica, 2025, 12(1): 99-113. doi: 10.1109/JAS.2024.124812
    [3]Qinchen Yang, Fukai Zhang, Cong Wang. Deterministic Learning-Based Neural PID Control for Nonlinear Robotic Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(5): 1227-1238. doi: 10.1109/JAS.2024.124224
    [4]Jin-Xi Zhang, Kai-Di Xu, Qing-Guo Wang. Prescribed Performance Tracking Control of Time-Delay Nonlinear Systems With Output Constraints[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(7): 1557-1565. doi: 10.1109/JAS.2023.123831
    [5]Mahdi Taheri, Khashayar Khorasani, Nader Meskin. On Zero Dynamics and Controllable Cyber-Attacks in Cyber-Physical Systems and Dynamic Coding Schemes as Their Countermeasures[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(11): 2191-2203. doi: 10.1109/JAS.2024.124692
    [6]Haijing Wang, Jinzhu Peng, Fangfang Zhang, Yaonan Wang. High-Order Control Barrier Function-Based Safety Control of Constrained Robotic Systems: An Augmented Dynamics Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(12): 2487-2496. doi: 10.1109/JAS.2024.124524
    [7]Zheng Wu, Yiyun Zhao, Fanbiao Li, Tao Yang, Yang Shi, Weihua Gui. Asynchronous Learning-Based Output Feedback Sliding Mode Control for Semi-Markov Jump Systems: A Descriptor Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(6): 1358-1369. doi: 10.1109/JAS.2024.124416
    [8]Wei Chen, Qinglei Hu. Sliding-Mode-Based Attitude Tracking Control of Spacecraft Under Reaction Wheel Uncertainties[J]. IEEE/CAA Journal of Automatica Sinica, 2023, 10(6): 1475-1487. doi: 10.1109/JAS.2022.105665
    [9]Zhiwen Deng, Chang Xu. Frequency Regulation of Power Systems With a Wind Farm by Sliding-Mode-Based Design[J]. IEEE/CAA Journal of Automatica Sinica, 2022, 9(11): 1980-1989. doi: 10.1109/JAS.2022.105407
    [10]Lei Liu, Lifeng Ma, Jie Zhang, Yuming Bo. Sliding Mode Control for Nonlinear Markovian Jump Systems Under Denial-of-Service Attacks[J]. IEEE/CAA Journal of Automatica Sinica, 2020, 7(6): 1638-1648. doi: 10.1109/JAS.2019.1911531
    [11]Shuping He, Jun Song. Finite-time Sliding Mode Control Design for a Class of Uncertain Conic Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2017, 4(4): 809-816. doi: 10.1109/JAS.2017.7510643
    [12]Longsheng Chen, Qi Wang. Adaptive Robust Control for a Class of Uncertain MIMO Non-affine Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(1): 105-112.
    [13]Shanwei Su. Output-feedback Dynamic Surface Control for a Class of Nonlinear Non-minimum Phase Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(1): 96-104.
    [14]Yanlong Zhou, Mou Chen, Changsheng Jiang. Robust Tracking Control of Uncertain MIMO Nonlinear Systems with Application to UAVs[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 25-32.
    [15]Xiaoli Liu, Shengchao Zhen, Kang Huang, Han Zhao, Ye-Hwa Chen, Ke Shao. A Systematic Approach for Designing Analytical Dynamics and Servo Control of Constrained Mechanical Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(4): 382-393.
    [16]Zixuan Liang, Zhang Ren, Xingyue Shao. Decoupling Trajectory Tracking for Gliding Reentry Vehicles[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 115-120.
    [17]Yuanqing Xia, Ning Zhou, Kunfeng Lu, Yong Li. Attitude Control of Multiple Rigid Bodies with Uncertainties and Disturbances[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 2-10.
    [18]Tengfei Liu, Zhongping Jiang. Distributed Control of Nonlinear Uncertain Systems: A Cyclic-small-gain Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(1): 46-53.
    [19]Aleksandra Faust, Peter Ruymgaart, Molly Salman, Rafael Fierro, Lydia Tapia. Continuous Action Reinforcement Learning for Control-Affine Systems with Unknown Dynamics[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(3): 323-336.
    [20]Huiyang Liu, Long Cheng, Min Tan, Zengguang Hou. Containment Control of General Linear Multi-agent Systems with Multiple Dynamic Leaders: a Fast Sliding Mode Based Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(2): 134-140.

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(12)

    Article Metrics

    Article views (1361) PDF downloads(79) Cited by()

    Highlights

    • This paper investigates the chaos suppression problem of a biological coronary artery system.
    • A scheme integrates the sliding mode controller and the disturbance observer.
    • The system stability is presented in the sense of Lyapunov.
    • Results are illustrated to support the scheme in comparison with a benchmark.

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return