CN110708267B - 频偏信息估计值确定方法 - Google Patents
频偏信息估计值确定方法 Download PDFInfo
- Publication number
- CN110708267B CN110708267B CN201910947297.5A CN201910947297A CN110708267B CN 110708267 B CN110708267 B CN 110708267B CN 201910947297 A CN201910947297 A CN 201910947297A CN 110708267 B CN110708267 B CN 110708267B
- Authority
- CN
- China
- Prior art keywords
- estimated
- frequency offset
- frequency
- carrier signal
- order
- 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
Links
Images
Abstract
本申请公开了一种频偏信息估计值确定方法。其中，该方法包括：采集包括待估计频偏信息的载波信号；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。本申请解决了现有技术中的粗估计算法分辨率较低的技术问题。
Description
技术领域
本申请涉及卫星通信领域，具体而言，涉及一种频偏信息估计值确定方法。
背景技术
在低轨卫星的高动态场景下，卫星具有一阶、二阶甚至更高阶的多普勒频偏变化率，这对卫星信号的捕获和跟踪造成了极大的困难。本申请以DVB-S2协议为例，图1为现有技术中的DVB-S2协议中同步流程示意图；图1中，接收信号经过基于星历的多普勒频偏预补偿12、定时同步14、帧同步16后进行载波同步18；图2为载波同步模块的结构示意图；载波同步模块22包括载波频率捕获模块24和载波相位跟踪模块26，载波频率捕获模块24包括粗估计模块28与细估计模块210。在载波频率捕获过程中，通常先采用“粗估计”，即对接收到的信号的频偏及其变化率进行粗略的估计。
在“粗估计”中，分辨率是制约算法的关键点。高分辨率能将频偏及其变化率缩小到更小的范围，是后面细估计的基础。现有的数据辅助类算法，如：Kay算法、Fitz算法、L&W算法、L&R算法、M&M算法和导频辅助检测的最大似然估计算法以及将M&M、L&R算法和细估计算法的级联算法等都视信号频偏为定值，不适用于频偏具体高阶变化率的高动态场景。基于时域延迟类的相关算法在低信噪比时算法性能较差。DCFT能很好的捕获频偏和变化率，但分辨率稍差。因此，现有的粗估计算法分辨率太低，对细估计的实施造成了影响。因而，需要设计出一种高分辨率的粗估计算法。
针对上述的问题，目前尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种频偏信息估计值确定方法，以至少解决现有技术中的粗估计算法分辨率较低的技术问题。
根据本申请实施例的一个方面，提供了一种频偏信息估计值确定方法，包括：采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
可选地，所述旋转因子的表达式为
可选地，所述载波信号的表达式为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率。
可选地，通过以下公式确定所述载波信号的频域信息：
其中，X(k，l)为所述频域信息。
可选地，通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
可选地，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果之后，所述方法还包括：基于所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果对所述待估计频偏以及所述待估计频偏的一阶变化率进行细估计。
根据本申请实施例的一个方面，提供了一种频偏信息估计值确定装置，包括：采集模块，用于采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；确定模块，用于根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
可选地，所述旋转因子的表达式为
其中：
k，1分别为对应所述频域信息的第一频率和所述第一频率的一阶变化率的点，MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数，nT所述载波信号的采样点对应的时间信息，N为所述载波信号的采样点数，T为所述载波信号的采样周期，(f0)max为所述载波信号的预设最大频偏，(f1)max为所述载波信号的频偏的预设最大一阶变化率，所述A为所述第一参数，所述B为所述第二参数。
可选地，所述载波信号的表达式为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率。
可选地，所述确定模块通过以下公式确定所述载波信号的频域信息：
其中，X(k，l)为所述频域信息。
可选地，所述确定模块通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
其中，所述
为所述待估计频偏的估计结果，所述
为所述待估计频偏的一阶变化率的估计结果，所述km为所述目标频率对应的点，所述lm为所述目标一阶变化率对应的点。
根据本申请实施例的一个方面，提供了一种存储介质，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行上述的频偏信息估计值确定方法。
根据本申请实施例的一个方面，提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行上述的频偏信息估计值确定方法。
在本申请实施例中，采用采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果的方式，通过引入第一参数与第二参数，达到了提高粗估计算法分辨率的技术效果，进而解决了现有技术中的粗估计算法分辨率较低的技术问题。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1为现有技术中的DVB-S2协议中同步流程示意图；
图2为载波同步模块的结构示意图；
图3是根据本申请实施例的一种可选的频偏信息估计值确定方法的流程示意图；
图4是根据本申请实施例的一种可选的频偏信息估计值确定装置的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
根据本申请实施例，提供了一种频偏信息估计值确定方法实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
在频率捕获过程中，主要是为了给出粗略的频偏值，可假设相位θ＝0。在DVB-S2这类系统中，辅助数据较短。载波捕获算法在设计时，可以忽略二级及其以上的高阶变化率频偏对信号的影响，故信号为：
式中，w(n)为信号相关高阶变化率频偏的信息，可以忽略不计，N是DFT变换区间长度，即载波信号的采样点数，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率，nT所述载波信号的采样点对应的时间信息，N为所述载波信号的第三采样点数，T为所述载波信号的采样周期，f0为nT时刻对应的待估计频偏，f1为nT时刻所述待估计频偏的一阶变化率，n为大于等于0，小于等于N-1的整数。
载波频率捕获算法可运用DCFT算法。DCFT的数学表达式如下：
式中，
为旋转因子，可将其简写成
k、l分别为DCFT运算中频率分量和一阶变化率分量对应的点，即k，1分别为所述频域信息中的第一频率和所述第一频率的一阶变化率的点数，X(k，l)为引入旋转因子以后，载波信号的频域信息得：
其中，X(k，1)为频域信息，当
且
时，|X(k，l)|最大，此时对应的k，l即为km和1m值，进而求出频偏f0及一阶变化率f1估计值
在DCFT算法中，因为k、1只能取整数值，直接采用此方法，f0分辨率为
频偏一阶变化率f1分辨率为
在实际的硬件实现中，FFT点数常取2n，N为1024。若低轨卫星通信系统符号速率为25Msps，则f0的分辨率为
Hz，f1的分辨率为
针对DCFT算法估计的精度不佳的问题，MDCFT(Modified Discrete ChirpFourier Transform)提出将旋转因子
替换为
MDCFT的定义如下：
式中，M表示频率一阶变化率点数，N表示频率点数，
同理，当
且
时，X(k，1)模值最大。通过搜索最大的模值，找到对应的km值和1m值，然后计算出f0和f1的估计值。
和DCFT算法中条件一样，N取1024，T取0.04us，f0的分辨率仍然为24414Hz。而f1分辨率和(f1)max、M值有关。若(f1)max取值5KHz/s，M取值500，则f1分辨率为20Hz/s。不难看出，MDCFT中旋转因子的改变改善了f1的分辨率，f0的分辨率仍然不够的缺点。
图3是根据本申请实施例的频偏信息估计值确定方法的流程示意图，如图3所示，该方法至少包括如下步骤：
步骤S302，采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；
在本申请的一些可选的实施例中，载波信号的表达式可以为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率，N为所述载波信号的第三采样点数，T为所述载波信号的采样周期，f0为nT时刻对应的待估计频偏，f1为nT时刻所述待估计频偏的一阶变化率，n为大于等于0，小于等于N-1的整数。
步骤S304，根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；
可选地，可通过以下公式确定所述载波信号的频域信息：
其中，X(k，l)为所述频域信息，可选地，
为所述旋转因子的表达式，其中：
其中：
k，1分别为对应所述频域信息中的第一频率和所述第一频率的一阶变化率的点数，MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数，nT所述载波信号的采样点对应的时间信息，N为所述载波信号的采样点数，T为所述载波信号的采样周期，(f0)max为所述载波信号的预设最大频偏，(f1)max为所述载波信号的频偏的预设最大一阶变化率，所述A为所述第一参数，所述B为所述第二参数。其中，k，1可以分别为计算中的对应所述频域信息的第一频率和所述第一频率的一阶变化率的点。
步骤S306，确定所述频域信息对应频域的幅值的最大值；
步骤S308，确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；
步骤S310，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
可选地，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果可以通过以下方式进行实现：
基于所述幅值最大值时，目标频率对应的点与所述目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
在本申请的一些可选的实施例中，引入旋转因子后频域信息的表达式为：
式中，
MB与MA分别表示频率和一阶变化率的点数，即MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数。
其中，当
且
时，|X(k，l)|最大，通过搜索x(k，1)最大模值，即可得到km和1m，从而计算法f0和f1的估计值
由上述可知，
分辨率和(f0)max、MB的取值有关，
分辨率和(f1)max、MA的取值有关。若以DVB-S2协议体制为例，频偏最大为符号速率的20％，本申请取符号速率25Msps，故(f0)max＝5MHz，(f1)max＝5kHz/s MA＝MB＝2048，则
分辨率为均为4882Hz，
分别按率为49Hz/s，分辨率高于未引入新的旋转因子前的算法。所述待估计频偏的估计结果为
与所述待估计频偏的一阶变化率的估计结果为
本申请粗估计算法如下：对信号引入新的旋转因子
后做DCFT运算，得到x(k，1)。其中
搜索x(k，1)最大模值，得到km和1m；
由
求出频偏和一阶变化率估计值
然后由
对信号进行补偿。
可选地，通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
可选地，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果之后，所述方法还包括：基于所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果对所述待估计频偏以及所述待估计频偏的一阶变化率进行细估计。
可选地，粗估计以后更有利于细估计的实施，将频偏估计到一定范围后再对信号进行“细估计”，跟踪信号，实现载波的精准捕获，获取捕获结果后，进行载波相位跟踪，以完成载波同步。
可选地，细估计以后，则可根据估计结果对频率进行补偿。
在本申请实施例中，采用采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果的方式，通过引入第一参数与第二参数，达到了提高粗估计算法分辨率的技术效果，进而解决了现有技术中的粗估计算法分辨率较低的技术问题。
根据本申请实施例，还提供了一种用于实施上述频偏信息估计值确定方法的频偏信息估计值确定装置，如图4所示，该装置包括：采集模块42、确定模块44；其中：
采集模块42，用于采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；
确定模块44，用于根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
可选地，所述旋转因子的表达式为
其中：
k，1分别为对应所述频域信息中的第一频率和所述第一频率的一阶变化率的点，MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数，nT所述载波信号的采样点对应的时间信息，N为所述载波信号的采样点数，T为所述载波信号的采样周期，(f0)max为所述载波信号的预设最大频偏，(f1)max为所述载波信号的频偏的预设最大一阶变化率。所述A为所述第一参数，所述B为所述第二参数。
可选地，所述载波信号的表达式为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率。
可选地，所述确定模块44通过以下公式确定所述载波信号的频域信息：
其中，X(k，l)为所述频域信息。
可选地，所述确定模块44通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
其中，所述
为所述待估计频偏的估计结果，所述
为所述待估计频偏的一阶变化率的估计结果，所述km为所述目标频率对应的点，所述lm为所述目标一阶变化率对应的点。
可选地，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果之后，所述装置还用于：基于所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果对所述待估计频偏以及所述待估计频偏的一阶变化率进行细估计。
需要说明的是，图4所示实施例的优选实施方式可以参见图3所示实施例的相关描述，此处不再赘述。
根据本申请实施例的另一个方面，还提供了一种存储介质，存储介质包括存储的程序，可选地，在本实施例中，存储介质被设置为存储用于执行以下步骤的程序代码：采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
根据本申请实施例的另一个方面，还提供了一种处理器，处理器用于运行程序，其中，程序运行时可以执行应用程序的频偏信息估计值确定方法中以下步骤的程序代码：采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；根据预设旋转因子与所述载波信号确定所述载波信号的频域信息，所述旋转因子为基于第一参数和第二参数共同确定的，所述第一参数为基于对应所述待估计频偏的一阶变化率的第一测试点数，以及所述载波信号的频偏的预设最大一阶变化率确定的，所述第二参数为基于对应所述待估计频偏的第二测试点数，以及所述载波信号的预设最大频偏确定的；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果。
上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。
在本申请的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上仅是本申请的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本申请原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本申请的保护范围。
Claims (11)
1.一种频偏信息估计值确定方法，其特征在于，包括：
采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；
根据预设旋转因子与所述载波信号确定所述载波信号的频域信息；
确定所述频域信息对应频域的幅值的最大值；
确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；
基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果；
其中，所述旋转因子的表达式为
式中，
k，l分别为对应所述频域信息的第一频率和所述第一频率的一阶变化率的点，MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数，N为所述载波信号的采样点数，n表示N个采样点中的第n个采样点，T为所述载波信号的采样周期，(f0)max为所述载波信号的预设最大频偏，(f1)max为所述载波信号的频偏的预设最大一阶变化率。
2.根据权利要求1所述的方法，其特征在于，所述载波信号的表达式为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率，nT为所述载波信号的第n个采样点对应的时间信息。
3.根据权利要求2所述的方法，其特征在于，通过以下公式确定所述载波信号的频域信息：
其中，X(k,l)为所述频域信息。
4.根据权利要求3所述的方法，其特征在于，通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
其中，所述
为所述待估计频偏的估计结果，所述
为所述待估计频偏的一阶变化率的估计结果，所述km为所述目标频率对应的点，所述lm为所述目标一阶变化率对应的点。
5.根据权利要求4所述的方法，其特征在于，基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果之后，所述方法还包括：
基于所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果对所述待估计频偏以及所述待估计频偏的一阶变化率进行细估计。
6.一种频偏信息估计值确定装置，其特征在于，包括：
采集模块，用于采集包括待估计频偏信息的载波信号，所述待估计频偏信息包括：待估计频偏与所述待估计频偏的一阶变化率；
确定模块，用于根据预设旋转因子与所述载波信号确定所述载波信号的频域信息；确定所述频域信息对应频域的幅值的最大值；确定所述最大值对应的目标频率对应的点和目标一阶变化率对应的点；基于所述目标频率对应的点和目标一阶变化率对应的点确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果；其中，所述旋转因子的表达式为
式中，
k，l分别为对应所述频域信息的第一频率和所述第一频率的一阶变化率的点，MA为对应所述待估计频偏的一阶变化率的第一测试点数，MB为对应所述待估计频偏的第二测试点数，N为所述载波信号的采样点数，n表示N个采样点中的第n个采样点，T为所述载波信号的采样周期，(f0)max为所述载波信号的预设最大频偏，(f1)max为所述载波信号的频偏的预设最大一阶变化率。
7.根据权利要求6所述的装置，其特征在于，所述载波信号的表达式为：
其中，x(n)为所述载波信号，f0为所述待估计频偏，f1为所述待估计频偏的一阶变化率，nT为所述载波信号的第n个采样点对应的时间信息。
8.根据权利要求7所述的装置，其特征在于，所述确定模块通过以下公式确定所述载波信号的频域信息：
其中，X(k,l)为所述频域信息。
9.根据权利要求8所述的装置，其特征在于，所述确定模块通过以下公式确定所述待估计频偏的估计结果与所述待估计频偏的一阶变化率的估计结果：
其中，所述
为所述待估计频偏的估计结果，所述
为所述待估计频偏的一阶变化率的估计结果，所述km为所述目标频率对应的点，所述lm为所述目标一阶变化率对应的点。
10.一种存储介质，其特征在于，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行权利要求1至5中任意一项所述的频偏信息估计值确定方法。
11.一种处理器，其特征在于，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1至5中任意一项所述的频偏信息估计值确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910947297.5A CN110708267B (zh) | 2019-09-30 | 2019-09-30 | 频偏信息估计值确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910947297.5A CN110708267B (zh) | 2019-09-30 | 2019-09-30 | 频偏信息估计值确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110708267A CN110708267A (zh) | 2020-01-17 |
| CN110708267B true CN110708267B (zh) | 2022-06-14 |
Family
ID=69196639
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910947297.5A Active CN110708267B (zh) | 2019-09-30 | 2019-09-30 | 频偏信息估计值确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110708267B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112911698B (zh) * | 2021-01-13 | 2023-04-11 | 北京中科晶上科技股份有限公司 | 通信系统中的定时同步方法及装置 |
| CN114615122B (zh) * | 2022-01-30 | 2023-11-03 | 北京邮电大学 | 通信信号的频偏确定方法及装置 |
| CN114978828B (zh) * | 2022-05-07 | 2023-04-28 | 广州海格通信集团股份有限公司 | 频偏估计方法、装置、电子设备及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6389040B1 (en) * | 1998-09-11 | 2002-05-14 | Lucent Technologies Inc. | Apparatus and method for generating a frequency offset estimate for communication systems having frequency selecting fading channels |
| US8407020B1 (en) * | 2010-09-28 | 2013-03-26 | The United States Of America As Represented By The Secretary Of The Navy | Fast method to search for linear frequency-modulated signals |
| CN104502898A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 将修正rft和mdcft相结合的机动目标参数估计方法 |
| CN109781153A (zh) * | 2018-12-17 | 2019-05-21 | 北京理工大学 | 物理参数估计方法、装置和电子设备 |
| CN110133598A (zh) * | 2019-05-09 | 2019-08-16 | 西安电子科技大学 | 基于FrFT的线性调频信号参数快速估计方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI487334B (zh) * | 2012-12-26 | 2015-06-01 | Mstar Semiconductor Inc | 載波頻率偏移補償裝置及方法 |
-
2019
- 2019-09-30 CN CN201910947297.5A patent/CN110708267B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6389040B1 (en) * | 1998-09-11 | 2002-05-14 | Lucent Technologies Inc. | Apparatus and method for generating a frequency offset estimate for communication systems having frequency selecting fading channels |
| US8407020B1 (en) * | 2010-09-28 | 2013-03-26 | The United States Of America As Represented By The Secretary Of The Navy | Fast method to search for linear frequency-modulated signals |
| CN104502898A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 将修正rft和mdcft相结合的机动目标参数估计方法 |
| CN109781153A (zh) * | 2018-12-17 | 2019-05-21 | 北京理工大学 | 物理参数估计方法、装置和电子设备 |
| CN110133598A (zh) * | 2019-05-09 | 2019-08-16 | 西安电子科技大学 | 基于FrFT的线性调频信号参数快速估计方法 |
Non-Patent Citations (4)
| Title |
|---|
| A modified discrete chirp-Fourier transform scheme;Pingyi Fan等;《WCC 2000 - ICSP 2000. 2000 5th International Conference on Signal Processing Proceedings. 16th World Computer Congress 2000》;20020806;全文 * |
| Doppler frequency parameters estimation for SAR imaging using a modified discrete chirp-Fourier transform;Li JunXian等;《IEEE International Conference Mechatronics and Automation, 2005》;20060508;全文 * |
| 一种改进MDCFT分辨率的方法;袁德生等;《科学技术与工程》;20140831;全文 * |
| 低轨卫星高动态通信链路同步技术研究;崔诵祺;《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》;20160215;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110708267A (zh) | 2020-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110708267B (zh) | 频偏信息估计值确定方法 | |
| EP1507378A1 (en) | Frame and frequency synchronization for OFDM | |
| CN101729461B (zh) | 消除单频干扰及多频干扰的系统及方法 | |
| KR102341875B1 (ko) | 송신기 및 수신기와 해당 방법들 | |
| CN111865865B (zh) | 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法 | |
| CN111123312A (zh) | 一种高动态低信噪比环境下的卫星弱信号捕获方法及系统 | |
| CN108802772B (zh) | Gnss卫星的捕获方法及装置、存储介质、终端 | |
| CN116366092A (zh) | 一种多普勒捕获方法、装置和存储介质 | |
| EP1079579A1 (en) | OFDM frame synchronisation | |
| CN110018500B (zh) | 一种基于圆周移位的北斗卫星信号捕获方法 | |
| CN117214928A (zh) | 一种伪码跟踪方法、装置、终端设备及存储介质 | |
| CN115296698A (zh) | 高动态卫星通信系统信号捕获方法 | |
| Svatoň et al. | Pre-and post-correlation method for acquisition of new GNSS signals with secondary code | |
| CN110445740A (zh) | 基于重复序列的频偏估计方法及系统 | |
| EP1442553B1 (en) | Apparatus and method for using training sequences to estimate timing error in a digital signal receiver | |
| CN114244674B (zh) | 一种超宽带基带接收机频偏估计方法及装置 | |
| CN113671547B (zh) | 一种改进的高动态捕获方法、装置、设备及存储介质 | |
| CN115643139A (zh) | 一种抗频偏的突发信号检测系统及检测方法 | |
| CN113315734A (zh) | 用于卫星的载波同步方法、装置、电子设备及可读介质 | |
| CN109302364A (zh) | 一种用于NB_IoT系统的频偏估计方法 | |
| CN113189542A (zh) | 窄带跳频信号tdoa定位方法、装置、系统 | |
| CN112764063A (zh) | 一种实现捕获处理的方法及接收机 | |
| US20060007571A1 (en) | Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive | |
| CN112987041B (zh) | 一种弱信号下的位同步方法、装置及计算机存储介质 | |
| CN115209135B (zh) | 一种电磁泄漏视频检测及还原系统 |
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 |