CN111131106B - Frequency offset estimation method, system, storage medium and receiving device of communication signal - Google Patents
Frequency offset estimation method, system, storage medium and receiving device of communication signal Download PDFInfo
- Publication number
- CN111131106B CN111131106B CN201811287703.1A CN201811287703A CN111131106B CN 111131106 B CN111131106 B CN 111131106B CN 201811287703 A CN201811287703 A CN 201811287703A CN 111131106 B CN111131106 B CN 111131106B
- Authority
- CN
- China
- Prior art keywords
- frequency offset
- communication signal
- phase difference
- offset estimation
- signal
- 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
- 238000004891 communication Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000005540 biological transmission Effects 0.000 claims abstract description 22
- 238000004590 computer program Methods 0.000 claims description 11
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 12
- 230000010363 phase shift Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 239000000969 carrier Substances 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
Images
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/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
-
- 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
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明提供一种通信信号的频偏估计方法、系统、存储介质及接收装置,包括以下步骤:接收通信信号发射端发送的N次长度为M的发射信号,所述发射信号经过信道传输后成为长度为M*N的接收信号;将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r;对所述N个频偏估计中间量r进行不同组合的共轭相乘运算,得到N*(N‑1)/2个相位差估计量;对所述N*(N‑1)/2个相位差估计量求和求平均得到相位差估计结果;对所述相位差估计结果进行计算得到频偏估计值。本发明的一种通信信号的频偏估计方法、系统、存储介质及接收装置,用于计算频偏估计值。
The present invention provides a method, system, storage medium and receiving device for estimating frequency offset of a communication signal, comprising the following steps: receiving a transmission signal of length M sent by a transmitting end of a communication signal N times, and the transmission signal becomes a transmission signal after channel transmission. A received signal with a length of M*N; a cross-correlation operation is performed on the received signal and the transmitted signal to obtain N frequency offset estimated intermediate quantities r; the N frequency offset estimated intermediate quantities r are combined in different combinations. yoke multiplication operation to obtain N*(N‑1)/2 phase difference estimators; sum and average the N*(N‑1)/2 phase difference estimators to obtain the phase difference estimation result; The estimated frequency offset is obtained by calculating the phase difference estimation result. A frequency offset estimation method, system, storage medium and receiving device of a communication signal of the present invention are used for calculating a frequency offset estimation value.
Description
技术领域technical field
本发明涉及通信技术领域,特别是涉及一种通信信号的频偏估计方法、系统、存储介质及接收装置。The present invention relates to the field of communication technologies, and in particular, to a frequency offset estimation method, system, storage medium and receiving device of a communication signal.
背景技术Background technique
在无线通信系统中,频率同步是数据正确传输的前提和基础。一般情况下接收信号的频率偏移主要由载波频偏和多普勒频偏两部分组成。载波频偏是因为接收机晶振和发射机晶振频率不同所造成的,而多普勒频偏是由发射机与接收机之间的相对移动形成的,这些频率偏移会使接收信号附加一个时变的相位旋转,进而导致接受信噪比的下降,影响通信质量。为了获得较高的通信质量,需要对频偏进行估计,并通过自动频率控制进行频率跟踪和补偿。In a wireless communication system, frequency synchronization is the premise and basis for correct data transmission. In general, the frequency offset of the received signal is mainly composed of carrier frequency offset and Doppler frequency offset. The carrier frequency offset is caused by the different frequencies of the receiver crystal oscillator and the transmitter crystal oscillator, and the Doppler frequency offset is formed by the relative movement between the transmitter and the receiver. These frequency offsets will add a time to the received signal. The variable phase rotation, which in turn leads to a decrease in the received signal-to-noise ratio, affects the communication quality. In order to obtain higher communication quality, it is necessary to estimate the frequency offset, and perform frequency tracking and compensation through automatic frequency control.
普遍的频偏估计方法都是通过插入在同步信号的前导、导频或循环前缀进行估计的。发射机在同步信号中插入用于频偏估计的训练序列,当接收机收到该同步信号时使用自相关法来提取前后两个相同序列的相位差,由此进行频偏估计。Common frequency offset estimation methods are estimated by inserting the preamble, pilot or cyclic prefix in the synchronization signal. The transmitter inserts a training sequence for frequency offset estimation into the synchronization signal, and when the receiver receives the synchronization signal, the autocorrelation method is used to extract the phase difference between the two identical sequences before and after, thereby performing frequency offset estimation.
在实现本发明的过程中,发明人发现现在有技术至少存在以下问题。In the process of realizing the present invention, the inventor found that the existing technology has at least the following problems.
在通信条件较差,信噪比较低的情况下,由于同步信号被噪声所淹没,对用于提取频偏的统计估计量影响很大,导致估计不准确,均方误差较大。使用估计结果对接收信号进行频偏补偿不仅不会产生任何增益,反而会使得接受性能下降。In the case of poor communication conditions and low signal-to-noise ratio, since the synchronization signal is submerged by noise, it has a great influence on the statistical estimator used to extract the frequency offset, resulting in inaccurate estimation and large mean square error. Using the estimation result to perform frequency offset compensation on the received signal will not produce any gain, but will degrade the receiving performance.
因此,希望能够解决如何更好地进行频偏估计,以避免因使用自相关方法计算得到的估计结果对接收信号进行频偏补偿而使性能下降的问题。Therefore, it is desirable to solve the problem of how to better estimate the frequency offset, so as to avoid the problem of performance degradation due to the use of the estimation result calculated by the autocorrelation method to perform frequency offset compensation on the received signal.
发明内容SUMMARY OF THE INVENTION
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种通信信号的频偏估计方法、系统、存储介质及接收装置,用于解决现有技术中如何更好地进行频偏估计,以避免因使用自相关方法计算得到的估计结果对接收信号进行频偏补偿而使性能下降的问题。In view of the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide a frequency offset estimation method, system, storage medium and receiving device for a communication signal, which are used to solve how to better perform frequency offset estimation in the prior art, In order to avoid the problem of performance degradation due to the use of the estimation result calculated by the autocorrelation method to perform frequency offset compensation on the received signal.
为实现上述目的及其他相关目的,本发明提供一种通信信号的频偏估计方法,包括以下步骤:接收通信信号发射端发送的N个长度为M的发射信号,所述发射信号经过信道传输后成为长度为M*N的接收信号;将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r;对所述N个频偏估计中间量r进行不同组合的共轭相乘运算,得到N*(N-1)/2个相位差估计量;对所述N*(N-1)/2个相位差估计量求和求平均得到相位差估计结果;对所述相位差估计结果进行计算得到频偏估计值。In order to achieve the above purpose and other related purposes, the present invention provides a method for estimating frequency offset of a communication signal, comprising the following steps: receiving N transmission signals with a length of M sent by a transmitting end of the communication signal, and after the transmission signal is transmitted through a channel become a received signal with a length of M*N; perform a cross-correlation operation on the received signal and the transmitted signal to obtain N frequency offset estimated intermediate quantities r; perform different combinations of the N frequency offset estimated intermediate quantities r Conjugate multiplication operation to obtain N*(N-1)/2 phase difference estimators; sum and average the N*(N-1)/2 phase difference estimators to obtain the phase difference estimation result; The phase difference estimation result is calculated to obtain a frequency offset estimation value.
于本发明的一实施例中,所述个数为M的发射信号表示为所述发射信号的功率已进行归一化为 In an embodiment of the present invention, the number of M transmission signals is expressed as The power of the transmitted signal has been normalized to
于本发明的一实施例中,所述长度为M*N的接收信号表示为所述N个频偏估计中间量r根据计算,1≤j≤M,1≤i≤N。In an embodiment of the present invention, the received signal whose length is M*N is expressed as The N frequency offset estimation intermediate quantities r are based on Calculation, 1≤j≤M, 1≤i≤N.
于本发明的一实施例中,所述N*(N-1)/2个相位差估计量为 In an embodiment of the present invention, the N*(N-1)/2 phase difference estimators are
于本发明的一实施例中,还包括根据计算频偏估计值fd,fs为采样率,θ为相位差估计结果。In an embodiment of the present invention, it also includes Calculate the frequency offset estimation value f d , where f s is the sampling rate, and θ is the phase difference estimation result.
为实现上述目的,本发明还提供一种通信信号的频偏估计系统,包括:接收模块、第一运算模块、第二运算模块、第三运算模块和第四运算模块;所述接收模块用于接收通信信号发射端发送的N次个数为M的发射信号,所述发射信号经过信道传输后成为M*N个接收信号;所述第一运算模块用于将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r;所述第二运算模块用于对所述N个频偏估计中间量r进行不同组合的共轭相乘运算,得到N*(N-1)/2个相位差估计量;所述第三运算模块用于对所述N*(N-1)/2个相位差估计量求和求平均得到相位差估计结果;所述第四运算模块用于对所述相位差估计结果进行计算得到频偏估计值。In order to achieve the above purpose, the present invention also provides a frequency offset estimation system for communication signals, comprising: a receiving module, a first computing module, a second computing module, a third computing module and a fourth computing module; the receiving module is used for Receive the N times the number of M transmission signals sent by the transmitting end of the communication signal, the transmission signals become M*N received signals after channel transmission; the first arithmetic module is used to combine the received signals with the transmitted signals The signal is subjected to cross-correlation operation to obtain N frequency offset estimated intermediate quantities r; the second operation module is used to perform different combinations of conjugate multiplication operations on the N frequency offset estimated intermediate quantities r to obtain N*(N -1)/2 phase difference estimators; the third operation module is used to sum and average the N*(N-1)/2 phase difference estimators to obtain a phase difference estimation result; the fourth The operation module is used for calculating the phase difference estimation result to obtain the frequency offset estimation value.
为实现上述目的,本发明还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现任一上述通信信号的频偏估计方法。In order to achieve the above object, the present invention also provides a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, implements any of the above methods for estimating frequency offset of a communication signal.
为实现上述目的,本发明还提供一种接收装置,包括:处理器和存储器;所述存储器用于存储计算机程序;所述处理器与所述存储器相连,用于执行所述存储器存储的计算机程序,以使所述接收装置执行任一上述的通信信号的频偏估计方法。In order to achieve the above object, the present invention also provides a receiving device, comprising: a processor and a memory; the memory is used for storing a computer program; the processor is connected to the memory and is used for executing the computer program stored in the memory , so that the receiving apparatus executes any one of the above methods for estimating frequency offsets of communication signals.
最后,本发明还提供一种通信信号的频偏估计系统,包括:包括上述的通信信号接收装置和通信信号发射装置;所述通信信号发射装置用于发送N个长度为M的发射信号。Finally, the present invention also provides a frequency offset estimation system for a communication signal, including the above-mentioned communication signal receiving device and communication signal transmitting device; the communication signal transmitting device is used for transmitting N transmission signals of length M.
于本发明的一实施例中,所述通信信号接收装置为天线、放大器、滤波器、混频器、模数转换器或处理器。In an embodiment of the present invention, the communication signal receiving device is an antenna, an amplifier, a filter, a mixer, an analog-to-digital converter or a processor.
如上所述,本发明的一种通信信号的频偏估计方法、系统、存储介质及接收装置,具有以下有益效果:基于频偏互相关的方法计算得到频偏估计值;所述频偏估计值应用于接收信号频偏补偿时能够提升性能、增加增益。另外,当信号条件较好的时候,还可以使用简化算法来更快的估计频偏,简化算法可以减少计算量,提升同步速度。As described above, the method, system, storage medium and receiving device for frequency offset estimation of a communication signal of the present invention have the following beneficial effects: a frequency offset estimated value is obtained by calculating a frequency offset cross-correlation method; the frequency offset estimated value When applied to the frequency offset compensation of the received signal, it can improve the performance and increase the gain. In addition, when the signal conditions are good, a simplified algorithm can also be used to estimate the frequency offset faster. The simplified algorithm can reduce the amount of calculation and improve the synchronization speed.
附图说明Description of drawings
图1显示为本发明的通信信号的频偏估计方法于一实施例中的流程图;FIG. 1 is a flowchart of a method for estimating frequency offset of a communication signal according to an embodiment of the present invention;
图2显示为本发明的发射信号与接收信号长度关系于一实施例中的示意图;FIG. 2 is a schematic diagram showing the relationship between the length of the transmitted signal and the received signal according to an embodiment of the present invention;
图3显示为本发明的将所述接收信号与所述发射信号进行互相关运算得到N个频偏估计中间量r于一实施例中的示意图;FIG. 3 is a schematic diagram showing N frequency offset estimated intermediate quantities r obtained by performing a cross-correlation operation on the received signal and the transmitted signal according to an embodiment of the present invention;
图4显示为本发明的对所述N个频偏估计中间量r进行不同组合的共轭相乘运算得到N*(N-1)/2个相位差估计量于一实施例中的示意图;4 is a schematic diagram showing N*(N-1)/2 phase difference estimators obtained by performing different combinations of conjugate multiplication operations on the N frequency offset estimation intermediate quantities r according to an embodiment of the present invention;
图5显示为本发明的1kHz频偏互相关方法残余频偏统计于一实施例中的示意图;5 is a schematic diagram of residual frequency offset statistics of the 1kHz frequency offset cross-correlation method of the present invention in an embodiment;
图6显示为本发明的1kHz频偏自相关方法残余频偏统计于一实施例中的示意图;FIG. 6 is a schematic diagram showing the residual frequency offset statistics of the 1kHz frequency offset autocorrelation method of the present invention in an embodiment;
图7显示为本发明的2kHz频偏互相关方法残余频偏统计于一实施例中的示意图;FIG. 7 is a schematic diagram showing the residual frequency offset statistics of the 2kHz frequency offset cross-correlation method of the present invention in an embodiment;
图8显示为本发明的2kHz频偏自相关方法残余频偏统计于一实施例中的示意图;FIG. 8 is a schematic diagram showing the residual frequency offset statistics of the 2kHz frequency offset autocorrelation method of the present invention in an embodiment;
图9显示为本发明的1kHz频偏互相关方法残余频偏统计于一实施例中的示意图;FIG. 9 is a schematic diagram showing the residual frequency offset statistics of the 1kHz frequency offset cross-correlation method of the present invention in an embodiment;
图10显示为本发明的1kHz频偏自相关方法残余频偏统计于一实施例中的示意图;FIG. 10 is a schematic diagram showing the residual frequency offset statistics of the 1kHz frequency offset autocorrelation method of the present invention in an embodiment;
图11显示为本发明的通信信号的频偏估计系统于一实施例中的结构示意图;FIG. 11 is a schematic structural diagram of a frequency offset estimation system for a communication signal according to an embodiment of the present invention;
图12显示为本发明的通信信号接收装置于一实施例中的结构示意图;FIG. 12 is a schematic structural diagram of a communication signal receiving apparatus according to an embodiment of the present invention;
图13显示为本发明的通信信号的频偏估计系统于又一实施例中的结构示意图。FIG. 13 is a schematic structural diagram of a frequency offset estimation system for a communication signal according to another embodiment of the present invention.
元件标号说明Component label description
111 接收模块111 Receiver module
112 第一运算模块112 The first operation module
113 第二运算模块113 Second operation module
114 第三运算模块114 The third operation module
115 第四运算模块115 Fourth arithmetic module
121 处理器121 processors
122 存储器122 memory
131 通信信号接收装置131 Communication signal receiving device
132 通信信号发射装置132 Communication signal transmitter
具体实施方式Detailed ways
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。The embodiments of the present invention are described below through specific specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments may be combined with each other under the condition of no conflict.
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。It should be noted that the drawings provided in the following embodiments are only used to illustrate the basic concept of the present invention in a schematic way, so the drawings only show the components related to the present invention rather than the number, shape and number of components in actual implementation. For dimension drawing, the type, quantity and proportion of each component can be changed at will in actual implementation, and the component layout may also be more complicated.
本发明的通信信号的频偏估计方法、系统、存储介质及接收装置,基于频偏互相关的方法计算得到频偏估计值;所述频偏估计值应用于接收信号频偏补偿时能够提升性能、增加增益。The frequency offset estimation method, system, storage medium and receiving device of the communication signal of the present invention calculate and obtain the frequency offset estimation value based on the frequency offset cross-correlation method; the frequency offset estimation value can improve the performance when applied to the frequency offset compensation of the received signal , increase the gain.
如图1所示,于一实施例中,本发明的通信信号的频偏估计方法,包括以下步骤:As shown in FIG. 1 , in an embodiment, the method for estimating frequency offset of a communication signal of the present invention includes the following steps:
步骤S11、接收通信信号发射端发送的N次个数为M的发射信号,所述发射信号经过信道传输后成为M*N个接收信号。Step S11: Receive N times of M transmission signals sent by the transmitting end of the communication signal, and the transmission signals become M*N received signals after channel transmission.
于本发明一实施例中,所述N次个数为M的发射信号为已归一化的圆周调制复信号,所述长度为M归一化的圆周调制复信号可表示为所述发射信号包括但不限于2PSK(2PSK是相移键控的最简单的一种形式,它用两个初相相隔为180的载波来传递二进制信息,所以也被称为BPSK)、QPSK(正交相移键控(Quadrature Phase ShiftKeyin,QPSK)是一种数字调制方式。它分为绝对相移和相对相移两种)等。所述发射信号是已归一的信号。具体地,Xk中的k可表示为M*N,如k=2*3时则M=2,N=3,表示为第3次发送的复信号中的x2,所述M*N个接收信号表示为 In an embodiment of the present invention, the N times the transmit signal with the number M is a normalized circularly modulated complex signal, and the normalized circularly modulated complex signal with a length of M can be expressed as The transmission signal includes but is not limited to 2PSK (2PSK is the simplest form of phase shift keying, it uses two carriers with an initial phase interval of 180 to transmit binary information, so it is also called BPSK), QPSK ( Quadrature Phase Shift Keying (Quadrature Phase Shift Keyin, QPSK) is a digital modulation method. It is divided into absolute phase shift and relative phase shift) and so on. The transmitted signal is a normalized signal. Specifically, k in X k can be expressed as M*N, for example, when k=2*3, M=2, N=3, which is expressed as the complex signal sent for the third time x 2 in , the M*N received signals are expressed as
如图2所示,于一实施例中,发送长度为M的发射信号N次,得到经过信道传输后长度为M*N的接收信号。As shown in FIG. 2 , in an embodiment, a transmit signal with a length of M is sent N times to obtain a received signal with a length of M*N after channel transmission.
步骤S12、将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r。Step S12: Perform a cross-correlation operation on the received signal and the transmitted signal to obtain N intermediate frequency offset estimates r.
于本发明一实施例中,所述长度为M*N的接收信号表示为所述N个频偏估计中间量r根据计算,1≤j≤M,1≤i≤N。所述xj *与xj为复共轭关系。In an embodiment of the present invention, the received signal whose length is M*N is expressed as The N frequency offset estimation intermediate quantities r are based on Calculation, 1≤j≤M, 1≤i≤N. The x j * and x j are in a complex conjugate relationship.
具体地,如图3所示, Specifically, as shown in Figure 3,
步骤S13、对所述N个频偏估计中间量r进行不同组合的共轭相乘运算,得到N*(N-1)/2个相位差估计量。Step S13: Perform different combinations of conjugate multiplication operations on the N intermediate frequency offset estimates r to obtain N*(N-1)/2 phase difference estimates.
于本发明一实施例中,所述N*(N-1)/2个相位差估计量为 其中所述rj *与rj为复共轭关系。In an embodiment of the present invention, the N*(N-1)/2 phase difference estimators are The r j * and r j are in a complex conjugate relationship.
具体地,如图4所示, Specifically, as shown in Figure 4,
一般情况下S13步骤需要计算N*(N-1)/2个ri。当SNR较高的时候(相比信号质量更差的情况)简化算法可以减少计算的ri的个数,比如可以只计算2个,所以可以减少计算量,提升同步速度。具体地,当SNR较高的时候所述简化算法可以在步骤S12将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r时只计算一个或两个r值,例如只计算r1和r2。后续步骤S13只需对r1和r2进行不同组合的共轭相乘运算得到一个相位差估计量。In general, step S13 needs to calculate N*(N-1)/2 r i . When the SNR is high (compared to the case of worse signal quality), the simplified algorithm can reduce the number of ri to be calculated, for example, only 2 ri can be calculated, so the amount of calculation can be reduced and the synchronization speed can be improved. Specifically, when the SNR is high, the simplified algorithm may perform a cross-correlation operation on the received signal and the transmitted signal in step S12, and only calculate one or two r values when N frequency offset estimation intermediate quantities r are obtained , for example to compute only r 1 and r 2 . Subsequent step S13 only needs to perform different combinations of conjugate multiplication operations on r 1 and r 2 to obtain a phase difference estimator.
步骤S14、对所述N*(N-1)/2个相位差估计量求和求平均得到相位差估计结果。Step S14, summing and averaging the N*(N-1)/2 phase difference estimates to obtain a phase difference estimation result.
于本发明一实施例中,对所述N*(N-1)/2个相位差估计量为求和,i<j,1≤j≤N,1≤i≤N,得到相位差估计结果。所述相位差估计结果即为所有的θi,j的求和平均值。In an embodiment of the present invention, the N*(N-1)/2 phase difference estimators are Sum, i<j, 1≤j≤N, 1≤i≤N, get the phase difference estimation result. The phase difference estimation result That is, the summed average of all θ i,j .
步骤S15、对所述相位差估计结果进行计算得到频偏估计值。Step S15: Calculate the phase difference estimation result to obtain a frequency offset estimation value.
于本发明一实施例中,还包括根据计算频偏估计值fd,fs为采样率,θ为相位差估计结果。In an embodiment of the present invention, it also includes Calculate the frequency offset estimation value f d , where f s is the sampling rate, and θ is the phase difference estimation result.
如图5所示,于一实施例中,1kHz频偏互相关方法残余频偏统图显示;如图6所示1kHz频偏自相关方法残余频偏统计图显示。由图可见,两种方法均能够估计出频偏。本发明采用的方法得到的残余频偏更小。As shown in FIG. 5 , in one embodiment, the residual frequency offset chart of the 1kHz frequency offset cross-correlation method is displayed; as shown in FIG. 6 , the residual frequency offset chart of the 1kHz frequency offset autocorrelation method is displayed. It can be seen from the figure that both methods can estimate the frequency offset. The residual frequency offset obtained by the method adopted in the present invention is smaller.
如图7所示,于一实施例中,2kHz频偏互相关方法残余频偏统计图显示;如图8所示2kHz频偏自相关方法残余频偏统计图。由图可见,自相关法已经无法估计出频偏,本发明的方法可以估计出频偏,并使得残余频偏下降。As shown in FIG. 7 , in an embodiment, the residual frequency offset statistics of the 2kHz frequency offset cross-correlation method are displayed; as shown in FIG. 8 , the residual frequency offset statistics of the 2kHz frequency offset autocorrelation method are displayed. It can be seen from the figure that the autocorrelation method has been unable to estimate the frequency offset, and the method of the present invention can estimate the frequency offset and reduce the residual frequency offset.
如图9所示,于一实施例中,1kHz频偏互相关方法残余频偏统计图显示;-如图10所示1kHz频偏自相关方法残余频偏统计图。由图可见,自相关法已经无法估计出频偏,本发明的方法可以估计出频偏,并使得残余频偏下降。图5,6与图9,10的区别在于SNR(信噪比)值不一样。图5,6的SNR为-10dB,图9,10的SNR为-20dB。本发明的技术效果是与传统的采用自相关方法的通信设备相比,在理论上使用互相关方法的通信设备的估计方差更小;在相同长度的训练符号情况下,理论方差缩小为原有的SNR:信噪比SIGNAL NOISERATIO。另外差分方式的互相关方法可以抵消信号的初始相位,提高估计精度。As shown in FIG. 9 , in one embodiment, the residual frequency offset statistics of the 1kHz frequency offset cross-correlation method is displayed; - the residual frequency offset statistics diagram of the 1kHz frequency offset autocorrelation method is shown in FIG. 10 . It can be seen from the figure that the autocorrelation method has been unable to estimate the frequency offset, and the method of the present invention can estimate the frequency offset and reduce the residual frequency offset. The difference between Figures 5 and 6 and Figures 9 and 10 is that the SNR (Signal to Noise Ratio) values are different. The SNR of Figures 5 and 6 is -10dB, and the SNR of Figures 9 and 10 is -20dB. The technical effect of the invention is that compared with the traditional communication equipment using the autocorrelation method, the estimated variance of the communication equipment using the cross-correlation method is theoretically smaller; in the case of training symbols of the same length, the theoretical variance is reduced to the original of SNR: Signal-to-noise ratio SIGNAL NOISERATIO. In addition, the cross-correlation method of the differential mode can cancel the initial phase of the signal and improve the estimation accuracy.
如图11所示,于一实施例中,本发明的通信信号的频偏估计系统,包括接收模块111、第一运算模块112、第二运算模块113、第三运算模块114和第四运算模块115。As shown in FIG. 11 , in an embodiment, the frequency offset estimation system of the communication signal of the present invention includes a receiving
所述接收模块111用于接收通信信号发射端发送的N个长度为M的发射信号,所述发射信号经过信道传输后成长度为M*N的接收信号。The receiving
于本发明一实施例中,所述N个长度为M的发射信号为已归一化的圆周调制复信号,所述长度为M归一化的圆周调制复信号可表示为所述发射信号包括但不限于QPSK(正交相移键控(Quadrature Phase Shift Keyin,QPSK)是一种数字调制方式。它分为绝对相移和相对相移两种)、2PSK(2PSK是相移键控的最简单的一种形式,它用两个初相相隔为180的载波来传递二进制信息,所以也被称为BPSK)等。具体地,Xk中的k可表示为M*N,如k=2*3时则M=2,N=3,表示为第3次发送的复信号中的x2,所述长度为M*N的接收信号表示为 In an embodiment of the present invention, the N transmit signals with a length of M are normalized circularly modulated complex signals, and the normalized circularly modulated complex signals with a length of M can be expressed as The transmission signal includes but is not limited to QPSK (Quadrature Phase Shift Keyin (QPSK) is a digital modulation method. It is divided into absolute phase shift and relative phase shift), 2PSK (2PSK is a phase shift The simplest form of shift keying, it uses two carriers with an initial phase interval of 180 to transmit binary information, so it is also called BPSK) and so on. Specifically, k in X k can be expressed as M*N, for example, when k=2*3, M=2, N=3, which is expressed as the complex signal sent for the third time x 2 in , the received signal of length M*N is expressed as
于一实施例中,发送长度为M的发射信号N次,得到经过信道传输后长度为M*N的接收信号。In an embodiment, a transmit signal with a length of M is sent N times, and a received signal with a length of M*N is obtained after channel transmission.
所述第一运算模块112用于将所述接收信号与所述发射信号进行互相关运算,得到N个频偏估计中间量r。The
于本发明一实施例中,所述长度为M*的个接收信号表示为所述N个频偏估计中间量r根据计算,1≤j≤M,1≤i≤N。所述xj *与xj为复共轭关系。In an embodiment of the present invention, the received signals whose length is M* are expressed as The N frequency offset estimation intermediate quantities r are based on Calculation, 1≤j≤M, 1≤i≤N. The x j * and x j are in a complex conjugate relationship.
具体地, specifically,
所述第二运算模块113用于对所述N个频偏估计中间量r进行不同组合的共轭相乘运算,得到N*(N-1)/2个相位差估计量。The
于本发明一实施例中,所述N*(N-1)/2个相位差估计量为 其中所述rj *与rj为复共轭关系。In an embodiment of the present invention, the N*(N-1)/2 phase difference estimators are The r j * and r j are in a complex conjugate relationship.
具体地, specifically,
所述第三运算模块114用于对所述N*(N-1)/2个相位差估计量求和求平均得到相位差估计结果。The
于本发明一实施例中,对所述N*(N-1)/2个相位差估计量为求和,i<j,1≤j≤N,1≤i≤N,得到相位差估计结果。所述相位差估计结果即为所有的θi,j的求和平均值。In an embodiment of the present invention, the N*(N-1)/2 phase difference estimators are Sum, i<j, 1≤j≤N, 1≤i≤N, get the phase difference estimation result. The phase difference estimation result That is, the summed average of all θ i,j .
所述第四运算模块115用于对所述相位差估计结果进行计算得到频偏估计值。The
于本发明一实施例中,还包括根据计算频偏估计值fd,fs为采样率,θ为相位差估计结果。In an embodiment of the present invention, it also includes Calculate the frequency offset estimation value f d , where f s is the sampling rate, and θ is the phase difference estimation result.
需要说明的是,应理解以上系统的各个模块的划分仅仅是一种逻辑功能的划分,实际实现时可以全部或部分集成到一个物理实体上,也可以物理上分开。且这些模块可以全部以软件通过处理元件调用的形式实现;也可以全部以硬件的形式实现;还可以部分模块通过处理元件调用软件的形式实现,部分模块通过硬件的形式实现。例如,x模块可以为单独设立的处理元件,也可以集成在上述装置的某一个芯片中实现,此外,也可以以程序代码的形式存储于上述装置的存储器中,由上述装置的某一个处理元件调用并执行以上x模块的功能。其它模块的实现与之类似。此外这些模块全部或部分可以集成在一起,也可以独立实现。这里所述的处理元件可以是一种集成电路,具有信号的处理能力。在实现过程中,上述方法的各步骤或以上各个模块可以通过处理器元件中的硬件的集成逻辑电路或者软件形式的指令完成。It should be noted that it should be understood that the division of each module of the above system is only a division of logical functions, and may be fully or partially integrated into a physical entity in actual implementation, or may be physically separated. And these modules can all be implemented in the form of software calling through processing elements; they can also all be implemented in hardware; some modules can also be implemented in the form of calling software through processing elements, and some modules can be implemented in hardware. For example, the x module may be a separately established processing element, or it may be integrated into a certain chip of the above-mentioned device to be implemented, in addition, it may also be stored in the memory of the above-mentioned device in the form of program code, and a certain processing element of the above-mentioned device Calls and executes the functions of the above x module. The implementation of other modules is similar. In addition, all or part of these modules can be integrated together, and can also be implemented independently. The processing element described here may be an integrated circuit with signal processing capability. In the implementation process, each step of the above-mentioned method or each of the above-mentioned modules can be completed by an integrated logic circuit of hardware in the processor element or an instruction in the form of software.
例如,以上这些模块可以是被配置成实施以上方法的一个或多个集成电路,例如:一个或多个特定集成电路(Application Specific Integrated Circuit,简称ASIC),或,一个或多个微处理器(Digital Singnal Processor,简称DSP),或,一个或者多个现场可编程门阵列(Field Programmable Gate Array,简称接收装置)等。再如,当以上某个模块通过处理元件调度程序代码的形式实现时,该处理元件可以是通用处理器,例如中央处理器(Central Processing Unit,简称CPU)或其它可以调用程序代码的处理器。再如,这些模块可以集成在一起,以片上系统(system-on-a-chip,简称SOC)的形式实现。For example, the above modules may be one or more integrated circuits configured to implement the above methods, such as: one or more specific integrated circuits (Application Specific Integrated Circuit, ASIC for short), or one or more microprocessors ( Digital Singnal Processor, DSP for short), or one or more Field Programmable Gate Arrays (Field Programmable Gate Array, receiving device for short), etc. For another example, when one of the above modules is implemented in the form of processing element scheduling program code, the processing element may be a general-purpose processor, such as a central processing unit (Central Processing Unit, CPU for short) or other processors that can call program codes. For another example, these modules can be integrated together and implemented in the form of a system-on-a-chip (SOC for short).
于本发明一实施例中,本发明还包括一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述任一通信信号的频偏估计方法。In an embodiment of the present invention, the present invention further includes a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, implements any of the above-mentioned frequency offset estimation methods for communication signals.
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过计算机程序相关的硬件来完成。前述的计算机程序可以存储于一计算机可读存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。Those of ordinary skill in the art can understand that all or part of the steps of implementing the above method embodiments may be completed by hardware related to computer programs. The aforementioned computer program may be stored in a computer-readable storage medium. When the program is executed, the steps including the above method embodiments are executed; and the foregoing storage medium includes: ROM, RAM, magnetic disk or optical disk and other media that can store program codes.
如图10所示,于一实施例中,本发明的接收装置包括:处理器121和存储器122;所述存储器122用于存储计算机程序;所述处理器121与所述存储器122相连,用于执行所述存储器122存储的计算机程序,以使所述接收装置执行任一所述的通信信号的频偏估计方法。As shown in FIG. 10, in an embodiment, the receiving apparatus of the present invention includes: a
具体地,所述存储器122包括:ROM、RAM、磁碟、U盘、存储卡或者光盘等各种可以存储程序代码的介质。Specifically, the
优选地,所述处理器121可以是通用处理器,包括中央处理器(CentralProcessing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital Signal Processor,简称DSP)、专用集成电路(ApplicationSpecific Integrated Circuit,简称ASIC)、现场可编程门阵列(Field ProgrammableGate Array,简称接收装置)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。Preferably, the
如图11所示,于一实施例中,本发明的通信信号的频偏估计系统,包括上述的通信信号接收装置131和通信信号发射装置132。所述通信信号发射装置132用于发送N次个数为M的发射信号。As shown in FIG. 11 , in an embodiment, the frequency offset estimation system of a communication signal of the present invention includes the above-mentioned communication
于本发明一实施例中,所述通信信号接收装置包括但不限于天线、放大器、滤波器、混频器、模数转换器或处理器。In an embodiment of the present invention, the communication signal receiving device includes, but is not limited to, an antenna, an amplifier, a filter, a mixer, an analog-to-digital converter, or a processor.
综上所述,本发明通信信号的频偏估计方法、系统、存储介质及接收装置,基于频偏互相关的方法计算得到频偏估计值;所述频偏估计值应用于接收信号频偏补偿时能够提升性能、增加增益。所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。To sum up, the frequency offset estimation method, system, storage medium and receiving device of the present invention are based on the frequency offset cross-correlation method to calculate and obtain the frequency offset estimation value; the frequency offset estimation value is applied to the frequency offset compensation of the received signal can improve performance and increase gain. Therefore, the present invention effectively overcomes various shortcomings in the prior art and has high industrial application value.
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。The above-mentioned embodiments merely illustrate the principles and effects of the present invention, but are not intended to limit the present invention. Anyone skilled in the art can make modifications or changes to the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or changes made by those with ordinary knowledge in the technical field without departing from the spirit and technical idea disclosed in the present invention should still be covered by the claims of the present invention.
Claims (7)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811287703.1A CN111131106B (en) | 2018-10-31 | 2018-10-31 | Frequency offset estimation method, system, storage medium and receiving device of communication signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811287703.1A CN111131106B (en) | 2018-10-31 | 2018-10-31 | Frequency offset estimation method, system, storage medium and receiving device of communication signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111131106A CN111131106A (en) | 2020-05-08 |
| CN111131106B true CN111131106B (en) | 2022-08-30 |
Family
ID=70485511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811287703.1A Active CN111131106B (en) | 2018-10-31 | 2018-10-31 | Frequency offset estimation method, system, storage medium and receiving device of communication signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111131106B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115801518B (en) * | 2022-11-30 | 2024-04-09 | 湖南智领通信科技有限公司 | Frequency offset estimation method and device based on probability distribution statistics and computer equipment |
| CN120182615B (en) * | 2025-05-22 | 2025-09-05 | 佛山职业技术学院 | Hyperspectral image noise detection and evaluation method and system |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1314033A (en) * | 1998-06-29 | 2001-09-19 | 艾利森电话股份有限公司 | Method, apparatus and system for determing location of frequency synchronization signal |
| CN101079857A (en) * | 2006-05-25 | 2007-11-28 | 北京泰美世纪科技有限公司 | A carrier residual frequency deviation tracking method based on OFDM system |
| CN101299737A (en) * | 2007-04-30 | 2008-11-05 | 华为技术有限公司 | Synchronization Estimation Method and System of Orthogonal Frequency Division Multiplexing Technology |
| CN101374134A (en) * | 2008-07-03 | 2009-02-25 | 天津大学 | A Precise Frequency Offset Estimation Method for Time Domain Synchronous Orthogonal Frequency Division Multiplexing Receiver |
| US7623467B1 (en) * | 2003-09-17 | 2009-11-24 | Atheros Communications, Inc. | Wireless channel estimation |
| CN101969329A (en) * | 2010-09-30 | 2011-02-09 | 中国科学院国家天文台 | Signal phase difference estimation device and method in antenna aligning array based on circulating cross-correlation |
| CN102137048A (en) * | 2010-01-22 | 2011-07-27 | 中兴通讯股份有限公司 | Method and system for frequency offset estimation |
| CN102170414A (en) * | 2011-05-17 | 2011-08-31 | 浙江瑞讯微电子有限公司 | Demodulation and timing synchronization combined method for GFSK (Gauss Frequency Shift Key) |
| CN103428147A (en) * | 2012-05-24 | 2013-12-04 | 普天信息技术研究院有限公司 | Frequency offset compensation method in TDD-LTE system |
| WO2017067601A1 (en) * | 2015-10-23 | 2017-04-27 | Huawei Technologies Co., Ltd. | Communication device and method for joint offset estimation |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100054211A1 (en) * | 2008-08-27 | 2010-03-04 | Qualcomm Incorporated | Frequency domain pn sequence |
| US10575272B2 (en) * | 2016-02-03 | 2020-02-25 | Lg Electronics Inc. | Method and apparatus for synchronizing frequency and time in a wireless communication system |
-
2018
- 2018-10-31 CN CN201811287703.1A patent/CN111131106B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1314033A (en) * | 1998-06-29 | 2001-09-19 | 艾利森电话股份有限公司 | Method, apparatus and system for determing location of frequency synchronization signal |
| US7623467B1 (en) * | 2003-09-17 | 2009-11-24 | Atheros Communications, Inc. | Wireless channel estimation |
| CN101079857A (en) * | 2006-05-25 | 2007-11-28 | 北京泰美世纪科技有限公司 | A carrier residual frequency deviation tracking method based on OFDM system |
| CN101299737A (en) * | 2007-04-30 | 2008-11-05 | 华为技术有限公司 | Synchronization Estimation Method and System of Orthogonal Frequency Division Multiplexing Technology |
| CN101374134A (en) * | 2008-07-03 | 2009-02-25 | 天津大学 | A Precise Frequency Offset Estimation Method for Time Domain Synchronous Orthogonal Frequency Division Multiplexing Receiver |
| CN102137048A (en) * | 2010-01-22 | 2011-07-27 | 中兴通讯股份有限公司 | Method and system for frequency offset estimation |
| CN101969329A (en) * | 2010-09-30 | 2011-02-09 | 中国科学院国家天文台 | Signal phase difference estimation device and method in antenna aligning array based on circulating cross-correlation |
| CN102170414A (en) * | 2011-05-17 | 2011-08-31 | 浙江瑞讯微电子有限公司 | Demodulation and timing synchronization combined method for GFSK (Gauss Frequency Shift Key) |
| CN103428147A (en) * | 2012-05-24 | 2013-12-04 | 普天信息技术研究院有限公司 | Frequency offset compensation method in TDD-LTE system |
| WO2017067601A1 (en) * | 2015-10-23 | 2017-04-27 | Huawei Technologies Co., Ltd. | Communication device and method for joint offset estimation |
Non-Patent Citations (2)
| Title |
|---|
| "Fast coarse frequency offset estimation for OFDM systems by using differentially modulated subcarriers";Bo-Seok Seo等;《IEEE Transactions on Consumer Electronics》;20021130;第48卷(第4期);全文 * |
| "延迟共轭相乘的弹道目标平动补偿";贺思三等;《雷达学报》;20141113(第5期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111131106A (en) | 2020-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109462566B (en) | An information transmission method and device | |
| CN102368756B (en) | Channel estimation for use on LTE terminals | |
| CN104092635B (en) | Carrier frequency offset estimation method, carrier compensation method and device | |
| US12081377B2 (en) | Demodulating modulated signals | |
| CN111131106B (en) | Frequency offset estimation method, system, storage medium and receiving device of communication signal | |
| CN106059984A (en) | Digital phase-modulated signal carrier phase offset estimation method | |
| US20040120248A1 (en) | Delta-predicted frequency offset compensation apparatus and method thereof | |
| US8139682B2 (en) | Channel estimation of orthogonal frequency division multiplexed systems | |
| CN115941406B (en) | Frequency offset estimation method and device, electronic equipment and readable storage medium | |
| US8842753B2 (en) | Orthogonal frequency division multiplexed (OFDM) demodulator imbalance estimation | |
| CN102118333B (en) | Method and device for canceling DC offsets | |
| US10841139B2 (en) | O-QPSK symbol timing estimation | |
| CN111355539B (en) | Method, terminal and computer-readable storage medium for obtaining channel estimation value | |
| CN114039829B (en) | Phase recovery method and system based on pilot frequency assistance | |
| US9876660B2 (en) | Apparatus and method for estimating carrier frequency offset | |
| CN104539564A (en) | Frequency offset estimation method and device for LTE (Long Term Evolution) system | |
| CN111147417B (en) | Method, device, terminal and storage medium for realizing carrier frequency offset estimation | |
| CN106464644B (en) | System and method for frequency synchronization between transmitter and receiver in a communication system | |
| CN109660481B (en) | Quadrature amplitude modulation and demodulation method and constellation diagram optimization method thereof | |
| CN113381949A (en) | Noise estimation method, device, receiver and storage medium of physical broadcast channel | |
| CN113141324A (en) | Channel estimation method and device | |
| CN105871433A (en) | Virtual array antenna signal generating method based on high-speed mobile communication | |
| KR20170050175A (en) | Pilot signal generating apparatus and generating method thereof, transmission apparatus | |
| CN117499191A (en) | Phase estimation method and device based on statistical average | |
| CN109672641A (en) | A kind of LTE down channel estimation method suitable under complex environment |
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 |