基于控制图的干扰源数量估计方法
【技术领域】
[0001] 本发明属于网络数据通信技术领域，特别是一种基于控制图的干扰源数量估计方 法。
【背景技术】
[0002] 干扰源定位是指多跳无线网络中的被干扰节点，通过主动测量和被动监听，得到 多个协议找层次的观察结果，利用无线信道传播特性和干扰区域几何知识，协作推断出干 扰源的相对或绝对位置。干扰源定位的典型算法包括质屯、定位算法、权重质屯、定位算法、虚 拟力迭代算法、凸壳定位方法等。运些定位算法大多针对单个干扰源提出。为了达到高效、 大范围的干扰，多个干扰源通常协同工作，W增强干扰效果，从而使得无线网络节点周围可 能存在多个干扰源。为此，如何有效的定位多个干扰源使得问题进一步复杂化。
[0003] 北京科技大学的程天被等人提出了一种称为X射线的干扰区域定位算法。该X射 线定位算法包括干扰区域映射、干扰区域骨架化和干扰源位置确定=个步骤。该算法通过 干扰区域骨架化，利用分叉点的位置信息反映干扰源的物理位置，由于区域边界受噪声影 响较大，因此X射线干扰区域定位算法定位精度偏低。
[0004] 为了准确定位多个干扰源，对干扰源数量的估计是基础。现有技术中未见关于多 个干扰源数量估计的方法，导致难W对多个干扰源进行准确定位。
【发明内容】
阳〇化]本发明的目的在于提供一种基于控制图的干扰源数量估计方法，能准确估计多个 干扰源的数量，从而有利于对多个干扰源进行准确定位。
[0006] 实现本发明目的的技术解决方案为：一种基于控制图的干扰源数量估计方法，包 括如下步骤：
[0007] (10)邻居监测：本节点接收来自邻居节点的报文，根据所述报文建立本节点的邻 居列表，记录每个邻居节点的信噪比，本节点与邻居节点之间通过交互报文获得彼此的邻 居列表；
[0008] (20)干扰判别：根据本节点的邻居列表及邻居节点的信噪比变化情况，判断本节 点是否受到干扰攻击，如否，则返回（10)邻居监测步骤；
[0009] (30)干扰估计：如本节点受到干扰攻击，则其向邻居节点发送静默请求，监测干 扰信号强度，形成干扰信号强度矩阵，采用控制图方法从干扰信号强度矩阵估计出干扰源 数量。
[0010] 本发明与现有技术相比，其显著优点为：
[0011] 1.开销小，实现简单。本发明提出的干扰源数量估计方法无需复杂的通信交互和 信号处理，仅需依赖廉价无线设备所具有的接收信号强度监测功能进行干扰源信号强度监 ，无需进行干扰信号分离，因而处理开销大大降低，便于在普通民用无线设备上实现。
[0012] 2.精度高，估计准确。本发明提出使用控制图方法对干扰信号强度矩阵进行变更 点检测，可W精确发现干扰信号强度变化较大的变更时刻，可W准确估计干扰源数量的变 化情况。
[0013] 下面结合附图和【具体实施方式】对本发明作进一步的详细描述。
【附图说明】
[0014] 图1是本发明基于控制图的干扰源数量估计方法的主流程图。
[0015] 图2是图1中干扰估计步骤的流程图。
[0016] 图3是存在干扰源的多跳无线网络场景示例图。
[0017] 图4是采用本发明基于控制图的干扰源数量估计方法对多个干扰源的估计的结 果不例。
【具体实施方式】
[0018] 如图1所示，本发明基于控制图的干扰源数量估计方法，包括如下步骤：
[0019] (10)邻居监测：本节点接收来自邻居节点的报文，根据运些报文建立本节点的邻 居列表，记录每个邻居节点的信噪比，然后，本节点与邻居节点之间通过交互报文获得彼此 的邻居列表；
[0020] 例如，本节点与邻居节点之间通过交互肥IGHBOR_LIST_REQ肥ST和肥IGHB0R_ LIST_REPLY报文，获得彼此的邻居列表；当邻居节点收到肥IGHBOR_LIST_REQUEST报文时， 邻居节点将自身的邻居列表放入肥IGHBOR_LIST_REPLY报文，回复本节点。
[0021] (20)干扰判别：根据本节点的邻居列表及邻居节点的信噪比变化情况，判断本节 点是否受到干扰攻击，如否，则返回（10)邻居监测步骤；
[0022] 所述（20)干扰判别步骤具体为：当本节点的邻居数量减少幅度大于等于30%，或 者尽管本节点邻居数量减少幅度小于30%，但有不少于30%的邻居节点的信噪比降幅大 于等于30%，则判定本节点受到干扰攻击。
[0023] (30)干扰估计：如本节点受到干扰攻击，则本节点向邻居节点发送静默请求，监 测干扰信号强度，形成干扰信号强度矩阵，采用控制图方法从干扰信号强度矩阵估计出干 扰源数量。
[0024] 如图2所示，所述（30)干扰估计步骤包括：
[00对 （31)静默请求发送：本节点向邻居节点发送报文，请求邻居节点在静默期内不向 本节点发送信息； 阳0%] 例如，本节点向邻居节点发送静默请求报文SILENCE_REQUEST。SILENCE_REQUEST 报文中规定了静默时长。在静默时长内，邻居节点将不会向节点发送信息。
[0027] (32)干扰信号强度监测：本节点在m个信道上，经过n个监测周期t，记录自身接 收到的干扰信号强度，形成干扰信号强度矩阵
[00巧]X= [Xi,而，…Xi,…，Xm]，
[0029]其中，Xi表示第i个信道上监测到的干扰信号强度时间序列， 阳030] Xi=Ix…枯，…，Xji,…，XnJ，
[00川其中，X，康示第i个信道上第j个监测周期内监测到的干扰信号强度；
[0032] (33)干扰信号数量估计：采用控制图方法，对干扰信号强度矩阵进行处理，估计 出干扰源数量；
[0033]所述（33)干扰信号数量估计步骤包括： 阳034] (331)构建n-1 行m列的矩阵Y= [yi，y2,…，y"J，使得y。= X(w)j-Xij; 阳0对 （33?对矩阵Y应用Shewhart充控制图方法，其上控制线为
[0037]下控制线为
[0039]其中，
[0041] 度为Y的各列极差的平均值，dz是依赖于参数n-1的一个常数；
[0042] (333)构建包含n-1个元素的向量 阳043] Z= {Zi,Z2,…，ZnJ， W44] 其中，
[0046] (334)设置第1个监测周期之前的干扰源数量为0 ;
[0047] (335)估计第i个监测周期之前的干扰源数量：从Zi开始到Z11为止，在保证干