[发明专利]一种基于Matlab2008的室内指纹定位系统的仿真方法

摘要

一种基于Matlab2008的室内指纹定位系统的仿真方法，涉及一种室内指纹定位系统的仿真方法，属于室内指纹定位技术领域。它解决了现有的室内指纹定位系统中指纹图的建立工作量大，验证指纹图采集密度和定位算法等因素对定位结果的影响时采集数据量大的问题。本发明基于Matlab2008实现，包括两部分的内容。第一部分是室内电磁波传播模型的建立，利用它可以模拟室内任意一点的信号强度；第二部分是定位方法的实现，以第一部分的模型为基础，进行环境设置、指纹图采集、定位和误差统计，从而达到对算法性能的仿真验证。本发明主要适用于基于WLAN的室内指纹定位技术领域。

主权项

一种基于Matlab2008的室内指纹定位系统的仿真方法，其特征在于所述仿真方法的实现是按照以下步骤进行的：步骤A、室内电磁波传播模型的建立，即待测点信号强度的仿真：设待测点坐标为（x，y），AP，即WIFI接入点的坐标为（x_ap，y_ap），具体过程如下：步骤A1、计算待测点与AP之间的距离d：$d=\sqrt{{(x-x_{\mathrm{ap}})}^2+{(y-y_{\mathrm{ap}})}^2}---\left(1\right)$步骤A2、根据环境中墙和门的分布，判断待测点与AP之间隔着的墙的数目l和门的数目h，利用下式计算等效的穿墙数目DWR=DAF/WAF   （2）$\tilde{l}=(l+\mathrm{DWR}\&CenterDot;h)---\left(3\right)$其中DAF代表信号穿过一扇门时引起的衰减，衰减范围为10～20dB，WAF代表信号穿过一堵墙时引起的衰减，衰减范围为1.5～5.5dB，DWR的物理意义为：一扇门对信号的衰减作用相当于DWR堵墙；步骤A3、根据的大小来判断信号的传播：时，P_l(d₀)=41.5dB，n=1.9dB；时，P_l(d₀)=37.7dB，n=3.3dB，其中P_l(d₀)为第一米路径衰减，n为路径衰减因子；步骤A4、将步骤A1中得到的d、步骤A2中得到的步骤A3中得到的P_l(d₀)和n带入下面的式子，计算路径衰减P_l：$P_l=\left\{\begin{array}{c}{P}_l\left(d_0\right)+10\&CenterDot;n\&CenterDot;\log (d/d_0)+\tilde{l}\&CenterDot;\mathrm{WAF}(\tilde{l}\&le;2)\\ P_l\left(d_0\right)+10\&CenterDot;n\log (d/d_0)+2\&CenterDot;\mathrm{WAF}(\tilde{l}>2)\end{array}\right.---\left(4\right)$步骤A5、利用Matlab2008自带函数norm(·)生成一个均值为0、标准差为3的随机数，代表阴影衰落P_s；步骤A6、将步骤A4得到的路径衰减P_l、步骤A5得到的阴影衰落P_s，带入下面的式子计算待测点的接收信号强度rss：rss=P_t‑P_l‑P_s   （5）其中P_t为AP发射功率，范围为‑40dBm到0dBm之间；步骤B、定位方法的实现，具体过程如下：步骤B1、设定仿真环境大小为L米×W米，L>0,W>0，在这个范围内设定墙和门的位置，以及M个AP的位置；步骤B2、在步骤B1设定的环境中，均匀选取N个点，间隔为0.5到1米，设第i个点的坐标为（x_i，y_i），其中i=1，2，…，N；步骤B3、利用式（5）在步骤B2中的每个点（x_i，y_i）处仿真信号强度，在每个点仿真T次取平均值，得到第i个点的接收信号强度向量RSS_i，RSS_i=[rss_i1，rss_i2，…，rss_iM]，rss_i1表示第1个AP在第i个点处的平均信号强度，依次类推；步骤B4、设定定位过程中仿真次数为T_s次，在步骤B1的仿真环境里随机选取T_s个（x_p，y_p）坐标作为测试点，其中p=1，2，…，T_s，设置迭代次数计数值cnt的初值为0；步骤B5、设置迭代次数计数值cnt加1，在步骤B4得到的T_s个坐标中随机选择第cnt个坐标（x_cnt，_ycnt）作为要定位的点，利用式（5）仿真得到该点信号强度RSS₀，RSS₀=[rss₀₁，rss₀₂，…，rss_0M]，rss₀₁表示第l个AP在该点处的平均信号强度，依次类推；对RSS₀与步骤B3中每一个RSS_i，分别计算欧氏距离d_i：$d_i={\left(\underset{j=1}{\overset{M}{\mathrm{\&Sigma;}}}{(r{\mathrm{ss}}_{\mathrm{ij}}-{\mathrm{rss}}_{0j})}^2\right)}^{1/2}---\left(6\right)$步骤B6、在步骤B5从小到大排序的欧氏距离d_i中选取前k个，重新从小到大排序后得到欧氏距离其中q=1，2，…，k；对应的坐标为步骤B7、利用步骤B6中得到的坐标和欧氏距离得到定位结果$\hat{x}=\underset{q=1}{\overset{k}{\mathrm{\&Sigma;}}}f\left({\tilde{d}}_q\right){\tilde{x}}_q---\left(7\right)$$\hat{y}=\underset{q=1}{\overset{k}{\mathrm{\&Sigma;}}}f\left({\tilde{d}}_q\right){\tilde{y}}_q---\left(8\right)$其中，f(·)为被测试的定位方法；步骤B8、比较步骤B5给出的真实位置（x_cnt，y_cnt），步骤B7给出的定位结果得到定位误差err_cnt：${\mathrm{err}}_{\mathrm{cnt}}={({(\hat{x}-x_{\mathrm{cnt}})}^2+{(\hat{y}-y_{\mathrm{cnt}})}^2)}^{1/2}---\left(9\right)$步骤B9、判断迭代次数的计数值cnt是否达到T_s，若cnt<T_s，则回到步骤B5进行下一次迭代，否则，就统计定位误差，其结果为[err₁,err₂,…,err_Ts]，其中，err₁代表第1次仿真的定位误差，依此类推，共进行了T_s次，统计误差方差var和平均误差mean，用来判断被测方法f(·)的定位误差的大小：$\mathrm{var}=\frac{1}{T_s}({{\mathrm{err}}_1}^2+{{\mathrm{err}}_1}^2+...+{\mathrm{err}}_{T_s}^2)---\left(10\right)$$\mathrm{mean}=\frac{1}{T_s}({\mathrm{err}}_1+{\mathrm{err}}_2+...+{\mathrm{err}}_{T_s})---\left(11\right).$