一种无接触感知定位方法与流程
本发明涉及感知定位技术，尤其涉及一种基于商用wifi设备的室内无接触感知定位方法。
背景技术：
接触感知定位指的是在定位目标不携带任何设备，也不与设备发生任何接触行为的情况下获得目标的位置。相比于接触式定位方式，无接触感知定位具有非侵扰性、方便、低成本的优点。而商用wifi设备目前广泛存在于我们的日常生活中，因此实现商用wifi设备上的无接触感知定位技术受到了广泛关注。
自2007年moustayoussef首次提出设备无关被动定位(device-freepassivelocalization,dfp)的概念以来，基于wifi的非接触式室内定位引起了很多学者的关注，涌现出一批优秀的工作。早期的工作主要利用wifirssi实现(youssefetal.,2007；zhangetal.,2007；seifeldinetal.,2013；xuetal.,2013)，尽管rssi可以在相对简单的空旷室内环境达到米级的定位精度，但是其性能并不可靠，研究显示，即使室内的空间布局不发生任何变化，在没有人的环境下，rssi依旧可以达到5db的波动(wuetal.,2012)。这种难以预测的不确定性波动，大大制约了基于rssi的室内定位的可靠性。因此大部分基于rssi可靠的定位方式也只能停留在房间级别的精度(yangetal.,2013)。自2011年danielhalperin首次在sigcomm公开基于intel5300网卡的csi测量方法和工具(halperinetal.,2011)，csi以其信号度量的精细化优势，基于wificsi的非接触式室内定位的精度得到前所未有的提升，如图1-5所示，从pilot(xiaoetal.,2013),monophy(abdel-nasseretal.,2013)到lifs(wangetal.,2016b),matrack(lietal.,2016b)，达到了米级以下的精度。lifs利用11个设备(4个ap配合7个wifi接收器)实现70m2环境中1.1m的定位中位误差，但是其位置的映射关系依然依赖指纹，且nlos定位精度较低；matrack利用多天线aoa技术直接采取几何映射的方案，在36m2环境中实现了0.6m的定位中位误差，但是其困难在于，一方面，基于静态路径融合虽然降低了需要识别的路径数量，但是也导致在融合角度附近造成感知盲区；另一方面，基于两维空间谱搜索算法带来了实时性难以保障的困难，同时累计误差依然是基于aoa定位的局限。因此，如何实现低成本(指纹库构建成本、设备数量的成本)、高精度的定位方案依然是室内定位领域尚未解决的问题。
菲涅尔区的相关概念：当射频信号由发送设备发送后，其在空间中通过多条路径达到接收端，也就是接收端接收到的信号是多个信号的叠加，根据波的干涉原理，当信号的两条路径长度相差半波长(λ/2)的奇数倍时产生减弱的信号，当信号的两条路径长度相差半波长的偶数倍时产生增强的信号。在本实例中，当发送设备和接收设备固定时，其直接路径(los)的长度是固定的，考虑上反射造成的相位偏转π，当反射路径与直接路径的长度差为半波长的奇数倍时产生增强的信号，为半波长的偶数倍时产生减弱的信号，这些位置就是菲涅尔区的边界。当反射路径与直接路径的长度差为i*λ/2时，即|tqi|+|qir|-|tr|＝iλ/2，其中，|tqi|是qi到发射设备t的距离，|qir|是qi到接收设备r的距离，|tr|为收发设备之间的距离，λ为载波的波长，所有满足该条件的点构成的椭圆为第i菲涅尔区的边界。由第1菲涅尔区边界围成的椭圆称为第1菲涅尔区，以此类推，由第i菲涅尔区边界和第i+1菲涅尔区边界围成的椭圆环为第i+1菲涅尔区。图1和图2分别示出了一维方向和二维方向上的菲涅尔区、菲涅尔边界以及交点区域。
当发送设备和接收设备的位置固定后，不同频率的射频信号根据自己的波长形成各自独立的菲涅尔区。
技术实现要素：
为了克服上述现有技术的不足，本发明提供一种室内无接触感知定位方法，基于多载波穿越模型，使用商用wifi设备，在目标不携带任何设备的情况下利用无线射频信号(rf)对目标在大范围内进行定位，并且不需要修改商用wifi设备的任何硬件信息，最少使用三个wifi接收设备就能获得室内环境中的一个运动物体的位置。
为了便于说明，本文约定：“rf”表示无线射频信号，csi表示信道状态信息。
本发明的原理是：rf信号分布空间为椭圆层次形状的空间菲涅尔区，且椭圆的大小由rf信号的波长决定；当给定两个载波λa和λb，动态反射路径长度为
在本发明的一个方面中，提供一种确定环境中的菲涅尔相位偏移的方法，其包括：接收设备r从发送设备t接收射频信号，所述射频信号包括不同波长的至少两个载波；以发送设备t的位置、接收设备r的位置为椭圆焦点，构建所述至少两个载波的k层菲涅尔区；在所述环境中布置反射目标，使得反射目标在所述环境中移动；针对所述至少两个载波中的至少一个载波对组合，确定所述反射目标引起的不同载波之间的理论菲涅尔相位差；针对所述至少两个载波中的所述至少一个载波对组合，确定所述反射目标和环境引起的不同载波之间的实际菲涅尔相位差；将实际菲涅尔相位差与相应的理论菲涅尔相位差进行比较，确定在所述环境中的所述至少一个载波对组合的菲涅尔相位偏移。在一个实施例中，在确定一个载波对的菲涅尔相位偏移时，确定理论菲涅尔相位差和实际菲涅尔相位差针对相同的载波对。
在一个实施例中，所述至少一个载波对组合包含所述至少两个载波中的所有载波对组合。
在一个优选实施例中，所述反射目标沿着垂直于发送设备t和接收设备r的方向移动。
在一个实施例中，所述反射目标为强反射目标，其使得到达其自身的射频信号的大部分被反射出去。
在一个实施例中，所述反射目标为金属。
在一个优选实施例中，所述反射目标为铁板。
在一个实施例中，所述反射目标为人体。
在一个优选的实施例中，确定在所述环境中的菲涅尔相位偏移的步骤包括：针对一段时间内接收到的信号，确定每对载波在所有固定大小的时间窗口内的菲涅尔相位偏移，取所有时间窗口的平均值或中位数作为所述至少一个载波对组合之间的菲涅尔相位偏移。
在本发明的一个方面中，提供一种利用以上所述的确定环境中的相位偏移的方法进行室内定位的方法，其包括：第一接收设备r1和第二接收设备r2同时从发送设备t接收射频信号，所述射频信号包含两个或更多个载波；针对所述两个或更多个载波中的至少一对载波a，以发送设备t的位置、第一接收设备r1的位置为椭圆焦点，构建m层菲涅尔区；针对所述两个或更多个载波中的至少一对载波b，以发送设备t的位置、第二接收设备r2的位置为椭圆焦点，构建n层菲涅尔区；在所述m层菲涅尔区中，根据实际菲涅尔相位差和菲涅尔相位偏移，确定目标位于的菲涅尔区f1；在所述n层菲涅尔区中，根据实际菲涅尔相位差和菲涅尔相位偏移，确定目标位于的菲涅尔区f2；确定菲涅尔区f1和f2的交点；根据菲涅尔区f1和f2的交点，确定目标的位置坐标。
在本发明的一个方面中，提供一种确定环境中的菲涅尔相位偏移的系统，其包括：接收设备r，其经配置执行射频信号接收；移动的反射目标，其经配置在所述环境中移动；确定菲涅尔相位偏移的设备，其经配置根据从发送设备t接收的包括不同波长的至少两个载波的射频信号，执行以下操作：以发送设备t的位置、接收设备r的位置为椭圆焦点，构建所述至少两个载波的k层菲涅尔区；针对所述至少两个载波中的至少一个载波对，确定所述反射目标引起的不同载波之间的理论菲涅尔相位差；针对所述至少两个载波中的所述至少一个载波对，确定所述反射目标和环境引起的不同载波之间的实际菲涅尔相位差；将实际菲涅尔相位差与相应的理论菲涅尔相位差进行比较，确定在所述环境中的菲涅尔相位偏移。在一个实施例中，在确定一个载波对的菲涅尔相位偏移时，确定理论菲涅尔相位差和实际菲涅尔相位差针对相同的载波对。
在一个实施例中，将将实际菲涅尔相位差与相应的理论菲涅尔相位差进行比较，确定在所述环境中的菲涅尔相位偏移的步骤进一步包括：将实际菲涅尔相位差与相应的理论菲涅尔相位差进行减法运算，将减法运算的差值作为在所述环境中的菲涅尔相位偏移。
在一个实施例中，所述确定相位偏移的系统选自包含所述接收设备r或任意能够执行所述操作的外部设备的组中的一个。
在一个实施例中，所述射频信号包括wi-fi信号。
根据本发明的一个方面，提供一种计算接可读存储介质，其存储有经配置由处理器执行的指令，该指令使得计算机执行上述确定相位偏移的方法所述的操作。
与现有技术相比，本发明的有益效果是：
本发明提供一种无接触感知定位方法，利用无线射频信号对移动目标定位，定位目标无需携带任何设备，不需要修改wifi无线网卡的任何硬件信息，仅利用商用wifi网卡就可以实现准确的室内无接触式动态目标定位，最少只需要一个wifi发射端和两个wifi接收端，就能够确定目标的位置坐标。
此外，由于本发明无需携带或接触任何设备，且射频信号(如wifi)已经广泛存在，本发明提供的技术方案具有非侵扰性、方便、低成本的优点。
附图说明
图1示出菲涅尔区、菲涅尔边界和同焦点椭圆的方程；
图2示出菲涅尔区、菲涅尔区交点区域以及目标位置所在的交点区域；
图3是本发明的一个实施例的确定菲涅尔相位偏移的系统，其中，(a)为发送设备(其支持包括wi-fi信号的射频信号发射，例如wi-fi设备，)；(b)为接收设备(其支持包括wi-fi信号的射频信号接收，例如wi-fi设备)；
图4是根据本发明的一个实施例的确定菲涅尔相位偏移的方法的流程框图。
图5是根据本发明的一个实施例的利用图4所示的确定菲涅尔相位偏移的方法进行定位的方法的流程框图。
图6是根据本发明的一个实施例的利用图4所示的确定菲涅尔相位偏移的方法进行定位的系统的框图。
图7以方框图的形式描述了一种示例数据处理系统。
具体实施方式
下面结合附图，通过实施例进一步描述本发明，本领域技术人员能够理解的是，该具体描述目的不在于以任何方式限制本发明的范围。
图3是根据本发明的一个实施例的确定菲涅尔相位偏移的系统的框图，其中，(a)为发送设备(其支持包括wi-fi信号的射频信号发射，例如wi-fi设备，)；(b)为接收设备(其支持包括wi-fi信号的射频信号接收，例如wi-fi设备)。应当理解的是，在一个实施例中，发送设备(a)可以不是确定菲涅尔相位偏移的系统的一部分。例如，发送设备可以是需要确定菲涅尔相位偏移的场所中所具有能够进行射频信号发送(例如wi-fi信号发送的设备)，其能够发送信号以协助确定菲涅尔相位偏移的系统确定菲涅尔相位偏移，但其本身不是确定菲涅尔相位偏移的系统的一部分。在一个实施例中，确定菲涅尔相位偏移的系统不包括发送设备(a)，但是包括经配置能够执行本文描述的确定菲涅尔相位偏移的方法的其他设备。在不包括发送设备(a)的实施例中，确定菲涅尔相位偏移的系统能够包括接收设备(b)或者经配置能够执行本文描述的确定菲涅尔相位偏移的方法的其他设备。在一个实施例中，确定菲涅尔相位偏移的系统不包括发送设备(b)，也不包括经配置能够执行本文描述的定位方法的其他设备，接收设备(b)经配置执行本文描述的确定菲涅尔相位偏移的方法。通过对本发明各个实施例的描述，本领域技术人员能够理解上述各种确定菲涅尔相位偏移的系统的变形和修改，其均落入本发明的保护范围内。
在图3中的确定菲涅尔相位偏移的系统以及本发明的其他实施例的确定菲涅尔相位偏移的系统中，取决于系统实时性、实现成本和测量精度等，收发设备的天线可适当增加或减少。例如，在一个实施例中，发送设备或接收设备的天线可选择为3根。在另一个实施例中，发送设备或接收设备的天线可选择为2根。应当理解，为了实施本发明中的定位方法，在发送设备或接收设备上需要分别安装至少一根天线，分别用于发送信号或接收信号。
在一个实施例中，确定菲涅尔相位偏移的系统中的发送设备和接收设备的天线，例如，可以是垂直极化的全向天线且垂直地面放置。
发送设备和接收设备，例如，可以是笔记本电脑、minipc、路由器以及任何支持rf(射频)信号收发的设备。
图4是根据本发明的一个实施例的确定菲涅尔相位偏移方法的流程框图，其适于使用图3中以及本发明中其他实施例中的确定菲涅尔相位偏移的系统进行确定菲涅尔相位偏移的方法。
图4中的方法开始于步骤410。需要注意的是，图4仅以示例的方式描述本发明的一个实施例，其中，用于执行图4中所示方法需要至少两台收发设备，一台发射设备(t)能提供射频信号例如wi-fi信号以提供信道状态信息，一台接收设备(r)能接收射频信号例如wi-fi信号以测量信道状态信息；移动的反射目标不携带任何设备。
在步骤410处，接收设备r从发送设备t接收射频信号，所述射频信号包括不同波长的至少两个载波。
在步骤420处，以发送设备t的位置、第一接收设备r1的位置为椭圆焦点，构建所述至少两个载波的k层菲涅尔区，其中，构建菲涅尔区的步骤具体如下：
以发送设备t的位置和接收设备r1的位置为椭圆焦点，构建k层菲涅尔区的步骤如下：确定波长为λa的第a个载波的k层菲涅尔区；令qi是第i层同心椭圆的轨迹，则qi的轨迹满足：
|tqi|+|qir1|-|tr1|＝iλa/2(式1)
其中，|tqi|是qi到发射设备t的距离，|qir1|是qi到接收设备r1的距离，|tr1|为收发设备之间的距离，λa为第a个载波的波长，m为菲涅尔区的最大层数，由内层到外层分别为第1菲涅尔区到第m菲涅尔区；取决于关注的范围，m越大，该层菲涅尔区距离收发设备越远。
例如，参照图2，获得通过目标反射的信号的总路径长度，相当于获得椭圆菲涅尔区的2an，由于已知发送设备和接收设备之间的距离d0，其相当于2c，因此，由an2-c2＝bn2可以获得目标所在的菲涅耳区的椭圆方程
定位目标穿过菲涅尔区时产生的信号为
其中，
其中，d0是los信号的长度，其中los为视距信号或者直接路径信号，即从发送设备t到接收设备的直接r1路径长度，
构建针对以发送设备t和接收设备r1为焦点的一维多载波k层菲涅尔区。例如，wi-fi设备可以提供多个载波。例如，5300网卡可以提供30个子载波，每个子载波之间的波长有轻微的差距，那么不同载波之间的菲涅尔相位具有一定的差值，即，菲涅尔相位差为不同载波之间菲涅尔相位的差值，表示为
其中，δf为两个载波λa和λb之间的频率差，c为光速，δεab为由室内静态路径引入的两个载波之间的菲涅尔相位偏移，不考虑设备因素，该值只与选择的载波和环境有关。为便于描述，当δεab为零时，将公式4中的
在步骤430处，在所述环境中布置反射目标，使得反射目标在所述环境中移动。在一个优选实施例中，所述反射目标沿着垂直于发送设备t和接收设备r的方向移动。在一个实施例中，所述反射目标为强反射目标，其使得到达其自身的射频信号的大部分被反射出去。在一个实施例中，所述反射目标为金属。在一个优选实施例中，所述反射目标为铁板。在一个实施例中，所述反射目标为人体。
在步骤440处，针对所述至少两个载波中的至少一个载波对组合，确定所述反射目标引起的不同载波之间的理论菲涅尔相位差。例如，根据发送设备t和接收设备r之间的距离、移动的反射目标的坐标位置，以及不同载波的频率或波长，确定所述反射目标引起的不同载波之间的理论菲涅尔相位差。
在步骤450处，针对所述至少两个载波中的所述至少一个载波对组合，确定所述反射目标和环境引起的不同载波之间的实际菲涅尔相位差。例如，通过接收设备r1，确定所述反射目标和环境引起的不同载波之间的实际菲涅尔相位差。
在步骤460处，将实际菲涅尔相位差与相应的理论菲涅尔相位差进行比较，确定在所述环境中的所述至少一个载波对组合的菲涅尔相位偏移。在一个优选的实施例中，确定在所述环境中的菲涅尔相位偏移的步骤包括：针对一段时间内接收到的信号，确定任意一对载波在所有固定时间窗口内的菲涅尔相位偏移，取所有固定时间窗口的菲涅尔相位偏移的平均值或中位数作为所述任意一对载波对之间的菲涅尔相位。由于菲涅尔相位偏移只跟两个载波有关系，和移动物体所在的位置无关，而且是一个固定值，所以需要针对每一对载波求得一个菲涅尔相位偏移，以便在定位时对每对载波的菲涅尔相位进行修正。在一个实施例中，接收设备根据接收到的信号，将其划分为多个固定大小的时间窗口，例如0.05s的时间窗口，针对每对载波在所有时间窗口内计算菲涅尔相位偏移，并将计算得到的菲涅尔相位偏移的平均值或中位数作为该对载波的菲涅尔相位偏移。
图5是根据本发明的一个实施例的利用图4所示的确定菲涅尔相位偏移的方法进行定位的方法的流程框图。图5中的方法开始于步骤510。需要注意的是，图5仅以示例的方式描述本发明的一个实施例，其中，用于执行图4中所示方法需要至少三台收发设备，一台发射设备(t1)能提供射频信号例如wi-fi信号以提供信道状态信息，两台接收设备(r1和r2)能接收射频信号例如wi-fi信号以测量信道状态信息；待定位的移动目标不携带任何设备。通过利用rf(射频)信号对定位目标实现无接触感知定位，获得目标的精确位置。
需要指出的是，相同的移动目标在相同时间位于定位环境中时，该移动目标引起的发送设备所发送的信号中的扰动在同时传达到不同接收设备时，接收设备r1和接收设备r2使用前面累积收到的一段时间信号形成的时间窗口对信号进行处理，以迭代滑动的方式在固定的时间窗口逐步进行。比如有5秒钟的信号，时间窗口是0.05秒，那么整个信号可以划分为100个时间窗口，使用每个时间窗口内的信号进行一次定位，也就是进行了100次定位；迭代的意思是每次定位为了使用更长时间的信号，例如，使用两个窗口定位，每次向后移动例如一个时间窗口。例如，先使用1、2时间窗口内的信号进行一次定位，下一次使用2、3时间窗口内的信号进行定位。本领域技术人员能够理解的是，为了使用更长时间的信号进行定位，可以选择其他长度的时间窗口进行定位，并且每次迭代可以移动其他数量的时间窗口。
在步骤510处，第一接收设备r1和第二接收设备r2同时从发送设备t接收射频信号，所述射频信号包含两个或更多个载波。本领域技术人员能够理解，第一接收设备r1和第二接收设备r2同时从发送设备t接收射频信号指的是，相同的移动目标在相同时间位于定位环境中时，该移动目标引起的发送设备所发送的信号中的扰动同时传达到不同接收设备r1和r2。接收设备r1和接收设备r2根据接收的累计信号对相同的固定的时间窗口内的信号进行处理，从而定位移动目标，如上文所述。
在步骤520处，针对所述两个或更多个载波中的至少一个载波，以发送设备t的位置、第一接收设备r1的位置为椭圆焦点，构建第一接收设备r1的m层菲涅尔区；在步骤530处，针对所述两个或更多个载波中的至少一个载波，以发送设备t的位置、第一接收设备r2的位置为椭圆焦点，构建n层菲涅尔区。其中，构建菲涅尔区的步骤具体如下：
步骤5200与步骤530中的方法类似，因此为简洁描述，本文以发送设备t的位置和接收设备r1的位置为例，具体描述构建第一接收设备r1的m层菲涅尔区的步骤。本领域技术人员参考该描述，将能够理解并实施针对所述两个或更多个载波中的至少一个载波，以发送设备t的位置和接收设备r2的位置，构建第二接收设备r2的n层菲涅尔区的步骤。
针对所述两个或更多个载波中的至少一个载波，以发送设备t的位置和第一接收设备r1的位置为椭圆焦点，构建第一接收设备r1的m层菲涅尔区的步骤如下：确定波长为λa的第a个载波的m层菲涅尔区；令qi是第i层同心椭圆的轨迹，则qi的轨迹满足：
|tqi|+|qir1|-|tr1|＝iλa/2(式1)
其中，|tqi|是qi到发射设备t的距离，|qir1|是qi到接收设备r1的距离，|tr1|为收发设备之间的距离，λa为第a个载波的波长，m为菲涅尔区的最大层数，由内层到外层分别为第1菲涅尔区到第m菲涅尔区；取决于关注的范围，m越大，该层菲涅尔区距离收发设备越远。
例如，参照图2，获得通过目标反射的信号的总路径长度，相当于获得椭圆菲涅尔区的2an，由于已知发送设备和接收设备之间的距离d0，其相当于2c，因此，由an2-c2＝bn2可以获得目标所在的菲涅耳区的椭圆方程
m可以根据实际情况设置取值。例如，为消除二义性，m为大于零的自然数且取值满足：在m层菲涅尔区中，任意两个载波之间的菲涅尔相位差小于2π。类似地，针对所述两个或更多个载波中的至少一个载波，在以发送设备t的位置和接收设备r2的位置，构建n层菲涅尔区时，n为大于零的自然数且取值满足：在n层菲涅尔区中，任意两个载波之间的菲涅尔相位差小于2π。
定位目标穿过菲涅尔区时产生的信号为
其中，
其中，d0是los信号的长度，其中los为视距信号或者直接路径信号，即从发送设备t到接收设备r1的直接路径长度，
构建针对以发送设备t和接收设备r1为焦点的一维多载波m层菲涅尔区。wi-fi设备通常可以提供多个载波。例如，5300网卡可以提供30个子载波，每个子载波之间的波长有轻微的差距，那么不同载波之间的菲涅尔相位具有一定的差值，即，两个载波之间菲涅尔相位的差值被称为菲涅尔相位差，表示为
其中，δf为两个载波λa和λb之间的频率差，c为光速，δεab为由室内静态路径引入的两个载波之间的菲涅尔相位偏移，不考虑设备因素，该值只与选择的载波和环境有关。
在步骤540中，确定针对以发送设备t和接收设备r1为焦点的一维方向上目标在m层菲涅尔区中所在的菲涅尔区f1。由式4可以知道，只要获得了两个载波之间的实际菲涅尔相位差和菲涅尔相位偏移，就能获得由目标反射的信号的总路径长度，也就获得目标所在的菲涅尔区。例如，参照图2，获得通过目标反射的信号的总路径长度，相当于获得椭圆菲涅尔区的2an，由于已知发送设备和接收设备之间的距离d0，其相当于2c，因此，由an2-c2＝bn2能够获得目标所在的菲涅尔区的椭圆方程
利用wi-fi设备，根据csi信息，获得目标所在的菲涅尔区的步骤包括：根据接收信号在滑动窗口内的周期数和时间差估计实际菲涅尔相位差，其中周期数计算方法包括：1、使用快速傅里叶变换转换到频域，其主频率为所要周期数；2、寻找信号的波峰处，根据其所占时间的大小估计周期数，具体步骤如下：
对csi信号进行平滑和标准化处理，使得信号具有可比性；
在一个时间窗口内计算两个载波之间的时延，设定时间窗口，使用互相关的方法求得相关性最大的时延为两个载波之间的时延；
估计时间窗口内的周期数，联合使用fft频域方法与数波峰方法计算精确的周期数；
根据时延和周期数计算两个载波之间的实际菲涅尔相位差；
菲涅尔相位偏移的获得方法如在图4中所述的确定菲涅尔相位偏移的方法。其中，根据图4中所述的方法，针对每组两个载波的多组载波计算菲涅尔相位偏移。该多个菲涅尔相位偏移可以存储在数据库或数据表中，供后续进行定位时使用。计算菲涅尔相位偏移的方法可以包括：固定时间窗口，在每个时间窗口内根据接收到的信号计算任意两个载波之间的实际菲涅尔相位差，根据目标所在位置计算理论上的菲涅尔相位差，根据实际菲涅尔相位差和理论菲涅尔相位差的差值，从而可得任意两个载波之间的菲涅尔相位偏移。
根据公式4可知，只要获得了两个载波之间的实际菲涅尔相位差和菲涅尔相位偏移，就能获得目标经过反射的信号的路径长度，也就获得目标所在的菲涅尔区。根据csi信息获得目标所在的菲涅尔区包括：根据菲涅尔相位偏移修正实际菲涅尔相位差，从而获得理论菲涅尔相位差，由理论菲涅尔相位差直接通过几何映射的方法获得目标所在的菲涅尔区。
由实际菲涅尔相位差和菲涅尔相位偏移获得理论菲涅尔相位差，进而获得目标所在的菲涅尔区，例如，直接通过几何映射的方法获得目标所在的菲涅尔区，其中几何映射是指通过查询数据库或数据表获得目标所在的菲涅尔区，其中数据库或数据表中预先存储有不同菲涅尔相位差对应的菲尔区。
在步骤550中，确定针对以发送设备t和接收设备r2为焦点的一维方向上目标在n层菲涅尔区中所在的菲涅尔区f2。该步骤类似于步骤540中的步骤，为简洁描述，不再赘述。本领域技术人员通过对步骤530的描述，能够理解并实施步骤550中的步骤。
在步骤560中，确定菲涅尔区f1和f2的交点。在步骤570中，根据菲涅尔区f1和f2的交点，确定目标的位置。也就是，确定目标在二维(即发送设备t、接收设备r1和接收设备r2所在的平面)上的位置坐标，在获得每一维方向上目标所在的菲涅尔区的基础上，使用解析或者映射的方法获得目标在二维上的位置坐标，其中解析方法是指通过计算两种菲涅尔区的椭圆方程的交点坐标获得位置坐标，映射方法是指通过查询数据库或数据表获得目标在二维方向上的位置坐标，其中数据库或数据表中预先存储有针对在每一维上的不同菲涅尔区对应的二维方向上的位置坐标。
需要指出的是，上述各种特征和步骤可以彼此独立使用，或可以以各种方式结合。所有可能的组合和子组合旨在落入本申请的范围内。此外，在一些实施方式中可以省略某些步骤或过程块。本文描述的方法和步骤也不限于任何特定顺序，并且能够以合适的其他顺序执行。例如，描述的步骤可以以不同于具体公开的顺序的顺序执行，或多个步骤可以在单个块中组合。示例步骤或过程块可以串行、并行或以某种其他方式执行。步骤或过程块可以被添加到所公开的示例实施例中或从所公开的示例实施例中移除。本文描述的示例系统和部件可以配置为不同于所描述的结构。例如，设备可以被添加到所公开的示例实施例、从所公开的示例实施例中移除或不同于所公开的示例实施例被重新布置。
为了更清楚地传达本发明的技术方案，下面更具体描述了一个实施例。以下实施例使用中心频率为5.745ghz的wifi信号作为rf信号；使用搭载8db全向天线、intel5300网卡的minipc作为rf收发设备，csi中包含30个载波，带宽为40mhz，则载波间隔为1.25mhz。
以6m×6m的房间为例，根据本发明提供的室内无接触感知定位方法，执行如下步骤：
将rf收发设备(minipc)上的收发天线垂直地面放置，rf发送设备放置在离任一墙角的对应墙壁0.5m的距离，rf接收设备分别距离发送设备5m，且互相垂直，发送设备的位置标记为t，接收设备的位置分别标记为r1，r2；
定位目标穿过菲涅尔区时产生的信号为
其中，
其中，d0是los信号的长度，当d0和λ固定时，菲涅尔相位随着反射路径而变化。本实施例中，d0＝5m，
获得室内多载波菲涅尔相位差，wifi设备可以提供多个载波，intel5300网卡提供了30个子载波，那么两个载波之间的菲涅尔相位差表示为
其中，δf为两个载波λa和λb之间的频率差，c为光速。本实施例中，δf＝(a-b)·1.25mhz，c＝3×108m/s，并且动态选择多组载波，以消除异常点的影响，使得结果更准确。δεab为由室内静态路径引入的两个载波之间的菲涅尔相位偏移，不考虑设备因素，该值只与选择的载波和环境有关。
由csi获得多载波菲涅尔相位差的步骤包括：
对csi信号进行平滑和标准化处理，使得信号具有可比性；本实施例中使用savitzky-golay降噪技术进行了10次平滑处理，并使用zscore方法进行标准化处理；
在一个时间窗口内计算两个载波之间的时延，设定时间窗口，使用互相关的方法求得相关性最大的时延为两个载波之间的时延；本实施例中设置信号的采样频率为500，使用0.05秒的时间窗口，即25个采样点，并重叠一个时间窗口计算两个载波之间的时延；
估计时间窗口内的周期数，联合使用fft频域方法与数波峰方法计算精确的周期数；本实施例中先用fft方法计算准确的周期数，然后使用数波峰方法计算更加精确的周期数，以两种方法中最优的周期数为该时间窗口内的周期数，即以fft方法计算的周期数为标准，若使用数波峰的方式计算的周期与fft方法的周期之间的误差在0.5的范围内，以数波峰的方式计算的周期为准，否则以fft方法的周期为准。例如，fft方法计算为2个周期，数波峰方法计算为2.7个周期，则最终结果为2个周期；若fft方法计算为2个周期，数波峰方法计算为2.1个周期，则最终结果为2.1个周期；
根据时延和周期数计算两个载波之间的菲涅尔相位差；本实施例中，两个载波之间的菲涅尔相位差为
计算室内菲涅尔相位偏移，因为在不考虑设备因素的情况下，两个载波之间的菲涅尔相位偏移只与选择的载波和环境有关，所以在室内为了获得准确的菲涅尔相位差，需要先计算出两个载波之间的菲涅尔相位偏移，根据两个载波之间获得的理论菲涅尔相位差和实际菲涅尔相位差获得菲涅尔相位偏移；本实例中在部署好设备后，首先使用铁板沿着每一维(发送设备t和接收设备r1的方向；以及发送设备t和接收设备r2的方向)los的中垂线(垂直于发送设备t和接收设备r1的方向；和垂直于发送设备t和接收设备r2的方向)移动，根据获得到的菲涅尔相位差和菲涅尔相位差的groundtruth计算得打每对载波之间的菲涅尔相位偏移；
确定一维方向上目标所在菲涅尔区，由式4可以知道，只要获得了两个载波之间的实际菲涅尔相位差和菲涅尔相位偏移，就能获得由目标反射的信号的路径长度，也就获得目标所在的菲涅尔区，根据csi信息获得目标所在的菲涅尔区包括
根据获得菲涅尔相位偏移对菲涅尔相位差进行修正
由修正后的菲涅尔相位差直接通过几何映射的方法获得目标所在的菲涅尔区；
确定目标在二维上的位置坐标，在获得每一维方向上目标所在的菲涅尔区的基础上，使用解析或者映射的方法获得目标在二维上的位置坐标，即求两个椭圆的交点为目标的实际位置，限定两维坐标都大于0，即可唯一确定一个交点，也就是目标的位置坐标。
图6是根据本发明的一个实施例的利用图4所示的确定菲涅尔相位偏移的方法进行定位的系统的框图，其可以与图5中的定位方法配合使用。其中，(a)为发送设备(其支持包括wi-fi信号的射频信号发射，例如wi-fi设备，)；(b)和(c)为接收设备(其支持包括wi-fi信号的射频信号接收，例如wi-fi设备)。应当理解的是，在一个实施例中，发送设备(a)可以不是定位系统的一部分。例如，发送设备可以是需要进行定位的场所中所具有能够进行射频信号发送(例如wi-fi信号发送的设备)，其能够发送信号以协助定位系统进行定位，但其本身不是定位系统的一部分。在一个实施例中，定位系统不包括发送设备(a)，但是包括经配置能够执行本文描述的定位方法的其他设备。在不包括发送设备(a)的实施例中，定位系统能够包括接收设备(b)或(c)或者经配置能够执行本文描述的定位方法的其他设备。在一个实施例中，定位系统不包括发送设备(a)，也不包括经配置能够执行本文描述的定位方法的其他设备，接收设备(b)和接收设备(c)中的一个经配置执行本文描述的定位方法。在一个实施例中，定位系统不包括发送设备(a)，也不包括经配置能够执行本文描述的定位方法的其他设备，接收设备(b)和接收设备(c)两者经配置协作执行本文描述的定位方法。通过对本发明各个实施例的描述，本领域技术人员能够理解上述各种定位系统的变形和修改，其均落入本发明的保护范围内。
在图6中的定位系统以及本发明的其他实施例的定位系统中，取决于系统实时性、实现成本和测量精度等，收发设备的天线可适当增加或减少。例如，在一个实施例中，发送设备或接收设备的天线可选择为3根。在另一个实施例中，发送设备或接收设备的天线可选择为2根。应当理解，为了实施本发明中的定位方法，在发送设备或接收设备上需要分别安装至少一根天线，分别用于发送信号或接收信号。
在一个实施例中，在定位系统中的发送设备和接收设备的天线，例如，可以是垂直极化的全向天线且垂直地面放置。在一个实施例中，从发送设备t到接收设备r1的向量方向，例如，可以与从发送设备t到接收设备r2的向量方向成直角。在其他实施例中，从发送设备t到接收设备r1的向量方向，例如，可以与从发送设备t到接收设备r2的向量方向成锐角或钝角。
在一个实施例中，发送设备t和接收设备r1和r2，例如，可以是笔记本电脑、minipc、路由器以及任何支持rf(射频)信号收发的设备。
图7以方框图的形式描述了数据处理系统。数据处理系统700可用于实施本发明中的各种系统。
如图所示，数据处理系统700包括通信架构702，其提供处理器单元704、存储装置706、通信单元708、输入/输出单元710和显示器之间的通信。在一些情形中，通信架构702可以被实施为总线系统。
处理器单元704被配置为执行软件的指令从而执行若干操作。处理器单元704根据实施可以包括若干处理器、多核处理器和/或一些其他类型的处理器。在一些情形中，处理器单元704可以采用硬件单元的形式，例如，电路系统、专用集成电路(asic)、可编程逻辑器件或一些其他合适类型的硬件单元。
处理器单元704运行的操作系统的指令、应用程序和/或程序可位于存储装置706中。存储装置706可以通过通信架构702与处理器单元704通信。如本文所用的，存储装置也被称为计算机可读存储装置，是任何能够将信息暂存和/或永久性存储的任何一件硬件。该信息可以包括，但不限于，数据、程序代码和/或其他信息。
存储器714和持久存储器716是存储装置706的实例。存储器714可以采用，例如随机存取存储器或一些易失性或非易失性存储装置的形式。持久存储器716可以包括任何数量的组件和装置。例如，持久存储器716可以包括硬件设备、闪存、可重写光盘、可重写磁带或上述其他组合。持久存储器716使用的媒介可以是或可以不是可拆卸的。
通信单元708允许数据处理系统700与其他数据处理系统和/或设备通信。通信单元708可以使用物理和/或无线通信链路提供通信。
输入/输出单元710允许从与数据处理系统700相连的其他设备处接收输入和发送输出。例如，输入/输出单元710可以允许用户通过键盘、鼠标和/或一些其他类型输入设备接收输入。作为另一个例子，输入/输出单元710可以允许发送输出至与数据处理系统700相连的打印机。
显示器712被配置为向用户显示信息。显示器712可以包括，例如，但不限于，监视器、触摸屏、激光显示器、全息显示器、虚拟显示设备和/或一些其他类型显示设备。
在这个示例性实例中，通过处理器单元704使用计算机实施指令可以执行不同示例性实施例的过程。这些指令可以被称为程序代码、计算机可用程序指令或计算机可读程序代码，而且处理器单元704中的一个或更多处理器可以读取和执行。
在这些实例中，程序代码718以功能性形式位于计算机可读介质720中，而且可以被加载到或被传输到数据处理系统700以被处理器单元执行，其中计算机可读介质720可以被选择性拆卸。程序代码718和计算机可读介质720共同形成计算机程序产品722。在这个示例性实例中，计算机可读介质720可以是计算机可读存储介质724或计算机可读信号介质726。
计算机可读存储介质724是用于存储程序代码718的物理的或有形存储设备而不是传播或传输程序代码718的媒介。计算机可读存储设备724可以是，例如，但不限于，与数据处理系统相连的光或磁盘或持续存储设备。
可替换地，使用计算机可读信号介质可以将程序代码718传输到数据处理系统700中。计算机可读信号介质726可以是，例如，包含程序代码718的传播的数据信号。该数据信号可以是电磁信号、光信号和/或一些其他能够通过物理和/或无线通信链路传输的信号。
图7所示的数据处理系统700并不意图提供可以实施这些示例性实施例的方式的结构限制。可以实施不同示例性实施例的数据处理系统包括附加组件或替代数据处理系统700所示的那些组件的组件。此外，图7中所示组件可以与所示实例中的不同。
本领域普通技术人员将理解，结合本文中所公开的示例所描述的各种说明性逻辑块、模块及过程可实施为电子硬件、计算机软件或两者的组合。此外，前述过程能够体现为计算机可读介质，其使处理器或计算机进行或执行某些功能。
为了清楚地说明硬件和软件的这种可交换性，上文通常就它们的功能对各种说明性的部件、框和模块进行了描述。该功能是被实现为硬件还是软件取决于特定的应用和施加于整个系统上的设计约束。本领域普通技术人员可以对每个特定的应用以不同方法来实现所描述的功能，但是这种实现决策不应解释为超出本公开的范围。
结合这里所公开的示例所描述的各种说明性的逻辑块、单元、模块和控制器可以利用通用处理器、数字信号处理器(dsp)、专用集成电路(asic)、现场可编程门阵列(fpga)或其他可编程逻辑器件、分立门或晶体管逻辑、分立硬件部件、或者设计为实现这里所描述的功能的任何组合来实现或者执行。通用处理器可以是微处理器，但是在替代方案中，处理器可以是任何常规处理器、控制器、微控制器或状态机。处理器也可以实现为计算装置的组合，例如，dsp与微处理器的组合、多个微处理器、一个或多个微处理器结合dsp核或者任何其他此类配置。
结合本文中公开的示例所描述的方法或过程的行为可直接在硬件中、在由处理器执行的软件模块中、或在这两者的组合中体现。方法或算法的步骤还可以与在示例中所提供的那些不同的替代顺序来执行。软件模块可以存在于ram存储器、闪速存储器、rom存储器、eprom存储器、eeprom存储器、寄存器、硬盘、可移除介质、光学介质或者本领域已知的任何其他形式的存储介质中。示例性存储介质耦合到处理器，使得该处理器能够从/向该存储介质读写信息。在替代方案中，存储介质可以与处理器集成。处理器和存储介质可以驻留在专用集成电路(asic)中。
需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。
- 还没有人留言评论。精彩留言会获得点赞！