干扰源的定位方法及装置与流程

本发明涉及通信领域，具体而言，涉及一种干扰源的定位方法及装置。

背景技术：

长期演进(Long Term Evolution，简称为LTE)系统最常遇到的干扰可以分为系统内干扰、系统外干扰几类。系统内干扰主要是同频干扰，包括如LTE时分双工(Time Division Duplexing简称为TDD)帧失步(GPS失锁)、TDD超远干扰、数据配置错误导致干扰、越区覆盖导致干扰等；系统外干扰主要是异系统非法使用LTE频段、以系统的杂散、阻塞或者互调干扰对本系统的影响。

外部干扰对于LTE网络性能影响脚大，对接入、掉话、速率都有较大影响。外部干扰主要包括码分多址(Code Division Multiple Access，简称为CDMA)干扰、小灵通干扰、广电干扰、干扰器干扰等。相关技术中排查干扰主要用扫频仪在楼顶使用八木天线判断干扰源方向，然后沿方向附近，根据经验找到附近另一个站点，再找出一个干扰源方向，最终根据两个方向交叉，找出干扰源的疑似位置，这种方法效率比较低。

针对相关技术中采用扫频仪找出干扰源的位置效率比较低的问题，目前还没有有效的解决方案。

技术实现要素：

本发明提供了一种干扰源的定位方法及装置，以至少解决相关技术中采用扫频仪找出干扰源的位置效率比较低的问题。

根据本发明的一个方面，提供了一种干扰源的定位方法，其特征在于，包括：获取以第一小区为中心预定范围的多个第二小区和所述第一小区组成的小区干扰组，其中，所述第一小区中有一个或多个资源块RB的接收干扰功率RIP值比所述第一小区内RB的RIP均值的高预定阈值，所述第一小区和多个所述第二小区分别有一段或多段连续的RB；在所述干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将所述路损均方差最小的点作为输出的干扰源位置，其中，所述经纬度信息包括：所述干扰小区组中的小区最大经纬度信息和最小经纬度信息。

进一步地，在获取以第一小区为中心预定范围的多个第二小区和所述第一小区组成的小区干扰组之前，所述方法还包括：获取网络内每一个小区中RB的RIP均值；将每一个小区中每个RB的RIP值与该小区的RIP均值进行取差运算；将差值大于所述预定 阈值的小区标记为干扰小区；记录所有所述干扰小区中RB值大于所述预定阈值的RB中具有连续标识ID的连续RB并为所述连续RB进行编号，以及记录所有所述连续RB的均值。

进一步地，获取以第一小区为中心预定范围的多个第二小区和所述第一小区组成的小区干扰组包括：获取以所述第一小区为中心预定范围内所有连续RB的编号；比较所述第一小区中连续RB编号和预定范围内其他小区的连续RB编号；将比较结果相同的所有连续RB以RIP均值的高低依次排序；将排序中排名靠前的多个小区标记为多个所述第二小区；获取以多个所述第二小区和所述第一小区组成的所述小区干扰组。

进一步地，获取每个栅格中心点的路损均方差包括：获取所述小区干扰组中各小区的参数信息，其中，所述参数信息包括：小区的RIP均值、小区天线高度、小区方位角；根据所述参数信息和预设规则计算得到每个栅格中心点的路损均方差。

进一步地，将所述路损均方差最小的点作为输出的干扰源位置之后，所述方法还包括：呈现所述第一小区中和多个所述第二小区中连续RB的RIP均值信息，以及所述干扰源位置信息。

根据本发明的另一个方面，提供了一种干扰源的定位装置，包括：第一获取模块，用于获取以第一小区为中心预定范围的多个第二小区和所述第一小区组成的小区干扰组，其中，所述第一小区中有一个或多个资源块RB的接收干扰功率RIP值比所述第一小区内RB的RIP均值的高预定阈值，所述第一小区和多个所述第二小区分别有一段或多段连续的RB；定位模块，用于在所述干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将所述路损均方差最小的点作为输出的干扰源位置，其中，所述经纬度信息包括：所述干扰小区组中的小区最大经纬度信息和最小经纬度信息。

