CN115996162A - 一种串行高效通信时域o&m定时同步方法 - Google Patents

一种串行高效通信时域o&m定时同步方法 Download PDF

Info

Publication number
CN115996162A
CN115996162A CN202211474592.1A CN202211474592A CN115996162A CN 115996162 A CN115996162 A CN 115996162A CN 202211474592 A CN202211474592 A CN 202211474592A CN 115996162 A CN115996162 A CN 115996162A
Authority
CN
China
Prior art keywords
timing
time domain
timing error
signal
filter
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
Application number
CN202211474592.1A
Other languages
English (en)
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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202211474592.1A priority Critical patent/CN115996162A/zh
Publication of CN115996162A publication Critical patent/CN115996162A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明公开了一种串行高效通信时域O&M定时同步方法，属于太赫兹毫米波高速通信领域。本方法使用ADC完成信号采样，数据捕获后提取复数序列到定时同步模块进行处理；计算该导频序列模方即包络幅值，将该结果分别送入IQ两路完成DFT运算；利用DFT运算后的实部和虚部，计算得到定时误差估计值；依据定时误差计算结果求出滤波器抽头位置，用来索引多相滤波器系数，多相滤波器在时域完成定时误差补偿和信号数据速率转换。本发明采用开环IQ两路并行和符号非整数倍采样设计，IQ各路串行处理样本数据，并在时域上进行定时误差补偿，避免DFT和IDFT运算的同时进一步降低了硬件资源的占用。

Description

