Research Article  Open Access
Adaptive FractionalOrder Operator Clock Synchronization Algorithm Based on Grey Prediction
Abstract
As a result of the influence of clock drift and uncertainty delay in synchronous message transmission, the clock synchronization model based on statistical distribution cannot accurately describe clock deviation. This model also requires a large number of timestamp samples that cause a storage occupation issue for wireless sensor nodes with limited resources. The modeling method based on grey prediction has advantages of low sample demand and simple modeling process. However, the accuracy of the existing clock synchronization models needs to be improved. Based on the grey prediction theory, this paper proposes an adaptive fractionalorder operator clock synchronization algorithm considering uncertainty delay. First, based on the clock model and clock offset model, the frequency offset between nodes is optimized by taking the mean on the clock frequencies. Second, a grey prediction algorithm based on a fractionalorder operator is proposed by estimating the uncertainty delay in message transmission to obtain the clock offset. Finally, the order of the fractionalorder accumulation is adjusted adaptively in the grey prediction model according to the collected timestamp sample values so that the estimation of the uncertainty delay is more accurate, thereby improving the accuracy of the clock offset. Compared with the firstorder accumulative grey prediction clock synchronization algorithms and timingsync protocol for sensor networks, the proposed scheme improved the synchronization accuracy by 29.18% and 44.01%, respectively, and reduced the variance of the clock offset by 48.66% and 64.89%. Thus, the proposed algorithm is characterized by improved stability.
1. Introduction
Clock synchronization between different nodes in wireless sensor networks (WSNs) is an important basis for coordinating network operations [1]. Data fusion, energy management, transmission scheduling, and other processes require sensor nodes to be clock synchronized [2, 3]. However, sensor nodes tend to have limited energy. Under the premise of ensuring a certain clock synchronization accuracy, reducing the energy consumption of the nodes during the clock synchronization process and extending the life cycle of the nodes are the primary considerations in designing the clock synchronization algorithm [4]. Clock synchronization usually relies on the way of message transmission between the nodes. However, the time delay is generated in the stage of message transmission between the nodes, and the delay can be usually divided into two parts: certainty delay and uncertainty delay. As the former can be directly measured, the latter is one of the key factors that affect the accuracy of synchronization. Therefore, the influence of uncertainty delay on synchronization accuracy in message transmission can be reduced by establishing an accurate prediction model [5, 6].
Existing clock synchronization methods can be generally classified into two types according to the transmission message interaction mode: senderreceiver synchronization (SRP) model and receiverreceiver synchronization (RRP) model. The SRP model can be divided into a twoway message exchange synchronization mechanism and a oneway broadcast synchronization mechanism according to the message exchange manner in the synchronization process. In the algorithm based on the twoway message exchange synchronization mechanism, two sensor nodes are synchronized with each other through a twoway message exchange. The message has a local timestamp value, and the node calculates the clock offset by four collected timestamps, which has the advantage of high synchronization precision. However, in this process, frequent message exchanges between two nodes cause an increase in energy consumption. In the algorithm based on the oneway broadcast synchronization mechanism, the reference node broadcasts a data packet including its timestamp. The node to be synchronized calculates the data communication delay and performs delay compensation as soon as they receive a timestamp message, and then completes the clock synchronization with the reference node. The advantage is that the algorithm has low complexity and energy consumption. However, the certainty delay is regarded as a part of the clock offset; thus, the synchronization precision is lower than the synchronization mechanism of the twoway message exchange. On the other hand, such clock synchronization algorithm depends on the broadcast node, and the intranetwork synchronization fails if only the broadcast node fails. In the timingsync protocol for sensor networks (TPSNs) based on the SRP model [7], the realization of clock synchronization with the help of a tree structure in which each leaf node can establish a clock synchronization with its parent node, and the average delay of the message, is calculated by using a twoway message exchange. The synchronization precision is relatively high, but the treebased structure easily results in limited robustness and scalability. Especially when a node dies or a new node needs to be added, the tree structure has to be rebuilt, thereby adding extra overhead. In the RRP model, the sensor node synchronizes with the reference node by listening to the twoway synchronization message. The reference broadcast synchronization protocol [8] is a typical method based on the RRP model. Generally, the reference node broadcasts the message, and the other nodes receive the message by comparing the two timestamps that two nodes received the same data packet, which can reduce the influence of the uncertainty delay in the message transmission and improve the synchronization accuracy. However, due to a large amount of message exchange, the largescale sensor network consumes a large amount of energy. Above all, the clock synchronization algorithm in this study is implemented based on the SRP model.
Furthermore, each sensor node has its clock while being equipped with a crystal oscillator. In the actual environment, the oscillator frequency is unstable due to the oscillator quality, node power battery, temperature, humidity, and device aging. These issues often cause clock frequency offset, and the clock synchronization deviation becomes increasingly larger over time. Maintaining accuracy requires frequent synchronization to correct the clock offset, which also causes energy waste. Thus, the effect of clock frequency offset needs to be considered in clock synchronization deviation. To reduce the energy consumption created by the clock synchronization process, extend the network life cycle, and improve the accuracy of clock synchronization, this study proposes an adaptive fractionalorder operator clock synchronization algorithm based on the grey prediction GM(1,1) model [9]. The main contributions are as follows: (i)Based on the GM(1,1) grey prediction model, the grey prediction clock synchronization algorithm using a fractional accumulative operator is proposed. The accuracy of clock synchronization is improved in comparison with the grey prediction clock synchronization algorithm of the first order(ii)The study uses the collected timestamp sample values to adaptively obtain the best fractionalorder cumulative generation order. When the timestamp received by the node is significantly different from the historical record, the proposed algorithm can be adaptive to adjust the cumulative generation order to fit the current timestamp sequence, thereby improving the accuracy of the estimation of the uncertainty delay(iii)Compared with the entire network synchronization based on the spanning tree algorithm, the proposed algorithm does not require a lot of energy to maintain the network structure. This method should be robust to the dynamic changes in the network topology; furthermore, the applicable scene is wider
The subsequent parts of this paper are organized as follows. Section 2 provides an indepth study of the proposed clock synchronization algorithm. Section 3 describes the series of models used in this paper. Section 4 details the grey prediction clock synchronization algorithm based on fractionalorder operators. Section 5 evaluates and compares the performance between our proposed algorithm and other existing algorithms through simulation experiments. Section 6 summarizes the entire study.
2. Related Work
Clock synchronization is essentially related to the prediction of clock parameters, and the accuracy of the synchronization algorithm depends on the statistical properties of the parameter estimation algorithm. At present, many clock parameter estimation algorithms based on a twoway message exchange mechanism have been proposed to solve the accuracy problem in clock synchronization [10–12]. In the algorithm of clock synchronization from the perspective of statistical signal processing, the uncertainty delay in the message transmission process is modeled as a probability density function model [13]. Under the assumption that the uncertainty delay obeys the exponential distribution, Jeske in [14] deduces the closed solution of the clock offset maximum likelihood estimator, and the derivation process assumes that the means of the certainty delay and exponential distribution are both unknown. In [15], under the premise of symmetric and asymmetric exponential delay distributions, the minimum variance undeviating estimator of the clock offset and the best linear undeviating estimator using the order statistic are derived. However, if only the clock offset is estimated, then the frequency offset of the clock becomes increasingly larger to maintain synchronization between the nodes. The clock offset needs to be corrected frequently, thereby increasing the energy consumption of the node. To extend the synchronization cycle, Chaudhari et al. in [16] deduced the joint maximum likelihood estimation of the clock offset and clock frequency offset under the exponential delay model. Noh et al. in [17] assumes that the certainty delay is known and the Gaussian delay model is applied, and then derives the joint maximum likelihood estimation of the clock offset and clock frequency offset, as well as the corresponding Cramer lower bound. In [18], under the case of unknown certainty delay and Gaussian delay model, three different clock offset and clock frequency offset estimation methods are proposed. In [19], the unknown certainty delay is regarded as a redundant parameter; based on the exponential delay model, the approximate Cramer’s lower bound of the joint maximum likelihood estimation is derived. Nevertheless, the performance of the maximum likelihood estimation depends on the type of random delay model, and its estimation accuracy on the Gaussian distribution model is relatively better than that of the exponential distribution model. At the same time, under the linear clock model, the clock parameters obtained by the maximum likelihood estimation are only valid for shortterm applications; in the nonlinear clock model, the search algorithm is required and the algorithm complexity is relatively high. In [20], an adaptive powerefficient time synchronization algorithm is proposed for mobile underwater sensor networks; the study takes into account the timevarying clock frequency offset and combines the DopplerEnhanced synchronization protocol (DESync) and the Kalman filter tracking the clock frequency offset to achieve time synchronization; the study reduces the energy consumption of the nodes and guarantees accurate time synchronization. In [21], the clock parameter estimation problem under the exponential model is transformed into a linear programming problem through which the geometric analysis of the feasible domain, a low power consumption estimator is proposed. Owing to the characteristic of slow variation in clock parameters, the clock parameter estimated value obtained in the previous synchronization cycle can be used for prediction in the next cycle; thus, the clock parameter tracking technique has recently received attention. Assuming that the clock frequency offset and clock offset can be directly observed under the influence of noise, and the uncertainty delay obeys the Gaussian distribution, Freris et al. in [22] uses the Kalman filter to track the clock frequency offset. In [23], the Kalman filter is also employed to track both the clock frequency offset and clock offset. In [24, 25], the factorgraph model is applied to derive the method of tracking clock offset under the hypothesis exponential distribution delay model. However, the probability density function of the uncertainty delay model is unknown in advance, and frequent changes in the application environment of the wireless sensor network result in no certain statistical distribution model that can simulate the uncertainty delay in the message transmission process [26–28]. The literature [29] implemented a time synchronization service for lowpower WSN using low frequency realtime clocks in each node; what is more, three algorithms based on different strategies have been proposed, the research fulfilled a tradeoff between the desired level of synchronization and the associated energy consumption. AbdelGhaffar in [30] also proves that the clock parameter estimation algorithm is highly sensitive to the choice of the delay model. Therefore, studying the clock parameter estimation method under any delay distribution model is necessary. The aforementioned algorithms do not study the clock synchronization under any delay distribution model.
The preceding analysis shows that only the clock offset correction without considering the clock frequency offset must be increasing the synchronization overhead, while the synchronization algorithm based on statistical distribution corrects the clock frequency offset. However, the delay distribution depends on the specific environment; thus, the adaptability is poor and requires a large number of timestamp samples. Energy consumption is large as well. Based on the treestructured clock synchronization protocol, when the network topology changes, additional energy is needed to maintain the network topology, especially when the reference node fails, the network structure needs to be reconstructed, and the energy consumption is relatively large. Wu et al. in [31] proposed a clock synchronization scheme based on grey prediction, which estimates the clock offset based on the grey prediction model of GM(1,1); the scheme requires a few timestamp samples and does not need to hypothesize that the uncertainty delay obeys the specific distribution model, which effectively solves the problem of relying on a specific statistical distribution model in the estimation of clock synchronization parameters. However, this scheme does not correct the clock frequency offset nor does it have a specific network structure to synchronize the entire network. In addition, modeling the timestamp data based on the firstorder grey accumulation sequence is not the optimal order for modeling. From the perspective of statistical signal processing, this paper proposes an adaptive grey prediction clock synchronization algorithm with fractionalorder operators, which can estimate the clock offset of a node without a specific delay distribution model. The proposed algorithm calculates the timestamp samples collected to determine the current best fractionalorder accumulator operator that is applied to the next clock offset estimation process, thereby improving the synchronization accuracy. Compared with the spanning treebased wholenetwork clock synchronization scheme [32], the proposed synchronization scheme in this paper based on the breadthfirst search of the graph does not need to consider the failure and rotation of the root node, to a certain extent, which prolongs the life cycle of the entire network.
3. System Models
3.1. Network Model
Based on the assumption that the sensor network consists of nodes, the entire network structure can be represented by an undirected graph , where represents the set of nodes in the network, and represents the edge set. If an edge exists, then communication between nodes and can be conducted directly. A node located within the broadcast communication range of the node is called a neighbor node of the node , and is represented by . Furthermore, the following assumptions are made: (i)The network communication link is reliable and no message packet loss occurs(ii)Each node is static and has a unique ID(iii)Each node has its own clock and is equipped with a crystal oscillator(iv)Any two nodes can communicate through a limited hop
3.2. Clock Model
Ideally, the clock of a sensor node is defined as , where represents the ideal time [33]. Usually, the clock function of the th node is
Among them, the parameter is the clock offset and is the clock frequency offset.
The clock model is shown in Figure 1. The three oblique lines indicate the fast clock of node , slow clock of node , and ideal clock. Equation (2) defines the clock relationship between nodes and .
Among the lines, and represent the clock offset and clock frequency offset between the two nodes, respectively. When and , the clocks of node and are synchronized.
3.2.1. Clock Offset Model
Assume that nodes and perform clock synchronization, and the clock offset model of twoway time message exchange is as shown in Figure 2.
As shown in Figure 2, during the timestamp exchange of the th round, node transmits its th timestamp at time , and then records the local time as . Node at time receives the timestamp sent from node at time and records the local time as , similarly. Node then sends its th timestamp at time , recording the local time as . At time , node receives the timestamp sent from node at time and records the local time as . Assume that node is the primary synchronization node and node is synchronized to node . Combined with formula (2), the local clocks of nodes and at time are as follows: where and are the clock frequency offsets of nodes and , respectively; and are the clock offsets of nodes and , respectively; and is the real time.
According to Figure 2, the relationship between the four time points of sending and receiving timestamps of nodes and and the transmission delay can be obtained as follows: where and are the propagation delays between nodes and , is the certainty delay and and are uncertainty delays. Combined with formula (3) and (4), the time when nodes and receive the data packet for the th time can be expressed as
4. Adaptive FractionalOrder Clock Synchronization Algorithm Based on Grey Prediction Theory
The clock synchronization algorithm in this study consists of four parts: graph structure establishment based on breadthfirst search, clock frequency offset correction, clock synchronization prediction based on fractionalorder operator, and adaptive synchronization process. In the establishment period of the graph structure, the sensor network is constructed as an undirected graph, all nodes are equally connected in a breadthfirst search manner in the network, and they communicate only with their neighbor nodes. In the clock frequency offset correction period, the nodes broadcast their respective clock frequencies, and correct the clock frequency offset among the nodes by taking the mean to the clock frequencies. In the clock offset synchronization period of the fractionalorder operator, the grey prediction of the fractionalorder operator is modeled for the clock offset so that the uncertainty delay is estimated. Finally, the estimated value of the clock offset is obtained. In the adaptive synchronization period, whether to optimize the fractionalorder operator is determined according to the clock offset. When an order needs to be optimized, the optimal order is calculated according to the previously received timestamp samples.
4.1. BreadthFirstWide NetworkWide Clock Synchronization
To improve the robustness of the clock synchronization network structure, this paper proposes a graphbased breadthfirst search synchronization algorithm. All nodes in the network are equal, do not need to maintain a specific network structure, and perform the same synchronization algorithm. Furthermore, no root node or gateway node exists. As all nodes only communicate with their neighbor nodes, the algorithm should be robust to the dynamic changes of the topology, and the scalability improves as the network scale increases. Compared with the spanning treebased clock synchronization protocol, the proposed algorithm consumes less energy and is more suitable for WSNs with frequent changes in topology. The network structure of the breadthfirst search based on the graph is shown in Figure 3.
sensor nodes are randomly laid on the monitored area, where each node has a different ID number and the same transmission range. The nodes in the same node transmission range are neighbor nodes. One node is randomly selected in the initial stage of the network as the origin node, which initiates the broadcast data packet and connects its neighbor nodes. In the process of connecting its neighbor nodes, the node with the small ID is connected before the node with the large ID. The process is repeated until all the nodes in the network are connected. All nodes in the network only communicate with their neighbor nodes, and nonneighbor nodes can communicate with each other through multihop transmission. As each node is peertopeer, all nodes can synchronize to the entire network by synchronizing with neighbor nodes.
4.2. Correction of Clock Frequency Offset
If the clock of the node is not corrected to its frequency offset, the clock offset increases over time. The clock frequency offset is a random variable that changes with time. It is not only affected by many factors but also difficult to model. Usually, in an extremely short period, the change in clock frequency offset is negligible, but the clock frequency offset needs to be corrected for longterm synchronization. Therefore, this study adopts a frequency offset optimization strategy to achieve frequency offset correction by exchanging and computing the clock frequency information between sensor nodes. Assume that in a WSN consisting of nodes, the clock frequencies of the nodes are, respectively, denoted as where . The specific frequency offset optimization algorithm is as follows:
Step 1. For each node in the connected network, , at time , , respectively, broadcast their own clock frequency information ω_{j}(n), .
Step 2. Node in the broadcast domain of node receives frequency information , , from its neighbor nodes, where is the total number of messages received. The node performs arithmetic mean on the received frequency information and then takes this mean as the frequency of the next time of the node , namely,
Step 3. Node broadcasts the updated frequency information to its neighbor node.
Step 4. The preceding process is repeated until the clock frequencies of all nodes in the network reach the same value, that is, the entire network reaches frequency synchronization.
The frequency offset optimization process is presented in Figure 4. Node 1 collects all clock frequency information from surrounding nodes for a period of time , computes the clock frequency of the entire group, and then sends the computed clock frequency value to its neighbor nodes 2, 3, and 4. Finally, the neighbor nodes update the clock frequency at the next moment according to the received frequency value.
4.3. Modeling of Clock Offset Based on FractionalOrder Operator
The grey prediction model [34] is based on the accumulation of small data of uncertainty on the offset of the cumulative generation of the original sequence. The application of the buffer operator data transformation technique [35] ensures the validity and rationality of the model prediction results. The modeling process is relatively simple and does not need to consider the probability distribution rule and statistical characteristics of the sample data. For the clock synchronization of wireless sensor network nodes with limited resources, the delay of information transmission is uncertain and the sample size is small, which is highly suitable for grey estimation theory for offset estimation.
As the frequency offset of the clock has been optimized in Section 2, when performing the clock synchronization prediction, this study only considers the offset of the clock, that is, the clock frequency offset is set to 1, which only estimates the clock offset. In this case, the clock offset between nodes and is . Thus, formula (5) can be rewritten as
The existing grey prediction model is mainly based on firstorder cumulative generation sequence modeling, and then the prediction value is obtained by firstorder cumulative reduction. However, the modeling based on the firstorder accumulation generation sequence is unsuitable for modeling data or experimental scenes of all features. Therefore, this study establishes a grey prediction model based on fractionalorder operator [36] without firstorder accumulation of sequences but adaptively finds the optimal order suitable for synchronization to improve the accuracy of grey prediction. At the same time, only a small sample is used to estimate the transmission delay of the clock information to achieve clock synchronization, thereby improving the application range of the algorithm. Thus, formula (7) can be rewritten as where , , is the certainty delay in message transmission, and and are the uncertainty delays in message transmission. Formula (8) shows that and are related to the delay of message transmission. Suppose that and are predicted values obtained by grey prediction, and and represent the minimum values of and , respectively.
From [14], the maximum likelihood estimation of the certainty delay and uncertainty delay , is where and represent the average of the samples and , respectively ( is the number of samples). The method of maximum likelihood estimation requires a large number of samples to obtain a meaningful mean, whereas the grey prediction method can be modeled with only a small sample, and can also be used for estimation. The method can be used in clock synchronization to reduce energy consumption. In this study, and are replaced by the estimated values obtained by grey prediction. The certainty delay and uncertainty delay , of the grey prediction can be obtained as follows: where is the number of samples used in calculation, and the grey model requires at least 4 samples. As increasing the number of samples does not necessarily improve the prediction accuracy [37], is set to 4 when considering the node energy consumption. The following is the step of estimating and using the grey model of the fractionalorder operator.
From [37], when is the order accumulation that generates operators of , where
From equation (8), , . are the difference values of the timestamp samples. The algorithm for estimating using fractionalorder operators is as follows:
Step 1. Initialize the order of the fractionalorder accumulative operator.
Step 2. Introduce the grey weakening buffer operator to reduce the randomness of the sample sequence. Let the original timestamp sample sequence be , and the disposed by the operator have a sequence that is recorded as , and be the average weakening buffer operator [38], among them Using the constructed average weakening buffer operator to work on the original data sequence smoothens the sequence and considers the principle of new information prioritization. The latest information remains unchanged under the action of the average weakening buffer operator, that is, when in formula (12), , the randomness of the sample sequence is reduced, thereby improving the accuracy of the prediction model [38].
Step 3. Perform order accumulation on the grey sequence of the sample sequence to obtain , among them
Step 4. Obtain the estimated value of using the fractionalorder accumulative operator grey model (where ) The parameters in equation (14) can be estimated by least squares as follows: where and are, respectively (where ), The expression for is
Step 5. Accumulate the fractionalorder ( order) to obtain the restored value of , that is, the estimated value of the th time of (where ): Similarly, the estimated value of the th time of (where ) is obtained as follows: When the estimated values of and are substituted into equation (10) with the deterministic delay and uncertainty delay , the estimated value of [ ].
Step 6. The ratio of the message transmission delay in the node communication link is According to formulas (8) and (20), the clock offsets of nodes and is obtained After computing the clock offset between nodes and , node can synchronize to node according to its current local clock time minus . As each node is peertopeer, the algorithms executed are also the same, and all nodes can achieve fullnetwork synchronization by synchronizing with neighboring nodes.
4.4. Synchronization Process of Adaptive Clock Offset
Due to changes in environmental factors, the timestamp value sent by the node after a period of time may change. Thus, setting a threshold of clock synchronization offset is necessary, where is equal to the mean of the clock offset obtained by the fractionalorder cumulative grey prediction synchronization algorithm. Based on the assumption that nodes and are clocksynchronized, the clock offset of their th synchronization is , where is the number of times that nodes and have been synchronized. The expression of threshold is
When the remaining energy of the node meets the requirement of one synchronization, we have to check whether the current clock offset is greater than a given threshold. If the remaining energy is greater than the threshold, then the node calculates an optimal fractional cumulative generation order according to the currently obtained timestamp value as a new fractionalorder cumulative generation order in the synchronization model. The clock offset is estimated using the new order when performing a new round of clock offset predictions.
A flowchart of the clock offset correction and adaptive clock offset synchronization process of the fractionalorder operator is shown in Figure 5.
5. Simulation Results
5.1. Simulation Setup
In this paper, the proposed fractionalorder grey prediction algorithm is simulated by MATLAB_R2017b and compared with the firstorder grey prediction algorithm [33] and TPSN [7] algorithm to analyze the performance. The simulation scene is set as follows: 100 nodes are randomly deployed in an area of . The nodes in the area are distributed randomly. In order to simulate the uncertainty of network environment, the value of the uncertain delay, the packet sending time and the channel access time are selected randomly in a fixed interval [4 μs, 55 μs], respectively. Set the number of synchronizations per node to 500, and each synchronization cycle lasts 10 s because this duration is sufficient to synchronize all nodes in the network [17]. Each message propagation of the three algorithms is conducted by broadcasting to ensure fairness, and the communication channel condition is supposed to be ideal. The details of the algorithm parameter settings are listed in Table 1. The three parameters that clock synchronization offset, clock offset variance, and distribution of the clock offset are evaluated for the three algorithms, and the performance differences between the firstorder accumulative operator and fractionalorder operator in the grey prediction are compared in this section.