进一步地，在获取以第一小区为中心预定范围的多个第二小区和所述第一小区组成的小区干扰组之前，所述装置还包括：第二获取模块，用于获取网络内每一个小区中RB的RIP均值；取差模块，用于将每一个小区中每个RB的RIP值与该小区的RIP均值进行取差运算；标记模块，用于将差值大于所述预定阈值的小区标记为干扰小区；记录模块，用于记录所有所述干扰小区中RB值大于所述预定阈值的RB中具有连续标识ID的连续RB并为所述连续RB进行编号，以及记录所有所述连续RB的均值。

进一步地，所述第一获取模块包括：第一获取单元，用于获取以所述第一小区为中心预定范围内所有连续RB的编号；比较单元，用于比较所述第一小区中连续RB编号和预定范围内其他小区的连续RB编号；排序单元，用于将比较结果相同的所有连续RB以RIP均值的高低依次排序；标记单元，用于将排序中排名靠前的多个小区标记为多个所述第二小区；第二获取单元，用于获取以多个所述第二小区和所述第一小区组成的所述小区干扰组。

进一步地，所述定位模块包括：第三获取单元，用于获取所述小区干扰组中各小区的参数信息，其中，所述参数信息包括：小区的RIP均值、小区天线高度、小区方位角；计算单元，用于根据所述参数信息和预设规则计算得到每个栅格中心点的路损均方差。

进一步地，将所述路损均方差最小的点作为输出的干扰源位置之后，所述装置还包括：呈现模块，用于呈现所述第一小区中和多个所述第二小区中连续RB的RIP均值信息，以及所述干扰源位置信息。

在本发明中，采用获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组，其中，该第一小区中有一个或多个资源块RB的接收干扰功率RIP值比第一小区内RB的RIP均值的高预定阈值，第一小区和多个第二小区分别有一段或多段连续的RB，在干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将路损均方差最小的点作为输出的干扰源位置，其中，该经纬度信息包括：干扰小区组中的小区最大经纬度信息和最小经纬度信息的方式，解决了相关技术中采用扫频仪找出干扰源的位置效率比较低的问题，到达了快速定位干扰源的效果。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是根据本发明实施例的干扰源的定位方法的流程图；

图2是根据本发明实施例的干扰源的定位装置的结构框图；

图3是根据本发明实施例的干扰源的定位装置的可选结构框图一；

图4是根据本发明实施例的干扰源的定位装置的可选结构框图二；

图5是根据本发明实施例的干扰源的定位装置的可选结构框图三；

图6是根据本发明实施例的干扰源的定位装置的可选结构框图四；

图7是根据本发明可选实施例的干扰排查方法的流程图。

具体实施方式

下文中将参考附图并结合实施例来详细说明本发明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。

在本实施例中提供了一种干扰源的定位方法，图1是根据本发明实施例的干扰源的定位方法的流程图，如图1所示，该流程包括如下步骤：

步骤S102：获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组，其中，第一小区中有一个或多个资源块(Resource Block简称为RB)的接收干扰功率接收干扰功率(Received Interference Power，简称为RIP)值比第一小区内RB的RIP均值的高预定阈值，第一小区和多个第二小区分别有一段或多段连续的RB；

步骤S104：在干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将路损均方差最小的点作为输出的干扰源位置，其中，经纬度信息包括：干扰小区组中的小区最大经纬度信息和最小经纬度信息。

在本实施例的步骤S102和步骤S104中，采用获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组，其中，该第一小区中有一个或多个资源块RB的接收干扰功率RIP值比第一小区内RB的RIP均值的高预定阈值，第一小区和多个第二小区分别有一段或多段连续的RB，在干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将路损均方差最小的点作为输出的干扰源位置，其中，该经纬度信息包括：干扰小区组中的小区最大经纬度信息和最小经纬度信息，也就是说在本实施例中采用首先计算出小区的平均RIP值，并依据该RIP均值和相关方案得到干扰小区组，进而根据该干扰小区组和预设方案找到干扰源位置的方式，解决了相关技术中采用扫频仪找出干扰源的位置效率比较低的问题，到达了快速定位干扰源的效果。

