CN102449940A

CN102449940A - 执行节点之间的时间同步的网络同步方法及装置

Info

Publication number: CN102449940A
Application number: CN2009801593976A
Authority: CN
Inventors: 杰弗里·M·加纳; 柳贤锡; 朴根柱; 李俊行
Original assignee: Samsung Electronics Co Ltd
Current assignee: Samsung Electronics Co Ltd
Priority date: 2009-05-18
Filing date: 2009-11-13
Publication date: 2012-05-09
Anticipated expiration: 2029-11-13
Also published as: EP2434673B1; KR20100124195A; US8467487B2; US20100290572A1; CN104468017B; KR101705592B1; WO2010134674A2; EP2434673A2; EP2434673A4; CN104468017A; CN102449940B; WO2010134674A3

Abstract

涉及执行节点之间的时间同步的网络同步方法及装置。提供如下的网络同步方法及装置，即，系统启动(start up)之后，若开始节点之间的同步，则通过逐渐增加用于测量传递时间(propagation time measurement)的窗口大小的方法，能够提高抖动(jitter)、漂移(wander)及时间同步的性能。而且，还提供如下的网络同步方法及装置，即，传递时间的测量执行至最大窗口大小之后，在传递时间的平均值计算中采用指数性的加权值，由此能够提高抖动、漂移及时间同步的性能。

Description

执行节点之间的时间同步的网络同步方法及装置

技术领域

技术领域涉及在网络中执行节点之间的时间同步的网络同步方法及装置。

背景技术

存在很多的用于利用时间戳(time stamps)传递时序(timing)的方案(schemes)。

在所有的这种方案中，构成网络的各个节点尝试与所述节点中的一个节点进行同步。此时，一个节点属于被称为主的节点，即主节点。而且，在这种方法中，时序能够在被称为最高主的一个节点，即最高主节点跟踪。

主节点向除了该主节点的主节点之外的附属在主节点的所有节点传送包含时间戳的消息。此时，时间戳可表示消息被传送的瞬间的时刻。各个从节点接收这种时间戳，并在时间戳上相加从节点与主节点之间的传递时间(propagation time)。即，从节点通过比较接收消息的时刻与相加传递时间的时间戳，即相加传递时间的主节点的传送时刻，计算主节点和从节点的偏移(offset)，从而能够进行同步。

发明内容

技术方案

提供如下网络同步方法及装置，即，系统启动(start up)之后，若开始节点之间的同步，则通过逐渐增加用于测量传递时间(propagation timemeasurement)的窗口大小的方法，能够提高抖动(jitter)、漂移(wander)及时间同步的性能。

而且，还提供如下的网络同步方法及装置，即，传递时间的测量执行至最大窗口大小之后，在传递时间的平均值计算中采用指数性的加权值，由此能够提高抖动、漂移及时间同步的性能。

本发明的一方面提供的网络同步方法包括如下步骤：在第k(k为自然数)次时间步长(time step)中，测量自从节点至主节点的传递时间；若所述k小于滑动窗口大小(window size)，则根据第一计算式计算在第k+1次时间步长下的传递时间平均值，若所述k大于或等于滑动窗口大小，则根据第二计算式计算在第k+1次时间步长下的传递时间平均值；以及，利用所述主节点上的消息传送时刻、所述从节点上的消息接收时刻以及所述传递时间平均值计算偏移。

此时，所述第一计算式及第二计算式可利用第k次时间步长下的传递时间平均值和在第k次时间步长下测量的传递时间计算出在所述第k+1时间步长下的传递时间平均值。

此时，所述第一计算式如下。

x_{k + 1} = \frac{{kx}_{k} + d_{k}}{k + 1}

其中，x_k是指在第k次时间步长下的传递时间平均值，d_k是指在第k次时间步长下测量的传递时间。

此时，所述第二计算式如下。

x_k+1＝ax_k+(1-a)d_k

其中，其中，x_k是指在第k次时间步长下的传递时间平均值，d_k是指在第k次时间步长下测量的传递时间，a是指一阶(first-order)数字滤波器的滤波器系数。

此时，随着所述滑动窗口大小增加，所述a可以指数性地减小。

此时，所述a可通过以下数学式计算。

a＝e^-1/M

在此，M表示滑动窗口大小。

此时，所述a可通过以下数学公式计算。

a＝e^-P/M

在此，M表示滑动窗口大小，P表示任意的常数。

技术效果

