上行干扰定位装置、方法及计算机可读存储介质与流程

1.本公开涉及上行干扰定位装置、上行干扰定位方法及计算机可读存储介质。更具体地，涉及基于mr(measurement report，测量报告)数据实现的上行干扰定位技术。


背景技术：

2.在无线通信网中，近年来都是2/3/4g基站共存的局面，lte(long term evolution，长期演进)基站容易受到上行干扰的影响，包括阻塞、互调和杂散干扰，此外还有其他无线电设备引起的外部同频干扰。
3.干扰问题严重将直接影响用户体验度，例如，在初始接入或切换过程中，由于上行干扰的影响，用户发出的上行信令可能丢失，从而导致接入、切换失败。另外，lte多为同频组网，处于小区边缘的用户容易受到相邻小区的干扰，可能由于服务小区给用户分配的资源位置与相邻小区用户的资源位置重叠，从而相互干扰而导致上行信道质量降低、上行业务速率降低；进而，如果上行控制信道受到强烈干扰，必然导致控制信道的解调能力下降，由于下行数据包的上行反馈的解调失败，将导致不必要的重传，从而造成链路重建掉话等网络问题。


技术实现要素：

4.目前，通常采用的排查lte系统中上行干扰源的方法，很大程度地依赖于网络优化人员的经验判断，所以对网络优化人员的专业水平以及工作经验要求较高。另外，对于一个受干扰小区来说，其干扰电平呈现的波形往往包含多种干扰源，而各种干扰源都具备其各自特有的规律以及特征，适用的排查方法也不相同，所以在由网络优化人员判断干扰源的类型以及定位干扰源位置时，准确性和工作效率上都受到较大的限制。
5.本公开的一个目的在于提供一种能够准确判定系统外干扰并且精确定位系统外干扰的干扰源位置的上行干扰定位装置、上行干扰定位方法及计算机可读存储介质。
6.在下文中给出了关于本公开的简要概述，以便提供关于本公开的一些方面的基本理解。但是，应当理解，这个概述并不是关于本公开的穷举性概述。它并不意图用来确定本公开的关键性部分或重要部分，也不是意图用来限定本公开的范围。其目的仅仅是以简化的形式给出关于本公开的某些概念，以此作为稍后给出的更详细描述的前序。
7.根据本公开的一个方面，提供一种上行干扰定位装置，其中，上述上行干扰定位装置包括处理电路，上述处理电路被配置为：获取来自用户终端装置的测量报告数据即mr数据；按照预定的统计周期，对与所述用户终端装置所在的lte小区对应的多个物理资源块带宽的上行接收干扰信号功率即prb粒度enb接收干扰功率进行统计，生成所述lte小区的数字化干扰波形；对生成的所述数字化干扰波形进行向量处理而生成干扰波形向量，基于预先存储的干扰波形向量模板，对从所述干扰波形向量减去特征已知的干扰的波形向量而剩余的干扰波形向量进行波形识别；在通过所述波形识别判定为所述lte小区存在系统外干扰的情况下，根据所述mr数据对所述系统外干扰的干扰源的位置进行可视化呈现。
8.根据本公开的另一方面，提供一种上行干扰定位方法，其中，上述上行干扰定位方法包括：获取来自用户终端装置的测量报告数据即mr数据；按照预定的统计周期，对与所述用户终端装置所在的lte小区对应的多个物理资源块带宽的上行接收干扰信号功率即prb粒度enb接收干扰功率进行统计，生成所述lte小区的数字化干扰波形；对生成的所述数字化干扰波形进行向量处理而生成干扰波形向量，基于预先存储的干扰波形向量模板，对从所述干扰波形向量减去特征已知的干扰的波形向量而得到的剩余干扰波形向量进行波形识别；在通过所述波形识别判定为所述lte小区存在系统外干扰的情况下，根据所述mr数据对所述系统外干扰的干扰源的位置进行可视化呈现。
9.根据本公开的另一方面，提供一种计算机可读存储介质，包括可执行指令，当上述可执行指令由计算机执行时，使上述计算机执行上述上行干扰定位方法。
10.根据本公开，能够准确判定干扰源类型并且精确定位干扰源位置。
附图说明
11.构成说明书的一部分的附图描述了本公开的实施例，并且连同说明书一起用于解释本公开的原理。参照附图，根据下面的详细描述，可以更清楚地理解本公开，其中：
12.图1是示意地示出本公开涉及的测量报告即mr数据的采集过程的一个例子的示意图；
13.图2是示意地示出存在于lte系统中的各种干扰类型及其对应的波形的一个例子的示意图；
14.图3是示出干扰向量波形的一个例子的图；
15.图4是示出波形向量匹配处理的一个例子的示意图；
16.图5是示意地示出本公开的实施例的上行干扰定位装置2000的构成的例子的框图；
17.图6是示意地示出本公开的实施例的上行干扰定位装置2000执行的处理的例示性的流程图；
18.图7是示出本公开的基于ripprb的数字化干扰波形的一个例子的示意图。
19.图8是示出本公开的实施例的上行干扰定位装置2000执行的对系统外干扰的干扰源的位置进行可视化呈现的处理的一个例子的流程图；
20.图9是示出receivedipower上行干扰信号强度区间颜色配置表的一个例子的示意图；
21.图10是示出本公开的栅格渲染的一个例子的示意图。
22.图11是示出一个具体实施例的系统外干扰源的定位时的栅格渲染的示意图。
23.图12示意地示出实现本公开的实施例的计算设备300的示例性的配置框图。
具体实施方式
24.在下文中，参照附图详细地描述本公开内容的优选实施例。需要注意的是，在本说明书和附图中，用相同的附图标记来表示具有基本相同的功能和结构的结构元件，并且省略对这些结构元件的重复说明。
25.首先，为便于理解，对本公开涉及的lte系统的频点频率的对应关系、干扰定位的
数据源、lte系统干扰的类型、基于快速模版匹的波形识别算法等进行说明。
26.1-1.《lte系统的频点频率的对应关系》
27.lte系统分为fdd(frequency-division duplex))和tdd(time-division duplex，时分双工)两种双工通信模式，fdd即上下行信道采用不同频段进行信号发送；tdd即上下行信道使用相同频率发送信号，但上下行信道的发送时间不同，也就是在同一频率不同时隙进行上下行信号发送。其中，例如，可以通过参照预先存储的“lte频点与频率对应表”，计算上下行频段的每个频点所对应的频率。
28.例如，可以参照如下公式(1)计算上行频点对应的频率：
29.ful＝ful_low 0.1(nul
–noffs-ul
)(公式1)
30.在公式1中，ful为上行频点的频率，ful_low为对应频段的最低上行频点的频率，nul为该载频的上行频点号，n
offs-ul
为对应频段的最低上行频点号。
31.在lte网络无线资源中，上行接收的干扰功率(received interfere power，rip)定义为一个物理资源块(physical radio block，prb)带宽上的干扰功率，包括热噪声。prb是lte系统中的通信载体，一个prb的频域宽度为180khz且时域上连续。在lte协议中，lte网络的上行无线带宽可以设置为1.4mhz到20mhz等多种不同类型，以20mhz为例，有100个prb，共18mhz，两边各留10mhz做保护带，id从prb0到prb99，中心频率对应于49号prb和50号prb的交界处dc子载波的中心位置，所以，50号prb起始频率就是小区中心频率加上7.5khz，49号rb起始频率就是小区中心频率减去187.5khz，其它prb的频率例如可以按照180khz一个prb进行计算。
32.1-2.《干扰定位的数据源》
33.另外，本公开在判定上行干扰类型以及对干扰源进行定位时，主要基于测量报告即mr数据。测量报告是lte系统的一项重要功能，可以用于系统中无线资源控制子层完成诸如小区选择/重选及切换等事件的触发，也可以用于系统操作维护，观察系统的运行状态。mr数据的每一条测量报告都包括测量到的发生业务过程中所占用的主服务小区和相邻小区的信息，包括receivedipower上行接收的干扰功率、enb天线到达角(aoa)、时间提前量(tadv)等关键测量信息。
34.以下，结合图1说明本公开涉及的测量报告(measurement report，mr)统计数据的采集过程。如图1的示意图所示，mr数据主要来自ue(user equipment，用户设备)和enodeb(evolved node b，即演进型node b，以下有时简称为“enb”，其表示lte中的基站)的物理层、rlc(radio link control，无线链路控制)层以及在无线资源管理过程中计算产生的测量报告，在业务信道上每480ms(信令信道上470ms)发送一次，可用于对网络中基站的干扰情况进行评估。原始测量数据可以以统计数据形式(例如，可以在enodeb或无线接入网网元管理系统(omc-r)上实现统计)报送到omc-r进行存储(例如，参照图1的(a))或者直接报送到omc-r以样本数据形式进行存储(例如，参照图1的(b))。
35.在本公开中，mr数据中至少包括以下表1所示的数据：
36.表1
37.名称含义主要作用a_gps部分用户所在地理位置，经纬度定位用户所在的位置receivedipower服务小区上行干扰功率表示服务小区的干扰情况
ripprbprb粒度enb接收干扰功率表示物理资源块的干扰情况aoa用户与服务小区的方位角定位用户所在的位置tadv时间提前量测量用户与服务小区距离enbid服务小区标记结合工参确定服务小区位置
38.其中，数据“a_gps”包含手机接入时的经度地理坐标“ac.longitude”和纬度“ac.latitude”，用于定位用户所在位置。
39.数据“receivedipower”表示服务小区的上行信号接收干扰功率，例如，其表示服务小区所有prb的上行接收干扰信号功率的平均值，在本公开中，数据“receivedipower”用于干扰源定位。
40.表2 receivedipower的取值范围
[0041][0042]
数据“ripprb”表示prb粒度enb接收干扰功率，是prb所在的180khz带宽内的平均干扰值。由于各prb的频率范围不一样，所以基站omc网管按照统计周期对lte小区的各prb的上行接收干扰信号功率的样本值进行统计，例如存储为表3所示的统计表，用于生成lte小区的数字化干扰波形。
[0043]
表3 ripprb统计表
[0044][0045]
数据“aoa”表示enb天线到达角，其定义了一个用户终端装置例如手机相对参考方向的估计角度。该测量数据可用于确定用户终端装置相对lte基站的方位，与时间提前量(time advance，tadv)、基站经纬度结合，可实现ue的位置定位。
[0046]
表4 aoa取值范围表
[0047][0048]
数据“tadv”为时间提前量(time advance)，其用于定义用户终端装置用于调整其主小区物理上行控制信道/物理上行共享信道/探测参考信号(pucch/pusch/srs)上行发送的时间。具体计算方法为：在随机接入过程，enodeb通过测量接收到导频信号来确定时间提前值，时间提前量取值范围为(0,1,2,...,1282)
×
16ts；在rrc连接状态下，enodeb基于测量对应ue的上行传输来确定每个ue的tadv调整值，这个调整值的范围为(0,1,2,...,63)
×
16ts。本次得到的最新的时间提前量即为上次记录的时间提前量与本次enodeb测量得到的调整值之和。
[0049]
利用mr数据来计算距离的方法：
[0050]
1ts＝1/(15000*2048)s
[0051]
对应的距离等于：
[0052]
(3*10^8*1/(15000*2048))/2＝4.89m
[0053]
1tadv＝16ts＝16*4.89m＝78.12m
[0054]
表5 tadv的取值范围表
[0055]
测量报告统计数据测量数据区间分布(单位ts)mr.tadv.00tadv《16mr.tadv.0116≤tadv《32
……
mr.tadv.11176≤tadv《192mr.tadv.12192≤tadv《224
……
mr.tadv.37992≤tadv《1024mr.tadv.381024≤tadv《1280
……
mr.tadv.411792≤tadv《2048mr.tadv.422048≤tadv《3072mr.tadv.433072≤tadv《4096mr.tadv.444096≤tadv
[0056]
数据“tadv”能够与数据aoa、lte基站经纬度相结合，实现ue的位置定位。
[0057]
数据“enbid”是用于识别lte小区的标识，与工参结合能够确定lte小区的基站位置(例如，经纬度)。
[0058]
另外，本公开中还涉及地里信息系统(geographic information system，gis)，gis是一种基于计算机仿真技术，能够对工件信息进行分析和处理。
[0059]
1-3.《lte系统干扰的类型》
[0060]
在lte系统中，干扰分为系统内干扰和系统外干扰。系统内干扰包括：无源器件(主要是室分站，包括合路器，功分器，耦合器，天馈干扰)和有源器件(例如rru(remote radio unit，射频拉远单元))干扰、因用户数增加导致底噪抬升形成的阻塞干扰、lte天馈系统因发射多个无线电频率产生的互调干扰、设备运行热噪声等。系统内容干扰主要通过优化系统布局、优化网络扩容方案以及提高电缆连接质量等方式来降低。另外，系统外干扰主要包括：杂散干扰、外部非法设立的无线电设备干扰、固定运行的无线电设备故障等。系统外干扰需要进行定位排查，以增强lte系统运行质量。本公开涉及的主要是对系统外干扰的判定以及定位技术。
[0061]
另外，各种干扰类型的波形都有其各自独特的特征，例如，可以根据各种干扰类型的波形特征预先创建波形向量模板，并存储到波形向量特征数据库。然后，通过基于计算机程序实现的快速模板匹配的波形识别算法，对输入的干扰波形进行类型识别。
[0062]
1-3-1.阻塞干扰
[0063]
在干扰源系统的发射频率与受害系统的接收频率距离较近、或者接收机滤波器的阻塞抑制能力较差时，可能产生图2(a-1)所示的阻塞干扰。如图2(a-2)所示，阻塞干扰的波形具有如下特征：
[0064]
①
干扰电平值呈现整体平直抬升的频谱特性；
[0065]
②
阻塞干扰与lte系统的流量吸收强度、prb利用率及网络负荷呈现正相关；
[0066]
③
空间特征分析：相同区域的lte各小区，具有较低程度的共同时频干扰特征。
[0067]
根据阻塞干扰波形的上述特征，可事先设计阻塞干扰的波形向量模板，建立阻塞干扰波的向量特征库x。
[0068]
1-3-2.谐波干扰
[0069]
由于发射机有源器件和无源器件的非线性，在其发射频率的整数倍频率上将产生较强的谐波产物。如图2(b-1)所示，当这些谐波产物正好落于受害系统接收机频段内，将导致受害接收机灵敏度损失，称之为谐波干扰。如图2(b-2)所示，谐波干扰波形具有如下特征：
[0070]
①
干扰电平值出现1-2个尖峰频谱特性，在prb的干扰特征中，受影响的prb组不连续。
[0071]
②
空间特征分析：相同区域的lte各小区，没有明显的同时频特征。
[0072]
根据谐波干扰波形的上述特征，可事先设计谐波干扰的波形向量模板，建立谐波干扰波的向量特征库x。
[0073]
1-3-3.互调干扰
[0074]
当两个或多个不同频率的发射信号通过非线性电路时，将在多个频率的线性组合频率上形成互调产物。如图2(c-1)所示，当这些互调产物与受害接收机的有用信号频率相同或相近时，将导致受害接收机灵敏度损失，称之为互调干扰。如图2(c-2)所示，互调干扰波形具有如下特征：
[0075]
①
干扰电平值出现多个尖峰的频谱特性；
[0076]
②
频域和时域比较固定；
[0077]
③
空间特征分析：相同区域的lte各小区，没有明显的同时频特征。根据互调干扰波形的上述特征，可事先设计互调干扰的波形向量模板，建立互调干扰波的向量特征库x。
[0078]
1-3-4.杂散干扰
[0079]
由于某些发射机中的功放、混频器和滤波器等稳定性问题，非线性器件在工作频带以外很宽的范围内产生辐射信号分量，图2(d-1)所示，包括谐波、寄生辐射、频率转换产物和互调产物等落入受害系统接收频段内，导致受害接收机的底噪抬升，造成灵敏度损失，称之为杂散干扰。此外，某些自然现象(如雷电)偶发性的发生瞬间干扰，也属于杂散干扰。如图2(d-2)所示，杂散干扰波形具有如下特征：
[0080]
①
时域上不连续，突发特征明显；
[0081]
②
频域上不固定，偶发性强；
[0082]
③
lte系统靠近干扰源发射频段的prb更容易受到干扰；
[0083]
④
空间特征分析：相同区域的lte各小区，具有较低程度的同时频干扰特征。根据杂散干扰波形的上述特征，可事先设计杂散干扰的波形向量模板，建立杂散干扰波的向量特征库x。
[0084]
1-3-5.外部干扰源干扰
[0085]
由于某些发射机非法设立或没有滤波设计，在受害lte系统工作频率带内产生的持续性的外部干扰源干扰。并且，外部干扰源干扰波形具有如下特征：
[0086]
①
时域上连续或具有一定的规律性，如24小时持续出现、周末定时出现、晚上定时出现等，主要与外部干扰源的加电工作时间相关，因而具有一定的时域特征；
[0087]
②
频域上稳定，对lte系统prb的干扰持续在固定的位置中；
[0088]
③
对lte各小区的干扰强度较为稳定；
[0089]
④
空间特征分析：当外部存在一个固定的干扰源时，该干扰源能对相同区域的lte各小区同时产生相同类别的干扰，具有较强程度的同时频干扰特征。根据外部干扰源干扰波形的上述特征，可事先设计外部干扰源干扰的波形向量模板，建立外部干扰源干扰波的向量特征库x。
[0090]
1-4.《基于快速模版匹的波形识别》
[0091]
本公开涉及的快速模板匹配的波形识别大致包括如下步骤：
[0092]
1-4-1.波形向量化
[0093]
在检测波形、字符检测波形等波形中，可以看到每种波形都有其特有的形状，只是在外部条件改变的情况下(如电压大小、采样频率等)，波形的幅值有所不同。
[0094]
对于每个数字化干扰波形w，可以用含有m个元素的向量x来表示：x＝[x1 x2
ꢀ…
xm-2 xm-1 xm]。
[0095]
其中，在本公开中，xm是数字化干扰波形在第m个prb采样时间t内的平均幅值，例如：xm可表示为：xm＝(采样时间t内，第m个prb上的平均rip幅值xm)。向量中的相邻两元素间的采样间隔δf相等。图3示出本公开的干扰波形向量的一个例子的示意图，图中的横坐标表示prb的编号，纵坐标表示每个prb对应的干扰值的平均值。另外，例如，20m带宽的系统中包括100个prb，在对数字化干扰波形进行向量化处理时，既可以将全部的100个prb的干扰平均值都抽样为为向量元素，即向量包括100个向量元素，也可以按照一定的间隔进行采样，例如间隔3个prb进行采样，即波形对应的向量最终包括25个向量元素。采样间隔的设定与模板向量的采样间隔一致即可。
[0096]
1-4-2.建立干扰波形特征向量模版库
[0097]
例如，在本公开中，根据上述各种干扰类型的特征，预先建立与各种干扰类型对应的干扰特征向量模板库xi，作为与输入信号进行匹配的模板。这里，需要注意的是，考虑到信号幅值呈正负分布，不对波形向量模板中采用的向量元素进行归一化处理。为了解决幅值一致性的问题，可以对输入波形y进行(公式2)所示的处理：
[0098][0099]
在公式2中，max(y)表示输入向量的最大值，max(xi)表示模板向量的最大值。
[0100]
1-4-3.向量波形匹配
[0101]
在向量波形匹配过程中，定义输入波形向量y与波形模板xi之间的距离di：
[0102][0103]
根据公式3，针对输入波形向量y的各元素yi与波形向量模板的各元素xi计算di值，对针对各模板结算出的di的值进行比较，取di的最小值，将输入波形y识别为di的最小值对应的模板xi的向量，输出识别结果。图4示出向量波形匹配的一个例子的示意图。图4的(a)示出的是杂散干扰对应的波形向量模板的模板1，图4的(b)示出的是外部干扰源干扰对应的波形向量模板的模板2，图4的(c)示出的是输入波形向量，图4的(d)示出波形匹配处理的结果。根据图4可知，由于输入波形向量与模板2之间的距离di最小，所以输入波形向量被识别为是对应于外部干扰源的向量。
[0104]
以下，结合附图详细说明本公开涉及的上行干扰定位处理。
[0105]
图5是示意地示出本公开的实施例的上行干扰定位装置2000的构成的例子的框图。
[0106]
在一些实施例中，上行干扰定位装置2000可以包括处理电路2010。上行干扰定位装置2000的处理电路2010提供上行干扰定位装置2000的各种功能。在一些实施例中，上行干扰定位装置2000的处理电路2010可以被配置为用于执行上行干扰定位装置2000中的上行干扰定位方法。另外，既可以在lte系统的网络管理侧的服务器中设置处理电路2010，即由网络管理侧的服务器实现上行干扰定位装置2000的上行干扰定位，也可以独立于网络管理侧的服务器而单独设置上行干扰定位装置2000。
[0107]
处理电路2010可以指在计算系统中执行功能的数字电路系统、模拟电路系统或混合信号(模拟和数字的组合)电路系统的各种实现。处理电路可以包括例如诸如集成电路(integrated circuit，ic)、专用集成电路(application specific integrated circuit，asic)这样的电路、单独处理器核心的部分或电路、整个处理器核心、单独的处理器、诸如现场可编程门阵列(field programmable gate array，fpga)的可编程硬件设备、和/或包括多个处理器的系统。
[0108]
在一些实施例中，处理电路2010可以包括：mr数据收集单元2020、ripprb统计单元2030、波形向量识别单元2040以及干扰位置可视化单元2050。其中，mr数据收集单元2020被配置为执行后述图6的流程图中的步骤s2002，ripprb统计单元2030被配置为执行后述图6的流程图中的步骤s2003，波形向量识别单元2040被配置为执行后述图6的流程图中的步骤s2004，干扰位置可视化单元2050被配置为执行后述图6的流程图中的步骤s2005以及图8的步骤s3001～s3003的处理。
[0109]
在一些实施例中，上行干扰定位装置2000还可以包括存储器(未图示)。上行干扰
定位装置2000的存储器可以存储由处理电路2010产生的信息以及用于上行干扰定位装置2000操作的程序和数据。存储器可以是易失性存储器和/或非易失性存储器。例如，存储器可以包括但不限于随机存取存储器(random access memory，ram)、动态随机存取存储器(dynamic random access memory，dram)、静态随机存取存储器(static random-access memory，sram)、只读存储器(read-only memory，rom)以及闪存存储器。
[0110]
另外，上行干扰定位装置2000可以以芯片级来实现，或者也可以通过包括其它外部部件而以设备级来实现。
[0111]
应当理解，上述mr数据收集单元2020、ripprb统计单元2030、波形向量识别单元2040以及干扰位置可视化单元2050仅仅是根据其所实现的具体功能所划分的逻辑模块，而不是用于限制具体的实现方式。在实际实现时，上述各个单元可被实现为独立的物理实体，或者也可由单个实体(例如，处理器(cpu或dsp等)、集成电路等)来实现。
[0112]
图6是示意地示出本公开的实施例的上行干扰定位装置2000所执行的检测处理的例示性的流程图。
[0113]
如图6的步骤s2001所示，该上行干扰定位处理被定期触发或者收到来自用户终端装置的投诉等时被触发。虽未图示，但应当理解的是，在定期触发上行干扰定位处理的情况下，在步骤s2001中反复判定预定的定时是否到来，直到判定为“是”而进入步骤s2002。
[0114]
在步骤s2002中，mr数据收集单元2020从用户终端装置收集mr数据。
[0115]
然后，在步骤s2003中，ripprb统计单元2030按照预定的统计周期，例如5分钟、15分钟等，根据mr数据，对用户终端装置ue所在的lte小区的对应的例如100个prb统计prb粒度enb接收干扰功率(ripprb)的样本值，生成如图7所示的数字化干扰波形。图7中示出基于ripprb的数字化干扰波形，横坐标为prb的编号，纵坐标为各个prb的上行干扰功率值。根据图7可知，该lte小区对应的三个扇区的每个扇区的各扇区对应的上行干扰功率各不相同。
[0116]
在基于快速模板匹配实现的波形识别中，针对单个波形的识别比较简单，但实际应用中，多个波形复合在一起的情况更加常见。如上所述，在lte系统化境中，上行干扰可能同时包括系统内干扰和系统外干扰，若能够迅速精确地定位排查出系统外干扰的干扰源，则能够大幅提高lte系统的运行稳定性，提高用户的体验度。
[0117]
为此，在本公开中，在步骤s2004中，波形向量识别单元2040首先将步骤2003中得到的数字化干扰波形向量化，得到复合波形向量，从复合波形向量中抽取出特征已知的波形向量(如热噪声波形向量、系统内干扰波形向量等)，然后，利用上述说明的模板匹配波形识别算法，例如，基于系统外干扰的波形向量模板，对抽取特征已知的波形向量后的剩余干扰波形向量进行波形匹配识别。
[0118]
在一些实施例中，例如，将let系统统计得到的波形向量设为xa(xa1，xa2,
…
xam)，将系统的热噪声波形向量设为xb(xb1，xb2,
…
xbm)，将阻塞干扰向量设为xc(xc1，xc2，
…
xcm)，计算xa-xb-xc＝xd，然后对xd与系统外干扰的波形向量模板(例如，外部干扰源干扰波形模板、杂散干扰波形模板)进行匹配识别，判定是否存在系统外干扰。
[0119]
接下来，在通过步骤s2004的波形识别判定为存在系统外干扰的情况下，在步骤s2005中，根据mr数据对系统外干扰的干扰源的位置进行可视化呈现。
[0120]
接下来，结合图8，具体说明对系统外干扰的干扰源的位置进行可视化呈现的处理的一个例子的流程图。
[0121]
如图8所示，在判定为存在系统外干扰的情况下，在步骤s3001中，判定mr数据中的“a_gps”是否为空。在“a_gps”非空的情况下，直接进入步骤s3003，根据“a_gps”的经纬度信息对地图栅格进行渲染。另一方面，在“a_gps”为空的情况下，进入步骤s3002的处理。
[0122]
在步骤s3002中，根据mr数据对用户终端装置的地理坐标进行计算，确定用户终端装置的经纬度坐标。
[0123]
在一些实施例中，首先确定mr数据对应的服务小区enbid，通过基站工参查询，得到基站的经纬度坐标，设经度为long0、纬度为lat0。然后，查询mr数据中的aoa查询表，取得对应的天线到达角aoa_angle平均值，设为a(度)。再通过mr中的tadv查询表，取的对应的天线距离设为d(米)，d表示基站天线到达用户终端ue的直线距离，由于基站一般建设在铁塔或楼顶等具有一定高度的场所，对于平面地图而言，需要计算的是平面投影的距离，这里可以取均值d/1.5。
[0124]
然后，根据公式4、公式5求出用户终端ue的地理坐标long1和lat1：
[0125][0126][0127]
接着，在步骤s3002中确定出用户终端装置的经纬度坐标之后，进入步骤s3003，对地图栅格进行渲染。
[0128]
在一些实施例中，在对地图栅格进行渲染时，将计算得到的用户终端装置的地理位置(经纬度)与地图栅格的经纬度进行映射，通过例如gis栅格仿真定位算法，将当前用户的位置信息(经纬度)定位到地图栅格上，从而得到一个待渲染的栅格。在本公开中，关于地图栅格化的方法没有特别限定，可以利用公知的方法进行。另外，由于mr数据是lte系统中的用户终端装置在执行业务过程中上报给无线网络管理侧的测量数据，其能够准确地反映出用户坐在区域的网络覆盖以及受到的干扰情况，因此，例如在以25m*25m的精度对地图进行栅格化的情况下，通过将海量用户级别的mr数据进行栅格聚类，能够直观地呈现网络的干扰情况。
[0129]
在一些实施例中，在将用户终端装置的地里位置与栅格位置进行映射后，结合当前mr测量记录的receivedipower平均值，生成栅格映射表，如下表6所示：
[0130]
表6上行干扰信号栅格映射表
[0131][0132]
然后，基于表6中的“栅格receivedipower平均值”，对用户终端装置的当前位置对应的栅格进行渲染。这里，所谓“对栅格进行渲染”是指以栅格为对象而呈现某个值。在本公开中，关于渲染处理的实现方法没有特别限定，可以通过公知的方法进行。在本公开中，将栅格receivedipower平均值渲染到栅格中。
[0133]
在一些实施例中，根据“栅格receivedipower平均值”，将上行干扰信号强度分为5
个区间，对应于5个干扰等级。例如，预先设定如图9所示的receivedipower上行干扰信号强度区间颜色配置表。在表中，1级干扰(《-116dbm)的栅格为空白显示，2级干扰(》＝-116dbm且《-110dbm)的栅格用蓝色(图10中示出为)渲染，3级干扰(》＝-110dbm且《-100dbm)的栅格用绿色(图10中示出为)渲染，4级干扰(》＝-100dbm且《-90dbm)的栅格用黄色渲染(图10中示出为)，5级干扰(》＝-90dbm且《-50dbm)的栅格用红色(图10中示出为)渲染。通过根据颜色配置表对相应的栅格进行颜色渲染，能够将上行干扰信号的栅格地理化为图10所示的栅格颜色渲染示意图。在图10中，直观地呈现出上行干扰源所处的位置(即3、4、5级干扰集中的区域，在地图栅格中表现为显眼的绿、黄、红色)。
[0134]
另外，在一些实施例中，并非通过不同的颜色，而是预先根据干扰等级设置不同的闪烁频率，例如，干扰等级越高，则设置为闪烁频率越快，干扰等级越低，则设置为闪烁频率越慢。这样，在对栅格进行了渲染时，同样能够根据不同栅格的闪烁频率迅速且准确地定位出干扰源所处的位置。
[0135]
进而，在一些实施中，并非通过不同的颜色，而是通过同一颜色的不同深浅对栅格进行渲染，例如，干扰等级越高则渲染为颜色越深。由此，同样能够根据渲染后的栅格的颜色深浅迅速准确地定位出干扰源所处的位置。
[0136]
关于渲染的表现方式，可以根据需要采用各种方式，只要能够直观地、能够明显区分地表现上行干扰的强度等级即可。
[0137]
根据本公开的上行干扰定位装置以及上行干扰定位方法，根据mr数据生成上行干扰数字波形，将数字波形向量化，并从复合在一起的数字波形向量中的特征已知的例如系统内干扰波形向量去除，实现了不同干扰类型特征的波形的分离，从而能够快速有效定位出受干扰小区的主要干扰类型。另外，通过mr数据的栅格仿真技术，对受干扰用户进行栅格聚类，并gis地理化，能够根据栅格颜色的变化、闪烁频率或栅格颜色深浅等，直观地呈现上行干扰源的位置，从而能够迅速准确地进行干扰排查，提升网络稳定性，提高了用户的体验度。
[0138]
另外，根据本公开上行干扰定位装置以及上行干扰定位方法，不仅根据mr数据中的receivedipower判定lte小区的上行干扰强度，还基于receivedipower的平均值对地图栅格进行了渲染，从而能够简单地将表示上行干扰强度的数据灵活运用到对于干扰源的定位中，提高了对系统外干扰的排查效率。
[0139]
以下，对一个具体的实施例进行说明。如图11所示，某运营商接到多起用户投诉白云区石头村区域的上网网速慢和volte掉话的情况，后台分析人员提取网管数据，发现682566_1_白云北二环俊发城下拉槽_拉远_电_电_c小区prb轮询曲线干扰噪声值均在-95dbm左右，经排查基站无故障，怀疑存在系统外干扰。然后，提取该站mr数据，进行仿真分析。根据栅格渲染图，很快发现石头村西北和西南方向两个位置(图11中用圆圈出的区域)存在较强的干扰值，结合用户投诉的位置信息，网络优化中心立即安排优化人员到上述两处怀疑区域进行现场排查。排查人员提供过对石头村的制高点位置进行扫频，很快定位出两组干扰源，均为用户装放大器导致的干扰，协调用户关闭放大器后，扫频仪干扰尖峰波形消失，后台网管指标亦回复正常。
[0140]
在上述干扰排查过程中，实际排查到两处干扰源位置与mr数据渲染位置大体一致，通过mr数据干扰定位识别方法，结合人工干扰排查，快速地定位出干扰位置信息，极大
地缩小了人工排查范围，有效提升了干扰排查效率，改善用户感知。
[0141]
接下来，图12示出了可以实现根据本公开的实施例的计算设备300的示例性配置。计算设备300是可以应用本公开的上述方面的硬件设备的实例。计算设备300可以是被配置为执行处理和/或计算的任何机器。计算设备300可以是但不限制于工作站、服务器、台式计算机、膝上型计算机、平板计算机、个人数据助手(pda)、智能电话、车载计算机或以上组合。
[0142]
如图12所示，计算设备300可以包括可以经由一个或多个接口与总线302连接或通信的一个或多个元件。总线302可以包括但不限于，工业标准架构(industry standard architecture，isa)总线、微通道架构(micro channel architecture，mca)总线、增强isa(eisa)总线、视频电子标准协会(vesa)局部总线、以及外设组件互连(pci)总线等。计算设备300可以包括例如一个或多个处理器304、一个或多个输入设备306、以及一个或多个输出设备308。一个或多个处理器304可以是任何种类的处理器，并且可以包括但不限于一个或多个通用处理器或专用处理器(诸如专用处理芯片)。各处理器304例如可以分别对应于图2中的处理电路2010，被配置为能够实现本公开的上行干扰定位装置2000的功能。输入设备306可以是能够向计算设备输入信息的任何类型的输入设备，并且可以包括但不限于鼠标、键盘、触摸屏、麦克风和/或远程控制器。输出设备308可以是能够呈现信息的任何类型的设备，并且可以包括但不限于显示器、扬声器、视频/音频输出终端、振动器和/或打印机。
[0143]
计算设备300还可以包括或被连接至非暂态存储设备314，该非暂态存储设备314可以是任何非暂态的并且可以实现数据存储的存储设备，并且可以包括但不限于盘驱动器、光存储设备、固态存储器、软盘、柔性盘、硬盘、磁带或任何其他磁性介质、压缩盘或任何其他光学介质、缓存存储器和/或任何其他存储芯片或模块、和/或计算机可以从其中读取数据、指令和/或代码的其他任何介质。计算设备300还可以包括随机存取存储器(ram)310和只读存储器(rom)312。rom 312可以以非易失性方式存储待执行的程序、实用程序或进程。ram 310可提供易失性数据存储，并存储与计算设备300的操作相关的指令。计算设备300还可包括耦接至数据链路318的网络/总线接口316。网络/总线接口316可以是能够启用与外部装置和/或网络通信的任何种类的设备或系统，并且可以包括但不限于调制解调器、网络卡、红外线通信设备、无线通信设备和/或芯片集(诸如蓝牙
tm
设备、802.11设备、wifi设备、wimax设备、蜂窝通信设施等)。
[0144]
应当理解，本说明书中“实施例”或类似表达方式的引用是指结合该实施例所述的特定特征、结构、或特性系包括在本公开的至少一具体实施例中。因此，在本说明书中，“在本公开的实施例中”及类似表达方式的用语的出现未必指相同的实施例。
[0145]
本领域技术人员应当知道，本公开被实施为一系统、装置、方法或作为计算机程序产品的计算机可读媒体(例如非瞬态存储介质)。因此，本公开可以实施为各种形式，例如完全的硬件实施例、完全的软件实施例(包括固件、常驻软件、微程序代码等)，或者也可实施为软件与硬件的实施形式，在以下会被称为“电路”、“模块”或“系统”。此外，本公开也可以任何有形的媒体形式实施为计算机程序产品，其具有计算机可使用程序代码存储于其上。
[0146]
本公开的相关叙述参照根据本公开具体实施例的系统、装置、方法及计算机程序产品的流程图和/或框图来进行说明。可以理解每一个流程图和/或框图中的每一个块，以及流程图和/或框图中的块的任何组合，可以使用计算机程序指令来实施。这些计算机程序指令可供通用型计算机或特殊计算机的处理器或其它可编程数据处理装置所组成的机器
来执行，而指令经由计算机或其它可编程数据处理装置处理以便实施流程图和/或框图中所说明的功能或操作。
[0147]
在附图中显示根据本公开各种实施例的系统、装置、方法及计算机程序产品可实施的架构、功能及操作的流程图及框图。因此，流程图或框图中的每个块可表示一模块、区段、或部分的程序代码，其包括一个或多个可执行指令，以实施指定的逻辑功能。另外应当注意，在某些其它的实施例中，块所述的功能可以不按图中所示的顺序进行。举例来说，两个图示相连接的块事实上也可以同时执行，或根据所涉及的功能在某些情况下也可以按图标相反的顺序执行。此外还需注意，每个框图和/或流程图的块，以及框图和/或流程图中块的组合，可藉由基于专用硬件的系统来实施，或者藉由专用硬件与计算机指令的组合，来执行特定的功能或操作。
[0148]
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。