Sliding-mode-based chaos suppression of coronary artery systems

Since coronary artery disease is the leading global cause of mortality and morbidity, this paper investigates the chaos suppression of coronary artery systems. The motivation of the paper is to discuss and analyze coronary artery disease in the field of dynamics. Firstly, the mathematic model of coronary artery systems is formulated and the properties of this model are illustrated by bifurcation diagram, information entropy analysis, phase plane trajectory, and Poincaré section. With regard to the uncertainties of coronary artery systems, the disturbance observer technique is adopted. Meanwhile, the smooth second-order sliding mode controller is designed to suppress the chaos phenomenon. In light of the combination of the controller and observer, the stability of such a closed-loop system is proven in the sense of Lyapunov. Finally, some numerical simulations demonstrate the feasibility and validity of such design.


Introduction
Coronary artery systems are two blood vessels that branch from the aorta.They are close to the point of departure from the heart.They take charge of carrying oxygen-rich blood to the myocardium.Blockage of any branch of the coronary arteries will cause death of a portion of the heart tissue because it is deprived of oxygen-rich blood.Such blockage leads from coronary artery disease. 1 It is reported that the coronary artery disease is the leading global cause of mortality and morbidity. 2rom the aspect of medical treatment, the coronary artery disease is characterized by an inadequate supply of oxygen-rich blood to the myocardium because of narrowing or blocking of a coronary artery by fatty plaques.If the oxygen depletion is extreme, the effect may be a myocardial infarction, that is, heart attack.If the deprivation is insufficient to cause the death of a section of heart muscle, the effect may be angina pectoris, that is, pain or discomfort in the chest.Both the effects can be fatal because they can lead to heart failure, ventricular fibrillation as far as sudden death. 3herefore, effective management and treatment of the coronary artery disease rely on risk assessment stratify.In order to assess the risk, mathematical modeling is powerful. 4hen heart failure or ventricular fibrillation takes place, it is characterized by an uncontrolled and uncoordinated contraction of the coronary artery systems.
From the aspect of mathematical modeling, such contraction can be treated as the phenomenon of chaos. 5ith regard to the chaotic phenomena, such abnormal and fatal chaos can induce the coronary artery disease because chaotic systems are very sensitive to perturbations.Unfortunately, there are just many perturbations in the coronary artery systems, that is, body temperature, blood pressure, blood viscosity, and so forth.This fact indicates that the effective management and treatment can be mathematically considered as the chaos suppression problem of the coronary artery systems. 6ince the coronary artery systems are inherently nonlinear, they are rather complicated. 7Such characteristics challenge the chaos suppression very much.So far, many methods have been presented to suppress the chaos of coronary artery systems, including sliding mode or variable structure control methods, 8 timedelay based control design, 9,10 adaptive fuzzy control approach, 11 to name but a few.The idea behind these methods is to suppress the fatal chaos so as to avoid serious health problems and illness development. 12For the purpose of effective treatment, the motion states in a vessel with pathological changes have to synchronize with the motion states of a normal vessel through the designed chaos-suppression methods.Therefore, two problems from the aspect of chaos suppression rise up for coronary artery systems.One is how to gain some insight into the chaotic behaviors.The other is how to develop the effective chaos-suppression methods.
As a powerful design tool, the methodology of sliding mode control (SMC) is advocated and investigated for the chaos-suppression problem of the coronary artery systems because the SMC methods exhibit better control performances than the conventional robust control methods. 13,14Although the SMC methods suffer from the chattering drawback, many sliding mode-based methods have been reported for the chaossuppression problem, that is, terminal sliding mode control, 15 higher order sliding mode adaptive control 16 and fractional-order sliding mode control. 17There is a variety of SMC design methods, wherein one is called second-order SMC.][20] Concerning the coronary artery systems, the achievement of the chaos-suppression problem indicates that the symptoms of the coronary artery disease are disappeared or relieved, which will be beneficial for patients a lot.Further, doctors argue that the chattering should be as less as possible during the process of chaos suppression, which will inspire and encourage us to develop more effective SMC methods with high performance.Motivated by the merits of the second-order SMC method, this paper investigates a smooth second-order SMC method for the chaos-suppression problem of the coronary artery systems.Once the abnormal and fatal chaotic behaviors take place in the coronary artery systems, the smooth second-order SMC method is able to reject them despite uncertainties.Meanwhile, the kind of SMC method is able to avoid the chattering.
Clinically, many risk factors are imperceptible to patients but pivotal to the coronary artery systems, such as, the perturbation of blood pressure, the concentration of serum cholesterol, the change of ambient temperature, even the fluctuation of emotion. 21These risk factors can hardly be measurable, predictive or known in advance.They can be treated as the disturbances of the coronary artery systems, which make the chaos suppression problem extremely tough.It is reported that the disturbance-observer technique is powerful to resist disturbances. 22,23This paper fuses the technique with the smooth second-order SMC method for this chaos-suppression problem.
Medically, the developed second-order SMC method means a certain clinical pharmacy and medical device for therapeutic purposes.The proposed strategy can promote the research and development because therapies can be treated as a kind of control input to suppress the chaotic vasospasm.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 also the potential application of the presented control design.
This paper faces with the chaos-suppression problem of the coronary artery systems.The remainder of this paper is organized as follows.Modeling and chaotic behavior analysis of the coronary artery systems are addressed in Section 2. Section 3 designs the smooth second-order SMC method in the sense of Lyapunov.In Section 4, the control design is implemented via a numerical platform and some comparisons are illustrated to support the control design.Finally, Section 5 concludes this paper.

