根升余弦脉冲成形滤波器FPGA实现

根升余弦脉冲成形滤波器FPGA实现

作者：赵林军

来源：《现代电子技术》2011年第01期

摘 要：提出了基于电路分割技术实现通信系统发送端根升余弦波形成形滤波器查表法的FPGA结构，节省了ROM单元，讨论了其ROM初始化时形波数据的组织方法，完成了该结构的VHDL实现，给出了该设计在Modelsim环境下的时序仿真结果。通过对仿真结果分析，表明所述的设计方法是可行的。该设计方案不随波形样本数目的增多而使电路系统变得更为复杂，它所实现的成形滤波器满足于高速成形的应用需求。 关键词：根升余弦； 成形滤波器； 查找表； FPGA 中图分类号：TN391.9-34文献标识码：A 文章编号：1004-373X(2011)01-0023-03

FPGA Implementation of Square Root Raised Cosine Pulse Shaping Filter ZHAO Lin-jun

(Department of Electronics & Information Engineering, Shaanxi University of Technology, Hanzhong 723003, China)

Abstract： The FPGA implementation of look-up table arithmetic-based square root raised cosine pulse shaping filter in communication system based on circuit segmentation is proposed, and how to calculate the initialization ROM data of a forming waveform is discussed. On this basis, the VHDL code and the Modelsim simulation results are elucidated. The design method is simple and easy, and is very suitable for high speed forming application.

Keywords： square root raised cosine; shaping filter; look-up table; FPGA 0 引 言

数字通信系统中，基带信号的频谱一般较宽，因此传递前需对信号进行成形处理，以改善其频谱特性，使得在消除码间干扰与达到最佳检测接收的前提下，提高信道的频带利用率。目前，数字系统中常使用的波形成形滤波器有平方根升余弦滤波器、高斯滤波器等。设计方法有卷积法或查表法［1-3］，其中：卷积法［4］的实现，需要消耗大量的乘法器与加法器，以构

龙源期刊网 http://www.qikan.com.cn

成具有一定延时的流水线结构。为降低硬件消耗，文献［5］提出了一种分布式算法(Distributed Arithmetic，DA)的滤波器设计结构。它将传统的乘、累加运算转化为移位、累加运算，当运算数据的字宽较小时，极大地降低了硬件电路的复杂度，提高了响应速度;当运算数据的字长较长时，因其需要更多的移位迭代运算而不适合高速处理的需求。为此，文献［2］提出了采用滤波器的多相结构与改进DA算法相结合的一种设计方法。当考虑ISI码元数目较多时，该设计所需要的ROM表个数就会增多，同时访问ROM所需的地址的产生电路就会变得更为复杂。文献［3,6］提出了采用存储器分割技术，可以降低ROM单元的数量，但是它是以增加系统的复杂性与响应时延、信号毛刺为代价的。文献［2,7］在滤波器设计时采用了CSD编码，虽然减少了乘法运算，但是需要设计CSD编解码电路。

文中论述的是二进制基带信号的连续查表法平方根升余弦波形成形滤波器(SRRC)的FPGA实现(滚降系数取0.22)，取冲击响应截断时间为8T，每T内样点数为8个，所用ROM单元数为

，每单元数据为

有符号整型数。查找ROM表所需11 b的地址由一个长

的数据移位寄存器与一个模8的采样时钟计数器链接而成。给出了设计在Modelsim 6.3

下的时域仿真波形，经与理论相比较，文中的设计方法是可行的，且当二进制码元的码间干扰数增多(码间样点增加)时，地址电路简单增长即可(不影响响应时间)，便于FPGA的实现。 1 二进制基带信号平方根升余弦成形原理

实际系统中，广义信道传递函数H(f)由发送滤波器器 (1)

根据乃奎斯特第一准则，当H(f)幅频特性满足的滚降系数为α升余弦滤波器特征时，可以实现无ISI传输时刻降低对采样时钟精度的要求，当信道噪声可以忽略时，取照接收滤波器的输出信噪比最大准则，有：

1,0≤f≤1－α2T

121+cosπ(2Tf－1+α)2α, 1－α2T＜f≤1+α2T 0,f＞1+α2T

，按

三部分共同构成，

即：

、信道

、接收滤波

龙源期刊网 http://www.qikan.com.cn

式中

为输入码元的周期;α为滚降系数。记

，由式(2)可推出滚降系数为α

平方根升余弦冲击响应为： h(t)=Tπt(T2－16t2T2){T2sin［ 4αtTcos［ (3)

其时域响应如图1所示。在数字化波形成形时，为确保h(t)采样后的h［n］保持第一类线性相位，舍去所示。

图1 SRRC冲击响应时域波形 设发送端传递的二进制数据是示为： y(t)=…+a－

－

－

－

－

(4)

图2 成形滤波器成形过程 可以看出，当前传递信息［57..64］［17..24］

－

［49..56］,h－［9..16］与

时刻对应的波形信号的上升沿y［1..8］分别由h－［41..48］,h－

［33..40］如式(5)所示: ［41］+a－

［33］

［25］

［25..32］

［1..8］线性表示，

［49］+a－［1］ ［50］+a－［2］ ［51］+a－h［3］

［43］+a－

［35］

［27］

［42］+a－

［34］

［26］

－

－

－

－

－

－

－

样点，同时对N(偶数)点h［n］右移N/2。文中采用的h(t)中t∈［－

于是，采样后根升余弦成形滤波器的64个归一化h［n］如表1

4T,4T］，采样间隔为T/8，

］}

－α)t］+

，则发送滤波器的输出如图2所示，该波形函数可表

y［0］=a－

［17］ y［1］=a－

［18］ y［2］=a－

［19］

［57］+a－［9］［58］+a－［10］［59］+a－［11］

龙源期刊网 http://www.qikan.com.cn

y［3］=a－

［20］ y［4］=a－

［21］ y［5］=a－

［22］ y［6］=a－

［23］ y［7］=a－

［24］ (5)

表1 8T×8点平方根升余弦滤波器的h［n］ h［n］值h［n］值h［n］值h［n］值

1,64-0.000 002 1359,560.001 143 27917,48-0.013 952 50525,400.088 917 492 2,63-0.000 260 45510,550.006 630 61218,47-0.046 308 47626,390.234 615 884 0 3,62-0.000 809 61111,540.013 880 64819,46-0.082 303 06927,380.403 618 856 4,61-0.001 658 14212,530.021 653 66020,45-0.114 696 20328,370.580 185 884 0 5,60-0.002 659 331 113,520.027 860 07021,44-0.134 201 22929,360.745 562 348 6,59-0.003 461 19814,510.029 792 44222,43-0.130 815 71430,350.880 695 468 7,58-0.003 514 48215,500.024 643 70023,42-0.095 621 35731,340.969 196 958 8,57-0.002 167 52116,490.010 266 82524,41-0.022 744 88032,331.000 000 000 注：用图1离散后的次大值归一化处理。

2 二进制基带信号平方根升余弦成形滤波器的FPGA实现

［60］+a－［12］［61］+a－［13］［62］+a－［14］［63］+a－［15］［64］+a－［16］

［52］+a－［4］ ［53］+a－［5］ ［54］+a－［6］ ［55］+a－［7］ ［56］+a－［8］

［48］+a－

［40］

［32］

［47］+a－

［39］

［31］

［46］+a－

［38］

［30］

［45］+a－

［37］

［29］

［44］+a－

［36］

［28］

龙源期刊网 http://www.qikan.com.cn

在分析文献［2-4,6，8］的基础上，文中波形成形滤波器的实现采用的查表法结构如图3所示。其中，ROM单元存储待成形的数据与成形滤波器的冲击响应的卷积结果。模8计数器的工作时钟速率是待成形数据速率的8倍。待成形数据从8位移位寄存器的低位移入后，选择ROM表中的数据块data i，同时模8计数器C从

～

计数，并用该计数结果C(j)可采用双极性码，而查找

选择输出data i中的y［j］。当计数器C计数归零时，新的待成形数据从低位移入8位移位寄存器。该设计的一个优点是：ROM表中的数据在计算时，码，即输入信息符号序列{0,1}映射为{+1,－1}，持续时间不变。 图3 查表法波形成形滤波器实现结构 2.1 连续查找表法的一种改进实现

