一种卫星通信载波同步方法、装置及通信设备与流程

1.本技术涉及卫星通信技术领域，具体而言，涉及一种卫星通信载波同步方法、装置，通信设备及计算机可读存储介质。


背景技术：

2.随着科技的飞速发展，越来越多的移动设备可接入到卫星网络中。其中，移动设备可包括中低速移动的地面终端设备(如高铁、汽车等)和高速移动的空中终端设备(如无人机、飞机等)。移动设备通常以突发模式随机接入卫星网络，以实现节省功耗和导频资源的目的。然而移动设备在接入卫星时面临的一个严重问题是多普勒效应引起的载波频率偏移。特别是考虑到高速移动设备以突发模式接入卫星系统时，接收到的观测数据比较短，由于链路预算导致的信噪比低，并且由于终端高速随机运动导致多普勒频移无法由固定的多普勒特性曲线补偿，这些原因使得该场景下的载波同步变得非常困难。
3.传统的载波同步方法中，锁相环需要很长的接收数据长度并且对信噪比要求很高；非数据辅助方法需要在高信噪比下才能正常工作；数据辅助方法依赖于导频，造成频谱资源的浪费，并且无法纠正大范围的频偏；编码辅助虽然可以在低信噪比工作并且不依赖导频，但是当频偏变大时，需要更长的接收数据，因此也无法适用于卫星与高移动设备的突发通信。
4.可见，目前的载波同步方法不能在低信噪比和短接收数据的场景下纠正大范围的多普勒频移，因此很难适用于卫星与高移动设备的突发通信。


技术实现要素：