需要说明的是，上述步骤S102和步骤S104中涉及到的预设阈值可选为5，当然该预设阈值也可以是小于或大于5的整数，即可以根据需要对该预设阈值进行相应的设置。此外，对于本实施例中涉及到的多个第二小区，该多个第二小区在本实施例的具体应用场景中可以是6个第二小区，当然这里仅仅是举例说明，具体的取值可以根据实际情况来确定。

在本实施例的一个可选实施方式中，在获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组之前，本实施例的方法还可以包括：

步骤S11：获取网络内每一个小区中RB的RIP均值；

步骤S12：将每一个小区中每个RB的RIP值与该小区的RIP均值进行取差运算；

步骤S13：将差值大于预定阈值的小区标记为干扰小区；

步骤S14：记录所有干扰小区中RB值大于预定阈值的RB中具有连续标识ID的连续RB并为连续RB进行编号，以及记录所有连续RB的均值。

在上述步骤S11至步骤S14中，通过获取小区RIP均值的方式，将小区中每个RB的RIP值与该小区的RIP均值进行取差运算，并将差值大于预定阈值的小区标记为干扰 小区，该干扰小区也就是本实施例中涉及到的第一小区，同时还记录所有干扰小区(第一小区)中RB值大于预定阈值的RB中具有连续标识ID的连续RB并为连续RB进行编号，以及记录连续RB的均值。

基于上述步骤S11和步骤S14本实施例涉及到的步骤S102获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组的方式，在本实施例的可选实施方式中，可以通过如下方式来实现，该方式的步骤包括：

步骤S21：获取以第一小区为中心预定范围内所有连续RB的编号；

步骤S22：比较第一小区中连续RB编号和预定范围内其他小区的连续RB编号；

步骤S23：将比较结果相同的所有连续RB以RIP均值的高低依次排序；

步骤S24：将排序中排名靠前的多个小区标记为多个第二小区；

步骤S25：获取以多个第二小区和第一小区组成的小区干扰组。

通过该步骤S21至步骤S25，得到了本实施例中涉及到的小区干扰组；需要说明的是，本实施例中涉及到的第一小区是干扰小区，而干扰小区是有多个的，即第一小区也是有多个的，通过上述方式确定的小区干扰组也可以是有多个的。

而对于本实施步骤S104中涉及到的获取每个栅格中心点的路损均方差的方式，在本实施例可以通过如下方式来实现包括：

步骤S31：获取小区干扰组中各小区的参数信息，其中，参数信息包括：小区的RIP均值、小区天线高度、小区方位角；

步骤S32：根据参数信息和预设规则计算得到每个栅格中心点的路损均方差。

对于步骤S31和步骤S32，在具体的应用场景可以是根据工程参数获取每个干扰组小区的经纬度信息，根据小区组内最大经纬度信息和最小经纬度信息，形成一个矩形，在矩形内划分栅格，再结合小区的RIP均值、小区天线高度、小区方位角计算出每个栅格中心点的路损均方差，取路损均方差最小的那个点，作为干扰位置。

此外，在本实施例的另一个可选实施方式中，在将路损均方差最小的点作为输出的干扰源位置之后，本实施例的方法还可以包括：呈现第一小区中和多个第二小区中连续RB的RIP均值信息，以及干扰源位置信息。

而对于本实施例中的呈现方式，在本实施例的一个应用场景中可以是：输出每个干扰源位置、主小区以及6个小区的连续RB的RIP均值信息记录。同时每条记录可以显示相关小区的RB1～100的RIP均值分布，每个小区一条线，并在地图上展示，展示的方式可以是：在装置地图上将全网符合RB的RBAvgDiff大于5的所有小区(整网中所有符合条件的小区)，显示以颜色渲染扇区，并以颜色区分该段连续RBID的RIP均值。 同时显示以疑似干扰点位置为圆心，半径为100米的实心圆，颜色用红色标识。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到根据上述实施例的方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，或者网络设备等)执行本发明各个实施例的方法。

