基于FPGA的北斗QPSK调制实现与解调验证

高 亮，宋茂忠

(南京航空航天大学电子信息工程学院，江苏南京 210016)

北斗卫星导航系统(BeiDou Navigation Satellite System)是我国正在实施的自主研发、完全独立运行的全球卫星导航系统，有着广泛的应用前景。北斗卫星导航系统信号采用正交相移键控(QPSK)调制，提高了数据传输速率，降低了信号间的相互干扰，改善了定位性能，成为目前全球卫星定位系统现代化发展的方向之一［1］。因此对北斗卫星导航接收机的需求日益增加，为了测试验收高性能卫星导航接收机的静态性能及动态性能，需要模拟产生北斗导航系统在各种环境下的真实卫星信号。

目前，针对北斗导航系统模拟信号源的研究较少，可参考GPS和GLONASS模拟信号源的研究，分析各个导航系统之间的差别，找到合适的研究方案。文献［2］分析了GLONASS信号的结构特点，研究了复杂环境下GLONASS导航信号的产生，文献［3］针对GPS信号模拟源的算法进行研究，并通过FPGA实现模拟源的产生，文献［4］分析了北斗卫星导航B1频段信号的结构，并用Simulink平台实现了信号模拟。

北斗导航系统已于2011年12月进入试运行阶段，并于2012年12月公布了空间信号接口控制文件［5］(Interface Control Document，ICD)。本文针对 ICD文件中北斗B1频段的QPSK调制信号进行结构分析，采用软件无线电的方法来研究QPSK信号的调制解调，分析QPSK信号调制解调的原理，利用FPGA仿真实现QPSK信号的调制解调，同时在射频输出端观察调制信号的功率密度谱。

1 北斗B1频段信号特征

北斗B1频段［6］信号由两个支路的“测距码+导航电文”正交调制在载波上构成。信号表达式如下

式中，j表示卫星编号;AI、AQ分别表示调制于B1频段载波同相、正交相支路的测距码振幅;CI、CQ分别表示同相、正交相支路的测距码;DI、DQ分别表示同相、正交相支路测距码上调制的数据码;f表示B1频段载波频率;φI、φQ分别表示B1频段载波同相、正交相支路的初相。

根据北斗B1频段信号QPSK调制方式的特点以及上式所示的信号结构，可知北斗B1频段每颗卫星均有一对独特的测距码，两者彼此不相关且正交，接收到同一颗卫星信号的两个支路受各种因素影响所导致的码延迟、多普勒频移等均相同。

2 北斗信号的QPSK调制

在“北斗二号”系统中，采用正交相移键控，QPSK，QPSK调制技术的核心思想是两个载波正交BPSK信号的合成，即将一个比特流的每一个四进制码元用两个二进制码元的组合来表示，两个二进制码元中的前一个码元用I表示，后一个码元用Q表示。由于在一个调制符号中传输2 Byte，QPSK调制比BPSK调制的带宽效率高一倍。载波的相位为4个间隔相等的值:π/4，3π/4，5π/4，7π/4。这种调制方式使同一载波每次可传输2 Byte信息，从而使载波的频带利用率比BPSK提高了1倍，同时抗干扰性更强。

图1 B1信号扩频码的自相关特性图

图2 QPSK调制后的自相关图

对QPSK调制后的数据进行相关运算，选取50位数据，通过仿真结果可以发现，加入QPSK调制技术后，原测距码的自相关曲线发生了变化［7］。在QPSK调制下的测距码，除了存在一个明显的高峰外，还具有两个比较明显的次高峰，分别位于靠近主峰的左右两侧。

3 北斗B1频段信号的模拟

北斗B1频段数字信号生成首先按照固定的采样时间间隔，也就是固定的采样频率(采样频率取1 561.098 MHz)，加入卫星的初始相位和多普勒频移，生成中频载波信号;然后将经过QPSK调制方式调制的测距码和北斗导航电文调制到载波上叠加。如果是MEO/IGSO的卫星需要在导航电文中在调制20位的NH码。每颗卫星的测距码和导航电文都不同，因此应该对每颗卫星的信号进行叠加，接着对多卫星信号的每颗卫星信号加入预设的强度，可设定的高斯白噪声，最后设定好模拟时长，将信号数据通过工作站储存模块导入到工作站保存。

图3 北斗B1频段信号模拟硬件功能框图

4 FPGA仿真实现QPSK调制解调

QPSK调制对输入数据流经内插、成形滤波和并串转换后，再经正交调制后即得到已调QPSK信号。内插有助于对基带信号进行波形成形，可以通过重复现存采样点或者插入零脉冲来实现，成形滤波器的作用是消除码间干扰和频谱扩散。QPSK解调对基带信号进行抽取、CIC滤波和Costas环载波同步后，再经过解调得到输出数据。采用多级CIC滤波器级联来实现较大的阻带衰减，Costas环来估计和校正解调过程中的多普勒频移［8－9］。在整个设计方案中，数字调制解调算法在FPGA中完成。FPGA硬件框图如图4所示，系统原理方案如图5所示。

图4 QPSK调制解调系统硬件框图

图5 QPSK调制解调原理框图

图5中CIC滤波器即级联积分器梳状(Cascade Integrator Comb，CIC)滤波器，滤波器的冲击响应如下形式

图6 CIC滤波器的幅频特性和相频特性(D=10)

图7 CIC滤波器的FPGA仿真

调制部分，数据通过串口输入，产生宽度为8 bit的并行数据流，然后分成I，Q两路数据流，经过补零内插，成形滤波和数据位截短后，通过并串转换，再通过DDS模块生成基带QPSK调制信号。调制输出数据的速率是128 kbit·s－1，而实现 DDS的器件 AD9857中设定输入数据速率是1 Mbit·s－1，为使两者速率匹配，设计中采用内插方法来实现。在内插模块的设计中，采用最简单的实现方法，即在数据之间插零。零的个数N由内插前后数据的速率决定，设计中N=7。内插模块通过数据锁存器和计数器实现。在内插模块设计中，由于采用了“插零”处理，导致码间干扰和带外辐射增大，为减小这些对信号解调的影响，设计采用成型滤波器。并串转换模块通过4级锁位寄存器实现，并行16位的输入数据按照并行4位的格式串行输入到DDS。

在解调部分，通过NCO进行数字下变频，经过抽取，CIC滤波器进行滤波，Costas环进行载波同步，最后解调处数据，通过串口输出显示。数字下变频模块主要由NCO和混频器组成。混频器为8位乘法器，采用Altera提供的宏功能模块LPM_MULT，将数据锁存模块输出的8位数据与分别NCO输出的正交载波进行相乘，乘积结果为一有符号的16位数据。为防止乘法器输出数据在后面模块运算处理中发生溢出，同样在设计中对数据进行了截短处理。CIC滤波器采用三级8倍抽取的积分滤波器与梳状滤波器并联。Costas环来估计和校正解调过程中的多普勒频移。

在测试部分，串口输入数据是1234FA42342343，串口输出12 34 FA 42 34 23 43，说明了QPSK调制解调系统的正确性。同时，通过频谱仪对设计的QPSK调制信号进行观察，FPGA开发板射频端口输出1 561.020 MHz的射频信号，证明生成的QPSK调制信号在频率特性、频谱特性上符合系统设计要求，从而证明了调制信号的正确性，功率密度谱如图8所示。

图8 功率密度谱图

5 结束语

北斗卫星导航系统，特别是B1频段信号，是未来民用的主要信号，因此，频段信号成为研究热点，而市场对北斗接收机的需求也会越来越多。为了测试北斗卫星导航接收机的性能，则需要模拟北斗卫星导航系统在各种环境下的真实卫星信号。本文针对北斗信号QPSK调制信号的结构特点，分析了QPSK调制解调的原理，并结合软件无线电的方法，通过FPGA验证了系统的可行性和正确性。

［1］《中国航天器》编辑部.中国的北斗导航卫星［J］.中国航天，2009(6):7－9.

［2］张鹏.复杂环境下GLONASS导航信号发生器的设计与实现［D］.南京:南京航空航天大学，2013.

［3］张敏虎，任章，华春红.惯性信息辅助的高动态弱GPS信号快速捕获［J］.系统工程与电子技术，2011，33(2):366－369.

［4］刘天旻.北斗卫星导航系统B1频段信号分析研究［D］.上海:上海交通大学，2013.

［5］中国卫星导航系统管理办公室.北斗卫星导航系统空间信号接口控制文件［S］.北京:中国卫星导航系统管理办公室，2012.

［6］黄隽祎，李荣冰，王翌，等.北斗B1 QPSK调制信号的高灵敏度捕获算法［J］.航空计算技术，2012，42(5):38－42.

［7］惠楠，黄海生，李鑫.一种北斗二号测距码生成器的设计［J］.西安邮电学院学报，2012，17(3):38 －47.

［8］姚彦鑫，杨东凯，丁凡，等.GPS码精密跟踪环路的研究［J］.武汉大学学报:信息科 学版，2010，35(11):1355－1359.

［9］崔冕，王宇.一种中高速突发QPSK数字解调系统捕获与载波同步算法的研究［J］.空间电子技术，2010，7(1):30－33，59.