CN114244673B - 一种基于数据辅助的频率跟踪方法、系统及存储介质 - Google Patents
一种基于数据辅助的频率跟踪方法、系统及存储介质 Download PDFInfo
- Publication number
- CN114244673B CN114244673B CN202111565220.5A CN202111565220A CN114244673B CN 114244673 B CN114244673 B CN 114244673B CN 202111565220 A CN202111565220 A CN 202111565220A CN 114244673 B CN114244673 B CN 114244673B
- Authority
- CN
- China
- Prior art keywords
- data
- value
- loop
- phase
- processing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012545 processing Methods 0.000 claims abstract description 65
- 238000012937 correction Methods 0.000 claims abstract description 30
- 238000000605 extraction Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 42
- 230000007480 spreading Effects 0.000 claims description 16
- 230000001360 synchronised effect Effects 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000004891 communication Methods 0.000 abstract description 4
- 230000008878 coupling Effects 0.000 abstract description 2
- 238000010168 coupling process Methods 0.000 abstract description 2
- 238000005859 coupling reaction Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000001186 cumulative effect Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
- H04L27/2659—Coarse or integer frequency offset determination and synchronisation
-
- 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/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)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明公开了一种基于数据辅助的频率跟踪方法、系统及存储介质，涉及无线通信领域，解决了现有技术中存在频率跟踪方法与调制方式耦合性强、跟踪速度和精度难以满足的缺陷，其技术方案要点是：对输入的数据进行处理，得到需要进行频率跟踪的目标数据，解扩该目标数据，对解扩后的目标数据进行频偏估计和校正处理，获得信号数据，然后对信号数据进行共轭相乘，对相乘的结果进行判决，生成期望数据，取期望数据的共轭将其与信号数据相乘，对相乘后的结果进行误差提取处理，最终对提取的误差进行环路反馈及环路输出。本发明提供了可运用于多种调制方式、跟踪速度快且精度高的频率跟踪方法。
Description
技术领域
本发明涉及一种无线通信领域中的频率跟踪领域，更具体地说，它涉及一种基于数据辅助的频率跟踪方法、系统及存储介质。
背景技术
低轨卫星移动通信系统有传输距离远、频率偏差大、传输信道多变的特点，为实现远距离信息传输，发送信号通常采用扩频加数字调制的方式进行无线传输，信号数字调制方式随传输信道质量变化自适应的进行调整，常用的调制方式有BPSK、QPSK、8PSK、16APSK、32APSK等，在接收端，为正确获取对方信息，在定时同步后通常采用载波跟踪同步后解扩频的接收方式。
现有技术中，先载波同步跟踪再进行解扩频接收方式，无法在载波跟踪中充分利用信号的扩频增益，载波跟踪性能受信号噪声影响，在低信噪比环境下，载波跟踪性能大大降低，甚至可能出现无法载波同步的情况。
现有技术中的频率跟踪算法实现方式大致可以分为两种：前馈方式和反馈方式。前馈方式通常根据调制方式采用对应的高阶谱方法，测频精度高，但适应频偏范围有限。在信噪比较低时，需要获取最佳采样点以将采样率降至最低，该方法有利于克服噪声的影响，高阶谱曲线显示出更加明显的谱线特征。反馈方式基于锁相环(PLL)，通常PLL直接对载波相位进行跟踪，PLL方式当环路稳定闭环时具有较高的跟踪精度，但在高动态环境下，采用PLL跟踪的高动态扩频接收机需要承受环路带宽与动态性能之间的折衷，即噪声引入跟踪误差随环路带宽降低而增加，难以满足跟踪精度与动态性能的要求。
发明内容
本发明的目的是提供一种基于数据辅助的频率跟踪方法、系统及存储介质，以解决现有技术中存在频率跟踪方法与调制方式耦合性强、跟踪速度和精度难以满足的缺陷。
本发明的上述技术目的是通过以下技术方案得以实现的：
第一方面，提供了一种基于数据辅助的频率跟踪方法，包括以下步骤：
步骤一，输入经匹配滤波处理后的数据，将所述输入的数据与本地序列的同步头进行相关运算，获得所述输入的数据的目标数据；
步骤二，解扩所述目标数据，获得解扩后的目标数据序列；
步骤三，对所述目标数据序列进行快速傅里叶变换，获得所述目标数据序列的频偏值，对所述频偏值进行频偏校正处理，获得频偏校正后的信号数据；
步骤四，输入初始值为零的环路误差值，对所述环路误差值进行处理，获得第一估计值，对所述信号数据与第一估计值进行处理，获得输出信号，对所述输出信号的相位值进行处理，获得期望数据；
步骤五，提取所述信号数据与期望数据的相位误差值，对所述相位误差值进行环路反馈处理，获得环路输出值；
步骤六，将所述环路输出值与所述初始值为零的环路误差值进行求和，更新所述环路输出值，并将更新后的所述环路输出值返回至步骤四进行循环迭代。
与现有技术比较而言，本发明频偏估计及频偏校正对输入的数据的粗略估计和校正，降低数据的误差值，由于对述输出信号的相位值的处理，只是对相位所处区域的相位值进行判决处理，其不受限于信号的调制方式，因此本发明所提及的不受限于调制方式的频率跟踪则可以实现，采用PLL锁相环方式，来跟踪输入的数据与理想解调信号之间的相位差，并不断跟踪该误差值，达到频率跟踪的动态平衡状态，从而避免了将噪声引入后误差随环路带宽降低而增加，满足了接收机的跟踪精度与动态性能。
进一步的，获取本地序列的同步头的前100个符号，将所述前100个符号进行移位和自相关处理，然后求和取模，对模值进行最大值判决，进行自适应处理查找获得所述数据的目标数据的数据起始位。
进一步的，在进行解扩处理之前，还需对所述目标数据进行扩频处理，扩频的方式为比特重复扩频。
进一步的，求取所述目标数据序列的数字模平方，然后在进行快速傅里叶变换并取变换结果的绝对值，获取所述目标数据序列的采样频率，根据所述采样频率和快速傅里叶变换的点数将所述绝对值的最大值索引号转换为频偏值，完成对所述目标数据序列的频偏估计。
进一步的，根据所述采样频率确认每个采样点对应的采样时间，根据所述采样时间和采样点计算每个采样点频偏值对应的相位补偿因子，将所述目标数据序列的每个信号数据点乘以所述相位补偿因子的共轭，完成频偏校正。
进一步的，获取所述输出信号的相位所处区域，判定所述相位所处区域的期望相位值，并形成所述期望数据。
进一步的，对所述期望数据和所述信号数据进行相乘处理，并对相乘后的数据进行误差提取处理，提取所述相位误差值的计算式为：err_temp＝Im(data_temp*conj(data_predict))，式中，Im表示取虚部，conj表示取共轭，data_temp表示信号数据，data_predict表示期望数据。
进一步的，对所述相位误差值进行环路反馈处理的计算式为：loop_out_temp＝g1*err_temp+g2*acc_err_temp，式中，g1表示相位系数，g2表示累计相位系数，err_temp表示本次相位误差值，acc_err_temp表示累计相位误差值。
第二方面，提供了一种基于数据辅助的频率跟踪系统，所述系统包括：
第一同步单元，用于输入经匹配滤波处理后的数据，将所述输入的数据与本地序列的同步头进行相关运算，获得所述输入的数据的目标数据；
解扩单元，用于解扩所述目标数据，获得解扩后的目标数据序列；
频偏单元，用于对所述目标数据序列进行快速傅里叶变换，获得目标数据序列的频偏值，对所述频偏值进行频偏校正处理，获得频偏校正后的信号数据；
第二同步单元，用于输入初始值为零的环路误差值，对所述环路误差值进行处理，获得第一估计值，对所述信号数据与第一估计值进行处理，获得输出信号，对所述输出信号的相位值进行处理，获得期望数据；
第三同步单元，用于提取所述信号数据与期望数据的相位误差值，对所述相位误差值进行环路反馈处理，获得环路输出值；
跟踪单元，用于将所述环路输出值与所述初始值为零的环路误差值进行求和，更新所述环路输出值，并将更新后的所述环路输出值返回至第二同步单元进行循环迭代。
第三方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现如第一方面所述的方法中所执行的操作。
与现有技术相比，本发明具有以下有益效果：
1.本发明频偏估计及频偏校正对输入的数据的粗略估计和校正，降低数据的误差值，由于对述输出信号的相位值的处理，只是对相位所处区域的相位值进行判决处理，其不受限于信号的调制方式，因此本发明所提及的不受限于调制方式的频率跟踪则可以实现，采用PLL锁相环方式，来跟踪输入的数据与理想解调信号之间的相位差，并不断跟踪该误差值，达到频率跟踪的动态平衡状态。
2.本发明先对已经扩频后的数据序列进行解扩，然后再对解扩后的数据进行频偏的估计及校正，最后再利用载波估计，即用于输入初始值为零的环路误差值，对环路误差值进行处理，获得第一估计值，对信号数据与第一估计值进行处理，这样即可避免无法在载波跟踪中充分利用信号的扩频增益，使得载波跟踪性能受信号噪声影响。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明实施例提供的方法流程图；
图2为本发明实施例提供的频率跟踪方法的数据走向图；
图3为本发明实施例提供的FFT相关峰图；
图4为本发明实施例提供的数据解扩前信号星座图；
图5为本发明实施例提供的数据解扩后信号星座图；
图6为本发明实施例提供的载波输出典型曲线图；
图7为本发明实施例提供的环路滤波输出典型曲线图；
图8为本发明实施例提供的经频率跟踪补偿后数据的星座图；
图9为本发明实施例提供的系统的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例一：
如图1所示，本申请实施例一提供了一种基于数据辅助的频率跟踪方法，可应用于低轨卫星通信系统；也可以应用于高速运动载荷之间通信，比如高铁，高速飞行器等。包括以下步骤：
步骤一，输入经匹配滤波处理后的数据，将输入的数据与本地序列的同步头进行相关运算，获得输入的数据的目标数据；
步骤二，解扩目标数据，获得解扩后的目标数据序列；
步骤三，对目标数据序列进行快速傅里叶变换，获得目标数据序列的频偏值，对频偏值进行频偏校正处理，获得频偏校正后的信号数据；
步骤四，输入初始值为零的环路误差值，对环路误差值进行处理，获得第一估计值，对信号数据与第一估计值进行处理，获得输出信号，对输出信号的相位值进行处理，获得期望数据；
步骤五，提取信号数据与期望数据的相位误差值，对相位误差值进行环路反馈处理，获得环路输出值；
步骤六，将环路输出值与初始值为零的环路误差值进行求和，更新环路输出值，并将更新后的环路输出值返回至步骤四进行循环迭代。
具体的，如图2所示，图2为频率跟踪方法的数据走向图，步骤一是相当于对输入数据进行粗定时同步，是获取经匹配滤波后的数据即输入的数据前面的同步头，获取到同步头也就得到数据的具体位置。粗定时同步处理过程如下：输入数据与本地序列的同步头实部符号位进行相关运算，获得数据的相关峰位置，相关峰位置即为输入的数据的同步头精确位置，同步头后面的是信号中所需要处理的数据。本发明所提及的数据辅助是指，利用同步头获取数据起始位，此过程便获取到后续处理步骤所需要的数据，即目标数据。
即对步骤一粗定时同步得到的数据进行解扩处理，解扩的数据是已被扩频的数据，然后进行频偏估计和频偏校正，主要是避免现有技术中先进行载波同步跟踪再进行解扩频接收方式。
频偏估计及频偏校正是对信号频偏的粗略估计和校正，信号中所残余的频偏需要用到精确频率补偿。
基于频偏估计及校正后得到的信号数据(即data_temp)，在步骤四中输入初始值为0的环路误差(即NCO_TEMP输入初值为0)，如图2所示，将其进行载波估计处理后，得到第一估计值，在与信号数据(即data_temp)进行相乘得到输出信号(即data_out)，对输出信号进行判决处理获得期望数据(即data_predict)，在将期望数据与信号数据进行共轭相乘处理，提取输入的数据与理想解调信号之间的相位误差值，然后进行环路反馈跟踪，循环迭代步骤四至步骤六，使得跟踪的相位误差达到动态平衡状态，并且将误差状态维持在一定水平。
本申请实施例一的又一个实施例中，获取本地序列的同步头的前100个符号，将所述前100个符号进行移位和自相关处理，然后求和取模，对模值进行最大值判决，进行自适应处理查找获得所述数据的目标数据的数据起始位。
具体的，相关运算是指选取本地同步头的前100个符号，将此序列进行移位后自相关，然后求和取模，对模值进行最大值判决，自适应查找相关峰位置，相关峰位置即为输入数据的同步头精确位置，同步头后面的是信号中所需要处理的数据。
本申请实施例一的又一个实施例中，在进行解扩处理之前，还需对目标数据进行扩频处理，扩频的方式为比特重复扩频。
具体的，解扩则根据扩频方式进行数据处理，本发明在数据发送到信道过程中采取扩频，扩频方式为比特重复扩频，扩频比为n，解扩的计算式如下式(1)：
despread(i)＝sum(din(4*(i-1)+1:4*i).*pnn4)(1)，式中，despread是解扩后数据序列，sum表示求和，din表示解扩前信号数据序列，pnn4表示解扩点乘矩阵。
本申请实施例一的又一个实施例中，求取目标数据序列的数字模平方，然后在进行快速傅里叶变换并取变换结果的绝对值，获取目标数据序列的采样频率，根据采样频率和快速傅里叶变换的点数将绝对值的最大值索引号转换为频偏值，完成对目标数据序列的频偏估计。
具体的，频偏估计是基于解扩后数据的流程，即对解扩后的数据进行数字模平方，然后在进行FFT变换再取绝对值，查找到最大值索引号，最后根据下式(2)将该索引号转换为频偏值。
上式中，f0是频偏，index是索引号，NFFT是FFT变换的点数，fs是信号采样频率。
本申请实施例一的又一个实施例中，根据采样频率确认每个采样点对应的采样时间，根据采样时间和采样点计算每个采样点频偏值对应的相位补偿因子，将目标数据序列的每个信号数据点乘以相位补偿因子的共轭，完成频偏校正。
具体的，频偏校正则是对所估计的频偏进行校正的过程。首先根据采样频率确认每个采样点对应的时间，如下式(3)所示。
上式中，n是采样点数，fs是信号采样频率。
然后计算每个采样点频偏对应的相位补偿因子，如下式(4)所示。
最后将信号数据每个点乘以相位补偿因子的共轭，即完成频偏校正，如下式(5)所示。
rx＝rx*conj(carr) (5)。
本申请实施例一的又一个实施例中，获取输出信号的相位所处区域，判定相位所处区域的期望相位值，并形成期望数据。
具体的，在获得输出信号之前还需进行的处理如下，如图2中的NCO_TEMP是环路误差输入，其初始值是0。
载波估计中，对初始值是0的NCO_TEMP取相位及共轭，如下式6所示。
conj(ej*NCO_TEMP) (6)
信号共轭相乘中，将信号数据即data_temp与载波估计得到的值进行相乘得到data_out，如下式7所示。
data_out＝data_temp*conj(ej*NCO_TEMP) (7)。
根据data_out信号的相位所处区域，判定其区域内的期望相位值，并形成期望数据即data_predict。例如：当信号相位处于[-pi/4，pi*3/4]区间时，判决其期望相位值为pi/4。由于算法中此处判决不受限于信号的调制方式，因此本发明的频率跟踪方法不受限于信号的调制方式。
本申请实施例一的又一个实施例中，对期望数据和信号数据进行相乘处理，并对相乘后的数据进行误差提取处理，提取相位误差值的计算式为：
err_temp＝Im(data_temp*conj(data_predict)) (8)，式中，Im表示取虚部，conj表示取共轭，data_temp表示信号数据，data_predict表示期望数据。
具体的，如图2所示，先去期望数据的共轭，然后再与信号数据相乘，在提取相乘后的信号数据与期望数据的相位误差值，即相当于提取输入的数据与理想解调信号之间的相位误差值。
本申请实施例一的又一个实施例中，对相位误差值进行环路反馈处理的计算式为：loop_out_temp＝g1*err_temp+g2*acc_err_temp(9)，式中，g1表示相位系数，g2表示累计相位系数，err_temp表示本次相位误差值，acc_err_temp表示累计相位误差值。
具体的，如图2所示，环路反馈中，对环路误差根据式(9)进行处理，获得环路输出值和环路输出值，将环路输出值。
本实施例一还提供了采用本申请的方法后的数据的星座图，请参考图3、4、5、6、7、和8，所采用的初始数据为经匹配滤波后得到的星座点数据，数据量为2400000个，形式为1x2400000的矩阵。本发明实例所用数据基本特征为：实例信号采用BPSK调试方式，符号速率30MHz，采样点数2048，频率偏移650KHz，采样率120MHz。
粗定时同步中，根据本地同步头与输入数据做相关运算后，匹配到输入数据中的同步头，形式为1x2048矩阵。
数据解扩中，式(1)中的点乘矩阵为(1 1 1 1)。对比解扩前后信号数据星座图，如图4和图5，可以明显看到解扩后信号数据星座图效果更好。
频偏估计及校正中，频偏估计相关峰曲线如图3所示。可以看到，相关峰较为明显，进行频率粗补偿即频率校正之后残余误差仅为24.4141Hz。
反馈精确频率同步中的式(9)中：g1取值为0.05333,g2取值为0.00142。图6和图7是载波输出典型曲线和环路滤波输出典型曲线，反应环路跟踪动态效果良好。图8为频率跟踪后的信号数据星座图，可以看出频率跟踪后星座图数据点位收敛效果较好。
实施例二：
基于同一构思，本申请实施例二在实施例一的基础上提供了一种基于数据辅助的频率跟踪系统，系统包括：
第一同步单元10，用于输入经匹配滤波处理后的数据，将输入的数据与本地序列的同步头进行相关运算，获得输入的数据的目标数据；
解扩单元20，用于解扩目标数据，获得解扩后的目标数据序列；
频偏单元30，用于对目标数据序列进行快速傅里叶变换，获得目标数据序列的频偏值，对频偏值进行频偏校正处理，获得频偏校正后的信号数据；
第二同步单元40，用于输入初始值为零的环路误差值，对环路误差值进行处理，获得第一估计值，对信号数据与第一估计值进行处理，获得输出信号，对输出信号的相位值进行处理，获得期望数据；
第三同步单元50，用于提取信号数据与期望数据的相位误差值，对相位误差值进行环路反馈处理，获得环路输出值；
跟踪单元60，用于将环路输出值与初始值为零的环路误差值进行求和，更新环路输出值，并将更新后的环路输出值返回至第二同步单元进行循环迭代。
本实施例二的又一个实施例中，系统还包括相关运算单元，用于获取本地序列的同步头的前100个符号，将所述前100个符号进行移位和自相关处理，然后求和取模，对模值进行最大值判决，进行自适应处理查找获得所述数据的目标数据的数据起始位。
本实施例二的又一个实施例中，系统还包括扩频单元，在进行解扩处理之前，用于对目标数据进行扩频处理，扩频的方式为比特重复扩频。
本实施例二的又一个实施例中，系统还包括频偏估计单元，用于求取目标数据序列的数字模平方，然后在进行快速傅里叶变换并取变换结果的绝对值，获取目标数据序列的采样频率，根据采样频率和快速傅里叶变换的点数将绝对值的最大值索引号转换为频偏值，完成对目标数据序列的频偏估计。
本实施例二的又一个实施例中，系统还包括频偏校正单元，用于根据采样频率确认每个采样点对应的采样时间，根据采样时间和采样点计算每个采样点频偏值对应的相位补偿因子，将目标数据序列的每个信号数据点乘以相位补偿因子的共轭，完成频偏校正。
本实施例二的又一个实施例中，系统还包括判定单元，用于获取输出信号的相位所处区域，判定相位所处区域的期望相位值，并形成期望数据。
本实施例二的又一个实施例中，系统还包括提取单元，用于对期望数据和信号数据进行相乘处理，并对相乘后的数据进行误差提取处理，提取相位误差值的计算式为：err_temp＝Im(data_temp*conj(data_predict))，式中，Im表示取虚部，conj表示取共轭，data_temp表示信号数据，data_predict表示期望数据。
本实施例二的又一个实施例中，系统还包括反馈单元，用于对相位误差值进行环路反馈处理的计算式为：loop_out_temp＝g1*err_temp+g2*acc_err_temp，式中，g1表示相位系数，g2表示累计相位系数，err_temp表示本次相位误差值，acc_err_temp表示累计相位误差值。
本实施例二中系统的各个单元所实现的均已在上述实施例一的方法中叙述清楚，因此，实施例二中不再一一叙述。
实施例三：
本申请实施例三还提供了一种计算机可读存储介质，计算机可读存储介质中存储有至少一条指令，该指令由处理器加载并执行时，使得计算机执行如上述实施例一中的方法所执行的操作。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于数据辅助的频率跟踪方法，其特征在于，包括以下步骤：
步骤一，输入经匹配滤波处理后的数据，将所述输入的数据与本地序列的同步头进行相关运算，获得所述输入的数据的目标数据；其中，结合经匹配滤波处理的输入的数据的同步头与本地序列的同步头的实部符号位进行相关运算，获得相关峰位置，以相关峰位置作为输入的数据的数据起始位，以数据起始位后的数据作为输入的数据的目标数据；
步骤二，解扩所述目标数据，获得解扩后的目标数据序列；
步骤三，对所述目标数据序列进行快速傅里叶变换，获得所述目标数据序列的频偏值，对所述频偏值进行频偏校正处理，获得频偏校正后的信号数据；
步骤四，输入初始值为零的环路误差值，对所述环路误差值进行处理，获得第一估计值，对所述信号数据与第一估计值进行处理，获得输出信号，对所述输出信号的相位值进行处理，获得期望数据；其中，对输入初始值为零的环路误差值进行载波估计处理后，得到第一估计值；将第一估计值与信号数据进行相乘得到输出信号；对输出信号进行判决处理获得期望数据，其中判决处理具体为根据输出信号的相位所处区域，判定相位所处区域的期望相位值，依据期望相位值形成所述期望数据；
步骤五，提取所述信号数据与期望数据的相位误差值，对所述相位误差值进行环路反馈处理，获得环路输出值；
步骤六，将所述环路输出值与所述初始值为零的环路误差值进行求和，更新所述环路输出值，并将更新后的所述环路输出值返回至步骤四进行循环迭代。
2.根据权利要求1所述的一种基于数据辅助的频率跟踪方法，其特征在于，获取本地序列的同步头的前100个符号，将所述前100个符号进行移位和自相关处理，然后求和取模，对模值进行最大值判决，进行自适应处理查找获得所述数据的目标数据的数据起始位。
3.根据权利要求1所述的一种基于数据辅助的频率跟踪方法，其特征在于，在进行解扩处理之前，还需对所述目标数据进行扩频处理，扩频的方式为比特重复扩频。
4.根据权利要求1所述的一种基于数据辅助的频率跟踪方法，其特征在于，求取所述目标数据序列的数字模平方，然后在进行快速傅里叶变换并取变换结果的绝对值，获取所述目标数据序列的采样频率，根据所述采样频率和快速傅里叶变换的点数将所述绝对值的最大值索引号转换为频偏值，完成对所述目标数据序列的频偏估计。
5.根据权利要求4所述的一种基于数据辅助的频率跟踪方法，其特征在于，根据所述采样频率确认每个采样点对应的采样时间，根据所述采样时间和采样点计算每个采样点频偏值对应的相位补偿因子，将所述目标数据序列的每个信号数据点乘以所述相位补偿因子的共轭，完成频偏校正。
6.根据权利要求1所述的一种基于数据辅助的频率跟踪方法，其特征在于，对所述期望数据和所述信号数据进行相乘处理，并对相乘后的数据进行误差提取处理，提取所述相位误差值的计算式为：
，式中，Im表示取虚部，conj表示取共轭，data_temp表示信号数据，data_predict表示期望数据。
7.根据权利要求6所述的一种基于数据辅助的频率跟踪方法，其特征在于，对所述相位误差值进行环路反馈处理的计算式为：
，式中，g1表示相位系数，g2表示累计相位系数，err_temp表示本次相位误差值，acc_err_temp表示累计相位误差值。
8.一种基于数据辅助的频率跟踪系统，其特征在于，所述系统包括：
第一同步单元，用于输入经匹配滤波处理后的数据，将所述输入的数据与本地序列的同步头进行相关运算，获得所述输入的数据的目标数据；其中，结合经匹配滤波处理的输入的数据的同步头与本地序列的同步头的实部符号位进行相关运算，获得相关峰位置，以相关峰位置作为输入的数据的数据起始位，以数据起始位后的数据作为输入的数据的目标数据；
解扩单元，用于解扩所述目标数据，获得解扩后的目标数据序列；
频偏单元，用于对所述目标数据序列进行快速傅里叶变换，获得目标数据序列的频偏值，对所述频偏值进行频偏校正处理，获得频偏校正后的信号数据；
第二同步单元，用于输入初始值为零的环路误差值，对所述环路误差值进行处理，获得第一估计值，对所述信号数据与第一估计值进行处理，获得输出信号，对所述输出信号的相位值进行处理，获得期望数据；其中，对输入初始值为零的环路误差值进行载波估计处理后，得到第一估计值；将第一估计值与信号数据进行相乘得到输出信号；对输出信号进行判决处理获得期望数据，其中判决处理具体为根据输出信号的相位所处区域，判定相位所处区域的期望相位值，依据期望相位值形成所述期望数据；
第三同步单元，用于提取所述信号数据与期望数据的相位误差值，对所述相位误差值进行环路反馈处理，获得环路输出值；
跟踪单元，用于将所述环路输出值与所述初始值为零的环路误差值进行求和，更新所述环路输出值，并将更新后的所述环路输出值返回至第二同步单元进行循环迭代。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现如权利要求1至7任一项所述的方法中所执行的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111565220.5A CN114244673B (zh) | 2021-12-20 | 2021-12-20 | 一种基于数据辅助的频率跟踪方法、系统及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111565220.5A CN114244673B (zh) | 2021-12-20 | 2021-12-20 | 一种基于数据辅助的频率跟踪方法、系统及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114244673A CN114244673A (zh) | 2022-03-25 |
| CN114244673B true CN114244673B (zh) | 2024-05-03 |
Family
ID=80759641
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111565220.5A Active CN114244673B (zh) | 2021-12-20 | 2021-12-20 | 一种基于数据辅助的频率跟踪方法、系统及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114244673B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109150233A (zh) * | 2018-09-13 | 2019-01-04 | 南京理工大学 | 一种直扩dpsk信号的调制解调方法 |
| CN112260977A (zh) * | 2020-10-19 | 2021-01-22 | 清华大学 | 一种动态卫星通信系统信号跟踪方法及装置 |
| CN112511182A (zh) * | 2020-11-25 | 2021-03-16 | 深圳航天东方红卫星有限公司 | 一种星载测控数传接收装置 |
| CN113238261A (zh) * | 2021-05-31 | 2021-08-10 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 低轨卫星扩频通信体制信号捕获跟踪系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106453187B (zh) * | 2016-09-26 | 2019-06-11 | 珠海全志科技股份有限公司 | 一种ofdm系统同步跟踪方法和装置 |
-
2021
- 2021-12-20 CN CN202111565220.5A patent/CN114244673B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109150233A (zh) * | 2018-09-13 | 2019-01-04 | 南京理工大学 | 一种直扩dpsk信号的调制解调方法 |
| CN112260977A (zh) * | 2020-10-19 | 2021-01-22 | 清华大学 | 一种动态卫星通信系统信号跟踪方法及装置 |
| CN112511182A (zh) * | 2020-11-25 | 2021-03-16 | 深圳航天东方红卫星有限公司 | 一种星载测控数传接收装置 |
| CN113238261A (zh) * | 2021-05-31 | 2021-08-10 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 低轨卫星扩频通信体制信号捕获跟踪系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114244673A (zh) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110545247B (zh) | 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法 | |
| US5875215A (en) | Carrier synchronizing unit | |
| CN107241286B (zh) | 一种结合频域和时域估计的载波同步方法 | |
| KR100355326B1 (ko) | Ｏｆｄｍ 통신 장치 및 전파로 추정 방법 | |
| US6628926B1 (en) | Method for automatic frequency control | |
| CN107135175B (zh) | 一种高速数传系统中同步均衡方法 | |
| KR100318952B1 (ko) | 무선통신시스템에서 채널 추정을 위한 장치 및 방법 | |
| WO2009016000A2 (en) | Iq gain imbalance suppression in presence of unknown phase imbalance | |
| WO1993022861A1 (en) | Receiver for digital communication system | |
| JP6559088B2 (ja) | 通信機および衛星通信システム | |
| CN113037671A (zh) | 一种低复杂度的高效soqpsk符号定时与相位联合同步算法 | |
| EP0838111B1 (en) | Carrier recovery for digitally phase modulated signals, using a known sequence | |
| CN114244673B (zh) | 一种基于数据辅助的频率跟踪方法、系统及存储介质 | |
| KR100840608B1 (ko) | 인터리브된 주파수 분할 다중 접속을 위한 주파수 오프셋보상 방법 및 장치 | |
| CN101116257A (zh) | 用于确定相关值的设备和方法 | |
| CN112953593B (zh) | LoRa高级接收器 | |
| US7428272B2 (en) | Method and apparatus for joint phase offset and frequency offset estimator for MPSK transmission | |
| CN110943949B (zh) | 一种基于sc-fde传输体制的宽带角跟踪方法 | |
| US7116727B2 (en) | Frequency offset estimation apparatus for intersymbol interference channels | |
| KR100585965B1 (ko) | 수신기 위상 보정장치 및 그 방법 | |
| CN100369386C (zh) | 一种结合频偏补偿的联合检测方法和装置 | |
| CN116016049B (zh) | 一种适合高动态低信噪比的频偏估计方法 | |
| CN114095317B (zh) | 一种适用于高阶apsk调制的载波参数估计解耦合方法和系统 | |
| KR20040036492A (ko) | 무선 통신시스템의 송수신기 성능을 개선하기 위한 방법 | |
| KR100569830B1 (ko) | 회전 큐피에스케이 신호를 이용한 위상편이 추정 및 보상방법 및 그 장치 |
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 |