基于fpga的fse-cma盲均衡器设计与实现

RESEARCHANDEXPLORATIONINLABORATORY

第38卷第10期Vol．38No．102019年10月Oct．2019

·专题研讨———虚拟仿真实验（74）·

CMA盲均衡器设计与实现基于FPGA的FSE-1，2

郭业才，李

11

峰，万逸儒，胡

峥

1

（1．南京信息工程大学电子与信息工程学院，南京210044；

2．南京信息工程大学滨江学院，江苏无锡214105）

要：针对传统常模盲均衡（CMA）算法相对复杂，太多乘法器导致占用过多逻

CMA盲均衡器。给出了盲均设计并实现了一种基于FPGA的FSE-辑资源的不足，摘

衡器的总体框架，设计了硬件电路。以1～20MHz的DPSK信号为输入信号，设计CMA的滤波器结构，了数字下变频模块和信号解调模块。改进了FSE-设计了硬件实现组成模块及抽头系数更新过程。搭建了硬件测试平台并进行级联了测试与

验证。

关键词：数字下变频；DPSK解调；滤波器结构改进；分数间隔均衡器中图分类号：TN911．7文献标志码：A文章编号：1006－7167（2019）10－0093－07

DesignandImplementationofFSE-CMABlindEqualizerBasedonFPGA

2

GUOYecai1，，LIFeng1，WANYiru1，HUZheng1

（1．CollegeofElectronicandInformationEngineering，NanjingUniversityofInformationScienceand

Technology，Nanjing210044，China；2．BingjiangCollege，NanjingUniversityofInformation

ScienceandTechnology，Wuxi214105，Jiangsu，China）

Abstract：InordertoovercometheshortcomingsoftraditionalCMAwhichhasrelativelycomplexstructureandtoomanymultipliers，andmayleadtotoomuchlogicresources，anFSE-CMAblindequalizerbasedonFPGAisdesignedandimplemented．Firstly，theoverallframeworkofblindequalizerisgivenandthehardwarecircuitisdesigned．Secondly，thedigitaldown-conversionmoduleandsignaldemodulationmodulearecreated，andittakestheDPSKsignalof1kHz～20MHzastheinputsignal．Thirdly，thefilterstructureofFSE-CMAisimproved，andthehardwaremoduleandtapcoefficientupdatingprocessareestablished．Finally，ahardwaretestplatformisbuiltforcascadetestingandverification．

Keywords：digitaldown-conversion；DPSKdemodulation；improvedfilterstructure；fractionallyspacedequalizer

0引言

近年来，随着无线通信技术的快速发展，盲自适应

［1］

均衡器的设计与实现研究成为重点关注的课题。

10-19收稿日期：2018-61371131），基金项目：国家自然科学基金（61673222，江苏省高校自然科学研究重大项目（13KJA510001）；江苏高校品牌专业建设项目（PPZY2015B134）；江苏省教育教学改革项目（2017JSJG168）

作者简介：郭业才（1962－），男，安徽安庆人，博士，教授，博导，主要研究方向：通信信号处理、水声信号处理，等。

Tel．：18861510036；E-mail：guo-yecai@163．com

2］文献［中设计的盲均衡器实现了Bussgang类算法中

G算法，3］文献［中设计的盲均衡器成功地用于的B-C，DVB-C）数字有线电视（DigitalVideoBroadcasting-全数字正交幅度调制（QuadratureAmplitudeModulation，QAM）接收机，4］中在深入分析文献［

Bussgang盲均衡算法及其特例基础上，给出了盲均衡系统的一种硬件并行实现结构并总结了硬件设计和实

5］中研究了分数间隔常模均文献［现中的重要问题，

ConstantModulus衡器（FractionallySpacedEqualization-Algorithm，FSE-CMA）的现场可编程门阵列（Field

PrograrnmableCateArray，FPGA）实现问题。上述盲均

94

衡器硬件实现成本高、占用逻辑资源多，难以用FPGA

直接进行硬件实现［6］

。为了克服所述不足，设计一款具有信号解调功能的盲均衡器硬件实现方法是十分必要的。

