CN110113272B

CN110113272B - 强窄带干扰的抑制方法、装置、通信设备和存储介质

Info

Publication number: CN110113272B
Application number: CN201910384075.7A
Authority: CN
Inventors: 彭岳峰; 周雄
Original assignee: Guangzhou Haige Communication Group Inc Co
Current assignee: Guangzhou Haige Communication Group Inc Co
Priority date: 2019-05-09
Filing date: 2019-05-09
Publication date: 2021-08-03
Anticipated expiration: 2039-05-09
Also published as: CN110113272A

Abstract

本申请涉及一种强窄带干扰的抑制方法、装置、通信设备和存储介质。所述方法包括：基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵；将P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量；根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵；利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对接收信号进行干扰抑制。利用该方法可以对接收信号进行快速瞬时的干扰抑制，从而可以有效地提高终端对接收信号的解码能力。

Description

强窄带干扰的抑制方法、装置、通信设备和存储介质

技术领域

本申请涉及计算机技术领域，特别是涉及一种强窄带干扰的抑制方法、装置、通信设备和存储介质。

背景技术

随着通信技术的不断发展，人们对通信数据的传输速率要求也越来越高，从刚开始的2G升级到3G，再从3G升级到4G，如今，随着第三代合作伙伴计划3GPP对5G标准的确立，从4G升级到5G的通信技术，以及5G的商用化也将会得到快速发展。

一般在4G通信体制中，正交频分复用***是基站和终端之间进行无线传输时一种比较经典的物理层空口传输技术，在专用通信领域(例如军用通信领域)的正交频分复用***中，有时会存在强窄带干扰(例如敌方恶意干扰)，当接收方不能及时感知干扰，并对干扰进行瞬时抑制，可能会导致接收方在接收数据时，解码失败，更严重的甚至会导致接收方和发送方之间的通信中断。

因此，在接收信号存在干扰时，如何对接收信号进行快速干扰抑制，以提高接收方的解码能力就成为亟待解决的技术问题。

发明内容

基于此，有必要针对上述技术问题，提供一种强窄带干扰的抑制方法、装置、通信设备和存储介质。

一种强窄带干扰的抑制方法，该方法包括：

基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵；

将P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量；

根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵；

利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对接收信号进行干扰抑制。

在其中一个实施例中，上述根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，包括：

利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值；

根据上述各组信道估计能量对应的干扰信息比特位值，确定各组信道在各个时刻是否存在干扰。

在其中一个实施例中，当上述门限包括初始阈值、第一预设阈值、第二预设阈值时，该第一预设阈值大于0，第二预设阈值小于0，上述利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值，包括：

将上述每组估计信道矩阵的信道估计能量分别与上述初始阈值、第一预设阈值、上述第二预设阈值进行对比，得到各组信道的对比结果；

根据上述各组信道的对比结果得到各组信道估计能量对应的干扰信息比特位值。

在其中一个实施例中，当上述门限包括一个绝对预设阈值时，上述利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值，包括：

将上述每组估计信道矩阵的信道估计能量分别与上述绝对预设阈值进行对比，得到各组信道的对比结果；

在其中一个实施例中，上述根据上述各组信道估计能量对应的干扰信息比特位值，确定各组信道在各个时刻是否存在干扰，包括：

若当前时刻存在为0的干扰信息比特位值，则确定上述为0的干扰信息比特位值对应的A组信道在当前时刻不存在干扰，以及确定为1的干扰信息比特位值对应的B组信道在当前时刻存在干扰。

在其中一个实施例中，上述基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵，包括：

基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，或者基于上述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻上述A组信道对应的本地信道矩阵；

将B组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述B组信道对应的本地信道矩阵。

在其中一个实施例中，上述基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，确定当前时刻上述A组信道对应的本地信道矩阵，包括：

根据公式(1)

确定当前时刻上述A组信道对应的本地信道矩阵

其中，α为预设的遗忘因子，

为当前时刻上述A组信道对应的估计信道矩阵，

为所述A组信道最近一个时刻不存在干扰时的本地信道矩阵，

的初始值为当前时刻上述A组信道对应的估计信道矩阵，t为当前时刻，t-1为最近一个时刻，k为第k个子载波，S为参考信号。

