ISRN Applied Mathematics

Volume 2014 (2014), Article ID 168247, 8 pages

http://dx.doi.org/10.1155/2014/168247

Research Article

A New Three-Oscillator Model for the Heart System in the Case of Time Delay and Designing Appropriate Controller for Its Synchronization

Siroos Nazari,¹ Aghileh Heydari,² Mahboobeh Tavakoli,³ and Javad Khaligh¹

¹Department of Mathematics, Payame Noor University, Tehran, Iran

²Department of Mathematics, Payame Noor University, Mashhad, Iran

³Department of Mathematics, Farhangian University, Mashhad, Iran

Received 6 October 2013; Accepted 12 November 2013; Published 5 June 2014

Academic Editors: A. El-Sayed and C.-H. Lien

Copyright © 2014 Siroos Nazari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

If the and oscillators are not synchronized, it may arise some kinds of blocking arrhythmias in the system of heart. In this paper, in order to examine the heart system more precisely, we apply the three-oscillator model of the heart system, and to prevent arrhythmias, perform the following steps. Firstly, we add a voltage with rang and frequency to node. Then, we use delay time factor in oscillators and finally the appropriate control is designed. In this paper, we have explained how simulating and curing these arrhythmias are possible by designing a three-oscillator system for heart in the state of delay and without delay and by applying an appropriate control. In the end, we present the simulation results.

1. Introduction

As we know, If the and oscillators are not synchronized, it may cause some kinds of blocking arrhythmias in the system of the heart. At present, there are devices called pacemaker using electric shocks to synchronize if there is an arrhythmia. [1–3]. The main defect of these devices is the error of diagnosis and not having accurate control. For example, they sometimes apply shocking of 600 to 700 voltages to the heart wrongly even when there is no blocking at all, which is too dangerous. Applied shock to the heart muscle constricts all fibers concurrently and consequently all fibers return to recovery state and hope to return to their normal rhythm after that and it may sometimes occur that the patient involves blocking and the device have not diagnosed it so the patient dies. Two popular types of arrhythmias blocking and ectopic are interesting for researchers [4]. In order to simulate ectopic arrhythmias, a voltage ranging and frequency is added to node.

2. Three-Oscillator Modelling of the Heart without Time Delay

As we know, cardiac normal rhythm is firstly produced by node (normal pacemaker) and then causes stimulating of node. Nevertheless, it has been observed that these two oscillators are not so accurate for producing signal. This is because, the signal of first oscillator is related to the activation of node and the right atrium and signal of second oscillator is only related to the depolarization of the left ventricular. Based on this hypothesis, it is possible to produce curve, but complex may not be produced, because this interval is mainly due to ventricular repolarization [5]. These observations and also the existence of blocking arrhythmias force us to put the third oscillator which represents the spread of the pulse through ventricles, which physiologically represents His-Purkinje complex. To create a general model, we assume that there is a mutual coupling asymmetry among all oscillators [6]. In addition, external stimulation is entered into the system, considering periodic stimulus sentence to each oscillator. This developed model can be shown by a system of differential equations as shown below, in which a Vander pol equation has been considered to model oscillator and another one to model the oscillator and the third one to model oscillator [7]:

in which pairs , , and show , , and oscillators, respectively.

To choose appropriate parameters from the Vander-pol system which are close to the system of the heart, we use trying and error test [8]. Regarding this test, appropriate parameters have been mentioned in Table 1. By choosing parameters as shown in Table 1, any oscillator system is synchronized with another system that has different parameters. In the case that two systems are not synchronized, with the help of changing parameters and coupling coefficients, arrhythmias will occur [9, 10].

Table 1: Appropriate parameters to synchronize three-oscillator system.

Parameter
Definition
Three-oscillators system

SA frequency
1

AV frequency
2

HP frequency
1

SA voltage range
5

AV voltage range
6

HP voltage range
4

SA damping Coefficient
5

AV damping Coefficient
6

HP damping Coefficient
7

Frequency

Coupling coefficients between and

1

1

1

2

5

3

SA membrance flow

SA membrance voltage

AV membrance flow

AV membrance voltage

HP membrance flow

HP membrance voltage

By changing , , and coefficients, we can produce different kinds of alternative responses. These coefficients influence the nonlinear sentence of the equation and cause stability of limit cycle in the phase plane . Having limit cycle with the heart behavior is adapted physiologically to the behavior of the heart [11].

In the discussion of synchronization, synchronizing of the slow oscillators with the fast oscillators is done more easily. In the case of heart, slower pacemakers such as and , should be synchronized with faster pacemakers such as . If it is not synchronized by changing the coupling coefficients, we try to synchronize oscillators by applying appropriate input controller to slow oscillators as shown below

The state variables of slow systems will be converged to state variables of fast system after a transient time. In fact, the second and third oscillators should follow the first oscillator which contains dominant frequency, so synchronization error is considered as follows [12]:

The purpose of the synchronization is to annihilate errors. So, selected control functions to vanish errors as , , are measured as follows:

So, by choosing control rules as follows:

matrix form of error equations is

By applying controllers and , the results of simulations show that diagrams related to synchronization are tracking each other and

