CN108989256B - 一种fsk/gfsk解调方法及装置 - Google Patents

Abstract

本发明提供了一种FSK/GFSK解调方法及装置，涉及通信技术，在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，确定接收到的FSK/GFSK基带信号

的接收相位向量

后，还要确定

与θ_l(i)之间的平均相位差β_l，并在去除平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，找出使相位匹配程度Q_l最大的l值，进而确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。由于在进行相位匹配时，去除了平均相位差的影响，增大了相位匹配的准确性，提高了相位域解调技术的性能。

Description

一种FSK/GFSK解调方法及装置

技术领域

本发明属于数据处理领域，具体涉及一种FSK/GFSK解调方法及装置。

背景技术

现有基于匹配滤波器组的FSK/GFSK解调方法主要是在复数域(I\Q平面)完成解调操作，需要进行大量的复数乘法、复数加减法运算，系统复杂度高。

由于FSK/GFSK信号具有恒定的包络，所传输的数字信息全部包含在信号的相位中，在相位域完成解调操作，能够将大量的复数运算转化为实数加减法运算，从而可以有效的降低系统复杂度。

然而，本发明的发明人发现，现有的相位域解调技术或者没有使用相位匹配滤波器组(即本发明中的相位匹配向量组)来进行解调操作，或者在使用相位匹配向量组时存在技术缺陷，使得现有的相位域解调技术的性能远低于复数域解调技术的性能。

发明内容

针对上面提到的相位域解调技术的性能较低的问题，提出了一种FSK/GFSK解调方法及装置，以提高相位域解调技术的性能。

依据本发明的第一方面，提供了一种FSK/GFSK解调方法，包括：

在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号；

确定接收到的FSK/GFSK基带信号

的接收相位向量

确定

与θ_l(i)之间的平均相位差β_l；

去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值；

确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。

进一步，所述去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值，具体包括：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

更进一步，所述f[x]具体为：

f[x]＝cos(x)；或者

或者

C为设定常数；或者

C为设定常数。

更进一步，所述确定

与θ_l(i)相位匹配程度时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，只使用其解调结果对应的θ_l(i)与

进行相位匹配。

进一步，所述在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，具体包括：

对于2FSK/2GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^{2M +1}个数字信息向量V_l(i),l＝1,2,…,2^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化；

对于4FSK/4GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。

进一步，所述确定

与θ_l(i)之间的平均相位差β_l，具体包括：

确定

与θ_l(i)之间的平均相位差为

依据本发明的第二方面，提供一种FSK/GFSK解调装置，包括：

相位匹配向量组确定单元，用于在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号；

接收相位向量确定单元，用于确定接收到的FSK/GFSK基带信号

的接收相位向量

平均相位差确定单元，用于确定

与θ_l(i)之间的平均相位差β_l；

匹配单元，用于去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定相位匹配程度Q_l最大的l值；

解调结果确定单元，用于确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。

进一步，所述匹配单元具体用于：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

更进一步，所述f[x]具体为：

f[x]＝cos(x)；或者

或者

C为设定常数；或者

C为设定常数。

更进一步，所述匹配单元计算

与θ_l(i)相位匹配程度Q_l时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，只使用其解调结果对应的θ_l(i)与

进行相位匹配。

进一步，相位匹配向量组确定单元具体用于：

对于2FSK/2GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^{2M +1}个数字信息向量V_l(i),l＝1,2,…,2^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化；

对于4FSK/4GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。

进一步，所述平均相位差确定单元具体用于：

确定

与θ_l(i)之间的平均相位差为

本发明实施例提供一种FSK/GFSK解调方法及装置，在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，确定接收到的FSK/GFSK基带信号

的接收相位向量

后，还要确定

与θ_l(i)之间的平均相位差β_l，并在去除平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值，进而确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。由于在进行相位匹配时，去除了平均相位差的影响，增大了相位匹配的准确性，提高了相位域解调技术的性能。

应当理解，上述说明仅是本发明技术方案的概述，以便能够更清楚地了解本发明的技术手段，从而可依照说明书的内容予以实施。为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举说明本发明的具体实施方式。

附图说明

通过阅读下文的示例性实施例的详细描述，本领域普通技术人员将明白本文所述的有点和益处以及其他优点和益处。附图仅用于示出示例性实施例的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的标号表示相同的部件。在附图中：

图1为本发明实施例提供的FSK/GFSK解调方法流程图；

图2为本发明实施例提供的FSK/GFSK解调装置结构示意图。

具体实施方式

下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。