一种串行高效通信时域O＆M定时同步方法
技术领域
本发明涉及太赫兹毫米波通信领域，具体涉及一种通信位同步的串行高效时域O&M定时同步方法。
背景技术
在软件无线电的领域里，信号经过信道时由于传输时延的存在以及收发端时钟偏移，都会造成接收到的信号产生偏差。信号接收端经过采样下变频处理后的数据是一串码元序列，只有调整接收端采样时钟使其与接收码元频率相位一致，才能在接收符号的最佳采样时刻对其进行采样判决。
位同步按有无数据辅助可分为数据辅助和非数据辅助两大类，后者不需要额外的辅助序列，相比于数据辅助的方法其频谱利用率更高，在工程中应用更广泛。非数据辅助位同步算法里比较有代表性的是Gardner算法和O&M算法。O&M算法是OerderM和Meyr H在1988年提出的利用平方律非线性变换对定时误差进行估计的一种算法。
但是，传统的O&M定时同步算法不满足高速率大带宽数据传输要求。
发明内容
本发明的目的在于提供一种串行高效通信时域O&M定时同步方法，可以从运算架构上压缩硬件逻辑资源的使用，在保证超高速调制速率传输的同时降低资源消耗，缩短定时同步时间。
本发明所采用的技术方案如下：
一种串行高效通信时域O&M定时同步方法，包括以下步骤：
步骤1，对信号进行采样，完成采样后数据的帧同步，对帧同步后的复数序列进行定时同步处理，得到复数定时同步序列；
步骤2，对复数定时同步序列进行I、Q两路串行输入定时处理，每路均以二进制符号形式按位相乘并累计求和，得到平方值，然后将I、Q两路得到的结果求和，得到接收信号模方；
步骤3，将信号包络值分别送入I、Q两路，每路信号分别与相应的旋转因子WnI和Wnq相乘并求和，完成DFT运算；
步骤4，利用信号DFT运算后的实部和虚部，进行定时误差估计，计算出频域相位偏移量，即定时误差估计结果；
步骤5，依据定时误差估计结果进行环路滤波平滑，计算出定时同步时域延迟，对应到滤波器抽头位置用于后续处理；
步骤6，不同定时误差估计结果对应的滤波器系数存储在ROM中，将步骤5输出的滤波器抽头位置作为ROM索引地址，动态控制ROM输出滤波器系数；
步骤7，将ROM输出的滤波器系数输入到多相滤波器模块，完成采样信号定时误差补偿。
进一步地，步骤1的具体方式为，使用与FPGA基带板卡相连的模数转换模块对信号进行采样，采样后的信号在基带板卡里完成帧同步，输出结果中包括从接收采样数据中提取出的导频序列，将该序列用于定时误差估计。
进一步地，步骤3的DFT运算中，采用满足奈奎斯特采样定理的采样速率，信号为复数形式，DFT运算表示为：
Figure BDA0003959091350000021
其中，xk为接收信号模方；Xm为xk的傅里叶变换，k为样本序号，L为符号累积长度，m为样本段号；I、Q路的傅里叶变换旋转因子分别表示为WnI＝Im[e-j2kπ/N]和Wnq＝Re[e-j2kπ/N]；N为DFT点数。
进一步地，步骤4中，Xm的实部对应采样点模方与旋转因子WnI的乘积，Xm的虚部对应采样点模方与旋转因子Wnq的乘积，符号累积长度为L，I、Q两路分别进行串行处理；Xm的实部、虚部分别表示为：
Figure BDA0003959091350000031
Figure BDA0003959091350000032
其中，Re[Xm]和Im[Xm]分别表示Xm的实部和虚部，Ik和Qk分别为接收样本信号的实部和虚部；
DFT运算由其接收信号的实部、虚部信号运算结果表示，L个符号的定时误差估计值为
Figure BDA0003959091350000033
其中，误差估计值
Figure BDA0003959091350000034
为ε的无偏估计。
进一步地，步骤5中，用于定时误差补偿的滤波器系数抽头位置的计算方式为
Figure BDA0003959091350000035
其中，
Figure BDA0003959091350000036
为定时误差估计值，fus为根生余弦滤波器符号上采样倍数。
进一步地，步骤7的具体方式为，将滤波器系数与信号样点在时域上进行卷积，滤波器采用多相结构设计，利用多相滤波器抽取和插值倍数的不同实现信号定时误差补偿，同时完成数据速率转换。
本发明的有益效果为：
1、在吉赫兹超大带宽、超高速调制解调***下，传统的平方律定时同步算法不再适用。本发明采用IQ路并行和符号非整数倍采样设计，在硬件采样速率受限下，既满足传输速率要求，又避免了符号运算点数过多的问题；开环的定时结构，相对于闭环结构具有定时同步时间短，算法实现简单等优点。
2、本发明解决了ADC无法在符号最佳采样时刻采样的问题，提升了接收信号的信噪比。相比于现有算法，具有控制信号反馈的结构势必会增加同步所需时间，收敛速度变慢，不适用于突发信号的解调。本发明仅需提取部分接收数据的导频序列用于定时同步算法处理，满足了高速通信的实时性和高效性；IQ路分别采用串行处理架构，在满足传输速率的同时也降低了硬件内部逻辑资源的使用，保证算法性能稳定。此外，本发明提出了一种新的误差校正方式，根据定时误差估计值和滤波器上采样倍数估计滤波器系数抽头位置，并在时域进行定时误差估计补偿，同时也省去了DFT和IDFT模块，有效减少了硬件资源的占用，缩短了定时同步算法所需时间。
附图说明
图1是本发明方法的流程图。
图2是本发明方法的原理示意图。
图3是本发明处理过程中的输入信号模方频谱图。
图4是本发明处理前后导频序列QPSK调制星座点对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例子，并且参照附图，对本发明进一步详细说明。
一种串行高效通信时域O&M定时同步方法，本方法使用ADC完成信号采样，数据捕获后提取复数序列到定时同步模块进行处理；计算该导频序列模方即包络幅值，将该结果分别送入IQ两路完成DFT运算；利用DFT运算后的实部和虚部，计算得到定时误差估计值；依据定时误差计算结果求出滤波器抽头位置，用来索引多相滤波器系数，多相滤波器在时域完成定时误差补偿和信号数据速率转换。
具体来说，本方法包括以下步骤：
步骤1，使用与FPGA基带板卡相连的ADC模数转换器完成信号采样，采样处理后的信号先进行数据捕获同步，捕获后的信号确定了帧头，提取复数导频序列并将其输出到定时同步模块进行定时误差估计。
步骤2，接收信号r(t)经过接收端N倍上采样后得到的样本为：
Figure BDA0003959091350000051
其中，an为发送数据符号，g(t)为成形滤波器、信道和接收滤波器的总响应，ε为定时误差。信号采样捕获后，提取的导频序列用正交复矢量可以表示为rk＝Ik+jQk，其中k＝0,1,2,...,N，提取的导频序列数据用于定时误差估计。
复数导频序列通过I、Q两路串行输入定时模块，各路采样二进制符号形式按位相乘并累计求和得到其平方值，将I路、Q路得到的结果求和即为接收信号模方，该包络值可用
Figure BDA0003959091350000052
表示。
步骤3，将得到的信号包络平方分别送入I、Q两路，每路信号分别与相应的旋转因子WnI和Wnq相乘并累积求和，完成DFT运算。本方法采样速率满足奈奎斯特采样定理，信号为复信号形式，其DFT变换可表示为：
Figure BDA0003959091350000061
其中WnI＝Im[e-j2kπ/N]和Wnq＝Re[e-j2kπ/N]分别为旋转因子，N为DFT点数。
步骤4，对信号模方xk求DFT运算得到Xm的实部和虚部，Xm的实部对应采样点模方与旋转因子WnI的乘积，Xm的虚部对应采样点模方与旋转因子Wnq的乘积，符号累积长度为L，I、Q各路串行处理，则Xm的实部虚部分别表示为
Figure BDA0003959091350000062
Figure BDA0003959091350000063
至此，DFT运算可由其接收信号的实部虚部信号运算结果表示，计算可得L个符号的定时误差估计值为
Figure BDA0003959091350000064
求出频域相位偏移量，得到定时误差估计结果。
步骤5，依据定时误差估计结果进行环路滤波平滑并计算得到定时同步时域延迟，用于定时误差补偿的滤波器系数抽头位置计算公式为
Figure BDA0003959091350000065
其中
Figure BDA0003959091350000066
为定时误差估计值，fus为根生余弦滤波器符号上采样倍数，计算得到的结果对应系数在ROM里的存储地址，索引输出的滤波器系数用于后续补偿处理。
步骤6，不同定时误差估计结果对应的滤波器系数存储在ROM中，将步骤5输出的滤波器抽头位置作为ROM索引地址，动态地控制ROM输出滤波器系数。
步骤7，将索引ROM得到的滤波器系数输入多相滤波器模块，与接收信号采样结果在时域上进行卷积，在完成采样信号定时误差补偿的同时实现数据速率转换。
本方法充分考虑了硬件采样速率的限制，在满足数据传输速率的前提下，采用开环IQ路并行和符号非整数倍采样设计，IQ各路串行处理样本的O&M算法有利于节约硬件内部资源，实现高效稳定的数据传输。O&M定时误差结果送到滤波器进行时域定时误差补偿，时域补偿相比频域实现省去了DFT和IDFT的步骤，进一步压缩了实现算法所需的资源数量。
下面为一个更具体的例子：
如图1所示，一种串行高效通信时域O&M定时同步方法，具体实现步骤如下：
步骤1，使用与FPGA基带板卡相连的ADC模数转换器完成信号采样，采样处理后的信号先进行数据捕获同步，捕获后的信号确定了帧头，提取复数导频序列并将其输出到定时同步模块进行定时误差估计。
步骤2，接收信号r(t)经过接收端N倍上采样后得到的样本为：
Figure BDA0003959091350000071
其中，an为发送数据符号，g(t)为成形滤波器、信道和接收滤波器的总响应，ε为定时误差。信号采样捕获后，导频序列用正交复矢量可以表示为rk＝Ik+jQk，其中k＝0,1,2,...,N，取N＝56，提取部分导频序列用于定时误差估计。
如图2所示，复数导频序列分别通过I、Q两路串行输入定时模块，各路以二进制符号形式按位相乘并累计求和得到其平方值，将I路、Q路得到的结果求和即为接收信号模方，其包络值可用
Figure BDA0003959091350000072
表示。
步骤3，将得到的信号包络值分别送入I、Q两路，每路信号分别与相应的旋转因子WnI和Wnq相乘并求和，完成DFT运算。
本方法采样速率满足奈奎斯特采样定理，信号为复信号形式，其中符号速率fB＝4.2GHz，信号带宽B＝2.1GHz，采样速率fs＝4.8GHz，接收端采样倍数为8/7倍采样，则其DFT变换可表示为：
Figure BDA0003959091350000081
其中WnI＝Im[e-j2kπ/N]和Wnq＝Re[e-j2kπ/N]分别为旋转因子，N为DFT点数。
步骤4，对信号模方xk求DFT运算得到Xm的实部和虚部，Xm的实部对应采样点模方与旋转因子WnI的乘积，Xm的虚部对应采样点模方与旋转因子Wnq的乘积，符号累积长度为L＝56，仿真中FFT点数为64点，其结果如图3功率谱所示；I、Q各路串行处理，则Xm的实部虚部分别表示为
Figure BDA0003959091350000082
Figure BDA0003959091350000083
至此，DFT运算可由其接收信号的实部虚部信号处理结果表示，计算可得L个符号的定时误差估计值为
Figure BDA0003959091350000084
计算出频域相位偏移量，得到定时误差估计结果。
步骤5，依据定时误差估计结果进行环路滤波平滑并计算得到定时同步时域延迟，定时误差补偿的滤波器系数抽头位置计算公式为
Figure BDA0003959091350000085
其中
Figure BDA0003959091350000086
为定时误差估计值，fus为根生余弦滤波器符号上采样倍数，计算得到的结果对应ROM存储地址，索引输出的滤波器系数用于后续补偿处理。
步骤6，不同定时误差估计值对应的滤波器系数存储在ROM中，该系数为个数N＝64，8位位宽的二进制数。将步骤5输出的滤波器抽头位置作为ROM索引地址，动态地控制ROM输出的滤波器系数。
步骤7，将索引ROM得到的滤波器系数输入多相滤波器模块，与接收信号采样结果在时域上进行卷积，在完成采样信号定时误差补偿的同时实现数据速率转换，定时误差补偿前后结果如图4所示。
下面为本方法的仿真结果：
1.条件
信号调制方式：QPSK
FFT点数：64
根升余弦滤波器滚降因子：0.25
根生余弦滤波符号上采样倍数：112
2.结果分析
本方法在上述条件下完成定时同步误差估计与补偿，结果如图4所示。左图为定时同步误差补偿前的星座图，右边为定时同步误差补偿后的星座图。可以明显看出，本方法在一定程度上降低了信号传输过程产生的误差，信号采样效果有了一定提升。
总之，本发明通过对高速通信传输数据进行定时误差估计，计算相位延迟误差，将其送入滤波器，从而完成时域定时误差的补偿。针对传统O&M定时同步算法在高速率大带宽传输***中硬件实现采样受限的问题，本发明采用开环IQ两路并行和符号非整数倍采样设计，IQ各路串行处理样本数据，并在时域上进行定时误差补偿，避免DFT和IDFT运算的同时进一步降低了硬件资源的占用。
上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种串行高效通信时域O&M定时同步方法，其特征在于，包括以下步骤：
步骤1，对信号进行采样，完成采样后数据的帧同步，对帧同步后的复数序列进行定时同步处理，得到复数定时同步序列；
步骤2，对复数定时同步序列进行I、Q两路串行输入定时处理，每路均以二进制符号形式按位相乘并累计求和，得到平方值，然后将I、Q两路得到的结果求和，得到接收信号模方；
步骤3，将信号包络值分别送入I、Q两路，每路信号分别与相应的旋转因子WnI和Wnq相乘并求和，完成DFT运算；
步骤4，利用信号DFT运算后的实部和虚部，进行定时误差估计，计算出频域相位偏移量，即定时误差估计结果；
步骤5，依据定时误差估计结果进行环路滤波平滑，计算出定时同步时域延迟，对应到滤波器抽头位置用于后续处理；
步骤6，不同定时误差估计结果对应的滤波器系数存储在ROM中，将步骤5输出的滤波器抽头位置作为ROM索引地址，动态控制ROM输出滤波器系数；
步骤7，将ROM输出的滤波器系数输入到多相滤波器模块，完成采样信号定时误差补偿。
2.根据权利要求1所述的一种串行高效通信时域O&M定时同步方法，其特征在于，步骤1的具体方式为，使用与FPGA基带板卡相连的模数转换模块对信号进行采样，采样后的信号在基带板卡里完成帧同步，输出结果中包括从接收采样数据中提取出的导频序列，将该序列用于定时误差估计。
3.根据权利要求2所述的一种串行高效通信时域O&M定时同步方法，其特征在于，步骤3的DFT运算中，采用满足奈奎斯特采样定理的采样速率，信号为复数形式，DFT运算表示为：
Figure FDA0003959091340000021
其中，xk为接收信号模方；Xm为xk的傅里叶变换，k为样本序号，L为符号累积长度，m为样本段号；I、Q路的傅里叶变换旋转因子分别表示为WnI＝Im[e-j2kπ/N]和Wnq＝Re[e-j2kπ/N]；N为DFT点数。
4.根据权利要求3所述的一种串行高效通信时域O&M定时同步方法，其特征在于，步骤4中，Xm的实部对应采样点模方与旋转因子WnI的乘积，Xm的虚部对应采样点模方与旋转因子Wnq的乘积，符号累积长度为L，I、Q两路分别进行串行处理；Xm的实部、虚部分别表示为：
Figure FDA0003959091340000022
Figure FDA0003959091340000023
其中，Re[Xm]和Im[Xm]分别表示Xm的实部和虚部，Ik和Qk分别为接收样本信号的实部和虚部；
DFT运算由其接收信号的实部、虚部信号运算结果表示，L个符号的定时误差估计值为
Figure FDA0003959091340000024
其中，误差估计值
Figure FDA0003959091340000025
为ε的无偏估计。
5.根据权利要求4所述的一种串行高效通信时域O&M定时同步方法，其特征在于，步骤5中，用于定时误差补偿的滤波器系数抽头位置的计算方式为
Figure FDA0003959091340000031
其中，
Figure FDA0003959091340000032
为定时误差估计值，fus为根生余弦滤波器符号上采样倍数。
6.根据权利要求5所述的一种串行高效通信时域O&M定时同步方法，其特征在于，步骤7的具体方式为，将滤波器系数与信号样点在时域上进行卷积，滤波器采用多相结构设计，利用多相滤波器抽取和插值倍数的不同实现信号定时误差补偿，同时完成数据速率转换。
CN202211474592.1A 2022-11-23 2022-11-23 一种串行高效通信时域o&m定时同步方法 Pending CN115996162A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211474592.1A CN115996162A (zh) 2022-11-23 2022-11-23 一种串行高效通信时域o&m定时同步方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211474592.1A CN115996162A (zh) 2022-11-23 2022-11-23 一种串行高效通信时域o&m定时同步方法