5.2. Performance Evaluation
The variance of the clock offset obtained by the firstorder cumulative grey prediction and fractionalorder cumulative grey prediction based on the same timestamp data is shown in Figure 6.
The green circle represents the variance of the firstorder accumulated clock synchronization offset, the blue solid line indicates the variance of the fractionalorder accumulated clock synchronization offset, and the red circle marks the position of the minimum offset variance obtained by the fractionalorder cumulative grey prediction. In the synchronization process, the deviation based on the 3.6order cumulative grey prediction clock synchronization offset is the smallest. Thus, this study selects the 3.6order fractionalorder to perform the grey accumulation prediction.
In addition, since the fractionalorder also contains the integerorder (the cumulative order is in the firstorder, the variance of the clock offset obtained by the fractional grey prediction and the firstorder grey prediction is the same), so the grey prediction is compared with the firstorder. The grey prediction of fractionalorder accumulation is more applicable. As shown in Figure 6, the result shows that the variance of the offset obtained by some fractional grey predictions are large, indicating that the accumulated order at that time is not suitable for grey prediction.
Figure 7(a) shows the clock offset change between two nodes with clock synchronization for 500 times. Since the fractionalorder operator is used to improve the adaptability of the algorithm and the clock frequency offset is optimized, the clock synchronization offset based on the fractionalorder grey prediction is closer to 0 and smaller than using the firstorder grey prediction algorithm and the TPSN algorithm. Since the TPSN algorithm does not estimate the uncertainty delay, the offset of the synchronization is large. Thus, it can be seen from the figure that the clock synchronization algorithm of the fractionalorder operator converges faster and can achieve synchronization quickly.
(a)
(b)
Figure 7(b) shows the clock offset change between the two nodes with clock synchronization for 100 times. It can be seen that the fractionalorder cumulative grey prediction algorithm obtains a smaller range of clock offset fluctuations and is more accurate. The average value of the clock offset is about 8.98 μs. The average value of the clock offset obtained by the firstorder cumulative grey prediction algorithm is about 12.68 μs, and the average value of the clock offset obtained by the TPSN algorithm is about 16.04 μs. Compared with the firstorder accumulative grey prediction algorithm, the fractionalorder cumulative grey prediction algorithm reduces the average offset by about 29.18%. Compared with the TPSN algorithm, the fractionalorder cumulative grey prediction algorithm reduces the average offset by about 44.01%.
Figure 8 shows the variance of the firstorder cumulative grey prediction algorithm, the fractionalorder accumulation grey prediction algorithm, and the TPSN algorithm at 500 times clock synchronizations. Because the algorithm proposed in this paper optimizes the clock frequency offset and adaptively adjusts the cumulative order of the grey prediction model, the obtained offset has smaller variance and better stability. Relative to the firstorder accumulative grey prediction algorithm, the fractionalorder cumulative grey prediction algorithm reduces the variance by about 48.66%. For the TPSN algorithm, the fractionalorder cumulative grey prediction algorithm reduces the variance by about 64.89%.
Figure 9 depicts the clock offset distribution for the three algorithms. Since the clock frequency offset is optimized and a reasonable fractional order is used, the clock offset based on the fractionalorder grey prediction clock synchronization algorithm is close to 0, the synchronization precision is higher, and the clock offset distribution range is more centralized, most of which are concentrated in between 8 μs and 12 μs, thereby suggesting that the algorithm is more stable than that of the other two algorithms. Meanwhile, the clock offset distribution of the firstorder cumulative grey prediction clock synchronization algorithm and the TPSN algorithm is more dispersed, and the stability of the algorithm is worse than that of the proposed algorithm.
5.3. Discussion
The stability of clock synchronization is an important criterion for evaluating a clock synchronization algorithm. Based on the simulation results in Figures 7(a) and 7(b), compared with the other two algorithms, the clock synchronization offset obtained by the algorithm proposed in this paper is the most stable. The average value of the clock synchronization offset is about 9 μs. In the process of clock synchronization proposed in this paper, when the clock synchronization offset is greater than the threshold, the adaptive mechanism is triggered to stabilize the clock offset in a small range. The experimental results in Figure 9 show that the synchronization accuracy of the proposed algorithm is the highest among the three algorithms. The main reasons for the stability of clock offset are that the proposed algorithm optimizes the clock frequency offset and selects the best fractionalorder accumulation order, which makes fluctuation of the clock offset small. The clock offset distribution of the firstorder grey prediction algorithm is relatively concentrated as the firstorder grey prediction algorithm does not adjust the accumulation order. The clock offset distribution of TPSN is scattered, and the offset range of TPSN is larger than that of fractionalorder cumulative grey prediction algorithm and firstorder cumulative grey prediction algorithm. Meanwhile, the offset of TPSN is very unstable as the clock offset estimating process is absent; thus, the clock correction of each node is not accurate, and the clock frequency offset is also not being optimized in TPSN.
6. Conclusions
To improve the clock synchronization accuracy, this study estimates the uncertainty delay in clock synchronization from the perspective of statistical signal processing and proposes an adaptive fractionalorder operator clock synchronization algorithm based on the grey prediction model. First, the frequency offset between nodes is optimized by taking the mean on the clock frequency so that the sensor nodes do not need to be frequently synchronized, thereby reducing the energy consumption of the nodes to a certain extent. Second, after frequency offset optimization, a fractionalorder operator clock model based on grey prediction theory is proposed to estimate the uncertainty delay in synchronous message transmission and then an accurate clock offset is obtained. In addition, the algorithm adaptively adjusts the order of the fractionalorder accumulation in the grey prediction model according to the collected timestamp sample values, thereby improving the accuracy of the uncertainty delay estimation. At the same time, the graphbased breadthfirst search is more suitable than the spanning tree to construct the clock synchronization network structure for the wireless sensor network with frequent changes in topology, which does not require a large amount of energy to maintain the topology network while achieving clock synchronization of the entire network. Compared with the existing firstorder cumulative grey prediction and TPSN algorithms, the proposed algorithm reduces the mean on the clock offset by approximately 29.18% and 44.01%, respectively. The proposed algorithm also reduces the variance of the offset by approximately 48.66% and 64.89%. The result of the distribution of the clock deviation shows that the proposed algorithm has higher synchronization accuracy and stability than the firstorder accumulated grey prediction and TPSN algorithms. In our further research work, we will testify the performance of the proposed algorithm in the real IoT platform and optimize the communication performance of the proposed clock synchronization algorithm; the characteristics of the communication channel will be considered simultaneously.
Data Availability
The test data, simulation data, and the proposed method used to support the findings of this study are available from the corresponding author upon request. The proposed algorithm is in the state of patentpending; thus, the access of the data used in the proposed algorithm is restricted now.
Conflicts of Interest
The authors declare no conflict of interest.
Acknowledgments
This research was funded in part by the National Key Research and Development Program (Grant No. 2017YFB1401800), the Jilin Province Foundation for Excellent Youth Talents, grant number 20190103052JH, and the Project of the Education Department of Jilin Province, grant number JJKH20200802KJ.
References
 M. Elsharief, M. A. Abd elGawad, and H. W. Kim, “LowPower scheduling for time synchronization protocols in a wireless sensor networks,” IEEE Sensors Letters, vol. 3, no. 4, article 7500104, 2019. View at: Publisher Site  Google Scholar
 S. Lu, Z. Wang, Z. Wang, and S. Zhou, “Throughput of underwater wireless ad hoc networks with random access: a physical layer perspective,” IEEE Transactions on Wireless Communications, vol. 14, no. 11, pp. 6257–6268, 2015. View at: Publisher Site  Google Scholar
 T. Qiu, Y. Zhang, D. Qiao, X. Zhang, M. L. Wymore, and A. K. Sangaiah, “A robust time synchronization scheme for industrial Internet of Things,” IEEE Transactions on Industrial Informatics, vol. 14, no. 8, pp. 3570–3580, 2018. View at: Publisher Site  Google Scholar
 F. Sivrikaya and B. Yener, “Time synchronization in sensor networks: a survey,” IEEE Network, vol. 18, no. 4, pp. 45–50, 2004. View at: Publisher Site  Google Scholar
 D. L. Mills, “Internet time synchronization: the network time protocol,” IEEE Transactions on Communications, vol. 39, no. 10, pp. 1482–1493, 1991. View at: Publisher Site  Google Scholar
 S. Ganeriwal, R. Kumar, and M. B. Srivastava, “Timingsync protocol for sensor networks,” in Proceedings of the first international conference on Embedded networked sensor systems  SenSys '03, pp. 138–149, Los Angeles, CA, USA, November 2003. View at: Publisher Site  Google Scholar
 A. B. Kulakli and K. Erciyes, “Time synchronization algorithms based on timingsync protocol in wireless sensor networks,” in 2008 23rd International Symposium on Computer and Information Sciencess, pp. 27–29, Istanbul, Turkey, October 2008. View at: Publisher Site  Google Scholar
 J. Elson, L. Girod, and D. Estrin, “Finegrained network time synchronization using reference broadcasts,” in Proceedings of the 5th symposium on Operating systems design and implementation  OSDI '02, pp. 147–163, Boston, MA, USA, December 2002. View at: Publisher Site  Google Scholar
 J. Guo, X. Xiao, J. Liu, and S. Mao, “Stability of GM(1,1) power model on vector transformation,” Journal of Systems Engineering and Electronics, vol. 26, no. 1, pp. 103–109, 2015. View at: Publisher Site  Google Scholar
 T. Qiu, X. Liu, M. Han, M. Li, and Y. Zhang, “SRTS : a selfrecoverable time synchronization for sensor networks of healthcare IoT,” Computer Networks, vol. 129, no. 2, pp. 481–492, 2017. View at: Publisher Site  Google Scholar
 T. Qiu, L. Chi, W. Guo, and Y. Zhang, “STETS: a novel energyefficient time synchronization scheme based on embedded networking devices,” Microprocessors and Microsystems, vol. 39, no. 8, pp. 1285–1295, 2015. View at: Publisher Site  Google Scholar
 K. S. Yildirim and A. Kantarci, “Time synchronization based on slowflooding in wireless sensor networks,” IEEE Transactions on Parallel and Distributed Systems, vol. 25, no. 1, pp. 244–253, 2014. View at: Publisher Site  Google Scholar
 K.L. Noh, E. Serpedin, and K. Qaraqe, “A New approach for time synchronization in wireless sensor networks: pairwise broadcast synchronization,” IEEE Transactions on Wireless Communications, vol. 7, no. 9, pp. 3318–3322, 2008. View at: Publisher Site  Google Scholar
 D. R. Jeske, “On maximumlikelihood estimation of clock offset,” IEEE Transactions on Communications, vol. 53, no. 1, pp. 5354, 2005. View at: Publisher Site  Google Scholar
 Q. M. Chaudhari, E. Serpedin, and K. Qaraqe, “On minimum variance unbiased estimation of clock offset in a twoway message exchange mechanism,” IEEE Transactions on Information Theory, vol. 56, no. 6, pp. 2893–2904, 2010. View at: Publisher Site  Google Scholar
 Q. M. Chaudhari, E. Serpedin, and K. Qaraqe, “On maximum likelihood estimation of clock offset and skew in networks with exponential delays,” IEEE Transactions on Signal Processing, vol. 56, no. 4, pp. 1685–1697, 2008. View at: Publisher Site  Google Scholar
 K. L. Noh, Q. M. Chaudhari, E. Serpedin, and B. W. Suter, “Novel clock phase offset and skew estimation using twoway timing message exchanges for wireless sensor networks,” IEEE Transactions on Communications, vol. 55, no. 4, pp. 766–777, 2007. View at: Publisher Site  Google Scholar
 M. Leng and Y.C. Wu, “On clock synchronization algorithms for wireless sensor networks under unknown delay,” IEEE Transactions on Vehicular Technology, vol. 59, no. 1, pp. 182–190, 2010. View at: Publisher Site  Google Scholar
 M. Leng and Y. C. Wu, “On joint synchronization of clock offset and skew for wireless sensor networks under exponential delay,” in Proceedings of 2010 IEEE International Symposium on Circuits and Systems, Paris, France, May 2010. View at: Publisher Site  Google Scholar
 F. Zhou, Q. Wang, G. Han, G. Qiao, Z. Sun, and A. Niaz, “APESync: an adaptive power efficient time synchronization for Mobile underwater sensor networks,” IEEE Access, vol. 7, pp. 52379–52389, 2019. View at: Publisher Site  Google Scholar
 M. Leng and Y. C. Wu, “Lowcomplexity maximumlikelihood estimator for clock synchronization of wireless sensor nodes under exponential delays,” IEEE Transactions on Signal Processing, vol. 59, no. 10, pp. 4860–4870, 2011. View at: Publisher Site  Google Scholar
 N. M. Freris, V. S. Borkar, and P. R. Kumar, “A modelbased approach to clock synchronization,” in Proceedings of the 48th IEEE Conference on Decision and Control held jointly with 2009 28th Chinese control conference, pp. 15–18, Shanghai, China, December 2009. View at: Google Scholar
 H. Kim, X. Ma, and B. R. Hamilton, “Tracking lowprecision clocks with timevarying drifts using Kalman filtering,” IEEE/ACM Transactions on Networking, vol. 20, no. 1, pp. 257–270, 2012. View at: Publisher Site  Google Scholar
 A. Ahmad, D. Zennaro, and E. Serpedin, “Timevarying clock offset estimation in twoway timing message exchange in wireless sensor networks using factor graphs,” in 2012 IEEE International Conference on Acoustics, Speech and Signal Processing, pp. 25–30, Kyoto, Japan, March 2012. View at: Google Scholar
 A. Ahmad, D. Zennaro, E. Serpedin, and L. Vangelista, “A factor graph approach to clock offset estimation in wireless sensor networks,” IEEE Transactions on Information Theory, vol. 58, no. 7, pp. 4244–4260, 2012. View at: Publisher Site  Google Scholar
 J. Heidemann, W. Ye, and J. Wills, “Research challenges and applications for underwater sensor networking,” in IEEE Wireless Communications and Networking Conference, pp. 03–06, Las Vegas, NV, USA, April 2006. View at: Google Scholar
 S. S. Kunniyur and S. Narasimhan, “Modelling the effect of network parameters on delay in wireless adhoc networks,” in 2005 Second annual IEEE communications society conference on sensor and ad hoc communications and networks, 2005. IEEE SECON 2005, pp. 26–29, Santa Clara, CA, USA, September 2005. View at: Publisher Site  Google Scholar
 A. A. Syed and J. Heidemann, “Time synchronization for high latency acoustic networks,” in Proceedings IEEE INFOCOM 2006. 25TH IEEE International Conference on Computer Communications, pp. 23–29, Barcelona, Spain, April 2006. View at: Publisher Site  Google Scholar
 L. T. Bruscato, T. Heimfarth, and E. P. de Freitas, “Enhancing time synchronization support in wireless sensor networks,” Sensors, vol. 17, no. 12, article 2956, 2017. View at: Publisher Site  Google Scholar
 H. S. AbdelGhaffar, “Analysis of synchronization algorithms with timeout control over networks with exponentially symmetric delays,” IEEE Transactions on Communications, vol. 50, no. 10, pp. 1652–1661, 2002. View at: Publisher Site  Google Scholar
 W. Chen, X. J. Gao, Z. H. Li, X. Jin, and L. Hong, “A fast clock synchronization method for wireless sensor networks based on grey prediction,” Adhoc & Sensor Wireless Networks, vol. 40, no. 1/2, pp. 145–161, 2018. View at: Google Scholar
 H. J. Li, B. Ding, C. M. Dong, and H. Zhou, “Clock synchronization network planning algorithm of electric power telecommunication system based on minimum spanning tree,” Journal of Electric Power, vol. 27, no. 2, pp. 136–139, 2012. View at: Google Scholar
 A. R. Swain and R. C. Hansdah, “A model for the classification and survey of clock synchronization protocols in WSNs,” Ad Hoc Networks, vol. 27, pp. 219–241, 2015. View at: Publisher Site  Google Scholar
 C.’. O.K. Chen and T.L. Tien, “The indirect measurement of tensile strength by the deterministic grey dynamic model DGDM(1, 1, 1),” International Journal of Systems Science, vol. 28, no. 7, pp. 683–690, 2007. View at: Publisher Site  Google Scholar
 C. I. Chen, H. L. Chen, and S. P. Chen, “Forecasting of foreign exchange rates of Taiwan’s major trading partners by novel nonlinear Grey Bernoulli model NGBM(1, 1),” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 6, pp. 1194–1204, 2008. View at: Publisher Site  Google Scholar
 W. Meng, D. Yang, and H. Huang, “Prediction of China’s sulfur dioxide emissions by discrete grey model with fractional order generation operators,” Complexity, vol. 2018, Article ID 8610679, 13 pages, 2018. View at: Publisher Site  Google Scholar
 H. C. Hsieh, H. R. Huang, and Y. Feng, “An approach to improve estimation performance of GM (1, 1) model,” Chinese Medicine, vol. 12, no. 1, pp. 249–253, 2008. View at: Google Scholar
 S. Liu, J. Forrest, and J. Yang, “A brief introduction to grey systems theory,” in Proceedings of 2011 IEEE International Conference on Grey Systems and Intelligent Services, pp. 15–18, Nanjing, China, September 2011. View at: Google Scholar
Copyright
Copyright © 2020 Hongwei Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.