HB滤波器在伽利略信号接收中应用及实现.docVIP

HB滤波器在伽利略信号接收中应用及实现 　　摘要：信号带宽的增加使得数字接收机对中频信号采样的速率越来越高。在工程实现常常要对采样后的信号进行抽样处理实现降速。当抽取因子为2时，可以在抽取后采用HB滤波器对信号进行处理。HB滤波器是一种特殊的FIR数字滤波器，它的的通带和阻带关于二分之一Nyquist频率对称，因而有近一半的滤波器系数为0，滤波器运算量小，实时性高，易于硬件实现，在多速率信号处理有着特别重要的意义。本文所描述的方案中，HB滤波器是通过MATLAB进行设计和仿真，获得参数后，在FPGA中硬件实现并且应用在伽利略信号的接收处理中。 　　关键词：HB滤波器 抽取 多相 MATLAB 数字接收机 　　中图分类号：TN713+.1 文献标识码：A 文章编号：1007-9416(2011)07-0100-02 　　 　　1、引言 　　目前工程中使用的数字接收机常常采用中频接口，根据接收信号的不同，中频信号的带宽有很大的区别。比如伽利略卫星E5采用了选择性BOC（AltBOC（15，10））调制方式，它引入了更大的带宽以及复杂的4个扩频码成份的交互作用，因此该信号的处理技术较传统的BPSK或普通的BOC调制而言具有更大的难度，。针对具有大带宽伽利略E5信号中频接收，为了不发生频谱混叠需要更高的采样速率。高采样率带来的问题就是采样后的数据流速率很高，导致后面的数字信号处理环节速度跟不上，从而不能实现对数据的实时处理。比如DDC过程中用到的乘法器就很难实现很高的速度。这样就需要对A/D后的数据进行抽取处理来实现降速的目的。 　　半带滤波器是一种基于抽取和内插原理的滤波器。它是一种特殊的低通滤波器――一种抽取因子为2的抽取滤波器。半带滤波器的通带和阻带相对于二分之一Nyquist频率对称，在设计时有近一半的滤波器系数为0。在实际的硬件实现中，系数为0的部分不消耗运算量，运算量减少带来的是计算效率的提高，实时性的提高，有利于滤波器的实时工作。同时运算量的减少在硬件上直接反应为所占用资源的减少，这在工程实现中也具有重要的意义，对芯片的选型，成本的控制都有好处。因此半带滤波器在数字接收机的设计中，在多速率信号处理中具有重要的有着重要的地位，并且得到了广泛的应用。 　　本文介绍了半带滤波器的性质和一种设计方法，然后以伽利略E5信号的接收为例，实现了该设计方法。 　　2、HB滤波器的性质和通用设计 　　2.1 HB滤波器的性质 　　图1为抽取速率为2的半带滤波器的框图。其中参数依次为:Hh（z）为抽取滤波器，fs为采样速率，fsl=1/2fs为输出采样速率。 　　为了保证抽取后的频谱不发生混叠，必须滤除fs/4以上的频率分量。如果不能滤除，则信号将会以fs/4为中心折叠刀有用的信号带内。 　　图2所展示的是一种特殊的FIR滤波器的频率特性，这就是HB滤波器，它具有如下的特性: 　　(1)通带波纹δp与阻带衰减δs相等; 　　(2)通带边频Fp与阻带边频Fa相对f s / 4 对称。 　　(3)HB滤波器近一半的系数精确为零，做2倍抽取时可以成倍地减少运算量，计算效率高、实时性强。 　　2.2 HB滤波器的设计 　　HB滤波器在设计的过程中可以采用窗函数法和频率采样法。由于频率采样法运算比较复杂，通常在工程实现中一般多采用窗函数法来完成HB滤波器的设计。用窗函数设计HB滤波器，窗函数可以有多种选择，例如:Bartlett窗，矩形窗，Hanning窗，Hamming窗，Blackman窗，Kaiser窗等。 　　采用窗函数法进行HB滤波器的设计的主要步骤有: 　　(1)计算相对过渡带宽; 　　(2)确定滤波器节数; 　　(3)计算窗函数; 　　(4)获得滤波器系数。 　　在MATLAB中可以采用的函数有remezord.m，remez.m，fir1.m，fir2.m等。可以根据实际需求采用不同的函数进行设计。 　　3、HB滤波器工程应用 　　伽利略E5信号经过DDC混频后，中心频率为1KHz，频谱带宽为52MHz，采样频率为165MHz，即，(＝0.31，＝0.69)。为了防止由于抽取引入的混叠，采用半带滤波器对信号进行滤波处理。 　　如果采用普通的滤波抽取结构，对器件的运算数率要求很高，这表现在抽取器模型中的低通滤波器位于抽取算子之前，也就是说低通滤波器是在降速之前实现的，其抗混叠数字滤波在高取样率条件下进行，这无疑大大提高了对运算速度的要求，对实时处理是极其不利的。 　　设数字滤波器的冲击响应为，则其Z变换定义为: 　　上式即为数字滤波器的多相滤波结构，其网络图如图3所示。将其应用于抽取器，并注意到抽取器的等效关系。即可得到抽取器的多相滤波结构如图4所示。 　　线性相位半带FIR