在其中一个实施例中，上述基于所述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻所述A组信道对应的本地信道矩阵，包括：

将A组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述A组信道对应的本地信道矩阵。

在其中一个实施例中，上述利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，包括：

根据公式(2)

计算各个时刻的强窄带干扰信号

其中，所述

为当前时刻所述B组信道对应的估计信道矩阵，所述

为所述B组信道当前时刻的本地信道矩阵，所述S为参考信号，所述t为当前时刻，所述m为所述B组信道中第m个子载波。

一种强窄带干扰的抑制装置，该装置包括：

估计模块，用于基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵；

计算模块，用于将上述P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量；

确定模块，用于根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵；

抑制模块，用于利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对所述接收信号进行干扰抑制。

一种通信设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：

一种可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

上述强窄带干扰的抑制方法、装置、通信设备和存储介质，首先基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵，接着将该P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量，然后根据每组信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵，最后利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并采用该强窄带干扰信号对接收信号进行干扰抑制。在本实施例中，由于可以计算得到各信道的估计信道矩阵、各组信道的信道估计能量、各组信道的本地信道矩阵，以及根据信道估计能量可以快速判断出是否存在干扰信息，并在存在干扰时，利用估计信道矩阵和本地信道矩阵可以快速计算出干扰信息，即可以快速感知到干扰信息，从而使得在一段时间内，可以利用计算出的干扰信息，对接收信号进行快速瞬时的干扰抑制，进而可以有效地提高对接收信号的解码能力。

附图说明

图1为一个实施例中强窄带干扰的抑制方法的应用环境图；

图2为一个实施例中强窄带干扰的抑制方法的流程示意图；

图3为另一个实施例中强窄带干扰的抑制方法的流程示意图；

图4为另一个实施例中强窄带干扰的抑制方法的流程示意图；

图5为另一个实施例中强窄带干扰的抑制方法的流程示意图；

图6为另一个实施例中强窄带干扰的抑制方法的流程示意图；

图7为另一个实施例中强窄带干扰的抑制方法的流程示意图；

图8为一个实施例中强窄带干扰的抑制装置的结构框图；

图9为另一个实施例中强窄带干扰的抑制装置的结构框图；

图10为一个实施例中通信设备的内部结构图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

本申请提供的强窄带干扰的瞬时抑制方法，可以应用于如图1所示的通信***中，该通信***可以是正交频分复用(orthogonal frequency division multiplexing，OFDM)***。其中，该***包括基站102和终端104，基站102和终端104之间可以通信，终端104可以是各种个人计算机、笔记本电脑、智能手机、平板电脑、计算机设备和便携式可穿戴设备等，本申请对终端104的具体形式并不做限定；可选的，基站102和终端104上均可以包括收发机，该收发机包括发射机和接收机。

需要说的是，本申请可以对上行信道进行估计，也可以对下行信道进行估计，上行信道指的是终端的发射机发射数据、同时基站的接收机接收数据时所在的信道，下行信道指的是基站的发射机发射数据、同时终端的接收机接收数据时所在的信道，本申请实施例是以发送下行参考信号，进行下行信道估计为例进行说明的。

在一个实施例中，如图2所示，提供了一种强窄带干扰的瞬时抑制方法，本实施例涉及的是终端如何基于参考信号对接收信号进行干扰抑制的具体过程，以该方法应用于图1中的终端为例进行说明，该方法包括以下步骤：

S202，基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵。

本实施例中，在OFDM***中，频域上一共包含M个子载波，基站可以将P个参考信号可以均匀的***到M个子载波中，之后基站可以通过这M个子载波将包含P个参考信号的数据发送给终端，可选的，基站发送的数据可以是时域数据；其中，参考信号也可以称为导频信号，可选的，P≤M，且M可以被P整除；终端在接收到带有参考信号的时域数据之后，可以将该时域数据进行快速傅里叶变换，之后就可以得到***导频信号的频域信号，该接收信号可以表示为公式(3)：

其中，

表示在t时刻、第k个子载波上对应参考信号的信道矩阵，k的范围为1到M，S表示参考信号，n₀表示高斯白噪声，