由图3不难实现采用本文的波形成形设计方案，共需要的ROM单元数目达2 048个。为此，可采用电路分割技术，将图3所示的8位移位寄存器输出的高8位地址同时给一个11位的中间寄存器，该中间寄存器的高8位又分为高4位与低4位，分别用于查找两个各具有16个单元的ROM表，之后再将各自的输出相加，此时消耗的ROM单元数共为256个。采用分割技术时，模8计数器、中间寄存器、ROM表三个部分的工作时钟相同。 2.2 根升余弦成形滤波器的VHDL实现

文中所述滤波器是在Modelsim 6.3d环境下采用VHDL实现的。Quarts环境中以文本方式调用LPM_ROM宏功能模块，定制ROM元件

－rom_16,元件的地址宽度分别是4 b，输

出数据字宽同为15 b。加法器的输出字宽16 b。实现的部分VHDL代码如下： Q1 Q2 Dout

L1:process(clk8x,cs) Ad1 Ad2 end process;

表地址产生电路使用单极性码。文中设计时，波形数据的计算采用了反逻辑、双极性、不归零

龙源期刊网 http://www.qikan.com.cn

u1: data_rom_16 port map(aclr=>clr, address=>Ad,clock=>clk8x,q=>Dout1); u1: data_rom_16 port map(aclr=>clr, address=>Ad,clock=>clk8x,q=>Dout2); 2.3 Modelsim时序仿真结果

Modelsim环境不能直接对mif格式的ROM初始化数据进行仿真，应在Quarts Ⅱ环境下先打开mif文件，再另存为hex格式，然后在Modelsim环境下编译后即可仿真。同时，如此操作又可将负值数据转为补码表示。Modelsim仿真结果如图4所示，其中clk的周期为160 ns，正好是一个din码元的宽度T，系统中地址产生电路的时钟周期是20 ns，以确保在一个码元持续时间内系统有8个样点输出。从图4中发现，一个码元成形后波形值延迟6T。 图4 根升余弦成形滤波器Modelsim仿真输出 3 结 论

文中所述的基于电路分割技术的查表法，实现通信系统发送端根升余弦滚降成形滤波器的FPGA实现方法简单可行，且当截断码元数目增多时或码内样点数目增加时，仅通过改变地址移位寄存器的长度或计数器的长度与ROM的长度即可，不至于使电路的复杂度成倍增加。 参 考 文 献

［1］KING Myung-Soon, CHUNG Jin-Gyun. Look-up table-based pulse-shaping filter ［J］. Electronics Letters, 2000, 36(17): 1505-1507.

［2］张会生,王效洪,耿光辉.基于FPGA实现根升余弦滤波器的研究

［3］郭继昌,李香萍,滕建辅.基于位串行分布式算法和FPGA实现FIR电路的研究［J］.电子测量与仪器学报，2001,15(2):15-21.

［4］杨贞斌,邬江兴.一种数字化基带成形实现新方法:基于码元间隔的定长滑动窗算法［J］.通信学报,1999,20(6):40-44.

［5］WHITE Stanley A. Application of distributed arithmetic to digital signal processing: a tutorial review ［J］. IEEE ASSP Magazine, 1989, 7: 4-18.

［6］王建新,蒋立平,吉训生,等.基带成形滤波器的FPGA实现［J］.电讯技术，2001,41(5):43-45.

龙源期刊网 http://www.qikan.com.cn

［7］张维良,郭兴波,潘长勇,等.平方根升余弦滚降FIR数字滤波器的设计［J］.电讯技术，2002,42(6):51-55.

［8］陈东华.升余弦滚降基带成型内插滤波器的FPGA实现［J］.华侨大学学报:自然科学版，2006,27(3):310-312.

［9］INRIA Alan Layec. Square raised root cosine filter block[EB/OL].[2007-07-25]. http://www-scicos.inria.fr/ScicosModNum/modnum_web/src/modnum_422/interf/scicos/help/eng/htm/SRRCF_c.htm.

［10］刘俊,刘会杰，尹增山.基于多速率DA的根升余弦滤波器的FPGA实现 注：本文中所涉及到的图表、注解、公式等内容请以PDF格式阅读原文