一种基于网格化频谱监测结果的无线干扰源定位方法技术

本发明专利技术涉及一种基于网格化频谱监测结果的无线干扰源定位方法，属于无线电信号源定位领域。本发明专利技术提出的无线电定位方法区别于传统的测向定位法，首先在被监测区域设有若干个矩阵排列的无线电监测器，通过对于无线电监测器接受功率的大小进行分析，测出其接受功率的大小，根据无线电监测器检测的数据，选出干扰源所在的一个小矩形区域，通过一定的算法，求得干扰源可能所在的区域。

【技术实现步骤摘要】

本专利技术涉及一种基于网格化频谱监测结果的无线干扰源定位方法，属于无线电信号源定位领域。
技术介绍
20世纪40年代开始出现双曲线无线电定位理论，人们开始就一些简单的侧向设备来围绕目标信号源进行多次测向，利用人工作图的方式确定出信号源的位置。到了20世纪80年代末，美国产生精确辐射源定位系统理论，在此基础上产生了时差定位理论、频差定位理论。目前主要采用的无源定位技术还是测向定位技术和时差定位技术。无源定位技术在我国发展较晚，从20世纪80年代初期才开始出现这方面的理论研究。目前我国也将无源定位技术的研究和使用作为电子领域的一项重要技术。近年来，随着无线电现代通信技术的进步和信息产业的兴起，频谱资源日益稀缺，所以频谱检测系统对现在的频谱频谱监管越来越重要，加强对于干扰源或非法电台的查出力度势在必行。传统的定位方法中大多采用多站测向定位、无源测向定位、基于方位角和网络拓扑的定位等，在传统的定位方法中，需要测出干扰源的方向，到达时间差等等，操作复杂，不够直观，而本专利技术提供了一种不需要测出干扰源具体方向，到达时间差等，只根据接收到的功率大小，就可以很直观的定位干扰源所在区域。
技术实现思路
本专利技术提供了一种基于网格化频谱监测结果的无线干扰源定位方法，以用于解决在不知道干扰源方向的情况下，就可以很直观的定位干扰源所在区域。本专利技术的技术方案是：一种基于网格化频谱监测结果的无线干扰源定位方法，在一片待测区域，设有若干个矩阵排列的无线电监测器，构成一个无线电监测网，用于监测该区域的干扰源的信号强度，测出其接收功率的大小，通过对于无线电监测器接收功率的大小进行分析，根据无线电监测器检测的数据，找出四个数据最大的监测点，即为干扰源所在的一个小矩形区域。通过对于数据的分析可知，干扰源所在的这个小矩形区域的四个顶点分别为a1、a2、a3、a4，它们实际对应的接收功率大小Pn1、Pn2、Pn3、Pn4，根据矩形区域的四个点的坐标及接收功率的大小，通过一定的算法便可推出干扰源位置所在的区域。所述基于网格化频谱监测结果的无线干扰源定位方法的具体步骤如下：Step1、首先确定电磁波在二维自由空间的传播模型如下： P r = A r 4 πd 2 P t G t ]]>式中Pt为发射功率，以圆辐射，λ为工作波长，Gt、Gr分别表示发射天线和接受天线增益，d为发射天线和接受天线间的距离；在二维自由空间中，接收功率Pr与发射天线和接受天线间的距离d2成反比，Gt、Gr、Ar均为已知常量，则电磁波在二维自由空间的传播模型可简化为如下： P r = k d 2 P t ]]>其中Step2、在一片待测区域，设有若干个矩阵排列的无线电监测器，构成一个无线电监测网，用于监测该区域的干扰源的信号强度，测出其接收功率的大小，通过对于无线电监测器接收功率的大小进行分析，根据无线电监测器检测的数据，找出相邻四个监测点数据之和最大的小矩形区域，即为干扰源所在的一个小矩形区域。Step3、设干扰源坐标为(X0，Y0)，发射功率为Pt，干扰源所在矩形区域的四个监测点坐标分别为a1(X1，Y1)、a2(X2，Y2)、a3(X3，Y3)、a4(X4，Y4)，在没有噪声干扰的理想环境下，对应监测点测得的接收功率为Pr1、Pr2、Pr3、Pr4，其中X1＝X4，X2＝X3，Y1＝Y2，Y3＝Y4，则X1<X0<X2，Y1<Y0<Y4，故根据电磁波在二维自由空间的传播模型可得： ( X 1 - X 0 ) 2 + ( Y 1 - Y 0 ) 2 = kP t P r 1 ]]> ( X 2 - X 0 ) 2 + ( Y 2 - Y 0 ) 2 = kP t P r 2 ]]> ( X 3 - X 0 ) 2 + ( Y 3 - Y 0 ) 2 = kP t P r 3 ]]> ( X 4 - X 0 ) 2 + ( 本文档来自技高网...

【技术保护点】
一种基于网格化频谱监测结果的无线干扰源定位方法，其特征在于：在待测区域，设有若干个矩阵排列的无线电监测器，通过对于无线电监测器接收功率的大小进行分析，测出其接收功率的大小，根据无线电监测器检测的数据，选取干扰源所在的一个小矩形区域，进一步计算求得干扰源可能所在的区域；具体步骤如下：Step1、在一片待测区域，设置若干个矩阵排列的无线电监测器，构成一个无线电监测网，根据无线电监测器检测的接收功率数据找出相邻四个监测点数据之和最大的矩形区域，即为干扰源所在的一个小矩形区域；Step2、设干扰源坐标为(X0，Y0)，发射功率为Pt，干扰源所在矩形区域的四个监测点坐标分别为a1(X1，Y1)、a2(X2，Y2)、a3(X3，Y3)、a4(X4，Y4)，在没有噪声干扰的理想环境下，对应监测点测得的接收功率为Pr1、Pr2、Pr3、Pr4，其中X1＝X4，X2＝X3，Y1＝Y2，Y3＝Y4，则X1<X0<X2，Y1<Y0<Y4；在没有噪声干扰的理想环境下，根据电磁波在二维自由空间的传播模型，可以得到对应监测点测得的接收功率Pr1、Pr2、Pr3、Pr4满足以下等式：(X1-X0)2+(Y1-Y0)2=kPtPr1]]>(X2-X0)2+(Y2-Y0)2=kPtPr2]]>(X3-X0)2+(Y3-Y0)2=kPtPr3]]>(X4-X0)2+(Y4-Y0)2=kPtPr4]]>进一步可得：Pr2Pr1[(X2-X0)2+(Y2-Y0)2]-(X1-X0)2-(Y1-Y0)2=0]]>Pr2Pr3[(X2-X0)2+(Y2-Y0)2]-(X3-X0)2-(Y3-Y0)2=0]]>Pr2Pr4[(X2-X0)2+(Y2-Y0)2]-(X4-X0)2-(Y4-Y0)2=0]]>Pr3Pr1[(X3-X0)2+(Y3-Y0)2]-(X1-X0)2-(Y1-Y0)2=0]]>Pr3Pr4[(X3-X0)2+(Y3-Y0)2]-(X4-X0)2-(Y4-Y0)2=0]]>Pr4Pr1[(X4-X0)2+(Y4-Y0)2]-(X1-X0)2-(Y1-Y0)2=0]]>Step3、设电波在空中传播受到高斯白噪声的影响，高斯白噪声的方差为c，均值为u，则高斯白噪声的数学模型可以表示为u+c*randn()，则在实际中对应监测点测得的接收功率分别为Pn1、Pn2、Pn3、Pn4，其中Pn1＝Pr1+u+c*randn()，Pn2＝Pr2+u+c*randn()，Pn3＝Pr3+u+c*randn()，Pn4＝Pr4+u+c*randn()。；故由电磁波在含有噪声的二维空间传播可得：|Pn2Pn1[(X2-X0)2+(Y2-Y0)2]-(X1-X0)2-(Y1-Y0)2|<d]]>|Pn2Pn3[(X2-X0)2+(Y2-Y0)2]-(X3-X0)2-(Y3-Y0)2|<d]]>|Pn2Pn4[(X2-X0)2+(Y2-Y0)2]-(X4-X0)2-(Y4-Y0)2|<d]]>|Pn3Pn1[(X3-X0)2+(Y3-Y0)2]-(X1-X0)2-(Y1-Y0)2|<d]]>|Pn3Pn4[(X3-X0)2+(Y3-Y0)2]-(X4-X0)2-(Y4-Y0)2|<d]]>|Pn4Pn1[(X4-X0)2+(Y4-Y0)2]-(X1-X0)2-(Y1-Y0)2|<d]]>其中d为每个式子的误差值；Step4、设定合适的误差值d，X0和Y0在小矩形区域内遍历取值，标记出能同时满足Step3中所有式子的坐标点，这些点所在的区域即为干扰源所在的区域。...