Bio-mathematical model
Coronary arteries extend over the surface of heart and branch into smaller capillaries.They move in a complex pattern during each cardiac cycle.Characteristically, such arteries are categorized as muscular blood vessels.The lumped-parameter model that is applicable to the interpretation of a coronary artery has the form of Here x 1 means the change of inside radius in the vessel, x 2 indicates the pressure change in the vessel, t is the time variable, b, c, and l denote the lumped parameters in this model, and E cos vt represents a periodical disturbance term.

Model-based analysis
No matter what coronary artery disease is, one of the clinic symptoms is vasospasm.Since the coronary artery spasm means an uncontrolled and uncoordinated contraction of the blood vessels that deliver blood to heart muscle, many serious consequences may be believed to be due to partial or complete occlusion of the blood vessels.( 1) is just a window that gains some insight into such spasm from the aspect of dynamics.
In (1), define the initial condition at t = 0 as ½x 1 (0) x 2 (0) T = ½0:2 0:2 T , consider the parameters in the disturbance term as E = 0:3 and v = 1, and set the time series as 0, 2p, 4p, 6p, 8p, 10p, and 12p where the period is 2p.Now, take the changes of the parameters l, b, and c into account.In order to visualize the effects of varying a certain parameter at a time and keeping others fixed, the bifurcation diagrams of the system (1) are employed.

Bifurcation diagram
The bifurcation diagrams of the coronary artery system are illustrated in Figure 1.In Figure 1 1 is equally divided into 4999 portions.Then, 5000 endpoints can be obtained in each closed interval.On each endpoint, x 1 and x 2 can be calculated according to the given time series.Finally, the curves can be plotted.

Information entropy analysis
According to its concept, information entropy refers to the degree of chaos in a system.As an indicator, the more orderly the system is, the lower value its information entropy has.On the contrary, the more chaotic the system is, the higher value its information entropy has.
Here the data in Figure 1 are employed.For each parameter, 5000 numbers in its closed interval are available.In order to calculate the information entropy of the coronary artery system, each 10 consecutive numbers in the closed interval are divided into a sub-interval so that 500 sub-intervals are gotten.Apparently, the maximum and minimum can be found in the 5000 numbers.
Define another closed interval that is between the extreme values.Then, divide this interval into 10 equal sub-intervals.Concerning each 10 consecutive numbers, count the quantity that they fall into each sub-interval.Finally, the value of information entropy H can be calculated by Here p i is the probability that the 10 consecutive numbers fall into the ith sub-interval, formulated by the quantity that the 10 consecutive numbers falls into the ith sub-interval over 10. Figure 2 illustrates the diagrams of information entropy of the coronary artery system.It can be calculated that the maximum of information entropy is log 2