表示在t时刻、第k个子载波上可能存在的强窄带干扰(就是说，该第k个子载波在t时刻，可能存在干扰，也可能不存在干扰)。在上面的公式中，接收信号的值是已知的，S参考信号也是已知的，n₀的值变化太快，其能量较弱，随机性也较大，一般在干扰抑制时不予考虑，那么上面的公式就剩下参考信号的信道矩阵和强窄带干扰信号未知，那么就要对这两个信号数据进行估计计算。

具体地，在利用参考信号进行信道估计时，可以采用公式(4)

来进行计算，以得到参考信号的估计信道矩阵；在该公式中，

表示在t时刻、第k个子载波上对应参考信号的估计信道矩阵，S^H表示参考信号S的共轭转置，也是已知数据，根据该公式就可以得到P个参考信号的估计信道矩阵，即P个估计信道矩阵。

另外，将公式(3)带入到公式(4)中，则上述公式(4)还可以展开成如下公式(5)的形式：

公式(5)

可见，估计出的估计信道矩阵中包含可能存在的强窄带干扰信号

其中，

也可以称为强窄带干扰信号。

S204，将P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量。

其中，终端在对P个估计信道矩阵进行分组时，可选的，可以将P个估计信道矩阵均分成L组，也可以是不均分的分成L组；为了便于分组统计，本实施例优选地是均分成L组，即每个分组里的估计信道矩阵个数相同，可选的，L≤P，P可以被L整除，P和L均为整数。

具体的，在计算每组估计信道矩阵的能量时，可以采用公式(6)来进行计算，公式(6)为

其中，i∈[1，L]，P/L表示的是每组信道中包含的估计信道矩阵的个数，||||表示的是范数，通过对每组信道中各个估计信道矩阵进行范数平方求和，就可得到每组估计信道矩阵的信道估计能量。

S206，根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵。

其中，在瑞利信道或莱斯信道中，估计信道矩阵在频域内，估计信道矩阵的信道估计能量具有一定的连续性和相关性，利用这种连续性和相关性，就可以估计出接收信号中是否存在强窄带干扰。

需要说明的是，在本申请中，OFDM***在初始状态时是不存在干扰的，在某一时刻之后，加入强窄带干扰；另外，终端确定的本地信道矩阵一般都是不存在干扰时的信道矩阵。

具体的，终端在计算出每组估计信道矩阵的信道估计能量之后，可以对每组信道估计能量进行干扰识别，从而得到干扰确定结果，在得到干扰确定结果后，可选的，终端可以根据该干扰确定结果，和预设的遗忘因子、最近一个时刻不存在干扰时的本地信道矩阵、以及至少一组估计信道矩阵，来计算得到每组信道对应的本地信道矩阵；可选的，终端可以利用预设的干扰识别方法来对每组信道估计能量进行干扰识别。在该步骤中，最近一个时刻不存在干扰时的本地信道矩阵，其的初始值是该组信道对应的当前时刻的估计信道矩阵，因此，该步骤中是需要至少一组估计信道矩阵。

S208，利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对接收信号进行干扰抑制。

其中，各个时刻包括某一时刻；也就是说，终端在计算各个时刻的强窄带干扰时，也可以分别计算出某个时刻的强窄带干扰。

具体的，终端在计算得到每组信道对应的本地信道矩阵之后，可以将估计信道矩阵和本地信道矩阵做差值，可以根据上述公式(5)等再对差值进行变换，就可以计算得到各个时刻强窄带干扰信号，在得到强窄带干扰信号后，终端可以将上面得到的本地信道矩阵作为接收信号中的信道矩阵，并将接收信号与本步骤中计算得到的强窄带干扰信号做差值，就可以对接收信号进行去噪，也就是对接收信号进行干扰抑制，从而就可以使终端在解调接收信号时，解调得到的信号更准确。需要说明的是，本申请在对接收信号进行干扰抑制时，不会重新调度时频资源，会利用当前的时频资源进行干扰抑制。

需要说明的是，本实施例中终端在进行干扰抑制时，是在原有的时频资源上进行干扰抑制，而不会重新调度时频资源，因此，不会增加空口的开销。

