一种面向室内Wi-Fi定位的信道状态信息特征提取方法

一种面向室内wi
‑
fi定位的信道状态信息特征提取方法
技术领域
1.本发明属于室内定位技术领域，涉及一种面向室内wi
‑
fi定位的信道状态信息特征提取方法。


背景技术：

2.随着移动互联网时代向物联网时代逐步发展，位置服务(location
‑
based service，lbs)在人们生活中的应用已经越来越广泛。目前较为流行的无线定位系统有全球定位系统(global position system，gps)、蜂窝定位系统、蓝牙定位系统、射频识别定位系统、zigbee定位系统以及wlan定位系统。其中，由于gps定位系统在室外定位领域具有较好的表现，因此最为常见。然而，在室内环境下由于各种障碍物和移动物体的遮挡，卫星信号衰减十分严重，使得其室内定位精度不尽人意。相比而言，wi
‑
fi网络具有部署成本较低，环境适应性强且通信范围广等优势，于是基于接收信号强度(received signal strength，rss)的wi
‑
fi定位已逐渐成为室内定位技术的主流。
3.信道状态信息(channel state information，csi)是一种可替代rss的新兴技术，可用于指纹识别。相较于基于rss的室内wi
‑
fi定位方法，csi包含了信号传输过程中更细粒度和多样化的物理层信息，这些信息代表了散射、衰落和功率随距离衰减的综合效应，描述了信号是如何从发射机传播到接收机的。csi比传统的rss具有更高的鲁棒性，故基于csi的室内wi
‑
fi定位方法通常具有更高的定位精度且其定位结果更为稳定。基于wi
‑
fi的室内定位方法主要是通过个人电脑来获取csi，目前使用智能手机获取csi用于定位的研究较少，schulz等人首次提出在移动设备上使用nexmon固件补丁，修改wi
‑
fi固件并构建wi
‑
fi测试平台来提取csi。


技术实现要素：