系统启动(start up)之后，若开始节点之间的同步，则通过逐渐增加用于测量传递时间(propagation time measurement)的窗口大小的方法，能够提高抖动(jitter)、漂移(wander)及时间同步的性能。

而且，传递时间的测量执行至最大窗口大小之后，在传递时间的平均值计算中采用指数性的加权值，由此能够提高抖动、漂移及时间同步的性能。

附图说明

图1为根据相关技术的计算传递时间的方法的一示例。

图2为根据相关技术的计算传递时间的方法的另一示例。

图3示出根据本发明的一实施例的传递时间测量方法。

图4为表示图3的330步骤的详细过程的图。

图5表示示出根据本发明的一实施例的网络装置的构成。

图6表示根据本发明的一实施例的用于提高测量传递时间的准确度的方法。

图7至图9表示根据本发明的一实施例的用于设定参数的信令。

具体实施方式

以下，参照附图来详细说明本发明的实施例。

图1为根据相关技术的计算传递时间的方法的一示例。

主节点101给从节点102传送消息。此时，作为传送所述消息的瞬间103的时刻的时间戳T1被记录在消息中。而且，当从节点接收消息时，从节点确认作为接收消息的瞬间104的时刻的接收时刻T2。若从节点102知道自主节点101的传递时间D1，则从节点102能够如同“偏移＝T2-T1-D1”计算出作为与主节点101之间的差的偏移，并可利用偏移执行与主节点101的同步。

但是，如同上述的假设，为了计算偏移，从节点102需要知道自主节点101的传递时间D1。即，从节点102需要计算传递时间D1。

图2为根据相关技术的计算传递时间的方法的另一示例。

分别附属于两个节点，即主节点201及从节点202的时钟(clock)之间的传递时间可由交换时间戳的所述时钟测量。作为所述时钟当中的一个的请求器(requester)在该请求器将第一消息传送至其他时钟时，在第一消息中记录该第一消息的传送时刻T1 203。作为所述另一时钟的应答器(responder)记录所述第一消息的接收时刻T2 204。

然后，所述应答器向所述请求器传送对应于所述第一消息的第二消息。此时，所述请求器在所述第二消息中记录所述第二消息的传送时刻T3 205而进行传送。最后，所述请求器确认所述第二消息的接收时刻T4 206。

所述请求器可利用如上确认的所述T1 203至所述T4 206计算如以下数学式2表示的所述传递时间。

【数学式1】

传递时间＝((T2-T1)+(T4-T3))/2

在此，假设测量所述传递时间的方法满足以下(1)至(4)的条件。

(1)所述传递时间是对称性的。即，两个方向的传递时间都相同。

(2)所述传递时间是稳定的。即，用于得出作为四个时间戳的所述T1至所述T4的时标(timescale)恒定。

(3)主节点和从节点之间的频率差小。

(4)所述时间戳的值准确。

但是，所述(1)及(2)的假设意味着所述请求器与所述应答器之间不能存在易变延迟(variable delay)的根源(source)。换言之，意味着所述主节点与从节点之间不能存在网络开关、网络桥等中继装置。任意这种装置的存在可导致抖动(Jitter)、漂移(Wander)及时间同步的不准确(inaccuracy)。

如上所述的主节点与从节点之间存在所述网络开关或所述网络桥，由此发生传递时间的变化及非对称性，作为用于调整传递时间的变化及非对称性的其他近似方法(approach)有在诸如所述网络开关或所述网络桥的中继装置中将各个消息的接收时刻及传送时刻记录在所述消息中的方法。这种所述方法要求在所述中继装置中包括被称为透明时钟(transparent clock)的时钟和记录时间的硬件(a clock and time stamping hardware)。

图3示出根据本发明的一实施例的传递时间测量方法。

根据本发明一实施例的传递时间测量方法计算测量值的平均值，并将计算出的平均值用作用于估计实际传递时间值的D_ave，k。据此，在本发明的一实施例中，可如同“偏移＝T2-T1-D_ave，k”计算出偏移。

图3的传递时间测量方法可在作为从设备运行的网络装置中执行。以下，将在网络中作为从设备运行的网络装置简单地称为从节点。相同地，将作为主设备运行的网络装置简单地称为主节点。

在步骤310中，从节点测量自从节点至主节点的传递时间(propagationtime)。此时，自从节点至主节点的传递时间可以是自从节点至主节点的传递时间T2-T1及自主节点至从节点的传递时间T4-T3的和或平均值。此时，对于传递时间的测量可以以预设定的时间单位测量k(k为自然数)次。因此，在第k次时间单位测量传递时间可应用为与在第k次时间步长(time step)下测量传递时间相同的意思。即，从节点在第k次时间步长下测量至主节点的传递时间。