本文采用1kHz～20MHz的差分移相键控（DifferrentialPhaseShiftKeying，DPSK）信号作为输入信号，通过加入信号解调模块、改进分数间隔盲均衡（FSE-CMA）结构，以利FPGA硬件实现。在盲均衡器的FPGA硬件实现过程中，通过Matlab和modelsim仿真验证可行性。最后，通过搭建硬件测试平台进行级

联测试，

验证均衡效果。1FSE-CMA盲均衡器的总体设计

FSE-CMA盲均衡器的实现方法有多种，如基于软

件平台的实现和基于硬件平台的实现。本文采用基于

FPGA硬件平台来实现FSE-CMA盲均衡器。基于该平台的系统框架结构，如图1所示。

图1FSE-CMA盲均衡器总体系统框架图

加噪后的AD模拟信号由AD9236采样，将得到12位的数字信号输入至FPGA中［7］；在FPGA芯片内部对信号进行数字下变频处理，通过解调后得到基带

信号；再由FSE-CMA盲均衡模块对基带信号进行运算，运算结果通过DA芯片AD9744

［8］

输出。NiosII嵌

入式平台由QSYS工具搭建，其外部扩展了256MB同

步动态原储器（Synchronousdynamicrandomaccessmemory，SDRAM）、MB串行FLASH芯片EPCS。SDRAM用来临时存放数据［9］，EPCS用来存储本设计的所有算法程序

［10］

。

2硬件电路设计

FSE-CMA盲均衡器的硬件电路主要由电源、

FPGA、A/D转换、D/A转换等模块等构成，如图2

所示。

考虑到价格因素

［11］

，选用Cyclone系列的

EP3C16Q240C8N芯片。该芯片共有15408个逻辑单

元，

0．5MB的RAM，160个I/O控制引脚，20个全局时钟域，4个锁相环，56个乘法器，具有低成本、低功耗、高性能的特点

［12］

。FPGA模块有以下配置：系统时

钟由50MHz的有源晶振提供；存储模块分为带电可

第38卷

图2硬件电路模块连接图

擦可编程只读存储器（ElectricallyErasable

ProgrammableFeadonlymemory，EEPROM）模块和SDRAM模块，EEPROM芯片采用的是MB的EPCSI16N串行Flash，用于存储程序；SDRAM芯片采用的是MB的MT48LC4M16A2，用于临时存储图片、

字符等数据，方便FPGA调用；FPGA程序下载接口采用JTAG+PS配置两种模式，与JTAG+AS配置模式相比，省去了一个程序下载接口，节约了PCB空间；程序下载接口与EPCSI16N上的专用引脚连接，SDRAM的时钟引脚与FPGA芯片的全局时钟引脚连

接，由PLL核产生时钟信号，其余引脚均与I/O引脚连接。

电源模块与其余各个模块息息相关。A/D芯片

AD9236采用3．3V电压供电；D/A芯片AD9744采用3．3V电压供电，参考值电压1．2V；FPGA芯片采用3．3、2．5、1．2V3种电压供电。因此，电源模块使用AMS1117-3．3、AMS1117-2．5、AMS1117-1．23块稳压

芯片分别输出3．3、

2．5和1．2V电压。3软件设计

该盲均衡器的程序由数字下变频、

DPSK信号解调、FSE-CMA盲均衡算法3部分组成。3．1

数字下变频

数字下变频采用级联多阶CIC滤波器的方式来处

理。根据CIC滤波器原理，其实现结构可分为无反馈有限脉冲响应（FiniteImpulseResponse，FIR）结构及有反馈的无限脉冲响应（InfiniteImpulseResponse，IIR）结构。采用Noble恒等式原理，变换多级CIC滤

波器结构需要采用IIR滤波器结构［13］

，如图3所示。

图3

多级CIC滤波器结构

第10期

CMA盲均衡器设计与实现郭业才，等：基于FPGA的FSE-

95

在FPGA上设计抽取倍数为5的抽取系统，采用5

阶3级CIC滤波器，并对抽取系统进行仿真测试，系统输入数据位宽为10bit，系统时钟与数据频率相同。在设计该滤波器之前，还需解决滤波器运算过程中的字长问题。估算运算过程字长为

WI=Win+Nlb（MD）

（1）

式中：WI为CIC滤波器的中间字长；Win为输入数据位宽；N为CIC滤波器的阶数；M为系统抽取因子；D为滤波器的级数。

