[发明专利]一种基于旋转RSS的单信标定位方法

权利要求书

1.一种基于旋转RSS的单信标定位方法，其特征在于，该方法具体按照以下步骤进行：

步骤一，场景设置：

步骤S11：通过一组发送节点和一个位置已知的接收节点搭建一个没有任何遮挡的传播环境，位置已知的接收节点作为信标，发送节点向接收节点发送无线信号，发送节点和接收节点之间形成视距链路，选择一个体格信息已知的人作为阻挡视距链路的载体；

步骤二，方向确定：

步骤S21：人围绕接收节点以相同的半径做匀速旋转，接收节点检测并记录各个发送节点的旋转RSS信息，并将该信息上传给服务器进行存储；

步骤S22：根据人的体格信息、人绕接收节点的旋转半径、无线发送节点和无线接收节点间的距离信息，分析并计算不同体格的人对无线信号的遮挡角度θ_B，具体计算过程如下所示：

b2′=d12+b2·|sin(θ-arctanbd1)|---(2)]]>

d1′′=d12+b2·|cos(θ-arctanbd1)|---(4)]]>

r1′′=λd1′′(d-d1′′)d---(9)]]>

B=(r1′′)2×arccosb2′r1′′-b2′×(r1′′)2-(b2′)2---(10)]]>

θ_m是满足下式的最小角度，即满足下式的最小θ；

Bπ(r1′′)2≤45%---(11)]]>

θ_B=2θ_m

式中：

2b表示人的宽度；

d₁表示未旋转时人与接收节点的距离，也是人围绕接收节点的旋转半径；

λ表示无线信号的波长；

d表示发送节点与接收节点的视距链路长度；

θ表示人的旋转角度；

2α表示未旋转时由人引起的遮挡角度；

b₂＇表示当人转动角度θ满足α＜θ＜π/2时，人左侧与视距链路的垂直距离；

d₁〞表示当人转动角度θ满足α＜θ＜π/2时，人左侧距接收节点的水平距离；

r₁〞表示当人转动角度θ满足α＜θ＜π/2时，在d₁〞处对应的第一菲涅尔区半径；

B表示当人转动角度θ满足α＜θ＜π/2时，遮挡第一菲涅尔区的面积；

步骤S23：根据步骤S22中计算出的遮挡角度θ_B选择移动窗口w，发送节点和接收节点的视距链路长度为d，d=d₁+d₂，d₁确定后，首先找到d的变化范围内大于并且最接近最大遮挡角度θ_B的偶数角度作为窗口角度，然后取这个窗口角度加1度得到的奇数作为移动窗口w，在0到360度的旋转范围内，计算移动窗口w对应角度范围内旋转RSS绝对值的累加和，将累加和最大时移动窗口的中心位置对应的角度作为第i个发送节点的方向θ_i，从而得到所有发送节点的方向集合Θ={θ_i}；

θi=argmaxθip(Σk=-w-12w-12|RSS(Di,θi(p+k))|),P∈(1,...,K)---(12)]]>

式中：

θ_i表示第i个发送节点相对于接收节点的方向；

θ_ip表示第i个发送节点旋转RSS的第p个角度；

D_i表示第i个发送节点相对于接收节点的距离；

RSS(D_i,θ_i(p+k))表示当人的旋转角度为θ_i(p+k)时，接收节点收到第i个发送节点的实际信号强度；

K表示旋转时记录的角度数目；

步骤三，距离确定：

步骤S31：根据无线通信中的绕射理论，建立简单的人遮挡无线信号的绕射模型，具体的绕射模型如下：

EdE0=F(v)=(1+j)2∫v∞exp((-jπt2)/2)dt---(13)]]>

v=b2(d1+d2)λd1d2---(14)]]>

G_d＝20log|F(v)|       (15)

式中：

E_d表示接收节点接收到的电场强度；

E₀表示自由空间的电场强度；

F(v)表示电场强度函数；

v表示Fresnel-Kirchoff绕射参数；

t表示积分变量；

d₂表示未旋转时人与发送节点的距离；

步骤S32：根据步骤S23中确定出的发送节点方向Θ={θ_i}和步骤S21中收集到的旋转RSS信息，计算得到无人遮挡时所有发送节点的实际信号强度集合R(d)={R(D_i)}和有人遮挡引起的实际绕射增益集合G_d={G_d(D_i,θ_i)}，具体计算公式如下所述：

R(Di)=1qΣk=1qRSS(Di,θk)---(20)]]>

Gd(Di,θi)=RSS(Di,θi)-1qΣk=1qRSS(Di,θk)---(21)]]>

式中：

R(D_i)表示无人遮挡时接收节点收到第i个发送节点的实际信号强度；

RSS(D_i,θ_k)表示当人的旋转角度为θ_k时，接收节点收到第i个发送节点的实际信号强度；

θ_k表示在角度区间(θ_i+π-π/4,θ_i+π+π/4)内的旋转角度，在此角度区间内人不遮挡发送节点和接收节点的通信链路；

q表示在角度区间(θ_i+π-π/4,θ_i+π+π/4)内，θ_k的总数目；

G_d(D_i,θ_i)表示人站在第i个发送节点与接收节点的视距链路上时引起的绕射增益；

RSS(D_i,θ_i)表示人站在第i个发送节点与接收节点的视距链路上时，接收节点收到第i个发送节点的实际信号强度；

步骤S33：在LDPL模型的两个参数P和n的取值范围内，根据每一组P_s、n_t和步骤S32中得到的实际信号强度集合R(d)={R(D_i)}计算出一组距离集合d_st={D_i}，然后利用此距离集合和步骤S31中建立的绕射模型计算出一组理论绕射增益集合G_dst={G_d(D_i)}，利用公式（17）得到距离集合，利用公式（15）得到理论绕射增益集合；

R(d)＝P-10nlogd+N    (17)

G_d=20log|F(v)|    (15)

步骤S34：将步骤S32中得到的实际绕射增益集合G_d={G_d(D_i,θ_i)}和步骤S33中得到的每一组理论绕射增益集合G_dst={G_d(D_i)}作差，得到相应的差异集合

DIF_st={|G_d(D_i,θ_i)-G_d(D_i)|}，将最小的一组差异集合DIF_mn对应的P_m和n_n作为LDPL模型的参数，并将此P_m和n_n对应的距离集合d_mn={D_i}作为所有发送节点相对于接收节点距离的最终结果；

步骤四，位置确定：

步骤S41：将步骤S23中所有发送节点相对于接收节点的方向结果和步骤S34中所有发送节点相对于接收节点的距离结果合并，得到每个发送节点相对于接收节点位置的极坐标表示，再根据接收节点的真实位置坐标，将发送节点位置相对于接收节点位置的极坐标转换为发送节点位置的笛卡尔坐标，即每个发送节点的位置坐标；

xi=x0+Di×cosθiyi=y0+Di×sinθi]]>

式中：

(x₀,y₀)是接收节点的位置坐标；

(x_i,y_i)是第i个发送节点的位置坐标；

D_i是第i个发送节点相对于接收节点的距离；

θi是第i个发送节点相对于接收节点的方向。