[发明专利]一种基于WiFi接收信号强度的混合室内定位方法在审

专利信息
申请号: 201410682677.8 申请日: 2015-08-04
公开(公告)号: CN104507159A 公开(公告)日: 2015-07-29
发明(设计)人: 杨东凯;修春娣;杨萌;刘源;罗智勇 申请(专利权)人: 北京航空航天大学
主分类号: H04W64/00 分类号: H04W64/00
代理公司: 北京慧泉知识产权代理有限公司 11232 代理人: 王顺荣;唐爱华
地址: 100191*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 一种基于WiFi接收信号强度的混合室内定位方法,它有五大步骤:步骤一:根据位置指纹室内定位技术,在离线阶段,根据选定的室内测试区域,依照对数距离路径损耗模型产生RSSI值,构成离线指纹数据库;步骤二:通过航位推测法在二维坐标下产生目标的运动路径;步骤三:通过模糊k-NN算出定位结果;步骤四:根据卡尔曼滤波方法的状态模型对运动目标的下一个位置进行预测;步骤五:将通过位置指纹技术得到的估计位置坐标加入到卡尔曼滤波方法的量测矩阵中,得到混合量测矩阵,通过混合量测模型对系统状态进行更新,得到滤波后的估计位置坐标;通过步骤四与步骤五不断的迭代更新,得到整个运动路径下的定位结果。
搜索关键词: 一种 基于 wifi 接收 信号 强度 混合 室内 定位 方法
【主权项】:
一种基于WiFi接收信号强度的混合室内定位方法,其特征在于:它包括以下步骤:步骤一:根据位置指纹室内定位技术,在离线阶段,根据选定的室内测试区域,依照对数距离路径损耗模型产生RSSI值,构成离线指纹数据库;其中,由于信号的长距离衰落特性服从对数正态分布,常用对数距离路径损耗模型表示,其中对数距离路径损耗模型的具体表示如下:<mrow><msub><mi>P</mi><mi>rx</mi></msub><mrow><mo>(</mo><mi>dB</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><msub><mi>d</mi><mi>o</mi></msub></msub><mrow><mo>(</mo><mi>dBm</mi><mo>)</mo></mrow><mo>-</mo><mn>10</mn><mi>l</mi><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><mfrac><mi>R</mi><msub><mi>R</mi><mi>o</mi></msub></mfrac><mo>)</mo></mrow><mo>-</mo><mi>wallLoss</mi></mrow>其中,Prx表示接收强度,是相对距离的接收强度或一米远的初始RSSI值,l是路径损耗指数,它能随传播环境的不同在2~6之间变化;Ro参考距离,R是接收设备与发射设备之间的距离;wallLoss是指由各面墙造成的损耗的和,这个因素决定于建筑布局,建筑材料,大量反射面,基础公共设施和移动物体;步骤二:通过航位推测法在二维坐标下产生目标的运动路径;xi=xi‑1+viΔtcos(θi)yi=yi‑1+viΔtsin(θi)其中,xi,yi表示当前物体的位置,xi‑1,yi‑1是笛卡尔坐标系统中前一个参考位置,vi和Δt分别表示速度和取样间隔;步骤三:通过模糊k‑NN算出定位结果;其中,离线阶段从各AP采集的接收信号强度RSSI向量为在线阶段,在特定位置采集到的RSSI值为:其中,N为AP的个数,n为测试区域栅格的个数;模糊k‑NN分类器采用一阶最近邻算法,则在位置γi的值μi(γ)表示为:<mrow><msub><mi>&mu;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>/</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><mi>&xi;</mi><mo>|</mo><mo>|</mo></mrow><mrow><mn>2</mn><mo>/</mo><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><mi>&xi;</mi><mo>|</mo><mo>|</mo></mrow><mrow><mn>2</mn><mo>/</mo><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow></mrow></mfrac></mrow>其中,其中K是最近邻的数量,模糊强度参数m用来决定计算各个邻居对函数值的贡献时距离的权重由多大,m的取值范围为(1,+∞);步骤四:根据卡尔曼滤波方法的状态模型对运动目标的下一个位置进行预测;其中,系统的状态方程是基于航位推测模型的,其中当前位置表示为x=[px py vx vy]T,则系统的状态方程的表达式如下:<mrow><msub><mi>A</mi><mi>k</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>T</mi><mi>k</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>T</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>由于对数路径损耗模型是非线性的,所以采用如下非线性模型进行预测:xk+1=fk(xk,k)+wkyk=hk(xk)+vk其中,wk为系统噪声,vk为量测噪声;系统噪声和量测噪声均为高斯白噪声,系统噪声的协方差矩阵为量测噪声的协方差矩阵为系统参数的线性一阶近似定义为雅可比矩阵<mrow><msup><mi>&Phi;</mi><mrow><mo>[</mo><mn>1</mn><mo>]</mo></mrow></msup><mrow><mo>(</mo><mover><mi>x</mi><mo>^</mo></mover><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><mo>&PartialD;</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub></mrow></msub><mo>,</mo><msup><mi>H</mi><mrow><mo>[</mo><mn>1</mn><mo>]</mo></mrow></msup><mrow><mo>(</mo><mover><mi>x</mi><mo>^</mo></mover><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><mo>&PartialD;</mo><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub></mrow></msub><mo>;</mo></mrow>则系统的预测过程如下:<mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>=</mo><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>=</mo><msub><mi>&Phi;</mi><mi>k</mi></msub><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>+</mo></msubsup><msubsup><mi>&Phi;</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub></mrow>步骤五:将通过位置指纹技术得到的估计位置坐标加入到卡尔曼滤波方法的量测矩阵中,得到混合量测矩阵,通过混合量测模型对系统状态进行更新,得到滤波后的估计位置坐标;通过步骤四与步骤五不断的迭代更新,得到整个运动路径下的定位结果;其中,步骤三中通过位置指纹方法得到的定位结果的坐标为则混合量测矩阵的表达式如下所示:<mrow><msub><mi>y</mi><mi>k</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>LF</mi><mi>x</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>P</mi><mi>LF</mi><mi>y</mi></msubsup></mtd></mtr><mtr><mtd><mi>y</mi><msub><mi>P</mi><msub><mi>rx</mi><mn>1</mn></msub></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mi>y</mi><msub><mi>P</mi><msub><mi>rx</mi><mn>1</mn></msub></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>LF</mi><mi>x</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>P</mi><mi>LF</mi><mi>y</mi></msubsup></mtd></mtr><mtr><mtd><msub><mi>P</mi><msub><mi>R</mi><mi>o</mi></msub></msub><mo>-</mo><mn>10</mn><mi>l</mi><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mi>o</mi></msub></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>P</mi><msub><mi>R</mi><mi>o</mi></msub></msub><mo>-</mo><mn>10</mn><mi>l</mi><msub><mi>log</mi><mn>10</mn></msub><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>N</mi></msub><msub><mi>R</mi><mi>o</mi></msub></mfrac><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow>通过卡尔曼滤波的量测模型进行更新,更新过程如下:<mrow><msub><mi>K</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>H</mi><mi>k</mi></msub><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow><mrow><msub><mi>I</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>y</mi><mi>k</mi></msub><mo>+</mo><msub><mi>H</mi><mi>k</mi></msub><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup></mrow><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>I</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>-</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>H</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup></mrow>其中,Kk为卡尔曼增益,Pk为误差协方差矩阵,Ik+1为真实值与量测值之间的误差,通过不断的迭代更新得到最终的定位结果。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410682677.8/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top