CN110691051A - 一种基于fft的gmsk信号频偏估计算法 - Google Patents
一种基于fft的gmsk信号频偏估计算法 Download PDFInfo
- Publication number
- CN110691051A CN110691051A CN201910934160.6A CN201910934160A CN110691051A CN 110691051 A CN110691051 A CN 110691051A CN 201910934160 A CN201910934160 A CN 201910934160A CN 110691051 A CN110691051 A CN 110691051A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- frequency offset
- fft
- gmsk
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
- H04L27/14—Demodulator circuits; Receiver circuits
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公基于FFT的GMSK信号频偏估计算法,解调时,去除GMSK调制信号的相位信息,利用FFT运算得到最大谱线与次谱线,计算出频移因子;设定判断偏移量是否处于中心区域的门限值,当根据门限值判定位于中心区域时,直接利用得出的频移因子计算频率偏移量;否则用频谱细化求得最大谱线左右两侧相邻谱线中心位置频率的幅值作为修正方向的判据,将待测信号进行频谱搬移,对频移因子进行修正之后再计算频率偏移量。本发明获得更高的估计精度且不增加不必要的计算量,从而使得在低信噪比环境和在待估频偏处于量化频率点附近的情况下,都能获得良好的估计精度。
Description
技术领域
本发明涉及信号处理技术领域,特别是涉及一种基于FFT的GMSK信号频偏估计算法。
背景技术
对GMSK信号解调的过程中,由于多普勒效应以及本振偏移等因素,造成接收的信号存在一定的频率偏移量。在相干解调过程中,解调性能将受到该频率偏差的影响,降低***的可靠性。现有的频偏估计算法可以分为时域类和频域类。时域类中,Kay算法具有实现简单的优点,但由于工作门限高而不适用于低信噪比条件。Fitz算法具有很好的低信噪比性能,但是频率估计范围有限。频域类中,快速傅里叶变换(FFT)法能直接检测频偏的绝对值,具有较快的速度,适用于实时性***。但是,FFT的栅栏效应将估计范围限制在了为采样频率的一半,且估计的精度与FFT的点数直接相关。
Rife算法是一种利用FFT实现的算法,Rife算法的估计精度随着量化频率所处的范围产生很大的波动。Quinn算法具有低的误判率,但量化频率点附近的性能恶化仍较严重。M-Rife算法是修正的Rife算法,该算法在低信噪比下不存在发散问题,但有时需进行二次移频。频域迭代插值(IIN)算法,通过不断迭代补偿来提高精度,该算法的计算量较大。
发明内容
本发明的目的是针对现有技术中存在的技术缺陷,而提供一种基于FFT的GMSK信号频偏估计算法。
为实现本发明的目的所采用的技术方案是:
一种基于FFT的GMSK信号频偏估计算法,包括步骤如下:
在解调时,先去除GMSK调制信号的相位信息,再利用反馈修正Rife算法进行频偏估计;利用反馈修正Rife算法进行频偏估计的步骤如下:
利用FFT运算得到最大谱线与次谱线,计算出频移因子;设定判断偏移量是否处于中心区域的门限值,当根据门限值判定位于中心区域时,直接利用得出的频移因子计算频率偏移量;否则,利用频谱细化求得最大谱线左右两侧相邻谱线中心位置频率的幅值作为修正方向的判据,将待测信号进行频谱搬移,对频移因子进行修正之后再计算频率偏移量。
本发明通过在GMSK的解调的中频段构造频率信息的辅助函数,去除调制相位信息的影响;利用阈值判定与反馈技术提出反馈修正Rife算法,解决Rife算法的精度在低信噪比环境下和频偏值位于量化频率附近时的恶化现象,获得更高的估计精度且不增加不必要的计算量,从而使得在低信噪比环境和在待估频偏处于量化频率点附近的情况下,都能获得良好的估计精度。
附图说明
图1是接收机中的数字解调块的结构原理图;
图2为本发明的基于FFT的GMSK信号频偏估计算法的流程图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明在去除GMSK调制信号的相位信息后,首先利用FFT得到最大谱线与次谱线,并利用比值计算出频移因子;设定判断偏移量是否处于中心区域的门限值,当根据门限值判定位于中心区域时,直接利用得出的频移因子计算频率偏移量,当判定不位于中心区域时,利用频谱细化求得最大谱线左右两侧相邻谱线中心位置频率的幅值作为修正方向的判据,将待测信号进行频谱搬移,对频移因子进行修正之后再计算频率偏移量,以获得更高的估计精度。
如图1所示,在接收机中,数字解调块由数字混频器、位同步、载波同步、解调、译码几个部分组成。载波同步通过频偏估计和频偏补偿实现,载波同步后的信号进行后续的解调及译码恢复出发送端的数字信号。
假设GMSK调制后的信号通过加行高斯白噪声信道(AWGN)传输并且完成定时恢复。依次通过混频器和带通滤波器后,接收到的信号表示为:
其中,f0和fd分别代表中间频率和多普勒频率偏移,θ是GMSK信号调制后的相位信息φ(t,α)。n(t)是高斯白噪声。GMSK调制信号可以表示为:
F(t)=ejθ(t) (2)
其中θ(t)是承载GMSK信号的相位信息。
为去除相位分量的影响,定义仅包含频率分量的辅助方程:
x(t)=rI(t)FI(t)+rQ(t)FQ(t) (3)
其中rI(t)和rQ(t)分别表示r(t)的同相分量和正交分量。同样,FI(t)和FQ(t)分别代表F(t)的同相分量和正交分量。
信号x(t)化简后可表示为:
x(t)=cos(2πfot+2πfdt)+w(t),w(t)=FI(t)nI(t)+FQ(t)nQ(t) (4)。
得到信号的简化模型后,利用提出的反馈修正Rife算法进行频偏估计。算法的流程图如图2所示。
首先通过步骤①对信号x(t)在[0,T]内进行周期为Ts的采样得到长度为N的序列x(n):
x(n)=cos(2πf0nTs+2fdnTs)+w(n),0≤n≤N-1 (5)
其中n∈[0,1,…,N]。通过步骤②对x(n)做N点FFT运算得到X(k)。w(n)为实部与虚部相互独立的,方差为2σ2的零均值复高斯白噪声,利用实信号FFT的对称性,只保留离散频谱的前N/2点,即只考虑频谱的正频率成分,可得:
当|X(k+1)|>|X(k-1)|时,r=1,否则,r=-1,
利用线性调频z变换(CZT)对k1-1与k1+1之间的频谱进行频谱细化:
由于|X(k1+0.5)|和|X(k1-0.5)|的值比|X(k1+1)|和|X(k1-1)|大,具有更强的抗噪声能力。判定不处于中心区域时r容易发生误判,因此通过比较|X(k1+0.5)|和|X(k1-0.5)|的大小再次得到r,如与第一次判别不同,则以此次为准。并返回步骤④再次通过式(7)计算的值,再通过步骤⑦,利用式(11)得到频率偏移量表示为:
当|X(k+0.5)|>|X(k-0.5)|时,r=1,否则,r=-1。
本发明通过构造频率辅助方程去除GMSK信号相位信息的影响,然后利用反馈修正Rife算法进行估计,通过频谱细化减小噪声环境对估计效果的影响,改善Rife算法在低信噪比环境下和部分频偏值的估计精度恶化的问题。且通过阈值判定,只将阈值外的点进行反馈重估,增加了精度的同时减少不必要的运算量。
以上所述仅是本发明的优选实施方式,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种基于FFT的GMSK信号频偏估计算法,其特征在于,包括步骤如下:
在解调时,先去除GMSK调制信号的相位信息,再利用反馈修正Rife算法进行频偏估计;利用反馈修正Rife算法进行频偏估计的步骤如下:
利用FFT运算得到最大谱线与次谱线,计算出频移因子;设定判断偏移量是否处于中心区域的门限值,当根据门限值判定位于中心区域时,直接利用得出的频移因子计算频率偏移量;否则,利用频谱细化求得最大谱线左右两侧相邻谱线中心位置频率的幅值作为修正方向的判据,将待测信号进行频谱搬移,对频移因子进行修正之后再计算频率偏移量。
2.根据权利要求1所述基于FFT的GMSK信号频偏估计算法,其特征在于,所述去除GMSK调制信号的相位信息后的信号,通过以下式实现:
x(t)=rI(t)FI(t)+rQ(t)FQ(t),F(t)=ejθ(t),
化简后为x(t)=cos(2πfot+2πfdt)+w(t),w(t)=FI(t)nI(t)+FQ(t)nQ(t),
其中rI(t)和rQ(t)分别表示接收到的信号r(t)的同相分量和正交分量,FI(t)和FQ(t)分别代表GMSK调制信号F(t)的同相分量和正交分量,f0和fd分别代表中间频率和多普勒频率偏移,θ是GMSK信号调制后的相位信息φ(t,α),n(t)是高斯白噪声,θ(t)是承载GMSK信号的相位信息。
3.根据权利要求2所述基于FFT的GMSK信号频偏估计算法,其特征在于,所述利用FFT运算得到最大谱线与次谱线,计算出频移因子的步骤如下:
对去除相位信息后的信号x(t)在[0,T]内进行周期为Ts的采样得到长度为N的序列x(n),x(n)=cos(2πf0nTs+2fdnTs)+w(n),0≤n≤N-1,其中n∈[0,1,…,N],w(n)为实部与虚部相互独立的、方差为2σ2的零均值复高斯白噪声;
对x(n)做N点FFT运算得到X(k),
其中,k=0,1,…N/2-1,W(k)是w(n)的离散傅立叶变换;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910934160.6A CN110691051A (zh) | 2019-09-29 | 2019-09-29 | 一种基于fft的gmsk信号频偏估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910934160.6A CN110691051A (zh) | 2019-09-29 | 2019-09-29 | 一种基于fft的gmsk信号频偏估计算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110691051A true CN110691051A (zh) | 2020-01-14 |
Family
ID=69110940
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910934160.6A Pending CN110691051A (zh) | 2019-09-29 | 2019-09-29 | 一种基于fft的gmsk信号频偏估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110691051A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458563A (zh) * | 2020-03-05 | 2020-07-28 | 熊军 | 一种双谱线幅度测量信号频偏的方法及装置 |
CN112014811A (zh) * | 2020-08-18 | 2020-12-01 | 南京航空航天大学 | 一种雷达载波频率的精细估计方法 |
CN114500188A (zh) * | 2021-12-24 | 2022-05-13 | 杭州电子科技大学 | 一种船舶自动识别***的频偏估计方法 |
CN114760178A (zh) * | 2022-04-08 | 2022-07-15 | 中国北方车辆研究所 | 一种基于多径能量的频偏估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102647382A (zh) * | 2011-02-16 | 2012-08-22 | 中兴通讯股份有限公司 | 整数倍频偏估计方法及装置 |
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
CN108206799A (zh) * | 2016-12-16 | 2018-06-26 | 中国航天科工飞航技术研究院 | 一种用于突发通信的载波频偏估计方法 |
-
2019
- 2019-09-29 CN CN201910934160.6A patent/CN110691051A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102647382A (zh) * | 2011-02-16 | 2012-08-22 | 中兴通讯股份有限公司 | 整数倍频偏估计方法及装置 |
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
CN108206799A (zh) * | 2016-12-16 | 2018-06-26 | 中国航天科工飞航技术研究院 | 一种用于突发通信的载波频偏估计方法 |
Non-Patent Citations (2)
Title |
---|
苑小华、罗武忠、罗来源: "MPSK信号载波频率盲估计", 《电子与信息学报》 * |
谢胜,陈航,于平: "基于FFT 并二次修正的Rife 频率估计算法", 《探测与控制学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458563A (zh) * | 2020-03-05 | 2020-07-28 | 熊军 | 一种双谱线幅度测量信号频偏的方法及装置 |
CN112014811A (zh) * | 2020-08-18 | 2020-12-01 | 南京航空航天大学 | 一种雷达载波频率的精细估计方法 |
CN112014811B (zh) * | 2020-08-18 | 2023-11-07 | 南京航空航天大学 | 一种雷达载波频率的精细估计方法 |
CN114500188A (zh) * | 2021-12-24 | 2022-05-13 | 杭州电子科技大学 | 一种船舶自动识别***的频偏估计方法 |
CN114500188B (zh) * | 2021-12-24 | 2024-01-26 | 杭州电子科技大学 | 一种船舶自动识别***的频偏估计方法 |
CN114760178A (zh) * | 2022-04-08 | 2022-07-15 | 中国北方车辆研究所 | 一种基于多径能量的频偏估计方法 |
CN114760178B (zh) * | 2022-04-08 | 2024-03-15 | 中国北方车辆研究所 | 一种基于多径能量的频偏估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110691051A (zh) | 一种基于fft的gmsk信号频偏估计算法 | |
CN108512791B (zh) | 基于定时频偏补偿的星载ais解调方法 | |
CN103095615B (zh) | 一种相位噪声估计、补偿的方法及装置 | |
CN104022981B (zh) | 一种正交幅度调制信号的盲载波频偏估计方法 | |
CN107342960B (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
CN112511477A (zh) | 一种基于星座图和深度学习的混合卫星通信调制识别方法及*** | |
EP3621259B1 (en) | Method and device for fsk/gfsk demodulation | |
US5093848A (en) | Method of controlling the frequency of a coherent radio receiver and apparatus for carrying out the method | |
CN109379310B (zh) | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 | |
CN103428153B (zh) | 卫星移动通信中gmsk信号接收方法 | |
CN103023831B (zh) | 一种适用于突发波形的载波频偏估计方法 | |
US9722845B2 (en) | Bluetooth low energy frequency offset and modulation index estimation | |
CN111935046B (zh) | 一种低复杂度的频移键控信号符号率估计方法 | |
CN113726716B (zh) | 一种基于判决反馈的载波相位误差鉴别方法 | |
CN104363194A (zh) | 基于波形变换的psk调制识别方法 | |
US20160323128A1 (en) | Circuits and methods for frequency offset estimation in fsk communications | |
CN104378317B (zh) | 一种应用于卫星电话的突发信号频偏校正方法 | |
JPH06205062A (ja) | 遅延検波回路 | |
CN109756435B (zh) | 一种对信号的频偏估计方法 | |
CN106534010B (zh) | 一种oqpsk信号解调的相位偏差估计与补偿方法 | |
WO2023078117A1 (zh) | Dpsk信号的解调方法、装置、设备及存储介质 | |
CN112671684B (zh) | 一种短时突发bpsk信号的自适应解调方法 | |
CN115118564A (zh) | 一种载波频率偏差估计方法及装置 | |
CN115632923A (zh) | 一种基于oqpsk的无人机与卫星超宽带通信方法及相关设备 | |
Xie et al. | An improved central frequency estimation method for frequency-hopping signal |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200114 |
|
WD01 | Invention patent application deemed withdrawn after publication |