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

一种fsk/gfsk解调方法及装置 Download PDF

Info

Publication number
CN108989256B
CN108989256B CN201811024709.XA CN201811024709A CN108989256B CN 108989256 B CN108989256 B CN 108989256B CN 201811024709 A CN201811024709 A CN 201811024709A CN 108989256 B CN108989256 B CN 108989256B
Authority
CN
China
Prior art keywords
phase matching
phase
determining
theta
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811024709.XA
Other languages
English (en)
Other versions
CN108989256A (zh
Inventor
聂宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TAILING MICROELECTRONICS (SHANGHAI) CO Ltd
Original Assignee
TAILING MICROELECTRONICS (SHANGHAI) CO Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by TAILING MICROELECTRONICS (SHANGHAI) CO Ltd filed Critical TAILING MICROELECTRONICS (SHANGHAI) CO Ltd
Priority to CN201811024709.XA priority Critical patent/CN108989256B/zh
Publication of CN108989256A publication Critical patent/CN108989256A/zh
Priority to US16/559,961 priority patent/US10763788B2/en
Priority to EP19195346.2A priority patent/EP3621259B1/en
Application granted granted Critical
Publication of CN108989256B publication Critical patent/CN108989256B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/10Frequency-modulated carrier systems, i.e. using frequency-shift keying
    • H04L27/14Demodulator circuits; Receiver circuits
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03DDEMODULATION OR TRANSFERENCE OF MODULATION FROM ONE CARRIER TO ANOTHER
    • H03D3/00Demodulation of angle-, frequency- or phase- modulated oscillations
    • H03D3/02Demodulation of angle-, frequency- or phase- modulated oscillations by detecting phase difference between two signals obtained from input signal
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/10Frequency-modulated carrier systems, i.e. using frequency-shift keying
    • H04L27/14Demodulator circuits; Receiver circuits
    • H04L27/144Demodulator circuits; Receiver circuits with demodulation using spectral properties of the received signal, e.g. by using frequency selective- or frequency sensitive elements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/20Modulator circuits; Transmitter circuits
    • H04L27/2003Modulator circuits; Transmitter circuits for continuous phase modulation
    • H04L27/2007Modulator circuits; Transmitter circuits for continuous phase modulation in which the phase change within each symbol period is constrained
    • H04L27/2017Modulator circuits; Transmitter circuits for continuous phase modulation in which the phase change within each symbol period is constrained in which the phase changes are non-linear, e.g. generalized and Gaussian minimum shift keying, tamed frequency modulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/22Demodulator circuits; Receiver circuits
    • H04L27/233Demodulator circuits; Receiver circuits using non-coherent demodulation
    • H04L27/2338Demodulator circuits; Receiver circuits using non-coherent demodulation using sampling

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Power Engineering (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明提供了一种FSK/GFSK解调方法及装置,涉及通信技术,在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},确定接收到的FSK/GFSK基带信号
Figure DDA0001788211710000011
的接收相位向量
Figure DDA0001788211710000012
后,还要确定
Figure DDA0001788211710000013
与θl(i)之间的平均相位差βl,并在去除平均相位差βl的影响后,计算
Figure DDA0001788211710000014
与θl(i)相位匹配程度Ql,找出使相位匹配程度Ql最大的l值,进而确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。由于在进行相位匹配时,去除了平均相位差的影响,增大了相位匹配的准确性,提高了相位域解调技术的性能。

Description

一种FSK/GFSK解调方法及装置
技术领域
本发明属于数据处理领域,具体涉及一种FSK/GFSK解调方法及装置。
背景技术
现有基于匹配滤波器组的FSK/GFSK解调方法主要是在复数域(I\Q平面)完成解调操作,需要进行大量的复数乘法、复数加减法运算,系统复杂度高。
由于FSK/GFSK信号具有恒定的包络,所传输的数字信息全部包含在信号的相位中,在相位域完成解调操作,能够将大量的复数运算转化为实数加减法运算,从而可以有效的降低系统复杂度。
然而,本发明的发明人发现,现有的相位域解调技术或者没有使用相位匹配滤波器组(即本发明中的相位匹配向量组)来进行解调操作,或者在使用相位匹配向量组时存在技术缺陷,使得现有的相位域解调技术的性能远低于复数域解调技术的性能。
发明内容
针对上面提到的相位域解调技术的性能较低的问题,提出了一种FSK/GFSK解调方法及装置,以提高相位域解调技术的性能。
依据本发明的第一方面,提供了一种FSK/GFSK解调方法,包括:
在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号;
确定接收到的FSK/GFSK基带信号
Figure BDA0001788211690000021
的接收相位向量
Figure BDA0001788211690000022
确定
Figure BDA0001788211690000023
与θl(i)之间的平均相位差βl
去除所述平均相位差βl的影响后,计算
Figure BDA0001788211690000024
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值;
确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。
进一步,所述去除所述平均相位差βl的影响后,计算
Figure BDA0001788211690000025
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值,具体包括:
确定相位匹配程度为
Figure BDA0001788211690000026
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
更进一步,所述f[x]具体为:
f[x]=cos(x);或者
Figure BDA0001788211690000027
或者
Figure BDA0001788211690000028
C为设定常数;或者
Figure BDA0001788211690000031
C为设定常数。
更进一步,所述确定
Figure BDA0001788211690000032
与θl(i)相位匹配程度时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],只使用其解调结果对应的θl(i)与
Figure BDA0001788211690000033
进行相位匹配。
进一步,所述在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},具体包括:
对于2FSK/2GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1,生成22M +1个数字信息向量Vl(i),l=1,2,…,22M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000034
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化;
对于4FSK/4GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000035
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。
进一步,所述确定
Figure BDA0001788211690000036
与θl(i)之间的平均相位差βl,具体包括:
确定
Figure BDA0001788211690000037
与θl(i)之间的平均相位差为
Figure BDA0001788211690000038
依据本发明的第二方面,提供一种FSK/GFSK解调装置,包括:
相位匹配向量组确定单元,用于在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号;
接收相位向量确定单元,用于确定接收到的FSK/GFSK基带信号
Figure BDA0001788211690000041
的接收相位向量
Figure BDA0001788211690000042
平均相位差确定单元,用于确定
Figure BDA0001788211690000043
与θl(i)之间的平均相位差βl
匹配单元,用于去除所述平均相位差βl的影响后,计算
Figure BDA0001788211690000044
与θl(i)相位匹配程度Ql,并确定相位匹配程度Ql最大的l值;
解调结果确定单元,用于确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。
进一步,所述匹配单元具体用于:
确定相位匹配程度为
Figure BDA0001788211690000045
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
更进一步,所述f[x]具体为:
f[x]=cos(x);或者
Figure BDA0001788211690000046
或者
Figure BDA0001788211690000047
C为设定常数;或者
Figure BDA0001788211690000048
C为设定常数。
更进一步,所述匹配单元计算
Figure BDA0001788211690000049
与θl(i)相位匹配程度Ql时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],只使用其解调结果对应的θl(i)与
Figure BDA0001788211690000051
进行相位匹配。
进一步,相位匹配向量组确定单元具体用于:
对于2FSK/2GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1,生成22M +1个数字信息向量Vl(i),l=1,2,…,22M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000052
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化;
对于4FSK/4GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000053
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。
进一步,所述平均相位差确定单元具体用于:
确定
Figure BDA0001788211690000054
与θl(i)之间的平均相位差为
Figure BDA0001788211690000055
本发明实施例提供一种FSK/GFSK解调方法及装置,在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},确定接收到的FSK/GFSK基带信号
Figure BDA0001788211690000056
的接收相位向量
Figure BDA0001788211690000057
后,还要确定
Figure BDA0001788211690000058
与θl(i)之间的平均相位差βl,并在去除平均相位差βl的影响后,计算
Figure BDA0001788211690000059
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值,进而确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。由于在进行相位匹配时,去除了平均相位差的影响,增大了相位匹配的准确性,提高了相位域解调技术的性能。
应当理解,上述说明仅是本发明技术方案的概述,以便能够更清楚地了解本发明的技术手段,从而可依照说明书的内容予以实施。为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举说明本发明的具体实施方式。
附图说明
通过阅读下文的示例性实施例的详细描述,本领域普通技术人员将明白本文所述的有点和益处以及其他优点和益处。附图仅用于示出示例性实施例的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的标号表示相同的部件。在附图中:
图1为本发明实施例提供的FSK/GFSK解调方法流程图;
图2为本发明实施例提供的FSK/GFSK解调装置结构示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
在本发明中,应理解,诸如“包括”或“具有”等术语旨在指示本说明书中所公开的特征、数字、步骤、行为、部件、部分或其组合的存在,并且不旨在排除一个或多个其他特征、数字、步骤、行为、部件、部分或其组合存在的可能性。
另外还需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
如图1所示,本发明实施例提供的FSK/GFSK解调方法,包括:
步骤S101、在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号;
步骤S102、确定接收到的FSK/GFSK基带信号
Figure BDA0001788211690000071
的接收相位向量
Figure BDA0001788211690000072
步骤S103、确定
Figure BDA0001788211690000073
与θl(i)之间的平均相位差βl
步骤S104、去除平均相位差βl的影响后,计算
Figure BDA0001788211690000074
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值;
步骤S105、确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。
由于在进行相位匹配时,去除了平均相位差的影响,增大了相位匹配的准确性,提高了相位域解调技术的性能。
由于在相位域解调FSK/GFSK信号,将大量的复数乘法、加减法运算转化为实数加减法运算,大大降低了解调的系统复杂度,并且,使用相位匹配向量组进行解调,在解调前消除了平均相位差对相位匹配的不良影响,增大了相位匹配的准确性,提高了相位域解调技术的性能。
进一步,步骤S104,去除平均相位差βl的影响后,计算
Figure BDA0001788211690000081
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值,具体包括:
确定相位匹配程度为
Figure BDA0001788211690000082
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
通过
Figure BDA0001788211690000083
确定相位匹配程度,可以较好的去除平均相位差的影响,从而提高相位匹配的准确性。
Figure BDA0001788211690000084
与θl(i)之间的平均相位差可以通过
Figure BDA0001788211690000085
Figure BDA0001788211690000086
确定。
N为T时间内的采样次数,N越大,一个周期T内的采样次数越多,性能越好,复杂度也越高,本领域技术人员可以根据实际需要设定N的值。
当f[x]是f[x]=-x2或f[x]=-|x|等常用函数时,解调性能较一般,特别是在低信噪比的情况下。为此,本发明实施例进一步根据复数域中欧几里德距离的计算方法通过理论分析,确定采用如下四种函数作为f[x]时,能够得到优秀的解调性能。
函数一:f[x]=cos(x),此函数是根据复数域中欧几里德距离的计算方法推导出的度量相位匹配程度的函数;
函数二:
Figure BDA0001788211690000091
此函数是函数一的改进型,能够消除当|x|≥π时,因cos(x)是周期函数而造成的性能下降。
函数三:
Figure BDA0001788211690000092
C为设定常数,此函数是函数二的简化型,能够消除cos(x)运算而降低系统复杂度。C建议值是π。
函数四:
Figure BDA0001788211690000093
此函数是函数二的另一种简化型,也能够消除cos(x)运算而降低系统复杂度。C建议值是π。
进一步,在计算
Figure BDA0001788211690000094
与θl(i)相位匹配程度Ql时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],可以只使用其解调结果对应的θl(i)与
Figure BDA0001788211690000095
进行相位匹配。
例如,对于M=1的情况,在解调a[2]时,需要将a[1]、a[2]、a[3]所有可能的取值组合所对应的θl(i)与
Figure BDA0001788211690000096
进行匹配。如果不使用a[1]的解调结果,若是2FSK/2GFSK的情况,a[1]、a[2]、a[3]一共有8种取值组合,因此
Figure BDA0001788211690000097
需要与8个相位匹配向量进行匹配。如果可以直接使用已经解调的a[1]的解调结果,由于a[1]已知,则a[1]、a[2]、a[3]一共有4种取值组合,
Figure BDA0001788211690000098
仅需与4个相位匹配向量进行匹配,减少了计算量。
由于a[k-M],a[k-M+1],…,a[k-1]的值已经知道,当调制方式是2FSK/2GFSK时,
Figure BDA0001788211690000099
仅需与22M+1个相位匹配向量中的2M+1个相位匹配向量进行匹配;当调制方式是4FSK/4GFSK时,
Figure BDA00017882116900000910
仅需与42M+1个相位匹配向量中的4M+1个相位匹配向量进行匹配。
对于2FSK/2GFSK的情况和4FSK/4GFSK的情况,可以通过如下方法确定a[k]的数字信息向量组{Vl(i)}及相位匹配向量组{θl(i)}:
对于2FSK/2GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1,生成22M +1个数字信息向量Vl(i),l=1,2,…,22M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000101
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化;
对于4FSK/4GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000102
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。
具体的,接收到的FSK/GFSK基带信号在复数域(I/Q平面)可以用以下复函数表示:
Figure BDA0001788211690000103
在公式(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基带信号的相位变化,
Figure BDA0001788211690000111
是接收到的FSK/GFSK基带信号的初始相位,n(t)是接收信号中的干扰和噪声。根据公式(1),在相位域,接收到的FSK/GFSK基带信号可以表示为:
Figure BDA0001788211690000112
在公式(2)中,
Figure BDA0001788211690000113
是干扰和噪声造成的相位噪声。
由于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,生成22M+1个数字信息向量Vl(i),l=1,2,…,22M+1。将每个数字信息向量分别代入下式:
Figure BDA0001788211690000114
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i)。
4FSK/4GFSK:让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1。将每个数字信息向量分别代入公式(3),再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i)。
解调操作:
在解调a[k]时,以t=(k-M)T+iT/N+ε0,i=1,2,…,(2M+1)N对接收到的FSK/GFSK基带信号
Figure BDA0001788211690000121
进行采样得到接收相位向量
Figure BDA0001788211690000122
其中,ε0是采样时间误差,其精确值不能确定,但常用的时钟同步算法能够使
Figure BDA0001788211690000123
为了降低系统复杂度,在解调a[k]时,使用a[k-M],a[k-M+1],…,a[k-1]的解调结果来减少需要匹配的相位向量。由于a[k-M],a[k-M+1],…,a[k-1]的值已经知道,当调制方式是2FSK/2GFSK时,
Figure BDA0001788211690000124
仅需与22M+1个相位匹配向量中的2M+1个相位匹配向量进行匹配;当调制方式是4FSK/4GFSK时,
Figure BDA0001788211690000125
仅需与42M+1个相位匹配向量中的4M+1个相位匹配向量进行匹配。
为了消除
Figure BDA0001788211690000126
和ε0造成的固定相位差对相位匹配的不良影响,确定
Figure BDA0001788211690000127
与θl(i)之间的平均相位差计算如下:
Figure BDA0001788211690000128
在得到平均相位差值βl之后,
Figure BDA0001788211690000129
与θl(i)的相位匹配度计算如下:
Figure BDA00017882116900001210
在上式中,f[x]是度量相位匹配程度的函数,当f[x]采用如下4种函数时,解调性能较佳。
f[x]=cos(x), (6)
Figure BDA0001788211690000131
Figure BDA0001788211690000132
Figure BDA0001788211690000133
在2M+1个(当调制方式是2FSK/2GFSK时)或4M+1个(当调制方式是4FSK/4GFSK时)相位匹配计算结果中找到使Ql最大的l,再根据l从相应的Vl(i)中找到a[k]对应的取值,此取值即为a[k]的解调结果。
本发明实施例在相位域解调FSK/GFSK信号,将大量的复数乘法、加减法运算转化为实数加减法运算,大大降低了解调的系统复杂度。并且,本发明实施例使用相位匹配向量组进行解调,在解调前消除了
Figure BDA0001788211690000134
和ε0等固定相位差对相位匹配的不良影响,并设计了简单有效的函数来度量相位匹配程度,从而能够达到与复数域解调技术相近的解调性能。由于在解调a[k]时,使用了a[k-M],a[k-M+1],…,a[k-1]的解调结果,从而可以大幅减少需要匹配的相位向量,进一步降低了使用相位匹配向量组进行解调的系统复杂度。
本发明实施例还相应提供一种FSK/GFSK解调装置,如图2所示,包括:
相位匹配向量组确定单元201,用于在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号;
接收相位向量确定单元202,用于确定接收到的FSK/GFSK基带信号
Figure BDA0001788211690000141
的接收相位向量
Figure BDA0001788211690000142
平均相位差确定单元203,用于确定
Figure BDA0001788211690000143
与θl(i)之间的平均相位差βl
匹配单元204,用于去除平均相位差βl的影响后,计算
Figure BDA0001788211690000144
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值;
解调结果确定单元205,用于确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果。
该装置可以通过软件实现,也可以通过硬件实现,例如集成电路,当通过硬件实现时,各个单元均通过数字电路实现。
例如,复数域到相位域的转化,可以将相位值存在rom里,以复数域输入作为地址,得到相应的相位值;
加减法、乘法、数值比较等运算,可以用硬件加法器,乘法器,比较器等来实现;
cos运算,可以将cos运算结果存在rom里,以相位输入作为地址,得到相应的cos值;
数据存储则可以用移位寄存器实现。
如何用硬件电路实现各单元的功能,是本领域技术人员所熟知,此处不再赘述。
进一步,匹配单元204具体用于:
确定相位匹配程度为
Figure BDA0001788211690000145
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
更进一步,f[x]具体为:
f[x]=cos(x);或者
Figure BDA0001788211690000151
或者
Figure BDA0001788211690000152
C为设定常数;或者
Figure BDA0001788211690000153
C为设定常数。
进一步,匹配单元204计算
Figure BDA0001788211690000154
与θl(i)相位匹配程度Ql时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],只使用其解调结果对应的θl(i)与
Figure BDA0001788211690000155
进行相位匹配。
进一步,相位匹配向量组确定单元201具体用于:
对于2FSK/2GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1,生成22M +1个数字信息向量Vl(i),l=1,2,…,22M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000156
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化;
对于4FSK/4GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1,将每个数字信息向量分别代入
Figure BDA0001788211690000157
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化。
进一步,平均相位差确定单元203具体用于:
确定
Figure BDA0001788211690000161
与θl(i)之间的平均相位差为
Figure BDA0001788211690000162
通过该装置,可以在相位域解调FSK/GFSK信号,将大量的复数乘法、加减法运算转化为实数加减法运算,大大降低了解调的系统复杂度。并且,使用相位匹配向量组进行解调,在解调前消除了
Figure BDA0001788211690000163
和ε0等固定相位差对相位匹配的不良影响,并设计了简单有效的函数来度量相位匹配程度,从而能够达到与复数域解调技术相近的解调性能。由于在解调a[k]时,使用了a[k-M],a[k-M+1],…,a[k-1]的解调结果,从而可以大幅减少需要匹配的相位向量,进一步降低了使用相位匹配向量组进行解调的系统复杂度。
附图中的流程图和框图,图示了按照本公开各种实施例的方法、装置和计算机可读存储介质的可能实现的体系架构、功能和操作。应当注意,流程图中的每个方框所表示的步骤未必按照标号所示的顺序进行,有时可以基本并行地执行,有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的硬件来实现,或者可以用硬件与计算机指令的组合来实现。
描述于本公开实施例中所涉及到的单元或模块可以通过软件的方式实现,也可以通过硬件的方式来实现。
通过以上对实施例的描述,本领域的技术人员可以清楚地了解到各实施例可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种FSK/GFSK解调方法,其特征在于,包括:
在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号,T是每个数字信息的周期,k为待解调码字编号;
确定接收到的FSK/GFSK基带信号
Figure FDA0002795356710000011
的接收相位向量
Figure FDA0002795356710000012
其中a是由a[k-M],a[k-M+1],…a[k-1],a[k],a[k+1],…,a[k+M-1],a[k+M]组成的向量;
确定
Figure FDA0002795356710000013
与θl(i)之间的平均相位差βl
去除所述平均相位差βl的影响后,计算
Figure FDA0002795356710000014
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值;
确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果,
其中,所述去除所述平均相位差βl的影响后,计算
Figure FDA0002795356710000015
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值,具体包括:
确定相位匹配程度为
Figure FDA0002795356710000016
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
2.如权利要求1所述的方法,其特征在于,所述f[x]具体为:
f[x]=cos(x);或者
Figure FDA0002795356710000017
或者
Figure FDA0002795356710000018
C为设定常数;或者
Figure FDA0002795356710000019
C为设定常数。
3.如权利要求1所述的方法,其特征在于,所述计算
Figure FDA00027953567100000110
与θl(i)相位匹配程度Ql时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],只使用其解调结果对应的θl(i)与
Figure FDA00027953567100000111
进行相位匹配。
4.如权利要求1所述的方法,其特征在于,所述在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},具体包括:
对于2FSK/2GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值1或-1,生成22M+1个数字信息向量Vl(i),l=1,2,…,22M+1,将每个数字信息向量分别代入
Figure FDA0002795356710000021
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成22M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化,h是频率调制指数;
对于4FSK/4GFSK的情况,让a[1],a[2],…,a[2M],a[2M+1]分别取值3,1,-1或-3,生成42M+1个数字信息向量Vl(i),l=1,2,…,42M+1,将每个数字信息向量分别代入
Figure FDA0002795356710000022
再以t=iT/N,i=1,2,…,(2M+1)N对θ(t,a)进行采样生成42M+1个相位匹配向量θl(i),其中,q(t)为每个数字信息在时间上造成的FSK/GFSK基带信号的相位变化,h是频率调制指数。
5.如权利要求1所述的方法,其特征在于,所述确定
Figure FDA0002795356710000023
与θl(i)之间的平均相位差βl,具体包括:
确定
Figure FDA0002795356710000024
与θl(i)之间的平均相位差为
Figure FDA0002795356710000025
6.一种FSK/GFSK解调装置,其特征在于,包括:
相位匹配向量组确定单元,用于在(2M+1)T的时间范围内,确定待解调码字a[k]的数字信息向量组{Vl(i)}及对应的相位匹配向量组{θl(i)},M为整数且M≥1,l为相位匹配向量编号,i为采样点编号,T是每个数字信息的周期,k为待解调码字编号;
接收相位向量确定单元,用于确定接收到的FSK/GFSK基带信号
Figure FDA0002795356710000028
的接收相位向量
Figure FDA0002795356710000026
其中a是由a[k-M],a[k-M+1],…a[k-1],a[k],a[k+1],…,a[k+M-1],a[k+M]组成的向量;
平均相位差确定单元,用于确定
Figure FDA0002795356710000027
与θl(i)之间的平均相位差βl
匹配单元,用于去除所述平均相位差βl的影响后,计算
Figure FDA0002795356710000031
与θl(i)相位匹配程度Ql,并确定使相位匹配程度Ql最大的l值;
解调结果确定单元,用于确定该l值在数字信息向量Vl(i)中对应的a[k]为解调结果,
其中,所述匹配单元具体用于:
确定相位匹配程度为
Figure FDA0002795356710000032
其中,N为T时间内的采样次数,f[x]是度量相位匹配程度的函数;
确定使相位匹配程度Ql最大的l值。
7.如权利要求6所述的装置,其特征在于,所述f[x]具体为:
f[x]=cos(x);或者
Figure FDA0002795356710000033
或者
Figure FDA0002795356710000034
C为设定常数;或者
Figure FDA0002795356710000035
C为设定常数。
8.如权利要求6所述的装置,其特征在于,所述匹配单元计算
Figure FDA0002795356710000036
与θl(i)相位匹配程度Ql时,对于已知解调结果的a[k-M],a[k-M+1],…,a[k-1],只使用其解调结果对应的θl(i)与
Figure FDA0002795356710000037
进行相位匹配。
CN201811024709.XA 2018-09-04 2018-09-04 一种fsk/gfsk解调方法及装置 Active CN108989256B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201811024709.XA CN108989256B (zh) 2018-09-04 2018-09-04 一种fsk/gfsk解调方法及装置
US16/559,961 US10763788B2 (en) 2018-09-04 2019-09-04 Method and device for FSK/GFSK demodulation
EP19195346.2A EP3621259B1 (en) 2018-09-04 2019-09-04 Method and device for fsk/gfsk demodulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811024709.XA CN108989256B (zh) 2018-09-04 2018-09-04 一种fsk/gfsk解调方法及装置

