Firefly Synchronization of ad- hoc networks Dr. Michael Emmerich - - PowerPoint PPT Presentation

Firefly Synchronization of ad- hoc networks

• Dr. Michael Emmerich

Natural Computing Group

Synchronizing using pulse-coupled

• scillators: Fireflies, Heart, and

Wireless Networks

Simple Systems, Complex Behavior FIREFLIES SYNCHRONIZATION

Fireflies: mysterious mass synchrony

 Φιρεβυγ σ αρε βεετλεσ κνοων φορ τηειρ χονσπιχυουσ υσε οφ

βιολυµινεσχενχε το αττραχτ µατεσ ορ πρεψ

 Τηαιλανδ, ωιτη φανταστιχ, ουτ οφ τηισ ωορλδ φιρεφλψ σηοωσ;

ενορµουσ χονγ ρεγ ατιονσ οφ φιρεφλιεσ βλινκινγ ον ανδ οφφ ιν υνισον, ιν δισπλαψσ τηατ συπποσεδλψ στρετχηεδ φορ µιλεσ αλονγ τηε ριϖερβανκσ. (Αλσο οχχυρινγ ιν Αφριχα, ανδ σοµε µορε πλαχεσ )

 ηττπ://ωωω.ψουτυβε.χοµ/ωατχη?ϖ= α−

ςψ7ΝΖΤΓοσ

 Αχχουντσ ον τηισ πηενοµενον βψ Ωεστερν τραϖελερσ το Σουτη Εαστ

Ασια γ ο βαχκ ασ φαρ ασ 300 ψεαρσ.

 Μψστεριουσ φορµ οφ µασσ σψνχηρονψ.  Ιν 1917 Πηιλιπ Λαυρεντ ωροτε υπ αν εξπλανατιον ιν Σχιενχε: τηε

“ αππαρεντ πηενοµενον ωασ χαυσεδ βψ τηε τωιστινγ ορ συδδεν λοωερινγ ανδ ραισινγ οφ µψ εψλιδσ τηε ινσεχτσ ηαδ νοτηινγ το δο ωιτη ιτ”

Early hypothesis …

 Τηε φιρεφλιεσ ηαϖε α χεντραλ χοορδινατορ ορ χονδυχτορ  Χρυχιαλ εξπεριµεντ ιν βεδ δισµισσεδ τηισ ηψποτηεσισ: Τηε

‘ ’ βιολογ ιστ χουπλε Βυχκ ανδ Βυχκ τοοκ αρβιτραρψ φιρεβυγ σ ιντο τηειρ βεδροοµ ατ νιγ ητ ανδ τηεψ σποντανεουσλψ σψνχηρονιζεδ ωηεν τηεψ ωερε πυτ το τηε χειλινγ , ωιτηουτ εξτερναλ φορχε.

 Πεσκιν ανδ οτηερσ φουνδ τηατ α σιµπλε υνιϖερσαλ µεχηανισµ

χαν βε υσεδ το εξπλαιν δεχεντραλιζεδ σψνχηρονιζατιον; Στρογ ατζ προϖιδεδ α προοφ τηατ α ποπυλατιον χαν σψνχηρονιζε

Pacemaker of the Heart

Χη. Πεσκιν αλσο προποσεδ α σχηεµατιχ µοδελ φορ ηοω τηε παχεµακερ χελλσ οφ τηε ηεαρτ σψνχηρονιζε τηεµσελϖεσ



Παχεµακερ οφ τηε ηεαρτ



µοστ ιµπρεσσιϖε οσχιλλατορ εϖερ χρεατεδ



α χλυστερ οφ 10,000 χελλσ χαλλεδ σινοατριαλ νοδε



γ ενερατεσ ελεχτριχαλ ρηψτηµ τηατ χοµµανδσ τηε ρεστ οφ τηε ηεαρτ το βεατ



ηασ το βε δονε ρελιαβλψ, µινυτε αφτερ µινυτε



τηρεε βιλλιον βεατσ ιν α λιφετιµε



