12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

518 Solutions to Exercises, Chap. 1019. Let χ = ϕ 2 − ϕ 1 ; then we have˙χ = ω − c s<strong>in</strong> χ .If |ω/c| ≤ 1, then there are two synchronized states, χ = arcs<strong>in</strong> (ω/c) and χ = π − arcs<strong>in</strong> (ω/c),one stable and the other unstable.20. From the theorem by Hoppensteadt and Izhikevich (1997) presented <strong>in</strong> Sect. 10.3.3 it followsthat Kuramoto’s model is a gradient system when c ij = c ji and ψ ij = −ψ ji . From Ermentrout’stheorem presented <strong>in</strong> the same section, it follows that the synchronized state ϕ i = ϕ j is stableif, e.g., all ψ ij = 0 and c ij > 0.21. S<strong>in</strong>ce the probability density function g(ω) is symmetric, the averaged frequency deviation ofthe network is zero, and rotat<strong>in</strong>g the coord<strong>in</strong>ate system, we can make the cluster phase ψ = 0.The network is split <strong>in</strong>to two populations: One oscillat<strong>in</strong>f with the cluster (|ω| < Kr), therebyform<strong>in</strong>g the cluster, and one drift<strong>in</strong>g <strong>in</strong> and out of the cluster. The latter does not contributeto the Kuramoto synchronization <strong>in</strong>dex, because contributions from different oscillators canceleach other on average. In the limit n → ∞, the sum (10.21) becomes the <strong>in</strong>tegral∫∫r = e iϕ(ω) g(ω)d ω ≈ e iϕ(ω) g(ω)d ω .|ω| 0 .(x)25. (Brown et al. 2004) The solution of x ′ = λx with x(0) = x 0 is x(t) = x 0 e λt . The period T =log(∆/x 0 )/λ is found from the condition x(T ) = ∆. Hence, Q(ϑ) = 1/(λx(ϑ)) = 1/(λx 0 e λϑ ) =e λ(T −ϑ) /(∆λ).26. Let us first consider the SNIC case ẋ = 1 + x 2 . The solution star<strong>in</strong>g with x(0) = x reset hasthe form (check by differentiat<strong>in</strong>g) x(t) = tan(t + t 0 ), where t 0 = atan x reset . The periodshould be found from the condition tan(T + t 0 ) = +∞, and it is T = π/2 − t 0 . Hence, x(t) =tan(t+π/2−T ) = − cot(t−T ). Now, Q(ϑ) = 1/(1+x(ϑ) 2 ) = 1/(1+cot 2 (ϑ−T )) = s<strong>in</strong> 2 (ϑ−T ).The homocl<strong>in</strong>ic case ẋ = −1 + x 2 is quite similar. The solution star<strong>in</strong>g with x(0) = x reset > 1has the form (check by differentiat<strong>in</strong>g) x(t) = − coth(t + t 0 ), where t 0 = acoth(−x reset ) = −acoth(x reset ). The period is found from the condition − coth(T +t 0 ) = +∞ result<strong>in</strong>g <strong>in</strong> T = −t 0 .Hence, x(t) = − coth(t − T ). F<strong>in</strong>ally, Q(ϑ) = 1/(−1 + x(ϑ) 2 ) = 1/(1 + coth 2 (ϑ − T )) =s<strong>in</strong>h 2 (ϑ − T ).

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

Saved successfully!

Ooh no, something went wrong!