CN110138491A - 一种基于广义加权最小二乘法的gps时钟补偿方法 - Google Patents
一种基于广义加权最小二乘法的gps时钟补偿方法 Download PDFInfo
- Publication number
- CN110138491A CN110138491A CN201910413217.8A CN201910413217A CN110138491A CN 110138491 A CN110138491 A CN 110138491A CN 201910413217 A CN201910413217 A CN 201910413217A CN 110138491 A CN110138491 A CN 110138491A
- Authority
- CN
- China
- Prior art keywords
- pps
- 1pps
- pulse per
- crystal oscillator
- gps
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L7/00—Arrangements for synchronising receiver with transmitter
- H04L7/0016—Arrangements for synchronising receiver with transmitter correction of synchronization errors
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
一种基于广义加权最小二乘法的GPS时钟补偿方法,该方法是通过时间差量测模块接收GPS接收机发送的GPS秒脉冲1PPS与基于本地高精度晶振产生的晶振秒脉冲PPS,量测1PPS与PPS之间的时间差,并将量测的时间差发送到误差补偿模块中的MCU;MCU接收1PPS与PPS之间的时间差,并建立基于广义加权最小二乘法的时钟补偿模型计算1PPS的随机误差和PPS的累积误差,对1PPS和PPS进行补偿。本发明考虑到晶振秒脉冲存在累积误差、GPS秒脉冲存在随机误差,利用两者互补的特点,建立本方法以输出精准秒脉冲。本方法能提高时间同步装置的输出时间的精度。
Description
技术领域
本发明涉及电力系统自动化技术领域,具体涉及一种基于广义加权最小二乘法的GPS时钟补偿方法。
背景技术
对于任何通信设备,都需要时钟为其提供工作频率,所以时钟性能是影响设备性能的一个重要方面。在电力自动化系统中,精度的同步时钟对电力系统广域测量、远程遥控等起着至关重要的作用。对于WAMS广域测控系统而言,其依赖于时间的高度一致,即使是0.12ms的时间偏差,折合成2.16°功角差,也会导致错误控制决策。
GPS接收机是电力系统卫星授时系统中常用设备,但是,GPS接收机存在秒脉冲误差,传统上常采用最小二乘法对GPS接收机进行时钟补偿,在此基础上,提出了普通一元一次线性最小二乘法和普通一元二次最小二乘法。但,普通最小二乘估计算法虽然在特定情况下仍具有无偏性和一致性,但其结果不再是有效估计,比如:均值、方差或自协方差如果受到趋势变动、季节效应和循环变化等影响而发生改变,使得GPS接收机时钟输出的精确度不高。
发明内容
本发明的目的是,针对现有技术中采用普通最小二乘法对GPS接收机进行时钟补偿后所存在的GPS接收机时钟输出精确度不高的问题,提出一种基于广义加权最小二乘法的GPS时钟补偿方法,进一步提高GPS接收机时钟输出的精确度,可为广域测量、故障测距等应用提供精确时标,并满足电网监控系统在线广域授时要求。
为达上述目的,本发明所采用的技术方案是:一种基于广义加权最小二乘法的GPS时钟补偿方法,该方法步骤如下:
步骤1:时间差量测模块接收GPS接收机发送的GPS秒脉冲1PPS与基于本地高精度晶振产生的晶振秒脉冲PPS,量测1PPS与PPS之间的时间差,并将量测的时间差发送到误差补偿模块的MCU;
上述提及的GPS秒脉冲1PPS由GPS接收机对接收的卫星信号解码和处理后得到。晶振秒脉冲PPS是由本地高精度晶振产生的高频振荡计数信号经过锁相环后分频得到(为现有技术)。1PPS与PPS之间的时间差的量测为现有技术,通过鉴相器识别相位差计算获取。
步骤2:误差补偿模块的MCU接收1PPS与PPS之间的时间差,并根据时间序列下基于世界标准时间的1PPS与PPS之间的时间差的计算公式建立基于广义加权最小二乘法的时钟补偿模型,计算1PPS的随机误差和PPS的累积误差,再根据1PPS的随机误差和PPS的累积误差对1PPS和PPS进行补偿;
其中,上述提及的在时间序列下基于世界标准时间的1PPS与PPS之间的时间差的计算公式为 为GPS接收机输出的第i个秒脉冲,εi为1PPS的随机误差;为本地高精度晶振输出的第i个秒脉冲,a为本地高精度晶振产生的秒脉冲序列中第0个秒脉冲与世界标准时间之间的初始偏差,b为时间间隔误差,a+bxi即为PPS的累积误差;和的计算公式中,xi为时间序列,x∈N,i∈N,x1=1,x2=2……xN=N。
上述提及的计算1PPS的随机误差和PPS的累积误差的过程如下:
在时间序列与时间差序列Yi下采用广义加权最小二乘法建模:
该模型中Wi为加权系数,W∈(0,1),W根据本地高精度晶振应用环境及自身物理特性确定;
引入相关系数ρ,相关系数ρ是指:εi=ρεi-1+μt;ρ采用柯克伦-奥克特迭代法求解或者通过德宾-沃森验证推算出;
将式(1)-式(2)×ρ,整理得下式(3):
(Yi-ρYi-1)=a(Wi-ρWi-1)+b(Wixi-ρxi-1*Wi-1)+μt (3)
利用下面式(4)-式(8)进行变量替换:
a*=a(Wi-ρWi-1); (6)
b*=b; (7)
μt=εi-ρεi-1; (8)
将式(3)化成
和能根据式(4)和式(5)计算出,然后根据最小二乘法能求解出a*和b*,然后根据式(6)和式(7)则能计算出a和b;从而,能计算得到GPS输出的第i个秒脉冲的随机误差εi和晶振输出的第i个秒脉冲的累积误差;
补偿过程为:藉由前述计算得到的GPS输出的第i个秒脉冲的随机误差εi和晶振输出的第i个秒脉冲的累积误差,根据公式和公式 计算GPS输出的第i个秒脉冲和晶振输出的第i个秒脉冲,从而对GPS输出的第i个秒脉冲和晶振输出的第i个秒脉冲进行补偿。
本发明考虑到晶振秒脉冲存在累积误差、GPS秒脉冲存在随机误差,利用两者互补的特点,通过分析电力系统卫星授时系统中现有误差补偿系统,根据广义加权最小二乘法的时钟补偿模型求解晶振秒脉冲的累计误差和GPS秒脉冲的随机误差,最后进行相应的补偿,输出精准秒脉冲。本方法能提高时间同步装置的输出时间的精度。
附图说明
图1是本发明的方法流程图。
图2是本发明方法补偿后的PPS的误差与1PPS的随机误差的比较图。
其中:(a)-实验开始;(b)-30min;(c)-60min。
具体实施方式
结合参见图1,本发明为一种基于广义加权最小二乘法的GPS时钟补偿方法,以提高时钟输出的精度,该方法步骤如下:
步骤1:时间差量测模块接收GPS接收机发送的GPS秒脉冲(1PPS)与基于本地高精度晶振(文中又以晶振表示)产生的晶振秒脉冲(PPS),量测1PPS与PPS之间的时间差,并将量测的时间差发送到误差补偿模块的MCU;
其中,GPS秒脉冲1PPS由GPS接收机对接收的卫星信号解码后得到。晶振秒脉冲PPS是由本地高精度晶振产生的高频振荡计数信号经过锁相环后分频得到(为现有技术)。1PPS与PPS之间的时间差的量测为现有技术,通过鉴相器识别相位差计算获取。
步骤2:误差补偿模块的MCU接收1PPS与PPS之间的时间差,并根据时间序列下基于世界标准时间的1PPS与PPS之间的时间差的计算公式建立基于广义加权最小二乘法的时钟补偿模型,计算和补偿GPS秒脉冲(1PPS)的随机误差和本地高精度晶振秒脉冲(PPS)的累积误差;具体过程如下:
在时间序列X∈N(x=1,2,3,4...N)下,基于世界标准时间,设定GPS接收机输出的第i个秒脉冲其中,εi为GPS秒脉冲的随机误差(即GPS秒脉冲与世界标准时间之间的误差);并设定本地高精度晶振产生的秒脉冲序列中第0个秒脉冲与世界标准时间之间的初始偏差为a,时间间隔误差为b,由于本地高精度晶振秒脉冲的随机误差远小于GPS秒脉冲的随机误差(如精度为10-9s的本地高精度晶振秒脉冲的随机误差<1ns),因此不考虑本地高精度晶振秒脉冲的随机误差,则本地高精度晶振输出的第i个秒脉冲a+bxi即为晶振秒脉冲的累积误差;前二式中,i∈N(i=1,2,3,4…N),i=x即x1=1,x2=2……xN=N。因此,晶振秒脉冲与GPS秒脉冲的偏差此处的偏差Yi即为步骤1中量测到的1PPS与PPS之间的时间差。
分析时间序列X与偏差序列Yi的相关特性,考虑到普通最小二乘估计算法虽然在特定情况下仍具有无偏性和一致性,但其结果不再是有效估计,比如:均值、方差或自协方差如果受到趋势变动、季节效应和循环变化等影响而发生改变,故对其进行广义加权最小二乘法估计。其思路是:根据高斯-马尔科夫定理,把不满足经典假设的误差项变成能满足的情况,然后采用广义加权最小二乘法建立模型:
该模型中Wi为加权系数,W∈(0,1),W根据本地高精度晶振应用环境及自身物理特性确定。
引入相关系数ρ,其中,相关系数ρ是指:εi=ρεi-1+μt。对于ρ,可采用柯克伦-奥克特迭代法求解,或者通过德宾-沃森验证推算出,为现有技术。
将(1)-(2)×ρ,整理得式(3):
(Yi-ρYi-1)=a(Wi-ρWi-1)+b(Wixi-ρxi-1*Wi-1)+μt (3)
利用下面式(4)-式(8)进行变量替换:
a*=a(Wi-ρWi-1); (6)
b*=b; (7)
μt=εi-ρεi-1; (8)
将式(3)化成
则:和是Yi和Xi的转换变量,此时把不满足经典假设的误差项变成能满足的情况。从而,在时间序列下,和能根据式(4)和式(5)求出,然后根据最小二乘法能求解出a*和b*(此为现有技术),根据式(6)和式(7)则可计算出a和b。
再根据前述的晶振秒脉冲与GPS秒脉冲的偏差Yi=a+bxi+εi,得到GPS秒脉冲的随机误差εi=Yi-a-bxi,而晶振秒脉冲的累积误差等于a+bxi,从而,藉由a和b的值,能计算得到随机误差εi和累积误差的值。
补偿时,藉由计算得到的GPS输出的第i个秒脉冲的随机误差εi值和晶振输出的第i个秒脉冲的累积误差值,根据公式和公式对GPS秒脉冲和晶振秒脉冲进行补偿,而得到补偿后的GPS和晶振的第i个秒脉冲。则,时间同步装置输出的秒脉冲为补偿后的秒脉冲,即晶振秒脉冲(a+bxi+xi)。
实施例1
利用某公司的时间同步装置采集1PPS(为解码卫星信号所产生,存在随机误差)和本发明方法补偿后输出的PPS,采用示波器进行观测,进行1h的仿真实验。实验开始时,将示波器的参考轴线调整到与补偿后的PPS的上升沿重合,如图2(a)所示;实验30min时1PPS与PPS的对比波形如图2(b)所示;实验60min时的对比波形如图2(c)所示。从图2可以看出,相对于参考轴线,根据本方法产生的PPS的误差小于1PPS的随机误差且在GPS存在随机误差的情况下长时间内仍能保持较高的时间精度,表明采用本方法能提高时间同步装置的输出时间的精度。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所的均等变化与修饰,皆应属本发明的涵盖范围内。
Claims (3)
1.一种基于广义加权最小二乘法的GPS时钟补偿方法,其特征在于,该方法步骤如下:
步骤1:时间差量测模块接收GPS接收机发送的GPS秒脉冲1PPS与基于本地高精度晶振产生的晶振秒脉冲PPS,量测1PPS与PPS之间的时间差,并将量测的时间差发送到误差补偿模块的MCU;
步骤2:误差补偿模块的MCU接收1PPS与PPS之间的时间差,并根据时间序列下基于世界标准时间的1PPS与PPS之间的时间差的计算公式建立基于广义加权最小二乘法的时钟补偿模型,计算1PPS的随机误差和PPS的累积误差,再根据1PPS的随机误差和PPS的累积误差对1PPS和PPS进行补偿;
其中,上述提及的在时间序列下基于世界标准时间的1PPS与PPS之间的时间差的计算公式为 为GPS接收机输出的第i个秒脉冲,εi为1PPS的随机误差,xi为时间序列,x∈N,i∈N,x1=1,x2=2……xN=N;为本地高精度晶振输出的第i个秒脉冲, a为本地高精度晶振产生的秒脉冲序列中第0个秒脉冲与世界标准时间之间的初始偏差,b为时间间隔误差,a+bxi即为PPS的累积误差;
上述提及的计算1PPS的随机误差和PPS的累积误差的过程如下:
在时间序列与时间差序列Yi下采用广义加权最小二乘法建模:
该模型中Wi为加权系数,W∈(0,1),W根据本地高精度晶振应用环境及自身物理特性确定;
引入相关系数ρ,相关系数ρ是指:εi=ρεi-1+μt;ρ采用柯克伦-奥克特迭代法求解或者通过德宾-沃森验证推算出;
将式(1)-式(2)×ρ,整理得下式(3):
(Yi-ρYi-1)=a(Wi-ρWi-1)+b(Wixi-ρxi-1*Wi-1)+μt (3)
利用式(4)-式(8)进行变量替换:
a*=a(Wi-ρWi-1); (6)
b*=b; (7)
μt=εi-ρεi-1; (8)
将式(3)化成
和能根据式(4)和式(5)计算出,然后根据最小二乘法能求解出a*和b*,然后根据式(6)和式(7)则能计算出a和b;从而,能计算得到GPS输出的第i个秒脉冲的随机误差εi和晶振输出的第i个秒脉冲的累积误差;
补偿过程为:藉由前述计算得到的GPS输出的第i个秒脉冲的随机误差εi和晶振输出的第i个秒脉冲的累积误差,根据公式和公式 计算GPS输出的第i个秒脉冲和晶振输出的第i个秒脉冲,从而对GPS输出的第i个秒脉冲和晶振输出的第i个秒脉冲进行补偿。
2.如权利要求1所述的一种基于广义加权最小二乘法的GPS时钟补偿方法,其特征在于,所述GPS秒脉冲1PPS由GPS接收机对接收的卫星信号解码和处理后得到。
3.如权利要求1所述的一种基于广义加权最小二乘法的GPS时钟补偿方法,其特征在于,所述晶振秒脉冲PPS是由本地高精度晶振产生的高频振荡计数信号经过锁相环分频后得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910413217.8A CN110138491B (zh) | 2019-05-17 | 2019-05-17 | 一种基于广义加权最小二乘法的gps时钟补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910413217.8A CN110138491B (zh) | 2019-05-17 | 2019-05-17 | 一种基于广义加权最小二乘法的gps时钟补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110138491A true CN110138491A (zh) | 2019-08-16 |
CN110138491B CN110138491B (zh) | 2021-04-16 |
Family
ID=67575036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910413217.8A Active CN110138491B (zh) | 2019-05-17 | 2019-05-17 | 一种基于广义加权最小二乘法的gps时钟补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110138491B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021077589A1 (zh) * | 2019-10-21 | 2021-04-29 | 泰斗微电子科技有限公司 | 用于测量秒脉冲信号时差的电路和装置 |
CN113419286A (zh) * | 2021-06-21 | 2021-09-21 | 中国科学院地质与地球物理研究所 | 一种卫星与恒温晶振互补的高精度频率合成补偿方法 |
CN114690614A (zh) * | 2022-04-13 | 2022-07-01 | 四川阿秒科技有限公司 | 一种基于时延控制的自主守时系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079687A (zh) * | 2006-05-25 | 2007-11-28 | 上海欣泰通信技术有限公司 | 一种基于最小二乘法模型的时钟调整算法 |
CN101931400A (zh) * | 2010-08-31 | 2010-12-29 | 新邮通信设备有限公司 | 一种基站时钟锁定方法和装置 |
CN103777518A (zh) * | 2014-01-22 | 2014-05-07 | 天津七一二通信广播有限公司 | 基于改进的最小二乘法的卫星导航授时系统及方法 |
CN106527113A (zh) * | 2016-10-31 | 2017-03-22 | 国网福建省电力有限公司 | 一种基于gps的高精度时钟在线校频与授时的方法 |
-
2019
- 2019-05-17 CN CN201910413217.8A patent/CN110138491B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079687A (zh) * | 2006-05-25 | 2007-11-28 | 上海欣泰通信技术有限公司 | 一种基于最小二乘法模型的时钟调整算法 |
CN101931400A (zh) * | 2010-08-31 | 2010-12-29 | 新邮通信设备有限公司 | 一种基站时钟锁定方法和装置 |
CN103777518A (zh) * | 2014-01-22 | 2014-05-07 | 天津七一二通信广播有限公司 | 基于改进的最小二乘法的卫星导航授时系统及方法 |
CN106527113A (zh) * | 2016-10-31 | 2017-03-22 | 国网福建省电力有限公司 | 一种基于gps的高精度时钟在线校频与授时的方法 |
Non-Patent Citations (1)
Title |
---|
菲天吾学: "第五章 自相关性", 《百度文库》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021077589A1 (zh) * | 2019-10-21 | 2021-04-29 | 泰斗微电子科技有限公司 | 用于测量秒脉冲信号时差的电路和装置 |
CN113419286A (zh) * | 2021-06-21 | 2021-09-21 | 中国科学院地质与地球物理研究所 | 一种卫星与恒温晶振互补的高精度频率合成补偿方法 |
CN113419286B (zh) * | 2021-06-21 | 2022-01-14 | 中国科学院地质与地球物理研究所 | 一种卫星与恒温晶振互补的高精度频率合成补偿方法 |
CN114690614A (zh) * | 2022-04-13 | 2022-07-01 | 四川阿秒科技有限公司 | 一种基于时延控制的自主守时系统及方法 |
CN114690614B (zh) * | 2022-04-13 | 2024-03-29 | 四川阿秒科技有限公司 | 一种基于时延控制的自主守时系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110138491B (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Panda et al. | Online estimation of system inertia in a power network utilizing synchrophasor measurements | |
CN110138491A (zh) | 一种基于广义加权最小二乘法的gps时钟补偿方法 | |
CN109299496B (zh) | 一种高精度同步时钟产生方法 | |
CN107345984B (zh) | 一种基于信号识别的自适应同步相量测量方法 | |
CN103616814A (zh) | 一种基于fpga的同步采样时钟闭环校正方法和系统 | |
CN102928803B (zh) | 一种基于同步脉冲输出电源的电子式互感器校验装置 | |
WO2012009365A1 (en) | Systems and methods for electrical power grid monitoring using loosely synchronized phasors | |
CN108896944B (zh) | 一种同步测量装置实验室校准仪及其同步相量测量方法 | |
CN109324228A (zh) | 一种基于多频率-泰勒动态模型的同步相量测量方法 | |
CN105403834A (zh) | 一种发电机动态状态估计方法 | |
CN104300969A (zh) | 一种基于全数字锁相环的高精度同步时钟实现方法 | |
CN110231554A (zh) | 一种gil/gis击穿故障四点声定位方法 | |
CN103970008A (zh) | 一种基于晶振误差补偿的守时方法 | |
CN109284933A (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
CN109870649A (zh) | 一种交流充电桩检定装置的校准系统及方法 | |
CN109560550A (zh) | 基于优化量测的电网谐波状态估计方法 | |
CN113569201A (zh) | 地磁Ap指数预报方法、装置及电子设备 | |
Zhao et al. | Recent development of frequency estimation methods for future smart grid | |
CN202939299U (zh) | 一种基于同步脉冲输出电源的电子式互感器校验装置 | |
Liu | UWB ranging error analysis based on TOA mode | |
CN112511056A (zh) | 一种基于相量测量的鲁棒发电机动态状态估计方法 | |
Albu et al. | Applicability of synchronised measurements in modern distribution grids. A discussion | |
CN103630744B (zh) | Pmu相角检测方法及系统 | |
Pirret | Testing and calibration of phasor measurement units | |
Macii et al. | Impact of wideband noise on synchrophasor, frequency and ROCOF estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |