## Precision of Pulse-Coupled Oscillator Synchronization on FPGA-Based Radios

G¨ unther Brandner ∗ , Johannes Klinglmayr , Udo Schilcher , Dominik Egarter , and Christian Bettstetter ∗ ∗ ∗ ∗§

∗ University of Klagenfurt, Mobile Systems Group, Inst. of Networked and Embedded Systems, Klagenfurt, Austria § Lakeside Labs GmbH, Klagenfurt, Austria

E-Mail: guenther.brandner@aau.at

Abstract -The precision of synchronization algorithms based on the theory of pulse-coupled oscillators is evaluated on FPGAbased radios for the first time. Measurements show that such algorithms can reach precision in the low microsecond range when being implemented in the physical layer. Furthermore, we propose an algorithm extension accounting for phase rate deviations of the hardware and show that an improved precision below one microsecond is possible with this extension in the given setup. The resulting algorithm can thus be applied in ad hoc wireless systems for fully distributed synchronization of transmission slots or sleep cycles, in particular, if centralized synchronization is impossible.

This paper intends to advance this direction of research. In particular, we evaluate the synchronization precision that algorithms based on pulse-coupled oscillators achieve in practice when being integrated into the physical layer on a programmable radio platform. Furthermore, as a result of our experimental research, we gain further insight into the behavior of pulse-coupling in real-world wireless scenarios, and thus propose an extension to the theory, which intends to correct phase rate deviations.

Our main contributions are as follows:

Index Terms -Synchronization, pulse-coupled oscillators, PCO, wireless systems, firefly synchronization, self-organization, programmable radio.

## I. INTRODUCTION AND MOTIVATION

The mathematical modeling of pulse-coupled biological oscillators, as proposed in [1] inspired by [2], offers a fully decentralized and scalable approach for time synchronization. There is a broad spectrum of work on pulse-coupled oscillators in physics, biology, neuroscience, and other disciplines (see, e.g., [3]-[11] and references therein). The communications engineering community has been interested to transfer these results to self-organizing synchronization of nodes in wireless systems [12] for purposes such as slot and frame synchronization, scheduling of cooperative transmissions and sleep cycles, and distributed sensing. A one-to-one transfer is, however, infeasible due to the differences between biological and radio communications. Several extensions and modifications are required with respect to delays, noise, multihop communications, and synchronization words, to mention a few (see [12]-[20] and references therein).

Despite the conceptional and theoretical advances in the design of pulse-coupled oscillator synchronization algorithms for wireless systems, real-world performance studies and experimental proofs of concepts are largely missing. There only exist a few implementations on low-cost wireless sensor platforms (see [15], [21], [22]), whose results are of interest, but their synchronization precision is limited by hardware capabilities. For example, the 'fifty percentile group spread' in a system with 24 MicaZ motes reported in [15] is in the order of 100 μ s, which is insufficient for certain applications, such as slot and frame synchronization.

- · Providing a proof of concept by implementing three pulse-coupled oscillator algorithms on fieldprogrammable gate array (FPGA)-based radios
- · Proposing an extension to the synchronization algorithms accounting for phase rate deviations
- · Comparing the synchronization precision of these algorithms by measurements
- · Showing by experiments that precisions below one μ s can be achieved

To the best of our knowledge, this is the first lower-layer implementation and real-world performance study of recent pulse-coupled oscillator algorithms on programmable radios.

## II. SYNCHRONIZATION ALGORITHMS

We evaluate the synchronization precision of the following recently proposed synchronization algorithms:

- · Synchronization by Pagliari &amp; Scaglione (PS) [22], [23],
- · Synchronization with inhibitory and excitatory coupling with stochastic pulse emission (IES) [11].
- · Synchronization with inhibitory coupling and selfadjustment (SISA) [24] and

We also propose and evaluate a modified version of the IES algorithm that applies phase rate correction (IES ∗ ). The objective of all four algorithms is to synchronize the phases of oscillators.

The general procedure is as follows: The oscillator's phase φ is increased from zero to one. When φ reaches one, φ is reset and a pulse is emitted, either always or with probability p &lt; 1 depending on the algorithm. When receiving a pulse from another oscillator, an oscillator adjusts its own phase according to an update function H φ ( ) .