在本发明中，应理解，诸如“包括”或“具有”等术语旨在指示本说明书中所公开的特征、数字、步骤、行为、部件、部分或其组合的存在，并且不旨在排除一个或多个其他特征、数字、步骤、行为、部件、部分或其组合存在的可能性。

另外还需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。

如图1所示，本发明实施例提供的FSK/GFSK解调方法，包括：

步骤S101、在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号；

步骤S102、确定接收到的FSK/GFSK基带信号

的接收相位向量

步骤S103、确定

与θ_l(i)之间的平均相位差β_l；

步骤S104、去除平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值；

步骤S105、确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。

由于在进行相位匹配时，去除了平均相位差的影响，增大了相位匹配的准确性，提高了相位域解调技术的性能。

由于在相位域解调FSK/GFSK信号，将大量的复数乘法、加减法运算转化为实数加减法运算，大大降低了解调的系统复杂度，并且，使用相位匹配向量组进行解调，在解调前消除了平均相位差对相位匹配的不良影响，增大了相位匹配的准确性，提高了相位域解调技术的性能。

进一步，步骤S104，去除平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值，具体包括：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

通过

确定相位匹配程度，可以较好的去除平均相位差的影响，从而提高相位匹配的准确性。

与θ_l(i)之间的平均相位差可以通过

确定。

N为T时间内的采样次数，N越大，一个周期T内的采样次数越多，性能越好，复杂度也越高，本领域技术人员可以根据实际需要设定N的值。

当f[x]是f[x]＝-x²或f[x]＝-|x|等常用函数时，解调性能较一般，特别是在低信噪比的情况下。为此，本发明实施例进一步根据复数域中欧几里德距离的计算方法通过理论分析，确定采用如下四种函数作为f[x]时，能够得到优秀的解调性能。

函数一：f[x]＝cos(x)，此函数是根据复数域中欧几里德距离的计算方法推导出的度量相位匹配程度的函数；

函数二：

此函数是函数一的改进型，能够消除当|x|≥π时，因cos(x)是周期函数而造成的性能下降。

函数三：

C为设定常数，此函数是函数二的简化型，能够消除cos(x)运算而降低系统复杂度。C建议值是π。

函数四：

此函数是函数二的另一种简化型，也能够消除cos(x)运算而降低系统复杂度。C建议值是π。

进一步，在计算

与θ_l(i)相位匹配程度Q_l时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，可以只使用其解调结果对应的θ_l(i)与

进行相位匹配。

例如，对于M＝1的情况，在解调a[2]时，需要将a[1]、a[2]、a[3]所有可能的取值组合所对应的θ_l(i)与

进行匹配。如果不使用a[1]的解调结果，若是2FSK/2GFSK的情况，a[1]、a[2]、a[3]一共有8种取值组合，因此

需要与8个相位匹配向量进行匹配。如果可以直接使用已经解调的a[1]的解调结果，由于a[1]已知，则a[1]、a[2]、a[3]一共有4种取值组合，

仅需与4个相位匹配向量进行匹配，减少了计算量。

由于a[k-M],a[k-M+1],…,a[k-1]的值已经知道，当调制方式是2FSK/2GFSK时，

仅需与2^2M+1个相位匹配向量中的2^M+1个相位匹配向量进行匹配；当调制方式是4FSK/4GFSK时，

仅需与4^2M+1个相位匹配向量中的4^M+1个相位匹配向量进行匹配。

对于2FSK/2GFSK的情况和4FSK/4GFSK的情况，可以通过如下方法确定a[k]的数字信息向量组{V_l(i)}及相位匹配向量组{θ_l(i)}：

对于2FSK/2GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^{2M +1}个数字信息向量V_l(i),l＝1,2,…,2^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化；

对于4FSK/4GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。

具体的，接收到的FSK/GFSK基带信号在复数域(I/Q平面)可以用以下复函数表示：

在公式(1)中，A是FSK/GFSK基带信号的幅度，h是频率调制指数，a[k]是所传输的数字信息：当调制方式是2FSK/2GFSK时，a[k]∈{1,-1}；当调制方式是4FSK/4FSK时，a[k]∈{3，1,-1，-3}。T是每个数字信息的周期，q(t)定义了每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化，

是接收到的FSK/GFSK基带信号的初始相位，n(t)是接收信号中的干扰和噪声。根据公式(1)，在相位域，接收到的FSK/GFSK基带信号可以表示为：

在公式(2)中，

是干扰和噪声造成的相位噪声。