Publications (2)

Publication Number Publication Date
CN108989256A CN108989256A (zh) 2018-12-11
CN108989256B true CN108989256B (zh) 2021-03-19

Family

ID=64548263

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811024709.XA Active CN108989256B (zh) 2018-09-04 2018-09-04 一种fsk/gfsk解调方法及装置

Country Status (3)

Country Link
US (1) US10763788B2 (zh)
EP (1) EP3621259B1 (zh)
CN (1) CN108989256B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2024309B1 (en) * 2019-11-26 2021-08-30 Univ Delft Tech Frequency Shift Keying signal demodulator system and method
CN113079120B (zh) * 2021-03-24 2022-05-17 湖南坤雷科技有限公司 一种2fsk调制信号的识别方法及装置
CN113098811B (zh) * 2021-04-01 2022-06-28 高拓讯达(北京)科技有限公司 一种gfsk信号的解调解码方法及解调解码装置
CN114640562B (zh) * 2022-03-16 2023-05-30 中山大学 一种cpfsk/gfsk信号非相干解调方法
CN115604062B (zh) * 2022-10-08 2024-04-12 杭州万高科技股份有限公司 Gfsk通信模式内的双比特组解调方法及解调器

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106161309A (zh) * 2015-03-27 2016-11-23 飞思卡尔半导体公司 解调频移键控调制的输入信号

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1643704B1 (en) * 2004-09-30 2007-10-31 Sony Deutschland GmbH Frequency-shift-keying demodulator and method of frequency-shift-keying
FR2932344B1 (fr) * 2008-06-06 2010-06-04 Thales Sa Procede de poursuite de la phase d'un signal module par modulation a phase continue et dispositif de synchronisation mettant en oeuvre le procede.
JP5347792B2 (ja) * 2009-07-16 2013-11-20 ソニー株式会社 信号処理装置、信号処理方法、及び、受信システム
CN101969329B (zh) * 2010-09-30 2014-05-21 中国科学院国家天文台 天线组阵中基于循环互相关的信号相位差估计装置与方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106161309A (zh) * 2015-03-27 2016-11-23 飞思卡尔半导体公司 解调频移键控调制的输入信号

Also Published As

Publication number Publication date
US10763788B2 (en) 2020-09-01
US20200076370A1 (en) 2020-03-05
EP3621259A1 (en) 2020-03-11
EP3621259B1 (en) 2021-01-27
CN108989256A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
CN108989256B (zh) 一种fsk/gfsk解调方法及装置
JPH0621982A (ja) Gmsk信号復調方法及びその装置
KR101828790B1 (ko) 주파수 편이 변조 신호의 수신 방법 및 장치
JPH0370231A (ja) 相関による同期方法
JP5354293B2 (ja) 位相同期装置および位相同期方法
CN116633456B (zh) 一种宽带射频系统信号的校准方法、装置及电子设备
CN107708145B (zh) 一种同步检测的方法及同步检测设备
JP3226561B2 (ja) Fsk信号受信回路
JPH05183592A (ja) 周波数変換回路、位相比較回路、およびこれらを備えた遅延検波復調装置
WO2020040889A1 (en) System-on-a-chip for reception of telemetry messages over a radio frequency channel
CN108353066B (zh) 用于载波频率偏移校正的设备和方法及其存储介质
CN112152758B (zh) 导航通信方法及装置
CN113765838A (zh) 一种dpsk信号的解调方法、装置、设备及存储介质
WO2021255108A1 (en) Demodulating frequency-modulated signals
CN107800656B (zh) 一种对于iq不平衡进行补偿的方法及设备
Webber et al. Implementing a/4 shift D-QPSK baseband modem using the TMS320C50
CN112448902A (zh) 一种新型发射机iq失衡估计和补偿方法及装置
CN114124632B (zh) 用于频移键控信号的自适应解调系统、自适应解调方法
JP2001177587A (ja) ディジタル変復調の同期方式
CN114184837B (zh) 一种基于Cordic算法的瞬时测频方法
JPH10308785A (ja) Tdmaデータ受信装置
CN116505978B (zh) 一种蓝牙信号的处理方法、装置、电子设备及存储介质
CN114363134B (zh) 一种信号解调方法、装置、电子设备及存储介质
CN111147417B (zh) 载波频偏估计的实现方法、装置、终端及存储介质
JP2739318B2 (ja) 最尤受信機

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: Building 3, no.1500 Zuchongzhi Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai, 201203

Applicant after: Tailing Microelectronics (Shanghai) Co.,Ltd.

Address before: Building 3, no.1500 Zuchongzhi Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai, 201203

Applicant before: TELINK SEMICONDUCTOR (SHANGHAI) Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant