一种面向GMSK突发通信系统的快速捕获与同步方法与流程

一种面向gmsk突发通信系统的快速捕获与同步方法
技术领域
1.本发明涉及无线通信技术领域，具体涉及一种面向gmsk突发通信系统的快速捕获与同步方法。


背景技术：

2.在通信过程中，传输时间越长，则越容易被敌方定位，信号的抗截获能力就越差，而gmsk突发通信能够尽可地避免这一问题。
3.gmsk突发通信的原理是在无线通信双方的信息收发过程中，发射机在极短的时间内，采用高斯滤波最小频移键控(gmsk，gaussian filtered minimum shift keying) 调制产生信号，并按照一定机制发送较短的数据帧，以此来提高通信系统的抗干扰性和信号的抗截获能力。在通信系统中，通讯方总是希望能将传输的信息隐藏起来，以达到不被侦测的目的。在常见的扩频和跳频系统中，虽然都可以提高通信系统的抗干扰性，但在这过程中都增加了信号的带宽。而gmsk突发通信只是减少了信号的传输时间，并没有增加信号的传输带宽，同时提高了系统的稳定性。因此，突发模式都得到了较为广泛的应用，如卫星通信系统、时分多址通信系统以及一些短波信道中的串行通信系统等。
4.gmsk突发通信的发送方式对接收机能够准确无误地恢复出原始数据提出了一定的要求：由于gmsk突发信号持续时间极短，因此接收机对捕获和同步的建立过程也要尽可能的短，否则会导致数据的缺失，最终无法恢复出原始数据信息。
5.传统的突发信号的捕获利用本地的训练序列与接收信号直接进行滑动相关检测，该方法的不足之处有：
6.1.计算量大，硬件运算复杂度高。由于逐点滑动相关处理，在降采样后需要对每个采样点进行滑动相关，计算量大，硬件运算复杂度高；
7.2.在低信噪比(如-8db)下相关检测无法正确检测出峰值，其初始同步位置处的相关值在低信噪比时已淹没在相关检测结果中；
8.3.无法抵抗较大频偏，较大频偏时容易出现较高的误捕或漏捕概率。


技术实现要素：