上述强窄带干扰的抑制方法中，首先基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵，接着将该P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量，然后根据每组信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵，最后利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用该强窄带干扰信号对接收信号进行瞬时干扰抑制。在本实施例中，由于可以计算得到各信道的估计信道矩阵、各组信道的信道估计能量、各组信道的本地信道矩阵，以及根据信道估计能量可以快速判断出是否存在干扰信息，并在存在干扰时，利用估计信道矩阵和本地信道矩阵可以快速计算出干扰信息，即可以快速感知到干扰信息，从而使得在一段时间内，可以利用计算出的干扰信息，对接收信号进行快速瞬时的干扰抑制，进而可以有效地提高对接收信号的解码能力。

在另一个实施例中，如图3所示，提供了另一种强窄带干扰的瞬时抑制方法，本实施例涉及的是如何利用每组信道估计能量确定各组信道在各个时刻是否存在干扰的具体过程，在上述实施例的基础上，上述S206中根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，可以包括以下步骤：

S302，利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值。

其中，预设的干扰识别方法可以包括相对比较法、绝对门限比较法等，其中，相对比较法包括两个门限，分别为初始阈值、第一预设阈值、第二预设阈值，第一预设阈值大于0，第二预设阈值小于0；绝对门限比较法包括一个门限，为一个绝对预设阈值；另外，干扰信息比特位值可以用B表示，其表征的是每组信道是否存在干扰的情况。

具体的，终端在得到每组信道估计能量之后，可以通过预设的干扰识别方法对每组信道估计能量进行门限比较，就可以得到各组信道估计能量对应的干扰信息比特位值。

可选的，在一种可能的实施方式中，在利用相对比较法来进行门限比较时，如图4所示，可以包括下列两个步骤：

S402，将上述每组估计信道矩阵的信道估计能量分别与上述初始阈值、上述第一预设阈值、上述第二预设阈值进行对比，得到各组信道的对比结果。

S404，根据上述各组信道的对比结果得到各组信道估计能量对应的干扰信息比特位值。

具体的，在得到每组信道估计能量后，可以采用公式(7)的相对比较法来对各个信道估计能量进行门限比较，公式(7)可以表示为：

当i＝1时，

i>1时：

当B_i-1＝0时，

当B_i-1＝1时，

在公式(7)中，

为初始阈值，

为第一预设阈值，

为第二预设阈值，B为干扰信息比特位值，i∈[1,L]，G_i表示第i组信道估计能量，G_i-1表示第i-1组信道估计能量，初始时刻判定是否存在干扰，如公式(7-1)；在公式(7-2)中，每组信道估计能量都需要同时与

和

进行对比，最终根据每组信道估计能量与两个预设阈值的对比结果，确定是否存在干扰。对于某一时刻、每一组B_i的判定，均与上一组B_i-1的值相关，并以此为前提条件。例如，B_i-1为0时，当满足

时，B_i为0；否则为1；B_i-1为1时，当满足

时，B_i为0；否则为1。

可选的，在另一种可能的实施方式中，在利用绝对门限比较法来进行门限比较时，如图5所示，可以包括下列两个步骤：

S502，将上述每组估计信道矩阵的信道估计能量分别与上述绝对预设阈值进行对比，得到各组信道的对比结果。

S504，根据上述各组信道的对比结果得到各组信道估计能量对应的干扰信息比特位值。

具体的，在得到每组信道估计能量后，可以采用公式(8)的绝对门限比较法来对各个信道估计能量进行门限比较，公式(8)可以表示为：

在公式(8)中，G_thre为一个绝对预设阈值，每组信道估计能量分别于该绝对预设阈值进行对比，当某一组信道估计能量小于绝对预设阈值时，则就可以确定该组没有干扰，也可以根据该一组的对比结果得到干扰信息比特位值；当另外一组信道估计能量大于等于绝对预设阈值时，则也可以确定该另一组存在干扰，同时也可以根据该另一组的对比结果确定干扰信息比特位值。

上述干扰识别的集合可以用

表示。

S304，根据上述各组信道估计能量对应的干扰信息比特位值，确定各组信道在各个时刻是否存在干扰。