Dynamic behaviors
Consider the coronary artery system described by (1)  with the selected representative parameters.Figure 3 illustrates the dynamic behaviors of the coronary artery system.Figure 3(a) in the phase plane demonstrates that the trajectory of this coronary artery system is very complex.It is well known that the technique of Poincare´section is a general method to explain complex dynamics because the technique can help to visualize the complexity.Here the Poincare´section analysis is carried out as insight into the complexity of the coronary artery system.Since the system (1) possesses two variables, an extended system should be defined in order to depict the system dynamics in 3 dimension.For this purpose, define x 3 = E cos vt and substitute x 3 into (1).Then, let x 3 = 0 as the section.Finally, the Poincare´section in three dimension is displayed in Figure 3(b), where the lime-green plane is the defined section x 3 = 0, the blue curve is the trajectory of the extended system in three dimension, and the red curve on the lime-green plane is just the Poincare´section.Figure 3(c) shows the Poincare´section in two dimension.Apparently, no matter what analysis tools are adopted, the dynamic behaviors of the coronary artery system revealed by such tools are very complex.In clinic, these dynamic behaviors referring to the vasospasm of the coronary artery system implicate a large number of coronary artery diseases.Unfortunately, such complexity challenges the chaos suppression problem of the coronary artery system very much.

Control design
According to the bio-mathematical model (1), a coronary artery system with uncertainties can have the form of where u(t) means the control input injected into the coronary artery system and it is generated by the designed controller.
Both Assumption 1 and Assumption 2 originate from the biological properties of the coronary artery system because violent disturbances may cause some irreversible damages of the coronary artery system.From the perspective, both the assumptions are mild.
The nominal system (4) with no control input can be drawn in light of the uncertain model (3).
x 2 T is the state vector of the nominal system.Apparently, (4) seems similar to (1) very much.In fact, they depict the coronary artery system from the two aspects.(1) formulates the bio-mathematical model and ( 4) is derived from the coronary artery system with uncertainties (3).Therefore, the different variables are adopted in ( 1) and ( 4).
According to the aforementioned model-based analysis about (1), not only the nominal system (4) exhibits the complex dynamic behaviors, but also it is sensitive to the system parameters.For the chaos suppression problem of the uncertain coronary artery system (3), the developed smooth second-order sliding mode controller carried out in (3) needs to possess the ability of forcing the error vector e in (5) to 0 231 once the uncertain coronary artery system becomes disordered.Here e = ½e 1 , e 2 T , e 1 = x 1 À x 1 and e 2 = x 2 À x 2 .Substituting (3) and ( 4) into the definition of e yields Let x1 = e 1 and x2 = _ e 1 .Then, differentiate x1 and x2 with respect to t. Substituting (6) into x : 2 yields Considering Assumption 2, there exists _ d 1 F0 in (7).Then, ( 7) can be written by (8) in the form of state space. Here Design of disturbance observer (9) indicates that d is also bounded concerning Assumption 1.Without loss of generality, the boundary of d is unknown such that d Ã 0 .0 is unknown.Unfortunately, the value of d Ã 0 is affiliated with the system stability and the control performance very much.In order to solve this issue concerning the unknown d Ã 0 , the disturbance-observer technique is adopted and a disturbance observer is designed to estimate d.Such a disturbance observer can be described by Here z 2 < 131 is the internal state variable of this observer, d 2 < 131 is the the estimate of d, L 2 < 132 is the observer gain vector and L is pre-defined by designers, the matrix A is written by 0 In (11), the derivative of the estimate error can be written by Compared with the dynamic characteristics of the observer (10), Assumption 2 indicates that the uncertainties in ( 2) are slowly time varying such that _ dF0 in (12).Finally, ( 13) can be obtained by substituting (10)  into (12).
The solution of ( 13) can be calculated by In ( 14), the vector L can be picked up by designers such that the eigenvalue of LB 2 is a positive constant.Such design means that the estimate error e d is able to exponentially converge as t !' if e d (0) is not infinite.Apparently, the condition that e d (0) is finite is not strict, that is, Assumption 3 is quit mild.Consequently, the output of this disturbance observer d can be treated as the estimate of d and d will be adopted to replace the unknown d in (8) in the following controller design.Inevitably, such replacement may affect the stability of the coronary artery system, which should be discussed in theory.

Controller design
In order to synchronize (8) by the smooth second-order SMC method, a sliding surface needs to be defined firstly.
Here the predefined constant b is positive.Accordingly, the derivative of s has the form of Substituting ( 8) into (16) yields In (17), the disturbance term d is unknown such that the output of the disturbance observer in ( 10) is employed to replace the term.By such replacement, ( 17) can be re-written by Since the coronary artery system (8) is biological, not only may the dramatic change badly damage the coronary artery system, but also such change may cause the coronary artery diseases.This indicates that the control input u in ( 8) should be as smooth as possible.For the purpose, the smooth second-order SMC law is defined as Here k 1 and k 2 are positive and they are predefined, sgn( Á ) is the sign function.

