CN114363133A - 一种频偏检测方法、系统、存储介质和电子设备 - Google Patents

一种频偏检测方法、系统、存储介质和电子设备 Download PDF

Info

Publication number
CN114363133A
CN114363133A CN202210067803.3A CN202210067803A CN114363133A CN 114363133 A CN114363133 A CN 114363133A CN 202210067803 A CN202210067803 A CN 202210067803A CN 114363133 A CN114363133 A CN 114363133A
Authority
CN
China
Prior art keywords
value
frequency
signal
frequency offset
detected
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.)
Granted
Application number
CN202210067803.3A
Other languages
English (en)
Other versions
CN114363133B (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.)
Wuhan Mengxin Technology Co ltd
Original Assignee
Wuhan Mengxin Technology 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 Wuhan Mengxin Technology Co ltd filed Critical Wuhan Mengxin Technology Co ltd
Priority to CN202210067803.3A priority Critical patent/CN114363133B/zh
Publication of CN114363133A publication Critical patent/CN114363133A/zh
Application granted granted Critical
Publication of CN114363133B publication Critical patent/CN114363133B/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/0014Carrier regulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2649Demodulators
    • H04L27/265Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • H04L27/2659Coarse or integer frequency offset determination and synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • H04L27/266Fine or fractional frequency offset determination and synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end
    • H04L2027/0026Correction of carrier offset
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0044Control loops for carrier regulation
    • H04L2027/0063Elements of loops
    • H04L2027/0065Frequency error detectors
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明涉及一种频偏检测方法、系统、存储介质和电子设备,所述方法包括:根据待检测信号的初始频偏值对待检测信号进行预纠偏,得到待检测信号的残留频偏范围,根据预设频偏搜索步长和预设时域数据长度,在残留频偏范围内搜索第一最大频点模值和第二最大频点模值,并根据第一最大频点模值的索引值和第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值。本发明的方法实现简单,快速高效,通过GNSS芯片的同步和定位过程,先确定一部分频偏值,并利用该频偏值进行预纠偏,然后仅需要在残留频偏范围内,进行DFT扫描,确认残留频偏值。不同波特率下,通过设置合理的预设参数就能保证LBAND信号的检测性能。

Description