在该步骤中，可选的，若当前时刻存在为0的干扰信息比特位值，则确定上述为0的干扰信息比特位值对应的A组信道在当前时刻不存在干扰，以及确定为1的干扰信息比特位值对应的B组信道在当前时刻存在干扰。

其中，不存在干扰的A组可以是一组，也可以是多组，只要不存在干扰均称为A组，存在干扰的B组可以是一组，也可以是多组，只要存在干扰均称为B组，本实施例对A组、B组的个数不做限定。

具体的，根据上述S404和S504可以得到干扰信息比特位值，当干扰信息比特位值为0时，则确定不存干扰，当干扰信息比特位值为1时，则确定存在干扰，下面就可以计算存在干扰时的强窄带干扰信号了。

本实施例提供的强窄带干扰的抑制方法，首先利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值，然后根据上述各组信道估计能量对应的干扰信息比特位值，确定各组信道在各个时刻是否存在干扰。在本实施例中，由于利用干扰信息比特位值确定各组信道是否存在干扰的方法比较简单，也比较快速，因此，利用该方法可以快速确定接收信号中是否存在干扰，从而在接收信号中存在干扰时，可以快速地对接收信号进行干扰抑制。

在另一个实施例中，如图6所示，提供了另一种强窄带干扰的瞬时抑制方法，本实施例涉及的是终端如何根据上述干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵的具体过程，在上述实施例的基础上，上述S206中基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵，可以包括以下步骤：

S602，基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，或者基于上述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻上述A组信道对应的本地信道矩阵。

本实施例中，在上述S304的描述基础上可知，A组信道是不存在干扰的，即干扰信息比特位值为0的都是A组信道，B组信道是存在干扰的，也即是干扰信息比特位值为1的都是B组，则该S602的步骤可以包括两种情况，第一种是基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，确定当前时刻上述A组信道对应的本地信道矩阵；第二种是基于上述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻上述A组信道对应的本地信道矩阵。

可选的，上述第一种情况是在B组存在干扰时，不替换不存在干扰的A组的本地信道矩阵的方案，A组可以采用遗忘记忆法保存当前的信道矩阵，该方案可以包括：根据公式(1)

确定当前时刻上述A组信道对应的本地信道矩阵

其中，α为预设的遗忘因子，

为当前时刻上述A组信道对应的估计信道矩阵，

为所述A组信道最近一个时刻不存在干扰时的本地信道矩阵，

的初始值为当前时刻上述A组信道对应的估计信道矩阵，t为当前时刻，t-1为最近一个时刻，k为第k个子载波，S为参考信号。在该步骤中，可选的，可以根据上述公式(1)确定当前时刻上述A组信道对应的本地信道矩阵

还可以通过对公式(1)进行可以简化的变形，确定当前时刻上述A组信道对应的本地信道矩阵

本实施例对此并不做限定。

可选的，上述第二种情况是在B组存在干扰时，同时替换不存在干扰的A组的本地信道矩阵的方案，该方案可以包括：将A组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述A组信道对应的本地信道矩阵。也就是说，A组信道当前时刻对应的本地信道矩阵可以被替换为最近一个时刻A组信道对应的本地信道矩阵。

S604，将B组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述B组信道对应的本地信道矩阵。

具体的，在确定B组信道存在干扰时，可以将B组信道当前时刻对应的本地信道矩阵替换为最近一个时刻不存在干扰时，B组信道对应的本地信道矩阵。

本实施例提供的强窄带干扰的抑制方法，首先基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，或者基于上述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻上述A组信道对应的本地信道矩阵，然后将B组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述B组信道对应的本地信道矩阵，得到各个时刻每组信道对应的本地信道矩阵。在本实施例中，由于本地信道矩阵都是不存在干扰时的信道矩阵，因此，利用该方法计算得到的每组信道的本地信道矩阵，在后续对接收信号进行干扰抑制时，得到的结果会更准确，也即会使得对接收信号的解码能力提高。

在另一个实施例中，提供了另一种强窄带干扰的抑制方法，本实施例涉及的终端如何计算各个时刻的强窄带干扰信号，并利用计算的强窄带干扰信号对接收信号进行干扰抑制的具体过程，在上述实施例的基础上，上述S208中的利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号可以包括：

根据公式(2)

