发明名称 一种基于非现场勘测的无线室内定位方法
摘要 本发明涉及一种基于非现场勘测的无线室内定位方法,属于无线定位技术领域。本方法首先利用智能手机自动采集指纹数据,形成数据集;基于数据集中无线信号的二阶差分值,利用K平均算法生成虚拟房间;根据传感器读数提取用户移动路径,从而获得虚拟房间的连通关系,构建逻辑平面图;将逻辑平面图映射到实际平面图上,得到指纹-位置关系;匹配用户提交的位置请求所在的虚拟房间,并检索指纹-位置关系数据库,判断用户所在位置。本发明方法无需对定位区域进行人工的现场勘测,无需无线接入点的位置、功率等信息,充分利用人的移动信息,实现房间级别的逻辑定位,给出用户当前位置,既提高了定位的精度与实时性,又节约了使用者的时间。
申请公布号 CN102724751A 申请公布日期 2012.10.10
申请号 CN201210181296.2 申请日期 2012.06.04
申请人 清华大学 发明人 吴陈沭;赵弋洋;杨铮;刘云浩
分类号 H04W64/00(2009.01)I 主分类号 H04W64/00(2009.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种基于非现场勘测的无线室内定位方法,其特征在于该方法包括以下各步骤:(1)通过手机采集无线网络的信号以及手机中的加速度传感器和方向传感器的读数,组成一个指纹数据集:D={F,S},其中F为采集的无线信号,F=[f<sub>1</sub>,f<sub>2</sub>,f<sub>3</sub>,......,f<sub>n</sub>],f<sub>i</sub>表示无线网络中第i个无线接入点的信号强度值,S为手机中的加速度传感器和方向传感器的读数,i=1,2,3,......,n,n为无线网络中无线接入点的个数;(2)设指纹数据集中的任意两条无线信号分别为:F=[f<sub>1</sub>,f<sub>2</sub>,f<sub>3</sub>,......,f<sub>n</sub>]和F’=[f’<sub>1</sub>,f’2,f’3,...,f’n],两个无线信号之间的区分度用无线信号强度的二阶差分值表示为:<maths num="0001"><![CDATA[<math><mrow><mi>&phi;</mi><mrow><mo>(</mo><mi>F</mi><mo>,</mo><msup><mi>F</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>=</mo><msqrt><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msup><mrow><mo>(</mo><mi>&omega;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><msub><msup><mi>f</mi><mo>&prime;</mo></msup><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>&omega;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>I</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>-</mo><msub><mi>f</mi><mi>j</mi></msub><mo>></mo><mn>0</mn><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>-</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,I为指示型函数;(3)根据步骤(2)中的二阶差分值,采用K平均聚类方法将指纹数据集中的无线信号划分为N类,即N个虚拟房间<img file="FDA00001722517400013.GIF" wi="478" he="48" />从K平均聚类结果中得到每个虚拟房间R的平均指纹F[R]和指纹阈值ξ,其中平均指纹F[R]为虚拟房间R所包含的所有无线信号指纹的平均值,指纹阈值ξ为该虚拟房间R所包含的所有无线信号中与平均指纹F[R]的最大的二阶差分值,<img file="FDA00001722517400014.GIF" wi="141" he="40" />(4)构建一个逻辑平面图G=(V,E),使逻辑平面图G=(V,E)中的每一个顶点R<sub>1</sub>∈V代表一个虚拟房间,每一条边e=(R<sub>1</sub>,R<sub>2</sub>)∈E表示虚拟房间R<sub>1</sub>与R<sub>2</sub>相连通,构建过程包括如下步骤:(4-1)根据步骤(1)中的连续采集的传感器读数A,形成手机在移动过程中的连续移动轨迹,移动轨迹上的无线信号强度依次记为<img file="FDA00001722517400015.GIF" wi="445" he="47" />与该连续无线信号强度相对应的步骤(3)的多个相连续的虚拟房间记为<img file="FDA00001722517400016.GIF" wi="479" he="48" />(4-2)根据步骤(4-1)中的连续虚拟房间序列<img file="FDA00001722517400017.GIF" wi="479" he="49" />将每个虚拟房间简化为一个点,在任意两个相邻的虚拟房间R<sub>i</sub>、R<sub>i+1</sub>之间连接出一条边;(4-3)遍历步骤(1)的指纹数据集中的所有连续采集的传感器读数,重复步骤(4-1)和(4-2),得到逻辑平面图G=(V,E);(5)将上述逻辑平面图映射到需无线室内定位区域的实际平面图,得到逻辑平面图中每一个顶点与实际平面图的对应关系,即指纹-位置关系数据库,包括如下步骤:(5-1)将需无线室内定位区域的实际平面图简化为一个无向图G′=(V′,E′),其中每个顶点A<sub>1</sub>∈V′代表需无线室内定位区域的实际平面图中的一个功能区域,每条边e=(A<sub>1</sub>,A<sub>2</sub>)∈E′表示两个功能区域A<sub>1</sub>和A<sub>2</sub>相连通;(5-2)根据介数中心性,建立上述逻辑平面图G=(V,E)和实际平面图G′=(V′,E′)中居于中心的顶点的相互映射关系,逻辑平面图G=(V,E)中顶点R的介数中心性C<sub>b</sub>(R)为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>C</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>s</mi><mo>&NotEqual;</mo><mi>t</mi><mo>&NotEqual;</mo><mi>R</mi><mo>&Element;</mo><mi>V</mi></mrow></munder><mfrac><mrow><msub><mi>&sigma;</mi><mi>st</mi></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow></mrow><msub><mi>&sigma;</mi><mi>st</mi></msub></mfrac></mrow></math>]]></maths>其中,σ<sub>st</sub>为逻辑平面图G=(V,E)中顶点s到顶点t的最短路径数,σ<sub>st</sub>(R)为穿过顶点R的顶点s到顶点t的最短路径数,按介数中心性从大到小排列得到逻辑平面图G=(V,E)的中心顶点序列{r<sub>1</sub>,r<sub>2</sub>,......,r<sub>k</sub>};实际平面图G′=(V′,E′)中顶点A的介数中心性C<sub>b</sub>(A)为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>C</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>s</mi><mo>&NotEqual;</mo><mi>t</mi><mo>&NotEqual;</mo><mi>A</mi><mo>&Element;</mo><msup><mi>V</mi><mo>&prime;</mo></msup></mrow></munder><mfrac><mrow><msub><mi>&sigma;</mi><mi>uv</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow></mrow><msub><mi>&sigma;</mi><mi>uv</mi></msub></mfrac></mrow></math>]]></maths>其中,σ<sub>uv</sub>为实际平面图G′=(V′,E′)中顶点u到顶点v的最短路径数,σ<sub>uv</sub>(A)为穿过顶点A的顶点u到顶点v的最短路径数,按介数中心性从大到小排列得到实际平面图G′=(V′,E′)的中心顶点序列{a<sub>1</sub>,a<sub>2</sub>,......,a<sub>k</sub>};将上述中心顶点序列{r<sub>1</sub>,r<sub>2</sub>,......,r<sub>k</sub>}与中心顶点序列{a<sub>1</sub>,a<sub>2</sub>,......,a<sub>k</sub>}按照序号一一对应,得到逻辑平面图G=(V,E)和实际平面图G′=(V′,E′)的中心顶点的映射关系;(5-3)设逻辑平面图G=(V,E)中非中心顶点v到所有其他顶点的最短路径长度的和为该顶点v的权值,即:w(v)=∑<sub>u∈V,u≠v</sub>d(v,u)其中,d(v,u)为顶点v到u的最短路径长度,重复本步骤,得到逻辑平面图G=(V,E)中所有非中心顶点的权值;实际平面图G′=(V′,E′)中非中心顶点v′到所有其他顶点的最短路径长度的和为该顶点v′的权值,即:w(v′)=∑<sub>u′∈V,u′≠v′</sub>d(v′,u′)其中,d(v′,u′)为顶点v′到u′的最短路径长度,重复本步骤,得到实际平面图G′=(V′,E′)中所有非中心顶点的权值;根据上述得到的所有非中心顶点的权值,利用KM算法求解带权最小二分图匹配模型,使得所有匹配顶点之间的权值之差最小,即使得W=∑<sub>v∈V</sub>|w(v)-w(v′)|最小,其中v′为v的映射,得到逻辑平面图和实际平面图中的非中心顶点的映射关系;(5-4)根据步骤(5-2)和(5-3)中得到的映射关系,实施以下步骤加以纠正:设逻辑平面图G=(V,E)中顶点p被映射到实际平面图G′=(V′,E′)中的顶点p′,定义顶点v的邻居集合N<sub>G</sub>(p)为逻辑平面图G=(V,E)中与顶点p直接相连的顶点集合,定义顶点p′的邻居集合N<sub>G′</sub>(p′)为实际平面图G′=(V′,E′)中与顶点p′直接相连的顶点集合,设N′<sub>G</sub>(p)为与顶点p的邻居集合N<sub>G</sub>(p)中的所有节点相映射的实际平面图G′=(V′,E′)中的节点集合,若N′<sub>G</sub>(p)与N<sub>G′</sub>(p′)中的相同节点数少于N<sub>G</sub>(p)所包含节点数的一半,则判定p到p′的映射不正确;遍历实际平面图G′=(V′,E′)除p′以外的所有顶点,寻找新顶点q,使得N′<sub>G</sub>(p)与q的邻居集合N<sub>G′</sub>(q)的相同顶点数最多,并确定顶点q与顶点p的映射关系;(5-5)重复步骤(5-4),得到逻辑平面图G=(V,E)与实际平面图G′=(V′,E′)之间的映射关系;(6)根据步骤(5)得到的指纹-位置关系数据库,实现对用户的定位,包括以下步骤:(6-1)根据用户提交的无线信号为F,若存在满足以下条件的虚拟房间R:<img file="FDA00001722517400031.GIF" wi="796" he="50" />与<maths num="0005"><![CDATA[<math><mrow><mi>&phi;</mi><mrow><mo>(</mo><mi>F</mi><mo>,</mo><mi>F</mi><mo>[</mo><mi>R</mi><mo>]</mo><mo>)</mo></mrow><mo>&lt;</mo><mi>&xi;</mi></mrow></math>]]></maths>其中<img file="FDA00001722517400033.GIF" wi="28" he="42" />为步骤(2)中定义的二阶差分值,R<sub>i</sub>为任意一个虚拟房间,<img file="FDA00001722517400034.GIF" wi="38" he="36" />为步骤(3)中的虚拟房间集合,F[R]为虚拟房间R的平均信号强度,ξ为步骤(2)中的与虚拟房间R对应的信号强度阈值,则判定用户所在虚拟房间为R;(6-2)根据上述判定的用户所在虚拟房间R,检索步骤(5)中指纹位置数据库,得到与上述虚拟房间R所对应的实际平面图中的功能区域,即为用户所在的位置。
地址 100084 北京市海淀区清华园1号