在本实施例中还提供了一种干扰源的定位装置，该装置用于实现上述实施例及可选实施方式，已经进行过说明的不再赘述。如以下所使用的，术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图2是根据本发明实施例的干扰源的定位装置的结构框图，如图2所示，该装置包括：第一获取模块22，用于获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组，其中，第一小区中有一个或多个资源块RB的接收干扰功率RIP值比第一小区内RB的RIP均值的高预定阈值，第一小区和多个第二小区分别有一段或多段连续的RB；定位模块24，与第一获取模块22耦合连接，用于在干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将路损均方差最小的点作为输出的干扰源位置，其中，经纬度信息包括：干扰小区组中的小区最大经纬度信息和最小经纬度信息。

图3是根据本发明实施例的干扰源的定位装置的可选结构框图一，如图3所示，在获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组之前，装置还包括：第二获取模块32，与取差模块34耦合连接，用于获取网络内每一个小区中RB的RIP均值；取差模块34，与标记模块36耦合连接，用于将每一个小区中每个RB的RIP值与该小区的RIP均值进行取差运算；标记模块36，与记录模块38耦合连接，用于将差值大于预定阈值的小区标记为干扰小区；记录模块38，与第一获取模块22耦合连接，用于记录所有干扰小区中RB值大于预定阈值的RB中具有连续标识ID的连续RB并为连续RB进行编号，以及记录连续RB的均值。

图4是根据本发明实施例的干扰源的定位装置的可选结构框图二，如图4所示，该第一获取模块22包括：第一获取单元402，用于获取以第一小区为中心预定范围内所有连续RB的编号；比较单元404，与第一获取单元402耦合连接，用于比较第一小区中连续RB编号和预定范围内其他小区的连续RB编号；排序单元406，与比较单元404耦合连接，用于将比较结果相同的所有连续RB以RIP均值的高低依次排序；标记单元408，用于将排序中排名靠前的多个小区标记为多个第二小区；第二获取单元410，用于获取以多个第二小区和第一小区组成的小区干扰组。

图5是根据本发明实施例的干扰源的定位装置的可选结构框图三，如图5所示，定位模块24包括：第三获取单元52，用于获取小区干扰组中各小区的参数信息，其中，参数信息包括：小区的RIP均值、小区天线高度、小区方位角；计算单元54，与第三获取单元52耦合连接，用于根据参数信息和预设规则计算得到每个栅格中心点的路损均方差。

图6是根据本发明实施例的干扰源的定位装置的可选结构框图四，如图6所示，将路损均方差最小的点作为输出的干扰源位置之后，装置还包括：呈现模块62，与定位模块24耦合连接，用于呈现第一小区中和多个第二小区中连续RB的RIP均值信息，以及干扰源位置信息。

下面结合本发明的可选实施例对本发明进行举例说明；

本可选实施例提供了一种LTE网络基于MR数据进行外部干扰排查的方法，该方法是一种基于测试报告(Measurement Report，简称为MR)数据，结合工参、无线数据以及天线损耗数据的网络干扰源确定以及定位方法，其中，该方法从MR中的接收干扰功率(Received Interference Power，简称为RIP)数据中计算出小区的平均RIP值，并依据具体判定算法得到干扰小区组，根据干扰小区组找到疑似干扰源位置在地图上进行定位分析，以辅助网优工程师快速查找和定位干扰源。

下面对本可选实施例中涉及到的分析流程包括数据导入、周期设置、数据处理、查找干扰组、定位干扰源、输出分析结果、地图展示等几个阶段进行相应的说明。

数据导入：将待分析的MR数据、工参数据、无线数据和天线损耗数据导入；

周期设置：由于干扰具有时间特征，某个干扰源有可能仅仅是在某个时间段开启，其它时段关闭，因此干扰排查需要进行时间设置，在装置上设置本次干扰排查的起始时间和结束时间以及排查时间粒度；

数据处理：对设定周期内的每个小区的RIP数据进行预处理分析，得到每个小区的连续资源块RB的RIP均值以及编号，需要注意的是一个小区有可能有好几段连续RB，每段连续RB都要记录一个RIP均值和RB编号；