The absolute time is called t . The period τ i j denotes the delay between oscillator i and oscillator j , i.e., the time it takes from the start of a pulse at i until it is processed at j . Let τ min , τ max , and τ denote the minimum, maximum, and mean values of all delays, respectively. Furthermore, φ t ( ) is an oscillator's phase at time t and φ t ( + ) its phase infinitely short after t . The term ν i is the phase rate deviation of oscillator i , i.e., the speed of oscillator i compared to a reference oscillator; ν max is the maximum phase rate deviation of all oscillators. Let h t ( ) denote the function which maps durations t in seconds to the corresponding phase, i.e. h t ( ) = t t c , where t c is the cycle length in seconds.

Figure 1 specifies the four algorithms. We use

$$H _ { \text{PS} } ( \phi ) = \min ( 1, a _ { 1 } \phi + a _ { 0 } )$$

for PS with parameters a 0 = 1 and a 1 = exp(1) (strong coupling). These parameters for PS are chosen as they optimize the convergence speed. Choosing different parameters does not influence the achieved precision. For SISA we apply

$$H _ { \text{SISA} } ( \phi ) = ( 1 + \alpha ) \, \phi \, \bmod \, 1,$$

with α = 0 5 . which is also applied in [24]. IES uses H IES ( φ ) = ˜ H IES ( φ -h τ ( min ) mod 1) + h τ ( min ) mod 1 . ˜ H IES is, mutatis mutandis , a function of the form [25]

$$\widetilde { H } _ { \text{IES} } ( \phi ) = \begin{cases} h _ { 1 } ( \phi ) & \text{if $\phi_{\text{ref} } < \phi \leq \frac{1}{2}$,} \\ h _ { 2 } ( \phi ) & \text{if $\frac{1}{2} < \phi \leq 1$,} \end{cases}$$

where h 1 and h 2 are continuous functions that satisfy the following requirements:

- · d h 1 d φ &gt; 0 , d h 2 d φ &gt; 0 ,
- · h 1 (0 5) . ≤ 0 25 . -τ max -τ min ,
- · h 1 ( τ max ) = τ max ,
- · h 2 (0 5 . + ) ≥ 0 75 + ( . τ max -τ min ) , and

In the following we use

- · h 2 (1) = 1 .

$$h _ { 1 } ( \phi ) = a \cdot [ \phi - h ( \tau _ { \max } ) ] + h ( \tau _ { \max } )$$

and

$$h _ { 2 } ( \phi ) = b \cdot [ \phi - 1 ] + 1,$$

with a = 1 4 -2 h τ ( max ) -h τ ( min ) 1 2 -h τ ( max ) and b = 1 2 + 2 h τ ( min ) -2 h τ ( max ) . These functions fulfill all requirements. IES ∗ uses H IES ∗ ( φ ) = ˜ H IES ( φ -h τ ( ) mod 1) + h τ ( ) mod 1 . For both IES and IES ∗ we apply p = 1 2 as the sending probability.

## III. IMPLEMENTATION ON PROGRAMMABLE RADIO

We implement all synchronization algorithms on WARP boards [26], which are FPGA-based radios. A custom singlecarrier physical layer is programmed with 5 MHz bandwidth and binary phase shift keying (BPSK). Boards operate at 2 4 . GHz and use a peak transmit power of about 20 dBm. The overall structure of the transceiver is shown in Figure 2(a). All components are implemented directly on the FPGA. On the transmitter side, the packetizer and modulator build the packet after they receive a trigger signal from the synchronization

- 1) Whenever φ t ( ) = 1 , the oscillator sends a pulse.
- 2) The refractory interval is

$$[ 0, \phi _ { \text{ref} } = 2 \, ( 1 + \nu _ { \max } ) \, h ( \tau _ { \max } ) ].$$

- 3) Upon reception of a pulse at time t ∗ :

$$\phi ( t _ { * } ^ { + } ) = \begin{cases} \phi ( t _ { * } ) & \text{if $\phi(t_{*}) \leq \phi_{ref}$}, \\ H _ { \text{PS} } ( \phi ( t _ { * } ) ) & \text{else}. \end{cases}$$

(a) PS

- 1) Whenever φ t ( ) = 1 , the oscillator adjusts its phase to φ t ( + ) = H SISA (1) and sends a pulse.
- 2) The refractory interval is

$$[ 0, \phi _ { r e f } = H _ { S I S A } ( 1 ) + 2 \left ( 1 + \nu _ { \max } \right ) h ( \tau _ { \max } ) ]. \quad \Big |$$

- 3) Upon reception of a pulse at time t ∗ :

$$\underline { \phi ( t _ { * } ^ { + } ) } = \begin{cases} \phi ( t _ { * } ) & \text{if $\phi(t_{*}) \leq \phi_{ref}$,} \\ H _ { S I S A } ( \phi ( t _ { * } ) ) & \text{else.} \end{cases}$$

(b) SISA

- 1) Whenever φ t ( ) = 1 , the oscillator sends a pulse with probability p &lt; 1 .
- 2) The refractory interval is

$$[ 0, \phi _ { \text{ref} } = ( 1 + \nu _ { \max } ) \, h ( \tau _ { \max } ) ].$$

- 3) Upon reception of a pulse at time t ∗ :

$$\phi ( t _ { * } ^ { + } ) = \begin{cases} \phi ( t _ { * } ) & \text{if $\phi(t_{*}) \leq \phi_{ref}$,} \\ H _ { \text{IES} } ( \phi ( t _ { * } ) ) & \text{else.} \end{cases}$$

(c) IES

- 1) Whenever φ t ( ) = 1 and no pulse has been received within the last τ -τ min seconds, the oscillator sends a pulse with probability p &lt; 1 .
- 2) The refractory interval is

$$[ 0, \phi _ { \text{ref} } = ( 1 + \nu _ { \max } ) \, h ( \tau _ { \max } ) ].$$

- 3) Upon reception of a pulse at time t ∗ :

$$\phi ( t _ { * } ^ { + } ) = \begin{cases} \phi ( t _ { * } ) & \text{if $\phi(t_{*}) \leq \phi_{ref}$,} \\ H _ { \text{IES} ^ { * } } ( \phi ( t _ { * } ) ) & \text{else.} \end{cases}$$

(d) IES ∗

Fig. 1: Synchronization algorithms logic. A modulated packet is fed into an interpolator and upconverter, and finally transmitted over the air. As we cannot send infinitely short pulses, as often assumed in theory, we send short packets instead. These packets have a length of 12 bytes, where the first 8 bytes are used for setting the receiving gains of the hardware (agc) and for mitigating carrier frequency offsets (cfo). The remaining 4 bytes represent a

Fig. 2: FPGA design

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-ac2ca088-88f8-4393-a44b-10ac2d22714c.jpg)

synchronization word consisting of pseudorandom bytes. The transmit duration for a packet is 19 2 . μ s.

c i in IES . ∗ These correction terms are determined for each board individually by manually measuring their phase rates with respect to a reference phase. The terms c i are then stored in the board's memory and applied during synchronization in the following way: at each clock cycle the correction factor c i is accumulated in a dedicated accumulator and the output of the accumulator is then added to the output of the counter.

The synchronization logic on the FPGA consists of implementations of the algorithms discussed above. The oscillator component (Figure 2(b)) replicates the oscillator on the board. The main part, generating the phase of the oscillator, is a 22 -bit wrap-around counter running at a clock frequency of 40 MHz. Thus, the cycle duration t c , i.e., the time it takes for the counter to increment from 0 to 2 22 -1 , is about 104 86 . ms. To get a value between zero and one we reinterpret the output as a fractional number by multiplying with 2 -22 . After a pulse is detected, the new phase H φ ( ) , determined by the synchronization algorithm component, is forwarded to the oscillator block and an accumulator is used to adjust to the new phase. Note that the accumulator only processes the value on its input if it is enabled, i.e. if a synchronization word has been detected. The output of the accumulator always reflects the current value, independent on whether or not the accumulator is enabled.

Note that phase rates depend on environmental factors, e.g., temperature. The purpose of applying these correction factors is to showcase the influence of phase rate deviations on the achievable synchronization precision. As future work we plan to propose a fully decentralized algorithm that not only synchronizes phases but also phase rates. Note that for all other algorithms, besides IES ∗ , we set c i to zero.

On the receiver side, the signals' inphase ( I ) and quadrature ( Q ) components are used to estimate and set the amplifier gains of the boards. The downconverter brings signals to the baseband. We implement a non-data aided algorithm [27] for removing carrier frequency offsets. To detect the synchronization word, a correlator implemented as an FIR filter is applied.

## IV. MEASUREMENT RESULTS

Six radios are setup to form a fully-connected network with six nodes. The distances between all nodes are a few meters. Thus, all nodes can generally receive packets of all other nodes, but packets might be lost due to, e.g., interference from co-located WLANs.

## A. Delays τ and Phase Rate Deviations ν

We measure the delays τ i j for various sender-receiver pairs and analyze the overall empirical probability density function (epdf). The epdf is derived by the method of kernel density estimation [28]. Figure 3(a) shows the epdf of τ based on

Due to manufacturing tolerances, boards exhibit phase rate deviations. These deviations limit the achievable synchronization precision. As a countermeasure we add correction terms six sender-receiver pairs and 10 000 transmissions each. Experiments show no significant difference between the senderreceiver pairs: The delay is always between τ min = 21 7 . μ s and τ max = 22 2 . μ s, and the average delay is τ = 21 92 . μ s; we use these values in the synchronization algorithms. The values presented for τ are accurate to ± 25 ns.

(a) Empirical pdf of delay τ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-3dd22488-7470-43aa-808f-33c7aa1db831.jpg)

(b) Phase rate deviations ν

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-c7509c60-d337-4ac4-a172-cd6f9318ef33.jpg)

Fig. 3: Empirical pdf of delay τ and phase rate deviations ν

Figure 3(b) shows the phase rate deviations ν in parts per million (ppm) from a reference phase for all six nodes. Only five lines are visible as two nodes have the same deviation of about 1 8 . ppm. The phases of all nodes run faster than the reference phase; the deviations range from 1 8 . to 6 ppm and remain constant over time. The values presented for ν are accurate to ± 0 25 . ppm.

## B. Synchronization Precision

The synchronization precision at time t in terms of the maximum phase difference between all nodes is evaluated as [11]

$$\Gamma ( t ) = t _ { c } \cdot \max _ { i, j } \left \{ \min [ | \phi _ { i } ( t ) - \phi _ { j } ( t ) |, 1 - | \phi _ { i } ( t ) - \phi _ { j } ( t ) | ] \right \}. \quad \underset { 1 4 9 } { [ 4 ] } \, \underbrace { L. \, 1 } { \text{net} }$$

Figure 4 shows the synchronization precision for n ∈ { 2 4 6 , , } nodes. Results are based on 100 synchronization runs, where prior to each run the phases of all nodes are randomly initiated. The x -axes show the cycle number. The cycle duration is 52 43 . ms for SISA and 104 86 . ms for all other algorithms. The shown values are accurate to ± 25 ns.

These measurement results can be interpreted as follows: PS converges very quickly to a synchronization precision of about 21 μ s. For SISA, the speed of convergence decreases with increasing n , which is due to the fact that packets lost over the wireless link can cause a deterioration of the synchronization precision: Whenever a node detects a synchronization word from some other node, and is not currently in refractory, it halves its current phase. In case that not all nodes detect the synchronization word, we have the situation that some nodes halve their phases, while others do not halve their phases. The convergence of IES is slower than that of PS with the given parameters, however, it converges to a mean precision of about 1 5 . μ s ( n = 2 ), 2 μ s ( 4 ) and 4 μ s ( 6 ). IES ∗ achieves a precision of about 200 ns ( n = 2 ), 400 ns ( 4 ) and 600 ns ( 6 ).

The fact that PS and SISA synchronize less precisely than IES in this setup is likely due to two reasons: (i) propagation delays are not considered in those algorithms and (ii) nodes cannot hear other nodes when sending. This result confirms that stochastic communication of synchronization words is an important design feature (see, e.g., [11]). This feature could in principle also be applied to PS.

## V. CONCLUSIONS

