[发明专利]一种基于RSSI测距的空气质量监测点网络定位的方法

摘要

本发明公开了一种基于RSSI测距的空气质量监测点网络定位的方法，步骤为Ⅰ、获取RSSI值；Ⅱ、建立RSSI测距模型；Ⅲ、RSSI测距初步估算；Ⅳ、基于RSSI的测量距离矫正；Ⅴ、目标监测节点的定位；其应用RF无线通信模型，无需求出路径衰减指数，推导出简便的距离估算方法，减少了环境因素对测距的不良影响，选取相对距离误差系数进行测距校正，获得精度较高的距离信息，定位估算部分采用递推最小二乘的定位算法估计未知监测节点的坐标。仿真试验表明，本算法将非线性观测方程线性化，从而采用线性方法求解，既降低了位置计算复杂性，充分考虑已有的网络先验信息的同时又达到了良好的定位效果。

主权项

一种基于RSSI测距的空气质量监测点网络定位的方法，在多个监测节点中对部分节点的位置进行标定，标定已知的节点为锚节点，其特征在于，该方法的具体步骤为：Ⅰ、获取RSSI值：通过无线通信获取各个监测节点的RSSI值；Ⅱ、建立RSSI测距模型：采用对数正态分布无线传播模型建立RSSI测距模型，其表达式为：Pr(d)=Pt-PL(d0)-10nlg(dd0)+Nσ---(1)]]>其中d0为信号传播参考距离，且d0＝10，单位为m；d为节点之间的测量距离，单位为m；Pr(d)是距离为d时节点接收到的信号强度，单位为dBm；n是与监测环境有关的路径损耗因子；Pt代表节点无线信号发射强度，单位为dBm；测量噪声Nσ为服从的高斯随机变量；PL(d0)为无线信号经传播距离d0后的信号强度损耗，单位为dB，可由信号传播的自由空间模型获得，即PL(d0)=-10lg[GtGrλ2(4π)2d02L]---(2)]]>其中Gt是节点发射天线增益，单位为dBi；Gr为节点接收天线增益，单位为dBi；L是系统损耗系数；λ是RF信号波长，单位为m；在网络中节点zi通信范围内，设其邻居节点zj到zi的测量距离为d，则zi处的无线信号接收强度为Pr(d)，令Pr(d0)＝Pt‑PL(d0)   (3)由式(1)可知Pr(d)=Pr(d0)-10nlg(dd0)+Nσ---(4)]]>其中，d>d0，Pr(d0)是参考距离d0对应的无线信号接收强度；Ⅲ、RSSI测距初步估算：RSSI测距的距离估算的精度直接决定其定位精度；由式(4)可写为Pr(d)＝PN‑10nlgd   (5)其中PN＝Pr(d0)‑Nσ   (6)则有10nlgd＝PN‑Pr(d)   (7)设定无线网络中的监测节点在监测区域内均匀随机部署，且监测节点具有统一相同的通信半径，考虑通信范围是规则的圆形区域；若某个监测节点的通信范围内均匀随机部署多个的其他监测节点，可依据Pr(d)与测量距离d成反比，认为监测节点接收的最小信号强度Pmin所对应的距离最大即dmax，因此有10nlgdmax＝PN‑Pmin   (8)先由式(8)求出n为n=PN-Pmin10lgdmax---(9)]]>再代入式(7)中，求出估算距离为d=dmaxPN-Pr(d)PN-Pmin≈rPN-Pr(d)PN-Pmin---(10)]]>由此可见，依据式(10)就可以计算出任一监测节点zi与其邻居节点间的测量距离；在未知节点获得的多个Pr(d)值中，取最小的Pr(d)值即为Pmin，对应dmax＝r，即为通信半径，从而可由式(10)获得未知节点到锚节点的d值；Ⅳ、基于RSSI的测量距离矫正：首先对锚节点间的测量距离与真实距离进行比较，获得由Pr(d)计算得到的测距误差，然后对监测节点与锚节点间的测量距离进行校正以减少随机测距误差；设网络中锚节点Ai(xi，yi)至其他锚节点的实际距离为rik，k＝1，2，…，n，k≠i；其中n为网络中锚节点个数；锚节点Ai至其他锚节点的测量距离为dik，k＝1，2，…，n，k≠i；则有锚节点Ai处实际距离与测量距离的距离平均相对误差系数μi为μi=1n-1Σk=1,k≠inrik-dikdik---(11)]]>其中n为监测区域内锚节点Ai的邻居锚节点个数；μi反映锚节点Ai处RSSI测距的精确程度，称为相对测量误差系数；事先设定好阈值τ，若|μi|≥τ，则测距误差过大，锚节点Ai被认作状况异常锚节点，不参与定位，否则参与定位计算；监测节点接收到锚节点的广播消息，即将接收到的RSSI值通过式(10)计算，得到监测节点与锚节点之间的测量距离d，再利用锚节点的相对误差系数，即依据式(12)对测量距离d进行测距校正duic=dui(1+μi)---(12)]]>其中dui是监测节点和锚节点Ai之间的测量距离，单位为m；是监测节点和锚节点Ai之间的校正距离，单位为m；μi为锚节点Ai的相对误差系数；Ⅴ、目标监测节点的定位：设网络中N个锚节点记为Ai(xi，yi)，i＝1，2，…，N；各个锚节点到未知位置的盲节点B(x1，y1)的测量距离为di，i＝1，2，…，N，即有N个测量方程为d^i=di+vi,i=1,2,...,N---(13)]]>其中di=(xi-x)2+(yi-y)2]]>vi为测量误差；式(13)两边平方，有：x2+y2-2xix-2yiy=(d^l-vi)2-(xi2+yi2)]]>采用批处理数据方式的最小二乘整理，可简写为h＝Gθ+v   (14)其中，G=-2x1-2y11.........-2xN-2yN1]]>h=d^12-(x12+y12),...,d^N2-(xN2+yN2)T]]>v=2d^1v1-v12,...,2d^NvN-vN2T]]>θ＝[x y R]TR＝x2+y2则LS解为：θ^ls=(GTG)-1GTh---(15)]]>则WLS解为：θ^wls=(GTWG)-1GTWh---(16)]]>其中W设置为单位阵；最终通过得出的估计值确定实时监测节点的位置。