一种MSK-Type信号的频率同步改进算法

迟永钢， 郑宇希， 杨 木， 陈三斌

(哈尔滨工业大学 通信技术研究所，哈尔滨 150001)

一种MSK-Type信号的频率同步改进算法

迟永钢， 郑宇希， 杨 木， 陈三斌

(哈尔滨工业大学 通信技术研究所，哈尔滨 150001)

MSK-Type信号的频率同步技术可以采用基于最大似然理论的估计方法. 但该算法结构复杂，计算量十分庞大，不利用工程上的实现，所以对MSK-Type的频率同步大多采用非线性同步的方式. 其算法主要有延迟相乘和2P-Power算法，但这种选用前馈结构的无数据辅助算法只适用于全响应的连续相位调制信号，且由于在非线性相乘过程中引入大量自噪声，从而导致同步性能下降. 为解决MSK-Type调制信号频率同步算法精度低、适用范围小及不利于工程上实现等问题，本文在延迟相乘频率同步算法的基础上，提出一种改进算法. 该算法通过对部分响应和全响应MSK-Type信号采用不同幂次非线性处理，消除了延迟相乘算法中的自噪声影响，与原算法相比，频偏估计精度可提高5dB左右. 另外，所提算法不仅将应用范围扩展到了部分响应MSK-Type信号，扩大了同步算法的适用范围，而且还具有便于工程实现的特点.

数字信号处理；MSK-Type调制；载波同步；非线性频率同步算法；频偏估计

连续相位调制(CPM)是一种非线性有记忆的频率调制方式，由于CPM信号相位连续,所以该信号具有良好的频谱特性以及带外辐射小、旁瓣衰落快、瞬时包络恒定等优点，在抗非线性失真方面有良好的性能，因此，在卫星通信、机载舰载通信等民用和军用短波电台以及战术数据链中获得大量应用.

通常将调制系数为0.5的CPM信号称为MSK-Type调制. 在MSK-Type的解调中，同步技术是各系统中迫切需要解决的技术难点. R. Melhan针对MSK调制提出MCM算法[1]，但该算法只适用于MSK调制，对于采用其他成型滤波器的MSK-Type信号算法性能恶化严重，无法同步. Andrea在文献[2]中提出基于最大似然原理采用反馈结构的无数据辅助(NDA)算法，该算法利用劳伦特分解，将信号分解为线性PAM信号叠加的形式，再进行处理得到频偏信息. 最初，此算法仅适用于MSK信号，后扩展到一般的CPM信号[3]. 文献[4]针对多进制和部分响应信号提出基于延迟相乘的频率同步算法，该方法采用前馈结构，频偏捕获范围可达到符号速率的150%，且适用范围广，复杂度低. 但是牺牲了同步的精度，信号的自乘引入了大量的自噪声，其频偏估计的MSE曲线在信噪比为5 dB时达到10-3后便不随着信噪比的增加而衰减，所以算法同步性能很差. 在文献[5]中，MENGALI和ANDREA 对之前的同步算法进行了总结. M. Morelli在MCM算法的基础上进行改进[6]，针对MSK-Type调制提出了载波频率和定时联合同步的算法，将原MCM算法在频率同步中的精度提高了多个数量级，且改进的MCM算法适用于GMSK信号. 之后提出的一些频率同步算法包括国内学者在此方面的研究成果都存在着同步效率低或适用范围小等问题[7-10]. 所以寻找一种精度高，适用范围广的频率同步算法成为迫切需要.

在上述算法中，延迟相乘法易于实现且适用范围最广，基本涵盖所有的MSK-Type信号，所以本文从如何消除延迟相乘法中的自噪声入手，在已有的延迟相乘频率同步算法的基础上提出改进算法. 改进算法在延迟相乘之前先对接收信号进行幂次非线性处理，这样可以将调制信息从信号中移除，仅留下频偏信息，进而消除了自噪声的影响. 改进算法也解决了原算法因自噪声影响无法同步的问题，并将算法的精度提高了大约5个数量级. 算法中需要对全响应与部分响应的信号进行不同的幂次处理，可以有效保证算法的精度和广适性.

1 MSK-Type信号模型

接收信号r(t)由信号和噪声构成，具体表达式如下：

r(t)=s(t)+wR(t)+jwI(t)

其中wR(t)和wI(t)为高斯白噪声的实部和虚部，两者相互独立，功率谱密度均定义为N0. 发送信号s(t)的表达式为

(1)

式(1)为MSK-Type的通用表达式，其中Es为每个符号的能量；T为符号周期；θ为载波相位；τ为定时误差.s(t)中ψ(α;t)为相位信息，具体表达式为

).

下面给出g(t)常用的几种脉冲形状的表达式：

LREC：

LRC：

GMSK：

L=1且h=0.5的矩形脉冲对应的调制称为MSK调制；L>1或任一脉冲对应的调制称为MSK-Type调制. 此外，将L=1时的信号称为全响应信号，将L>1的信号称为部分响应信号. 本文只讨论g(t)是长度为L的矩形(LREC)和升余弦形(LRC)两种情况. 此时q(t)=1/2-q(LT-t)，满足几何对称的特性.

现假设处理信号的抗混叠滤波器为矩形，采样率为符号周期的N倍，抗混叠滤波器的带宽BAAF=N/2T. 则信号经过抗混叠滤波器后进行采样得到

(2)

2 频率同步算法

2.1延迟相乘法

接收信号可由调制信和高斯白噪声构成：

x(kTs)=s(kTs)+n(kTs).

利用定时同步算法中自相关的思想，忽略噪声的影响，对接收信号进行延迟相乘得

z(kTs)=x(kTs)x*[(k-D)Ts]=

其中D表示延迟的采样点数. 对于数据符号取期望得到

(3)

由于实际中期望不可获得，因此约等于

(4)

延迟相乘算法框图如图1所示.

图1 延迟相乘发算法框图

2.2基于延迟相乘法的改进算法

对信号进行延迟相乘可以消除信源和相位对频率估计的影响，但同时也引入了大量的自噪声，导致同步性能很差. 不论是全响应信号还是部分响应信号频偏估计误差的方差与克拉美罗界相差甚远. 所以为了消除自噪声，考虑将调制信息从信号中移除，只留下频偏信息. 而通过对指数信号进行幂次处理可以将调制信息最终变为±π的整数倍，进而消除掉调制信息. 本文采用先对信号进行幂次非线性处理，后进行延迟相乘的方式对信号进行同步. 由于L参数的影响，需要根据L的不同而采取不同的非线性处理方式，以使得算法适用于部分相应信号. 具体方法如下：

对于调制指数h=1/2和脉冲长度为L的MSK-Type信号，假设定时信息已经全部获得，采样率为符号速率，信号形式为

先考虑全响应信号的情况，即脉冲长度L=1时，对信号进行平方非线性处理得到

z(kT)=[x(kT)x*(kT-T)]2=

(5)

式中

因为此时脉冲长度L=1，所以在一个码元周期内，相位增加了π/2，因此上式结果为

2[ψ(kT,α)-ψ(kT-T,α)]=±π,

因此式(5)中最后一个指数项为-1，变为

下面再考虑部分响应信号的情况，例如当脉冲长度L=2时，由于L的影响，要对信号进行2L也就是4次方非线性处理便可得到

z(kT)=[x(kT)x*(kT-T)]4=

(6)

因为此时在一个码元周期内，相位增加了±2π，所以得到

4[ψ(kT,α)-ψ(kT-T,α)]=0,

因此式(5)中最后一个指数项为1，则式(6)可以写为

这样，综合公式(4)可得估计频偏的计算式为

其中

z(kT)=[x(kT)x*(kT-T)]2L.

改进后的算法框图如图2所示.

图2 改进算法的原理框图