在步骤320中，从节点计算在第k+1次时间步长下的传递时间平均值。在步骤320中，若所述k小于滑动窗口大小(window size)，则从节点依据第一计算式计算在第k+1次时间步长下的传递时间平均值，若所述k大于或等于滑动窗口大小，则从节点依据第二计算式计算在第k+1次时间步长下的传递时间平均值。此时，“滑动窗口”可以简单地表示为“窗口”。此时，滑动窗口大小M(M为自然数)可以是测量的传递时间被保存的次数。滑动窗口大小M与配备与从节点的缓存器的大小相关。

此时，如图4所示，步骤320可包括步骤401、步骤403、步骤405而执行。

在步骤401中，从节点可判断作为之前的时间步长的k是否小于滑动窗口大小M。此时，若k＜M，则从节点执行步骤403，若k大于或等于M，则执行步骤405。

在步骤403中，从节点可依据第一计算式计算第k+1次时间步长下的传递时间平均值。

在步骤405中，从节点可依据第二计算式计算第k+1次时间步长下的传递时间平均值。所述第一计算式及第二计算式利用第k次时间步长下的传递时间平均值及在第k次时间步长下测量的传递时间计算所述第k+1次时间步长下的传递时间平均值。

第一计算式与“系统启动(star up)之后，节点之间的同步开始进行时，逐渐增加用于测量传递时间(propagation time measurement)的窗口大小”有关。第二计算式与“传递时间的测量执行至最大窗口大小之后，在传递时间平均值中应用指数性的加权值”有关。

第一计算式可如以下的数学式2定义。

【数学式2】

x_{k + 1} = \frac{{kx}_{k} + d_{k}}{k + 1}

在此，x_k表示在第k次时间步长下的传递时间平均值，d_k表示第k次时间步长下测量的传递时间。此时，为了数学式2的计算，从节点可以配备保存(save)在第k时间步长下的传递时间平均值x_k或kx_k的单元。

第二计算式可如以下的数学式3定义。

【数学式3】

x_k+1＝ax_k+(1-a)d_k

在此，x_k表示在第k次时间步长下的传递时间平均值，d_k表示第k次时间步长下测量的传递时间，a表示一阶(first-order)数字滤波器的滤波器系数。为了数学式3的计算，从节点可以配备用于保存(save)在第k时间步长下的传递时间平均值x_k的单元。所述数学式3可依据一阶数字滤波器计算，此时数字滤波器可以是无限脉冲响应(IIR，infinite impulse response)滤波器。

在数学式3中，a可被设定为随着滑动窗口大小的增加而指数性地减少。因此，a可通过以下数学式4或数学式5计算。

【数学式4】

a＝e^-1/M

【数学式5】

a＝e^-P/M

在此，P表示任意的常数。所述数学式5表示将数学式4一般化的公式。

滑动窗口大小M应选择为大于改变在第k次时间步长下测量的传递时间d_k的样品数的值。而且，滑动窗口大小M可选择为大于消息的实际传递时间(actual propagation time)的变量(variation)的时标(time scale)的值。在某种状况下，M可被设定为1000，但并不局限于此。

在步骤330中，从节点利用主节点的消息传送时刻、所述从节点的消息接收时刻以及所述传递时间平均值计算偏移。偏移可如同“偏移＝T2-T1-D_ave，k”计算。此时，从节点和主节点分别计算传递时间平均值并维持。

图5示出根据本发明的一实施例的网络装置的构成。

图5所示的网络装置500可作为从节点而运行。

网络装置500包括传递时间测量单元510、传递时间平均值计算单元520及偏移计算单元530。

传递时间测量单元510在第k(k为自然数)次时间步长下测量至主节点的传递时间。传递时间测量单元510可利用时间戳测量至主节点的传递时间。传递时间测量单元510为了测量当前时间步长下的传递时间，可将作为第一消息的传送时刻的第一时间戳记录于所述第一消息并传送给所述主节点。此时，所述主节点可将所述第一时间戳、作为所述第一消息的接收时刻的所述第二时间戳及作为对应于所述第一消息的第二消息的传送时刻的第三时间戳记录于所述第二消息而传送给网络装置500。传递时间测量单元510确认作为所接收的所述第二消息的接收时刻的第四时间戳。传递时间测量单元510可利用所述第一时间戳及所述第四时间戳测量当前时间步长下的传递时间。