计算各个时刻的强窄带干扰信号

其中，所述

为当前时刻所述B组信道对应的估计信道矩阵，所述

具体的，在上述确定B组信道存在干扰时，假设只有一个B组存在干扰，上述可以根据干扰信息比特位值得到如下存在干扰的公式，如公式(9)所示，

其中，在公式(9)中，l为存在干扰的一组，可以根据公式(9)计算出存在干扰的l组中每个信道的序号m，m可以表示为：

根据上述序号m，可以利用公式(2)估计出强窄带干扰信号，可选的，可以根据上述公式(2)计算各个时刻的强窄带干扰信号，还可以通过对公式(2)进行可以简化的变形，计算各个时刻的强窄带干扰信号，本实施例对此并不做限定；在本步骤中，

也称为估计强窄带干扰信号。

在公式(2)中，

可以通过上述公式(4)计算得到，在这里为已知值，

可以用最近一个时刻B组的本地信道矩阵

代替，而

也可以根据上述S602和S604确定，在这里为已知值，因此根据公式(2)就可以计算得到估计强窄带干扰信号；在这里，根据公式(4)，上述公式(2)还可以展开如下公式(10)所示：

在公式(10)中，

为强窄带干扰估计误差，假设终端移动速度为中低时速时，则在一段时间间隔内(毫秒级或者数十毫秒级)，快衰落信道的多径效应、多普勒效应可以认为不变，即存在较强的相关性，那么

忽略高斯白噪声n0，那么上面公式(10)中的估计强窄带干扰信号就可以认为是强窄带干扰信号。

需要说明的是，上述S602和S604是在存在强窄带干扰时，对信道进行重新拟合的过程，包括替换全部子载波的信道矩阵和仅替换存在强窄带干扰信号的信道的信道矩阵，上述A组的第一种情况和B组的结合属于仅替换存在强窄带干扰信号的信道的信道矩阵，上述A组的第二种情况和B组的结合属于替换全部子载波的信道矩阵。

上述替换全部子载波的信道矩阵表示为如下：

上述仅替换存在强窄带干扰信号的信道的信道矩阵可以表示为如下：

其中

和

均表示拟合后的信道矩阵。

在上述S602和S604得到P个信道的本地信道矩阵之后，即得到P个信道不存在干扰时的信道矩阵之后，如果上述S202中的全部信道数量M大于P，那么可以用拟合后的信道矩阵进行插值，得到全频段，即M个信道对应的信道矩阵，以便后续用于数据解调。

在上述计算得到全部信道的信道矩阵和强窄带干扰信号之后，接收信号可以用公式(11)表示，如下所示：

其中，

d表示实际数据传输过程中的信号，这时可以利用上面公式(10)得到的强窄带干扰信号对接收信号进行去噪，如下公式(12)所示：

其中，

表示在滤除强窄带干扰时产生的误差。在接收信号去噪后，对接收信号进行频域均衡，就可以得到均衡后的数据，即就可以解调出基站发送的数据。

本实施例提供的强窄带干扰的抑制方法，可以根据每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算得到各个时刻的强窄带干扰信号，并可以利用该强窄带干扰信号对接收信号进行干扰抑制，也就是去噪。在本实施例中，由于可以根据公式快速计算出强窄带干扰信号，并利用该强窄带干扰信号对接收信号进行去噪，因此，可以有效地提高对接收信号的瞬间解码能力，提高通信的质量和通信的稳定性。

为了便于本领域技术人员的理解，以下对本发明提供的强窄带干扰的抑制方法进行详细介绍，如图7所示，该方法可以包括：

S701，基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵。

S702，将P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量。

S703，基于每组信道估计能量，计算每组信道的干扰信息比特位值。

S704，判断各组信道的干扰信息比特位值是否为0，若是，则执行S705，若否，则执行S706。

S705，确定为0的干扰信息比特位值对应的A组信道不存在干扰，并基于当前时刻A组信道对应的估计信道矩阵和预设的遗忘因子，或者基于A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻A组信道对应的本地信道矩阵。

S706，确定为1的干扰信息比特位值对应的B组信道存在干扰，并将B组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻B组信道对应的本地信道矩阵。