Stability analysis
In (19), the output of the disturbance observer is involved in the smooth second-order sliding mode controller to achieve the chaos suppression of this coronary artery system.Due to the combination of the controller and the observer, the stability of this chaos suppression problem should be taken into consideration.Replacing u in ( 8) by (19) gives Then, the stability of the coronary artery system (8) by the smooth second-order sliding mode controller (19)  with the disturbance observer ( 10) is equivalent to the stability of (20).
Theorem 1: Take Assumptions 1, 2, and 3 into account; consider the bio-mathematical model of the coronary artery system (8); design the disturbance observer (10); formulate the sliding surface (15); the smooth secondorder sliding mode controller has the form of (19).If k 1 .0 is as large as possible and k 2 .0 in (19), then the closed-loop system is asymptotically stable in the presence of uncertainties.
Proof: Select the following Lyapunov candidate V as The derivative of V with respect to time t in (21) has the form of Substituting ( 20) into ( 22) yields Now, the inequation k 1 jsj5e Ã d s will be discussed in the following five cases.is also apparently held true such that dV dt = 0. Case 3: For 0 \ s \ 1, there exists jsj \ s 1 3 .In order to have dV dt 40, not only the value of k 1 is larger than e Ã d , but also k 1 has to be as large as possible to assure k 1 jsj5e Ã d s Case 4: For s = 1, there exists jsj = s 1 3 .Similarly, the value of k 1 needs to be larger than e Ã d to assure k 1 .e Ã d in order to have dV dt \ 0. Case 5: For s .1, there exists jsj .s In (23), dV dt 40 will be held true if k 1 .0, k 2 .0, and , that is, the closed-loop control system is of asymptotical stability in the sense of Lyapunov.Ä Remark: From the aforementioned discussions, for any s, dV dt 40 can be held true if k 1 is large enough and k 1 .e Ã d .However, e Ã d .0 is unknown according to Assumption 3, which means that the value of k 1 does not have any reference at all.Fortunately, (14) indicates that e d is exponentially convergent to zero as t !' and that the maximum of e d in one dynamic process is e d (0), which means that e Ã d must be larger than or equal to e d (0).Namely, the value of k 1 must not be smaller than e d (0) from the aspect of the system stability.

Simulation results
For the purpose of verifying the feasibility, the designed control method in this section will be carried out by a benchmark.Take the bio-mathematical model of a coronary artery system with uncertainties into consideration, where the system parameters in (3) are determined by l = À 0:65, b = 0:15, c = À 1:7, E = 0:3, and v = 1.These parameters are kept unchanged from the model presented by Zhao et al. 10,11,16 The initial states of the system in (3) are depicted by ½ x 1 (0) x 2 (0) T = ½ 0:2 0:2 T .Concerning the disturbance observer, the observer gain vector is assigned to L = ½ 5 5 T .With regard to the smooth second-order sliding mode controller, its parameter in the sliding surface ( 15) is picked up as b = 70 and the parameters in the control law (19) are given by k 1 = 150 and k 2 = 100.
In (3), both the equations contain uncertainties.But their properties of the uncertainties are inherently different.Apparently, the term d 2 enters the coronary artery system via the control channel such that it belongs to the matched uncertainties.Pointed out by Utkin, 14 the invariance means a sliding mode control system is insensitive to such matched uncertainties when the sliding mode is reached.On the other hand, the mismatched term d 1 is completely in the coronary artery system with no direct effects by the control channel, which is just indirectly effected by the control channel according to the dynamics (3).
In order to illustrate the validity of the designed control method, the mismatched term is set to d 1 = 2 sin t due to its challenges of.Accordingly, the matched term is set to d 2 = 0 due to the invariance.Figure 4 demonstrates the phase plane trajectories of the uncertain system (3) under the effects of mismatched term, the error system (7) and the nominal system (4).From Figure 4, the uncertain coronary artery system (3) in Figure 4(a) can synchronize with the nominal system (4) in Figure 4(c) by the designed control method as proven in Theorem 1.Meanwhile, the tracking errors of the control system in Figure 4(b) can be asymptotically convergent to zero. Figure 5 displays that the designed control method can suppress the chaotic phenomena and overcome the mismatched uncertainties in the coronary artery system.The error of the disturbance observer in (11) is illustrated in Figure 6.The sliding surface in (15) and the control input in (19) are shown in Figure 7.
The results in Figures 4 to 7 indicate that the designed control method possesses the ability of dealing with the adverse effects of the uncertainties in the coronary artery system.The method is able to synchronize the uncertain coronary artery system with the nominal one.In practice, therapy can be considered as a control input for any disease.Such method is potential to develop a machine or medicine for therapeutic purposes in order to synchronize chaotic blood vessels with a nominal one.