4.有鉴于此，本发明的目的在于提供一种以多维度分析为主线的张量图像处理框架，分别实现了张量图像的降噪和特征提取，提升了数据处理分析的能力。
5.为达到上述目的，本发明提供如下技术方案：
6.一种面向室内wi
‑
fi定位的信道状态信息特征提取方法，以多维主成分分析为基础进行csi维度约减，并在尽可能保持原始csi数据变化信息的前提下，于张量各维度上提取最显著的信号成分，引入交替迭代最小二乘算法(alternating iterative least squares，ails)来求解投影矩阵的最优解，并得到最终的约减张量，最后，将重构得到的张量作为csi位置指纹实现特征提取。
7.进一步，本方法具体包括以下步骤：
8.s1：在n
f
个参考点(reference point，rp)处分别采集到来自同一个无线接入点(access point，ap)的信道状态信息(channel state information，csi)；
9.s2：对采集到的csi进行解析得到csi幅值，并表示为n
f
个x轴为子载波、y轴为数据包且z轴为csi幅值的三维图像，将所述三维图像视为n
f
个二元三维张量其
中，w1、w2和w3分别为子载波、数据包、幅值这三个维度的维度大小；
10.s3：将第n个参考点所对应的张量θ
n
划分为m个各维度大小相等且互不重叠的子张量
11.s4：对第n个参考点的第m个子张量a
(m,n)
进行去中心化得到其中，
12.s5：将去中心化子张量在第d(d＝1,...,3)维上展开为矩阵对协方差矩阵进行特征值分解后得到从大到小排序的特征值，初始化第d维投影矩阵为前d
d
个最显著的特征值所对应的特征向量的组合；
13.s6：利用交替迭代最小二乘算法ails来求解各维度投影矩阵的最优解和最终的约减张量
14.s7：重复步骤s3
‑
s6，得到第n个参考点的m个约减张量再将这m个约减张量按原位置重构得到张量并将其作为第n个参考点的csi位置指纹，其中，v1、v2和v3分别为子载波、数据包、幅值这三个维度的维度大小；
15.s8：重复步骤s3
‑
s7，得到n
f
个参考点的csi位置指纹实现特征提取。
16.进一步，步骤s6具体包括以下步骤：
17.s61：初始化子张量的约减张量为其能量为
18.s62：当第k次迭代时，先将张量在第d维(d＝1,...,3)上展开为矩阵再对协方差矩阵进行特征值分解后得到从大到小排序的特征值，然后更新第d维投影矩阵为前d
d
个最显著的特征值所对应的特征向量的组合，利用更新后的更新
19.s63：重复步骤s62依次更新第1维至第3维的投影矩阵后即完成一次迭代；
20.s64：获取第k次迭代的约减张量以及
该约减张量的能量
21.s65：如果相邻两次迭代所得能量的差值(η
(k,n)
‑
η
((k
‑
1),n)
)小于预先设定的阈值，则满足收敛要求，输出各维度的投影矩阵并将第k次迭代的约减张量作为最终的约减张量反之，重复步骤s62至步骤s64直到满足收敛要求。
22.本发明的有益效果在于：本发明充分利用了高阶张量可以描述数据信息和结构的特性，采用张量形式表达复杂数据，最大的优势是无需改变数据原有的存在形式，这也使得以张量形式表达的数据能够最大限度地保留数据内在的结构信息，提升了数据处理分析的能力。
23.本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
24.为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作优选的详细描述，其中：
25.图1为本发明所述面向室内wi
‑
fi定位的信道状态信息特征提取方法流程图；
26.图2为子张量在各维度展开为矩阵的示意图。
具体实施方式
27.以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
28.其中，附图仅用于示例性说明，表示的仅是示意图，而非实物图，不能理解为对本发明的限制；为了更好地说明本发明的实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
29.本发明实施例的附图中相同或相似的标号对应相同或相似的部件；在本发明的描述中，需要理解的是，若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此附图中描述位置关系的用语仅用于示例性说明，不能理解为对本发明的限制，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
30.请参阅图1～图2，如图1所示的一种面向室内wi
‑
fi定位的信道状态信息特征提取
方法，具体包含以下步骤：
31.步骤一、在n
f
个参考点(reference point，rp)处分别采集到来自同一个无线接入点(access point，ap)的信道状态信息(channel state information，csi)，再对采集到的csi进行解析得到csi幅值，并表示为n
f
个x轴为子载波、y轴为数据包且z轴为csi幅值的三维图像；
32.步骤二、将步骤一得到的三维图像视为n
f
个二元三维张量其中，w1、w2和w3分别为子载波、数据包、幅值这三个维度的维度大小。
33.步骤三、将第n(n＝1,...,n
f
)个参考点所对应的张量θ
n
划分为m个各维度大小相等且互不重叠的子张量，记为
34.步骤四、对第n个参考点的第m个子张量a
(m,n)
进行去中心化得到其中，
35.步骤五、将子张量在第d(d＝1,...,3)维上展开为矩阵对协方差矩阵进行特征值分解后得到从大到小排序的特征值，初始化第d维投影矩阵为前d
d
个最显著的特征值所对应的特征向量的组合；
36.步骤六、为了保留子张量的特征变化信息，利用ails算法来求解各维度投影矩阵的最优解和最终的约减张量具体包括以下步骤：
37.步骤六(一)、初始化子张量的约减张量为其能量为
38.步骤六(二)、当第k次迭代时，先将张量在第d维(d＝1,...,3)上展开为矩阵再对协方差矩阵进行特征值分解后得到从大到小排序的特征值，然后更新第d维投影矩阵为前d
d
个最显著的特征值所对应的特征向量的组合，利用更新后的更新
39.步骤六(三)、重复步骤六(二)依次更新第1维至第3维的投影矩阵后即完成一次迭代；
40.步骤六(四)、获取第k次迭代的约减张量以及该约减张量的能量如果相邻两次迭代所得能量的差值(η
(k,n)
‑
η
((k
‑
1),n)
)小于预先设定的阈值，则满足收敛要求，输出各维度的投影矩阵并将第k次迭代的约减张量作为最终的约减张量
41.步骤六(五)、反之，重复步骤六(二)至步骤六(四)直到满足收敛要求；
42.步骤七、重复步骤三至步骤六，可得到第n个参考点的m个约减张量再将这m个约减张量按原位置重构得到张量并将其作为第n个参考点的csi位置指纹，其中，v1、v2和v3分别为子载波、数据包、幅值这三个维度的维度大小；
43.步骤八、重复步骤三至步骤七，可得到n
f
个参考点的csi位置指纹
44.最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。