Simulation diagrams related to synchronization in the case of three oscillators without time delay have been shown below. See Figures 1, 2, and 3 (the horizontal axis time (s) and the vertical axis millivolts).

(a)

(b)

Figure 1: Synchronization of states , without time delay befor control (a) and after control (b).

(a)

(b)

Figure 2: Synchronization of states , without time delay befor control (a) and after control (b).

(a)

(b)

Figure 3: Error states , and , without time delay after control ((a), (b)).

3. Three-Oscillator Modelling of the Heart with Time Delay

Since even small delays may change the dynamics of the system, including of time delay in differential equations can cause drastic changes and chaos in a system that has regular behavior [7, 13]. Thus, system (1) is changed as below in which represents time delay and :

By substitution, we will have

Now, by applying the appropriate controller on it which is as follows [14]:

the state variables of slow systems will be converged to state variables of fast system, after a transient time. In fact, the second and third oscillators should follow the first oscillator containing dominant frequency, so synchronization error is considered as follows [11]:

The purpose of the synchronization is to annihilate errors. So, chose control functions to vanish errors as , , that are measured as follows:

In this case,

By substitution, we will have

Also,

And control functions are as follows:

Also,

By substitution, we will have

And control functions is as follows:

and matrix form of error equations is

In this case, the third oscillator acts as a pacemaker and its advantage is that instead of a sudden shock and a voltage of , it applies sinusoidal voltage for treatment.

Simulation diagrams related to synchronization in the case of three oscillators with time delay have been shown below. See Figures 4, 5, and 6 (the horizontal axis time (s) and the vertical axis millivolts).

(a)

(b)

Figure 4: Synchronization of states , with time delay befor control (a) and after control (b).

(a)

(b)

Figure 5: Synchronization of states , with time delay befor control (a) and after control (b).

(a)

(b)

Figure 6: Error states , and , with time delay after control ((a), (b)).

4. Conclusion

In this paper, it has been explained that how simulating and curing some arrhythmias are possible, that one example of curing method, is modelling three-oscillator heart system in the state of delay and without delay and by applying an appropriate control.

In this paper after applying control to the model of three-oscillator heart system in case of time delay and without time delay and also observing of simulated diagrams corresponding to synchronization of this model, we observe that, the time of synchronization , and , and also the time when the error converges to zero in , and , after applying designed controller; in case of without time delay is less than 2(s), and in case of time delay less than 1.5(s).

These results indicate that the model of three-oscillator of the heart system in case of time delay and designed controller corresponding to synchronization of this model, is more appropriate and accurate than the state of without time delay.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

A. M. dos Santos, S. R. Lopes, and R. L. Viana, “Rhythm synchronization and chaotic modulation of coupled Van der Pol oscillators in a model for the heartbeat,” Physica A, vol. 338, no. 3-4, pp. 335–355, 2004.

D. Di Bernardo, M. G. Signorini, and S. Cerutti, “A model of two nonlinear coupled oscillators for the study of heartbeat dynamics,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 8, no. 10, pp. 1975–1985, 1998.

B. B. Ferreira, A. S. de Paula, and M. A. Savi, “Chaos control applied to heart rhythm dynamics,” Chaos, Solitons and Fractals, vol. 44, no. 8, pp. 587–599, 2011.

L. Glass, “Synchronization and rhythmic processes in physiology,” Nature, vol. 410, no. 6825, pp. 277–284, 2001.

W. Thanom and R. N. K. Loh, “Observer-based nonlinear feedback controls for heartbeat ECG tracking systems,” Intelligent Control and Automation, vol. 3, pp. 251–261, 2012.

B. van der Pol and J. van der Mark, “The heartbeat considered as a relaxation oscillator and an electrical model of the heart,” Philosophical Magazine and Journal of Science Series, vol. 7, supplement 6, pp. 763–775, 1928.

S. R. F. S. M. Gois and M. A. Savi, “An analysis of heart rhythm dynamics using a three-coupled oscillator model,” Chaos, Solitons and Fractals, vol. 41, no. 5, pp. 2553–2565, 2009.

G. Gottwald and L. Melbourne, “A new test for chaos in deterministic systems,” Nonlinear Dynamics, vol. 69, no. 3, pp. 1255–1262, 2012.

C. Nguyen, Van Der Pol Oscillators Synchronization: Methods and Applications, 2009.

L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990.

S. Sato, S. Doi, and T. Nomura, “Bonhoeffer-van der Pol Oscillator model of the sino-atrial node: a possible mechanism of heart rate regulation,” Methods of Information in Medicine, vol. 33, no. 1, pp. 116–119, 1994.

H. Gholizade-Narm, A. Azemi, M. Khademi, and M. Karimi-Ghartemani, “A state observer and a synchronization method for heart pacemakers,” Journal of Applied Sciences, vol. 8, no. 18, pp. 3175–3182, 2008.

K. Rompala, R. Rand, and H. Howland, “Dynamics of three coupled van der Pol oscillators with application to circadian rhythms,” Communications in Nonlinear Science and Numerical Simulation, vol. 12, no. 5, pp. 794–803, 2007.

W. Thanom and R. N. K. Loh, “Nonlinear control of heartbeat models,” Journal on Systemics, Cybernetics and Informatics, vol. 9, no. 1, pp. 21–27, 2009.