Conclusions
Abnormal and fatal chaos of coronary artery systems is able to induce the coronary artery disease.This paper has addressed the chaos-suppression problem by the sliding-mode-based design.The mathematic model of the coronary artery systems is described by dynamics.Some insight into the model has been illustrated by bifurcation diagram, information entropy analysis, phase plane trajectory, and Poincare´section.Due to the uncertainties the coronary artery systems tolerate, the disturbance-observer technique is employed to estimate the uncertainties.The smooth second-order sliding mode controller is also explored to achieve the chaos-suppression problem with the assistance of this observer.The stability of the closed-loop control system with the controller and the observer has been proven via the direct Lyapunov method.Finally, some numerical results have been displayed to verify the feasibility and validity of the presented control design.Concerning the therapeutic purposes, such design is potential and instructive to develop a clinical pharmacy or to invent a medical device.Patients are also able to benefit from such design a lot to avoid the risks of coronary artery disease.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

ORCID iD
Dianwei Qian https://orcid.org/0000-0001-7153-7997 (a), the diagram shows the changes of x 1 and x 2 with respect to the change of l.Here l is considered in the closed interval ½À1, 0 and the other two parameters are fixed as b = 0:15 and c = À 1:7.Similarly, the diagram in Figure 1(b) reveals the changes of x 1 and x 2 with respect to the change of b.Here b is considered in the closed interval ½0, 1 and the other two parameters are fixed as l = À 0:65 and c = À 1:7.The diagram in Figure 1(c) displays the changes of x 1 and x 2 with respect to the change of c.Here c is considered in the closed interval ½À2, À 1:5 and the other two parameters are fixed as l = À 0:65 and b = 0:15.The closed interval of each parameter in Figure

Figure 1 .
Figure 1.Bifurcation diagrams: (a) obtained by the change of l from À1 to 0, (b) obtained by the change of b from 0 to 1, and (c) obtained by the change of c from À2 to À1.

Figure 2 .
Figure 2. Diagrams of information entropy: (a) calculated by the change of l from À1 to 0, (b) calculated by the change of b from 0 to 1, and (c) calculated by the change of c from À2 to À1.

Figure 3 .
Figure 3. Diagrams of dynamic behaviors: (a) phase plane trajectory, (b) Poincare ´section in three dimension, and (c) Poincare ´section.

.
Define the estimate error of this disturbance observer e d as e d = d À d ð11Þ Assumption 3: e d is bounded, that is, je d j4e Ã d , where e Ã d .0 is unknown.

1 3 .
In order to have dV dt \ 0, the value of k 1 should be picked up to guarantee k 1 .e Ã d such that k 1 jsj .e Ã

Figure 4 .
Figure 4. Phase plane trajectories of the uncertain system (3) under the effects of mismatched uncertainties, the error system (7) and the nominal system (4).(a) Uncertain system, (b) error system, and (c) nominal system.

Figure 5 .
Figure 5. Error states in (7) by the designed control method.

Figure 7 .
Figure 7. Sliding surface and control input by the designed control method: (a) sliding surface and (b) control input.
1and d 2 indicate the uncertain factors that may affect the inside radius and pressure in the vessel, that is, body temperature, blood pressure, blood viscosity, and environmental temperature.In (3), define x = ½x 1 , x 2 T as the state vector and d = ½d 1 , d 2 T as the disturbance vector.With regard to d, the following two assumptions are considered.
Assumption 1: d is bounded, depicted by jjdjj ' 4d Ã where d Ã is a positive constant.Assumption 2: d is slowly time varying, formulated by _ d ' 0 231 where