基于数字滤波的低频时码秒脉冲抖动平滑方法研究.pdf

第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 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自适应滤波算法和卡尔曼滤波算法对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 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 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 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 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 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 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 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 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 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 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年建成的商丘低频时码授时台,采用幅度调制体制,调制速率为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。为提高定时精度,需对秒脉冲的抖动进行平滑。本文将利用接收到的我国商丘低频时码发播台发播的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之间存在一个随机的相位差,该相位差是由于设备上电瞬间本地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的相位差波动情况进行量化反映。对该抖动的相位差进行滤波平滑处理后,结合已标定的时延对本地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 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 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所示自适应滤波器原理,输入信号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 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。卡尔曼滤波可对平稳和非平稳随机信号进行较高精度的估计,由于系统的当前状态只与前一时刻的系统状态以及当前测量值有关,所以无需过去所有的测量值,只需前一时刻的系统状态与当前时刻的测量值 就可以对当 前状态 实 现 最 优 估计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为状态观测矩阵;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为卡尔曼增益;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仿真分析。该接收设备的本地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阶时,相邻采样点之间的波动在1m s左右,甚至在部分时间段波动超出1m s。当阶数增加至1 6阶时,平滑效果稍有改善,随着阶数继续增加,滤波结果的稳定性也明显改善。分别统计8阶、1 6阶、3 2阶和6 4阶LM S自适应滤波器滤波结果标准差以及原始数据标准差如表1所示,原始数据标准差结果为2.8 9 4m s,经LM S自适应滤波处理后,标准差可以提升到0.5 6 9m s。如图1 2所示的LM S自适应滤波器的阶数同滤波结果标准差之间的变化趋势可以发现,标准差的变化趋势随着阶数变大而慢慢趋于平稳,最终稳定在0.5m s附近。统计滤波前后均值可以发现,与原始数据均值相差最大的是采用8阶滤波,相差0.2 0 1m s,与原始数据均值相差最小的是采用3 2阶滤波,相差仅0.0 0 2m s。所以,在实际应用中采用LM S自适应滤波应当选择合适的滤波阶数,以达到较好的滤波效果。虽然随着滤波阶图4 8阶LM S自适应滤波仿真结果F i g.4 S i m u l a t i o nr e s u l t so f8-o r d e rLM Sa d a p t i v e f i l t e r i n g图5 1 6阶LM S自适应滤波仿真结果F i g.5 S i m u l a t i o nr e s u l t so f1 6-o r d e rLM Sa d a p t i v e f i l t e r i n g图6 3 2阶LM S自适应滤波仿真结果F i g.6 S i m u l a t i o nr e s u l t so f3 2-o r d e rLM Sa d a p t i v e f i l t e r i n g26 第3期基于数字滤波的低频时码秒脉冲抖动平滑方法研究图7 6 4阶LM S自适应滤波仿真结果F i g.7 S i m u l a t i o nr e s u l t so f 6 4-o r d e rLM Sa d a p t i v e f i l t e r i n g图8 8阶LM S自适应滤波仿真结果一阶差分F i g.8 F i r s t-o r d e rd i f f e r e n c eo f 8-o r d e rLM Sa d a p t i v ef i l t e r i n g s i m u l a t i o nr e s u l t s图9 1 6阶LM S自适应滤波仿真结果一阶差分F i g.9 F i r s t-o r d e rd i f f e r e n c eo f 1 6-o r d e rLM Sa d a p t i v ef i l t e r i n g s i m u l a t i o nr e s u l t s图1 0 3 2阶LM S自适应滤波仿真结果一阶差分F i g.1 0 F i r s t-o r d e rd i f f e r e n c eo f 3 2-o r d e rLM Sa d a p t i v ef i l t e r i n g s i m u l a t i o nr e s u l t s图1 1 6 4阶LM S自适应滤波仿真结果一阶差分F i g.1 1 F i r s t-o r d e rd i f f e r e n c eo f 3 2-o r d e rLM Sa d a p t i v ef i l t e r i n gs i m u l a t i o nr e s u l t s表1 LM S自适应滤波结果的标准差及均值统计T a b.1 S t a n d a r dd e v i a t i o na n dm e a ns t a t i s t i c so fLM Sa d a p t i v e f i l t e r i n gr e s u l t sm s原始观测数据8阶LM S滤波结果1 6阶LM S滤波结果3 2阶LM S滤波结果6 4阶LM S滤波结果均值6 5.2 4 76 5.0 4 66 5.1 9 46 5.2 4 56 5.2 7 1标准差2.8 9 41.2 4 00.9 2 50.7 1 20.5 6 9图1 2 LM S自适应滤波器阶数与滤波结果标准差关系图F i g.1 2 T h e r e l a t i o n s h i pb e t w e e nt h eo r d e ro fLM Sa d a p t i v ef i l t e ra n dt h e s t a n d a r dd e v i a t i o no f f i l t e r i n gr e s u l t s数的增加,会获得更好的滤波效果,但是在电路实现中,多阶数滤波对应消耗更多的电路逻辑资源。综上所述,LM S自适应滤波可以有效地平滑低频时码秒脉冲的抖动,但是在实际应用中应根据滤波性能和资源使用情况合理选择滤波阶数。2.2 卡尔曼滤波仿真结果分析卡尔曼滤波仿真结果如图1 3所示,可以发现卡尔曼滤波在平滑秒脉冲抖动应用上具有非常明显的效果。图中红线为滤波结果,收敛后的结果接近于一条直线。卡尔曼滤波结果的一阶差分结果如图1 4所示,在滤波的初始阶段大约有近7 0 0s的收敛时间,但收敛后的结果表现十分稳定。计算已收敛时间段内的均值和标准差如表2所示,标准差由36 导航定位与授时2 0 2 3年5月2.8 9 4m s提升至0.0 8 4m s,均值稍有0.3 4 2m s的偏差。图1 3 卡尔曼滤波仿真结果F i g.1 3 K a l m a nf i l t e r s i m u l a t i o nr e s u l t s图1 4 卡尔曼滤波仿真结果一阶差分F i g.1 4 F i r s t-o r d e rd i f f e r e n c e r e s u l t s o fK a l m a n f i l t e r s i m u l a t i o n表2 卡尔曼滤波结果的标准差及均值统计T a b.2 S t a n d a r dd e v i a t i o na n dm e a ns t a t i s t i c so fK a l m a nf i l t e r r e s u l t sm s项目原始观测数据卡尔曼滤波结果均值6 5.2 4 76 4.9 0 5标准差2.8 9 40.0 8 42.3 不同滤波方法对比分析对比L M S自适应滤波和卡尔曼滤波在秒脉冲平滑处理的结果,发现卡尔曼滤波的效果要优于L M S自适应滤波,但是卡尔曼滤波存在一定时间的收敛过程,所以滤波响应时间上,L M S自适应滤波要优于卡尔曼滤波。如图1 5所示,利用A l l a n方差对比两种滤波结果的频率稳定度可以发现,卡尔曼滤波的结果明显优于L M S自适应滤波。所以在实际应用中,可根据需求合理选择滤波策略进行秒脉冲抖动处理,若对定时稳定性及精度要求较高可以选择卡尔曼滤波,若对系统响应时间要求较高,而对定时精度没有过高需求可以选择L M S自适应滤波。图1 5 原始数据与不同滤波方法的A l l a n方差F i g.1 5 A l l a nv a r i a n c eo fo r i g i n a l d a t aa n dd i f f e r e n tf i l t e r i n gm e t h o d s3 结论本文通过研究B P C1 P P S抖动平滑原理,并利用B P C接收机实际采集的B P C1 P P S和本地1 P P S相位差数据,对LM S自适应滤波和卡尔曼滤波秒脉冲抖动平滑仿真分析,得到以下结论:1)B P C信号在传播过程中受环境等因素的影响而发生的秒脉冲抖动可以通过B P C1 P P S和本地1 P P S的相位差反映,对该相位差抖动平滑后用于校准本地1 P P S可以有效提高B P C定时精度和稳定性。2)卡尔曼滤波和LM S自适应滤波对于低频时码秒脉冲的抖动平滑都具有明显效果,卡尔曼滤波的平滑效果和频率稳定度表现更优,收敛后滤波结果标准差为0.0 8 4s。LM S自适应滤波随着滤波器阶数的增加,滤波效果也不断改善,最终滤波结果标准差稳定在0.5m s附近。3)上述两种滤波方法有各自的优势和缺陷,LM S自适应滤波的滤波结果响应时间快,但滤波性能受阶数的影响,而高阶滤波器在电路实现中会消耗更多资源;卡尔曼滤波具有较好的滤波性能,但需要一定的时间进行收敛。所以,在秒脉冲抖动平滑的实际应用中,应根据实际需求合理选择滤波器。参考文献1 刘军,吴贵臣,许林生.低频时码接收机数字化改进研究J.无线电工程,2 0 0 2,3 2(5):5.L i uJ u n,W uG u i c h e n,X uL i n s h e n g.R e s e a r c ho nd i g i t a li m p r o v e m e n to f l o wf r e q u e n c yt i m ec o d er e c e i-v e rJ.R a d i oE n g i n e e r i n g,2 0 0 2,3 2(5):5(i nC h i n e s e).46 第3期基于数字滤波的低频时码秒脉冲抖动平滑方法研究2 刘小花.基于F P G A的便携式B P C定时接收机设计与实现D.西安:中国科学院研究生院(国家授时中心),2 0 1 2.L i uX i a o h u a.D e s i g n a n d i m p l e m e n t a t i o no f p o r t a b l eB P Ct i m i n gr e c e i v e rb a s e do nF P G AD.X i-a n:G r a d u a t eS c h o o l o fC h i n e s eA c a d e m yo fS c i e n c e s(N a t i o n a lT i m eS e r v i c eC e n t e r),2 0 1 2(i nC h i n e s e).3 冯平.低频时码系统附加扩频授时研究D.西安:中国科学院研究生院(国家授时中心),2 0 0 8.F e n gP i n g.R e s e a r c ho na d d i t i o n a ls p r e a ds p e c t r u mt i m i n go f l o wf r e q u e n c y t i m e c o d e s y s t e mD.X i a n:G r a d u a t e S c h o o lo f C h i n e s e A c a d e m y o f S c i e n c e s(N a t i o n a lT i m eS e r v i c eC e n t e r),2 0 0 8(i nC h i n e s e).4 顾卓璟.低频时码接收系统的研究设计D.南京:南京大学,2 0 1 2.G uZ h u o j i n g.R e s e a r c ha n dd e s i g no fl o wf r e q u e n c yt i m ec o d er e c e i v i n gs y s t e mD.N a n j i n g:N a n j i n gU-n i v e r s i t y,2 0 1 2(i nC h i n e s e).5 李国栋.基于F P G A的低频时码数字接收机设计与实现D.西安:中国科学院研究生院(国家授时中心),2 0 0 7.L iG u o d o n g.D e s i g na n di m p l e m e n t a t i o no fl o w-f r e q u-e n c yt i m e-c o d ed i g i t a lr e c e i v e rb a s e do n F P G AD.X i a n:G r a d u a t eS c h o o lo fC h i n e s eA c a d e m yo fS c i e n c e s(N a t i o n a lT i m eS e r v i c eC e n t e r),2 0 0 7(i nC h i n e s e).6 刘小花,许林生,华宇.基于F P G A的便携式B P C定时接收机设计J.时间频率学报,2 0 1 2,3 5(2):8 8-9 5.D O I:1 0.1 3 8 7 5/j.i s s n.1 6 7 4-0 6 3 7.2 0 1 2.0 2.0 0 5.L i uX i a o h u a,X uL i n s h e n g,H u aY u.D e s i g no f p o r t a b l eB P Ct i m i n gr e c e i v e rb a s e do nF P G AJ.J o u r n a lo fT i m ea n dF r e q u e n c y,2 0 1 2,3 5(2):8 8-9 5.D O I:1 0.1 3 8 7 5/j.i s s n.1 6 7 4-0 6 3 7.2 0 1 2.0 2.0 0 5(i nC h i n e s e).7 刘 雄 飞,高 金 定,齐 海 兵.LM S自 适 应 滤 波 器F P G A实现 的 新 方 法 J.压 电 与 声 光,2 0 0 7,2 9(1):3.L i uX i o n g f e i,G a oJ i n d i n g,Q iH a i b i n g.N e w m e t h o df o rF P G Ai m p l e m e n t a t i o no fLM Sa d a p t i v e f i l t e rJ.P i e z o e l e c t r i ca n dA c o u s t o-o p t i c,2 0 0 7,2 9(1):3(i nC h i n e s e).8 夏晓.自适应滤波器中LM S算法的研究及应用D.北京:北京邮电大学,2 0 1 3.X i aX i a o.R e s e a r c ha n da p p l i c a t i o no fLM Sa l g o r i t h mi na d a p t i v ef i l t e rD.B e i j i n g:B e i j i n g U n i v e r s i t yo fP o s t sa n dT e l e c o mm u n i c a t i o n s,2 0 1 3(i nC h i n e s e).9 马伟富.自适应LM S算法的研究及应用D.成都:四川大学,2 0 0 5.M aW e i f u.R e s e a r c ha n da p p l i c a t i o no f a d a p t i v eLM Sa l g o r i t h mD.C h e n g d u:S i c h u a n U n i v e r s i t y,2 0 0 5(i nC h i n e s e).1 0 彭丁聪.卡尔曼滤波的基本原理及应用J.软件导刊,2 0 0 9,8(1 1):3 2-3 4.P e n gD i n g c o n g.B a s i cp r i n c i p l e sa n da p p l i c a t i o n so fK a l m a nf i l t e r i n gJ.S o f t w a r e,2 0 0 9,8(1 1):3 2-3 4(i nC h i n e s e).1 1 严恭敏,邓瑀.传统组合导航中的实用K a l m a n滤波技术评 述 J.导 航 定 位 与 授 时,2 0 2 0,7(2):5 0-6 4.D 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 0.0 2.0 0 7.Y a nG o n g m i n,D e n gY u.R e v i e wo f p r a c t i c a lK a l m a n f i l-t e r i n gt e c h n i q u e s i nt r a d i t i o n a l i n t e g r a t e dn a v i g a t i o nJ.N a v i g a t i o n,P o s i t i o n i n ga n dT i m i n g,2 0 2 0,7(2):5 0-6 4.D 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 0.0 2.0 0 7(i nC h i n e s e).1 2 杨丹.卡尔曼滤波器设计及其应用研究D.湘潭:湘潭大学,2 0 1 4.Y a n gD a n,.D e s i g na n da p p l i c a t i o no fK a l m a n f i l t e rD.X i a n g t a n:X i a n g t a nU n i v e r s i t y,2 0 1 4(i nC h i n e s e).1 3 王勋,左启耀,洪诗聘,等.一种基于卡尔曼滤波的定位解算性能评估新方法J.导航定位与授时,2 0 1 8,5(5):7 3-8 2.D O I:1 0.1 9 3 0 6/j.c n k i.2 0 9 5-8 1 1 0.2 0 1 8.0 5.0 1 3.W a n gX u n,Z u oQ i y a o,H