由于q(t)的变化往往是在超过一个T的时间内完成，由此而产生的码间串扰使得只在一个T的时间范围内解调a[k]的技术(无论该技术是在复数域还是在相位域)都性能不佳。本发明实施例中根据q(t)变化的时间长度，在相位域中使用相位匹配向量组在多个T的时间范围内解调a[k]。

具体方法描述如下：

根据q(t)变化的时间长度，确定在(2M+1)T的时间范围内解调a[k]。M越长性能越好，但解调的复杂度越高，若通过硬件实现，硬件成本也相应提高，一般M常用1或者2。

生成相位匹配向量组：

2FSK/2GFSK：让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^2M+1个数字信息向量V_l(i),l＝1,2,…,2^2M+1。将每个数字信息向量分别代入下式：

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)。

4FSK/4GFSK：让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1。将每个数字信息向量分别代入公式(3)，再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)。

解调操作：

在解调a[k]时，以t＝(k-M)T+iT/N+ε₀,i＝1,2,…,(2M+1)N对接收到的FSK/GFSK基带信号

进行采样得到接收相位向量

其中，ε₀是采样时间误差，其精确值不能确定，但常用的时钟同步算法能够使

为了降低系统复杂度，在解调a[k]时，使用a[k-M],a[k-M+1],…,a[k-1]的解调结果来减少需要匹配的相位向量。由于a[k-M],a[k-M+1],…,a[k-1]的值已经知道，当调制方式是2FSK/2GFSK时，

仅需与2^2M+1个相位匹配向量中的2^M+1个相位匹配向量进行匹配；当调制方式是4FSK/4GFSK时，

仅需与4^2M+1个相位匹配向量中的4^M+1个相位匹配向量进行匹配。

为了消除

和ε₀造成的固定相位差对相位匹配的不良影响，确定

与θ_l(i)之间的平均相位差计算如下：

在得到平均相位差值β_l之后，

与θ_l(i)的相位匹配度计算如下：

在上式中，f[x]是度量相位匹配程度的函数，当f[x]采用如下4种函数时，解调性能较佳。

f[x]＝cos(x)， (6)

在2^M+1个(当调制方式是2FSK/2GFSK时)或4^M+1个(当调制方式是4FSK/4GFSK时)相位匹配计算结果中找到使Q_l最大的l，再根据l从相应的V_l(i)中找到a[k]对应的取值，此取值即为a[k]的解调结果。

本发明实施例在相位域解调FSK/GFSK信号，将大量的复数乘法、加减法运算转化为实数加减法运算，大大降低了解调的系统复杂度。并且，本发明实施例使用相位匹配向量组进行解调，在解调前消除了

和ε₀等固定相位差对相位匹配的不良影响，并设计了简单有效的函数来度量相位匹配程度，从而能够达到与复数域解调技术相近的解调性能。由于在解调a[k]时，使用了a[k-M],a[k-M+1],…,a[k-1]的解调结果，从而可以大幅减少需要匹配的相位向量，进一步降低了使用相位匹配向量组进行解调的系统复杂度。

本发明实施例还相应提供一种FSK/GFSK解调装置，如图2所示，包括：

相位匹配向量组确定单元201，用于在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号；

接收相位向量确定单元202，用于确定接收到的FSK/GFSK基带信号

的接收相位向量

平均相位差确定单元203，用于确定

与θ_l(i)之间的平均相位差β_l；

匹配单元204，用于去除平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值；

解调结果确定单元205，用于确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果。

该装置可以通过软件实现，也可以通过硬件实现，例如集成电路，当通过硬件实现时，各个单元均通过数字电路实现。

例如，复数域到相位域的转化，可以将相位值存在rom里，以复数域输入作为地址，得到相应的相位值；

加减法、乘法、数值比较等运算，可以用硬件加法器，乘法器，比较器等来实现；

cos运算，可以将cos运算结果存在rom里，以相位输入作为地址，得到相应的cos值；

数据存储则可以用移位寄存器实现。

如何用硬件电路实现各单元的功能，是本领域技术人员所熟知，此处不再赘述。

进一步，匹配单元204具体用于：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

更进一步，f[x]具体为：

f[x]＝cos(x)；或者

或者

C为设定常数；或者

C为设定常数。

进一步，匹配单元204计算

与θ_l(i)相位匹配程度Q_l时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，只使用其解调结果对应的θ_l(i)与

进行相位匹配。

进一步，相位匹配向量组确定单元201具体用于：

对于2FSK/2GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^{2M +1}个数字信息向量V_l(i),l＝1,2,…,2^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化；

对于4FSK/4GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。

进一步，平均相位差确定单元203具体用于：

确定