【技术特征摘要】
1.一种基于网格化频谱监测结果的无线干扰源定位方法，其特征在于：在待测区域，设有若干个矩阵排列的无线电监测器，通过对于无线电监测器接收功率的大小进行分析，测出其接收功率的大小，根据无线电监测器检测的数据，选取干扰源所在的一个小矩形区域，进一步计算求得干扰源可能所在的区域；具体步骤如下：Step1、在一片待测区域，设置若干个矩阵排列的无线电监测器，构成一个无线电监测网，根据无线电监测器检测的接收功率数据找出相邻四个监测点数据之和最大的矩形区域，即为干扰源所在的一个小矩形区域；Step2、设干扰源坐标为(X0，Y0)，发射功率为Pt，干扰源所在矩形区域的四个监测点坐标分别为a1(X1，Y1)、a2(X2，Y2)、a3(X3，Y3)、a4(X4，Y4)，在没有噪声干扰的理想环境下，对应监测点测得的接收功率为Pr1、Pr2、Pr3、Pr4，其中X1＝X4，X2＝X3，Y1＝Y2，Y3＝Y4，则X1<X0<X2，Y1<Y0<Y4；在没有噪声干扰的理想环境下，根据电磁波在二维自由空间的传播模型，可以得到对应监测点测得的接收功率Pr1、Pr2、Pr3、Pr4满足以下等式： ( X 1 - X 0 ) 2 + ( Y 1 - Y 0 ) 2 = kP t P r 1 ]]> ( X 2 - X 0 ) 2 + ( Y 2 - Y 0 ) 2 = kP t P r 2 ]]> ( X 3 - X 0 ) 2 + ( Y 3 - Y 0 ) 2 = kP t P r 3 ]]> ( X 4 - X 0 ) 2 + ( Y 4 - Y 0 ) 2 = kP t P r 4 ]]>进一步可得： P r 2 P r 1 [ ( X 2 - X 0 ) 2 + ( Y 2 - Y 0 ) 2 ] - ( X 1 - X 0 ) 2 - ( Y 1 - Y 0 ) 2 = 0 ]]> P r 2 P r 3 [ ( X 2 - X 0 ) 2 + ( Y 2 - Y 0 ) 2 ] - ( X 3 - X 0 ) 2 - ( Y 3 - Y 0 ) 2 = 0 ]]> P r 2 P r 4 [ ( X 2 - X 0 ) 2 + ( Y 2 - Y 0 ) 2 ] - ( X 4 - X 0 ) 2 - ( Y 4 - Y 0 ) 2 = 0 ]]> P r 3 P r 1 [ ( X 3 - X 0 ) 2 + ( Y 3 - Y 0 ) 2 ] - ( X 1 - X 0 ) 2 - ( Y 1 - Y 0 ) 2 = 0 ]]> P r 3 P r 4 [ ( X 3 - X 0 ) 2 + ( Y 3 - Y 0 ) 2 ] - ( X 4 - X 0 ) 2 - ( Y 4 - Y 0 ) 2 = 0 ]]> P r 4 P r 1 [ ( X 4 - X 0 ) 2 + ( Y 4 - Y 0 ) 2 ] - ( X 1 - X 0 ) 2 - ( Y 1 - Y 0 ) 2 = 0 ]]>Step3、设电波在空中传播受到高斯白噪声的影响，高斯白噪声的方差为c，均值为u，则高斯白噪声的数学模型可以表示为u+c*randn()，则在实际中对应监测点测得的接收功率分别为Pn1、Pn2、Pn3、Pn4，其中Pn1＝Pr1+u+c*randn()，Pn2＝Pr2+u+c*randn()，Pn3＝Pr3+u+c*randn()，Pn4＝Pr4+u+c*randn()。；故由电磁波在含有噪声的二维空间传播可得： | P n 2 P n 1 [ ( X 2 - ...

【专利技术属性】
技术研发人员：沈韬，闫真真，李金山，
申请(专利权)人：昆明理工大学，
类型：发明
国别省市：云南;53