由于该算法对信号实施了非线性变换处理，从算法推导中也可以看出，非线性处理弱化了噪声的影响，改进算法的同步性能更好. 但由于非线性变换的存在该算法的频偏估计范围必然会减小，非线性次数越高，估计范围越小.

3 性能分析

下面给出算法的性能分析，以及与现有算法的比较，频率估计的克拉美罗界为

3.1改进算法精度分析

所有仿真中观测长度均设为100个符号，蒙特卡罗循环为2 000次. 图3和图4分别针对多调制阶数和部分响应信号对延迟相乘法和改进算法进行了性能仿真，从图3和图4可以看出改进算法的精度比延迟相乘法高出很多. 由于延迟相乘法的自噪声非常大，在信噪比5～30 dB范围内，频偏估计误差的MSE没有随着信噪比的增大而衰减，而改进算法随着信噪比的增加，均方误差曲线衰减且精度越来越高，在20dB时的精度比延迟相乘法提高了4个数量级，在25 dB时的精度比延迟相乘法提高了5个数量级. 此外，从仿真中也可以看出改进算法也适用于部分响应的MSK-Type信号，进而证明了改进算法精度更高，涵盖范围更广. 具体仿真结果如下.

1)仿真条件：进制数M=2、4；脉冲长度L=1；脉冲形状REC；归一化频偏vT为10%T(T为符号速率)；过采样倍数N=4；抗混叠滤波器为8阶低通巴特沃斯滤波器，带宽Bw=0.75/T；Eb/N0=5～35 dB.

2)仿真条件：进制数M=2，4；脉冲长度L=1，2，脉冲形状为RC；其余条件同上.

图3 不同进制下频偏估计的方差与克拉美罗界对比

Fig.3 The frequency offset estimation variance compared with MCRB in differentM

图4不同进制及不同脉冲长度下频偏估计的方差与克拉美罗界对比

Fig.4 The frequency offset estimation variance compared with MCRB in differentMandL

3.2改进算法频偏捕获范围分析

仿真条件：信噪比Eb/N0=20 dB，脉冲长度L=1，脉冲形状为REC，其余条件不变. 频偏估计结果如图5所示.

图5 不同频偏下的估计范围对比

Fig.5 The comparison of estimation range under different frequency offsets

在频偏捕获范围方面，改进算法的估计范围为符号速率的25%，而延迟相乘法的估计范围为180%. 这是由于引入了幂次非线性变换，所以改进算法的频偏估计范围必然会减小，非线性次数越高，估计范围越小，仿真结果与理论一致.

3.3改进算法与改进的MCM算法比较

现有的MSK-Type信号频率同步算法中比较成熟的是MCM算法和改进的MCM算法. MCM算法仅适用于MSK调制，而改进的MCM算法虽然适用于非矩形脉冲成形的MSK-Type信号，但是对于L>1的部分响应信号无法同步. 针对这两种算法进行了比较.

从图6可以看出，在L=1时即对MSK的同步中，改进的MCM算法同步性能优于本文的改进算法，但是当L=2时，即在部分响应信号中，改进的MCM算法同步性能急剧恶化，无法同步,而本文的改进算法在信噪比为25dB时依旧可以衰减到10-8. 进一步证明了本文提出的改进算法的广适性和良好的同步性能.

图6 改进算法与改进的MCM算法的性能比较

Fig.6 The performance comparison of the improved algorithm and the improved MCM algorithm

4 结 语

在MSK-Type信号的频率同步算法中，基于非线性的延迟相乘法在相乘过程中引入大量自噪声导致同步性能不佳，估计方差不随信噪比的增大而衰减，且曲线距离克拉美罗界较远. 为了解决上述问题，本文在延迟相乘算法的基础上提出了改进算法. 根据信号频率脉冲归一化长度L的不同采取不同的幂次非线性处理，将调制信息从信号中移除，仅剩下频偏信息. 本算法消除了自噪声，解决了现有算法同步性能不佳的问题，将算法精度提高了约5个数量级. 同时改进后的算法在部分响应和多进制信号的同步中性能依旧良好.

[1] POLLOK A, MCkILLIAM R. Modified Cramér-Rao bounds for continuous-phase modulated signals[J]. IEEE Transactions on Communications, 2014, 62(5):1681-1690.

[2] MEHLAN R, CHEN Y E, MEYR H. A fully digital feedforward MSK demodulator with joint frequency offset and symbol timing estimation for burst mode mobile radio[J]. IEEE Transactions on Vehicular Technology, 1993, 42(4):434-443.

[3] D'Andrea A.N., GINESI A, MENGALI U. Frequency detectors for CPM signals[J]. IEEE Transactions on Communications, 1995, 43(2): 1828-1837.

[4] D'Andrea A.N., GINESI A, MENGALI U.Digital carrier frequency estimation for multilevel CPM signals[C]// IEEE International Conference on Communications. Seattle: IEEE, 1995:1041-1045.

[5] MENGALI U, D'Andrea A.N.. Synchronization techniques for digital receivers[M]. New York: Plenum Press, 1997.

[6] MICHELE M, UNBERTO M. Joint frequency and timing recovery for MSK-type modulation[J]. IEEE Transactions on Communications, 1999, 47(6): 938-946.

[7] GUDOVSKIY D A, CHU L, LEE S. A novel non-data aided synchronization algorithm for MSK-type-modulated signals[J]. Communications Letters IEEE, 2015, 19:1552-1555.

[8] 刘毅. CPM信号载波同步研究[D]. 西安：西安电子科技大学, 2008.

LIU Yi. Research on carrier synchronization of CPM signal [D]. Xi'an: Xi'an Electronic and Science University, 2008.

[9] 乔慧. 连续相位调制信号载波同步研究[D]. 西安：西安电子科技大学, 2012.

QIAO Hui. Research on carrier synchronization of continuous phase modulated signal [D]. Xi'an: Xi'an Electronic and Science University, 2012.

[10] 樊康敏. CPM收发机的关键技术研究[D]. 杭州：浙江大学, 2015.

FAN Kangmin. Key technology research on CPM transceiver [D]. Hangzhou: Zhejiang University, 2015.

AnimprovedfrequencysynchronizationalgorithmforMSK-Typesignals

CHI Yonggang, ZHENG Yuxi, YANG Mu, CHEN Sanbin

(Communication Research Center, Harbin Institute of Technology, Harbin 150001, China)

Frequency estimation algorithm for Minimum Shift Keying-Type(MSK-Type) signals is usually based on maximum likelihood theory. However, the algorithm is complex, and it is not practically realized in engineering. So the frequency synchronization for MSK-Type signal mostly adopts nonlinear methods such as delay-multiplication and 2P-Power algorithms. But these feedforward structure of non-data aided algorithm applies only to full response MSK-Type signal, and due to the nonlinear process of multiplying, these algorithms introduce large amounts of self-noise, so that it leads to poor synchronization performance. In order to solve the problem of lower accuracy and narrower application range and unpractical realization in engineering for the existing nonlinear algorithms, an improved algorithm based on delay-multiplication frequency synchronization algorithm is proposed in this paper. After different power nonlinear processing on the MSK-Type signal with partial response and full response, the improved algorithm eliminates the effects of self-noise in the delay-multiplication algorithm. Compared with the original algorithm, the estimated accuracy for the frequency offset is improved about 5dB, and at the same time, the application range is extended to the partial response MSK-Type signal. The improved algorithm enlarges the application range and it is easy to be implemented in engineering.

digital signal processing; MSK-Type modulation; carrier synchronization; nonlinear frequency synchronization algorithm; frequency offset estimation

10.11918/j.issn.0367-6234.201608021

TN911.72

A

0367-6234(2017)11-0036-05

2016-08-02

科技部重大科学仪器专项(MGZFSBQ13029)

迟永钢(1972—)，男，副教授

迟永钢，chiyg@hit.edu.cn

(编辑王小唯， 苗秀芝)