主权项 |
一种基于测距信息融合的室内定位方法,其特征在于所述方法包括如下步骤:在WLAN网络中加入精准测距基站BS,在精准测距基站BS位置(x<sub>1</sub>,y<sub>1</sub>)已知情况下,以精准测距基站BS为圆心,以测得的与移动站MS的距离d<sub>1</sub>为半径画圆,将所估计的移动站MS坐标定位于圆上,圆上各点坐标(x<sub>Mi</sub>,y<sub>Mi</sub>)表示为:x<sub>Mi</sub>=x<sub>1</sub>+d<sub>1</sub>*cos(i*θ) (1)y<sub>Mi</sub>=y<sub>1</sub>+d<sub>1</sub>*sin(i*θ)其中i=1,2,...,N,θ是角位移的单位转动量,θ=2π/N;同时,WLAN网络中的各个接入点AP还会扫描到移动站MS的发射信号强度,产生对应移动站MS的接收信号强度RSS值,通过对数正态模型,在接收信号强度RSS估计中测量值P<sub>ij</sub>是节点i接收到由节点j发出的信号强度,其对数形式P<sub>ij</sub>(dBm)=10log<sub>10</sub>P<sub>ij</sub>服从高斯正态分布<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>dBm</mi><mo>)</mo></mrow><mo>~</mo><mi>N</mi><mrow><mo>(</mo><msub><mover><mi>P</mi><mo>‾</mo></mover><mi>ij</mi></msub><mrow><mo>(</mo><mi>dBm</mi><mo>)</mo></mrow><mo>,</mo><msubsup><mi>σ</mi><mi>dB</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000011.GIF" wi="1623" he="114" /></maths>其中均值<img file="FDA0000715618420000012.GIF" wi="60" he="87" />可以表示为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>P</mi><mo>‾</mo></mover><mi>ij</mi></msub><mrow><mo>(</mo><mi>dBm</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>dBm</mi><mo>)</mo></mrow><mo>-</mo><mn>10</mn><mi>β</mi><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><msub><mi>d</mi><mi>ij</mi></msub><mo>/</mo><msub><mi>d</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000013.GIF" wi="823" he="146" /></maths>其中:P<sub>0</sub>(dBm)是参考距离d<sub>0</sub>的接收信号强度分贝,根据自由空间路径损耗模型计算得到,路径损耗系数β是环境的函数,节点i接收到由节点j发出的信号强度测量值p<sub>ij</sub>的概率密度函数表示为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>ij</mi></msub><mo>;</mo><msub><mi>γ</mi><mi>i</mi></msub><mo>,</mo><msub><mi>γ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>10</mn><mo>/</mo><mi>log</mi><mn>10</mn></mrow><msqrt><mn>2</mn><mi>π</mi><msubsup><mi>σ</mi><mi>dB</mi><mn>2</mn></msubsup></msqrt></mfrac><mfrac><mn>1</mn><msub><mi>P</mi><mi>ij</mi></msub></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><mi>b</mi><mn>8</mn></mfrac><msup><mrow><mo>(</mo><mi>log</mi><mfrac><msubsup><mi>d</mi><mi>ij</mi><mn>2</mn></msubsup><msubsup><mover><mi>d</mi><mo>~</mo></mover><mi>ij</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000014.GIF" wi="1624" he="222" /></maths>其中<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><mn>10</mn><mi>β</mi></mrow><mrow><msub><mi>σ</mi><mi>dB</mi></msub><mi>log</mi><mn>10</mn></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000021.GIF" wi="408" he="265" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>=</mo><msub><mi>d</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>/</mo><msub><mi>p</mi><mi>ij</mi></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>β</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000022.GIF" wi="445" he="189" /></maths><img file="FDA0000715618420000023.GIF" wi="66" he="89" />是在给定接收信号强度p<sub>ij</sub>时,距离d<sub>ij</sub>的最大似然估计值,最终生成各个接入点AP与移动站MS的距离估计;提取接收信号强度RSS值最强的两个接入点AP与移动站MS的距离估计值并以d<sub>2</sub>,d<sub>3</sub>表示,在定位过程中,首先计算所画圆上各点与接收信号强度RSS值最强两个接入点AP的距离,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mn>2</mn><mi>θi</mi></mrow></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>Mi</mi></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>Mi</mi></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow>]]></math><img file="FDA0000715618420000024.GIF" wi="732" he="118" /></maths> (7)<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mn>3</mn><mi>θi</mi></mrow></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>Mi</mi></msub><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>Mi</mi></msub><mo>-</mo><msub><mi>y</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow>]]></math><img file="FDA0000715618420000025.GIF" wi="727" he="120" /></maths>并与距离估计做差,Err<sub>1i</sub>=|d<sub>2</sub>‑d<sub>2θi</sub>| (8)Err<sub>2i</sub>=|d<sub>3</sub>‑d<sub>3θi</sub>|代价函数最终表示为,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><munder><mi>Min</mi><mi>i</mi></munder><msqrt><msup><msub><mi>Err</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Err</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000715618420000026.GIF" wi="1558" he="130" /></maths>搜索整个圆上的代价函数,选取代价函数最小值所对应的移动站MS坐标值,作为移动站MS的位置估计值。 |