Measurement results of pulse-coupled oscillator synchronization implemented on FPGA-based radios show that the synchronization precision can reach values below one μ s. Key factors for reaching this precision are the explicit consideration of propagation and processing delays, the stochastic nature in communications of synchronization words, and a phase rate correction. The latter mitigates precision limitations caused by phase rate deviations of the hardware.

## ACKNOWLEDGEMENTS

This work was supported by Lakeside Labs with funding from ERDF, KWF, BABEG under grant 20214/23794/35530. The work of Johannes Klinglmayr was partially supported by the Linz Center of Mechatronics (LCM) in the framework of the Austrian COMET-K2 program.

## REFERENCES

- [1] R. E. Mirollo and S. H. Strogatz, 'Synchronization of pulse-coupled biological oscillators,' SIAM J. Appl. Math. , vol. 50, no. 6, pp. 16451662, 1990.
- [3] J. Buck, E. Buck, J. Case, and F. Hanson, 'Control of flashing in fireflies. V . Pacemaker synchronization in pteroptyx cribellata ,' J. Comp. Physiology A , vol. 144, pp. 630-633, Sept. 1981.
- [2] C. S. Peskin, Mathematical Aspects of Heart Physiology , pp. 268- 278. Courant Institute of Mathematical Sciences, 1975.
- [4] L. F. Abbott and C. A. van Vreeswijk, 'Asynchronous states in neural networks of pulse-coupled oscillators,' Phys. Rev. E , vol. 48, pp. 14831490, 1993.
- [5] C. A. van Vreeswijk, L. F. Abbott, and G. B. Ermentrout, 'When inhibition, not excitation, synchronizes neural firing,' J. Comp. Neurosci. , vol. 1, pp. 313-321, Dec. 1994.

(i) n = 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-e67b0f15-fc29-4ec0-92c5-95d792aa72f6.jpg)

(ii)

n

= 4

(iii) n = 6

(a) PS

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-40c45138-3118-4c94-9600-e6948ca67e03.jpg)

(i) n = 2

(ii) n = 4

(b) SISA

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-ec2eb0a3-6058-4995-99b1-099a36986f5e.jpg)

(i) n = 2

(ii)

n

= 4

(iii) n = 6

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-d15da023-e622-4343-8581-1658427ea8aa.jpg)

(iii) n = 6

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-5e26c22e-3842-417e-a05c-a96c73134ecd.jpg)

(c) IES

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0652v2-figure-c5482d10-2fb2-460a-b29b-6a544ea82075.jpg)

(d) IES ∗

Fig. 4: Synchronization precision

