发明名称 一种增强手机基站定位精度的方法
摘要 本发明属于手机基站定位领域,涉及一种利用隐马尔可夫地图匹配算法增强手机基站定位精度的方法。利用隐马尔可夫地图匹配算法来提高手机基站定位的精度,算法将定位点匹配到候选匹配道路上的概率称作初始匹配概率,将候选匹配道路发生转换的概率称作道路转移概率,算法使用动态编程来快速找到初始匹配概率和道路转移概率乘积最大的路径。考虑到车辆在实际行驶过程中不可能频繁地转换道路,我们在算法计算定位点道路转移概率时引入了约束限制进行改进,从而提高了算法性能。
申请公布号 CN104900059B 申请公布日期 2017.04.12
申请号 CN201510273848.6 申请日期 2015.05.26
申请人 大连理工大学 发明人 申彦明;张鹏飞
分类号 G08G1/01(2006.01)I;H04W64/00(2009.01)I;H04W4/02(2009.01)I;G01C21/30(2006.01)I 主分类号 G08G1/01(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 李宝元;梅洪玉
主权项 一种增强手机基站定位精度的方法,其特征在于如下步骤,(1)初始匹配概率点z<sub>t</sub>与每一条临近的道路r<sub>i</sub>存在一个匹配概率p(z<sub>t</sub>|r<sub>i</sub>);x<sub>t,i</sub>表示点z<sub>t</sub>匹配到道路r<sub>i</sub>上的点;点z<sub>1</sub>匹配到候选匹配道路r<sub>1</sub>,r<sub>2</sub>,r<sub>3</sub>上的候选匹配点为x<sub>1,1</sub>,x<sub>1,2</sub>,x<sub>1,3</sub>;坐标点和候选匹配点的大圆距离为‖z<sub>t</sub>‑x<sub>t,i</sub>‖<sub>great circle</sub>;相对于t+1时的点z<sub>t+1</sub>来说,z<sub>t+1</sub>匹配到道路r<sub>j</sub>上的点为x<sub>t+1,j</sub>;点x<sub>t,i</sub>和点x<sub>t+1,j</sub>之间车辆行驶的距离称为“路径距离”,记做‖x<sub>t,i</sub>‑x<sub>t+1,j</sub>‖<sub>route</sub>;模拟的定位误差为零均值的高斯误差,即:<maths num="0001"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>t</mi></msub><mo>|</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt><msub><mi>&sigma;</mi><mi>z</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mn>0.5</mn><msup><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>|</mo><msub><mo>|</mo><mrow><mi>g</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>t</mi><mi>c</mi><mi>i</mi><mi>r</mi><mi>c</mi><mi>l</mi><mi>e</mi></mrow></msub><mo>/</mo><msub><mi>&sigma;</mi><mi>z</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194364350000011.GIF" wi="1261" he="220" /></maths>σ<sub>z</sub>在式中代表定位测量值的标准偏差;初始匹配概率π<sub>i</sub>,i=1,…,N<sub>r</sub>,表示在地图匹配开始运行前,从所有路段中找出车辆所在路段的概率;当部分HMM的构想值均匀分布于π<sub>i</sub>上时,假设没有测量值被采用,则在第一个测量值上开始,有π<sub>i</sub>=p(z<sub>t</sub>|r<sub>i</sub>),使用第一测量值z<sub>1</sub>;(2)道路转移概率转移概率是指在这两个时间t到t+1内,车辆在候选匹配道路之间的转移概率;对于测量值z<sub>t</sub>和其候选匹配道路段r<sub>i</sub>来说,在r<sub>i</sub>上最可能匹配点是x<sub>t,i</sub>;测量值z<sub>t+1</sub>在其候选匹配道路段r<sub>j</sub>上的候选匹配点是x<sub>t+1,j</sub>;测量点之间的“路径距离”||x<sub>t,i</sub>‑x<sub>t+1,j</sub>||<sub>route</sub>和大圆距离||z<sub>t</sub>‑z<sub>t+1</sub>||<sub>greatcircle</sub>,正确匹配时这两个值相当接近,通过比较这两个值判断道路是否发生转换;通过对正确匹配的大圆距离和路径距离差的绝对值进行分析,得出指数概率分布:<maths num="0002"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>d</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>&beta;</mi></mfrac><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>d</mi><mi>t</mi></msub><mo>/</mo><mi>&beta;</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194364350000012.GIF" wi="1342" he="126" /></maths>式中:<maths num="0003"><math><![CDATA[<mrow><msub><mi>d</mi><mi>t</mi></msub><mo>=</mo><mo>|</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>z</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>|</mo><mo>|</mo></mrow><mi>great circle</mi></msub><mo>-</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><msup><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow><mo>*</mo></msup></msub><mo>-</mo><msub><mi>x</mi><msup><mrow><mi>t</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mo>*</mo></msup></msub><mo>|</mo><mo>|</mo></mrow><mi>route</mi></msub><mo>|</mo></mrow>]]></math><img file="FDA0001194364350000013.GIF" wi="1032" he="118" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>&beta;</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac><msub><mi>median</mi><mi>t</mi></msub><mrow><mo>(</mo><mo>|</mo><mrow><mo>|</mo><mo>|</mo><msub><mi>z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>z</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>|</mo><msub><mo>|</mo><mrow><mi>g</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>t</mi><mi>c</mi><mi>i</mi><mi>r</mi><mi>c</mi><mi>l</mi><mi>e</mi></mrow></msub><mo>-</mo><mo>|</mo><mo>|</mo><msub><mi>x</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi><mo>*</mo></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>*</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mrow><mi>r</mi><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi></mrow></msub></mrow><mo>|</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194364350000021.GIF" wi="1124" he="126" /></maths>其中,i<sup>*</sup>和j<sup>*</sup>指的是车辆的真实行驶道路段;(3)优化引入一个参数Δs来对这一行为进行约束;当定位点的候选匹配道路发生改变,候选匹配点离开原行驶的道路,其路径距离就远离大圆距离50米进行惩罚,就是当路径距离大于大圆距离时,路径距离加上50米,路径距离小于大圆距离时,路径距离减去50米;这样不正确匹配点的路径距离就会远离大圆距离,更容易找出正确匹配点;(4)输出最佳匹配路径通过用公式(1)中的匹配概率和公式(2)中的转换概率,使用维特比算法在HMM点阵中来计算出最佳匹配路径;维特比算法使用动态编程来快速找到匹配概率和转移概率乘积最大的路径,从而推断出坐标点的正确匹配道路段;如果对于从t时刻的点z<sub>t</sub>匹配到道路r<sub>i</sub>上的点x<sub>t,i</sub>到达t+1时刻的点z<sub>t+1</sub>匹配到r<sub>j</sub>道路上的点x<sub>t+1,j</sub>发生道路转换时,要考虑点z<sub>t</sub>是通过哪条道路段r<sub>x</sub>到达点z<sub>t+1</sub>,不能只单独考虑道路匹配概率p(z<sub>t</sub>|r<sub>i</sub>)和p(z<sub>t+1</sub>|r<sub>j</sub>),要考虑两点之间发生的道路转换概率,即匹配概率乘以相应的道路转换概率,从中选出最大值即为此时点的最佳匹配道路段,并记下该点在此道路段上的最佳匹配值是由哪条道路段转换而来,以此类推;算法对定位点的处理区分为对第一个定位点的处理和对普通定位点的处理;具体如下:对于一个点z<sub>t</sub>,它匹配到每条候选道路段上的状态概率为p(z<sub>t</sub>|r<sub>i</sub>);如果该点是车辆定位点的第一个点,则不需要对此概率进行修正,如果该点不是该车辆定位点的第一个点,仅为分组后该组点的第一个点,则要对状态概率p(z<sub>t</sub>|r<sub>i</sub>)进行修正;需要将状态概率乘以上时刻点的匹配道路段到本时刻点的匹配路段的转移概率,即p(z<sub>t</sub>|r<sub>i</sub>)=p(z<sub>t</sub>|r<sub>i</sub>)×p(d<sub>t‑1</sub>);对于其余点的处理,只需考虑每个点z<sub>t</sub>匹配到每条可行弧r<sub>i</sub>的状态概率;设p(z<sub>t</sub>|r<sub>i</sub>)表示第t个点z<sub>t</sub>匹配到弧r<sub>i</sub>的最大概率,设从弧r<sub>i</sub>转移到弧r<sub>j</sub>的转移概率为p(r<sub>i</sub>|r<sub>j</sub>),第t+1个点定位到弧r<sub>j</sub>的状态概率为p(z<sub>t+1</sub>|r<sub>j</sub>),则其转移方程为p(z<sub>t+1</sub>|r<sub>j</sub>)=max(p(z<sub>t+1</sub>|r<sub>j</sub>),p(z<sub>t</sub>|r<sub>i</sub>)×p(z<sub>t+1</sub>|r<sub>j</sub>)×p(r<sub>i</sub>|r<sub>j</sub>));如果p(z<sub>t+1</sub>|r<sub>j</sub>)取值为p(z<sub>t</sub>|r<sub>i</sub>)×p(z<sub>t+1</sub>|r<sub>j</sub>)×p(r<sub>i</sub>|r<sub>j</sub>),则需要记录转换路径,在输出结果时以此反推最佳匹配道路段;当计算完全部点的状态概率后,确定最后一个点的最佳匹配道路段,并以此道路段反向推出从第一个点到最后一个点的车辆实际行驶轨迹的道路段;整个行驶轨迹计算完毕后,最终输出车辆的最佳匹配路径。
地址 116024 辽宁省大连市甘井子区凌工路2号