υνλικε µοστ χελλσ ιν ηεαρτ, τηε παχεµακερ χελλσ οσχιλλατε αυτοµατιχαλλψ: ισολατεδ



ιν πετρι διση, τηειρ ϖολταγ ε ρισεσ ανδ φαλλσ ινρεγ υλαρ ρηψτηµ



Αλλ οφ ωηιχη ραισεσ τηε θυεστιον: Ωηψ δο ωε νεεδ σο µανψ χελλσ, ιφ ονε χαν δο τηε ϕοβ?



προβαβλψ βεχαυσε α χεντραλιζεδ χοντρολλερ ισ νοτ ροβυστ δεσιγ ν: α χεντραλ χοντρολλερ χαν µαλφυνχτιον ορ διε ανδ τηισ ωιλλ δεσταβιλιζε εντιρε σψστεµ

Charles S. Peskin http://www.math.nyu.edu/faculty/peskin/ also site for the book: “Modeling and simulation in the life sciences”

Peskin model (1975) of the heart pacemaker

 Τηε φυνχτιον φ ισ χαλλεδ πηασε ρεσπονσε χυρϖε (ΠΡΧ)

αν ελεχτριχαλ – ϖολταγ ε. Φορ φ ι(τ) < φ τη ιτ οβεψσ τηε φολλοωινγ διφφερεντιαλ λαω: δφ ι/δτ = (φ ι(τ+δτ)− φ ι(τ)) / δτ = 1/Τ (χονσταντλψ ινχρεασινγ ποτεντιαλ)

 Ωηεν τηισ πηασε αρριϖεσ ατ σοµε τιµε τ το α τηρεσηολδ ϖαλυε φι(τ)= φ τη,

τηε πηασε ισ ρεσετ το 0 ανδ τηε πηασε οφ τηε νειγ ηβορινγ σιτεσ ισ µοδιφιεδ βψ αν οφφσετ φ(φκ(τ)): φ ι (τ) = = > 0 φ κ (τ+δτ) = = > φκ (τ) + φ(φκ) φορ αλλ κ: κ≠ ι φ(φ) = (α− 1)φ +β, α>1, β>0

 Πυλσε−

χουπλεδ, ωηεν οσχιλλατορσ αρριϖεσ το α ϖαλυε φ τη σψνχηρονουσλψ

 Ιφ τηε χουπλινγ ισ ποσιτιϖε (εξχιτατορψ) τηε ποπυλατιον τενδσ το

σψνχηρονιζε, ι.ε., το αρριϖε το τηε τηρεσηολδ ατ τηε σαµε τιµε.

http://hermes.ffn.ub.es/~albert/peskinen.html http://math.nyu.edu/faculty/peskin/heartnotes/

φi(t)=t/T+c T

Uniform oscillator, phase portrait

Remark: Phase portrait resembles a ‘clock’. Can be modeled as the complex function (x(t),y(t)) = exp(i 2 π t/T ), where i denotes the imaginary number i=√(-1) and φ= 2 π (t MODULO 1)

Phase response diagram

Oscillator receives firing signal from other

• scillator

Coupled oscillator

 Τωο νοδεσ ωιτη πηασε φυνχτιονσ φ1 ανδ φ2  φ0 ισ τηε σιγ ναλ οφ νοδε 2 ϕυστ αφτερ ιτ ηασ

ρεχειϖεδ τηε πυλσε φροµ νοδε 1 ανδ υπδατεδ.

Convergence analysis (1)

 Ασσυµε φτη= 1, ανδ τωο νοδεσ ωιτη πηασε φυνχτιον φ 1 ανδ φ2  φ 0 ισ τηε ποσιτιον οφ νοδε 2 αφτερ νοδε 1 ηασ φιρεδ  Φιρστ λινεαρ εϖολυτιον: (φ 1, φ 2)= (0, φ 0) το (1−

φ 0,0)

 Νοω νοδε 2 ισ ρεσετ το φ 2= 0 ανδ νοδε 1 ϕυµπσ το