本设计中，滤波器阶数N=5，系统抽取因子M=5，滤波器级数D=3，则中间字长需要39bit。3．2

DPSK信号解调

本设计采用Costas环解调DPSK信号。Costas环

解调DPSK信号原理，

如图4所示，它是由输入信号分别乘以同相和正交两路载波信号而得名。输入信号的

上、下支路分别乘以同相和正交载波后，分别通过低通滤波器；两个低通滤波器输出再相乘，就完成了鉴相功能；最后，环路滤波器输出控制本地振荡器的误差电压。

图4Costas环解调DPSK信号原理

在设计Costas环之前，需要设计低通滤波器、数字鉴相器和环路滤波器。3．2．1

低通滤波器设计

低通滤波器的设计是载波同步环中的一个重点，关键在于如何确定滤波器的过渡带及截止频率。过渡带带宽的选择原则：必须确保滤除相邻的A/D镜像频

率成分；需要滤除数字下变频引入的倍频分量［14］

。根据采样定理，

A/D镜像频率的最小间隔公式为Δfad=min［

2f1－kfs，（k+1）fs－2fh］（2）

式中：f1为输入信号的下边缘频率；fh为输入信号的上边缘频率；fs为采样频率；min［］表示取小操作。本设计中，Δfad=4．8MHz。

数字下变频引入倍频分量的最低频率为fcddc=

min［－2f0+（m+1）fs，2f0－mfs］－Bf/2

（3）

式中：f0为载波频率，Bf为输入信号处理带宽。本设计中，

fcddc=8．4MHz。低通滤波器的截止频率为

fc=min［fcddc，Bf/2+Δfad］（4）

本设计中，式（4）中Bf=7．2MHz，故fc=8．4MHz。

由Matlab提供的凯撒窗函数进行设计，得到低通滤波器的频率特性，如图5所示。

图5低通滤波器频率特性图

3．2．2数字鉴相器设计

数字鉴相器实际上是同相支路与正交支路的相乘运算，在FPGA实现过程中，乘法运算需要消耗大量的硬件资源，而且运算速度也会受。若只取同相支路的符号位作为过零检测脉冲，并与正交支路进行异或运算，则会大大减少对硬件资源的消耗并提高运算速度。其实现如图6所示。

图6数字鉴相器设计图

3．2．3

环路滤波器设计

环路滤波器在Costas环中，不仅能低通滤波，而且

对环路参数调整起着决定作用［15］

。其结构如图7

所示。

图7

环路滤波器结构

由图7可得，数字化系统函数为：

C－1

F（z）=C2z1+

1－z－1

（5）

对于载波同步环路来说，环路中最重要的参数有：自然角频率ωn、阻尼系数ξ、环路总增益K。工程设计

96

第38卷

K通常取接近于1的值，中，ξ通常取0．707，它们之间

的关系为：

8ξ（S/N）iBiΔωL

＜ωn＜2ξ（S/N）L

（6）

0．7854。环路滤波器中C1、C2分别为

C1=2ξωnT/K

（8）

C2=（ωnT）2/K（9）

代入已知量到式（8）、（9）中，计算得C1

=0．002=2－9，C2=0．0625=2－4。3．2．4

Costas环的FPGA实现

Costas环的RTL图，如图8所示。图8表明，本设

2个乘法器、2个低通滤波器和一计Costas环有：DDS、

个鉴相及环路滤波器等模块组成。DDS模块由NCO

核生成；乘法器模块由乘法器核直接生成，两路输入分别为8bit和10bit有符号数，输出保留所有有效位；低通滤波器由FIR核生成，其参数为上述设计的低通滤波器参数。

（S/N）L＞6，i为输入。且要求ωn/fs＜0．1、

fs=16得，70．72rad·kHz＜设计中，由Bi=12，ωn＜1053rad·kHz。又根据ωn/fs＜0．1，得，ωn＜1．6×10rad/s=254．65rad·kHz。为兼顾稳态相差及快捕带宽，本设计中ωn=150rad·kHz、滤波器长度为

38、系数绝对值之和为768，这时，滤波器输出需增加13bit位宽。而锁相环路总增益为

2πfs

K=NTdds（2Bloop－2）

2

6

（7）

可容易计算出，符合本设计的环路总增益为