查找干扰组：在所有小区中按照连续RB的RIP均值最高的小区从大到小排列。如果一个小区有多段连续RB，那么该小区排多次。按照顺序对小区进行干扰分析，得到干扰小区组，每个小区组在本可选实施例中最多有6个小区；

定位干扰源：根据工程参数获取每个干扰组小区的经纬度信息，根据小区组内最大经纬度信息和最小经纬度信息，形成一个矩形，在矩形内划分栅格，再结合小区的RIP均值、小区天线高度、小区方位角计算出每个栅格中心点的路损均方差，取路损均方差最小的那个点，作为疑似干扰位置；

输出分析结果：输出每个干扰源位置、主小区以及6个小区的连续RB的RIP均值 信息记录。同时每条记录可以显示相关小区的RB1～100的RIP均值分布，每个小区一条线；

地图展示：通过地图将全网符合RB的RBAvgDiff大于5的所有小区(整网中所有符合条件的小区)，显示以颜色渲染扇区，并以颜色区分该段连续RBID的RIP均值。同时显示以疑似干扰点位置为圆心，半径为100米的实心圆，颜色用红色标识。

下面结合附图和本可选实施例的实施例对本可选实施例上述分析流程进行详细说明；

图7是根据本发明可选实施例的干扰排查方法的流程图，如图7所示，该方法的步骤包括：

步骤S702：MR数据、工参、无线数据、天线数据入库；

步骤S704：设置处理数据周期；

步骤S706：查询周期内数据；

步骤S708：找到N组干扰小区；

步骤S710：分析干扰源位置；

步骤S712：输出分析结果。

对于步骤S702至步骤S712涉及到的过程下面将进行详细说明；

首先，以小区为单位数据处理：(1)计算每个小区的平均RIP值，按照每个小区的100个RB的RIP做平均，得到AvgRIP；

(2)对每个小区的每个RB的RIP与该小区的RIP均值进行比较，得到100个RB与均值的差值(RB1Diff～RB100Diff)。

(3)从RB1Diff到RB100Diff中，分别计算RBDiff大于0的所有RB的均值，以及RBDiff小于0的所有RB的均值，然后用前面均值减去后面均值，记为RBAvgDiff；

(4)如果RBAvgDiff大于5，说明有部分RB的干扰比较高，将此小区记为有干扰的小区，其中，RBAvgDiff大于5是本可选实施例的举例说明，并不构成对本可选实施例中其他方式的限定；

(5)判断RB中哪些RB是连续RB：将RBDiff大于5的所有RB以RBID从小到大排序，然后用第二个RB的ID减第一个RB的ID，判断是否为1，是则第一个RBID为连续RB的起始ID，以此循环处理，直到第N个RB的ID减第N-1个RB的ID不等于1，则第N-1个RB的ID为连续RB的结束ID，然后得到连续RB的个数，当个数大于等于2时，这段RB为连续RB，否则从第N+1个RB重新计算起始RBID的位置；

(6)如果RBAvgDiff小于等于5，判断AvgRIP是否大于100dBm，若是则将此小区的全带宽(1～N RB)号记录为干扰RB位置；

(7)记录这些连续RB的RIP均值，并记录连续RB的编号，便于后续与周边小区进行比较。注意一个小区有可能有好几段连续RB，每段连续RB都要记录一个RIP均值和RB编号。

其次，查找干扰组：在所有小区中按照连续RB的RIP均值最高的小区从大到小排列。如果一个小区有多段连续RB，那么该小区排多次，以连续RB的起始中心频率+RB个数为索引。

连续RB的起始中心频点计算：读取网管无线参数表中的中心频率earfcnul，bandwidthul分别得到绝对频率和上行带宽，然后按照一个RB 180k进行转换，将每个RBID转换为绝对频率，具体算法如表1：

表1

(1)从第一个小区开始，将其设置为主小区，并以其为圆心，找到距离3公里以内(可设置)的所有非共站同频小区，并查询出这些小区的连续RB的编号；

(2)在找到的同频小区中依次与主小区的连续RB编号进行对比，只要有一段连续的RB编号相同，则将这些小区记录下来；