一种频偏检测方法、系统、存储介质和电子设备
技术领域
本发明涉及卫星通信领域技术领域,尤其涉及一种频偏检测方法、系统、存储介质和电子设备。
背景技术
随着通信业务由海上用户向陆地移动和航空用户扩展,1994年改名为国际移动卫星通信组织(International Mobile Satellite Organization,IMSO)但其缩写仍然是Inmarsat。Inmarsat系统能覆盖地球南北纬76度之内的所有区域,能够提供多种通信服务,承担了目前海上绝大部分的通信业务。
在卫星移动通信系统中,由于卫星与移动终端之间的相对运动,会导致多普勒效应。因此必须分析在卫星可视范围内移动终端的多普勒频偏和变化规律,从而补偿多普勒效应,提高通信质量。同时基于成本考虑,绝大部分接收机自身的晶振精度不够高,也会产生时钟漂移频率偏差的问题。例如石英晶体的频率偏差率一般在0.25ppm,以L1频率的振荡器而言,0.25ppm相当于约400Hz的频率误差。该频率误差在基带信号处理时也必须要被考虑进来。
发明内容
为解决上述技术问题,本发明提供了一种频偏检测方法、系统、存储介质和电子设备。
本发明的一种频偏检测方法的技术方案如下:
S1、根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
S2、根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
S3、根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
本发明的一种频偏检测方法的有益效果如下:
本发明的方法根据待检测信号的初始频偏值对待检测信号进行预纠偏,得到待检测信号的残留频偏范围,根据预设频偏搜索步长和预设时域数据长度,在残留频偏范围内搜索第一最大频点模值和第二最大频点模值,并根据第一最大频点模值的索引值和第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值。因此,本发明的方法实现简单,快速高效,通过GNSS芯片的同步和定位过程,先确定一部分频偏值,并利用该频偏值进行预纠偏,然后仅需要在残留频偏范围内,进行DFT扫描,确认残留频偏值。不同波特率下,通过设置合理的预设参数就能保证LBAND信号的检测性能。
在上述方案的基础上,本发明的一种频偏检测方法还可以做如下改进。
进一步,在所述S1之前,还包括:
S01、通过GNSS接收机获取所述待检测信号的初始频偏值;
所述S1具体包括:
S11、根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure BDA0003480864190000021
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
S12、将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
S13、根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure BDA0003480864190000031
n=0,1,2,...,Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引,j为待检测信号的系数;
S14、根据所述零中频基带值确定所述待检测信号的残留频偏范围。
进一步,在所述S2之前,还包括:根据预设函数公式确定所述预设频偏搜索步长的取值范围,并将所述频偏搜索步长的取值范围内的合理取值确定为所述预设频偏搜索步长,其中,所述预设函数公式为:
Figure BDA0003480864190000032
Tcoh为相干积分长度,fe为频率误差,sinc(feTcoh)为所述频偏搜索步长的取值范围,用于获取所述频偏搜索步长的取值范围内的合理取值。
进一步,所述S2具体包括:
S21、获取待检测信号的符号速率Symbol_Rate;
S22、根据降采样公式,对所述零中频基带值和所述符号速率进行处理,得到降采样值,其中,所述降采样公式为:
Figure BDA0003480864190000033
M=Fs/Fw,Fs为所述采样频率,Fw=8×Symbol_Rate,k为所述降采样公式中的降采样值的离散数据点的索引,i为所述降采样公式中的零中频基带值的离散数据点的索引;
S23、根据低通滤波公式,对所述降采样值进行滤波处理,得到滤波值,其中,所述低通滤波公式为:
Figure BDA0003480864190000041
u=0,1,2…,c(n)u为第u次滤波,LLPF为低通滤波器阶数,hLPF(n)为滤波器系数,p(u)为第u次滤波后的滤波值;
S24、根据所述修正公式,对所述滤波值进行修正,得到修正值,其中,所述修正公式为:x(n)=(I(n)2-Q(n)2)+j*2*I(n)*Q(n)0≤n<NLen,I(n)和Q(n)分别为p(n)的实部和虚部,n表示离散的数据点的索引;
S25、基于离散傅里叶变换公式,在所述残留频偏范围内根据所述修正值、所述预设频偏搜索步长和所述预设时域数据长度得到所述待检测信号的第一最大频点模值,其中,所述离散傅里叶变换公式为:
Figure BDA0003480864190000042
X(m)为多个搜索频点,N1为第一搜索频点数,所述第一搜索频点数由预设频偏搜索步长和所述残留频偏范围确定,Δf为所述预设频偏搜索步长,NInt为预设的时域数据长度所对应的时域积分样点个数;
S26、在所述多个搜索频点中,获取所述第一最大频点模值的索引值idx1,idx1∈[1,...N1];
S27、在预设搜索条件下,根据所述第一最大频点模值得到所述第二最大频点模值的索引值idx2,idx2∈[1,...N2];
S28、基于预设转换公式,根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值得到第三最大频点模值的索引值,其中,预设转换公式为:
Figure BDA0003480864190000043
idx3为所述第三最大频点模值的索引值,N2为预设搜索频点个数。
进一步,所述根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,具体包括:
判断所述第一最大频点模值的索引值和所述第三最大频点模值的索引值之间的差值是否大于1,若是,则所述待检测信号中不存在所述残留频偏值,若否,则所述待检测信号中存在所述残留频偏值。
进一步,还包括:根据残留频偏公式得到所述残留频偏值,其中,所述残留频偏公式为:
Figure BDA0003480864190000051
本发明的一种频偏检测系统的技术方案如下:
包括:预处理模块、获取模块、判断模块;
所述预处理模块用于:根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
所述获取模块用于:根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
所述判断模块用于:根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
本发明的一种频偏检测系统的有益效果如下:
本发明的系统根据待检测信号的初始频偏值对待检测信号进行预纠偏,得到待检测信号的残留频偏范围,根据预设频偏搜索步长和预设时域数据长度,在残留频偏范围内搜索第一最大频点模值和第二最大频点模值,并根据第一最大频点模值的索引值和第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值。因此,本发明的系统实现简单,快速高效,通过GNSS芯片的同步和定位过程,先确定一部分频偏值,并利用该频偏值进行预纠偏,然后仅需要在残留频偏范围内,进行DFT扫描,确认残留频偏值。不同波特率下,通过设置合理的预设参数就能保证LBAND信号的检测性能。
在上述方案的基础上,本发明的一种频偏检测系统还可以做如下改进。
进一步,在所述预处理模块之前,还包括:接收模块;
所述接收模块用于:通过GNSS接收机获取所述待检测信号的初始频偏值;
所述预处理模块具体用于:
根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure BDA0003480864190000061
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure BDA0003480864190000062
n=0,1,2,...,Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引,j为待检测信号的系数;
根据所述零中频基带值确定所述待检测信号的残留频偏范围。
本发明的一种存储介质的技术方案如下:
存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如本发明的一种频偏检测方法的步骤。
本发明的一种电子设备的技术方案如下:
包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时,使所述计算机执行如本发明的一种频偏检测方法的步骤。
附图说明
图1为本发明实施例的一种频偏检测方法的流程示意图;
图2为本发明实施例的一种频偏检测方法的LBAND信号的高精度定位芯片的基带结构图;
图3为本发明实施例的一种频偏检测系统的结构示意图。
具体实施方式
如图1所示,本发明实施例的一种频偏检测方法,包括如下步骤:
S1、根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
其中,在卫星移动通信系统中,由于卫星与移动终端之间的相对运动,会导致多普勒效应。因此必须分析在卫星可视范围内移动终端的多普勒频偏和变化规律,从而补偿多普勒效应,提高通信质量。
其中,待检测信号为LBAND信号。本实施例的技术方案基于LBAND信号的高精度定位芯片的基带如图2所示。GNSS基带和LBAND信号基带在同一芯片内部,共用同一个晶振,因此完全可以先利用GNSS基带的同步和定位,确定终端钟漂和多普勒运行带来的频偏,对信号进行预纠偏。这样在LBAND信号检测阶段,需要搜索的残留频偏范围就会大大减小。同时基于成本考虑,绝大部分接收机自身的晶振精度不够高,也会产生时钟漂移频率偏差的问题。例如石英晶体的频率偏差率一般在0.25ppm,以L1频率的振荡器而言,0.25ppm相当于约400Hz的频率误差。该频率误差在基带信号处理时也必须要被考虑进来。
S2、根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
具体地,采样离散傅里叶变换(DFT)计算方法,在预先定义的频偏搜索步长和时域数据长度,做第一次DFT得到待检测信号的第一最大频点模值,其中,第一最大频点模值为第一次DFT后,得到的所有频点模值中的最大值(最大功率所对应的频偏值),然后再以第一最大频点模值为中心,在预设搜索频点个数的搜索频点范围内做第二次DFT,得到待检测信号的第二最大频点模值,其中,第二最大频点模值为第二次DFT后,得到的所有频点模值中的最大值(最大功率所对应的频偏值)。
S3、根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
具体地,根据两次最大值的位置差值确定LBAND信号的残留频偏,并确认LBAND信号是否存在。
较优地,在所述S1之前,还包括:
S01、通过GNSS接收机获取所述待检测信号的初始频偏值;
其中,GNSS接收机通过与GPS/BD/GAL等卫星实现同步和定位,得到初始频偏值。
所述S1具体包括:
S11、根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure BDA0003480864190000081
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
其中,RFFreq_new为32位寄存器的值,也就是补偿频率对应的32位的频率控制字。
S12、将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
其中,射频前端持续输出得到的d(0)、d(1)、d(2)...,最终得到基带信号值:d=[d(0),d(1),…]。
S13、根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure BDA0003480864190000091
n=0,1,2,...,Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引,j为待检测信号的系数;
S14、根据所述零中频基带值确定所述待检测信号的残留频偏范围。
其中,将32位的频率控制字传给LBAND基带后,经过预纠偏,信号的残留频偏能保证在±1KHz范围之内。用频率控制字生成频率为对应频率的载波,对LBAND信号进行下变频处理,即预纠偏。预纠偏后信号的残留频偏在±1kHz范围内。
较优地,在所述S2之前,还包括:根据预设函数公式确定所述预设频偏搜索步长的取值范围,并将所述频偏搜索步长的取值范围内的合理取值确定为所述预设频偏搜索步长,其中,所述预设函数公式为:
Figure BDA0003480864190000092
Tcoh为相干积分长度,fe为频率误差,sinc(feTcoh)为所述频偏搜索步长的取值范围,用于获取所述频偏搜索步长的取值范围内的合理取值。
其中,如表1中的内容所示,信号具有4种不同的波特率。时域上相同数量的离散点对应的相干积分长度Tcoh并不相同,因此在DFT计算之前,为了保证因频率差异引起的相干结果幅值衰减始终在3db范围之内,不同的波特率下应该采用不同的频偏搜索步长Δf。相干积分结果的幅值随频率误差Fe的变化曲线为函数|sinc(Fe*Tcoh)|,其中函数
Figure BDA0003480864190000101
sin为正弦函数,该曲线能指导频偏搜索步长的合理取值范围。
表1 LBAND信号各波特率下参数细节表
Figure BDA0003480864190000102
较优地,所述S2具体包括:
S21、获取待检测信号的符号速率Symbol_Rate;
S22、根据降采样公式,对所述零中频基带值和所述符号速率进行处理,得到降采样值,其中,所述降采样公式为:
Figure BDA0003480864190000103
M=Fs/Fw,Fs为所述采样频率,Fw=8×Symbol_Rate,k为所述降采样公式中的降采样值的离散数据点的索引,i为所述降采样公式中的零中频基带值的离散数据点的索引;
S23、根据低通滤波公式,对所述降采样值进行滤波处理,得到滤波值,其中,所述低通滤波公式为:
Figure BDA0003480864190000111
u=0,1,2…,c(n)u为第u次滤波,LLPF为低通滤波器阶数,hLPF(n)为滤波器系数,p(u)为第u次滤波后的滤波值;
其中,由于降采样前后采样率差距过大,会在降采样后的通带内引入混叠噪声。因此为了提高频偏估计的准确性,需要先过一级低通滤波器。滤波器的通带由符号速率加1KHz最大频偏确定。
S24、根据所述修正公式,对所述滤波值进行修正,得到修正值,其中,所述修正公式为:x(n)=(I(n)2-Q(n)2)+j*2*I(n)*Q(n)0≤n<NLen,I(n)和Q(n)分别为p(n)的实部和虚部,n表示离散的数据点的索引;
其中,需要说明的是,如果输入数据存在频偏fd,经过修正后,频偏变成fnew=2*fd。后续的频偏估计实际上都是对fnew进行估计,从而折算出fd
S25、基于离散傅里叶变换公式,在所述残留频偏范围内根据所述修正值、所述预设频偏搜索步长和所述预设时域数据长度得到所述待检测信号的第一最大频点模值,其中,所述离散傅里叶变换公式为:
Figure BDA0003480864190000112
X(m)为多个搜索频点,N1为第一搜索频点数,所述第一搜索频点数由预设频偏搜索步长和所述残留频偏范围确定,Δf为所述预设频偏搜索步长,NInt为预设的时域数据长度所对应的时域积分样点个数;
S26、在所述多个搜索频点中,获取所述第一最大频点模值的索引值idx1,idx1∈[1,...N1];
S27、在预设搜索条件下,根据所述第一最大频点模值得到所述第二最大频点模值的索引值idx2,idx2∈[1,...N2];
S28、基于预设转换公式,根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值得到第三最大频点模值的索引值,其中,预设转换公式为:
Figure BDA0003480864190000121
idx3为所述第三最大频点模值的索引值,N2为预设搜索频点个数。
较优地,所述根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,具体包括:
判断所述第一最大频点模值的索引值和所述第三最大频点模值的索引值之间的差值是否大于1,若是,则所述待检测信号中不存在所述残留频偏值,若否,则所述待检测信号中存在所述残留频偏值。
较优地,还包括:根据残留频偏公式得到所述残留频偏值,其中,所述残留频偏公式为:
Figure BDA0003480864190000122
其中,若存在残留频偏值,则最终频偏值为残留频偏值和初始频偏值之和。
需要说明的是,本实施例中的预设参数可以根据算法要求和场景,动态做各种各样的修改,本发明并不局限该参数的设定。
本实施例的方法根据待检测信号的初始频偏值对待检测信号进行预纠偏,得到待检测信号的残留频偏范围,根据预设频偏搜索步长和预设时域数据长度,在残留频偏范围内搜索第一最大频点模值和第二最大频点模值,并根据第一最大频点模值的位置和第二最大频点模值的位置之间的差值判断所述待检测信号中是否存在残留频偏值。因此,本实施例的方法实现简单,快速高效,通过GNSS芯片的同步和定位过程,先确定一部分频偏值,并利用该频偏值进行预纠偏,然后仅需要在残留频偏范围内,进行DFT扫描,确认残留频偏值。不同波特率下,通过设置合理的预设参数就能保证LBAND信号的检测性能。
如图3所示,本发明实施例的一种频偏检测系统200,包括:预处理模块210、获取模块220、判断模块230;
所述预处理模块210用于:根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
所述获取模块220用于:根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
所述判断模块230用于:根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
较优地,在所述预处理模块210之前,还包括:接收模块;
所述接收模块用于:通过GNSS接收机获取所述待检测信号的初始频偏值;
所述预处理模块210具体用于:
根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure BDA0003480864190000131
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure BDA0003480864190000141
n=0,1,2,...,Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引;
根据所述零中频基带值确定所述待检测信号的残留频偏范围。
本实施例的系统根据待检测信号的初始频偏值对待检测信号进行预纠偏,得到待检测信号的残留频偏范围,根据预设频偏搜索步长和预设时域数据长度,在残留频偏范围内搜索第一最大频点模值和第二最大频点模值,并根据第一最大频点模值的位置和第二最大频点模值的位置之间的差值判断所述待检测信号中是否存在残留频偏值。因此,本实施例的系统实现简单,快速高效,通过GNSS芯片的同步和定位过程,先确定一部分频偏值,并利用该频偏值进行预纠偏,然后仅需要在残留频偏范围内,进行DFT扫描,确认残留频偏值。不同波特率下,通过设置合理的预设参数就能保证LBAND信号的检测性能。
上述关于本实施例的一种频偏检测系统200中的各参数和各个模块实现相应功能的步骤,可参考上文中关于一种频偏检测方法的实施例中的各参数和步骤,在此不做赘述。
本发明实施例提供的一种存储介质,包括:存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如一种频偏检测方法的步骤,具体可参考上文中一种频偏检测方法的实施例中的各参数和步骤,在此不做赘述。
计算机存储介质例如:优盘、移动硬盘等。
本发明实施例提供的一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时,使所述计算机执行如一种频偏检测方法的步骤,具体可参考上文中一种频偏检测方法的实施例中的各参数和步骤,在此不做赘述。
所属技术领域的技术人员知道,本发明可以实现为方法、系统、存储介质和电子设备。
因此,本发明可以具体实现为以下形式,即:可以是完全的硬件、也可以是完全的软件(包括固件、驻留软件、微代码等),还可以是硬件和软件结合的形式,本文一般称为“电路”、“模块”或“系统”。此外,在一些实施例中,本发明还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可读的程序代码。可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM),只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (10)

1.一种频偏检测方法,其特征在于,包括:
S1、根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
S2、根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
S3、根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
2.根据权利要求1所述的一种频偏检测方法,其特征在于,在所述S1之前,还包括:
S01、通过GNSS接收机获取所述待检测信号的初始频偏值;
所述S1具体包括:
S11、根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure FDA0003480864180000011
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
S12、将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
S13、根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure FDA0003480864180000021
Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引,j为待检测信号的系数;
S14、根据所述零中频基带值确定所述待检测信号的残留频偏范围。
3.根据权利要求2所述的一种频偏检测方法,其特征在于,在所述S2之前,还包括:根据预设函数公式确定所述预设频偏搜索步长的取值范围,并将所述频偏搜索步长的取值范围内的合理取值确定为所述预设频偏搜索步长,其中,所述预设函数公式为:
Figure FDA0003480864180000022
Tcoh为相干积分长度,fe为频率误差,sinc(feTcoh)为所述频偏搜索步长的取值范围,用于获取所述频偏搜索步长的取值范围内的合理取值。
4.根据权利要求3所述的一种频偏检测方法,其特征在于,所述S2具体包括:
S21、获取待检测信号的符号速率Symbol_Rate;
S22、根据降采样公式,对所述零中频基带值和所述符号速率进行处理,得到降采样值,其中,所述降采样公式为:
Figure FDA0003480864180000023
M=Fs/Fw,Fs为所述采样频率,Fw=8×Symbol_Rate,k为所述降采样公式中的降采样值的离散数据点的索引,i为所述降采样公式中的零中频基带值的离散数据点的索引;
S23、根据低通滤波公式,对所述降采样值进行滤波处理,得到滤波值,其中,所述低通滤波公式为:
Figure FDA0003480864180000024
c(n)u为第u次滤波,LLPF为低通滤波器阶数,hLPF(n)为滤波器系数,p(u)为第u次滤波后的滤波值;
S24、根据所述修正公式,对所述滤波值进行修正,得到修正值,其中,所述修正公式为:x(n)=(I(n)2-Q(n)2)+j*2*I(n)*Q(n)0≤n<NLen,I(n)和Q(n)分别为p(n)的实部和虚部,n表示离散的数据点的索引;
S25、基于离散傅里叶变换公式,在所述残留频偏范围内根据所述修正值、所述预设频偏搜索步长和所述预设时域数据长度得到所述待检测信号的第一最大频点模值,其中,所述离散傅里叶变换公式为:
Figure FDA0003480864180000031
X(m)为多个搜索频点,N1为第一搜索频点数,所述第一搜索频点数由预设频偏搜索步长和所述残留频偏范围确定,Δf为所述预设频偏搜索步长,NInt为预设的时域数据长度所对应的时域积分样点个数;
S26、在所述多个搜索频点中,获取所述第一最大频点模值的索引值idx1,idx1∈[1,...N1];
S27、在预设搜索条件下,根据所述第一最大频点模值得到所述第二最大频点模值的索引值idx2,idx2∈[1,...N2];
S28、基于预设转换公式,根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值得到第三最大频点模值的索引值,其中,预设转换公式为:
Figure FDA0003480864180000032
idx3为所述第三最大频点模值的索引值,N2为预设搜索频点个数。
5.根据权利要求4所述的一种频偏检测方法,其特征在于,所述根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,具体包括:
判断所述第一最大频点模值的索引值和所述第三最大频点模值的索引值之间的差值是否大于1,若是,则所述待检测信号中不存在所述残留频偏值,若否,则所述待检测信号中存在所述残留频偏值。
6.根据权利要求5所述的一种频偏检测方法,其特征在于,还包括:根据残留频偏公式得到所述残留频偏值,其中,所述残留频偏公式为:
Figure FDA0003480864180000041
7.一种频偏检测系统,其特征在于,包括:预处理模块、获取模块、判断模块;
所述预处理模块用于:根据待检测信号的初始频偏值对所述待检测信号进行预纠偏,得到所述待检测信号的残留频偏范围;
所述获取模块用于:根据预设频偏搜索步长和预设时域数据长度,在所述残留频偏范围内搜索并得到所述待检测信号的第一最大频点模值,并根据预设搜索频点个数和所述第一最大频点模值得到所述待检测信号的第二最大频点模值;
所述判断模块用于:根据所述第一最大频点模值的索引值和所述第二最大频点模值的索引值之间的差值判断所述待检测信号中是否存在残留频偏值,若是,则根据所述初始频偏值和所述残留频偏值确定所述待检测信号的最终频偏值,若否,则将所述初始频偏值确定为所述待检测信号的最终频偏值。
8.根据权利要求7所述的一种频偏检测系统,其特征在于,在所述预处理模块之前,还包括:接收模块;
所述接收模块用于:通过GNSS接收机获取所述待检测信号的初始频偏值;
所述预处理模块具体用于:
根据频偏值转化公式,将所述初始频偏值转化为频率控制字,其中,所述频偏值转化公式为:
Figure FDA0003480864180000042
fd0为所述初始频偏值,RFFreq_new为所述频率控制字,RFfrequency为所述待检测信号的理论中频频率,Fs为所述待检测信号的采样频率;
将所述频率控制字发送至LBAND基带,并经所述LBAND基带的射频前端持续输出得到基带信号值,其中,所述基带信号值为d=[d(0),d(1),…];
根据下变频公式对所述基带信号值进行处理,得到零中频基带值,其中,所述下变频公式为:
Figure FDA0003480864180000051
Fs为所述待检测信号的采样频率,Fc为所述待检测信号的中频频率,n表示所述下变频公式中的离散数据点的索引,j为待检测信号的系数;
根据所述零中频基带值确定所述待检测信号的残留频偏范围。
9.一种存储介质,其特征在于,所述存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如权利要求1至6中任一项所述的一种频偏检测方法。
10.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时,使所述计算机执行如权利要求1至6中任一项所述的一种频偏检测方法。
CN202210067803.3A 2022-01-20 2022-01-20 一种频偏检测方法、系统、存储介质和电子设备 Active CN114363133B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210067803.3A CN114363133B (zh) 2022-01-20 2022-01-20 一种频偏检测方法、系统、存储介质和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210067803.3A CN114363133B (zh) 2022-01-20 2022-01-20 一种频偏检测方法、系统、存储介质和电子设备

Publications (2)

Publication Number Publication Date
CN114363133A true CN114363133A (zh) 2022-04-15
CN114363133B CN114363133B (zh) 2023-10-03

Family

ID=81091388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210067803.3A Active CN114363133B (zh) 2022-01-20 2022-01-20 一种频偏检测方法、系统、存储介质和电子设备

Country Status (1)

Country Link
CN (1) CN114363133B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117176529A (zh) * 2023-11-02 2023-12-05 北京融为科技有限公司 Fsk信号解调方法及装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020089770A1 (en) * 2001-01-10 2002-07-11 Keirn Zachary A. Channel quality monitor (CQM) for digital peak detection (DPD)
US20040240593A1 (en) * 2003-06-02 2004-12-02 Hsiao-Lan Su Frequency offset compensation estimation system and method for a wireless local area network
US20060230769A1 (en) * 2005-01-11 2006-10-19 Kabushiki Kaisha Toshiba Superconducting magnet device-monitoring system, method of monitoring superconducting magnet device and MRI device
CN1960235A (zh) * 2005-11-03 2007-05-09 北京信威通信技术股份有限公司 一种基于正交频分复用的混合自动请求重传的方法
CN102202021A (zh) * 2011-05-04 2011-09-28 福州瑞芯微电子有限公司 一种广播定位系统的粗略频偏估计方法
WO2014135785A2 (fr) * 2013-03-05 2014-09-12 Electricite De France Procédé de détection d'un défaut de court-circuit dans des bobinages d'un rotor d'une machine électrique tournante
WO2018054058A1 (zh) * 2016-09-23 2018-03-29 深圳市中兴微电子技术有限公司 一种光传输信道中频偏估计方法、装置及计算机存储介质
JP2018107763A (ja) * 2016-12-28 2018-07-05 ラピスセミコンダクタ株式会社 半導体装置及び振幅検出方法
CN109889460A (zh) * 2019-01-25 2019-06-14 武汉虹信通信技术有限责任公司 上行频偏跟踪补偿方法及装置
CN110855590A (zh) * 2019-11-18 2020-02-28 紫光展锐(重庆)科技有限公司 频偏补偿方法、系统、电子设备及计算机可读存储介质
WO2020052683A1 (zh) * 2018-09-13 2020-03-19 中兴通讯股份有限公司 频偏估计补偿方法、装置、通信设备及计算机可读存储介质
CN112367283A (zh) * 2020-10-30 2021-02-12 青岛鼎信通讯股份有限公司 一种用于gfsk系统的载波频偏估计方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020089770A1 (en) * 2001-01-10 2002-07-11 Keirn Zachary A. Channel quality monitor (CQM) for digital peak detection (DPD)
US20040240593A1 (en) * 2003-06-02 2004-12-02 Hsiao-Lan Su Frequency offset compensation estimation system and method for a wireless local area network
US20060230769A1 (en) * 2005-01-11 2006-10-19 Kabushiki Kaisha Toshiba Superconducting magnet device-monitoring system, method of monitoring superconducting magnet device and MRI device
CN1960235A (zh) * 2005-11-03 2007-05-09 北京信威通信技术股份有限公司 一种基于正交频分复用的混合自动请求重传的方法
CN102202021A (zh) * 2011-05-04 2011-09-28 福州瑞芯微电子有限公司 一种广播定位系统的粗略频偏估计方法
WO2014135785A2 (fr) * 2013-03-05 2014-09-12 Electricite De France Procédé de détection d'un défaut de court-circuit dans des bobinages d'un rotor d'une machine électrique tournante
WO2018054058A1 (zh) * 2016-09-23 2018-03-29 深圳市中兴微电子技术有限公司 一种光传输信道中频偏估计方法、装置及计算机存储介质
JP2018107763A (ja) * 2016-12-28 2018-07-05 ラピスセミコンダクタ株式会社 半導体装置及び振幅検出方法
WO2020052683A1 (zh) * 2018-09-13 2020-03-19 中兴通讯股份有限公司 频偏估计补偿方法、装置、通信设备及计算机可读存储介质
CN109889460A (zh) * 2019-01-25 2019-06-14 武汉虹信通信技术有限责任公司 上行频偏跟踪补偿方法及装置
CN110855590A (zh) * 2019-11-18 2020-02-28 紫光展锐(重庆)科技有限公司 频偏补偿方法、系统、电子设备及计算机可读存储介质
CN112367283A (zh) * 2020-10-30 2021-02-12 青岛鼎信通讯股份有限公司 一种用于gfsk系统的载波频偏估计方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A, KHANNA;R, KHURANA;A, KYRIACOU;R, DAVIES;D W, RAY: "Management of adrenocortical insufficiency with continuous subcutaneous hydrocortisone infusion: long-term experience in three patients.", ENDOCRINOLOGY, DIABETES & METABOLISM CASE REPORTS, vol. 2015 *
WARR, P. A.: "Amplitude Offset Estimation by Phase Comparison in Suppression Loops", IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, vol. 58, no. 7, XP011310071 *
杨勇;张冬玲;彭华;涂世龙;: "单通道高阶调制线性混合信号的频偏估计", 电子学报, no. 01 *
秦文兵;罗来源;: "CDMA2000移动台的SSOA估计技术研究", 电子学报, no. 04 *
陈晓挺;刘会杰;梁旭文;: "扩频系统中多普勒频偏快速估计算法的性能分析", 宇航学报, no. 02 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117176529A (zh) * 2023-11-02 2023-12-05 北京融为科技有限公司 Fsk信号解调方法及装置
CN117176529B (zh) * 2023-11-02 2024-04-02 北京融为科技有限公司 Fsk信号解调方法及装置

Also Published As

Publication number Publication date
CN114363133B (zh) 2023-10-03

Similar Documents

Publication Publication Date Title
CN102540221B (zh) 信号处理装置与信号处理方法
TWI358546B (en) Method for acquisition of gps signals and gps rece
CN110531393B (zh) 基于自适应卡尔曼滤波的弱信号跟踪方法
CN114363133B (zh) 一种频偏检测方法、系统、存储介质和电子设备
CN103873079A (zh) 接收装置
US8792591B1 (en) Systems and methods for I/Q imbalance correction and calibration of variable modulus signals
CN108027442B (zh) 具有定点西格玛罗滤波器的卫星导航接收机
JP5029190B2 (ja) 受信回路、電子機器及びノイズキャンセル方法
CN111510409B (zh) 利用bpsk数据估计天基机会信号多普勒的方法和系统
EP1847838B1 (en) Method and apparatus for frequency estimation
WO2020040889A1 (en) System-on-a-chip for reception of telemetry messages over a radio frequency channel
JP2006217601A (ja) 地理的局所化システムの位置決め信号を捕捉するための方法、地理的局所化システム用受信機及び上記の方法を実行するためのプログラム命令を含んでいるコンピュータデータ媒体
NZ201807A (en) Modem clock extraction:double band pass filter and square
JP3809432B2 (ja) 受信信号処理装置及び衛星測位システム
CN112014864B (zh) 一种基于最速上升法的高动态卫星跟踪方法
CN116636151A (zh) 用于卫星导航接收器的自适应窄频带和宽频带干扰抑制
CN110554413B (zh) 基于相位-幅值锁定环路的连续波干扰检测与辨识方法
CN109633713B (zh) 一种北斗导航弱信号接收的频率跟踪方法
KR100826376B1 (ko) Cdma 수신기의 매핑을 이용한 상관 방법 및 신호 처리방법
WO2020171735A1 (en) Method of energy signal detection
CN116505978B (zh) 一种蓝牙信号的处理方法、装置、电子设备及存储介质
GB2412738A (en) Enhanced sonde recognition
CN113721271B (zh) 一种北斗信号接收机接收信号快速并行捕获方法和系统
CN116046022B (zh) 光纤陀螺的仿真测试方法及系统
US20230058231A1 (en) Methods of and apparatus for digital filtering

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
GR01 Patent grant
GR01 Patent grant