图8Costas环的RTL图

3．33．3．1

FSE-CMA的FPGA实现滤波器结构改进

［16］

常规型滤波器结构如图9所示。图9中，粗线

表示该滤波器的延迟路径，延迟路径上的运算单元包括1个乘法器和N－1个加法器，如8阶滤波器有1个乘法器和7个加法器。可见，常规滤波器结构的计算

这严重影响均衡器延迟随滤波器阶数的增加而增加，效果。

图10

图9

常规型滤波器结构

改进后的8阶改进型滤波器结构

3．3．2FSE-CMA的组成模块

将常规型滤波器结构改造为具有流水线式高速结

构，就除去了关键路径上的单位延时。改进后的8阶滤波器结构，如图10所示。图10表明，改进的8阶滤波器的延迟路径上有1个乘法器，但只有4个加法器；与8阶常规滤波器结构相比，减少了3个加法器。也就是说，对于N阶滤波器而言，改进后加法器数量为logN，大大缩短了计算所消耗的时间。

FSE-CMA的组成模块，如图11所示。

（1）总控模块。该模块控制系统的工作状态，是

系统的核心，主要为各模块提供全局控制信号、初始化各个模块并启动均衡器。

（2）输入信号延时模块。该模块对输入信号进行延时存储；在输入时钟上升沿将输入信号向量输入到输出计算模块和误差计算模块中。

第10期

CMA盲均衡器设计与实现郭业才，等：基于FPGA的FSE-97

图11FSE-CMA的组成

（3）输出计算模块。该模块由乘法器、加法器、平

方器和截位器组成，由其对输入信号向量与相应子均衡器的各个抽头系数进行乘法、加法、平方和截位运算，以得到均衡器输出信号z（n）。（4）误差计算模块。该模块用于计算FSE-CMA

更新方程f2

n+1=fn－μz（n）r（n）（|z（n）|－1）中误差

信号向量μz（n）r（n）（|z（n）|2

－1）。计算过程如下：均衡器输出z（n）由平方器计算|z（n）|2

并减1，

得到|z（n）|2

－1，同时均衡器输入信号向量r（n）通过乘法器与z（n）相乘得到向量z（n）r（n）；再将|z（n）|

2－1与向量z（n）r（n）相乘，得到z（n）r（n）（|z（n）|2

－

1）后乘以步长因子μ；最终，得到误差信号向量

μz（n）r（n）（|z（n）|2

－1）。该计算无须消耗时钟，简

化了设计，

提高了系统的处理速度。（5）抽头系数调整模块。该模块完成均衡器抽头

系数的更新运算。均衡器的抽头系数向量fn存储在寄存器中，误差计算模块输出误差向量信号

μz（n）r（n）（|z（n）|2

－1）与寄存器中当前抽头系数向量fn对应相减，得到新的均衡器抽头向量系数fn+1。

待到下一个输入时钟上升沿被触发，fn+1与新的输入信号向量rn+1相乘，进行下一次迭代过程。

3．3．3FSE-CMA的抽头系数更新过程

T/2分数间隔均衡器的实现框图，如图12所示。

图12

T/2间隔的FSE-CMA实现框图

计算步骤如下：

步骤1两路向量r（0）（1）

n和rn分别与子均衡器抽

头系数向量f（0）

n和f（1）

n的各个分量对应相乘、累加输

出z（n）；

步骤2

均衡器输出z（n）平方后和常数1相减，

得中间变量A。与此同时，均衡器输入向量rn的每一个分量与z（n）相乘，

得中间变量B；步骤3A、B相乘后，再乘步长因子μ，得中间变

量C，这里取2的整数次方；步骤4

C与抽头系数fn的每一个分量对应相减，得新抽头向量fn+1；步骤5新抽头系数向量再与下一个均衡器输入

向量重复进行步骤1～4操作，

如此循环就实现了均衡过程。

3．3．4FSE-CMA盲均衡器软件总体系统

该盲均衡器软件总体系统接口，如图13所示。图中，Clk表示模块时钟信号接口，与PLL提供的时钟信号相连接；Reset表示模块复位信号接口；Enable表示模块使能信号接口，因为该模块不间断工作，因此该接口维持高电平以持续工作；Data_in［13．．0］表示输入信号接口，与DPSK信号解调模块的输出端相连；Data_out［13．．0］表示输出信号接口，与D/A模块相连。

图13

盲均衡器总体系统接口图

4

测试平台与测试分析

4．1

测试平台

该盲均衡器的硬件平台，如图14所示。

图14盲均衡器硬件平台

测试平台搭建效果，如图15所示。该平台由高阶

调制信号发生器、卫星信道模拟器、盲均衡器和示波器

组成。其中，高阶调制信号发生器、卫星信道模拟器均由本团队研发，示波器采用DS4052。4．2

测试与结果

该盲均衡器编译报表，如图16所示。

98

第38卷

FPGA完成FSE-CMA盲均衡器的计图16表明，

算需要消耗4342个逻辑资源、占用2877个寄存器、占用73088bit存储容量。

DPSK信号解调模块仿真，df如图17所示。图中，

di为同相支路输出结果，dq为正交支表示收敛情况，

路输出结果。

FSE-CMA盲均衡器模块仿真结果如图18所示。

Xin为输入信号，clk为FPGA整体时钟信号，clk_data

图15

测试平台搭建

CMA盲均衡模块时钟信号，rst为复位信号，为FSE-Error为误差收敛信号。图19显示了FPGA仿真和Matlab仿真的误差收敛情况。该图表明，两种仿真的误差收敛速度基本相同。也就是说，基于FPGA的FSE-CMA实现结果能满足要求。实际测试结果，如图20所示。使用信号发生器产

如图20（a）所示，将此信号生600kHz的DPSK信号，

输入卫星信道模拟器。模拟器产生50dB的高斯白噪

声，如图20（b）所示，随后将信号输出至盲均衡器。盲均衡器经过解调和均衡后，得到输出波形，如图20（c）所示。

图16盲均衡器编译报表

图17DPSK信号解调模块仿真结果

图18FSE-CMA模块仿真结果

图19FPGA仿真与Matlab仿真误差收敛情况

（a）输入DPSK信号（b）加噪后的DPSK信号（c）解调、均衡后的DPSK信号

图20输入DPSK信号实际测试结果

第10期

CMA盲均衡器设计与实现郭业才，等：基于FPGA的FSE-

99

5结语

［6］任全会，杨保海．基于FPGA弱小信号测量系统研究与实现［J］．实验室研究与探索，2018，37（3）：92-96．

焱，周圣泽，罗

军，等．基于FPGA的并串转换电路硬件

CMA盲均衡器的硬件实现问题。本文研究了FSE-CMA盲在研究中给出了DPSK信号解调模块和FSE-均衡算法模块的设计思路。以1kHz～20MHz的

DPSK信号作为输入信号，设计了盲均衡器硬件结构和PCB结构；利用VerilogHDL语言进行了实现，编译成功后得到了完成该实现所需占用的逻辑资源数和储存空间；将Matlab和modelsim仿结合，给出了两种仿真结果；通过高阶调制信号发生器、卫星信道模拟器以及示波器级联测试，验证了均衡效果。参考文献（References）：

［1］李玉龙．无线通信系统中自适应均衡技术的研究［D］．贵阳：贵

2008．州大学，

［2］张益萍．基于FPGA的盲均衡器的设计［D］．西安：西北工业大

2004．学，

［3］许［4］夏［5］李

玲，蒋文军．一种用于数字QAM接收机的盲均衡器实现．成都：电威．Bussgang盲均衡算法研究及FPGA实现［D］

松．盲均衡器的FPGA实现［D］．郑州：信息工程大［J］．电视技术，2003，258（12）：15-17．2005．子科技大学，2004．学，

［7］刘

2017，43（12）：21-24．实现［J］．电子技术应用，

［8］ReederR．

Transformer-CoupledfrontendforwidebandA/D

converters［J］．NumberAnalogDevices，2007（114）：1–4．

［9］郭业才，万逸儒，彭舒，等．卫星信道模拟器触摸屏参数控制系

J］．实验室研究与探索，2018，37（2）：66-70．统的设计［

［10］朱宏光，司静．一种用于高速数据采集的SDRAM控制器［J］．

2010，40（4）：62-．无线电工程，

［11］郭业才．群体智能与计算智能优化的盲均衡算法［M］．北京：清

2018．华大学出版社，

［12］张思齐．基于三种不同FPGA芯片的数字中频接收机设计及实

D］．成都：西南交通大学，2009．现［

［13］张怀武．现代印刷电路原理与工艺［M］．2版．北京：机械工业出

2015．版社，

［14］谢海霞，孙志雄．多级CIC滤波器的FPGA实现［J］．电子设计

2016，24（23）：183-185．工程，

［15］杜

J］．广勇，刘帝英．基于FPGA的相乘微分型AFC环设计［2013（1）：59-62．东通信技术，

［16］张安安，杜勇，韩方景．全数字Costas环在FPGA上的设计与实

J］．信息化研究，2006，32（1）：18-20．现［

［17］荣雅君，杨秋霞．自适应滤波器的设计及其应用［J］．河北大学

2003，23（3）：307-311．学报（自然科学版），

檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿（上接第49页）参考文献（References）：

［1］JiangS，CaoD，LiY，etal．Low-THDfast-transientandcost-effectivesynchronous-framerepetitivecontrollerforthree-phaseUPSinverters［J］．IEEETransactionsonPowerElectronics，2012，27（6）：2994-3005．

［2］ZhangK，KangY，XiongJ，etal．DirectrepetitivecontrolofSPWM

inverterforUPSpurpose［J］．Electronics，2003，18（3）：784-792．

［3］HasanKomurcugil．Combineduseofdouble-bandhysteresiscurrent

andproportionalresonantcontrolmethodsforsingle-phaseUPSinverters［C］//IECON201440thAnnualConferenceoftheIEEEIndustrialElectronicsSociety．1305-1311．

［4］DongD，TimothyThacker，IgorCvetkovic，etal．Modesofoperation

andsystem-levelcontrolofsingle-phasebidirectionalPWMconverter．IEEETransactionsonSmartGrid，2012，3formicrogridsystems［J］（1）：93-104．

［5］BommegowdaK．B，VishwasK，SuryanarayanaK，etal．Single

phaseinvertercontrolwithcapacitorcurrentfeedback［C］//2015InternationalConferenceonPowerandAdvancedControlEngineering（ICPACE）．Bangalore，India：［s．n．］，2015：418-422．

［6］RyanMJ，LorenzRD．Ahighperformancesinewaveinverter

controllerwithcapacitorcurrentfeedbackand“back-EMF”decoupling［C］//26thAnnualIEEEPowerElectronicsSpecialists

Dallas，USA：［s．n．］，2014：IEEETransactionsonPower

，1995：507-513．Conference．Atlanta，USA：［s．n．］

［7］QiY，PengL，ChenML，etal．Loaddisturbancesuppressionfor

voltage-controlledthree-phasevoltagesourceinverterwithtransformer［J］．IETPowerElectronics，2014，7（12）：3147–3158．［8］林茜．基于GaN器件的单相逆变器数字控制技术研究［D］．哈

2017．尔滨：哈尔滨工程大学，

［9］TurnerRobert，WaltonSimon，DukeRichard．

Robusthigh-performanceinvertercontrolusingdiscretedirect-designpoleplacement［J］．IEEETransactionsonIndustrialElectronics，2011，58（1）：348-357．

［10］邱关源．电路［M］．北京：高等教育出版，2006．

［11］沈定坤，王卫星，徐军，等．基于坐标变换的单相光伏并网逆变器

J］．电力电子技术，2013，47（1）：55-57．研制［

［12］MonfaredM，GolestanS，GuerreroJM．Analysis，design，and

experimentalverificationofasynchronousreferenceframevoltagecontrolforsingle-phaseinverters［J］．258-269．

［13］李翠艳，．电机与控制学报，张东纯，庄显义．重复控制综述［J］

2005，9（1）：37-44．

［14］WuM，XuB，CaoW，etal．Aperiodicdisturbancerejectionin

repetitive-controlsystems［J］．IEEETransactionsonControlSystems2014，22（3）：1044-1051．Technology，

［15］马海啸，龚春英．负载电流前馈双闭环控制逆变器的研究［J］．

2012，32（3）：30-35．南京邮电大学学报（自然科学版），

IEEETransactionson

IndustrialElectronicsinIndustrialElectronics，2014，61（1）：