- [6] U. Ernst, K. Pawelzik, and T. Geisel, 'Synchronization induced by temporal delays in pulse-coupled oscillators,' Phys. Rev. Lett. , vol. 74, pp. 1570-1573, Feb. 1995.
- [8] U. Ernst, K. Pawelzik, and T. Geisel, 'Delay-induced multistable synchronization of biological oscillators,' Phys. Rev. E , vol. 57, pp. 21502162, Feb. 1998.
- [7] W. Gerstner, 'Rapid phase locking in systems of pulse-coupled oscillators with delays,' Phys. Rev. Lett. , vol. 76, pp. 1755-1758, Mar. 1996.
- [9] D. Golomb and G. B. Ermentrout, 'Bistability in pulse propagation in networks of excitatory and inhibitory populations,' Phys. Rev. Lett. , vol. 86, pp. 4179-4182, Apr. 2001.
- [11] J. Klinglmayr, C. Kirst, C. Bettstetter, and M. Timme, 'Guaranteeing global synchronization in networks with stochastic interactions,' New J. Phys. , vol. 14, July 2012. 073031.
- [10] J. Nishimura and E. J. Friedman, 'Robust convergence in pulse-coupled oscillators with delays,' Phys. Rev. Lett. , vol. 106, May 2011. 194101.
- [12] R. Mathar and J. Mattfeldt, 'Pulse-coupled decentral synchronization,' SIAM J. Appl. Math. , vol. 56, no. 4, pp. 1094-1106, 1996.
- [14] D. Lucarelli and I.-J. Wang, 'Decentralized synchronization protocols with nearest neighbor communication,' in Proc. ACM Conf. Embedded Networked Sensor Systems (SenSys) , (Baltimore, MD, USA), Nov. 2004.
- [13] Y.-W. Hong and A. Scaglione, 'A scalable synchronization protocol for large scale sensor networks and its applications,' IEEE J. Sel. Areas Commun. , vol. 23, no. 5, pp. 1085-1099, 2005.
- [15] G. Werner-Allen, G. Tewari, A. Patel, M. Welsh, and R. Nagpal, 'Fireflyinspired sensor network synchronicity with realistic radio effects,' in Proc. ACM Conf. Embedded Networked Sensor Systems (SenSys) , (San Diego, CA, USA), pp. 142-153, Nov. 2005.
- [17] O. Simeone, U. Spagnolini, Y. Bar-Ness, and S. H. Strogatz, 'Distributed synchronization in wireless networks,' IEEE Signal Process. Mag. , vol. 25, no. 5, pp. 81-97, 2008.
- [16] A. Tyrrell, G. Auer, and C. Bettstetter, 'Fireflies as role models for synchronization in ad hoc networks,' in Proc. Intern. Conf. on BioInspired Models of Network, Information, and Computing Systems (BIONETICS) , (Cavalese, Italy), Dec. 2006.
- [18] A. Tyrrell, G. Auer, and C. Bettstetter, 'Emergent slot synchronization in wireless networks,' IEEE Trans. Mobile Comput. , vol. 9, pp. 719-732, May 2010.
- [20] Y. Wang, F. N´ u˜ez, and F. J. Doyle III, 'Increasing sync rate of pulsen coupled oscillators via phase response function design: Theory and application to wireless networks,' IEEE Trans. Control Syst. Technol. , vol. 21, pp. 1455-1462, July 2013.
- [19] Y. Wang, F. N´ u˜ez, and F. J. Doyle III, 'Statistical analysis of the pulsen coupled synchronization strategy for wireless sensor networks,' IEEE Trans. Signal Process. , vol. 61, pp. 5193-5204, Nov. 2013.
- [21] R. Leidenfrost and W. Elmenreich, 'Firefly clock synchronization in an 802.15.4 wireless network,' EURASIP J. Embed. Syst. , vol. 2009, 2009. 186406.
- [23] R. Pagliari and A. Scaglione, 'Correction to 'scalable network synchronization with pulse-coupled oscillators',' IEEE Trans. Mobile Comput. , vol. 10, no. 5, 2011.
- [22] R. Pagliari and A. Scaglione, 'Scalable network synchronization with pulse-coupled oscillators,' IEEE Trans. Mobile Comput. , vol. 10, no. 3, pp. 392-405, 2011.
- [24] J. Klinglmayr and C. Bettstetter, 'Self-organizing synchronization with inhibitory-coupled oscillators: Convergence and robustness,' ACM Trans. Auton. Adapt. Syst. , vol. 7, pp. 30:1-30:23, Oct. 2012.
- [26] 'WARP Project: Wireless Open-Access Research Platform, Rice University.' http://warp.rice.edu; accessed December 14, 2014.
- [25] J. Klinglmayr, Self-organizing network synchronization: Convergence and robustness in pulse-coupled oscillator systems . PhD thesis, AlpenAdria-Universit¨ at, Klagenfurt, Austria, 2013.
- [27] J. R. Barry, E. A. Lee, and D. G. Messerschmitt, Digital Communications . Springer, 2004.
- [28] E. Parzen, 'On estimation of a probability density function and mode,' Ann. Math. Statist. , vol. 33, no. 3, pp. 1065-1076, 1962.
Günther Brandner Johannes Klinglmayr Udo Schilcher Dominik Egarter Christian Bettstetter
Precision of Pulse-Coupled Oscillator Synchronization on FPGA-Based Radios
The precision of synchronization algorithms based on the theory of
pulse-coupled oscillators is evaluated on FPGA-based radios for the first time.
Measurements show that such algorithms can reach precision in the low
microsecond range when being implemented in the physical layer. Furthermore, we
propose an algorithm extension accounting for phase rate deviations of the
hardware and show that an improved precision below one microsecond is possible
with this extension in the given setup. The resulting algorithm can thus be
applied in ad hoc wireless systems for fully distributed synchronization of
transmission slots or sleep cycles, in particular, if centralized
synchronization is impossible.