Publications (1)

Publication Number Publication Date
CN115996162A true CN115996162A (zh) 2023-04-21

Family

ID=85989597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211474592.1A Pending CN115996162A (zh) 2022-11-23 2022-11-23 一种串行高效通信时域o&m定时同步方法

Country Status (1)

Country Link
CN (1) CN115996162A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117792338A (zh) * 2024-02-27 2024-03-29 南京朗立微集成电路有限公司 一种滤波器及其设计方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117792338A (zh) * 2024-02-27 2024-03-29 南京朗立微集成电路有限公司 一种滤波器及其设计方法

Similar Documents

Publication Publication Date Title
CN111212007A (zh) 一种通用600Mbps中速解调器实现方法及调制解调器
CN111490866B (zh) 一种高速并行解调结构中的定时恢复方法、数字通信***
CN109274623B (zh) 一种基于大点数fft修正载波频偏的实现方法
CN110300079B (zh) 一种msk信号相干解调方法及***
CN101136890A (zh) 一种优化的多载波信号削波装置及其方法
CN111194077B (zh) 一种低采样率下的定时同步方法
CN113037671B (zh) 一种低复杂度的高效soqpsk符号定时与相位联合同步算法
CN101237435A (zh) 一种降低多载波信号峰均比的方法和装置
CN115996162A (zh) 一种串行高效通信时域o&m定时同步方法
CN107682294B (zh) 一种基于FPGA的高速16apsk信号的相位模糊校正方法
CN110752870A (zh) 滚降系数可变宽带卫星传输***的定时恢复方法及装置
CN107135175A (zh) 一种高速数传***中同步均衡方法
CN101483627B (zh) 减小多载波信号峰均比的方法
CN114938255A (zh) 一种基于o&m算法的高速时域并行定时同步方法
TWI603603B (zh) In the downlink transmission system to solve the imperfect joint estimation of pre-compensation method
CN1327621C (zh) 一种新的超宽带信号的接收方法
CN102801674A (zh) 连续相位调制信号的非相干检测方法及装置
CN113923083A (zh) 一种基于伪随机导频的等效时间采样太赫兹信道估计方法
CN107332632B (zh) MSK调制方式下Alpha稳定分布噪声参数的估计方法
CN110324269A (zh) 一种高速解调器的符号同步***及实现方法
CN107612860A (zh) 适用于802.11ac接收机的同步与下采样估计方法
CN115694549A (zh) 用于uqpsk-dsss信号的定时同步方法及***
CN109714144B (zh) 一种定时同步恢复方法及***
CN114337979B (zh) 一种基于aprx架构的前馈型频域定时同步方法
CN114338304B (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