S707，利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对所述接收信号进行干扰抑制。

应该理解的是，虽然图2-7的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2-7中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

在一个实施例中，如图8所示，提供了一种强窄带干扰的抑制装置，包括：估计模块10、计算模块11、确定模块12和抑制模块13，其中：

估计模块10，用于基于P个参考信号对接收信号进行信道估计，得到P个估计信道矩阵；

计算模块11，用于将上述P个估计信道矩阵分成L组，并计算每组估计信道矩阵的信道估计能量；

确定模块12，用于根据每组估计信道矩阵的信道估计能量确定各组信道在各个时刻是否存在干扰，并基于干扰确定结果和至少一组估计信道矩阵，确定各个时刻每组信道对应的本地信道矩阵；

抑制模块13，用于利用每组信道对应的本地信道矩阵以及每组信道对应的估计信道矩阵，计算各个时刻的强窄带干扰信号，并利用强窄带干扰信号对所述接收信号进行干扰抑制。

本实施例提供的强窄带干扰的抑制装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

在另一个实施例中，如图9所示，提供了另一种强窄带干扰的抑制装置，在上述实施例的基础上，上述确定模块12包括比较单元121、第一确定单元122，其中，

比较单元121，用于利用预设的干扰识别方法对上述每组估计信道矩阵的信道估计能量进行门限比较，得到各组信道估计能量对应的干扰信息比特位值；

第一确定单元122，用于根据上述各组信道估计能量对应的干扰信息比特位值，确定各组信道在各个时刻是否存在干扰。

可选的，当上述门限包括初始阈值、第一预设阈值、第二预设阈值时，该第一预设阈值大于0，第二预设阈值小于0，上述比较单元121还用于将上述每组估计信道矩阵的信道估计能量分别与上述初始阈值、上述第一预设阈值、上述第二预设阈值进行对比，得到各组信道的对比结果；根据上述各组信道的对比结果得到各组信道估计能量对应的干扰信息比特位值。

可选的，当上述门限包括一个绝对预设阈值时，上述比较单元121还用于将上述每组估计信道矩阵的信道估计能量分别与上述绝对预设阈值进行对比，得到各组信道的对比结果；根据上述各组信道的对比结果得到各组信道估计能量对应的干扰信息比特位值。

可选的，上述第一确定单元122还用于若当前时刻存在为0的干扰信息比特位值，则确定上述为0的干扰信息比特位值对应的A组信道在当前时刻不存在干扰，以及确定为1的干扰信息比特位值对应的B组信道在当前时刻存在干扰。

在另一个实施例中，请继续参见图9所示，上述确定模块12可以包括：第二确定单元123、第三确定单元124，其中，

第二确定单元123，用于基于当前时刻上述A组信道对应的估计信道矩阵和预设的遗忘因子，或者基于上述A组信道最近一个时刻不存在干扰时的本地信道矩阵，确定当前时刻上述A组信道对应的本地信道矩阵；

第三确定单元124，用于将B组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述B组信道对应的本地信道矩阵。

可选的，上述第二确定单元123还用于根据公式(1)

确定当前时刻上述A组信道对应的本地信道矩阵

其中，α为预设的遗忘因子，

为当前时刻上述A组信道对应的估计信道矩阵，

为所述A组信道最近一个时刻不存在干扰时的本地信道矩阵，

可选的，上述第二确定单元123还用于将A组信道最近一个时刻不存在干扰时的本地信道矩阵确定为当前时刻上述A组信道对应的本地信道矩阵。

在另一个实施例中，上述抑制模块13还用于根据公式(2)

计算各个时刻的强窄带干扰信号

其中，所述

为当前时刻所述B组信道对应的估计信道矩阵，所述

在一个实施例中，提供了一种通信设备，该通信设备可以是终端或者基站，其内部结构图可以如图10所示。该通信设备包括连接的收发机、处理器、存储器、网络接口、显示屏和输入装置。其中，该通信设备的处理器用于提供计算和控制能力。该通信设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***和计算机程序。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该通信设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种强窄带干扰的抑制方法。

本领域技术人员可以理解，图10中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的通信设备的限定，具体的通信设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种通信设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：

在一个实施例中，提供了一种可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。