φ1 = ηφ(φ 0)= α φ 1 + β = − α φ 0 + (α+β)

 Ονλψ ϖαλιδ ιφ νοδε 2 δοεσ ρεαχη φ τηνοτ φορχε νοδε 1 το φιρε

ιµµεδιατελψ, ιν ωηιχη χασε βοτη πυλσεσ αρε σψνχηρονιζεδ. Τηισ µεανσ φ 0 ∈ ]1− φ λ,1[ ωιτη φ λ := (1− β)/α (φ λ ισ χαλλεδ χηαραχτεριστιχ ηοριζον)

Convergence analysis (2)

 Αφτερ φιρινγ οφ νοδε 2 τηε σψστεµ ισ ιν στατε

(φ1,φ2)= (ηφ(φ0),0).

 Τηε νεξτ νοδε το φιρε ισ νοδε 1.  Νοω, φ1 = 0, ανδ φ1 ισ οβταινεδ βψ

ηΡ(φ0)= ηφ (ηφ(φ0))= α(1− ηφ(φ0))+β= α2 φ0+(1− α)(α+β)

 ωηιχη ισ ονλψ ϖαλιδ ιφ τηε ιντιαλ πηασε οφ νοδε 2 ισ ιν

ιντερϖαλ ]1− φλ,1− φλ(1− 1/α)[, (οτηερωισε σψνχηρονιζατιον ηασ αλρεαδψ τακεν πλαχε)

 Ωε χαν νοω στυδψ τηε ρετυρν µαπ ηΡ(φ0) ωηετηερ ιτ ηασ

σταβλε φιξ ποιντσ φ0 = ηΡ(φ0)., σοµετηινγ λικε ηΡ (ηΡ(φ0)), ετχ.

Convergence analysis (3)

 Ρετυρν µαπσ:

 Study dynamics of a recurrent system

xi+1 = f(xi), xi ∈[0,1] (automorphism)

 Use Verhulst diagram (‘cobweb’):

xi+1 f(xi) xi+1 f(xi) x2 x3 x3

Convergence Analysis (4)

 Τηε ρετυρν µαπ ατ α φιξ ποιντ ξφιξ ωιτη ξφιξ=

φ(ξφιξ)

 Φιξποιντσ αρε ιντερσεχτιον ποιντσ ωιτη βισεχτριξ  Αφτερ α σµαλλ περτυρβατιον οφ ξφιξ :

 for stable fixpoint the system bounces back to point;  for instable fixpoint it moves away from fixpoint

 Σλοπε |φ (ξ

’

φιξ)| < 1 = = > σταβλε;

Σλοπε | φ (ξ ’

φιξ)| > 1 = = > ινσταβλε

xi f(xi) xi+1 f(xi) x fix1+ε xfix2 xfix1 xi+1 xi f(xi) xi+1 x fix2+ε

Verhulst diagram animation

http://en.wikipedia.org/wiki/Cobweb_diagram

Convergence Analysis (5)

 Ρετυρν µαπσ οφ ηΡ(φ0) διφφερεντ α ανδ β  Τηε τωο φιξποιντσ

ατ λεφτ ανδ ριγ ητ βουνδαρψ αρε σταβλε: φ (ξ) = 0<1 ’

 Τηε φιξ ποιντ ιν τηε

µιδδλε ισ αλωαψσ ινσταβλε: φ (ξ)= ’ α2>1

 Ηενχε: Σψστεµ ωιλλ

πραχτιχαλλψ αλωαψσ χονϖεργ ε το φ0= 1 ορ φ 1= 1, ωηερε σψνχηρονψ ισ ρεαχηεδ.

More than two oscillators

 Μορε διφφιχυλτ το

αναλψζε. Προοφ ιν [1].

 Πηασε πορτραιτ υσεφυλ

τοολ φορ ϖισυαλιζατιον

 Φιρστ σµαλλ χλυστερσ

σψνχηρονιζε (γ ρεψ ποιντσ).

 Τηεν χλυστερσ

σψνχηρονιζε

[1] R. Mirollo and S. Strogatz, “Synchronization of pulse-coupled Biological oscillators,” SIAM J. APPL. MATH, vol. 50, no. 6, pp. 1645–1662, Dec. 1990. Phase portrait for N=20 nodes with

• ne node firing; grey nodes will synchronize

with the black nodes in next step.

Convergence and time to synchronicity

 ιτ ωασ σηοων ιν [1] τηατ ιφ τηε νετωορκ ισ

φυλλψ µεσηεδ ανδ ‘ ’ ε> 0 ανδ β> 0, τηε σψστεµ αλωαψσ χονϖεργ εσ

 αλλ οσχιλλατορσ ωιλλ φιρε ασ ονε

ινδεπενδεντλψ οφ ινιτιαλ χονδιτιονσ.

 τηε τιµε το σψνχηρονψ ισ ινϖερσελψ

προπορτιοναλ το τηε προδυχτ βε ιν

[1] R. Mirollo and S. Strogatz, “Synchronization of pulse-coupled Biological oscillators,” SIAM J. APPL. MATH, vol. 50, no. 6, pp. 1645–1662,

• Dec. 1990.

0.2 0.4 0.6 0.8 1 200 400 600 800 1000 1200 1400 1600 1800 2000 0.2 0.4 0.6 0.8 1 200 400 600 800 1000 1200 1400 1600 1800 2000

clear; t=1; TMAX=20000; PHITH=2000; PHI1(1)=1000; PHI2(1)=400; d=1; epsilon=0.1; b =0.8; alpha = exp(b*epsilon); beta = (exp(b*epsilon)-1) / (exp(b)-1); for (i=1:TMAX-1) PHI1(t+1) = PHI1(t) + d; PHI2(t+1) = PHI2(t) + d; if (PHI1(t+1) > PHITH) PHI1(t+1) = 0; PHI2(t+1) = … min(alpha*PHI2(t) + beta, PHITH ); end if (PHI2(t+1) > PHITH) PHI2(t+1) = 0; PHI1(t+1) = … min(alpha*PHI1(t) + beta, PHITH ); end t=t+1 end figure(6); plot(PHI1, 'blue'); figure(7); plot(PHI2, 'red');

Experiment in MATLAB

 Τηε πηασε φυνχτιον σψνχηρονιζεσ χα. ατ τιµε στεπ 1.3, 5

τη ιτερατιον

Wireless sensor networks

 ∆εφ.: Α ωιρελεσσ σενσορ

νετωορκ (ΩΣΝ) ισ α ωιρελεσσ νετωορκ χονσιστινγ οφ σπατιαλλψ διστριβυτεδ αυτονοµουσ δεϖιχεσ υσινγ σενσορσ το χοοπερατιϖελψ µονιτορ πηψσιχαλ ορ ενϖιρονµενταλ χονδιτιονσ, συχη ασ τεµπερατυρε, σουνδ, ϖιβρατιον, πρεσσυρε, µοτιον ορ πολλυταντσ.

 Αδ−

ηοχ νετωορκσ: διστριβυτεδ ρανδοµλψ ανδ χοµµυνιχατε το νειγ ηβορσ; οφτεν αχχυ− ενεργ ψ ισ λιµιτεδ.

Source: Wikipedia Observer

Application to ad-hoc networks

 Τηε ΠΧΟ σψνχηρονιζατιον σχηεµε δεσχριβεδ πρεϖιουσλψ

χαν βε αππλιεδ το ωιρελεσσ σψστεµσ.

 Εξχηανγ ε οφ ινφορµατιον ρεθυιρεσ ενεργ ψ ανδ χαν βε δονε

ονλψ ιν χερταιν σψνχηρονιζεδ τιµε σλοτσ. Νοδεσ χαν βε ον λοω ενεργ ψ λεϖελ (ηιλβερνατε) φορ τηε ρεµαινινγ τιµε.

 Τηισ σαϖεσ ενεργ ψ ανδ ινχρεασεσ οπερατιον τιµε ανδ

ρανγ ε οφ αυτονοµουσ βαττερψ δριϖεν ωιρελεσσ αδ− ηοχ νετωορκσ

 Σψνχηρονιζατιον µυστ τακε ιντο αχχουντ πραχτιχαλ

προβλεµσ, µαινλψ χαυσεδ βψ δελαψσ …

Firefly synchronizations in wireless networks literature

 Α. Τψρρελλ, Γ. Αυερ, ανδ Χ. Βεττστεττερ,

“Φιρεφλψ σψνχηρονιζατιον ιν αδ ηοχ νετωορκσ,” ιν Προχ. ΜιΝΕΜΑ Ωορκσηοπ 2006, Φεβ. 2006. (βασιχ ιδεασ)

 Ψ.−

Ω. Ηονγ ανδ Α. Σχαγ λιονε, Α σχαλαβλε “ σψνχηρονιζατιον προτοχολ φορ λαργ ε σχαλε σενσορ νετωορκσ ανδ ιτσ αππλιχατιονσ, ” ΙΕΕΕ ϑουρναλ ον Σελεχτεδ Αρεασ ιν Χοµµυνιχατιονσ, ππ. 1085 1099, Μαψ 2005. – (προοφσ φορ δελαψ τρεατµεντ)

 Αλεξανδερ Τψρελλ, Φιρεφλψ σψνχηρονιζατιον ιν ωιρελεσσ

“ νετωορκσ”, ∆ισσερτατιον, Υνιϖερσιτ τ Κλαγ ενφυρτ, 2009 (χοµπρεηενσιϖε οϖερϖιεω, ρανγ ε οφ αππλιχατιονσ)

Propagation delays

 Ιφ α προπαγ ατιον δελαψ Τ0 οχχυρσ βετωεεν τωο πυλσε

χουπλεδ οσχιλλατορσ, τηε σψστεµ χαν βεχοµε ινσταβλε

 Τηε πυλσε οφ ονε οσχιλλατορ χουλδ χαυσε τηε οτηερ

οσχιλλατορ το τρανσµιτ αφτερ Τ0, ανδ τηισ τρανσµιττεδ πυλσε χαυσεσ τηε φιρστ οσχιλλατορ το φιρε αγ αιν αφτερ Τ0, ανδ σο ον.

 Το αϖοιδ τηισ αϖαλανχηε εφφεχτ α ρεφραχτορψ περιοδ

οφ δυρατιον Τρεφρ νεεδσ το βε αδδεδ αφτερ τρανσµισσιον.

 ∆υρινγ τηισ περιοδ, τηε πηασε φυνχτιον οφ α νοδε σταψσ

εθυαλ το 0 ανδ ισ νοτ µοδιφιεδ ιφ ρεχειϖινγ α πυλσε

 Σταβιλιτψ ισ µαινταινεδ ονλψ ιφ εχηοεσ αρε νοτ

ρεχειϖεδ, ωηιχη τρανσλατεσ το α χονδιτιον Τρεφρ > 2 • Τ0

• U. Ernst, K. Pawelzik, and T. Geisel, “Synchronization induced by

temporal delays in pulse-coupled oscillators,” Physical Review Letters,

• vol. 74, no. 9, pp. 1570–1573, Feb. 1995.

Multiple delays

 Τ0: Προπαγ ατιον δελαψ: τιµε το προπαγ ατε φροµ αν

εµιττινγ νοδε το α ρεχειϖινγ νοδε. Τηισ τιµε ισ προπορτιοναλ το τηε διστανχε βετωεεν τωο νοδεσ.

 Τξ: Τρανσµιττινγ δελαψ: λενγ τη οφ τηε βυρστ. Ωηιλε

τρανσµιττινγ , α νοδε ισ ιν α τρανσµιτ στατε ανδ χαννοτ λιστεν το οτηερ σψνχηρονιζατιον µεσσαγ εσ.

 Τδεχ ∆εχοδινγ δελαψ: τιµε ρεθυιρεδ βψ τηε ρεχειϖερ

το δεχοδε α σψνχηρονιζατιον µεσσαγ ε.

 Τρεφρ Ρεφραχτορψ δελαψ: τιµε νεχεσσαρψ αφτερ

τρανσµιττινγ το µαινταιν σταβιλιτψ. Α νοδε ισ ιν ρεφραχτορψ στατε δυρινγ τηισ περιοδ.

Synchronization with multiple delays (1)

 Το χοµβατ τηε λοσσ οφ αχχυραχψ

τηε τρανσµιττερ ισ δελαψεδ ιν ιτσ τρανσµισσιον φορ α χερταιν τιµε Τωαιτ εθυαλ το: Τωαιτ = Τ − (Τξ + Τδεχ) ωηερε Τ δενοτεσ τηε σψνχηρονιζατιον περιοδ.

 Τηισ σχηεµε µοδιφιεσ τηε νατυραλ

οσχιλλατορψ περιοδ οφ αν οσχιλλατορ, ωηιχη ισ νοω εθυαλ το 2 •Τ.

 Τηε τιµε δυρινγ ωηιχη τηε πηασε

φυνχτιον ωιλλ ινχρεµεντ ισ ρεδυχεδ βψ τηε ωαιτινγ , τρανσµιττινγ ανδ ρεφραχτορψ δελαψσ. Ιτ ισ νοω εθυαλ το ΤΡξ= 2Τ− Τωαιτ− Τξ− Τρεφρ

Y.-W. Hong and A. Scaglione (2009)

Synchronization with multiple delays

 Ατ ινσταντ 0, οσχιλλατορ 1

ρεαχηεσ φτη. Ιτ ωαιτσ υντιλ τ1 = Τωαιτ βεφορε σταρτινγ το τρανσµιτ α σψνχηρονιζατιον βυρστ.

 Ατ τ2 = Τωαιτ + Τξ + Τδεχ = Τ,

οσχιλλατορ 2 ηασ συχχεσσφυλλψ ρεχειϖεδ ανδ δεχοδεδ τηε βυρστ.

 Ασ τηε τωο οσχιλλατορσ αρε

αλρεαδψ σψνχηρονιζεδ, ιτ ωιλλ φολλοω τηε σαµε σχηεµε ασ οσχιλλατορ 1 ανδ ωαιτ υντιλ τ3 = Τ + Τωαιτ βεφορε τρανσµιττινγ .

t2

Other types of synchronization and periodicity in nature

 Χοντρολλεδ βψ χεντραλ σιγ ναλ (παχεµακερ)

θυαρτζ χλοχκσ, – σεασονσ, µοντηλψ ρηψτηµσ (µοον), δαψ− νιγ ητ ρηψτηµσ (συν), σεασοναλ ρηψτηµ, ΧΠΥσ χλοχκ, χλασσιχαλ µυσιχ ορχηεστρα

 Αντιχιπατιον−

βασεδ, σενσε οφ ρηψτηµ ιν ηυµανσ ‘ ’ – ιµπροϖισατιον ιν ϕαζζ− ενσεµβλε, σινγ ινγ σοχχερ συππορτερσ ‘ ’

 Σουρχεσ οφ περιοδιχ σιγ ναλσ ιν νατυρε:

 Waves = periodic limit cycles in dynamical systems; fix points of f(f(x)), f(f(f(x)))

and so on.

 Intermittence/pseudo-periodicity: Regular peaks in chaotic systems followed

by periods of deterministic chaos, often periods are multiples of 3; caused by near tangential return maps of f(f(x)), f(f(x)), f(f(f(x))), and so on.

Summary

 Χουπλεδ Οσχιλλατορσ χαν εξπλαιν σποντανεουσ

σψνχηρονιζατιον οφ φιρεφλιεσ

 Τηε σαµε µοδελ χαν βε υσεδ φορ µοδελινγ

ηεαρτβεατ (10000 οσχιλλατορσ), βραιν χελλ σψστεµσ, εαρτηθυακεσ.

 Χοµπυτερ σχιενχε αππλιχατιον το ωιρελεσσ

σενσορ νετωορκ σψνχηρονιζατιον

 ∆ελαψ σ (νο ινσταντανεουσ τρανσµισσιον)

’ δεµανδ φορ αδαπτατιονσ συχη ασ ρεφραχτορψ τιµεσ