5.本技术实施例的目的在于提供一种卫星通信载波同步方法、装置，通信设备及计算机可读存储介质，以解决在高移动设备接入的场景下，多普勒效应对卫星通信的影响。
6.本发明是这样实现的：
7.第一方面，本技术实施例提供一种卫星通信载波同步方法，包括：步骤一：获取基于移动设备的突变信号生成的频偏粗估计值；步骤二：将所述频偏粗估计值作为高斯过程模型的初始观测值，构建训练数据集，并基于所述高斯过程模型构建目标函数；其中，所述训练数据集包括观测值及与每个所述观测值对应的目标函数观测向量；所述目标函数观测向量通过sccpm系统解码输出的均方软输出得到；所述观测值包括所述初始观测值；所述目标函数为通过所述高斯过程模型构造所述均方软输出得到；步骤三：基于所述训练数据集，计算所述目标函数的均值和所述目标函数的方差；步骤四：基于复合积分法则和高斯累积分布函数，通过所述目标函数的均值和所述目标函数的方差从多普勒频偏的搜索范围中估计所述目标函数的极大值；步骤五：基于高斯估计函数，通过所述目标函数的均值、所述目标函数的方差和所述目标函数的极大值，得到新的频偏估计值；步骤六：将所述新的频偏估计值作为一个观测值添加到所述训练数据集中，重复步骤一至步骤五，直至得到第t个频偏估计值；其中，t为预设的大于零的自然数；所述第t个频偏估计值为频偏精估计值；步骤七：
基于所述频偏精估计值补偿频偏，得到载波同步信号。
8.本技术实施例所提供的卫星通信载波同步方法，通过高斯过程模型所构建的目标函数来估计基于串行级联连续相位调制(sccpm)的突发信号的大范围多普勒频偏。借助复合积分法则和高斯累积分布函数，得到了高斯估计模型中目标函数的极大值。然后通过高斯估计函数，得到新的频偏估计值，该函数能够实现更加准确地估计频偏并避免手动调节参数。相比较于传统算法，本技术实施例提供的卫星通信载波同步方法可以在观测数据较少的情况下估计大范围的频偏并且具有很低的最优参数搜索复杂度。
9.结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述获取基于移动设备的突变信号生成的频偏粗估计值，包括：获取卫星接收到的所述移动设备发送的突发信号；基于期望最大化载波同步算法对所述突发信号进行载波粗同步，得到所述频偏粗估计值。
10.结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述高斯估计函数的表达式为：
[0011][0012]
其中，δf
t
表示所述新的频偏估计值；表示所述目标函数的极大值，表示所述目标函数的均值；表示所述目标函数的方差。
[0013]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，在所述得到新的频偏估计值之后，所述方法还包括：获取所述频偏粗估计值的误差长度；基于所述新的频偏估计值和所述频偏粗估计值的差值与所述误差长度的大小关系，确定所述目标函数的均值和所述目标函数的方差中核函数的长度尺度参数的取值。
[0014]
在本技术实施例中，通过新的频偏估计值和频偏粗估计值的差值与误差长度的大小关系，来实现对长度尺度参数的自适应调节，进而使得后续输出的频偏估计值更加准确。
[0015]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述基于所述新的频偏估计值和所述频偏粗估计值的差值与所述误差长度的大小关系，确定所述目标函数的均值和所述目标函数的方差中核函数的长度尺度参数的取值，包括：当所述新的频偏估计值和所述频偏粗估计值的差值大于所述误差长度时，将所述长度尺度参数减去第一预设增量；当所述新的频偏估计值和所述频偏粗估计值的差值小于所述误差长度时，将所述长度尺度参数加上所述第一预设增量。
[0016]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述目标函数为凹函数，所述方法还包括：基于所述频偏粗估计值对应的目标函数观测向量与预设值对应的目标函数观测向量确定所述多普勒频偏的搜索范围；其中，所述预设值为所述频偏粗估计值加上第二增量。
[0017]
在本技术实施例中，通过凹函数的函数特性，基于频偏粗估计值对应的目标函数观测向量与预设值对应的目标函数观测向量确定所述多普勒频偏的搜索范围，进而缩小多普勒频偏的搜索范围，提高了通信设备的处理效率。
[0018]
结合上述第一方面提供的技术方案，在一些可能的实现方式中，所述基于所述频
偏粗估计值对应的目标函数观测向量与预设值对应的目标函数观测向量确定所述多普勒频偏的搜索范围，包括：当所述频偏粗估计值对应的目标函数观测向量大于所述预设值对应的目标函数观测向量时，所述多普勒频偏的搜索范围为频偏最小值至所述频偏粗估计值；当所述频偏粗估计值对应的目标函数观测向量小于所述预设值对应的目标函数观测向量时，所述多普勒频偏的搜索范围为所述频偏粗估计值至频偏最大值。
[0019]
第二方面，本技术实施例提供一种卫星通信载波同步装置，包括：获取模块，用于执行步骤一：获取基于移动设备的突变信号生成的频偏粗估计值；构建模块，用于执行步骤二：将所述频偏粗估计值作为高斯过程模型的初始观测值，构建训练数据集，并基于所述高斯过程模型构建目标函数；其中，所述训练数据集包括观测值及与每个所述观测值对应的目标函数观测向量；所述目标函数观测向量通过sccpm系统解码输出的均方软输出得到；所述观测值包括所述初始观测值；所述目标函数为通过所述高斯过程模型构造所述均方软输出得到；计算模块，用于执行步骤三：基于所述训练数据集，计算所述目标函数的均值和所述目标函数的方差；估计模块，用于执行步骤四：基于复合积分法则和高斯累积分布函数，通过所述目标函数的均值和所述目标函数的方差从多普勒频偏的搜索范围中估计所述目标函数的极大值；生成模块，用于执行步骤五：基于高斯估计函数，通过所述目标函数的均值、所述目标函数的方差和所述目标函数的极大值，得到新的频偏估计值；处理模块，用于执行步骤六：将所述新的频偏估计值作为一个观测值添加到所述训练数据集中，重复步骤一至步骤五，直至得到第t个频偏估计值；其中，t为预设的大于零的自然数；所述第t个频偏估计值为频偏精估计值；补偿模块，用于执行步骤七：基于所述频偏精估计值补偿频偏，得到载波同步信号。
[0020]
第三方面，本技术实施例提供一种通信设备，包括：处理器和存储器，所述处理器和所述存储器连接；所述存储器用于存储程序；所述处理器用于调用存储在所述存储器中的程序，执行如上述第一方面实施例和/或结合上述第一方面实施例的一些可能的实现方式提供的方法。
[0021]
第四方面，本技术实施例提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序在被处理器运行时执行如上述第一方面实施例和/或结合上述第一方面实施例的一些可能的实现方式提供的方法。
附图说明
[0022]
为了更清楚地说明本技术实施例的技术方案，下面将对本技术实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0023]
图1为本技术实施例提供的一种通信设备的结构示意图。
[0024]
图2为本技术实施例提供的一种卫星通信载波同步方法的步骤流程图。
[0025]
图3为本技术实施例提供的一种卫星通信载波同步装置的模块框图。
[0026]
图标：100
‑
通信设备；110
‑
处理器；120
‑
存储器；200
‑
卫星通信载波同步装置；210
‑
获取模块；220
‑
构建模块；230
‑
计算模块；240
‑
估计模块；250
‑
生成模块；260
‑
处理模块；270
‑
补偿模块。
具体实施方式
[0027]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行描述。
[0028]
请参阅图1，本技术实施例提供的一种应用卫星通信载波同步方法及装置的通信设备100的结构示意图。
[0029]
本技术实施例中，通信设备100可以是卫星。通信设备100也可以是地面工作站上的任意电子设备(如计算机、服务器)，电子设备与卫星之间通信连接。
[0030]
当通信设备100为卫星时，卫星接收移动设备所发送的突变信号，并基于突变信号执行本技术实施例所提供的卫星通信载波同步方法。而当通信设备100为地面工作站上的任意电子设备时，当卫星接收到移动设备所发送的突变信号后，将突变信号发送给地面工作站上的电子设备，由电子设备执行本技术实施例所提供的卫星通信载波同步方法。
[0031]
在一些实施例中，卫星在接收到移动设备所发送的突变信号后，可以先基于突变信号生成的频偏粗估计值，然后将生成的频偏粗估计值发送给地面工作站上的电子设备，然后再由电子设备执行本技术实施例所提供的卫星通信载波同步方法。
[0032]
上述的移动设备可包括中低速移动的地面终端设备(如高铁、汽车等)和高速移动的空中终端设备(如无人机、飞机等)。
[0033]
在结构上，通信设备100可以包括处理器110和存储器120。
[0034]
处理器110与存储器120直接或间接地电性连接，以实现数据的传输或交互，例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。卫星通信载波同步装置包括至少一个可以软件或固件(firmware)的形式存储在存储器120中或固化在通信设备100的操作系统(operating system，os)中的软件模块。处理器110用于执行存储器120中存储的可执行模块，例如，卫星通信载波同步装置所包括的软件功能模块及计算机程序等，以实现卫星通信载波同步方法。处理器110可以在接收到执行指令后，执行计算机程序。
[0035]
其中，处理器110可以是一种集成电路芯片，具有信号处理能力。处理器110也可以是通用处理器，例如，可以是中央处理器(central processing unit，cpu)、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、分立门或晶体管逻辑器件、分立硬件组件，可以实现或者执行本技术实施例中的公开的各方法、步骤及逻辑框图。此外，通用处理器可以是微处理器或者任何常规处理器等。
[0036]
存储器120可以是，但不限于，随机存取存储器(random access memory，ram)、只读存储器(read only memory，rom)、可编程只读存储器(programmable read
‑
only memory，prom)、可擦可编程序只读存储器(erasable programmable read
‑
only memory，eprom)，以及电可擦编程只读存储器(electric erasable programmable read
‑
only memory，eeprom)。存储器120用于存储程序，处理器110在接收到执行指令后，执行该程序。
[0037]
需要说明的是，图1所示的结构仅为示意，本技术实施例提供的通信设备100还可以具有比图1更少或更多的组件，或是具有与图1所示不同的配置。此外，图1所示的各组件可以通过软件、硬件或其组合实现。
[0038]
请参阅图2，图2为本技术实施例提供的卫星通信载波同步方法的步骤流程图，该方法应用于图1所示的通信设备100中。需要说明的是，本技术实施例提供的卫星通信载波同步方法方法不以图2及以下所示的顺序为限制，该方法包括：步骤一
‑
步骤七。
[0039]
步骤一(s101)：获取基于移动设备的突变信号生成的频偏粗估计值。
[0040]
于本技术实施例中，移动设备的突变信号指的是移动设备在高速移动过程中以突发模式所接入到卫星的信号。其中，移动设备的高速的具体数值可以根据情况而定，本技术不作限定。
[0041]
可选地，步骤一可具体包括：获取卫星接收到的移动设备发送的突发信号；基于期望最大化载波同步算法对突发信号进行载波粗同步，得到频偏粗估计值。
[0042]
上述的突发信号由移动设备发出，卫星在接收到突发信号后采用期望最大化载波同步算法进行粗同步，以得到频偏粗估计。当然，也可以是卫星接收到突发信号将其发送至地面工作站的电子设备中，以使电子设备采用期望最大化载波同步算法进行粗同步，以得到频偏粗估计。(为了便于说明，后文中均以卫星作为执行主体对本技术实施例的卫星通信载波同步方法的进行说明)。
[0043]
具体的，移动设备产生突变信息序列a＝[a0,a1,...,a
k
‑1]；然后经过卷积编码器被编码为序列c＝[c0,c1,...,c
l
‑1]；信道编码的编码率为r＝k/l；k为信息比特长度；l为编码后的序列长度。
[0044]
然后编码信号通过连续相位调制器以产生基带信号s＝[s0,s1,...,s
l
‑1]，之后该信号经过上变频发射出去。
[0045]
当移动设备在高速移动时，卫星接收到的信号将遭受到较大的多普勒频率偏移。考虑到噪声和多普勒频移的存在，卫星接收的信号可以表示为：
[0046]
r(t)＝s(t)exp(j2π(f
c
f
d
(t))t) ω(t)
ꢀꢀꢀꢀꢀꢀ
(1)
[0047][0048]
在公式(1)和公式(2)中，是cpm(continue phase modulation，连续相位调制)倾斜相位，f
c
是载波频率，f
d
(t)＝f
c
ν
a
(t)/c，ν
a
(t)是移动设备的速度，而ω(t)是零均值和方差为n0/2的复高斯白噪声；e表示发射信号能量，t为发射信号持续时间，c为编码比特，j为复数虚部符号。f
d
(t)在t0＝0附近的泰勒级数展开由下式给出：
[0049]
f
d
(t)＝f
d
(0) f
′
d
(0)t o(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0050]
公式(3)中，f
d
(0)是多普勒频移，f
′
d
(0)是多普勒速率，o(t)表示t的高阶无穷小。将(3)代入(1)得到：
[0051]
r(t)＝s(t)exp(j2π(f
d
(0)t f
′
d
(0)t2)) ω(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0052]
考虑到卫星突发通信中帧的时间通常小于10毫秒，f
′
d
(0)t2对每个的影响符号可以忽略不计。
[0053]
然后再将公式(4)经过下变频和离散，得到的接收信号可以表示为：
[0054]
r(n)＝s(n)exp(j2πnf
d
t
s
) ω(n),n＝1,2,
…
,kr
s
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0055]
其中r
s
是每个符号k的采样点数，t
s
是符号周期，f
d
t
s
被称为归一化频率偏移(nfo)。
[0056]
然后卫星再根据后验概率思想，进行粗略的解调与译码，为后续载波同步产生必要的信息。
[0057]
需要说明的是，由于前述的所有过程均为本领域所熟知，因此，本技术不作详细说
明。
[0058]
在完成粗略的解调与译码后，基于期望最大化载波同步算法对突发信号进行载波粗同步。期望最大化载波同步算法的公式包括：
[0059][0060][0061][0062][0063]
在公式(6)～公式(9)中，为接收信号的先验对数概率，为sccpm(串行级联连续相位调制)系统发送符号s的后验概率，可以通过最大后验概率思想进行计算，r
s
为采样倍数，s
n
＝(v
n
,u
n
)表示连续相位调制发送的符号，v
n
表示n时刻的状态，u
n
表示传输的信息序列，r
l
表示卫星接收的离散信号，(v的表示相同)，f
d
表示频偏，t
s
表示系统信息发送时间间隔，r表示卫星接收的离散信号向量，e
s
表示符号能量，n0表示噪声能量。通过公式(6)即可计算得到频偏粗估计值。
[0064]
步骤二(s102)：将频偏粗估计值作为高斯过程模型的初始观测值，构建训练数据集，并基于高斯过程模型构建目标函数。
[0065]
其中，训练数据集包括观测值及与每个观测值对应的目标函数观测向量；目标函数观测向量通过sccpm系统解码输出的均方软输出得到；观测值包括初始观测值；目标函数为通过高斯过程模型构造均方软输出得到。
[0066]
具体的，根据步骤一中的载波粗同步算法，得到频偏的粗略估计值δf
coarse
。选择观测频偏向量f0＝[δf
i
；δf
coarse
；δf
j
]，其中，δf
i
和δf
j
为δf
coarse
相邻的频偏数值。当然，观测频偏向量f0可以包括更多个数的频偏数值，比如五个、七个等等，本技术不作限定。
[0067]
然后基于观测频偏向量f0＝[δf
i
；δf
coarse
；δf
j
]计算对应的目标函数观测向量ψ＝[ψ(δf
i
)；ψ(δf
coarse
)；ψ(δf
j
)]，进而得到训练数据集初始化长度尺度参数l1。
[0068]
需要说明的是，当频率偏移为零时，从sccpm解码输出的对数似然比获得的均方软输出(maso)具有全局最大值，因此它可以作为目标函数来反映sccpm系统频率偏移校正的程度。本技术实施例中，通过均方软输出来得到目标函数观测向量。该均方软输出表示为：
[0069][0070]
公式(10)中，λ(n)是sccpm系统解码输出的对数似然比，其具有对数正态分布的结构。
[0071]
由于对数似然比，其具有对数正态分布的结构，因此其可以通过高斯分布来近似，也即，本技术实施例中采用高斯过程(gp)模型来构建目标函数。
[0072]
真实随机变量的高斯过程(gp)模型是一组随机变量，服从联合高斯分布。在实际系统中，可以合理地认为目标函数由与频偏相关的函数和噪声构成。因此，其表达式为：
[0073][0074]
ψ(δf)＝gp(m(f),κ(f,f'))
ꢀꢀꢀꢀꢀꢀꢀ
(12)
[0075]
在上述公式(11)～公式(12)中，是预测目标函数ψ(δf)的平均值，ε是高斯噪声服从m(f)为gp模型的均值，κ(f,f')为gp模型的协方差。
[0076]
步骤三(s103)：基于训练数据集，计算目标函数的均值和目标函数的方差。
[0077]
在得到训练数据集后，即可基于训练数据集计算上述构建的目标函数的真实均值和真实方差。
[0078]
当给出不同的多普勒频率偏移点f＝[δf1,δf2,...,δf
q
]和目标函数的观测值ψ＝[ψ(δf1),ψ(δf2),
…
,ψ(δf
q
)]作为输出，目标函数基于{f,ψ}的后验分布服从高斯分布，其平均值和方差的表达式如下：
[0079][0080][0081]
在公式(13)～公式(14)，δf表示待测试的多普勒频偏点，i表示单位矩阵，κ(f,δf)＝κ(δf,f)
t
，κ(δf,f)＝[κ(δf,δf1),
…
,κ(δf,δf
q
)],
[0082]
需要说明的是，κ(
·
)称为核函数，通常选择平方指数协方差函数：
[0083][0084]
在公式(15)中，η表示幅度，l表示长度尺度参数。
[0085]
步骤四(s104)：基于复合积分法则和高斯累积分布函数，通过目标函数的均值和目标函数的方差从多普勒频偏的搜索范围中估计目标函数的极大值。
[0086]
在求得目标函数的平均值和方差后，再利用复合积分法则和高斯累积分布函数从多普勒频偏的搜索范围中估计目标函数的极大值
[0087]
具体的，定义目标函数的极大值的表达式为：
[0088][0089]
公式(16)中，表示的均值，是观测频率偏移点和目标函数观测向量的集合，是目标函数的最大值，可以看作是一个高斯变量，y可以看作高斯变量。本技术实施例中，y大于零，目标函数的极大值的计算方法如下：
[0090][0091]
其中，为多普勒频偏的搜索范围。
[0092]
设置m0＝max
δf∈f
(ψ(δf))，则：
[0093][0094]
当y∈[0,m0]时，的概率为1，则有
[0095][0096][0097]
在公式(19)～公式(20)中，φ(x)为高斯累积函数。公式(19)～公式(20)提供了一种计算目标函数最大值的自适应方法。然而，在实际模拟中，发现由上式计算的最大值小于实际值，导致置信界的提升不足。这是因为上式假设观测过程中的噪声可以忽略，而目标函数ψ(δf)由预测的平均值和噪声ε组成。因此，需要得到的上界，以便灵活地设置其值，为自适应高斯置信界改进算法提供依据。上式可由复合积分规则计算，复合积分规则为：
[0098][0099]
在公式(21)中，x
k
＝a kh'，n为积分点数，h、h
′
为积分的步长，k为变量，a和b为积分的上下限。则将公式(21)代入公式(19)得：
[0100][0101]
其中h
′
为积分步长，假设当k＝n0时，
[0102]
其意味着因此，m
up
＝m0 n0h'为的上界。为了更好地估计目标函数的最大值，需要更多关于目标函数的先验信息。因此，此处模拟了均方软输出在不同e
b
/n0下的全局极大值，发现均方软输出的值随着e
b
/n0的增加而增加。因此，可以根据e
b
/n0灵活地设置以更好地估计频偏。其中，e
b
/n0也称为信噪比，表示比特信号能量与噪声能量的比值，e
b
为信号发射功率，n0为噪声功率。
[0103]
综上可得，目标函数的极大值为：m
up
为的上界，m
up
＝m0 n0h'，m0＝max
δf∈f
(ψ(δf))，h
′
为积分步长，n0满足，当k＝n0时，其中，0＜k≤1。
[0104]
步骤五(s105)：基于高斯估计函数，通过目标函数的均值、目标函数的方差和目标函数的极大值，得到新的频偏估计值。
[0105]
由高斯估计函数(gp
‑
est函数)得到新的频偏估计值δf
t
；高斯估计函数的表达式为：
[0106][0107]
公式(23)中，δf
t
表示新的频偏估计值；表示目标函数的极大值，表示目标函数的均值；表示目标函数的方差。
[0108]
步骤六(s106)：将新的频偏估计值作为一个观测值添加到训练数据集中，重复步骤一至步骤五，直至得到第t个频偏估计值；其中，t为预设的大于零的自然数；第t个频偏估计值为频偏精估计值。
[0109]
然后，将新的频偏估计值δf
t
(及与δf
t
对应的目标函数观测向量)添加到训练数据集中，重复执行上述步骤一至步骤五，直至输出第t个频偏估计值。然后，将第t个频偏估计值确定为频偏精估计值。其中，t可以根据模拟情况设定的具体数值，如10、20，本技术不作限定。
[0110]
可选地，在每次得到新的频偏估计值之后，还会对目标函数的均值和方差中的核函数的长度尺度参数进行更新，具体的，更新方式包括：获取频偏粗估计值的误差长度；基于新的频偏估计值和频偏粗估计值的差值与误差长度的大小关系，确定目标函数的均值和目标函数的方差中核函数的长度尺度参数的取值。
[0111]
确定方式具体包括：当新的频偏估计值和频偏粗估计值的差值大于误差长度时，将长度尺度参数减去第一预设增量；当新的频偏估计值和频偏粗估计值的差值小于误差长度时，将长度尺度参数加上第一预设增量。
[0112]
高斯过程模型的核函数长度标度l的设置是提高置信界，改进函数估计性能的关键。从具有平稳核的高斯过程模型中抽取的典型样本在整个输入空间中的变化率趋于相似，减少l将增加预测函数的基函数类，从而增大置信区间，以找到凸函数的最大点。设定频偏粗估计值的误差长度δυ；其中，误差长度δυ为误差总长度的二分之一。
[0113]
当|δf
t
‑
δf
coarse
|＞δυ时，表示高斯过程模型估计的频偏值离移动设备频偏点较远，因此需要减小l以增加移动设备频偏点周围的探测。相反，l则需要增加。更新长度尺度参数l
t 1
的表达式为：
[0114]
[0115]
公式(24)中，l
t
当前的长度尺度参数；l
t 1
表示更新的长度尺度参数；δ2表示第一预设增量；else表示其他情况，该情况下，长度尺度参数数值不变，δf
t
表示新的频偏估计值。
[0116]
在本技术实施例中，通过新的频偏估计值和频偏粗估计值的差值与误差长度的大小关系，来实现对长度尺度参数的自适应调节，进而使得后续输出的频偏估计值更加准确。
[0117]
步骤七(s107)：基于频偏精估计值补偿频偏，得到载波同步信号。
[0118]
最后，在得到频偏精估计值后，即可根据频偏精估计值对频偏进行补偿，进而得到载波同步信息。
[0119]
综上，本技术实施例所提供的卫星通信载波同步方法，通过高斯过程模型所构建的目标函数来估计基于串行级联连续相位调制(sccpm)的突发信号的大范围多普勒频偏。借助复合积分法则和高斯累积分布函数，得到了高斯估计模型中目标函数的最大估计值的上界。然后通过高斯估计函数，得到新的频偏估计值，该函数能够实现更加准确地估计频偏并避免手动调节参数。相比较于传统算法，本技术实施例提供的卫星通信载波同步方法可以在观测数据较少的情况下估计大范围的频偏并且具有很低的最优参数搜索复杂度。
[0120]
此外，由于目标函数为凹函数，因此，本技术实施例还提供一种缩小多普勒频偏的搜索范围的方式，以提高通信设备的处理效率。该方法还包括：基于频偏粗估计值对应的目标函数观测向量与预设值对应的目标函数观测向量确定多普勒频偏的搜索范围；其中，预设值为频偏粗估计值加上第二增量。
[0121]
具体的，当频偏粗估计值对应的目标函数观测向量大于预设值对应的目标函数观测向量时，多普勒频偏的搜索范围为频偏最小值至频偏粗估计值；当频偏粗估计值对应的目标函数观测向量小于预设值对应的目标函数观测向量时，多普勒频偏的搜索范围为频偏粗估计值至频偏最大值。
[0122]
示例性的，目标函数的全局最大值在δf
d
，因此载波粗同步算法估计的粗频偏点δf
coarse
可能在最大点的左侧或右侧。
[0123]
如果频率偏移点在目标函数的递增范围ψ(δf
coarse
δ)＞ψ(δf
coarse
)，则δf
coarse
＜δf
d
。相反，若ψ(δf
coarse
δ)＜ψ(δf
coarse
)，δf
coarse
＞δf
d
，δ表示增量。
[0124]
因此，多普勒频偏的搜索范围缩小方法如下：当ψ(δf
coarse
δ1)＞ψ(δf
coarse
)，则否则其中，为多普勒频偏的搜索范围，δ1为第二增量，为频偏最大值，为频偏最小值。
[0125]
请参阅图3，基于同一发明构思，本技术实施例还提供一种卫星通信载波同步装置200，该装置包括：
[0126]
获取模块210，用于执行步骤一：获取基于移动设备的突变信号生成的频偏粗估计值。
[0127]
构建模块220，用于执行步骤二：将所述频偏粗估计值作为高斯过程模型的初始观测值，构建训练数据集，并基于所述高斯过程模型构建目标函数；其中，所述训练数据集包括观测值及与每个所述观测值对应的目标函数观测向量；所述目标函数观测向量通过sccpm系统解码输出的均方软输出得到；所述观测值包括所述初始观测值；所述目标函数为通过所述高斯过程模型构造所述均方软输出得到。
[0128]
计算模块230，用于执行步骤三：基于所述训练数据集，计算所述目标函数的均值和所述目标函数的方差。
[0129]
估计模块240，用于执行步骤四：基于复合积分法则和高斯累积分布函数，通过所述目标函数的均值和所述目标函数的方差从多普勒频偏的搜索范围中估计所述目标函数的极大值。
[0130]
生成模块250，用于执行步骤五：基于高斯估计函数，通过所述目标函数的均值、所述目标函数的方差和所述目标函数的极大值，得到新的频偏估计值。
[0131]
处理模块260，用于执行步骤六：将所述新的频偏估计值作为一个观测值添加到所述训练数据集中，重复步骤一至步骤五，直至得到第t个频偏估计值；其中，t为预设的大于零的自然数；所述第t个频偏估计值为频偏精估计值。
[0132]
补偿模块270，用于执行步骤七：基于所述频偏精估计值补偿频偏，得到载波同步信号。
[0133]
可选地，获取模块210具体用于获取卫星接收到的所述移动设备发送的突发信号；基于期望最大化载波同步算法对所述突发信号进行载波粗同步，得到所述频偏粗估计值。
[0134]
可选地，该装置还包括长度尺度参数确定模块。长度尺度参数确定模块用于获取所述频偏粗估计值的误差长度；基于所述新的频偏估计值和所述频偏粗估计值的差值与所述误差长度的大小关系，确定所述目标函数的均值和所述目标函数的方差中核函数的长度尺度参数的取值。
[0135]
可选地，长度尺度参数确定模块具体用于当所述新的频偏估计值和所述频偏粗估计值的差值大于所述误差长度时，将所述长度尺度参数减去第一预设增量；当所述新的频偏估计值和所述频偏粗估计值的差值小于所述误差长度时，将所述长度尺度参数加上所述第一预设增量。
[0136]
可选地，该装置还包括搜索范围确定模块。搜索范围确定模块用于基于所述频偏粗估计值对应的目标函数观测向量与预设值对应的目标函数观测向量确定所述多普勒频偏的搜索范围；其中，所述预设值为所述频偏粗估计值加上第二增量。
[0137]
可选地，搜索范围确定模块具体用于当所述频偏粗估计值对应的目标函数观测向量大于所述预设值对应的目标函数观测向量时，所述多普勒频偏的搜索范围为频偏最小值至所述频偏粗估计值；当所述频偏粗估计值对应的目标函数观测向量小于所述预设值对应的目标函数观测向量时，所述多普勒频偏的搜索范围为所述频偏粗估计值至频偏最大值。
[0138]
需要说明的是，由于所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0139]
基于同一发明构思，本技术实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序在被运行时执行上述实施例中提供的方法。
[0140]
该存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，dvd)、或者半导体介质(例如固态硬盘solid state disk(ssd))等。
[0141]
在本技术所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方
式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0142]
另外，作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0143]
再者，在本技术各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
[0144]
在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
[0145]
以上所述仅为本技术的实施例而已，并不用于限制本技术的保护范围，对于本领域的技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。