发明名称 一种无线传感器节点随机部署的被动式目标定位方法
摘要 本发明公开了一种无线传感器节点随机部署被被动式目标定位方法,具体包括如下步骤:步骤1:在监测区域内部署传感器节点和基站;步骤2:监测前无干扰情况下的RSS值采集;步骤3:监测阶段的RSS值采集及分类;步骤4:由于监测区域内节点是随机部署的,要想准确定位到目标位置,首先需要确定监测区域内被目标所干扰的链路的位置,即构成这些被干扰链路的节点的位置;步骤5:监测区域目标位置的确定。本发明的方法无需先验知识即可进行,另外,该方法无需网格划分,从而规避了基于网格划分的定位方法中网格大小影响定位精度的缺点。
申请公布号 CN103634907B 申请公布日期 2017.01.25
申请号 CN201310651614.1 申请日期 2013.12.04
申请人 西北大学 发明人 房鼎益;王举;陈晓江;常俪琼;邢天璋;聂卫科;刘晨;肖云
分类号 H04W64/00(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W64/00(2009.01)I
代理机构 西安恒泰知识产权代理事务所 61216 代理人 李婷
主权项 一种无线传感器节点随机部署被被动式目标定位方法,其特征在于,具体包括如下步骤:步骤1:在监测区域内部署传感器节点和基站:1)节点部署:在监测区域内随机分散部署N个传感器节点,保证节点两两之间能够通信;2)基站部署:在监测区域之外部署一个基站;步骤2:监测前无干扰情况下的RSS值采集,具体操作如下:定义l<sub>ij</sub>表示由发送节点i和接收节点j构成的链路,F<sub>ij</sub>表示链路l<sub>ij</sub>在未受干扰时的RSS值,该值由接收节点j发送给基站,监测区域中共有N(N‑1)条链路;目标未进入监测区域时,监测区域内的N(N‑1)条链路都未受干扰;基站采集每条链路l<sub>ij</sub>的RSS值F<sub>ij</sub>,其中1<i<N,1<j<N;基站将采集好所有链路未受干扰的RSS值发送给PC机,得到N(N‑1)个链路在未受干扰时的RSS值;步骤3:监测阶段的RSS值采集及分类:开始监测后,基站每隔一定时间收集监测区域内N(N‑1)条链路中每条链路l<sub>ij</sub>的RSS值y<sub>ij</sub>,即收集到N(N‑1)个链路的RSS值,然后将这些RSS值发送给PC机;由于目标的进入会对监测区域内多条链路的RSS值产生干扰;我们知道,目标在LOS造成的绕射会使链路l<sub>ij</sub>的RSS值比没有目标时的链路的RSS值F<sub>ij</sub>低;而在NLOS造成的散射会使链路l<sub>ij</sub>的RSS值比没有目标时的链路的RSS值F<sub>ij</sub>高;PC机将收集到的所有链路的RSS值分类为绕射干扰RSS集合Y<sub>d</sub>、散射干扰RSS集合Y<sub>s</sub>和未受干扰RSS集合Y<sub>f</sub>;若y<sub>ij</sub>属于集合Y<sub>d</sub>或集合Y<sub>s</sub>,则说明链路l<sub>ij</sub>受到了目标的干扰,认为监测区域内有目标进入;若y<sub>ij</sub>属于集合Y<sub>f</sub>,则说明链路l<sub>ij</sub>没有受到目标的干扰,认为监测区域内没有目标进入;步骤4:由于监测区域内节点是随机部署的,要想准确定位到目标位置,首先需要确定监测区域内被目标所干扰的链路的位置,即构成这些被干扰链路的节点的位置;步骤5:监测区域目标位置的确定;所述步骤4中被干扰链路的节点位置的获取包括如下步骤:第一步:通过人工测量或其他方法得到监测区域内任意5个节点的位置C<sub>i</sub>、路径损耗因子γ<sub>i</sub>和节点发射功率P<sub>i</sub>,其中C<sub>i</sub>=[x<sub>i</sub>,y<sub>i</sub>]是二维向量,x<sub>i</sub>,y<sub>i</sub>代表节点i的二维坐标;对于监测区域内其余未知的N‑5个节点,分别随机给定每个节点的位置C<sub>i</sub>、路径损耗因子γ<sub>i</sub>和节点发射功率P<sub>i</sub>做为初始值;C<sub>i</sub>的各维度的取值在0~max(a,b)范围内;γ<sub>i</sub>取2~6;节点的发射功率P<sub>i</sub>为0~‑20dBm;a、b分别为监测区域的长、宽;第二步:将位置C<sub>i</sub>,路径损耗因子γ<sub>i</sub>和节点发射功率P<sub>i</sub>均已知的节点的集合定义为集合C<sub>done</sub>;定义一个集合O,并令其初始为空集;第三步:定位组成被干扰链路的节点的位置:针对集合Y<sub>d</sub>和集合Y<sub>s</sub>中每一个链路,找到所有位置C<sub>i</sub>未知的节点;第四步:针对位置C<sub>i</sub>未知的节点i,即<img file="FDA0001066804890000021.GIF" wi="182" he="55" />遍历监测区域内所有与节点i相连的链路l<sub>ij</sub>,如果j∈C<sub>done</sub>,则令O=O∪{F<sub>ij</sub>};当集合O的元素个数不少于5时,求解式(2)的最优解,即求得一组能满足式(2)最小的解{P<sub>i</sub>,C<sub>i</sub>,γ<sub>i</sub>};<maths num="0001"><math><![CDATA[<mrow><mi>m</mi><mi>i</mi><mi>n</mi><mfrac><mn>1</mn><mrow><mo>|</mo><mi>O</mi><mo>|</mo></mrow></mfrac><mi>&Sigma;</mi><mo>|</mo><msub><mi>F</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mi>i</mi></msub><mo>+</mo><mn>10</mn><msub><mi>&gamma;</mi><mi>i</mi></msub><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000022.GIF" wi="982" he="151" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000023.GIF" wi="1125" he="93" /></maths>其中,|O|表示集合的大小;F<sub>ij</sub>∈O;C<sub>i</sub>=[x<sub>i</sub>,y<sub>i</sub>],x<sub>i</sub>,y<sub>i</sub>代表节点i的二维坐标;C<sub>j</sub>=[x<sub>j</sub>,y<sub>j</sub>],x<sub>j</sub>,y<sub>j</sub>代表节点j的二维坐标;第五步:重复执行第三、四步,定位出所有组成被干扰链路的节点的位置;所述步骤5监测区域目标位置的确定具体包括如下步骤:第一步:随机给定目标的位置X<sub>t</sub>,X<sub>t</sub>=[x<sub>t</sub>,y<sub>t</sub>]是二维向量,x<sub>t</sub>,y<sub>t</sub>代表目标的二维坐标;一般情况下,监测区域近似为a×b的矩形,则X<sub>t</sub>的各维度取值在0~max(a,b)范围内;第二步:求满足式(4)最小的解X<sub>t</sub>,并将X<sub>t</sub>作为目标t的真实位置;<maths num="0003"><math><![CDATA[<mrow><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><mfrac><mn>1</mn><mrow><mo>|</mo><msub><mi>Y</mi><mi>d</mi></msub><mo>|</mo></mrow></mfrac><mi>&Sigma;</mi><mo>|</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>P</mi><mi>i</mi></msub><mo>+</mo><mn>10</mn><msub><mi>&gamma;</mi><mi>i</mi></msub><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>D</mi><mrow><mi>i</mi><mi>j</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>|</mo><mo>+</mo></mrow>]]></math><img file="FDA0001066804890000024.GIF" wi="1509" he="190" /></maths><maths num="0004"><math><![CDATA[<mrow><mfrac><mn>1</mn><mrow><mo>|</mo><msub><mi>Y</mi><mi>s</mi></msub><mo>|</mo></mrow></mfrac><mi>&Sigma;</mi><mo>|</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>P</mi><mi>i</mi></msub><mo>+</mo><mn>10</mn><msub><mi>&gamma;</mi><mi>i</mi></msub><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>i</mi><mi>j</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>|</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000031.GIF" wi="1420" he="182" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>i</mi><mi>j</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>20</mn><mi>l</mi><mi>o</mi><mi>g</mi><mo>&lsqb;</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><mo>&CenterDot;</mo><mo>|</mo><msubsup><mo>&Integral;</mo><mi>v</mi><mi>&infin;</mi></msubsup><mi>exp</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msup><mi>j&pi;t</mi><mn>2</mn></msup></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>|</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000032.GIF" wi="885" he="164" /></maths>式(5)中‑jπt<sup>2</sup>中的j是虚数单位,t是积分变量;<maths num="0006"><math><![CDATA[<mrow><mi>v</mi><mo>=</mo><mi>h</mi><msqrt><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mi>&lambda;</mi><mo>&CenterDot;</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow></msub></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000033.GIF" wi="602" he="165" /></maths>式(6)中h的值等于目标高度减去节点距离地面的高度H;<maths num="0007"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>i</mi><mi>j</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>10</mn><mi>l</mi><mi>o</mi><mi>g</mi><mfrac><mrow><msub><mi>P</mi><mi>i</mi></msub><msub><mi>G</mi><mi>i</mi></msub><msub><mi>G</mi><mi>j</mi></msub><msup><mi>&lambda;</mi><mn>2</mn></msup><mi>&sigma;</mi></mrow><mrow><msup><mrow><mo>(</mo><mn>4</mn><mi>&pi;</mi><mo>)</mo></mrow><mn>3</mn></msup><mo>&CenterDot;</mo><msubsup><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow><mn>2</mn></msubsup><mo>&CenterDot;</mo><msubsup><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000034.GIF" wi="612" he="150" /></maths><maths num="0008"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000035.GIF" wi="1142" he="94" /></maths><maths num="0009"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>t</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>t</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000036.GIF" wi="1117" he="87" /></maths><maths num="0010"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>j</mi><mi>t</mi></mrow></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>j</mi></msub><mo>-</mo><msub><mi>C</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>C</mi><mi>j</mi></msub><mo>-</mo><msub><mi>C</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>t</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>-</mo><msub><mi>y</mi><mi>t</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001066804890000037.GIF" wi="1174" he="92" /></maths>其中,y′<sub>ij</sub>是集合Y<sub>d</sub>内的y<sub>ij</sub>的取值;y″<sub>ij</sub>是集合Y<sub>s</sub>内的y<sub>ij</sub>的取值;γ<sub>i</sub>为链路l<sub>ij</sub>的路径损耗因子;P<sub>i</sub>为节点i的发射功率;d<sub>ij</sub>为节点i和节点j之间的距离;Y<sub>d</sub>为绕射干扰的RSS值集合,Y<sub>s</sub>为散射干扰的RSS值集合;|Y<sub>d</sub>|和|Y<sub>s</sub>|分别是集合Y<sub>d</sub>和Y<sub>s</sub>的大小;λ是节点发出的无线信号波长,λ=0.125米;G<sub>i</sub>,G<sub>j</sub>分别为节点i和节点j的天线增益,G<sub>i</sub>=G<sub>j</sub>=1;σ是目标的RCS,取1;d<sub>it</sub>是节点i和目标t之间的距离,d<sub>jt</sub>是节点j和目标t之间的距离;C<sub>t</sub>=X<sub>t</sub>。
地址 710069 陕西省西安市太白北路229号