基于Gardner位定时同步算法

1.1位同步算法在软件无线电接收机中，要正确的恢复出发送端所携带的信号，接收端必须知道每个码元的起止时刻，以便在每个码元的中间时刻进行周期性的采样判决恢复出二进制信号[43]。信号在传播过程中的延时一般是未知的，而且由于传输过程中噪声、多径效应等影响，造成接收到的信号与本地时钟信号不同步，这就需要位同步算法，恢复出与接收码元同频同相的时钟信号。正确的同步时钟是接收端正确判断的基础，也是影响系统误码率的重要因素；没有准确的位同步算法，就不可能进行可靠的数据传输，位同步性能的好坏直接影响整个通信系统的性能[44]。实现位同步算法的种类很多，按照处理方式的不同可分为模拟方式、半数字方式和全数字方式如图3-10所示。a)b)c)图3-10位同步算法模型Fig.3-10BitSynchronousAlgorithmModel图3-10(a)模型为全模拟位同步实现技术，通过在模拟域计算出输入信号的位同步定时控制信号去控制本地时钟，对信号进行同步采样。图3-10(b)模型为半模拟同步模型，该模型的主要思想是通过将采样后的信号经过一系列的数字化处理，提取出输入信号与本地时钟的偏差值，通过这个偏差来改变本地时钟的相位达到位同步。(a)(b)两种方式都需要适时改变本地时钟的相位，不利于高速数字信号的实现且集成化程度较低。图3-10(c)为全数字方式的位同步是目前比较常用方法，全数字方式的位同步算法十分适用于软件无线电的实现。该方法通过一个固定的本地时钟对输入的模拟信号进行采样，将采样后的信号经过全数字化的处理实现同步；采用此种方法，实现简单，且便于数字化实现，对本地时钟的要求大大降低。本次设计主要分析了基于内插方式的Gardner定时恢复算法。1.1.1Gardner定时恢复算法原理Gardner定时恢复算法是基于内插的位同步方式，全数字方式的位同步算法模型中，固定的本地采样时钟不能保证能在信号的极值点处实现采样，所以需要通过改变重采样时钟或输入信号来实现极值处采样[45-46]。Gardner定时恢复算法就是通过改变输入信号的方式实现，利用内插滤波器恢复出信号的最大值再进行重采样，算法原理如图3-11所示。图3-11Gardner定时恢复算法原理Fig.3-11GardnerTimingRecoveryTheory输入信号为离散信号x(mTs)，采样率为Ts，符号周期为T，重采样时钟为Ti，这里的重采样时钟周期Ti=n*T(n为一小整数)。Gardner定时恢复算法的基本思想就是，输入信号x(mTs)经过一个D/A器件和一个模拟滤波器h(t)，将数字信号恢复为模拟信号y(t)进行重采样，得到同步的输出信号y(kTi)。插值滤波器模型中包含了虚拟的D/A变换和模拟滤波器，但是只要具备下面三个条件，则内插完全可以通过数字方式实现。1输入采样序列x(mTs)2内插滤波器脉冲响应h(t)3输入采样时间Ts和输出采样时间Ti也就是说，图中的D/A以及模拟滤波器都可以通过设计数字内插滤波器的方式实现。这里Ts和Ti为固定的两个变量，Ts/Ti不一定为整数，为表示出它们之间的变换过程，通过换算得到Ti和Ts的关系如公式(3-4)所示skksisiTumTTkTkT)()((3-4)mk为比值的整数部分，可看做一个基本指针，表示了本地重采样时钟Ti对采样率为Ts的输入信号的整数倍重采样时刻，而uk为比值的分数部分，指示了滤波器对输入信号的插值时刻。一种典型的Gardner定时恢复算法结构框图如图3-12所示。图3-12Gardner定时恢复算法模型Fig.3-12GardnerTimingRecoveryModel符号速率为T的模拟输入信号x(t)经过本地固定时钟周期Ts采样后变为离散信号x(mTs)（Ts与T满足奈奎斯特基本采样定律）。经过插值滤波器得出的值送入定时误差检测器得出输入信号与本地时钟的相位误差τ(n)，再通过一个环路滤波器滤除其中的噪声及高频成分，将得到的值e(n)送入数控振荡器计算出整数采样时刻mk和插值滤波器插值点位置uk从而得到定时输出y(kTi)。从图3-12中可以看出一个完整的定时恢复算法主要由定时误差检测器、环路滤波器、数控振荡器和插值滤波器组成。其中环路滤波器与前一章中载波同步算法的环路滤波器设计方法相同。这里主要介绍其他模块的设计方法。1.1.2定时误差检测器定时误差检测器采用一种非数据辅助的误差检测算法(Gardner定时误差检测算法)，内插后的信号每个符号内需要两个重采样点，一个点对应信号的最佳采样点；另一个为最佳采样点中间时刻的内插值。定时误差计算公式为：)]1()([)2/1()(nynynyn(3-5)式中，τ(n)为定时误差检测值；y(n)为信号的采样值；n为第n个符号，输出信号的周期为Ti。由(3-5)式可以看出，Gardner算法只需要每个符号周期内的两个采样值，因此取Ti=T/2即可满足算法要求。定时误差检测算法示意图如图3-13所示。图3-13定时误差检测Fig.3-13TimingErrorDetecter该算法具有明显的物理含义。在没有定时误差时，如果有符号转换，则平均的中间采样点应该为零。反之，中间采样点的值不为零，其大小取决于定时误差的大小，或者说中间采样点的值表示了定时误差的大小，但它不能表示定时误差的方向（超前或滞后）。为了表示定时误差的方向，算法考虑中间采样点两边判决点的差值。如果有符号转换，则该差值的符号就表示了定时误差的方向。这样两者的乘积就完全确定了定时误差的大小和方向。如果没有符号转换，则两边采样点的差为零，此时不能获取定时信息。图3-13(a)中表示了当本地采样时钟与插值滤波器输出值同步时，定时误差检测器的采样值；同步时，两个极值采样点均为最大值，中间采样点的值为0，这时环路滤波器的输出值为0表示本地时钟已经与信号同步。(b)图中，表示本地时钟超前的情况，本地时钟超前，则在中间采样点的值为正，表示本地时钟比信号超前，需要内插滤波器向后插值。(c)图中，表示本地时钟比信号滞后，滞后的结果是中间时刻采样点的值为负，需要内插滤波器向前进行插值处理。1.1.3NCO模块设计对于数控振荡器NCO的设计，由于NCO只是用于计算插值点的有效位置，也就不需要采用在ROM表中预存输出波形的采样值。可以根据输入信号来实时产生输出信号脉冲和差值点。NCO计算原理如下图所示图3-14NCO原理图Fig.3-14NCOSchematicDiagramNCO用于对以Ts为采样时钟的输入信号进行抽样。因而NCO的工作时钟与输入信号的工作时钟一致也为Ts，而生成的重采样周期应该与输入信号的符号率同步为Ti。每次NCO寄存器溢出一次则表示要执行一次重采样操作。每次NCO寄存器过零点的时刻(mk+1)Ts便是内插滤波器进行一次运算的时刻（总是位于内插估计点位置的后一个Ts整点采样时刻）。NCO寄存器深度为1，假设当前样点mkTs时刻NCO寄存器的值为η(mk)，环路滤波器输出的控制字为W(mk)，表示每次递减的步进为W(mk)，用差分公式可表示为]1mod[)]()([)1(kkkmWmm当η(mk)W(mk)时，就表示下一个符号周期即将到来，NCO也将产生一次过零点，寄存器的值模1后的值设为下一个符号周期NCO的初始值。从图3-14经过几何分析不难得出：)1(1)1()(kkkkmumu从而得到分数倍插值位置uk为：)()()()1(1)(kkkkkkmWmmmmu通过精确的除法运算，就可以实时的得到分数间隔值uk，这样，内插滤波器的控制参数也就通过NCO完全提取出来。1.1.4插值滤波器设计Gardner定时恢复算法中的插值滤波器主要作用就是通过输入信号x(mTs)与采样点mk与分数插值点uk来实时生成与本地时钟相位相同的信号。插值滤波器输入信号x(mTs)与输出信号y(kTi)的关系可表示为：])[(])[(])[()(21skIIIiskskkiTuihTimxTumykTy(3-6)式中，I1、I2决定插值滤波器的抽头系数，hI为插值滤波器的冲激响应。mk、uk由数控振荡器(NCO)提供，mk决定内插器的整数倍插值位置，它以重采样时钟触发方式体现。uk控制小数倍插值位置直接送给插值滤波器,控制插值点的位置。输出的定时恢复信号的性能主要与插值滤波器的设计方式有很大的关系，下面就来具体分析插值滤波器的实现方法。插值滤波器的实质是对信号经过低通滤波器后再重采样的过程。考虑理想插值情况，根据Shannon定理，采用理想插值可以由带限的输入信号x(t)的抽样值x(mTs)精确得到x(t)在任意时刻的值，即ksssTmTtcmTxtx]/)([sin)()((3-7)其中sssTtTtTtc/)/sin()/(sin(3-8)它的频域表达式为sssTfTfTfH21,021,)((3-9)因而，内插后的序列x(kTi)可表示为：kssisiTmTkTcmTxkTx]/)[(sin)()((3-10)由于理想的内插滤波器是非因果系统，它需要无穷多个信号样值点，物理上具有不可实现性。因而，将理想插值滤波器的脉冲响应进行截断，并根据最优化准则逼近最佳性能。内插滤波器可以通过不同的截断函数得出无穷多种内插函数，但都必须遵守线性相位的条件，即参与插值的采样点数为偶数。常用的内插滤波器包括：两点线性内插滤波器、立方内插滤波器、分段抛物线内插滤波器。这里主要讨论立方插值滤波器的原理和实现结构立方插值滤波器是多项式的插值滤波器的一种，它是基于4点样值的拉格朗日函数2112)()(NNiiiNNxCty(3-11)这里jijNijNjittttC21,，这里N=4，那么N1=N/2=2，N2=N/2-1=-1从而得出立方插值滤波器的时域表达式为：其他,02,16/)/(11)/(6/)/(0,12/)/()/(2/)/(0,12/)/()/(2/)/(2,16/)/(11)/(6/)/()(23232323ssssssssssssssssssTtTTtTtTtTtTtTtTttTTtTtTtTtTTtTtTtth(3-12)归一化令t=(i+u)Ts,则可得h(t)的系数Ci(u)为：;3121616)2)(1()(;121212)2)(1)(1()(;21212)2)(1()(;61616)1)(1()(23123023132uuuuuuuCuuuuuuuCuuuuuuuCuuuuuuC(3-13)对多项式滤波器的实现结构，可采用Farrow结构实现[47]。该结构不必实时计算抽头系数，只需要根据当前时偏u，经过如公式(3-13)的少量计算，就可以得到内插滤波器的系数实现内插。表3-2给出了立方插值滤波器的系数表。图3-15立方插值滤波器Farrow结构实现框图Fig.3-15CubicInterpolationFilterFarrowRealizeStructure表3-2立方插值滤波器Farrow结构实现系数Tab.4-1CubicInterpolationFilterFarrowRealizeStructureCoefficientil0123-20-1/601/6-1011/2-1/201-1/2-11/210-1/3-1/2-1/6通过表3-2可以看出，Farrow结构的插值滤波器实现，每计算一个内插值只需要传送一个变量，即内插估计点值u，并通过简单的计算直接求出内插点的值，而不需要计算中间滤波器的系数(系数为固定值)。图3-15给出了立方插值滤波器的Farrow结构实现框图。1.1.5Simulink算法仿真及性能分析经过以上分析，在matlab中构建了一个BPSK信号的定时恢复模型，符号速率为2MHz，固定采样时钟为