9.本发明的目的在于提供一种面向gmsk突发通信系统的快速捕获与同步方法，用以解决现有gmsk突发通信系统中存在较大频偏难估计、估计精度较低，且gmsk 突发无线通信系统捕获效率有待提高的问题。
10.为实现上述目的，本发明的技术方案为：
11.一种面向gmsk突发通信系统的快速捕获与同步方法，包括以下步骤：
12.步骤1：利用地面终端设备接收来自发射端发射后的突发信号，并进行低通滤波以及下采样；
13.步骤2：将经过低通滤波以及下采样后的突发信号分别向后取若干个样点，分别得到两段信号，分别是第一段信号、第二段信号，将两段信号进行共轭相关计算，并得到共轭
值；
14.步骤3：将步骤2得到的共轭值与相关门限值进行比较，若共轭值超过相关门限值，则标记当前相关次数，并进行步骤4，否则，继续进行步骤2的相关计算，并再次与相关门限值进行比较；
15.步骤4：利用超过相关门限值时的信号进行频偏估计，并得到频偏估计结果，然后将频偏估计结果补偿到第二段信号；
16.步骤5：将频偏补偿后的信号与本地导频信号分别进行fft计算，并将分别的计算结果进行共轭相乘；
17.步骤6：将步骤5得到的共轭相乘结果进行ifft计算，并取模值平方，求取该模值平方的最大值，并将该最大值与ifft门限进行比较，当最大值超过ifft门限时，标记当前索引值，并进行步骤7，否则重复步骤2-6；
18.步骤7：利用步骤3中标记的相关次数与步骤6得到的索引值，计算同步头的起始位置，完成同步捕获。
19.进一步地，作为优选技术方案，还包括同步头结构设计步骤：将无线突发通信系统同步头设计为三段训练序列，每段导频序列长度为l
sym
个符号，每段导频由n
bit
比特的m序列经过高斯滤波最小频移键控调制产生，然后级联为一帧，同步头内的三段导频符号完全相同，在进行无线突发通信同步捕获时，接收端存储一段导频符号，该段导频符号表示为：s
loc
(m)，1≤m≤l
sym
。
20.进一步地，作为优选技术方案，所述步骤1中，终端实际的接收信号为 r(m),1≤m≤lr，其中lr为接收到的信号的总长度。
21.进一步地，作为优选技术方案，所述步骤2的具体过程为：
22.步骤2-1：共轭相关计算每次滑动步进的步长为d
step
个采样点的数据 (1≤d
step
≤l
sym
)，且在第n次共轭相关计算时，已滑动距首个样点(n-1)
·dstep
长度的距离；
23.步骤2-2：将经过低通滤波以及下采样后的突发信号从m＝(n-1)
·dstep
+1开始向后取l
sym
个样点的信号，得到第一段信号，该段信号表示为r
first,n
(m),1≤m≤l
sym
；
24.步骤2-3：将经过低通滤波以及下采样后的突发信号从m＝(n-1)
·dstep
+l
sym
+1开始向后取l
sym
个样点的信号，得到第二段信号，该段信号表示为 r
second,n
(m),1≤m≤l
sym
；
25.步骤2-4：将第一段信号与第二段信号做共轭相关计算，具体为：
[0026][0027]
其中，(
·
)
*
表示求取共轭值。
[0028]
进一步地，作为优选技术方案，所述步骤3中，当共轭值超过相关门限值时，标记当前相关次数n。
[0029]
进一步地，作为优选技术方案，所述步骤4的具体过程为：
[0030]
步骤4-1：利用当前超过相关门限值时的信号r
first,n
(m),1≤m≤l
sym
、 r
second,n
(m),1≤m≤l
sym
进行频偏估计，频偏估计具体表示为：
[0031]
[0032][0033]
其中，(
·
)
*
表示求取共轭值，fs为降采样后的信号频率；arg{x}为求得x的相位角，为频偏估计的结果；
[0034]
步骤4-2：将频偏估计的结果补偿到第二段信号r
second,n
(m),1≤m≤l
sym
，得到如下信号：
[0035][0036]
进一步地，作为优选技术方案，所述步骤5的具体过程为：
[0037]
步骤5-1：频偏补偿后的信号与本地保存的导频信号 s
loc
(m)，1≤m≤l
sym
分别做l
sym
点的fft计算，表示如下：
[0038][0039]
s(m)＝fft(s
loc
(m))，1≤m≤l
sym
；
[0040]
步骤5-2：将步骤5-1得到的fft计算结果进行共轭相乘，表示如下：
[0041]cfft
(m)＝r(m)
·s*
(m),1≤m≤l
sym
[0042]
其中，fft(
·
)表示为fft计算。
[0043]
进一步地，作为优选技术方案，所述步骤6的具体过程为：
[0044]
步骤6-1：将步骤5-1得到的共轭相乘结果c
fft
进行l
sym
点的ifft计算，并取模值平方，表示如下：
[0045]cifft
(m)＝|ifft(c
fft
(m))2,1≤m≤l
sym
[0046]
其中，ifft(
·
)表示为ifft计算，|
·
|为取模值；
[0047]
步骤6-2：对c
ifft
(m),1≤m≤l
sym
取出最大值p
max
以及最大值对应索引m
loca
，表示如下：
[0048][0049]
其中，为使得x最大的最大值p
max
、索引值m；
[0050]
步骤6-3：将最大值p
max
与ifft门限比较，当最大值p
max
超过ifft门限时，标记当前索引值m
p
，并进行同步头位置计算，其位置索引为(n-1)
·dstep
+m
p
，否则循环进行步骤2-6。
[0051]
进一步地，作为优选技术方案，所述步骤7的具体过程为：根据当前相关计算的次数n以及ifft相关最大值过门限的索引，计算得到同步头的具体位置，即同步头位置为(n-1)
·dstep
+m
p
，完成同步捕获。
[0052]
本发明相对于现有技术，具有如下有益效果：
[0053]
(1)本发明通过将同步头设计为三段相同的m序列分别经过gmsk调制后级联组成，且级联后的三段导频符号完全相同，采用前后段信号共轭相关计算，能够以步进小于等于一段导频符号的长度进行滑动相关计算，实现了非逐点滑动检测的大步进前后滑动检测，使接收终端设备在低运算量下进行正常捕获帧头位置，提高gmsk 突发无线通信系统捕获的效率，能够满足在低信噪比和大频偏环境下的无线突发通信系统的快速捕获与同步。
[0054]
(2)本发明通过采用前后段数据信号频偏估计，得到频偏估计结果，再将频偏估计值补偿到第二段信号中，进行频偏纠正，能够提高频偏估计结果的准确性，使得纠正后的信号能够进行后续的信号处理，从而解决gmsk突发通信系统中存在较大频偏难估计、估计精度较低的问题。
附图说明
[0055]
图1为本发明的流程图；
[0056]
图2为本发明的gmsk突发帧的同步头结构示意图；
[0057]
图3为本发明的共轭相关计算框图。
具体实施方式
[0058]
以下实施例用于说明本发明，但不用来限制本发明的范围。
[0059]
实施例
[0060]
如图1-图3所示，本实施例所述的一种面向gmsk突发通信系统的快速捕获与同步方法，具体如下：
[0061]
首先，如图2所示，本实施例将无线突发通信系统同步头设计为三段训练序列，每段导频序列长度为l
sym
个符号，每段导频由n
bit
比特的m序列经过高斯滤波最小频移键控(gmsk，gaussian filtered minimum shift keying)调制产生，然后级联为一帧，同步头内的三段导频符号完全相同，在进行无线突发通信同步捕获时，接收端存储一段导频符号，该段导频符号表示为：s
loc
(m)，1≤m≤l
sym
。具有三段相同导频的同步头，能够实现大步进的滑动共轭相关，为快速实现信号的捕获提供了条件。
[0062]
步骤1：利用地面终端设备接收来自发射端发射后的突发信号，并进行低通滤波以及下采样，假设终端实际的接收信号为r(m),1≤m≤lr，其中lr为接收到的信号的总长度。
[0063]
步骤2：如图3所示为两段数据共轭相关计算的框图，共轭相关计算每次滑动步进的步长为d
step
个采样点的数据(1≤d
step
≤l
sym
)，且在第n次共轭相关计算时，已滑动距首个样点(n-1)
·dstep
长度的距离，将经过低通滤波以及下采样后的突发信号从 m＝(n-1)
·dstep
+1开始向后取l
sym
个样点的信号，该段信号表示为 r
first,n
(m),1≤m≤l
sym
；从m＝(n-1
·
)d
espt
+l sym
+1开始向后取l
sym
个样点的信号，该段信号表示为r
second,n
(m),1≤m≤l
sym
，将两段信号做共轭相关计算，计算方法为：
[0064][0065]
其中，(
·
)
*
表示求取共轭值。
[0066]
步骤3：将步骤2得到的共轭值与相关门限值进行比较，若共轭值超过相关门限值，则标记当前相关次数n，并进行步骤4，否则，继续进行步骤2的相关计算，并再次与相关门限值进行比较。
[0067]
步骤4：利用当前超过相关门限值时的信号r
first,n
(m),1≤m≤l、 r
second,n
(m),1≤m≤l
sym
进行频偏估计，频偏估计具体表示为：
[0068][0069][0070]
其中，(
·
)
*
表示求取共轭值，fs为降采样后的信号频率；arg{x}为求得x的相位角，为频偏估计的结果；
[0071]
将频偏估计的结果补偿到第二段信号r
second,n
(m),1≤m≤l
sym
，得到如下信号：
[0072][0073]
步骤5：将频偏补偿后的信号与本地保存的导频信号 s
loc
(m)，1≤m≤l
sym
分别做l
sym
点的fft计算，表示如下：
[0074][0075]
s(m)＝fft(s
loc
(m))，1≤m≤l
sym
；
[0076]
将步骤5-1得到的fft计算结果进行共轭相乘，表示如下：
[0077]cfft
(m)＝r(m)
·s*
(m),1≤m≤l
sym
[0078]
其中，fft(
·
)表示为fft计算。
[0079]
步骤6：步骤5得到的共轭相乘结果c
fft
进行l
sym
点的ifft计算，并取模值平方，表示如下：
[0080]cifft
(m)＝|ifft(c
fft
(m))2,1≤m≤l
sym
[0081]
其中，ifft(
·
)表示为ifft计算，|
·
|为取模值；
[0082]
对c
ifft
(m),1≤m≤l
sym
取出最大值p
max
以及最大值对应索引m
loca
，表示如下：
[0083][0084]
其中，为使得x最大的最大值p
max
、索引值m；
[0085]
将最大值p
max
与ifft门限比较，当最大值p
max
超过ifft门限时，标记当前索引值m
p
，并进行同步头位置计算，其位置索引为(n-1)
·dstep
+m
p
，否则循环进行步骤 2-6。
[0086]
步骤7：利用步骤3中标记的相关次数与步骤6得到的索引值，计算同步头的起始位置，即同步头位置为(n-1)
·dstep
+m
p
，完成同步捕获。
[0087]
本发明通过对接收信号进行滤波和下采样、前后段数据共轭相关计算、相关门限比较、频偏估计与补偿、快速傅立叶变换(fft，fast fourier transform)与共轭计算、快速傅里叶逆变换(ifft，inverse fast fourier transform)计算与ifft门限比较，最后计算同步头的位置，能够满足在低信噪比和大频偏环境下的无线突发通信系统的快速捕获与同步，使接收终端设备在低运算量下进行正常捕获帧头位置，提高 gmsk突发无线通信系统捕获的效率。
[0088]
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。