(3)对(2)中找到的所有小区，根据RIP均值从高到低，并取出与主小区非共站的最多6个(可设置)RIP均值最高的小区(不含主小区)。6个小区RIP均值都要小于主小区的RIP均值，若大于则不处理此干扰组，进行下一个循环；

(4)得到6个小区的连续RB的RIP均值，若不足N个则结束，进行下一个循环；

再次，干扰源定位：

(1)从工参中确定6个小区的最大经度，最大纬度，最小经度和最小纬度，形成一个矩形，然后在矩形内分割为50*50栅格的经纬度集合，每个栅格都取该栅格中心点位置的经纬度，形成Point1～Pointn的集合；

(2)取出6个小区的RIP值，得到RIP1～RIP6；

(3)取出经纬度集合中的第一个经纬度，记为Point1，计算出该点的位置和6个小区的距离，并转换为公里的单位，分别以D1～D6表示；

(4)在工参中取出6个小区的天线高度H1～H6，单位为米；

(5)在工参中分别计算该点与6个小区所在站点的方位角(即起始点为小区，该点为假设的天线法线方向，得到一个方位角，方位角正北为0，顺时针逐渐增加)，得到Azimuth1～Azimuth6。

(6)分别用Azimuth1～Azimuth6减去每个小区的方位角，如果结果小于0则加上360，得到Azimuth1’～Azimuth6’。

(7)将Azimuth1’～Azimuth6’分别代入到天线损耗数据文件中的HorizontalDegree，计算出HorizontalLoss1～HorizontalLoss6。

(8)根据以上已经得到的信息，分别计算每个小区的pathloss：

pathloss＝-RIP-44.9*log(D)-13.82*Log(H)-6.55*Log(H)*Log(D)-HorizontalLoss，这样每个点得到6个pathloss的值，分别为pathloss1～pathloss6。

计算(pathloss1，pathloss2，pathloss3，pathloss4，pathloss5，pathloss6)的均方差，记为第一个点的StandDev1，以此方法循环得到所有点的StandDev；

(9)取StandDev最小的那个点，做为疑似干扰位置；

(10)当干扰点位置到主小区的距离大于1000米(可设置)，忽略此条干扰记录；当干扰点位置到辅助定位小区的距离大于1000米，则忽略此辅助小区；

(11)得到所有小区的干扰位置数据，最后结果可能会有重复的情况，因为一个干扰源会干扰多个小区，因此最后结果，对于同一个3公里以内的相同的连续干扰RB，只显示一行记录。干扰结果展示如表2：

表2

同时每条记录可以显示相关小区的RB1～100的RIP值分布图，每个小区一条线。在地图上找出全网所有符合该条RB的RBAvgDiff大于5的所有小区(整网中所有符合条件的小区)，显示以颜色渲染扇区，并以颜色区分该段连续RBID的RIP均值，颜色和区段个数均可设置。同时显示以疑似干扰点位置为圆心，半径为100米的实心圆，颜色用红色标识。

需要说明的是，上述各个模块是可以通过软件或硬件来实现的，对于后者，可以通过以下方式实现，但不限于此：上述模块均位于同一处理器中；或者，上述模块分别位于多个处理器中。

本发明的实施例还提供了一种存储介质。可选地，在本实施例中，上述存储介质可以被设置为存储用于执行以下步骤的程序代码：

S1，获取以第一小区为中心预定范围的多个第二小区和第一小区组成的小区干扰组，其中，第一小区中有一个或多个资源块RB的接收干扰功率接收干扰功率RIP值比第一小区内RB的RIP均值的高预定阈值，第一小区和多个第二小区分别有一段或多段连续的RB；

S2，在干扰小区组中的小区的经纬度信息所形成的矩形内划分预定个数的栅格，并获取每个栅格中心点的路损均方差，将路损均方差最小的点作为输出的干扰源位置，其中，经纬度信息包括：干扰小区组中的小区最大经纬度信息和最小经纬度信息。

可选地，在本实施例中，上述存储介质可以包括但不限于：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。

可选地，本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例，本实施例在此不再赘述。

显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。

以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员 来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。