Research Article
Distributed synchronization in large-scale wireless sensor networks using group consensus protocol
International Journal of Distributed Sensor Networks 2017, Vol. 13(11) Ó The Author(s) 2017 DOI: 10.1177/1550147717718114 journals.sagepub.com/home/ijdsn
Yanliang Jin, Xiaoshuai Zhang and Bin Yao
Abstract This article analyzes the clock model of sensor nodes and the basic principle of time synchronization. Then, it systematically introduces the existing time synchronization mechanisms and compares their advantages and disadvantages. In order to solve the problems of consensus-based time synchronization, such as propagation delay and convergence speed, Group Consensus Time Synchronization algorithm is proposed. In this protocol, minimum data exchange is achieved for both frequency and phase offset compensation, and at the same time, the propagation delay is considered and compensated for. The simulation results show that the protocol can adapt to dynamic topology and be suitable to large-scale wireless sensor networks. Keywords Distributed synchronization, group consensus, large scale, wireless sensor network, dynamic topology
Date received: 27 September 2016; accepted: 15 May 2017 Handling Editor: Luis Orozco-Barbosa
Introduction Time synchronization is an important issue, which must be addressed in all distributed systems, especially in wireless sensor networks (WSNs). It is crucial for WSNs when performing a number of fundamental operations,1 such as low-power medium access control (MAC) protocol operation, reliable data fusion,2 target tracking, and precise positioning. The research on time synchronization has a long history. In the past few decades, many algorithms and mechanisms have been proposed and applied. However, the existing synchronization mechanism is not suitable for WSNs due to the special nature of WSNs. In WSNs, time synchronization refers to exchanging and adjusting clock information of all nodes through single-hop or multi-hop manner to make all nodes running with the same clock. In order to achieve time synchronization in the WSN, the following challenges need to be overcome:3
1.
The uncertainty of communication time delay. The majority of clock synchronization algorithms are required to obtain clock information of other nodes through information communication to achieve synchronization. However, a variety of communication delay affects the transmission of information between the sensor nodes in WSNs. The communication delay
Key laboratory of Specialty Fiber Optics and Optical Access Networks, Joint International Research Laboratory of Specialty Fiber Optics and Advanced Communication, Shanghai Institute for Advanced Communication and Data Science, Shanghai University, Shanghai, China Corresponding author: Yanliang Jin, Key laboratory of Specialty Fiber Optics and Optical Access Networks, Joint International Research Laboratory of Specialty Fiber Optics and Advanced Communication, Shanghai Institute for Advanced Communication and Data Science, Shanghai University, Room 1211, Building Xingjian, No.149, Yanchang Road, Zhabei District, Shanghai 200444, China. Email:
[email protected]
Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (http://www.uk.sagepub.com/aboutus/ openaccess.htm).
2
International Journal of Distributed Sensor Networks
2.
3.
4.
brings some difficulties and challenges to time synchronization in WSNs. During the transmission process, the main communication delay includes send time, access time, transmission time, propagation time, reception time, and received time. The largest communication delay is the access time, which is incurred in the MAC layer waiting for access to the transmit channel. Scalability. In sensor network applications, both the size of network and the density of nodes are varied, so time synchronization mechanism must be able to adapt to these changes. Robustness. All sorts of reasons may cause the sensor node failure and the environment affects the communication quality of wireless link, so the time synchronization mechanism must have good robustness. Convergence Speed. It is known to all that the sensor network topology is dynamic and sensor nodes are energy limited. Therefore, time synchronization establishing time should be short enough.
Related work So far, many algorithms can solve some of the problems described above. All of these agreements have some common characteristics: the simple packet transport protocol, exchanging clock information between nodes, reducing the effect of communication delay, and the implementation of the corresponding different schemes and algorithms. They are classified into two types: centralized and distributed time synchronization. The centralized time synchronization algorithms, such as Reference Broadcasts Synchronization (RBS),4 Timing-sync Protocol for Sensor Networks (TPSN),5 and Flooding Time Synchronization Protocol (FTSP),6 have many advantages including fast convergence speed, low complexity, and high convergence precision. However, these synchronization algorithms are based on structure type and mostly depend on the reference node. So, their robustness as well as scalability is poor. For example, if a node dies or a new node is added, it is necessary to rebuild the tree or cluster, which requires additional implementation overhead and may lead to worse synchronous precision. Distributed time synchronization algorithms, such as Reachback Firefly Algorithm (RFA)7 and Distributed Time Synchronization Protocol (DTSP),8 achieve time synchronization using the time information of neighbor nodes. In order to overcome the limitations of centralized time synchronization algorithm and the distributed algorithm, the concept of consistency is introduced in sensor networks. Average Time Synchronization Protocol (ATS)9 takes clock phase shift and frequency shift into
consideration. At the same time, each node in the network complies with the same consensus protocol. Consensus Clock Synchronization (CCS)10 also uses consensus protocol to compensate for clock phase offset. The main idea is that weighted coefficient is introduced into the protocol and each node can get own clock offset, ultimately correct itself through the information. In a sense, the CCS is an enhanced version of ATS, but the drawback is that it still has a slow convergence rate. Time Synchronization using Max and Average Consensus Protocol (TSMA),3 based on ATS and CCS, achieves time synchronization by the way of the combination of the maximum consensus and average consensus, which can make the clock converge to a virtual global clock. Wu et al.11 use two-order consistency algorithm to improve the convergence speed of the algorithm, but it also increases the complexity of the algorithm. Consistency algorithms discussed above have ignored a major problem: propagation delay. However, it is significant in the actual time synchronization. The algorithm in the work by Leng and Wu12 describes two cases of fixed and random distribution of propagation delay. Du and Wu13 put forward a kind of asynchronous synchronization algorithm. Different nodes can update their estimated value based on different frequencies. However, the time and space complexity of these two algorithms are relatively high, which need more storage and computing resources. Due to these reasons, both of them are not suitable for large-scale WSNs. So, this article presents a fully novel DTSP in WSN called Group Consensus Time Synchronization Protocol (GCTS). It has simple calculation, fast convergence rate, and low energy consumption.
System model The algorithm based on GCTS mainly includes two parts: the compensation for the clock frequency offset and phase offset. The compensations are realized through the exchange and calculation of time information between sensor nodes. Finally, all nodes in the network will converge to a virtual clock. In the process of information exchange, the propagation delay may affect the accuracy of the time synchronization. So, this article proposes the two-way information exchange to attain the time synchronization.
The propagation delay model As shown in Figure 1, node 1 requires the time information of node 2. Before node 1 sends a request instruction, it records its instantaneous timestamp t1 first. Node 2 receives a request instruction from node 1 and records its timestamp at once. Then, node 2 transmits time synchronization information packet back to node
Jin et al.
3
Figure 1. Two-way message exchange.
1 including transmission time t2 and receiving time t3 . When receiving the packet, node 1 records its timestamp t4 . In this way, the propagation delay between two nodes can be estimated using these four timestamps as shown in equation (1) d=
(t4 t1 ) (t3 t2 ) 2
ð1Þ Figure 2. Point-to-point message exchange.
Then, the estimated propagation delay can be used to compensate for the clock frequency offset and phase offset.
The correction of the clock frequency offset The clock frequency offset correction for the sensor nodes is designed to allow all nodes to run at the same frequency because if two sensor nodes are not running at the same frequency, the network even gets synchronization at the present time and the next moment will lose synchronization. Moreover, with the passage of time, the time offset between the two nodes will become bigger. Therefore, there would be a need for more frequent time synchronization information exchanges to realize time synchronization. A method14 that can use two-way synchronization information package exchange to estimate the frequency of the node is proposed. In addition, the frequency drift of the most of sensor nodes can be found in the specification. In the initial phase of the network, it is easy to obtain the clock frequency through the programming operation. In the WSN composed of N nodes, the clock frequency of node i(i = 1, 2, . . . , N ) in state kT is denoted as wi (k). Then, the clock frequency in state (k + 1)T is denoted as wi (k + 1). The clock frequency adjustment of node is affected by itself and surrounding nodes with information exchange. The collection of these nodes is defined as neighbor set named Ni . The specific frequency offset correction algorithm is as follows: Step 1. Node i periodically broadcasts a request packet to all neighbor nodes which makes all nodes on the broadcast range of node i transmit their own clock frequency wk ,k 2 Ni .
Step 2. Neighbor node returns a data packet with its clock frequency information. Step 3. Compensate for the propagation delay in this process. Step 4. The node performs the operation based on GCTS according to the clock frequency of itself and surrounding nodes. It means that the node calculates its new frequency offset value by averaging it with the estimates of its neighbors. Step 5. Update the clock frequency. Step 6. The node broadcasts a compensation table including the address and the compensation for the clock frequency of the corresponding node. Step 7. Neighbor nodes update their own clock frequency according to the received compensation package. Step 8. Repeat the process above. In the end, w of all nodes may tend to be consistent. The threshold refers to synchronization accuracy. In other words, clock synchronization error maintains the fluctuation at a particular value. This value is the threshold. In this article, the threshold is set to 25 ms. In the process of frequency compensation, GCTS algorithm not only uses the group information exchange mechanism to replace the original point-to-point mechanism but also considers the compensation for the propagation delay. As is shown in Figures 2 and 3, this article compares the group exchange mechanism with the point-to-point exchange mechanism. In Figure 2, through broadcasting synchronization packets, the consistent clock frequency between two nodes is realized. Node 2, node 3, and node 4 can
4
International Journal of Distributed Sensor Networks operation in the WSNs means that the local time of all or part of the nodes is updated. Assuming the starting time t = 0, the expression of all the nodes at time t after several operations is shown below 8 ^ t + e1 (t) c1 (t) = w > > > < c2 (t) = w ^ t + e2 (t) .. > . > > : ^ t + eN (t) cN (t) = w
ð3Þ
ei (t) represents the time offset of node i ^ t, ei (t) = ci (t) w
Figure 3. Group message exchange.
achieve the clock frequency through operation after receiving the synchronization information packet of the node 1. In this synchronization mode, the propagation delay is ignored. After a round of synchronization, the clock frequencies of the four nodes are still different but their standard deviation is reduced. This is a typical consensus algorithm called consensus time synchronization algorithm based on two nodes. As is shown in Figure 3, it is just the time synchronization method used in this article. Node 1 waits for a period of time to collect all the clock information from the surrounding nodes and then calculates the propagation delay and the consistent clock information of the entire group. This method can realize accurate and fast time synchronization.
The correction of the clock phase offset According to the previous analysis, after correcting the clock frequency offset, all local virtual clock estimators ^ and will run at the will eventually have the same drift w same speed. At this point, it is only necessary to compensate for possible phase offset. The local time of node i is expressed as ^ t + fi + m, ci (t) = w
i = 1, . . . , N
ð4Þ
Whether the sensor network synchronization is completed or not, the local time of all network nodes is always increasing with the passage of time. From equations (3) and (4) above, when the network node clock achieves synchronization, variable ei (t) is constant. Therefore, it is very convenient to analyze the evolution of the node clock synchronization using ei (t) instead of ci (t). In order to obtain the WSN synchronization, it is necessary to discrete node clock. After k round operations in WSNs, node’s local clock can be expressed as 8 ^ t(k) + e1 (k) c1 (k) = w > > > < ^ t(k) + e2 (k) c2 (k) = w ð5Þ .. > . > > : ^ t(k) + eN (k) cN (k) = w where ci (k) is the local time of the node i after k round operations. t(k) is the discrete time expression. ei (k) is the discrete time offset of the node i after k round operations ^ t(k), ei (k) = ci (k) w
i = 1, . . . , N
ð6Þ
Therefore, it is convenient to study ei (k) to determine whether the synchronization of WSNs is achieved or not. In the WSN composed of N nodes, the clock offset of node i(i = 1, 2, . . . , N) in state kT is denoted as ei (k). Then, the clock offset in state (k + 1)T is denoted as ei (k + 1). The purpose of discrete protocol control is to establish the relationship between the present state and the upper state of a system so that the relationship between the final convergence state and the initial state can be obtained by the recurrence relation. Equation (7) represents the group consensus protocol discussed in this article
ð2Þ
where ci (t) is the local time of node i at the moment of ^ is the crystal oscillator frequency of node i. t. w fi = ci (0) represents the initial time (t = 0). m is the random noise. In order to express more clearly, a single
i = 1, . . . , N
ei (k + 1) = ei (k) + a
X
aij (ej (k) ei (k))
ð7Þ
j2Ni
where a 2 (0, 1) is the initial tuning parameter. aij is related to adjacency matrix which describes the
Jin et al.
5
neighborhood relationship between nodes. The matrix A is an N 3 N matrix where N represents the number of sensor nodes. A simple commonly used definition is as follows: if the node i is adjacent to j, and node i can send information to j, then aij = 1, otherwise aij = 0. This protocol is used to correct the current state of the sensor node with the deviation of the last state. The specific clock offset correction algorithm is as the same as the correction of the clock frequency offset algorithm. The difference is transmitting the clock offset information ei (k) instead of the clock frequency offset.
jsi j = j1 ali j 1
It can be obtained as below by solving the above inequality 0\a\ max (di )
lim (1)1 Ek = xyT
k!‘
Convergence proof of GCTS
ð8Þ
T
where E(k) = (e1 (k), e2 (k), . . . , eN (k)) , I is the unit matrix, L is the Laplace matrix corresponding to the topological map, and P can be regarded as a state transition matrix. In order to prove the convergence of GCTS, it is necessary to prove the equation lim (E(k)) = lim (Pk E(0)), k!‘
k!‘
so the problem is transformed into the existence of the limit lim (Pk ). From equation (8), it can be known that k!‘
1 T 3 1N 3 1 3 (1N 3 1 ) , k 2 N N
ð15Þ
Then lim (E(k)) =
k!‘
E(k + 1) = (I aL)E(k) = PE(k)
ð14Þ
Within this range, 1 = s1 .s2 sN .0. Therefore, it can be inferred that only s1 = 1 and the rest is 0. According to Perron theorem,15 equation (15) can be deduced
=
The control protocol of the discrete system composed of N nodes is written as a matrix expression, as shown in equation (8)
ð13Þ
1 3 1N 3 1 3 (1N 3 1 )T E(0) N
ð16Þ
So, the final value converges to the mean of the time offset of the initial state.
Simulations In order to verify the performance of the algorithm based on GCTS, this article uses MATLAB simulation software to simulate the algorithm.
Simulation under different network topology types P = I aL
ð9Þ
where the characteristic value of L can be recorded as 0 = l1 \l2 lN . Then, matrix L can be transformed as below JL = SL LSL1
ð10Þ
Similarly, the matrix P can be changed into equation (11) JP = SL PSL1 = SL (I aL)SL1 = I aJL
ð11Þ
According to equations (9)–(11), it can be obtained k
lim (P ) = lim
k!‘
k!‘
k (SL JP SL1 )
= lim (SL JPk SL1 )
First, GCTS is performed in four kinds of small networks composed of eight nodes shown in Figure 4. They are linear network (L8), ring network (R8), star network (S8), and full-connection Network (F8). These small networks are really part of a large network. At the same time, in order to measure the degree of synchronization of the whole network, a synchronization criterion called Kuramoto is introduced.14 Kuramoto is shown as follows r=
N X 1 3 3 exp ( j2pf) exp ( j2pfi ) N i=1
is the instantaneous average phase of all nodes where f and fi is the instantaneous phase of node i. r = 1 when ð12Þ
k!‘
SL and SL1 are fixed matrices, as long as lim (JPk ) can k!‘ converge to a value. The diagonal elements of JL are 0 = l1 \l2 lN \ max (di ). Then, the corresponding diagonal elements of JP are si = 1 ali . The next step is to ensure that lim (ski ) is convergent. k!‘
Equation (13) could be obtained as follows
ð17Þ
Figure 4. Four different network topologies.
6
International Journal of Distributed Sensor Networks
Figure 6. Distribution of nodes in the simulation. Figure 5. Kuramoto parameters of nodes in different network topologies based on GCTS.
Table 1. Simulation parameters. Network scale Number Transmission radius Simulation times Frequency offset distribution Phase offset distribution Algorithm
100 m 3 100 m 100 10 m 100 N (1, 105 ) U (1, 1000) ATS, CCS, TSMA, GCTS
ATS: Average Time Synchronization Protocol; CCS: Consensus Clock Synchronization; TSMA: Time Synchronization using Max and Average Consensus Protocol; GCTS: Group Consensus Time Synchronization Protocol.
Specific simulation conditions are shown in Table 1. 100 nodes are distributed in 100 m 3 100 m region (Figure 6). The transmission radius of the node is 10 m and the number of simulation operations is 100. The node’s clock frequency of crystal oscillator is 32.768 kHz. Therefore 1 tick = 1/32,768 Hz = 30.5 ms, and the initial frequency shift of the node obeys normal distribution with mean 1 and the variance 105 . The initial clock offset of the node obeys the uniform distribution from 0 to 1000 ticks.
Simulation of a round iteration
all nodes are in the same phase. At this point, the network can achieve synchronization. As shown in Figure 5, all nodes under different network topology types can achieve synchronization but the convergence speed is different. In these four networks, the number of iterations required to achieve steady-state synchronization from fast to slow is all-connected, star-shaped, ring-shaped, and linear networks. Because in a fully connected network each node can exchange information with all other nodes, convergence speed is the fastest. In the star network, information can also be exchanged with neighbor nodes through the central node, so the convergence speed is secondary. In the ring network and linear network, information can only be passed through multi-hop to the whole network; meanwhile, there is additional time information transmission process, so the convergence rate is slow.
Before and after running the operation based on group consistency algorithm, average time error area distribution of each node in the network is shown in Figure 7. According to Figure 7(a), it can be seen that the maximum relative time error of the nodes in the network can be up to 500 ticks before the group consistency algorithm is implemented. From Figure 7(b), the maximum relative error of nodes in the network is reduced to less than 100 ticks. Moreover, the time error distribution after running the group consistency algorithm is more stable. A statistical histogram of the average error of all nodes before and after the operation is shown in Figure 8. From Figure 8(a), the standard deviation of the time error is 293.33 before implementing group consistency algorithm, while the standard deviation is 43.69 after the operation, as shown in Figure 8(b). So, after running the group consistency algorithm, the time synchronization performance of the network is improved obviously.
Simulation in large-scale WSNs
Performance simulation of several algorithms
This article analyzes time synchronization in large-scale WSNs.
Convergence rate is an important factor to measure the performance of the distributed time synchronization
Jin et al.
7
Figure 7. Synchronization error (a) before and (b) after one round of group consensus algorithm.
Figure 8. Histogram of the synchronization errors (a) before and (b) after one round of group consensus algorithm.
algorithm. In order to study the convergence rate of GCTS algorithm, it is simulated under the same environment with other classical time synchronization algorithms like ATS, CCS, and TSMA. The relationship between rounds and Kuramoto parameter of all nodes’ phase through the corresponding iterations is shown in Figure 9.
The Kuramoto parameter represents the convergence degree of the current state of the whole network. The closer to 1 Kuramoto parameter of the node phase, the higher the degree of convergence under this condition. From Figure 9, it can be seen that the performance of GCTS is better than ATS, CCS, and TSMA. After each round operation, Kuramoto parameter of GCTS
8
Figure 9. Kuramoto parameters of nodes after running ATS, CCS, TSMA, and GCTS.
is closer to 1 than ATS, CCS, and TSMA. This is because GCTS exchanges information between that node with all the neighbor nodes, while ATS, CCS, and TSMA algorithms are carried out between two nodes. In order to analyze the change in the phase of each node in the network more clearly, the relationship
International Journal of Distributed Sensor Networks between the phase error of each node and the number of iterations for each node after running ATS, CCS, TSMA, and GCTS algorithms is shown in Figure 10. It can be seen from the figure that the convergence rate of GCTS is the fastest. In summary, GCTS has better convergence performance. The plot in Figure 11 clearly shows the relationship between the scale of network and the convergence rate. The synchronization speed of the GCTS, CCS, and TSMA algorithms can still maintain a steady state with the size of the network increasing. For distributed algorithms, although the number of nodes in the network is increasing, the number of neighbor nodes increases as well as the range of iteration when the propagation distance is constant. However, the synchronization speed of GCTS is significantly faster than that of CCS and TSMA since GCTS algorithm employs group mechanism instead of point-to-point communication mechanism. Therefore, in large-scale WSNs, the performance of GCTS is better. The condition above is for network within a fixed region. However, if the average node density and transmission range remain same, the partial synchronization is still the same as before. As the network size increases, the number of iterations increases and the convergence speed of each algorithm will decrease correspondingly.
Figure 10. Clock errors of 100 simulated nodes running ATS, CCS, TSMA, and GCTS.
Jin et al.
9 radio, the convergence of the whole network is affected little. In stage D, after replacing these dead nodes, due to no clock correction, these nodes cannot be synchronized. When these nodes begin to perform the GCTS algorithm, all the nodes are synchronized soon.
Conclusion
Figure 11. Average number of rounds of CCS, TSMA, and GCTS in networks of different sizes.
In this article, a new algorithm for clock synchronization of the large-scale WSNs is proposed based on the group consensus protocol. In the protocol, the main idea is to use the group information exchange mechanism to replace the original point-to-point mechanism. At the same time, it takes propagation delay into consideration and compensates for the clock frequency shift together with the clock phase shift. The advantages of the protocol are fully distributed, scalable, and it has simple calculation, fast convergence rate, and low energy consumption. Finally, a set of simulations is presented to show the good performance of the proposed protocol. Declaration of conflicting interests
Figure 12. The performance of the algorithms in dynamic topology network.
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
Simulation in the dynamic topology network This article also analyzes convergence performance under the dynamic topology network. When the node dies and is replaced, the convergence of the network and the stability of the clock error after implementing the GCTS algorithm are observed in Figure 12. The experimental process is divided into four stages: stage A as the beginning phase, stage B as the stable phase, stage C as the death phase, and stage D as the recovery phase. In stage A, as a result of just running GCTS algorithm, the clock phase offset of each node is different initially. So, through about 30 rounds of iterations, the entire network enters the synchronization phase. In stage B, 3% of nodes closed randomly, and once the node is open, it will begin running GCTS algorithm. Because most of the nodes have been in the synchronization phase, new nodes added have little influence on the whole event adjustment. At the same time, it can also be seen from Figure 12 that nodes can quickly enter the synchronization state. In stage C, 5% of nodes are dead due to the longtime work. When a small number of nodes turn off their
The author(s) received no financial support for the research, authorship, and/or publication of this article.
References 1. Sichitiu ML. Cross-layer scheduling for power efficiency in wireless sensor networks. In: Proceedings of the twenty-third annual joint conference of the IEEE computer and communications societies, Hong Kong, China, 7–11 March 2004, pp.1740–1750. New York: IEEE. 2. Cadell P and Upcroft B. Distributed data acquisition unit with microsecond-accurate wireless clock synchronisation. In: Proceedings of the 2013 IEEE eighth international conference on intelligent sensors, sensor networks and information processing, Melbourne, VIC, Australia, 2–5 April 2013, pp.117–122. New York: IEEE. 3. Dengchang Z, Zhulin A and Yongjun X. Time synchronization in wireless sensor networks using max and average consensus protocol. Int J Distrib Sens N 2013; 8: 388–391. 4. Cena G, Scanzio S, Valenzano A, et al. Implementation and evaluation of the reference broadcast infrastructure synchronization protocol. IEEE T Ind Inform 2015; 11: 801–811. 5. Seth G and Harisha A. Energy efficient timing-sync protocol for sensor network. In: Proceedings of the 2015
10
6.
7.
8.
9.
International Journal of Distributed Sensor Networks international conference on computing and network communications (CoCoNet), Trivandrum, India, 16–19 December 2015, pp.912–916. New York: IEEE. Huang DJ, You KJ and Teng WC. Secured flooding time synchronization protocol. In: Proceedings of the 2011 IEEE eighth international conference on mobile adhoc and sensor systems, Valencia, 17–22 October 2011, pp.620–625. New York: IEEE. Sun Y, Jiang Q and Zhang K. A clustering scheme for Reachback Firefly Synchronicity in wireless sensor networks. In: Proceedings of the 3rd IEEE international conference on network infrastructure and digital content, Beijing, China, 21–23 September 2012, pp.27–31. New York: IEEE. Solis R, Borkar VS and Kumar PR. A new distributed time synchronization protocol for multihop wireless networks. In: Proceedings of the 45th IEEE conference on decision and control, San Diego, CA, 13–15 December 2006, pp.2734–2739. New York: IEEE. Schenato L and Fiorentin F. Average TimeSynch: a consensus-based protocol for clock synchronization in wireless sensor networks. Automatica 2011; 47(9): 1878–1886.
10. Maggs MK, O’Keefe SG and Thiel DV. Consensus clock synchronization for wireless sensor networks. IEEE Sens J 2012; 12(6): 2269–2277. 11. Wu J, Bai Y and Zhang L. Distributed time synchronization in wireless sensor networks via second-order consensus algorithms. T Tianjin Univ 2015; 21: 113–121. 12. Leng M and Wu YC. Distributed clock synchronization for wireless sensor networks using belief propagation. IEEE T Signal Proces 2011; 59(11): 5404–5414. 13. Du J and Wu YC. Distributed clock skew and offset estimation in wireless sensor networks: asynchronous algorithm and convergence analysis. IEEE T Wirel Commun 2013; 12(11): 5908–5917. 14. Noh KL, Chaudhari QM, Serpedin E, et al. Novel clock phase offset and skew estimation using two-way timing message exchanges for wireless sensor networks. IEEE T Commun 2007; 55(4): 766–777. 15. Matcovschi MH and Pastravanu O. Perron-Frobenius theorem and invariant sets in linear systems dynamics. In: Proceedings of the 2007 Mediterranean conference on control and automation, Athens, 27–29 June 2007, pp.1–6. New York: IEEE.
