一种跳频MSK信号检测算法及FPGA 实现

      引言

　　采用MSK 调制的跳频通信具有主瓣能量集中、旁瓣衰落滚降快、频谱利用率高和抗干扰能力强等优点，在军事通信中应用广泛。如美军现役的联合战术信息分发系统采用的通信信号，工作带宽969~1 206 MHz，跳频速率为70000 多跳/ s， 单个频点驻留时间约为13 s，信号持续时间* s， 总共有51个间隔为3 MHz 的信道，码速率为5 MHz。已知在该工作频段内主要还存在单频、窄带调幅和线性调频等信号。为了准确截获并识别目标信号，针对此信号环境设计了一种MSK 信号检测识别方法，并使用FPGA 进行了设计实现。

　　1 算法设计

　　1.1 宽带跳频信号实时检测算法

　　用现代技术来实现宽带数字化接收的一个实用的方法是通过信道化技术，实现信道化通常的方法是采用快速傅里叶变换（FFT）。利用FFT技术比用单个滤波器设计法更容易实现，因为FFT所需要的运算量更少。

　　某个由FFT运算输出的频率分量，可以看成输入信号与某个脉冲函数的卷积。因此可以把FFT的每个输出看成滤波器的脉冲响应函数与输入信号的卷积。为了处理一个连续的输入信号。必须在不同时刻对各段数据进行FFT处理。通常，起始点记为n = 0，数据段可以滑动M点，相应的FFT可以写成：

　 M 的值必须随着输入信号连续变化，这种运算也叫作短时傅里叶变换（STFT）。

　　FFT的长度和重叠点数是非常重要的参数，这些参数与最小脉宽和频率分辨率有关，它们决定了接收机的灵敏度。若FFT的长度为N，信号的采样频率为f s，那么经FFT计算后，信号的频率分辨率为：

　　数据重叠点数决定了时间分辨率和处理的最短脉宽，数据重叠率越高，则时间分辨率越高。

　　本设计中使用STFT 的方法实现一个粗测频引导数字接收机，为覆盖整个跳频带宽，采用700 MHz采样率对目标信号进行采样，粗测频引导精度在1 MHz以内，因此FFT长度选择为1 024 点，数据重叠率50%，保证时间分辨率在1 s 之内。

　　利用粗测频接收单元的检测结果去引导一个精测单元，可以对目标信号进行更为精确的测量和识别。精测单元采用数字正交下变频的实现方法，数字本振频率随粗测引导结果而设置，变频后的带宽依据目标信号而确定。

　　1.2 MSK 信号识别算法

　　MSK 信号可以写成：

　　式中，Tb 为码元周期，θk 是第k 个码元的相位常数，取值为nπ， Pk 为二进制双极性码元，取值为 1，所以MSK 信号相位分段线性变化，每个码元周期内相对前*元载波相位上升或下降π/2

　　根据三角函数展开，设θk 起始参考值为0，得：

　　又根据Ik= 1, Qk= 1, 令f L= f c-14Tb, f H=f c+14Tb,MSK 信号经过平方环可得：

　　可知，MSK 信号经过平方运算后，含有2f L 和2fH两个离散频率分量，反映在FFT幅度谱上，在这2 个频点上存在2 个明显的谱峰，2 个谱峰的距离为1 个码元速率，并且距离2 倍载频处均为码元速率的一半。这些特征是此频段内其他信号不具备的，可以根据这些特征对MSK 进行有效的识别。因此，将经过正交下变频的信号进行平方运算，确定信号出现后对其平方的结果进行FFT处理计算幅度谱，再对谱峰点进行分析即可完成MSK 信号的识别。

　　2 算法的FPGA 实现

　　2.1 算法实现

　　跳频MSK 信号检测识别的实现框图如图1 所示，采用FPGA 实现，包括短时傅里叶变换（STFT） 粗测频引导、数字正交下变频、平方运算和幅度谱分析等主要模块。STFT 粗测频引导在宽带条件下进行实时的信号检测和频率粗测，测量的结果引导数字正交下变频模块，对信号进行变频、滤波和抽取，得到低采样速率的零中频数据，平方运算模块对零中频数据进行平方处理，在确定存在信号后，对平方运算模块的输出进行FFT运算得到信号的幅度谱，通过幅度谱分析模块得到最终的识别结果。