传递时间平均值计算单元520计算在第k+1次时间步长下的传递时间平均值。传递时间平均值计算单元520可包括用于维持M个传递时间测量值的缓存器。若所述k小于滑动窗口大小M，则传递时间平均值计算单元520依据第一计算式计算在第k+1次时间步长下的传递时间平均值。若所述k大于或等于滑动窗口大小M，则传递时间平均值计算单元520依据第二计算式计算在所述第k+1次时间步长下的传递时间平均值。

偏移计算单元530利用所述主节点的消息传送时刻、所述消息的接收时刻以及所述传递时间平均值计算偏移。偏移可如同“偏移＝T2-T1-D_ave，k”计算。

<传递时间测量的精度提高方法>

根据本发明的一实施例，可利用传递时间平均值提高测量传递时间的精度。此时，用于网络同步的时间精度可根据在相关节点中使用的应用中的要求最高时间精度的应用所决定。例如，要求图像和音响之间的时间同步的应用要求比决定用于移动通信的频率值的应用更低的时间精度。

图6示出根据本发明的一实施例的用于提高测量传递时间的准确度的方法。

图6所示的例中，从最高主设备610至从节点640的中继段(hop)数量为N-1个，全部的节点数为N个。此时，在第一中继节点620及第二中继节点630之间可存在其他中继节点，第一中继节点620及第二中继节点630可以是时间-感知系统(time-aware system)。

在S601步骤中，第二中继节点630自从节点640接收第一消息。此时，第一消息为请求与时间精度相关的性能参数(performance parameter)信息的消息。第一消息可包含表示中继段的数量的字段n，此时表示中继段的数量的字段n的值为1。

由于第二中继节点630不是最高主设备，因此在S602步骤中向第一中继节点620传送第二消息。此时，第二消息是请求与时间精度相关的性能参数信息的消息。第二消息可包含表示中继段的数量的字段n，此时表示中继段的数量的字段n的值为1。

由于第一中继节点620不是最高主设备，因此在S603步骤中向最高主设备610传送第三消息。第三消息可包含表示中继段的数量的字段n′，此时表示中继段的数量的字段n′的值是n′＝n+1。

第一消息、第二消息及第三消息可分别包含相关节点的性能参数信息。此时，相关节点的性能参数信息可包括数字滤波器类型、数字滤波器大小、滑动窗口大小及本地时钟(local clock)性能中的至少某一个。此时，第二消息及第三消息还可以包含之前节点的性能参数。即，第二消息可包含从节点640的性能参数信息及第二中继节点630的性能参数信息。同样地，第三消息可包含从节点640的性能参数信息、第二中继节点630的性能参数信息及第一中继节点620的性能参数信息。

由此，将S601步骤至S603步骤一般化为在任意的中继节点执行的动作而进行说明时如下。将任意的中继节点假设为第N-2中继节点时，第N-2中继节点从第N-1节点接收性能参数请求信息，并判断第N-2中继节点本身是否为最高主节点。若第N-2次中继节点判断为本身不是最高节点时，将从第N-1节点接收的性能参数请求消息传递给第N-3节点。此时，从第N-1节点接收的性能参数请求消息中包含第N-1节点的性能参数信息。此时，传递到第N-3节点的性能参数请求消息包含第N-1节点的次性能参数信息及第N-2节点的性能参数信息。

最高主节点610接收到第三消息时，生成性能表格。此时，性能表格中可包含关于N-2个中继节点的性能参数信息及从节点640的性能参数信息。最高节点610利用所生成的性能表格生成包含所有节点的性能参数信息的性能信息包，并返还性能信息包。性能信息包通过S604步骤至S606步骤中被传递到从节点640。

此时，性能信息包可在各个节点被更新之后传递到下一个节点。例如，第一中继节点620可对包含于性能信息包中的自身的信息更新之后传递至第二中继节点630。

另外，从节点640或最高主设备650可参照性能信息包或性能表格向各个中继节点请求性能参数变更。即，从节点640或最高主设备650可向各个中继节点请求增加滑动窗口大小或改变滤波器常数。例如，若当前正在使用的应用是对时间精度要求高的应用，则从节点640可向中继节点请求增加滑动窗口大小或变更滤波器常数。

S601步骤至S606步骤可在图3的330步骤之后执行。而且，各个节点还可以在测量传递时间之前执行S601步骤至S606步骤，或可以周期性或非周期性地执行S601步骤至S606步骤。各个节点通过周期性或非周期性地执行S601步骤至S606步骤，能够提高测量传递时间的精度。

