30.06.2013 Views

multiple time scale dynamics with two fast variables and one slow ...

multiple time scale dynamics with two fast variables and one slow ...

multiple time scale dynamics with two fast variables and one slow ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

papers. Geometric singular perturbation theory has been used successfully to<br />

analyze (3.5). In [69] the <strong>fast</strong> pulse is constructed using the exchange lemma<br />

[72, 68, 14]. The exchange lemma has also been used to prove the existence of a<br />

codimension <strong>two</strong> connection between <strong>fast</strong> <strong>and</strong> <strong>slow</strong> waves in (s,ǫ, a)-parameter<br />

space [82]. An extension of Fenichel’s theorem <strong>and</strong> Melnikov’s method can<br />

be employed to prove the existence of heteroclinic connections for parameter<br />

regimes of (3.5) <strong>with</strong> <strong>two</strong> fixed points [103]. The general theory of relaxation<br />

oscillations in <strong>fast</strong>-<strong>slow</strong> systems applies to (3.5) (see e.g. [93, 51]) as does - at<br />

least partially - the theory of canards (see e.g. [104, 35, 37, 86]).<br />

The equations (3.5) have been analyzed numerically by Champneys, Kirk,<br />

Knobloch, Oldeman <strong>and</strong> Sneyd [18] using the numerical bifurcation software<br />

AUTO [32, 33]. They considered the following parameter values:<br />

γ=1, a= 1<br />

, δ=5<br />

10<br />

We shall fix those values to allow comparison of our results <strong>with</strong> theirs. Hence<br />

we also write f1/10(u)= f (u). Changing from the <strong>fast</strong> <strong>time</strong> t to the <strong>slow</strong> <strong>time</strong>τ<br />

<strong>and</strong> relabeling <strong>variables</strong> x1= u, x2= v <strong>and</strong> y=w we get:<br />

ǫ ˙x1 = x2<br />

ǫ ˙x2 = 1<br />

5 (sx2−x1(x1− 1)( 1<br />

10 − x1)+y− p)= 1<br />

5 (sx2− f (x1)+y− p) (3.6)<br />

˙y = 1<br />

s (x1− y)<br />

From now on we refer to (3.6) as “the” FitzHugh-Nagumo equation. Investi-<br />

gating bifurcations in the (p, s) parameter space <strong>one</strong> finds C-shaped curves of<br />

homoclinic orbits <strong>and</strong> a U-shaped curve of Hopf bifurcations; see Figure 3.1.<br />

Only part of the bifurcation diagram is shown in Figure 3.1. There is another<br />

64

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!