一种基于马鞍地形山区无线传感器网络节点的定位方法
【专利摘要】本发明公开了一种基于马鞍地形山区无线传感器网络节点的定位方法，本发明无需测距的定位机制不需测量节点间的绝对距离或方位，降低了对节点硬件的成本，体积和能量消耗，更适合于大规模传感器网络。但非测距的节点定位方法的精度较低，在理想环境下其定位精度仍可满足应用的需求，但在山区复杂地形上的应用时其定位精度难于满足实际应用的要求。本发明是在山区马鞍形地形上的无线传感器网络非测距节点定位方法，大幅度提高了无线传感器网络节点的定位精度，完全能满足山区地形下无线传感器网络节点定位的要求。
【专利说明】
-种基于马較地形山区无线传感器网络节点的定位方法
技术领域
[0001] 本发明设及一种基于马鞍地形山区无线传感器网络节点的定位方法，属于无线传 感器网络定位技术领域。
【背景技术】
[0002] 我国山区地形(包括丘陵和高原）占69.1 %，森林覆盖率高，具有丰富的矿产资源， 但矿产的开采易造成环境污染。我国地质条件复杂，构造活动频繁，滑坡、泥石流与地裂缝 等灾害隐患多、分布广，且突发性和破坏性强，防范难度大，是世界上地质灾害较为严重、受 威胁人口较多的国家之一。因此，开展环境监测、森林防火、山体滑坡等监测与预警研究具 有十分重要的意义。
[0003] 无线传感器网络(Wireless Sensor化twork，WSN)的出现引起了全世界范围的广 泛关注。1999年，著名的美国商业周刊将无线传感器网络列为21世纪最具影响的21项技术 之一;2003年，Μ口技术评论(Technology Review)在预测未来技术发展的报告中，将其列为 改变世界的10大新技术之一。环境监测与预报、森林防火、地质灾害监测与预警等是无线传 感器网络的重要应用领域。
[0004] 节点定位是无线传感器网络应用的重要基础。在无线传感器网络的应用中，位置 信息对传感器网络的监测活动至关重要。例如：在大面积环境监测中需要知道污染源发生 的地点；森林火灾灾情监测中，需要知道火灾发生的地域;在地质灾害监测预警中，必需要 知道发生险情的时间与地域，W便迅速采取有效行动。
[0005] 目前现有的大多数无线传感器网络定位都是假设应用在理想环境下提出的。例如 二维、Ξ维的定位算法与路由算法适合Ξ维自由空间内的随机节点分布情况，即假定信号 传输模型为理想的球体。无线传感器网络实际应用环境往往是山区，此时，把非平面的网络 用平面的或Ξ维自由空间定位来确定节点位置，往往无法达到理想的性能要求，节点定位 误差较大，难于满足实际应用的要求。
[0006] 根据是否测量距离分为基于测距定位算法和无需测距的定位算法。目前常用的测 距技术有RSSI，T0A，TD0A和A0A。无需测距的定位方法主要有质屯、算法、DV-化P算法、 Amorphous算法、AP 口算法等，其中影响最大、应用最广泛的当属DV-Hop算法。无需测距的定 位机制不需测量节点间的绝对距离或方位，降低了对节点硬件的成本，体积和能量消耗，更 适合于大规模传感器网络。但非测距的节点定位方法的精度较低，在理想环境下其定位精 度仍可满足应用的需求，但在山区复杂地形上的应用时其定位精度难于满足实际应用的要 求。
【发明内容】
[0007] 本发明的目的在于提供一种基于马鞍地形山区无线传感器网络节点的定位方法。 本发明是在山区马鞍形地形上的无线传感器网络非测距节点定位方法，大幅度提高了无线 传感器网络节点的定位精度，完全能满足山区地形下无线传感器网络节点定位的要求。
[0008] 为实现上述目的，本发明采用如下技术方案：
[0009] 本发明中所述的一种基于马鞍地形山区无线传感器网络节点的定位方法，包括如 下步骤：
[0010] (1)在山区马鞍地形上用Ξ维DV-H0P定位算法求出无线传感器网络节点的初始位 置；
[0011] (2)建立山区马鞍地形方程，采用双曲抛物面作为马鞍地形的拟合函数，并结合最 小二乘的曲面拟合算法，在山区电子地形图上取一组离散点的坐标来求拟合参数，求得马 鞍地形方程；
[0012] (3)将步骤(1)所求出节点的初始位置垂直投影到马鞍地形表面上；
[0013] (31)求得过马鞍地形双曲抛物面任一点的法线方程；
[0014] (32)通过在马鞍地形表面上的投影点的法线方程可得2个关于投影点坐标的方 程，另外联立该投影点满足所述马鞍地形方程可求出关于投影点的坐标；
[0015] (33)消去步骤(32)所求得投影点的任意两个坐标，得到关于投影点的坐标的任一 坐标的一元五次方程；
[0016] (34)所述节点初始位置在马鞍地形表面上的投影点有可能有多个，取距离节点初 始位置更近的投影点作为未知节点的修正坐标；
[0017] (4)解关于投影点的坐标的任一坐标的一元五次方程，所述一元五次方程的解法 采用方程的数值解法，步骤如下：
[0018] 步骤一：自动捜索方程实数根的有根区间；
[0019] 步骤二:用二分法求方程的数值解，
[0020] a、如果方程只有唯一的实数根，该一元五次方程的解就是所求的初始位置在马鞍 地形上的投影，即所求未知节点的坐标估计值，最终确定未知节点的坐标；
[0021] b、如果方程有η个实数根，找出η个根中距离所求出的未知节点初始位置最近的那 个根，即为未知节点的最终坐标；
[0022] C、如果方程无实数根，则所求的初始位置就是未知节点的最终坐标。
[0023] 所述定位方法具体包括如下步骤：
[0024] (1)在山区马鞍地形上用Ξ维DV-H0P定位算法求出无线传感器网络节点的初始位 置，将未知节点到信标节点之间的距离用平均每跳距离和两者之间跳段数的乘积表示，使 用四边测量法或极大似然估计法获得未知节点位置信息，具体如下：
[0025] (11)计算未知节点与每个信标节点的最小跳数；
[0026] (12)计算未知节点与信标节点的实际跳段距离；
[0027] 每个信标节点根据记录的其他信标节点的位置信息和相距最小跳数，利用公式S1 计算平均每跳距离hopsize:
[002引
[0029] 其中（Xi，yi，zi)，（Xj，yj，zj)是信标节点i，j的坐标，hj为信标节点i与j(i声j)之间 的跳段数；
[0030] 未知节点i与信标节点j的实际跳段距离为：
[0031] Sdij = hminij Xhopsizei
[0032] (13)未知节点利用上述记录的到各个信标节点的跳段距离，利用四边测量法或极 大似然估计法计算自身坐标；
[0033] (131)所述四边测量法具体如下：
[0034] 已知A、B、C、D四个节点的坐标分别为（Xa，ya，Za)，（xb，yb，zb)，（Xc，yc，Zc)，（xd，yd， Zd)，W及它们到未知节点E的距离分别为da,山，dc，dd;
[0035] 根据两点间距离公式得方程组：
[0036]
[0037] 解方程组：
[003引 X=A-ib
[0039]其中；
[0042] (133)所述极大似然估计法具体如下：
[0043] 若1个未知节点（x，y，z)有n(n〉4)个错节点（（xi，yi，zi)，（X2，y2，Z2)，...（xn，yn, Zn))，可得到矛盾方程组:Αχ = b，
[0046] 则矛盾方程组的解为:x=(ATA厂lATb;
[0047] (2)建立山区马鞍地形方程，采用双曲抛物面方程-=麽作为马鞍地形的拟 合函数，a,b参数可通过实际地形数据拟合马鞍形地形而求出，并结合最小二乘的曲面拟合 算法，在山区电子地形图上取一组离散点的坐标来求拟合参数，求得马鞍地形方程;对双曲 抛物面方程进行等价变换:两边同等乘于a2b2得:b2x2-a 2y2 = a2b2z，式中：a、b为正常数；
[0048] (3)将步骤(1)所求出节点的初始位置垂直投影到马鞍地形表面上；
[0049] (31)求得过马鞍地形双曲抛物面任一点的法线方程；
[0050] 设双曲抛物面上任意一点(X，y，Z)，则 [0051 ] F(x,y ,ζ) =b\^-aV^~a^^z = 〇
[0化2]
[0化3] 法向量为：
[0054] (2b2x，-2a2y，-a2b2)
[0055] 过Po(m，n，k)的法向量方程为：
[0化6]
[0057] (32)通过在马鞍地形表面上的投影点的法线方程可得2个关于投影点坐标的方 程，另外联立该投影点满足所述马鞍地形方程可求出关于投影点的坐标；
[005引（321)由过P0(m，n，k)的法向量方程有
[0062] (322)由过P〇(m，n，k)的法向量方程有：
[00化]整理得：（2z-a2)m-2址+曰2义=0 S3
[0066] (323)投影点P〇(m，n，k)坐标满足双曲抛物面方程
[0067] b2m2-a2n2-a2b2k = 0 S4
[0068] (33)联立方程S1、S2及S3,并消去步骤(32)所求得投影点的任意两个坐标，得到关 于投影点的坐标的任一坐标(运里是关于m)的一元五次方程，并整理得：
[0073] (34)所述节点初始位置在马鞍地形表面上的投影点有可能有多个，取距离节点初 始位置更近的投影点作为未知节点的修正坐标；
[0074] (4)解关于投影点的坐标的任一坐标的一元五次方程，所述一元五次方程的解法 采用方程的数值解法，步骤如下：
[0075] 步骤一：自动捜索方程实数根的有根区间；
[0076] 将公式S5转换为标准方程(m5的系数为1):
[0085]则该方程根的绝对值(即根模)的上下界有如下结论：
[0086] 若a=max{ |ai|，|日2|，…，|an| }，贝防程的根的绝对值小于日+1，对方程S6有：
[0087] aa=max{ I ai I , I 日21 ,…,I an I }
[008引方程S6的根的绝对值小于(aa+1)，则根的捜索区间为：
[0089] [-(aa+1),(aa+1)]
[0090] 设无线传感器网络节点分布范围：[a0，b0]。因为，我们只需要求出范围[a0，b0]内 的根。所W，如果[-(aa+l)，（aa+l)]范围超出[a0，bO]，则从节点的分布范围[a0，b0]捜索方 程的有根区间；
[0091] 步骤二:用二分法求方程的数值解，
[0092] a、如果方程只有唯一的实数根，该一元五次方程的解就是所求的初始位置在马鞍 地形上的投影，即所求未知节点的坐标估计值，最终确定未知节点的坐标；
[0093] al、若方程S5(即f(m)=0)在区间[an，bn]内有唯一一个根，分别将an、bn及 代入方程S5中并计算；
[0094] 曰2、判断，如果
均根；
[0095] 若f (^) · f (孤）巧，则方程的根位于[an, 内，将^赋给bn，则得 到新有根区间为[an,bn];
[0096] 若f .:f化η) <0,则方程的根位于[SSp，地]内，将赋给an,则新 有根区间为[an,bn];
[0097] 若|bn-an|<e(e为精度要求），计算终止，此时否则转曰2;
[0098] b、如果方程有η个实数根，找出η个根中距离所求出的未知节点初始位置最近的那 个根，即为未知节点的最终坐标；
[0099] C、如果方程无实数根，则所求的初始位置就是未知节点的最终坐标。
[0100] 优选地，所述节点的坐标估计值与节点的实际坐标值之间的差距是：
[0101] 其中，（χ/ i，/ i，i)为节点的坐标估计值，（Xi，yi，zi)为节点的实际坐标值。
[0102] 优选地，传感器网络中η个未知节点的平均定位误差Δ表示为：
[0103]
[0104] 其中R为无线通讯距离。
[0105] 本发明的有益效果在于，本发明采样非测距定位机制的节点硬件成本低，体积小 和能量消耗低，更适合于大规模传感器网络。应用传统的非测距无线传感器节点定位技术 在山区地形中的定位误差较大，而采用本发明一山区马鞍地形无线传感器网络节点定位方 法，在山区马鞍地形中应用的定位精度有大幅度提高，能够满足实际应用的要求。本发明无 需测距的定位机制不需测量节点间的绝对距离或方位，降低了对节点硬件的成本，体积和 能量消耗，更适合于大规模传感器网络。但非测距的节点定位方法的精度较低，在理想环境 下其定位精度仍可满足应用的需求，但在山区复杂地形上的应用时其定位精度难于满足实 际应用的要求。本发明是在山区马鞍形地形上的无线传感器网络非测距节点定位方法，大 幅度提高了无线传感器网络节点的定位精度，完全能满足山区地形下无线传感器网络节点 定位的要求。
【附图说明】
[0106] 图1是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法中山 区马鞍形地形示意图；
[0107] 图2是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法中无 线传感器网络节点在山区马鞍形地形随机分布示意图；
[0108] 图3是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法中所 求节点的初始位置垂直投影到马鞍地形表面示意图，其中P(x，y，z)为所求节点初始位置， Po(m,n,k)与Pi(ml ,η? ,kl)是P在马鞍地形表面上的投影点；
[0109] 图4是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法中所 述四边测量法示意图；图中四个错节点A、B、C、D的坐标分别为(Xa，ya，Za)，（Xb，yb，Zb)，（Xc， yc，Zc)，（xd，yd，zd)，W及它们到未知节点E的距离分别为da，db，dc，dd;
[0110] 图5是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法仿真 实验中不同错节点比例的平均相对误差曲线；
[0111] 图6是本发明所述的一种基于马鞍地形山区无线传感器网络节点的定位方法仿真 实验中不同通信半径的平均相对误差曲线。
【具体实施方式】
[0112] 结合W下具体实例对本发明的原理和特征进行描述，所举实例只用于解释本发 明，并非用于限定本发明的范围。
[0113] -种基于马鞍地形山区无线传感器网络节点的定位方法，所述定位方法具体包括 如下步骤：
[0114] (1)在山区马鞍地形上用Ξ维DV-HOP定位算法求出无线传感器网络节点的初始位 置，参阅图2,将未知节点到信标节点之间的距离用平均每跳距离和两者之间跳段数的乘积 表示，使用四边测量法或极大似然估计法获得未知节点位置信息，具体如下：
[0115] (11)计算未知节点与每个信标节点的最小跳数；
[0116] (12)计算未知节点与信标节点的实际跳段距离；
[0117] 每个信标节点根据记录的其他信标节点的位置信息和相距最小跳数，利用公式S1 计算平均每跳距离hopsize:
[011 引
[0119]其中（Xi，yi，zi)，（Xj，yj，zj)是信标节点i，j的坐标，hj为信标节点i与j(i声j)之间 的跳段数；
[0120] 未知节点i与信标节点j的实际跳段距离为：
[0121] Sdij = hminij Xhopsizei
[0122] (13)未知节点利用上述记录的到各个信标节点的跳段距离，利用四边测量法或极 大似然估计法计算自身坐标；
[0123] (131)参阅图4，所述四边测量法具体如下：
[0124] 已知A、B、C、D四个节点的坐标分别为（Xa，ya，Za)，（xb，yb，zb)，（Xc，yc，Zc)，（xd，yd， Zd)，W及它们到未知节点E的距离分别为da,山，dc，dd;
[0125] 根据两点间距离公式得方程组：
[0126]
[0127]解方程组：
[012 引 X=A-ib
[0129]其中；
[0132] (134)所述极大似然估计法具体如下：
[0133] 若1个未知节点（X ,y, Z)有n(n〉4)个错节点（（XI ,yi ,ζι), (X2，y2，Z2)，…（xn，yn, Zn))，可得到矛盾方程组:Αχ = b，
[0136] 则矛盾方程组的解为:x= (ATA厂Vb;
[0137] (2)建立山区马鞍地形方程，因双曲抛物面能很好的表示山区马鞍地形，如图1所 示，因此，采用双曲抛物面方程
作为马鞍地形的拟合函数，a, b参数可通过实际 地形数据拟合马鞍形地形而求出，并结合最小二乘的曲面拟合算法，在山区电子地形图上 取一组离散点的坐标来求拟合参数，求得马鞍地形方程;对双曲抛物面方程进行等价变换： 两边同等乘于a2b2得:b2x2-aY = a2b2z，式中：a、b为正常数；
[0138] (3)将步骤(1)所求出节点的初始位置垂直投影到马鞍地形表面上，参阅图3;
[0139] (31)求得过马鞍地形双曲抛物面任一点的法线方程；
[0140] 设双曲抛物面上任意一点(X，y，Z)，则
[0143] 法向量为：
[0144] (2b\,-2aV.-aV)
[0145] 过P〇(m，n，k)的法向量方程为：
[0146]
[0147] (32)通过在马鞍地形表面上的投影点的法线方程可得2个关于投影点坐标的方 程，另外联立该投影点满足所述马鞍地形方程可求出关于投影点的坐标；
[014引（321)由过P0(m，n，k)的法向量方程有
[0152] (322)由过Po(m，n，k)的法向量方程有：
[0156] (323)投影点Po(m，n，k)坐标满足双曲抛物面方程
[0157] b2m2-a2n2-a2b2k=0 S4
[0158] (33)联立方程S1、S2及S3,并消去步骤(32)所求得投影点的任意两个坐标，得到关 于投影点的坐标的任一坐标的一元五次方程，并整理得：
[0162] (34)所述节点初始位置在马鞍地形表面上的投影点有可能有多个，取距离节点初 始位置更近的投影点作为未知节点的最终坐标；
[0163] (4)解关于投影点的坐标的任一坐标的一元五次方程，所述一元五次方程的解法 采用方程的数值解法，步骤如下：
[0164] 步骤一：自动捜索方程实数根的有根区间；
[0165] 将公式S5转换为标准方程(m5的系数为1):
[0174] 则该方程根的绝对值(即根模)的上下界有如下结论：
[0175] 若a=max{ I ai I，I日21，…，I an I }，则方程的根的绝对值小于a+1，对方程S6有：
[0176] aa=max{ |ai| , |日2| ,…,|an| }
[0177] 方程S6的根的绝对值小于(aa+1)，则根的捜索区间为：
[017引[-(aa+1),(aa+1)]
[0179] 设无线传感器网络节点分布范围：[a0，b0]。因为，我们只需要求出范围[a0，b0]内 的根。所W，如果[-(aa+1)，（aa+1)]范围超出[a0，bO]，则从节点的分布范围[a0，b0]捜索方 程的有根区间；
[0180] 步骤二:用二分法求方程的数值解，
[0181] a、如果方程只有唯一的实数根，该一元五次方程的解就是所求的初始位置在马鞍 地形上的投影，即所求未知节点的坐标估计值，最终确定未知节点的坐标；
[0182] 曰1、若方程55(即^111)=0)在区间[日11，扣]内有唯一一个根，分别将 代入方程S5中并计算；
[0183] 曰2、判断，如果
拍(m)=0的根；
[0184] 老
武给bn,则得 到新有根区间为[an,bn];
[0185] 老
武给an,则新 有根区间为[an,bn];
[0186] 若|bn-an|<e(e为精度要求），计算终止，此时，否则转a2;
[0187] b、如果方程有η个实数根，找出η个根中距离所求出的未知节点初始位置最近的那 个根，即为未知节点的最终坐标；
[0188] C、如果方程无实数根，则所求的初始位置就是未知节点的最终坐标。
[0189] 值得注意的是，所述节点的坐标估计值与节点的实际坐标值之间的差距是：
[0190] 其中，（χ/l，/l，z/l)为节点的坐标估计值，（Xl，yl，Zl)为节点的实际坐标值。传感 器网络中η个未知节点的平均定位误差Δ表示为：
[0191]
[0192] 其中R为无线通讯距离。
[0193] 仿真实验：
[0194] 为了检验算法的性能，在山区马鞍地形环境中，对经典Ξ维DV-Hop定位算法、改进 的Ξ维DV-Hop定位算法和马鞍地形节点定位算法在Matlab平台上进行了仿真对比分析。研 究了错节点不同比例与定位误差的关系W及节点不同通信半径R与定位误差的关系。每次 仿真都是在马鞍地形中均匀随机布放200个节点，马鞍地形在水平面的投影范围是lOOmX 100m。图5是一次仿真实验的无线传感器网络节点均匀随机分布图。
[01M] (1)不同错节点比例的节点定位算法仿真
[0196] 通信半径R = 30m，分别对错节点占15%，25%，35%，45%和50%时进行了模拟实 验。对不同错节点比例的仿真都次都进行了 50次，然后取相对均方误差的平均值，每次实验 中未知节点和错节点的坐标都是随机产生的，实验结果如图5所示。从图可知，马鞍地形节 点定位算法精度明显高于传统定位算法的精度，错节点比例越大，定位精度越高。经典Ξ维 DV-H0P定位算法与改进的Ξ维DV-H0P定位算法在山区定位误差太大，难于满足实际定位需 要，而马鞍地形节点定位算法精度高，完全能满足实际应用的需求。
[0197] (2)不同通信半径的节点定位算法仿真
[019引无线传感器节点总数200,错节点50,错节点比例25%，分别对通信半径R = 30m， 35m，40m，45m，50m，55m，60m时进行了模拟实验。对每种通信半径仿真都进行50次，然后取相 对均方误差的平均值，每次实验中未知节点和错节点的坐标都是随机产生的，实验结果如 图6所示。从图可知，马鞍地形节点定位算法精度明显高于传统定位算法的精度，定位精度 与通信半径关系不大。而马鞍地形节点定位算法精度高，完全能满足实际应用的需求。
[0199] 基于上述，本发明采样非测距定位机制的节点硬件成本低，体积小和能量消耗低， 更适合于大规模传感器网络。应用传统的非测距无线传感器节点定位技术在山区地形中的 定位误差较大，而采用本发明一山区马鞍地形无线传感器网络节点定位方法，在山区马鞍 地形中应用的定位精度有大幅度提高，能够满足实际应用的要求。
[0200] 尽管本发明是参照具体实施例来描述，但运种描述并不意味着对本发明构成限 审IJ。参照本发明的描述，所公开的实施例的其他变化，对于本领域技术人员都是可W预料 的，运种的变化应属于所属权利要求所限定的范围内。
【主权项】
1. 一种基于马鞍地形山区无线传感器网络节点的定位方法，其特征在于，包括如下步 骤： (1) 在山区马鞍地形上用三维DV-HOP定位算法求出无线传感器网络节点的初始位置； (2) 建立山区马鞍地形方程，采用双曲抛物面作为马鞍地形的拟合函数，并结合最小二 乘的曲面拟合算法，在山区电子地形图上取一组离散点的坐标来求拟合参数，求得马鞍地 形方程； (3) 将步骤（1)所求出节点的初始位置垂直投影到马鞍地形表面上； (31) 求得过马鞍地形双曲抛物面任一点的法线方程； (32) 通过在马鞍地形表面上的投影点的法线方程可得2个关于投影点坐标的方程，另 外联立该投影点满足所述马鞍地形方程可求出关于投影点的坐标； (33) 消去步骤(32)所求得投影点的任意两个坐标，得到关于投影点的坐标的任一坐标 的一元五次方程； (4) 解关于投影点的坐标的任一坐标的一元五次方程，所述一元五次方程的解法采用 方程的数值解法，其步骤如下： 步骤一：自动搜索方程实数根的有根区间； 步骤二:用二分法求方程的数值解， 方程的数值解有如下三种情况： a、 如果方程只有唯一的实数根，该一元五次方程的解就是所求的初始位置在马鞍地形 上的投影，即所求未知节点的坐标估计值，最终确定未知节点的坐标； b、 如果方程有η个实数根，找出η个根中距离所求出的未知节点初始位置最近的那个 根，即为未知节点的最终坐标； c、 如果方程无实数根，则所求的初始位置就是未知节点的最终坐标。2. 如权利要求1所述的一种基于马鞍地形山区无线传感器网络节点的定位方法，其特 征在于，所述定位方法具体包括如下步骤： (I) 在山区马鞍地形上用三维DV-HOP定位算法求出无线传感器网络节点的初始位置， 将未知节点到信标节点之间的距离用平均每跳距离和两者之间跳段数的乘积表示，使用四 边测量法或极大似然估计法获得未知节点位置信息，具体如下： (II) 计算未知节点与每个信标节点的最小跳数hmin; (12) 计算未知节点与信标节点的实际跳段距离Sd; 每个信标节点根据记录的其他信标节点的位置信息和相距最小跳数，利用公式S1计算 平均每跳距离hopsize:其中（Xi，yi，zi)，（Xj，yj，zj)是信标节点i，j的坐标，hj为信标节点i与j(i辛j)之间的跳 段数； 未知节点i与信标节点j的实际跳段距离为： Sdi j - hmini j X h〇pS iΖΘ? (13) 未知节点利用上述记录的到各个信标节点的跳段距离，利用四边测量法或极大似 然估计法计算自身坐标； (131) 所述四边测量法具体如下： 已知A、B、C、D四个节点的坐标分别为(xa，ya，za)，（Xb，yb，zb)，（Xc，y c，zc)，（Xd，yd，zd)，以 及它们到未知节点E的距离分别为c^dhde^dd; 根据两点间距离公式得方程组：解方程组： X=A_1b 其中：(132) 所述极大似然估计法具体如下： 若1 个未知节点（x，y，z)有n(n>4)个锚节点（（xi，yi，zi)，（X2，y2，Z2)，'"（xn，y n，zn))，可 得到矛盾方程组:Ax = b，则矛盾方程组的解为:X = (αΤαγΙ% ; (2)建立山区马鞍地形方程，采用双曲抛物面方程作为马鞍地形的拟合函 数，a，b参数可通过实际地形数据拟合马鞍形地形而求出，并结合最小二乘的曲面拟合算 法，在山区电子地形图上取一组离散点的坐标来求拟合参数，求得马鞍地形方程;对双曲抛 物面方程进行等价变换:两边同等乘于a 2b2得:b2x2-a2y 2 = a2b2z，式中:a、b为正常数； (3) 将步骤（1)所求出节点的初始位置垂直投影到马鞍地形表面上； (31) 求得过马鞍地形双曲抛物面任一点的法线方程； 设双曲抛物面上任意一点(x，y，z)，则 F(x,yz) =b2x2-a2y2-a2b 2z = 0法向量为： (2b2x，_2a2y，_a2b 2) 过Ρο (m，η，k)的法向量方程为：(32) 通过在马鞍地形表面上的投影点的法线方程可得2个关于投影点坐标的方程，另 外联立该投影点满足所述马鞍地形方程可求出关于投影点的坐标； (321) 由过P〇 (m，η，k)的法向量方程有a2n(x_m) =-b2m(y_n) 整理得：b2ym-(a2+b2)mn+a2xn = 0 S2 (322) 由过P〇 (m，η，k)的法向量方程有：a2(x_m) =-2(z_k)m 整理得：（2z-a2)m-2mk+a2x0 S3 (323) 投影点P〇(m，n，k)坐标满足双曲抛物面方程 b2m2-a2n2-a2b 2k = 0 S4 (33) 联立方程S1、S2及S3,并消去步骤(32)所求得投影点的任意两个坐标，得到关于投 影点的坐标的任一坐标的一元五次方程，并整理得： 2 (a2+b2) V-4a2 (a2+b2) xm4+a2 [ (a2-2z) (a2+b2)2+2a2x2 -2b2y2 ]m3_a4(a2+b2) (3a2+b2_4z) xm2+a6 (3a2+2b2 -2z)x2m_a8x3 = 0 S5 (34) 所述节点初始位置在马鞍地形表面上的投影点有可能有多个，取距离节点初始位 置更近的投影点作为未知节点的修正坐标，； (4) 解关于投影点的坐标的任一坐标的一元五次方程，所述一元五次方程的解法采用 方程的数值解法，其步骤如下： 步骤一：自动搜索方程实数根的有根区间； 将公式S5转换为标准方程(m5的系数为1): m5+a2m4+b2m3+C2m2+d2m+e2 = 0 S6若有方程： f (X) = xn+aixn_1+a2Xn-2+· · · +an-ix+an=Ο 则该方程根的绝对值(即根模)的上下界有如下结论： 若a=max{ |ai|，|a〗|，···，|an| }，则方程的根的绝对值小于a+Ι，对方程S6有： aa=max{ |ai|，|a〗|，···，|an| } 方程S6的根的绝对值小于(aa+1)，则根的搜索区间为： [_(aa+l)，（aa+Ι)] 设无线传感器网络节点分布范围：[aO，bO]，因为，我们只需要求出范围[aO，bO]内的 根，所以，如果[_(aa+l)，（aa+l)]范围超出[aO，bO]，则从节点的分布范围[aO，bO]搜索方程 的有根区间； 步骤二:用二分法求方程的数值解， a、 如果方程只有唯一的实数根，该一元五次方程的解就是所求的初始位置在马鞍地形 上的投影，即所求未知节点的坐标估计值，最终确定未知节点的坐标； al、若方程S5在区间[an，bn]内有唯 个根，分别将an、bn及代入方程S5中并计 算；，则方程的根位于赋给bn，则得到新有根区间为[an，bn];，则方程的根位于 赋给an，则新有根 区间为[an，bn]; 若| bn-an | <ε (ε为精度要求），计算终止：b、 如果方程有η个实数根，找出η个根中距离所求出的未知节点初始位置最近的那个 根，即为未知节点的最终坐标； c、 如果方程无实数根，则所求的初始位置就是未知节点的最终坐标。3.如权利要求1-2任一项所述的一种基于马鞍地形山区无线传感器网络节点的定 位方法，其特征在于，所述节点的坐标估计值与节点的实际坐标值之间的差距是：其中，（X7 iY iY i)为节点的坐标估计值，^#，2：1)为节点的实际坐标值。4.如权利要求3所述的一种基于马鞍地形山区无线传感器网络节点的定位方法，其特 征在于，传感器网络中η个未知节点的平均定位误差△表示为：其中R为无线通讯距离。
【文档编号】H04W64/00GK106060927SQ201610668135
【公开日】2016年10月26日
【申请日】2016年8月15日 公开号201610668135.4, CN 106060927 A, CN 106060927A, CN 201610668135, CN-A-106060927, CN106060927 A, CN106060927A, CN201610668135, CN201610668135.4
【发明人】胡中栋, 曾珽, 王振东
【申请人】江西理工大学