CN105721095B - 一种变电站设备时钟同步的改进方法 - Google Patents
一种变电站设备时钟同步的改进方法 Download PDFInfo
- Publication number
- CN105721095B CN105721095B CN201610109863.1A CN201610109863A CN105721095B CN 105721095 B CN105721095 B CN 105721095B CN 201610109863 A CN201610109863 A CN 201610109863A CN 105721095 B CN105721095 B CN 105721095B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- clock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 11
- 230000005540 biological transmission Effects 0.000 claims abstract description 19
- 239000013078 crystal Substances 0.000 claims abstract description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000003111 delayed effect Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 abstract description 3
- 238000007619 statistical method Methods 0.000 abstract description 2
- 108700009949 PTP protocol Proteins 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0635—Clock or time synchronisation in a network
- H04J3/0638—Clock or time synchronisation among nodes; Internode synchronisation
- H04J3/0658—Clock or time synchronisation among packet nodes
- H04J3/0661—Clock or time synchronisation among packet nodes using timestamps
- H04J3/0667—Bidirectional timestamps, e.g. NTP or PTP for compensation of clock drift and for compensation of propagation delays
Abstract
本发明公开了一种变电站设备时钟同步的改进方法，步骤一：本地时钟的频偏计算，解得本地时钟频偏η的最大似然估计值；步骤二：本地时钟的相偏计算，计算N次对时相偏Φk的平均值作为相偏值Φ；步骤三：对本地时钟晶振固有频率进行补偿；步骤四：对本地时钟进行补偿。针对智能变电站网络负荷较重时，网络传输往返延时不对称，传统的精确时钟同步协议(Precision Time Protocol，PTP)单次对时在本地时钟的校正上误差较大，本发明提供的一种变电站设备时钟同步的改进方法，通过统计学方法在多次对时后计算出本地时钟频偏与相偏误差的最优解进而进行补偿。
Description
技术领域
本发明涉及一种变电站设备时钟同步的改进方法，属于变电站网络对时技术领域。
背景技术
现代化变电站以全站信息数字化、通信平台网络化、信息共享标准化为基本要求，自动完成信息采集、测量、控制、保护、监测和计量等基本功能。但无论测控装置、保护保护装置，还是计量装置的计算处理都要求采样数据应在同一个时间点上采集，以免相位和幅值产生误差。对于过电流保护等保护，因为微机保护的本地晶振时钟的短期稳定性非常高，不会对保护的动作精度造成影响。但对于差动保护和计量，由于微机保护本地晶振时钟并不十分准确，经过长时间的误差累积，会造成跨间隔间不同采样数据的相位误差和幅值误差的逐渐扩大，导致差动保护的误动作和计量的严重误差，因此需要对时技术对本地时钟进行校正。
IEEE1588定义了一个能在测量和控制系统中实现高精度时钟同步的协议——精确时钟同步协议(Precision Time Protocol，PTP)。PTP协议的基本原理是主时钟与本地时钟之间进行同步信息包的发送，对信息包的发出时间和接收时间信息进行记录，并且在每一条信息包上打上时间戳标签。PTP协议为了使用和描述时间信息，定义了四种多点传送的信息包，即同步信息包Sync、同步信息包之后的跟随信息包Follow_Up、延时测量信息包Delay_Req和延时测量应答信息包Delay_Resp。
同步信息包传递的机制成为“延时—请求响应机制”，如图1所示，主时钟周期性发送包含时钟质量的Sync信息包，紧接着发送Follow_up信息包告知本地时钟上个信息包的实际发送时间Tm1；本地时钟记录Sync信息包的到达时间Ts1，随后在Ts3时刻发送Delay_Req信息包；主时钟记录Delay_Req信息包的到达时间Tm3，并发送信息包Delay_Resp把Tm3告知本地时钟。本地时钟根据这四个时间信息计算出时钟的偏移和传输延迟。假设在时钟的信息通信中认为信息包的往返延迟是相等的，则图1中Tm1、Ts1、Ts3、Tm3四个点的连线是一个等腰梯形。
本地时钟便可计算出自身与主时钟之间的信息交换延时Tdelay为
Tdelay＝[(Tm3﹣Tm1)﹣(Ts3﹣Ts1)]/2 (1.1)
本地时钟与主时钟的时钟相位偏移Toffset为
Toffset＝[(Ts1﹣Tm1)+(Ts3﹣Tm3)]/2 (1.2)
本地时钟根据计算出来的相位偏移修改本地时间，从而达到与主时钟同步。
PTP协议通过这种“延时—请求响应机制”消除了物理层处理延时的影响，从而进一步提高了对时精度。虽然IEEE1588通过PTP协议消除了上层处理延时与物理层处理延时造成的影响，但它与所有网络对时算法一样存在局限性，PTP协议为计算时钟相偏所进行的四次测量都以对时信息的传输延迟在传输方向上相同为基础，但在实际使用中，不可能绝对满足这个前提。虽然IEEE1588协议引入了透明时钟和边界时钟，但需要交换机支持IEEE1588时钟机制，同时网络传输延时测量的随机性和较低频率，但当网络负载突发性变化时，造成信息包在传输方向上的延时差异变大，因此公式(1.1)与公式(1.2)成立前提的假设——信息包交换的往返延迟相等便不成立，图1中Tm1、Ts1、Ts3、Tm3四个点的连线将是一个不等腰梯形，若简单地通过公式(1.2)对本地时钟相偏进行计算，便会对本地时钟的校正带来较大误差。
该方法只对本地时钟的相位偏移进行了补偿，并未对本地时钟的频率偏移进行补偿。而造成相位偏移的根本原因是本地晶振时钟频率并不十分准确，与主时钟频率相比存在误差，经过长时间的误差累积而造成相位偏差。
发明内容
目的：为了克服现有技术中存在的不足，本发明提供一种变电站设备时钟同步的改进方法。
技术方案：为解决上述技术问题，本发明采用的技术方案为：
一种变电站设备时钟同步的改进方法，包括步骤如下：
步骤一：本地时钟的频偏计算，假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型；一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间；忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk；选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻；则 与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏；本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：
有
同理有T3,k＝(T4,k-d-Yk)(1+η)+Φ
(2.2)
由式(2.1)，用T2,N减去T2,1可得：
T2,N-T2,1＝(T1,N-T1,1+XN-X1)(1+η) (2.3)
同理，根据式(2.2)可得：
T3,N-T3,1＝(T4,N-T4,1+Y1-YN)(1+η) (2.4)
将第N个与第1个时间戳的差值分别定义为D1＝T1,N-T1,1，D2＝T2,N-T2,1，D3＝T3,N-T3,1，D4＝T4,N-T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN-X1与DY＝YN-Y1，则式(2.3)与(2.4)可分别改写为
D2＝(D1+DX)(1+η) (2.5)
D3＝(D4-DY)(1+η) (2.6)
假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2；因此DX与DY的联合概率密度函数为：
式(2.5)与(2.6)描述的简化观测模型的似然函数为：
令代入式(2.8)得：
通过使式(2.8)的似然函数L(η′,λ)取得最大值，则可以得到η′的最大似然估计值，为
式中为的从小到大顺序统计量，Ki为与δ(i)对应的D2或D3，当δ(1)＝δ1时，K1＝D2，K2＝D3；当δ(1)＝δ2时，K1＝D3，K2＝D2。
解得本地时钟频偏η的最大似然估计值为：
步骤二：本地时钟的相偏计算，相偏计算时需要先对本地时钟的观测时间T2,k、T3,k(k＝1…N)进行频率补偿，得到修正后的时间T’2,k、T’3,k：
T′2,k＝T2,k-(T1,k+d+Xk)ηMLE (3.1)
T′3,k＝T3,k-(T4,k-d-Yk)ηMLE (3.2)
由式(3.3)计算出N次对时每次的相偏Φk
Φk＝T′2,k+T′3,k-T1,k-T4,k (3.3)
计算N次对时相偏Φk的平均值作为相偏值Φ
步骤三：设本地时钟晶振固有频率为f，对其进行补偿后，实际频率f’为
f’＝f×(1+ηMLE) (2.12)
步骤四：设本地时钟此时输出时间为Ci(t)，对其进行补偿后，实际输出时间C’i(t)为
C’i(t)＝Ci(t)+Φ (3.5)
有益效果：针对智能变电站网络负荷较重时，网络传输往返延时不对称，传统的精确时钟同步协议(Precision Time Protocol，PTP)单次对时在本地时钟的校正上误差较大，本发明提供的一种变电站设备时钟同步的改进方法，通过统计学方法在多次对时后计算出本地时钟频偏与相偏误差的最优解进而进行补偿。
附图说明
图1为延时—请求响应机制原理示意图；
图2为时钟相偏与频偏的对时信息交换的时钟模型示意图；
图3为本地时钟频偏与相偏补偿流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
一、本地时钟的频偏计算
假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型，如图2所示。一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间。忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk。选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻。则与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏。本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：
有
同理有T3,k＝(T4,k-d-Yk)(1+η)+Φ
(2.2)
从式(2.1)与(2.2)中可以看出，受到频偏的影响，本地时钟的相对观测时间分别会受到式中(T1,k+d+Xk)η项与(T4,k-d-Yk)η项的影响。
从图2中可以看出，基于线性时钟模型由于受到频偏的影响，主时钟与本地时钟间的时间偏移不断增大，第N次信息交换与第1次交换时间戳的时间差异最大，因此可以通过使用第N个和第1个信息交换时的观测值，得到对于本地时钟频偏的类最大似然估计。
由式(2.1)，用T2,N减去T2,1可得：
T2,N-T2,1＝(T1,N-T1,1+XN-X1)(1+η) (2.3)
同理，根据式(2.2)可得：
T3,N-T3,1＝(T4,N-T4,1+Y1-YN)(1+η) (2.4)
将第N个与第1个时间戳的差值分别定义为D1＝T1,N-T1,1，D2＝T2,N-T2,1，D3＝T3,N-T3,1，D4＝T4,N-T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN-X1与DY＝YN-Y1，则式(2.3)与(2.4)可分别改写为
D2＝(D1+DX)(1+η) (2.5)
D3＝(D4-DY)(1+η) (2.6)
假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2。因此DX与DY的联合概率密度函数为：
式(2.5)与(2.6)描述的简化观测模型的似然函数为：
令代入式(2.8)得：
通过使式(2.8)的似然函数L(η′,λ)取得最大值，则可以得到η′的最大似然估计值，为
式中为的从小到大顺序统计量，Ki为与δ(i)对应的D2或D3，当δ(1)＝δ1时，K1＝D2，K2＝D3；当δ(1)＝δ2时，K1＝D3，K2＝D2。
解得本地时钟频偏η的最大似然估计值为：
设本地时钟晶振固有频率为f，对其进行补偿后，实际频率f’为
f’＝f×(1+ηMLE) (2.12)
该方法不需要考虑相偏的影响，只需通过第N次与第1次主时钟与本地时钟间信息交换的观测值即可快速算出，可简单快捷地对频率进行补偿频率补偿。
二、本地时钟的相偏计算
相偏计算时需要先对本地时钟的观测时间T2,k、T3,k(k＝1…N)进行频率补偿，得到修正后的时间T’2,k、T’3,k：
T′2,k＝T2,k-(T1,k+d+Xk)ηMLE (3.1)
T′3,k＝T3,k-(T4,k-d-Yk)ηMLE (3.2)
由式(3.3)计算出N次对时每次的相偏Φk
Φk＝T′2,k+T′3,k-T1,k-T4,k (3.3)
计算N次对时相偏Φk的平均值作为相偏值Φ
如图3所示，设本地时钟此时输出时间为Ci(t)，对其进行补偿后，实际输出时间C’i(t)为
C’i(t)＝Ci(t)+Φ (3.5)
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种变电站设备时钟同步的改进方法，其特征在于：包括步骤如下：
步骤一：本地时钟的频偏计算，假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型；一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间；忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk；选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻；则 与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏；本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>+</mo>
<mi>d</mi>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</msubsup>
<mo>+</mo>
<mi>d</mi>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&eta;</mi>
<mo>+</mo>
<msup>
<mi>&Phi;</mi>
<mi>r</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>d</mi>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&Phi;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.1</mn>
<mo>)</mo>
</mrow>
</mrow>
有
同理有T3,k＝(T4,k-d-Yk)(1+η)+Φ
(2.2)
由式(2.1)，用T2,N减去T2,1可得：
T2,N-T2,1＝(T1,N-T1,1+XN-X1)(1+η) (2.3)
同理，根据式(2.2)可得：
T3,N-T3,1＝(T4,N-T4,1+Y1-YN)(1+η) (2.4)
将第N个与第1个时间戳的差值分别定义为D1＝T1,N-T1,1，D2＝T2,N-T2,1，D3＝T3,N-T3,1，D4＝T4,N-T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN-X1与DY＝YN-Y1，则式(2.3)与(2.4)可分别改写为
D2＝(D1+DX)(1+η) (2.5)
D3＝(D4-DY)(1+η) (2.6)
假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2；因此DX与DY的联合概率密度函数为：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>X</mi>
</msub>
<mo>,</mo>
<msub>
<mi>D</mi>
<mi>Y</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>X</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>Y</mi>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.7</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2.5)与(2.6)描述的简化观测模型的似然函数为：
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.8</mn>
<mo>)</mo>
</mrow>
</mrow>
令代入式(2.8)得：
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<mo>|</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.9</mn>
<mo>)</mo>
</mrow>
</mrow>
通过使式(2.8)的似然函数L(η′,λ)取得最大值，则可以得到η′的最大似然估计值，为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<mo>|</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
</munder>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
</munder>
<mi>h</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中 为的从小到大顺序统计量，Ki为与δ(i)对应的D2或D3，当δ(1)＝δ1时，K1＝D2，K2＝D3；当δ(1)＝δ2时，K1＝D3，K2＝D2；
解得本地时钟频偏η的最大似然估计值为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>M</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>&eta;</mi>
<mo>&prime;</mo>
</msup>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<msub>
<mi>D</mi>
<mn>4</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>4</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mn>4</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2.11</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤二：本地时钟的相偏计算，相偏计算时需要先对本地时钟的观测时间T2,k、T3,k(k＝1…N)进行频率补偿，得到修正后的时间T’2,k、T’3,k：
T′2,k＝T2,k-(T1,k+d+Xk)ηMLE (3.1)
T′3,k＝T3,k-(T4,k-d-Yk)ηMLE (3.2)
由式(3.3)计算出N次对时每次的相偏Φk
Φk＝T′2,k+T′3,k-T1,k-T4,k (3.3)
计算N次对时相偏Φk的平均值作为相偏值Φ
<mrow>
<mi>&Phi;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>N</mi>
</munderover>
<msub>
<mi>&Phi;</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3.4</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤三：设本地时钟晶振固有频率为f，对其进行补偿后，实际频率f’为
f’＝f×(1+ηMLE) (2.12)
步骤四：设本地时钟此时输出时间为Ci(t)，对其进行补偿后，实际输出时间C’i(t)为
C’i(t)＝Ci(t)+Φ (3.5)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610109863.1A CN105721095B (zh) | 2016-02-26 | 2016-02-26 | 一种变电站设备时钟同步的改进方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610109863.1A CN105721095B (zh) | 2016-02-26 | 2016-02-26 | 一种变电站设备时钟同步的改进方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105721095A CN105721095A (zh) | 2016-06-29 |
| CN105721095B true CN105721095B (zh) | 2018-03-02 |
Family
ID=56156213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610109863.1A Active CN105721095B (zh) | 2016-02-26 | 2016-02-26 | 一种变电站设备时钟同步的改进方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105721095B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106330378B (zh) * | 2016-09-05 | 2019-03-29 | 华北电力大学 | 一种智能变电站频域时延优化的时间同步方法 |
| CN109343333A (zh) * | 2018-12-11 | 2019-02-15 | 齐鲁工业大学 | 一种利用电压工频信号的配电终端对时方法和系统 |
| CN110505683B (zh) * | 2019-09-23 | 2021-06-08 | 重庆邮电大学 | 一种面向平均一致性时钟同步的频率偏移估计方法 |
| CN111698049B (zh) * | 2020-05-19 | 2023-05-02 | 东南大学 | 一种考虑往返时延不对称的高精度时钟同步方法 |
| CN111830463B (zh) * | 2020-06-02 | 2023-02-21 | 广州南方卫星导航仪器有限公司 | 一种基于tdoa的定位系统、方法及存储介质 |
| CN113687686A (zh) * | 2021-08-10 | 2021-11-23 | 北京小米移动软件有限公司 | 时钟同步方法、装置、电子设备和存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103812595A (zh) * | 2014-03-11 | 2014-05-21 | 重庆邮电大学 | 一种基于ieee1588同步机制的tps时间同步改进算法 |
| CN103888237A (zh) * | 2014-04-04 | 2014-06-25 | 瑞斯康达科技发展股份有限公司 | 一种实现时钟时间同步的方法及装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2429105B1 (en) * | 2010-09-13 | 2013-03-27 | Ntt Docomo, Inc. | Node in a wireless system with time and clock frequency synchronizing and corresponding method |
| CN103368721A (zh) * | 2013-07-23 | 2013-10-23 | 电子科技大学 | 时间触发以太网中透明时钟的计算方法 |
-
2016
- 2016-02-26 CN CN201610109863.1A patent/CN105721095B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103812595A (zh) * | 2014-03-11 | 2014-05-21 | 重庆邮电大学 | 一种基于ieee1588同步机制的tps时间同步改进算法 |
| CN103888237A (zh) * | 2014-04-04 | 2014-06-25 | 瑞斯康达科技发展股份有限公司 | 一种实现时钟时间同步的方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105721095A (zh) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105721095B (zh) | 一种变电站设备时钟同步的改进方法 | |
| CN102104475B (zh) | 基于ieee1588的同步系统及其同步方法 | |
| CN111585683B (zh) | 一种面向时间敏感网络的高可靠时钟同步系统及方法 | |
| CN104836630B (zh) | Ieee1588时钟同步系统及其实现方法 | |
| CN105743598B (zh) | 一种工业以太网时钟同步方法及系统 | |
| CN104737490B (zh) | 通信装置、通信系统和时刻同步方法 | |
| US9178637B2 (en) | Method and devices for synchronization using linear programming | |
| EP2254267B1 (en) | Multi input timing recovery over packet networks | |
| US7835366B2 (en) | Three-way message exchange clock synchronization | |
| CN100561904C (zh) | 一种实现以太网链状网络节点间同步的装置和方法 | |
| CN103959688B (zh) | 一种多时钟域的时钟同步方法、线卡及以太网设备 | |
| CN103812592A (zh) | 基于链状工业以太网的时间同步协议系统及同步方法 | |
| Garner et al. | Synchronization of audio/video bridging networks using IEEE 802.1 AS | |
| CN103888237A (zh) | 一种实现时钟时间同步的方法及装置 | |
| CN102983927B (zh) | 一种基于ieee 1588协议的主从时钟对时的时间补偿方法 | |
| CN106160914B (zh) | 一种基于干扰观测反馈控制技术的ieee1588时钟同步方法 | |
| CN102546142A (zh) | 透明时钟的频率同步方法以及同步报文的存储转发方法 | |
| CN105376043B (zh) | 一种双板卡系统的时间同步方法 | |
| CN105306159A (zh) | 一种时钟的时间戳补偿方法及装置 | |
| CN106603183A (zh) | 一种时间戳过滤方法及装置 | |
| CN109495203A (zh) | 一种ptp从钟的恢复系统 | |
| CN101977105B (zh) | 一种时间延迟不对称差值自动均衡锁相补偿方法 | |
| CN103117829A (zh) | 一种不对称网络间的时间同步与补偿方法或者装置 | |
| CN106656385A (zh) | 中继系统的空口时间同步方法、设备 | |
| WO2021004005A1 (zh) | 时间戳抖动补偿方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |