0 引言
如何对射频(RF)信号进行数字化是软件无线电(SDR)面临的首要问题，即如何对感兴趣的信号进行采样[1]。带通采样可以使用比RF信号最高频率低得多的采样率，因此被广泛用于SDR接收机中。随着无线通信的发展，多频段接收机[2-4]应用更加广泛。但是当接收到多频信号时，多频信号之间容易发生混叠。大多数研究人员通过选择采样频率来避免混叠问题，这限制了采样频率的选择[5-7]。二阶带通采样(BPS)可以用来降低采样率和处理混叠。对于二阶BPS系统，大部分研究都基于正交信号处理来进行抗混叠。例如文献[8]采样后采用希尔伯特变换，而文献[9]采用90°混合交叉信号进行采样。在文献[10]中，两个采样装置产生90°的相位差。这些研究允许采样系统使用较小的采样频率，但只能消除信号自身图像的混叠，并不能处理两个信号之间的混叠。改进的带通采样方法也可以实现对混叠信号的处理，但相应的重构算法较为复杂[11-13]。
在以前的工作中，二阶BPS可调用延时来实现两个重叠信号的接收[14]。本文提出延时可调的三阶带通采样结构用于处理3路带通信号混叠问题。通过对三阶带通采样信号频谱的分析，根据3路采样信号间相位差设计抗混叠滤波器，消除混叠信号，并通过仿真验证系统性能。
1 三阶BPS的结构
假设带通信号R(f)的带宽限制为B，采样率为fs=2B。定义信号位置索引n由式(1)表示，在索引位置n处的任意信号经带通采样后都可以恢复到第一奈奎斯特区中，其中第一奈奎斯区域|f|<B是索引为零的频率区。带宽B意味着处理带宽，即在该频带中可以存在多个的信号。
对于射频信号经带通采样后在同一频域内会引起三信号的混叠现象，设置三阶抗混叠滤波器来处理。
三阶BPS采样系统的结构如图1所示。
如图2所示，具有频谱R1(f)、R2(f)和R3(f)的3个RF信号S1、S2和S3，使用图1的结构进行三阶BPS采样。3个采样通道分别表示为通道A、通道B和通道C，将通道A中的频谱定义为：
在采样通道B中，频谱为：
在采样通道C中，频谱为：
其中，
频谱如图5所示。
2 抗混叠滤波器设计
在3个通道中分别定义抗混叠滤波器SA(f)、SB(f)和SC(f)，恢复信号的频谱为：
负频谱和正频谱是对称的，这里只讨论正频谱。使用式(3)和式(4)，式(6)和式(7)可以简化为：
将式(8)更改为矩阵形式：
基于式(10)～式(12)，可以通过对其时域数据进行采样作为FIR滤波器的系数来设计抗混叠滤波器。
由式(10)～式(12)可以看出，如果固定采样频率和三阶带通采样间采样延时，抗混叠滤波器只与信号位置索引n有关，因此对于不同位置的射频信号，只需要调整滤波器中系数就可以实现抗混叠效果。
3 仿真结果
基于抗混叠算法的分析，在MATLAB中使用FDA工具设计了抗混叠滤波器。测试输入3个信号S1、S2和S3，中心频率分别为fc1=1.320 GHz、fc2=2.025 GHz和fc3=1.630 GHz。采样频率选取fs=100 MHz。通过信号位置索引的定义可以得到n1=13，n2=20，n3=16。3个信号的带宽都为10 MHz。
在三阶BPS之后，3个信号同时被接收，频谱如图6所示。在BPS之后，3个信号的中心频率分别为20 MHz、25 MHz和30 MHz。从图6可以看出，这3个带通信号在采样后相互重叠。
根据式(10)～式(12)，应用三阶BPS在采样通道B和采样通道C中添加时间延迟TΔ1=0.4×10-9 s、TΔ2=1.2×10-9 s来设计抗混叠滤波器，信号S1的频谱恢复如图7所示。即由信号S2和信号S3引起的混叠被抑制。
使用相同的方法设计抗混叠滤波器来恢复S2和S3，结果如图8和图9所示。从图7～图9可以看出，混叠均被抑制超过40 dBW。
4 重构性能
对抗混叠处理后的信号进行下变频及解调，图10为中心频率为1.32 GHz 时重构信号的星座图，可算出信噪比为34.399 dB。
为了进一步验证其他频率范围内的重构性能，对频率偏置位于0～50 MHz（0～fs/2）范围内的信号进行测试，以n=20为例，根据一系列实验仿真结果得到相应的抗混叠性能分析如图11所示。结果表明接收信号信噪比在28 dB以上。
5 结论
本文基于BPS，给出了三阶BPS的结构。通过设计抗混叠滤波器，可以分离出3个重叠的信号。仿真结果表明，设计的抗混叠滤波器可以有效消除其他两种信号引起的混叠现象。从信号重建分析出其抗混叠性能很好。
参考文献
[1] MACHADO R G，WYGLINSKI A M.Software-defined radio：bridging the analog-digital divide[J].Proceedings of the IEEE，2015，103(3)：409-423.
[2] 李艳霞.射频直接数字化接收机[J].电子技术与软件工程，2017(16)：92-92.
[3] 殷楠，李玉峰，冯宪周.DME应答接收机技术研究[J].电子技术应用，2018，44(5)：27-30.
[4] 唐然，吴虹，程树军，等.带通采样星载AIS非相干接收机的FPGA实现[J].电子技术应用，2018，44(1)：33-36.
[5] MESSAOUD M A，GHANNOUCHI F M，BARRAK R，et al.SDR based multi-band subsampling receivers for GNSS applications[C].1st International Conference on Advanced Technologies for Signal and Image Processing，Sousse，TUNISIA，2014.
[6] THABET J，BARRAK R，GHAZEL A.A new reconfigurable architecture for Multistandard GNSS subsampling receiver[C].4th International Conference on Communications and Networking，Hammamet，TUNISIA，2014.
[7] THABET J，BARRAK R，GHAZEL A.Enhancement of bandpass sampling efficiency in direct RF subsampling receivers：application to multiband GPS subsampling receiver[C].International Conference on Multimedia Computing and Systems，Marrakech，MOROCCO，2014.
[8] AL-ABOODI M，LAMI I A.OCBPSR：orthogonal complex bandpass sampling receiver[C].World Congress on Computer Applications and Information Systems，Hammamet，TUNISIA，2014.
[9] Cao Minghua，Li Jianqiang，Dai Yitang，et al.Multiband phase-modulated RoF link with coherent detection and bandpass sampling[J].IEEE Photonics Technology Letters，2015，27(21)：2308-2311.
[10] JE-IN B.Design of second-order BPS systems for the cancellation of multiple aliasing[J].Journal of the Institute of Electronics and Information Engineers，2015，52(3)：162-170.
[11] 张睿，李维英，李建东.带通采样技术在软件接收机中的应用[J].西安电子科技大学学报，2000(3):326-329，367.
[12] 顾嘉辉，刘洋，李振华.高速采样信号并行数字下变频及FIR滤波方法[J].空间电子技术，2017，14(4)：22-27，38.
[13] Liu Gaohui，GaoYong，Yu Ningmei.Research on non uniform bandpass sampling[J].Journal of Beijing University of Posts And Telecom，2005，28(4)：66-68.
[14] Wang Hongmei，Yao Chong，Wang Faguang，et al.Multistandard receiver design for telemedicine monitoring system[J].Journal of Sensors，2018(4)：1-8.
作者信息:
王洪梅，姚 冲，王法广，李世银，宋金玲
(中国矿业大学 信息与控制学院，江苏 徐州221200)
