发明名称 基于卫星导航和Kalman滤波的高精度授时方法
摘要 本发明设计了一种基于卫星导航和Kalman滤波的高精度授时系统及方法,该系统包含数据解算与授时滤波两部分。数据解算部分采用传统的最小二乘法解算出北斗二代接收机的位置信息以及钟差信息,之后利用Kalman滤波器对钟差信息进行滤波得到秒信号发生器的相位误差信息以及频率误差信息,最后利用滤波得到秒信号发生器的相位误差信息以及频率误差信息对秒信号发生器的初始相位以及频率计数字进行调整,实现高精度授时的功能。本发明解决了现有技术条件下授时精度不高,以及秒信号发生器易受晶振频率误差影响的问题。
申请公布号 CN103176190B 申请公布日期 2015.01.28
申请号 CN201310069980.6 申请日期 2013.03.06
申请人 西北工业大学 发明人 何伟;廉保旺;唐成凯;雷浩;廉伟林;杨琼
分类号 G01S19/27(2010.01)I 主分类号 G01S19/27(2010.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种基于卫星导航和Kalman滤波的高精度授时方法,采用以下步骤:步骤1:卫星接收模块接受卫星位置信息与伪距信息,根据卫星位置信息与伪距信息,采用最小二乘法解算卫星接收模块的钟差信息δt<sub>g</sub>;步骤2:对卫星接收模块的钟差信息进行Kalman滤波,对于第k次滤波具体步骤如下:步骤2.1:Kalman滤波的时间更新过程:对状态向量X<sub>k‑1</sub>=[δt<sub>k‑1</sub> δf<sub>k‑1</sub>]进行估计,得到状态估计值<img file="FDA0000541362580000011.GIF" wi="93" he="97" />其中δt<sub>k‑1</sub>为秒信号发生器的钟差信息、δf<sub>k‑1</sub>为秒信号发生器频差信息,估计式为<img file="FDA0000541362580000012.GIF" wi="291" he="97" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000541362580000013.GIF" wi="258" he="166" /></maths>当k=1时,X<sub>k‑1</sub>取值为[0 0];根据式<img file="FDA0000541362580000014.GIF" wi="350" he="84" />得到协方差阵P<sub>k‑1</sub>的预测矩阵<img file="FDA0000541362580000015.GIF" wi="80" he="97" />其中Q为过程噪声协方差阵;步骤2.2:Kalman滤波的量测更新过程:根据式<img file="FDA0000541362580000016.GIF" wi="497" he="83" />得到增益矩阵K<sub>k</sub>,其中C=[1 0],R为量测噪声协方差阵;根据式<img file="FDA0000541362580000017.GIF" wi="576" he="89" />得到滤波后的状态向量矩阵X<sub>k</sub>,根据式<img file="FDA0000541362580000018.GIF" wi="397" he="83" />得到更新后的协方差矩阵P<sub>k</sub>,其中I为二维单位矩阵;步骤3:利用每一次滤波得到的秒信号发生器的钟差信息乘以晶振的频率f<sub>0</sub>得到秒信号发生器的相位控制字N<sub>p</sub>的调整量ΔN<sub>p</sub>,用N<sub>p</sub>减去调整量ΔN<sub>p</sub>得到精确的秒信号发生器的相位控制字;利用每一次滤波得到的秒信号发生器的频差信息乘以晶振的频率f<sub>0</sub>得到秒信号发生器的频率控制字N<sub>f</sub>的调整量ΔN<sub>f</sub>,用N<sub>f</sub>减去调整量ΔN<sub>f</sub>得到精确的秒信号发生器的频率控制字;根据得到的精确的秒信号发生器的相位控制字和频率控制字,实现高精度授时。
地址 710072 陕西省西安市友谊西路127号