基于数字滤波的低频时码秒脉冲抖动平滑方法研究.pdf
《基于数字滤波的低频时码秒脉冲抖动平滑方法研究.pdf》由会员分享,可在线阅读,更多相关《基于数字滤波的低频时码秒脉冲抖动平滑方法研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、第1 0卷 第3期2 0 2 3年5月导航定位与授时N a v i g a t i o nP o s i t i o n i n ga n dT i m i n gV o l.1 0N o.3M a y2 0 2 3d o i:1 0.1 9 3 0 6/j.c n k i.2 0 9 5-8 1 1 0.2 0 2 3.0 3.0 0 6基于数字滤波的低频时码秒脉冲抖动平滑方法研究齐 震1,2,3,程浪浪1,2,3,刘 强1,2,王 昕1,2,3,赵 凡1,2,3,高玉平1,3,冯 平1,2,3(1.中国科学院国家授时中心,西安7 1 0 6 0 0;2.封开低频时码授时台,广东 肇庆5 2
2、 6 5 0 0;3.中国科学院大学,北京1 0 0 0 4 9)摘 要:低频时码授时信号在接收时,由于环境等因素的影响会导致解调得到的低频时码秒脉冲出现抖动。为减小该抖动现象对低频时码定时精度的影响,基于我国低频时码授时系统(B P C),研究并设计基于数字滤波的秒脉冲抖动平滑方法,通过仿真实验,验证分析了最小均方误差(LM S)自适应滤波算法和卡尔曼滤波算法的有效性和抖动平滑性能。结果表明,低频时码秒脉冲的抖动现象可以通过B P C1 P P S(o n ep u l s ep e rs e c o n d)和本地1 P P S之间的相位差波动情况反映;LM S自适应滤波算法和卡尔曼滤波算
3、法对B P C1 P P S的抖动平滑处理均有明显效果,卡尔曼滤波算法的抖动平滑效果更优,但存在一定的收敛时间,LM S自适应滤波算法的滤波结果响应速度快,但抖动平滑性能受滤波器阶数的影响。所以,在实际应用中,应根据实际需求选择合适的滤波器及相关参数。关键词:抖动平滑;数字滤波;低频时码;秒脉冲中图分类号:V 1 1 文献标志码:A 文章编号:2 0 9 5-8 1 1 0(2 0 2 3)0 3-0 0 5 9-0 7R e s e a r c ho nj i t t e r s m o o t h i n gm e t h o do f l o w-f r e q u e n c yt i
4、 m e-c o d e1 P P Sb a s e do nd i g i t a l f i l t e r i n gQ IZ h e n1,2,3,CHE N GL a n g l a n g1,2,3,L I U Q i a n g1,2,WANGX i n1,2,3,Z HAOF a n1,2,3,G AOY u p i n g1,3,F E N GP i n g1,2,3(1.N a t i o n a lT i m eS e r v i c eC e n t e r,C h i n e s eA c a d e m yo fS c i e n c e s,X i a n7 1
5、0 6 0 0,C h i n a;2.F e n g k a iL o w-F r e q u e n c yT i m e-C o d eT i m eS e r v i c eS t a t i o n,Z h a o q i n g,G u a n g d o n g5 2 6 5 0 0,C h i n a;3.U n i v e r s i t yo fC h i n e s eA c a d e m yo fS c i e n c e s,B e i j i n g1 0 0 0 4 9,C h i n a)A b s t r a c t:Wh e n t h e l o w-f
6、 r e q u e n c y t i m e-c o d e s i g n a l i s r e c e i v e d,d u e t o t h e i n f l u e n c eo f e n v i r o n m e n-t a l f a c t o r s,t h ed e m o d u l a t e dl o w-f r e q u e n c yt i m e-c o d e1 P P S(o n ep u l s ep e rs e c o n d)w i l l j i t t e r.I no r d e r t or e d u c e t h e i
7、 n f l u e n c eo f l o w-f r e q u e n c yt i m e-c o d e1 P P Sj i t t e ro nt i m i n ga c c u r a c y,b a s e do nl o w-f r e q u e n c yt i m e-c o d e t i m i n gs y s t e m(B P C)i nC h i n a,a n1 P P S j i t t e r s m o o t h i n gm e t h o d i s s t u d-i e da n d d e s i g n e d o n t h e
8、b a s i s o f d i g i t a lf i l t e r i n g.T h r o u g h s i m u l a t i o n e x p e r i m e n t s,t h ee f f e c t i v e n e s sa n d j i t t e r s m o o t h i n gp e r f o r m a n c eo f l e a s tm e a ns q u a r e(LM S)a d a p t i v e f i l t e r a n dK a l-m a nf i l t e ra r ev e r i f i e da
9、 n da n a l y z e d.T h e r e s u l t s s h o wt h a t t h e j i t t e ro fB P C1 P P Sc a nb e r e f l e c t e db yt h ep h a s ed i f f e r e n c ef l u c t u a t i o nb e t w e e nB P C1 P P Sa n dl o c a l1 P P S.B o t hLM Sa d a p t i v ef i l t e ra n dK a l m a nf i l t e rh a v eo b v i o u
10、ss m o o t h i n ge f f e c t so nt h ej i t t e ro f1 P P S,a m o n gw h i c h K a l m a n 收稿日期:2 0 2 2-0 9-3 0;修订日期:2 0 2 3-0 1-0 5作者简介:齐震(1 9 9 7-),男,博士研究生,主要从事无线电通信及测量技术方面的研究。通信作者:冯平(1 9 7 5-),男,研究员,主要从事授时技术及信号处理研究、无线电通信及测量技术方面的研究。导航定位与授时2 0 2 3年5月f i l t e rh a sb e t t e r j i t t e rs m o o t
11、 h i n gr e s u l t sb u tw i t hac o n v e r g e n c e t i m e.T h e f i l t e r i n gr e s u l to fLM Sa d a p t i v e f i l t e rh a s a s h o r t r e s p o n s e t i m e,b u t i t s j i t t e r s m o o t h i n gp e r f o r m a n c e i s a f f e c t e db y t h eo r-d e ro f f i l t e r.T h e r e
12、f o r e,i np r a c t i c a l a p p l i c a t i o n s,a p p r o p r i a t e f i l t e r s a n dr e l a t e dp a r a m e t e r s s h o u l db es e l e c t e da c c o r d i n gt oa c t u a ln e e d s.K e yw o r d s:J i t t e r s m o o t h i n g;D i g i t a l f i l t e r i n g;L o w-f r e q u e n c y t
13、i m e-c o d e;O n ep u l s ep e r s e c o n d0 引言低频时码授时技术属于长波授时技术的一种扩展,是国际电信联盟(I n t e r n a t i o n a lT e l e c o mm u n i c a t i o nU n i o n,I T U)一直提倡的一种授时技术1。目前,已有多国建有低频时码授时系统,如中国B P C、美国WWV B及日本J J Y等2。低频时码授时的优势在于借助低频载波通过天波和地波进行时间信号传输,具有较高的相位稳定性,传播损耗小,可以远距离传输1。我国于2 0 0 7年建成的商丘低频时码授时台,采用幅度调制体
14、制,调制速率为1b i t/s,通过4种不同脉宽传递时间信息,在每秒开始时刻使幅度下跌9 0%标志U T C(N T S C)秒的起始时刻3-4。其中,低频时码授时系统的系统时间与U T C(N T S C)的时间同步通过卫星共视实现。随着数字电子技术的不断发展,低频时码接收机也朝着数字化、小型化的趋势发展,所采用的定时手段也趋向数字化解决方案5。低频时码授时信号从发播台至用户接收终端的过程中,会受到各种环境的干扰以及信号传播链路时延的影响,因此接收端往往接收到的低频时码秒脉冲会存在随机抖动现象,在一定程度上影响定时精度6。为提高定时精度,需对秒脉冲的抖动进行平滑。本文将利用接收到的我国商丘低
15、频时码发播台发播的B P C授时 信 号,分 别 对 最 小 均 方 误 差(l e a s t m e a ns q u a r e,LM S)自适应滤波和卡尔曼滤波这两种数字滤波方法在低频时码秒脉冲抖动平滑处理中的平滑性能进行仿真验证分析。1 原理及方法1.1 低频时码秒脉冲抖动平滑方法B P C授时信号到达接收终端后通过解调得到的1 P P S与国家授时中心标准时间UT C(NT S C)1 P P S之间存在一个由于传播路径及设备时延等因素产生的固定相位差,该相位差可通过理论模型计算和接收终端设备标定获得。同时,接收终端的内部工作时钟源产生的本地1 P P S与B P C1 P P S
16、之间存在一个随机的相位差,该相位差是由于设备上电瞬间本地1 P P S初始相位随机产生导致,可通过接收终端内部测量获得,每次断电重启后会发生变化。而在实际B P C信号接收时,受接收环境干扰等影响,解调得到的B P C1 P P S存在抖动现象。因此B P C1 P P S与UT C(NT S C)1 P P S之间的相位差会对应 地 出 现 幅 度 抖 动,同 理,B P C1 P P S与 本 地1 P P S之间的相位差会出现同样幅度变化的抖动。所以,接收终端解调得到的B P C1 P P S的抖动现象可以通过接收终端内部测量B P C1 P P S与本地1 P P S的相位差波动情况进
17、行量化反映。对该抖动的相位差进行滤波平滑处理后,结合已标定的时延对本地1 P P S进行修正,使其与UT C(NT S C)1 P P S保持同步,接收终端即可输出标准稳定的1 P P S。原理如图1所示。图1 B P C接收终端秒脉冲同步修正原理图F i g.1 B P Cr e c e i v e r1 P P Ss y n c h r o n i z a t i o nc o r r e c t i o ns c h e m a t i cd i a g r a m B P C接收终端秒信号修正原理如图2所示。天线接 收B P C授 时 信 号 并 解 调,然 后 得 到B P C1 P
18、 P S(o n ep u l s ep e rs e c o n d),同时本地高稳晶振提供的工作频率信号经分频倍频得到本地1 P P S,再对06 第3期基于数字滤波的低频时码秒脉冲抖动平滑方法研究上述两个1 P P S的脉冲边沿进行捕获,利用时钟计数器得到两个脉冲之间的时差值,该时差值即为B P C1 P P S和本地1 P P S之间的相位差。对获得的时差值进行数字滤波,将滤波结果作为本地1 P P S进行第一次移相修正,然后再附加时延修正,即可输出与UT C(NT S C)1 P P S保持同步的1 P P S。图2 B P C接收终端秒信号修正原理图F i g.2 B P Cr e
19、 c e i v e r s e c o n ds i g n a l c o r r e c t i o ns c h e m a t i cd i a g r a m1.2 数字滤波原理图3 LM S自适应滤波器原理图F i g.3 P r i n c i p l eo fLM Sa d a p t i v e f i l t e r1.2.1 LM S自适应滤波本文采用基于最小均方误差(l e a s tm e a ns q u a r e,L M S)算法的自适应滤波方法,是一种改进的最陡下降算法,不需要对信号提取先验统计,具有计算量小、易于实现的特点7。如图3所示自适应滤波器原理,输
20、入信号x(n)经过可调节参数滤波器后得到输出信号y(n),再与期望信号d(n)对比后得到误差信号e(n),然后通过自适应算法对滤波器参数进行调整,最终使e(n)的均方值最小8。自适应滤波利用前一时刻已得的滤波器参数的结果,自动调节当前时刻的滤波器参数,以适应信号和噪声随时间变化的统计特性,从而实现最优滤波9。所以,LM S自适应滤波算法主要分为两个部分:滤波和和自适应调整。其算法过程如下7y(n)=Ni=0wi(n)xi(n)(1)e(n)=d(n)-y(n)(2)w(n+1)=w(n)+e(n)x(n)(3)式中,N为滤波器的阶数;为步长因子,取值需要满足公式(4),即01/m a x,m
21、a x为输入信号自相关矩阵的最大特征值。1-2m a x1(4)根据LM S自适应滤波原理,在对秒脉冲进行滤波处理时需设置4个输入信号,分别为时差信号、时差期望信号、步长因子信号和滤波器阶数信号。时差信号为B P C1 P P S和本地1 P P S之间的相位差;时差期望信号为一定时间内时差信号的均值,同时利用该时间段内的时差信号计算自相关矩阵的最大特征值,从而确定步长因子;滤波器阶数可根据硬件资源情况确定,一般为8阶,而本文通过仿真进一步分析了不同阶数对滤波性能的影响。1.2.2 卡尔曼滤波卡尔曼滤波利用递归算法解决时域离散数据的滤波问题,是基于最小误差准则的最优线性估计1 0-1 1。卡尔
22、曼滤波可对平稳和非平稳随机信号进行较高精度的估计,由于系统的当前状态只与前一时刻的系统状态以及当前测量值有关,所以无需过去所有的测量值,只需前一时刻的系统状态与当前时刻的测量值 就可以对当 前状态 实 现 最 优 估计1 2-1 3。卡尔曼滤波是通过预测和更新两个阶段实现对系统的最优估计,基本模型如下1 4。首先构建一个离散线性系统模型,见式(5)、式(6)1 5xk=Axk-1+Buk+k-1(5)zk=Hxk+vk(6)其中,xk为k时刻系统的状态矩阵;xk-1为k-1时刻系统的状态矩阵;A为状态转移矩阵;B为控制输入矩阵;uk为系统控制输入;k-1为过程噪声;zk为k时刻系统的观测量;H
23、为状态观测矩阵;vk为系统噪声。基于该离散系统模型,假设过程噪声与系统噪声均为高斯白噪声,即P()N(0,Q),P(v)N(0,R),其中Q和R分别为和v的协方差。可得16 导航定位与授时2 0 2 3年5月预测方程为1 5x-k=Axk-1+Buk-1(7)p-k=Apk-1AT+Q(8)更新方程为K=p-kHT(Hp-kHT+R)-1(9)xk=x-k+K(zk-Hx-k)(1 0)pk=(E-KH)p-k(1 1)其中,x-k为当前状态的预测值;xk-1为前一时刻的最优估计值;xk为当前状态的最优估计值;p-k为系统k时刻先验误差协方差;pk-1为系统k-1时刻误差的协方差;K为卡尔曼增
24、益;pk为系统k时刻误差的协方差;E为单位矩阵;其他同上。在秒脉冲卡尔曼滤波处理中,由于信号在高斯信道中传输,当前时刻系统的预测值为上一时刻系统的输出值,所以令状态转移矩阵A设为1,控制输入矩阵B设为0,状态观测矩阵H设为1,单位矩阵E设为1。故式(7)(1 1)可更新为x-k=xk-1(1 2)p-k=pk-1+Q(1 3)K=p-k(p-k+R)-1(1 4)xk=x-k+K(zk-x-k)(1 5)pk=(1-K)p-k(1 6)2 仿真验证分析本文利用已有的B P C接收设备所采集的B P C1 P P S和本地1 P P S时差作为滤波模块的输入信号进行MA T L A B仿真分析。
25、该接收设备的本地1 P P S由J K O C 3 6 A-1 0MH z型号的恒温晶振提供,其频率准确度为0.0 31 0-6,日老化率为0.51 0-9。数据选用2 3点后采集的1h数据,共36 0 0个采样点。2.1 LM S自适应滤波仿真结果分析LM S自适应滤波的仿真结果如图4图7所示。为研究滤波器不同阶数对滤波结果的影响,本文分别从8阶、1 6阶、3 2阶和6 4阶分析LM S自适应滤波的性能。根据仿真结果发现,LM S自适应滤波算法可以有效平滑秒脉冲的抖动现象,并且随着滤波器阶数的增加,平滑效果越好。对滤波结果一次差分如图8图1 1所示,在滤波器阶数为8阶时,相邻采样点之间的波动
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 数字 滤波 低频 时码秒 脉冲 抖动 平滑 方法 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。