CN112235002B - 谐波消除方法及装置、存储介质、终端 - Google Patents
谐波消除方法及装置、存储介质、终端 Download PDFInfo
- Publication number
- CN112235002B CN112235002B CN202011461607.1A CN202011461607A CN112235002B CN 112235002 B CN112235002 B CN 112235002B CN 202011461607 A CN202011461607 A CN 202011461607A CN 112235002 B CN112235002 B CN 112235002B
- Authority
- CN
- China
- Prior art keywords
- data block
- input data
- harmonic
- current
- frequency point
- 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 41
- 230000008030 elimination Effects 0.000 title claims abstract description 30
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 30
- 238000004590 computer program Methods 0.000 claims description 11
- 238000004891 communication Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010295 mobile communication Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
- H04B2001/1072—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal by tuning the receiver frequency
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
一种谐波消除方法及装置、存储介质、终端，谐波消除方法包括：确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。本发明技术方案能够在不影响系统稳定性的情况下实现对接收机谐波干扰的消除。
Description
技术领域
本发明涉及通信技术领域，尤其涉及一种谐波消除方法及装置、存储介质、终端。
背景技术
在第五代移动通信中，接收机工作在某些特定的频点上，手机主板上的某些系统模块的时钟信号的谐波分量会耦合进入接收机，导致接收机在弱场情况下，灵敏度会受到比较严重的影响，可能导致掉话或者通话信号质量下降。
目前消除谐波干扰的主要方案是从干扰源头处理。现有主流方案有三种：方案一：找到时钟干扰源头，降低时钟幅度，从而降低时钟谐波分量；方案二：找到时钟干扰源头，给时钟加上模拟滤波器，比如低通滤波器，从而衰减高频分量，降低时钟谐波分量；方案三:找到时钟干扰源头，改变时钟的频率，从而改变时钟谐波的频率，时钟谐波分量避开接收机的接收频点。
但是，方案一降低时钟幅度会造成系统稳定性下降；方案二需要增加硬件成本，而且往往受限于面积，难以增加，如果是空间耦合造成的干扰，即使增加滤波器，也无法有效抑制干扰；方案三时钟频率往往是固定的，不能随意修改，另外即使修改，如果接收机的接收频带比较宽，其时钟谐波干扰也很难避开。
发明内容
本发明解决的技术问题是如何在不影响系统稳定性的情况下实现对接收机谐波干扰的消除。
为解决上述技术问题，本发明实施例提供一种谐波消除方法，谐波消除方法包括：确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
可选的，所述根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分之后还包括：利用所述当前输入数据块的谐波分量确定下一搜索周期的初始频点，并按照预设搜索间隔确定预设数量个下一待搜索频点，所述下一待搜索频点包括所述下一搜索周期的初始频点，所述下一待搜索频点用于计算所述下一输入数据块的谐波分量。
可选的，所述利用所述当前输入数据块的谐波分量确定下一搜索周期的初始频点包括：如果所述当前搜索周期的初始频点与预设值之和小于所述当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点增加所述预设搜索间隔，以作为所述下一搜索周期的初始频点；如果所述当前搜索周期的初始频点与预设值之和大于所述当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点减小所述预设搜索间隔，以作为所述下一搜索周期的初始频点。
可选的，采用以下公式计算所述预设值：(SearchNum/2)×F_search_step，其中，SearchNum为所述预设数量，F_search_step为搜索预设搜索间隔。
可选的，所述根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量包括：计算当前输入数据块与各个当前待搜索频点的搜索相关值；确定大于预设门限的所述搜索相关值对应的当前待搜索频点作为所述当前搜索周期的输入数据块的谐波分量。
可选的，所述根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除包括：计算所述谐波相关值和所述当前输入数据块的谐波分量的乘积；计算所述下一输入数据块与所述乘积的差值，以作为谐波消除后的数据块。
可选的，所述确定当前搜索周期的初始频点包括：如果所述当前搜索周期为首个搜索周期，则确定所述当前搜索周期的初始频点为预设的初始频点；如果所述当前搜索周期非首个搜索周期，则利用上一输入数据块的谐波分量确定所述当前搜索周期的初始频点。
可选的，所述预设的初始频点是通过对多个输入数据块进行谐波分量测试得到的。
可选的，采用以下公式计算所述谐波相关值：
；其中，
表示第m-1个输入数据块与第m个输入数据块的谐波相关值，N表示码片的
数量或者符号数量，Input(m-1)表示第m-1个输入数据块，F_init表示所述初始频点，
表示所述预设搜索间隔。
可选的，所述谐波分量的搜索周期小于等于输入数据块的数据周期。
为解决上述技术问题，本发明实施例还公开了一种谐波消除装置，谐波消除装置包括：初始频点确定模块，用于确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；谐波分量估计模块，用于根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；谐波相关值计算模块，用于至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；谐波消除模块，用于根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
本发明实施例还公开了一种存储介质，其上存储有计算机程序，所述计算机程序被处理器运行时执行所述谐波消除方法的步骤。
本发明实施例还公开了一种终端，包括存储器和处理器，所述存储器上存储有可在所述处理器上运行的计算机程序，所述处理器运行所述计算机程序时执行所述谐波消除方法的步骤。
与现有技术相比，本发明实施例的技术方案具有以下有益效果：
本发明技术方案通过计算待搜索频点与输入数据块的相关性可以准确地估计出输入数据块中的谐波分量；同时，在确定当前数据周期的待搜索频点时，初始频点是根据上一输入数据块的谐波分量实时更新的，能够动态跟踪并估算出时钟谐波的频点和能量；进而基于谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除，大幅提高通话信号质量，减少掉话概率。
附图说明
图1是本发明实施例一种谐波消除方法的流程图；
图2是本发明实施例一种具体应用场景的示意图；
图3是本发明实施例一种谐波消除装置的结构示意图。
具体实施方式
如背景技术中所述，方案一降低时钟幅度会造成系统稳定性下降；方案二需要增加硬件成本，而且往往受限于面积，难以增加，如果是空间耦合造成的干扰，即使增加滤波器，也无法有效抑制干扰；方案三时钟频率往往是固定的，不能随意修改，另外即使修改，如果接收机的接收频带比较宽，其时钟谐波干扰也很难避开。
本发明技术方案通过计算待搜索频点与输入数据块的相关性可以准确地估计出输入数据块中的谐波分量；同时，在确定当前数据周期的待搜索频点时，初始频点是根据上一输入数据块的谐波分量实时更新的，能够动态跟踪并估算出时钟谐波的频点和能量；进而基于谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除，大幅提高通话信号质量，减少掉话概率。
相对于现有技术中对产生干扰的时钟信号进行修改，而修改时钟信号会导致使用该时钟信号的系统模块性能不稳定，本发明技术方案则不会对产生干扰的时钟信号产生影响，也就不会影响系统模块性能稳定性。
本发明技术方案可以应用于传统移动通信系统，例如宽带码分多址（WidebandCode Division Multiple Access, WCDMA）、长期演进(Long Term Evolution，LTE)系统，也可以用于新无线（New Radio，NR）系统，如第五代移动通信技术(5th generation mobilenetworks or 5th generation wireless systems, 5G)，或者更新的通信系统，本发明实施例对此不作限制。
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施例做详细的说明。
图1是本发明实施例一种谐波消除方法的流程图。
图1所示谐波消除方法可以用于终端设备侧，具体可以是配置有接收机的终端设备。图1所示方法能够对接收机的接收信号进行谐波消除，以提升接收信号的信噪比，提升信号质量。
具体地，图1所示谐波消除方法可以包括以下步骤：
步骤S101：确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；
步骤S102：根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；
步骤S103：至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；
步骤S104：根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
需要指出的是，本实施例中各个步骤的序号并不代表对各个步骤的执行顺序的限定。
本实施例中，输入数据块是指接收机接收到的数据块，输入数据块是周期性数据，其传输周期称为数据周期。
在步骤S101和步骤S102的具体实施中，对当前输入数据块进行谐波估计，也即估计出当前输入数据块的谐波分量的大小。为了准确地估计出谐波分量的大小，首先确定多个待搜索频点，具体是以初始频点为中心/起始频点/结束频点/任意中间频点，按照预设搜索间隔来选取；再计算当前输入数据块与各个当前待搜索频点的相关性。输入数据块与待搜索频点的相关性越高，表示输入数据块与待搜索频点余额接近，从而可以依据相关性的大小确定与输入数据块最接近的待搜索频点，以作为当前输入数据块的谐波分量。
需要说明的是，可以采用现有技术中任意可实施的算法计算相关性，本发明实施例对此不作限制。
在一个非限制性的例子中，初始频点为2140MHz，以初始频点为中心，预设搜索间隔为5MHz，预设数量为13，则各个待搜索频点分别为2110、2115、2120、2125、2130、2135、2140、2145、2150、2155、2160、2165、2170（单位为MHz）。
可以理解的是，预设搜索间隔、预设数量的具体数值可以根据实际应用场景的需求来设置，如预设搜索间隔越小、预设数量越多，估计的准确性越高。计算的复杂度也相应增加，本发明实施例对此不作限制。
在步骤S103和步骤S104的具体实施中，实现对下一输入数据块的谐波消除。由于谐波是相对稳态的，不会突变，同时实时对下一周期进行谐波估计和消除在实现上比较困难，因此本发明实施例是利用当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
具体而言，Output(m)=Input(m)-Correlation（m-1）×ref_value(m-1)，其中，Output(m)是第m个数据块经过谐波消除后的输出，Input(m)是第m个输入数据块，Correlation（m-1）是第m-1个数据块的谐波分量与第m个数据块的谐波分量的相关值，ref_value(m-1)是第m-1个数据块的谐波分量。
具体而言，谐波相关值表示所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的相似程度。当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值可以基于当前输入数据块的谐波分量进行估计，具体可以根据当前输入数据块的谐波分量、当前输入数据块的初始频点以及预设搜索间隔进行估计。
在一个非限制性的实施例中，谐波相关值可以采用以下公式计算：
；其中，N表示
通信系统中码片的数量（如WCDMA中码片的数量）或者符号数量（如LTE/ NR中符号数量），
Input(m-1)表示第m-1个输入数据块，F_init表示所述初始频点，
表示所述
预设搜索间隔。
具体地，N是码片数量或符号数量取决于终端所使用的通信系统，如果终端使用的通信系统是WCDMA，则N为码片(chip)数量；如果终端使用的是4G/5G通信系统，由于其隶属于LTE/NR，其调制方式是正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)，则N为符号（symbol）数量。
本发明技术方案避免对时钟信号的源头做额外处理，只需要通过软件控制数字基带完成谐波的估算、动态跟踪和消除，减少了系统风险，并有效提升灵敏度性能。
在一个非限制性的实施例中，可以对下一输入数据块进行谐波估计，也即估计出下一输入数据块的谐波分量的大小。具体的计算方式与前述步骤类似，也即首先确定多个待搜索频点，具体是以初始频点为中心/起始频点/结束频点/任意中间频点，按照预设搜索间隔来选取；再计算下一输入数据块与各个下一待搜索频点的相关性。输入数据块与待搜索频点的相关性越高，表示输入数据块与待搜索频点余额接近，从而可以依据相关性的大小确定与输入数据块最接近的待搜索频点，以作为下一输入数据块的谐波分量。
其中，在确定下一待搜索频点时，下一搜索周期的初始频点与当前搜索周期的初始频点可以相同，也可以不同。下一搜索周期的初始频点具体可以根据当前输入数据块的谐波分量与当前搜索周期的初始频点的相对位置来确定。
具体地，当前搜索周期的初始频点与预设值之和小于当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点增加所述预设搜索间隔；反之，当前搜索周期的初始频点与预设值之差大于当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点减小所述预设搜索间隔，以作为所述下一搜索周期的初始频点。
在一个非限制性的例子中，预设值可以是(SearchNum/2)×F_search_step，其中，SearchNum为预设数量，F_search_step为预设搜索间隔。
通过对每个数据周期进行谐波估计所采用的初始频点进行更新，可以实现对谐波频点的动态追踪，实现谐波估计的准确性。
在一个非限制性的实施例中，图1所示步骤S102可以包括以下步骤：计算当前输入数据块与各个当前待搜索频点的搜索相关值；确定大于预设门限的所述搜索相关值对应的当前待搜索频点作为所述当前搜索周期的输入数据块的谐波分量。
如前所述，输入数据块与待搜索频点的相关性越高，表示输入数据块与待搜索频点余额接近，故而本实施例选取搜索相关值的最大值对应的当前待搜索频点作为所述当前搜索周期的输入数据块的谐波分量。
在一个非限制性的实施例中，图1所示步骤S104可以包括以下步骤：计算所述谐波相关值和所述当前输入数据块的谐波分量的乘积；计算所述下一输入数据块与所述乘积的差值，以作为谐波消除后的数据块。
谐波相关值表示当前输入数据块与下一输入数据块的相似度，在此基础上，可以利用当前输入数据块的谐波分量对下一输入数据块进行谐波消除，也即将相应的谐波分量从输入数据块中剔除即可。
在一个非限制性的实施例中，图1所示步骤S101可以包括以下步骤：如果所述当前搜索周期为首个搜索周期，则确定所述当前搜索周期的初始频点为预设的初始频点；如果所述当前搜索周期非首个搜索周期，则利用所述上一输入数据块的谐波分量确定所述当前搜索周期的初始频点。
每一搜索周期所使用的初始频点是动态更新的。具体是根据上一搜索周期计算得到的谐波分量确定的，如果不存在上一搜索周期，则直接确定为预设的初始频点。
在根据上一搜索周期计算得到的谐波分量确定当前搜索周期的初始频点时，上一搜索周期的初始频点与预设值之和小于谐波分量，则增大上一搜索周期的初始频点，作为当前搜索周期的初始频点；反之，上一搜索周期的初始频点与预设值之和大于谐波分量，则减小上一搜索周期的初始频点，作为当前搜索周期的初始频点。
进一步而言，所述预设的初始频点是通过对多个输入数据块进行谐波分量测试得到的。
本实施例中，时钟信号对输入数据块产生的干扰谐波分量通常是稳定的，因此可以通过预先对多个输入数据块进行测试，获得一个确定谐波分量的频点，作为预设的初始频点。
在一个非限制性的实施例中，所述谐波分量的搜索周期小于等于输入数据块的数据周期。
在一个具体的应用场景中，结合具体例子对初始频点的更新过程进行说明。以宽带码分多址（Wideband Code Division Multiple Access, WCDMA）为例，宽带的下行频率范围是2110M-2170MHz。
1、设置预设的初始频点F_init=2140MHz，搜索10M范围内的频点，即搜索频率范围为2139M-2149M。
2、以100KHz为预设搜索间隔搜索谐波(spur)，即F_search_step=100KHz。
3、移动间隔F_step=100KHz。
4、共搜索SearchNum=（2145M-2135M）/0.1M=100个频点。
5、如果谐波的位置在2145.1M，计算F_init+(SearchNum/2)×F_search_step=2140+(100/2)*0.1=2140+5=2145M；F_init+(SearchNum/2) ×F_search_step<Ftone，即2145M<2145.1M。
6、更新初始频点F_init = F_init +F_step =2140+0.1M=2140.1M，继续迭代。
7、如果谐波的位置在2139.9M，计算F_init+(SearchNum/2)×F_search_step=2140.1+(100/2)*0.1=2140+5=2145.1M，F_init+(SearchNum/2) ×F_search_step>Ftone，即2145.1M>2139.9M。
8、更新初始频点F_init = F_init -F_step =2140.1-0.1M=2140M，继续迭代。
在另一个具体的应用场景中，请参照图2，对第m-1个数据块进行谐波估计，得到第m-1个数据块的谐波参考值。利用第m-1个数据块的谐波参考值对第m个数据块进行谐波消除。
以此类推，对第m个数据块进行谐波估计，得到第m个数据块的谐波参考值。利用第m个数据块的谐波参考值对第m+1个数据块进行谐波消除。
具体的谐波估计和谐波消除的实施方式可参照前述实施例，此处不再赘述。
图3是本发明实施例一种谐波消除装置的结构示意图。谐波消除装置30可以包括：
初始频点确定模块301，用于确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；
谐波分量估计模块302，用于根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；
谐波相关值计算模块303，用于至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；
谐波消除模块304，用于根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
本发明实施例通过计算待搜索频点与输入数据块的相关性可以准确地估计出输入数据块中的谐波分量；同时，在确定当前数据周期的待搜索频点时，初始频点是根据上一输入数据块的谐波分量实时更新的，能够动态跟踪并估算出时钟谐波的频点和能量；进而基于谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除，大幅提高通话信号质量，减少掉话概率。
关于所述谐波消除装置30的工作原理、工作方式的更多内容，可以参照图1至图2中的相关描述，这里不再赘述。
本发明实施例还公开了一种存储介质，所述存储介质为计算机可读存储介质，其上存储有计算机程序，所述计算机程序运行时可以执行图1中所示的谐波消除方法的步骤。所述存储介质可以包括ROM、RAM、磁盘或光盘等。所述存储介质还可以包括非挥发性存储器（non-volatile）或者非瞬态（non-transitory）存储器等。
本发明实施例还公开了一种终端，所述终端可以包括存储器和处理器，所述存储器上存储有可在所述处理器上运行的计算机程序。所述处理器运行所述计算机程序时可以执行图1中所示的谐波消除方法的步骤。所述终端包括但不限于手机、计算机、平板电脑等终端设备。
应理解，上述的处理器可以是通用处理器、数字信号处理器(digital signalprocessor，DSP)、专用集成电路(application specific integrated circuit，ASIC)、现成可编程门阵列(field programmable gate array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件，还可以是系统芯片(system on chip，SoC)，还可以是中央处理器(central processor unit，CPU)，还可以是网络处理器(networkprocessor，NP)，还可以是数字信号处理电路(digital signal processor，DSP)，还可以是微控制器(micro controller unit，MCU)，还可以是可编程控制器(programmable logicdevice，PLD)或其他集成芯片。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
还应理解，本发明实施例中提及的存储器可以是易失性存储器或非易失性存储器，或可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(read-only memory，ROM)、可编程只读存储器(programmable ROM，PROM)、可擦除可编程只读存储器(erasable PROM，EPROM)、电可擦除可编程只读存储器(electrically EPROM，EEPROM)或闪存。易失性存储器可以是随机存取存储器(random access memory，RAM)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(static RAM，SRAM)、动态随机存取存储器(dynamic RAM，DRAM)、同步动态随机存取存储器(synchronous DRAM，SDRAM)、双倍数据速率同步动态随机存取存储器(double datarate SDRAM，DDR SDRAM)、增强型同步动态随机存取存储器(enhanced SDRAM，ESDRAM)、同步连接动态随机存取存储器(synchlink DRAM，SLDRAM)和直接内存总线随机存取存储器(direct rambus RAM，DR RAM)。应注意，本文描述的系统和方法的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
需要说明的是，当处理器为通用处理器、DSP、ASIC、FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件时，存储器(存储模块)集成在处理器中。应注意，本文描述的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(read-only memory，ROM)、随机存取存储器(random access memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。
Claims (13)
1.一种谐波消除方法，其特征在于，包括：
确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；
根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；
至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；
根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
2.根据权利要求1所述的谐波消除方法，其特征在于，所述根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量之后还包括：
利用所述当前输入数据块的谐波分量确定下一搜索周期的初始频点，并按照预设搜索间隔确定预设数量个下一待搜索频点，所述下一待搜索频点包括所述下一搜索周期的初始频点，所述下一待搜索频点用于计算所述下一输入数据块的谐波分量。
3.根据权利要求2所述的谐波消除方法，其特征在于，所述利用所述当前输入数据块的谐波分量确定下一搜索周期的初始频点包括：
如果所述当前搜索周期的初始频点与预设值之和小于所述当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点增加所述预设搜索间隔，以作为所述下一搜索周期的初始频点；
如果所述当前搜索周期的初始频点与预设值之和大于所述当前输入数据块的谐波分量，则将所述当前搜索周期的初始频点减小所述预设搜索间隔，以作为所述下一搜索周期的初始频点。
4.根据权利要求3所述的谐波消除方法，其特征在于，采用以下公式计算所述预设值：(SearchNum/2)×F_search_step，其中，SearchNum为所述预设数量，F_search_step为搜索预设搜索间隔。
5.根据权利要求1所述的谐波消除方法，其特征在于，所述根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量包括：
计算当前输入数据块与各个当前待搜索频点的搜索相关值；
确定大于预设门限的搜索相关值对应的当前待搜索频点作为所述当前搜索周期的输入数据块的谐波分量。
6.根据权利要求1所述的谐波消除方法，其特征在于，所述根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除包括：
计算所述谐波相关值和所述当前输入数据块的谐波分量的乘积；
计算所述下一输入数据块与所述乘积的差值，以作为谐波消除后的数据块，所述当前输入数据块与所述下一输入数据块为接收机的接收信号。
7.根据权利要求1所述的谐波消除方法，其特征在于，所述确定当前搜索周期的初始频点包括：
如果所述当前搜索周期为首个搜索周期，则确定所述当前搜索周期的初始频点为预设的初始频点；
如果所述当前搜索周期非首个搜索周期，则利用上一输入数据块的谐波分量确定所述当前搜索周期的初始频点。
8.根据权利要求7所述的谐波消除方法，其特征在于，所述预设的初始频点是通过对多个输入数据块进行谐波分量测试得到的。
9.根据权利要求1所述的谐波消除方法，其特征在于，采用以下公式计算所述谐波相关值：
10.根据权利要求1所述的谐波消除方法，其特征在于，所述谐波分量的搜索周期小于等于输入数据块的数据周期。
11.一种谐波消除装置，其特征在于，包括：
初始频点确定模块，用于确定当前搜索周期的初始频点，并按照预设搜索间隔确定预设数量个当前待搜索频点，所述当前待搜索频点包括所述初始频点；
谐波分量估计模块，用于根据当前输入数据块与各个当前待搜索频点的相关性确定所述当前输入数据块的谐波分量；
谐波相关值计算模块，用于至少根据所述当前输入数据块的谐波分量估计所述当前输入数据块的谐波分量与下一输入数据块的谐波分量的谐波相关值；
谐波消除模块，用于根据所述谐波相关值和所述当前输入数据块的谐波分量对所述下一输入数据块进行谐波消除。
12.一种存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器运行时执行权利要求1至10中任一项所述谐波消除方法的步骤。
13.一种通信终端，包括存储器和处理器，所述存储器上存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器运行所述计算机程序时执行权利要求1至10中任一项所述谐波消除方法的步骤。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011461607.1A CN112235002B (zh) | 2020-12-14 | 2020-12-14 | 谐波消除方法及装置、存储介质、终端 |
| PCT/CN2021/130950 WO2022127476A1 (zh) | 2020-12-14 | 2021-11-16 | 谐波消除方法及装置、存储介质、终端 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011461607.1A CN112235002B (zh) | 2020-12-14 | 2020-12-14 | 谐波消除方法及装置、存储介质、终端 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112235002A CN112235002A (zh) | 2021-01-15 |
| CN112235002B true CN112235002B (zh) | 2021-03-09 |
Family
ID=74124631
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011461607.1A Active CN112235002B (zh) | 2020-12-14 | 2020-12-14 | 谐波消除方法及装置、存储介质、终端 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112235002B (zh) |
| WO (1) | WO2022127476A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112235002B (zh) * | 2020-12-14 | 2021-03-09 | 展讯通信（上海）有限公司 | 谐波消除方法及装置、存储介质、终端 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104868956A (zh) * | 2015-04-14 | 2015-08-26 | 陈景竑 | 基于声波信道的数据通信方法 |
| GB2517847B (en) * | 2012-05-21 | 2015-08-26 | Aceaxis Ltd | Detection of intermodulation products |
| CN111313918A (zh) * | 2018-12-29 | 2020-06-19 | 上海司南卫星导航技术股份有限公司 | 一种接收装置、终端装置和计算机可读存储介质 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101233711B (zh) * | 2005-10-24 | 2012-11-21 | 松下电器产业株式会社 | 干扰信号特征量存储方法和装置、干扰信号特征量获取方法和装置、干扰信号抑制方法和装置 |
| KR100770839B1 (ko) * | 2006-04-04 | 2007-10-26 | 삼성전자주식회사 | 음성 신호의 하모닉 정보 및 스펙트럼 포락선 정보,유성음화 비율 추정 방법 및 장치 |
| JP4869455B2 (ja) * | 2009-08-28 | 2012-02-08 | 三菱電機株式会社 | ラジオ受信機 |
| GB2511865B (en) * | 2013-03-15 | 2015-07-01 | Aceaxis Ltd | Processing interference in a wireless network |
| CN107222226B (zh) * | 2016-03-21 | 2019-08-27 | 华为技术有限公司 | 消除谐波干扰的方法和装置 |
| CN108345039B (zh) * | 2018-01-12 | 2019-07-23 | 吉林大学 | 一种消除地面核磁共振数据中邻频谐波干扰的方法 |
| DE102019135901A1 (de) * | 2019-02-25 | 2020-08-27 | Samsung Electronics Co., Ltd. | Drahtloskommunikationsvorrichtung und kanalschätzverfahren davon |
| CN110931035B (zh) * | 2019-12-09 | 2023-10-10 | 广州酷狗计算机科技有限公司 | 音频处理方法、装置、设备及存储介质 |
| CN112235002B (zh) * | 2020-12-14 | 2021-03-09 | 展讯通信（上海）有限公司 | 谐波消除方法及装置、存储介质、终端 |
-
2020
- 2020-12-14 CN CN202011461607.1A patent/CN112235002B/zh active Active
-
2021
- 2021-11-16 WO PCT/CN2021/130950 patent/WO2022127476A1/zh active Application Filing
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2517847B (en) * | 2012-05-21 | 2015-08-26 | Aceaxis Ltd | Detection of intermodulation products |
| CN104868956A (zh) * | 2015-04-14 | 2015-08-26 | 陈景竑 | 基于声波信道的数据通信方法 |
| CN111313918A (zh) * | 2018-12-29 | 2020-06-19 | 上海司南卫星导航技术股份有限公司 | 一种接收装置、终端装置和计算机可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112235002A (zh) | 2021-01-15 |
| WO2022127476A1 (zh) | 2022-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111404853B (zh) | 一种载波频偏估计方法、装置及计算机存储介质 | |
| KR20130100557A (ko) | 무선통신 시스템에서 메모리 클럭 주파수를 제어하는 장치 및 방법 | |
| JP5765120B2 (ja) | 無線通信端末及び自動利得制御方法 | |
| CN113050030B (zh) | 基于到达角度测距的定位方法及装置 | |
| CN112235002B (zh) | 谐波消除方法及装置、存储介质、终端 | |
| CN109729032B (zh) | 一种频偏估计值的校正方法、装置及计算机可读存储介质 | |
| CN115865594A (zh) | 一种基于fdd系统的无源交调校正方法及系统 | |
| CN101959299B (zh) | 自动频率控制方法及接收端 | |
| KR20160121267A (ko) | 이동 통신 시스템에서 기지국의 채널 추정 방법 및 장치 | |
| WO2017080359A1 (zh) | 一种干扰消除方法、装置及基站 | |
| CN114755629B (zh) | 提高室内定位精度的方法及装置 | |
| CN108111445B (zh) | 一种数字前端均衡的方法和装置 | |
| CN113359161B (zh) | 高动态卫星通信信号捕获方法、装置、介质及计算设备 | |
| CN115396055A (zh) | 信道预测方法及装置 | |
| CN106896382B (zh) | 一种导航信号捕获方法及装置 | |
| CN110602740B (zh) | 确定同频干扰的方法、装置、存储介质及用户设备 | |
| JP2003510885A (ja) | 信号の時間追跡方法およびその装置 | |
| CN112769725A (zh) | 基于全相位频谱纠正的Costas序列时频联合同步方法 | |
| Safapourhajari et al. | On the low complexity implementation of the DFT-based BFSK demodulator for ultra-narrowband communications | |
| CN116455719B (zh) | 频偏估计方法、装置、通信设备和可读存储介质 | |
| US9148319B2 (en) | Dynamic task scheduling for multi-receive-path equalizer | |
| CN114726689B (zh) | 一种信号估计的方法及装置 | |
| CN115052331B (zh) | 整周跳变误差的处理方法、装置及存储介质 | |
| CN116347329B (zh) | 定位信号捕获方法及装置、计算设备及计算机存储介质 | |
| CN112433191B (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 |