与θ_l(i)之间的平均相位差为

通过该装置，可以在相位域解调FSK/GFSK信号，将大量的复数乘法、加减法运算转化为实数加减法运算，大大降低了解调的系统复杂度。并且，使用相位匹配向量组进行解调，在解调前消除了

和ε₀等固定相位差对相位匹配的不良影响，并设计了简单有效的函数来度量相位匹配程度，从而能够达到与复数域解调技术相近的解调性能。由于在解调a[k]时，使用了a[k-M],a[k-M+1],…,a[k-1]的解调结果，从而可以大幅减少需要匹配的相位向量，进一步降低了使用相位匹配向量组进行解调的系统复杂度。

附图中的流程图和框图，图示了按照本公开各种实施例的方法、装置和计算机可读存储介质的可能实现的体系架构、功能和操作。应当注意，流程图中的每个方框所表示的步骤未必按照标号所示的顺序进行，有时可以基本并行地执行，有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的硬件来实现，或者可以用硬件与计算机指令的组合来实现。

描述于本公开实施例中所涉及到的单元或模块可以通过软件的方式实现，也可以通过硬件的方式来实现。

通过以上对实施例的描述，本领域的技术人员可以清楚地了解到各实施例可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种FSK/GFSK解调方法，其特征在于，包括：

在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号，T是每个数字信息的周期，k为待解调码字编号；

确定接收到的FSK/GFSK基带信号

的接收相位向量

其中a是由a[k-M],a[k-M+1],…a[k-1],a[k],a[k+1],…,a[k+M-1],a[k+M]组成的向量；

确定

与θ_l(i)之间的平均相位差β_l；

去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值；

确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果，

其中，所述去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值，具体包括：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

2.如权利要求1所述的方法，其特征在于，所述f[x]具体为：

f[x]＝cos(x)；或者

或者

C为设定常数；或者

C为设定常数。

3.如权利要求1所述的方法，其特征在于，所述计算

与θ_l(i)相位匹配程度Q_l时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，只使用其解调结果对应的θ_l(i)与

进行相位匹配。

4.如权利要求1所述的方法，其特征在于，所述在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，具体包括：

对于2FSK/2GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1，生成2^2M+1个数字信息向量V_l(i),l＝1,2,…,2^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成2^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化，h是频率调制指数；

对于4FSK/4GFSK的情况，让a[1],a[2],…,a[2M],a[2M+1]分别取值3，1，-1或-3，生成4^2M+1个数字信息向量V_l(i),l＝1,2,…,4^2M+1，将每个数字信息向量分别代入

再以t＝iT/N,i＝1,2,…,(2M+1)N对θ(t,a)进行采样生成4^2M+1个相位匹配向量θ_l(i)，其中，q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化，h是频率调制指数。

5.如权利要求1所述的方法，其特征在于，所述确定

与θ_l(i)之间的平均相位差β_l，具体包括：

确定

与θ_l(i)之间的平均相位差为

6.一种FSK/GFSK解调装置，其特征在于，包括：

相位匹配向量组确定单元，用于在(2M+1)T的时间范围内，确定待解调码字a[k]的数字信息向量组{V_l(i)}及对应的相位匹配向量组{θ_l(i)}，M为整数且M≥1，l为相位匹配向量编号，i为采样点编号，T是每个数字信息的周期，k为待解调码字编号；

接收相位向量确定单元，用于确定接收到的FSK/GFSK基带信号

的接收相位向量

其中a是由a[k-M],a[k-M+1],…a[k-1],a[k],a[k+1],…,a[k+M-1],a[k+M]组成的向量；

平均相位差确定单元，用于确定

与θ_l(i)之间的平均相位差β_l；

匹配单元，用于去除所述平均相位差β_l的影响后，计算

与θ_l(i)相位匹配程度Q_l，并确定使相位匹配程度Q_l最大的l值；

解调结果确定单元，用于确定该l值在数字信息向量V_l(i)中对应的a[k]为解调结果，

其中，所述匹配单元具体用于：

确定相位匹配程度为

其中，N为T时间内的采样次数，f[x]是度量相位匹配程度的函数；

确定使相位匹配程度Q_l最大的l值。

7.如权利要求6所述的装置，其特征在于，所述f[x]具体为：

f[x]＝cos(x)；或者

或者

C为设定常数；或者

C为设定常数。

8.如权利要求6所述的装置，其特征在于，所述匹配单元计算

与θ_l(i)相位匹配程度Q_l时，对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1]，只使用其解调结果对应的θ_l(i)与

进行相位匹配。

Images