<用于设定参数的信令>

图7至图9示出根据本发明的一实施例的用户设定参数的信令。图7至图9所示的信令可应用于从节点640或最高主设备650向各个中继节点请求变更性能参数。而且，图7至图9所示的信令还可以应用于设定用于计算传递时间平均值的参数。此时，用于计算传递时间平均值的参数例如可以是窗口大小。另外，图7至图9所示的信令还可以应用于设定网络所需的任意的参数。图7至图9中，时间-感知桥(time-aware bridge)2，、3、4、5及6可执行中继节点的功能。

此时，信令可通过传送信令消息而被执行。信令消息传递时间-感知系统之间的信息、请求(requests)及/或命令(commands)。此时，信令消息可依据一个或一个以上的类型长度值(TLV：type-length-value)传递信息、请求及/或命令。

图7所示的例表示在时间-感知终端站(time-aware end-station)1中请求设定参数的情形。此时，时间-感知终端站1可执行从节点的功能。时间-感知终端站1将具有请求TLV的第一信令消息传送至作为相邻节点的时间-感知桥5。此时，第一信令消息可以是请求利用包含于请求TLV的值设定用于计算传递时间平均值的参数。

时间-感知桥5通过主端口(master port)接收第一信令消息，在保存关于请求(request)TLV的消息之后，通过第二信令消息转发请求(request)TLV。此时，关于请求TLV的消息例如可以是TLV类型及第一信令消息的序列号。此时，保存关于请求(request)TLV的原因在于，当设定参数的请求被拒绝时，使请求-拒绝(request-denied)TLV仅通过转发请求TLV的路径的逆方向传递。

如图7所示，请求TLV可通过时间-感知桥(time-aware bridge)2被转发至最高主设备时间-感知桥(grandmaster time-aware bridge)1。此时，最高主设备时间-感知桥1可执行与最高主设备650相同的功能。最高主设备时间-感知桥1同样保存关于请求TLV的信息。最高主设备时间-感知桥1可以对参数设定请求决定允许(grant)或拒绝(denial)。此时，对于参数设定请求的允许或拒绝可基于性能参数信息决定。

图8所示的例表示最高主设备时间-感知桥1拒绝参数设定请求的情形。图8中，请求-拒绝TLV是表示参数设定请求被拒绝的信令消息。

参照图8，当时间-感知桥2接收请求-拒绝TLV时，删除(remove)预先保存的关于请求TLV的信息，并转发请求-拒绝TLV。如图8所示，请求-拒绝TLV被传送至时间-感知终端站1。

图9所示的例表示最高主设备时间-感知桥1允许参数设定请求的情形。

当参数设定请求被允许时，最高主设备时间-感知桥1可利用包含于请求TLV的值设定用于计算传递时间平均值的参数。此时，最高主设备时间-感知桥1可删除预先保存的关于请求TLV的信息，并转发请求-允许(request-granted)TLV。此时，请求-允许TLV表示参数设定请求被允许。与请求-拒绝TLV不同，请求-允许TLV可通过最高主设备时间-感知桥1的所有主端口传递。由此，时间-感知桥2、3及4可接收请求-允许TLV。

时间-感知桥2、3及4利用包含于请求-允许TLV的值设定用于计算传递时间平均值的参数，并可通过主端口传递请求-允许TLV。此时，时间-感知桥2可删除预先保存的关于请求TLV的信息，并转发请求-允许TLV。当时间-感知终端站2接收到请求-允许TLV时，利用包含于请求-允许TLV的值设定用于计算传递时间平均值的参数。时间-感知终端站2不转发请求-允许TLV。

根据图7至图9中说明的信令，用于设定参数的信令消息可被传递至最高主设备。而且，参数设定请求被承认时，包含于参数设定请求的参数可被传递至网络内的所有节点。网络内的所有节点可设定包含于参数设定请求中的参数。如果参数设定请求被拒绝，则参数设定拒绝消息只能通过传递参数设定请求的路径(path)传递。

根据本发明实施例的方法可被实现为通过多种计算机手段执行的程序命令形态，从而记录于计算机可读介质中。所述计算机可读介质可单独或组合包含程序命令、数据文件、数据结构等。记录于所述介质的程序命令可以是专为本发明特别设计和构成的程序命令，或可以是计算机软件领域的技术人员所公知而能够使用的程序命令。

如上说明的本发明虽然仅通过有限的实施例和附图进行了说明，但本发明并不局限于上述实施例，本发明所属技术领域的具有普通知识的技术人员可通过这种记载进行多种修改及变更。

因此，本发明的范围并不局限于已说明的实施例，应根据权利要求书及与该权利要求书等同的内容定义。

Claims

1.一种网络同步方法，包括以下步骤：

在第k次时间步长中，测量自从节点至主节点的传递时间，其中，k为自然数；

若所述k小于滑动窗口大小，则根据第一计算式计算在第k+1次时间步长下的传递时间平均值，若所述k大于或等于滑动窗口大小，则根据第二计算式计算在第k+1次时间步长下的传递时间平均值；以及

利用所述主节点上的消息传送时刻、所述从节点上的消息接收时刻以及所述传递时间平均值计算偏移。

2.根据权利要求1所述的网络同步方法，其中，所述第一计算式及所述第二计算式利用第k次时间步长下的传递时间平均值及第k次时间步长下测量的传递时间计算所述第k+1次时间步长下的传递时间平均值。

3.根据权利要求2所述的网络同步方法，其中，所述第一计算式为

x_{k + 1} = \frac{{kx}_{k} + d_{k}}{k + 1}

其中，x_k为第k次时间步长下的传递时间平均值，d_k为第k次时间步长下测量的传递时间。

4.根据权利要求2所述的网络同步方法，其中，所述第二计算式为

x_k+1＝ax_k+(1-a)d_k

其中，x_k为第k次时间步长下的传递时间平均值，d_k为第k次时间步长下测量的传递时间，a为1阶数字滤波器的滤波器系数。

5.根据权利要求4所述的网络同步方法，其中，随着所述滑动窗口大小增加，所述a指数性地减小。

6.根据权利要求4所述的网络同步方法，其中，所述a通过如下的计算式计算：

a＝e^-1/M

其中，M为滑动窗口大小。

7.根据权利要求4所述的网络同步方法，其中，所述a通过如下的计算式计算：

a＝e^-P/M

其中，M为滑动窗口大小，P为任意的常数。

8.根据权利要求1所述的网络同步方法，还包括步骤：向中继节点传送性能参数请求消息，并接收对应于所述性能参数请求消息的性能信息包。

9.根据权利要求8所述的网络同步方法，还包括步骤：基于所述性能信息包向所述中继节点请求变更性能参数。

10.根据权利要求8所述的网络同步方法，所述性能信息包包括数字滤波器类型、数字滤波器大小、滑动窗口大小及本地时钟性能中的至少一个。

11.根据权利要求1所述的网络同步方法，包括步骤：发出用于计算所述传递时间平均值的参数设定请求。

12.根据权利要求11所述的网络同步方法，还包括步骤：接收针对所述参数设定请求的允许或拒绝。

13.根据权利要求11所述的网络同步方法，若所述参数设定请求被允许，则所述参数设定请求传递至网络内的所有的节点。

14.一种网络同步装置，包括：

传递时间测量单元，在第k次时间步长中，测量自从节点至主节点的传递时间，其中，k为自然数；

传递时间平均值计算单元，若所述k小于滑动窗口大小，则根据第一计算式计算在第k+1次时间步长下的传递时间平均值，若所述k大于或等于滑动窗口大小，则根据第二计算式计算在第k+1次时间步长下的传递时间平均值；以及

偏移计算单元，利用所述主节点上的消息传送时刻、所述从节点上的消息接收时刻以及所述传递时间平均值计算偏移。

15.根据权利要求14所述的网络同步装置，其中，所述第一计算式及所述第二计算式利用第k次时间步长下的传递时间平均值及第k次时间步长下测量的传递时间计算所述第k+1次时间步长下的传递时间平均值。

16.根据权利要求14所述的网络同步装置，其中，所述第一计算式为

x_{k + 1} = \frac{{kx}_{k} + d_{k}}{k + 1}

17.根据权利要求14所述的网络同步装置，其中，所述第二计算式为

x_k+1＝ax_k+(1-a)d_k

18.根据权利要求17所述的网络同步装置，其中，随着所述滑动窗口大小增加，所述a指数性地减小。

19.根据权利要求17所述的网络同步装置，其中，所述a通过如下的计算式计算：

a＝e^-1/M

其中，M为滑动窗口大小。

20.根据权利要求17所述的网络同步装置，其中，所述a通过如下的计算式计算：

a＝e^-P/M

其中，M为滑动窗口大小，P为任意的常数。