一种细粒度免校验的室内被动人体检测方法

1.一种细粒度免校验的室内被动人体检测方法，其特征在于，包括以下几个步骤：
步骤一：构建标准信号特征文件；
步骤二：根据标准信号特征文件计算接收机敏感度因子；
步骤三：根据接收机敏感度因子估计入侵检测阈值；
步骤四：接收机实时收集当前检测信号的信道状态信息，并从中提取检测信号特征，计算检测信号特征与标准信号特征相关性；
步骤五：判断检测信号特征与标准信号特征相关性是否小于检测阈值，如果小于阈值，则有入侵人体出现，否则无人出现；
所述的构建标准信号特征文件Hstd的方法为：在初始时刻，无人情况下，选取m个数据包，提取每一个数据包中的信号特征，形成标准信号特征文件Hstd，其中信号特征为从信道状态信息中提取频域上所有子载波的振幅向量；
所述的信号特征求取过程为：
建立无线通信模型：

其中，t是时间，y是接收信号，x是发送信号，h是信道状态信息，z为高斯白噪音，计算每个子载波信道状态信息的振幅向量作为信号特征，
在频率域中，OFDM系统在OFDM子载波粒度上信道频率响应为：
H＝{H(1),H(2),...,H(N)}
每个子载波的CFR信息为：
H(f)＝|H(f)|exp(jsin(∠H(f)))
其中H(f)表示子载波的振幅响应，而∠H(f)表示子载波的相位响应；
所述的计算接收机敏感度因子的方法为，
1)将每个子载波的振幅归一化到中心频率上：

其中H(fk)是第k个子载波原始振幅值，Hnorm(fk)第k个子载波归一化后的振幅值，f0是中心频率值，fk是第k个子载波的频率值，
2)求取归一化后的振幅向量的变异系数：

其中，std(Hnorm)和mean(Hnorm)分别是归一化CSI振幅Hnorm的标准差与均值，
3)求取接收机敏感度因子变量Ks：

其中x为变异系数ν的测量值，μ为变异系数ν的均值，σ为变异系数ν的标准差值。
2.根据权利要求1所述的一种细粒度免校验的室内被动人体检测方法，其特征在于：所述的根据接收机敏感度因子估计入侵检测阈值C为：
C＝a·Ks+b
其中，a和b为系数。
3.根据权利要求1所述的一种细粒度免校验的室内被动人体检测方法，其特征在于：所述的检测信号为入侵检测阈值估计完成后，接收机实时采集到的信号，从检测信号的信道状态信息中提取的信号特征为检测信号特征Hdetect，检测信号特征Hdetect与标准信号特征Hstd相关性为：

一种细粒度免校验的室内被动人体检测方法\n技术领域\n[0001] 本发明属于一种无线定位技术，尤其能够自动预测检测阈值的，一种细粒度免校验的室内被动人体检测方法。\n背景技术\n[0002] 随着科学技术的不断发展，室内入侵技术也在不断的提升。基于压力感知、红外感知、图像感知的传统入侵检测逐渐被破解。无线入侵检测系统继而得到了极大的关注。无线入侵检测系统通过分析接收到的无线信号特征来捕获环境变化，及时发现入侵行为。无线信号所具有全向覆盖、环境变化感知敏感的特点帮助其在入侵检测领域发挥出优异的表现。近几年，室内WLAN技术的快速发展给室内入侵检测带来了新的发展浪潮。价格便宜的，易于部署的WLAN设备目前已被广泛部署在工厂、楼宇、家庭环境中。这些WLAN设备不仅仅是数据传输的路由设备，更是廉价的无线信号发射机。若他们被应用到室内无线入侵检测系统中，将使室内入侵检测系统具有更强的普适性，能够应用到不同场景下。因此，许多学者在基于WLAN设备的室内无线入侵检测系统中付诸了大量精力进行研究。\n[0003] RSSI作为一个来自链路层的信号特征，其易于获取、易于处理的特点使其首先作为无线信号特征被应用到室内无线入侵检测系统中。其中，室内无线入侵检测，也称为室内被动入侵检测，被认为是入侵者不携带任何无线设备，系统通过分析监测区域内链路RSSI值变化来实现入侵检测。为了增大监测区域和提高检测精度，基于RSSI的室内入侵检测系统采用多链路架构，然而这种架构增加了系统花费和部署开销。另一方面，在室内环境下由于信号多径效应的影响，RSSI值具有时间上的不稳定性。这个缺点限制了室内入侵检测系统精度的提高。因此，许多学者探索利用细粒度的物理层信息来克服多径效应，以提高室内入侵检测系统的精度。\n[0004] 物理层信息中最常用的特征是信道响应。室内信号在传播过程中遭受了不同的能量衰减、信号反射、散射，以分簇形式先后达到接收机。为了从时域上区分多径信号，无线信道被建模为一个时间线性过滤器，即信道时域响应。然而我们发现在室内环境下，信道时域响应对于入侵行为的敏感度较低。值得注意的是，正交频分复用系统提供了一个简单的信道频域响应，将信道划分为多个不同频率的子载波，每个子载波具有独立的振幅与相位。室内环境下，信道频域响应具有频域选择性衰减。不同频率的子载波在经历相似的穿透、反射后，其振幅将发生不同的变化。因而其对入侵行为具有较高的敏感度。幸运地是，信道频域响应能够从普通商业网卡中提取并以信道状态信息的形式提交给用户态。大量的基于信道状态信息的细粒度室内入侵检测系统继而问世,并表现出优异的性能。\n[0005] 在值得庆幸的同时，我们注意到目前已有的细粒度室内无线入侵检测系统，包括两个部分：离线训练和在线检测。离线训练过程需要采集无人情况下和有入侵者情况下的信号特征信息来确定检测阈值。这个过程也是无法避免的，因为在现实环境中，由于环境结构不同，信道频域响应不同。当入侵行为发生时，入侵者对多径信号产生的干扰导致的频域选择性衰减不同，因此需要现场勘测来确定检测阈值。但是离线训练是一个耗时的、劳动密集型的过程，并且需要专业人士进行现场勘测和参数调控，这无疑增加了系统开销和部署难度。为此，我们在进行了大量的实验和深入的分析后发现，检测阈值大小与室内多径信号变化表现出的敏感度相关。在保证不丢失大量的检测精度的前提下，通过探测室内多径信号成分的多少可以预测检测阈值大小。估算多径信号成分与视距路径比重常用的方法是计算莱斯因子。然而在室内环境下，当视距路径信号能量小时，莱斯因子对多径变化表现较低的敏感度。同时分析室内信号多径传播往往需要利用专业设备，如VNA，SDN等，这些设备价格昂贵，携带不便，操作繁琐。\n发明内容\n[0006] 本发明的目的是提供部署快捷方便并且能够预测检测阈值的，细粒度免校验的室内被动人体检测方法。\n[0007] 本发明是通过以下技术方案实现的：\n[0008] 一种细粒度免校验的室内被动人体检测方法，包括以下几个步骤：\n[0009] 步骤一：构建标准信号特征文件；\n[0010] 步骤二：根据标准信号特征文件计算接收机敏感度因子；\n[0011] 步骤三：根据接收机敏感度因子估计入侵检测阈值；\n[0012] 步骤四：接收机实时收集当前检测信号的信道状态信息，并从中提取检测信号特征，计算检测信号特征与标准信号特征相关性；\n[0013] 步骤五：判断检测信号特征与标准信号特征相关性是否小于检测阈值，如果小于阈值，则有入侵人体出现，否则无人出现。\n[0014] 本发明一种细粒度免校验的室内被动人体检测方法，还可以包括：\n[0015] 1、构建标准信号特征文件Hstd的方法为：在初始时刻，无人情况下，选取m个数据包，提取每一个数据包中的信号特征，形成标准信号特征文件Hstd，其中信号特征为从信道状态信息中提取频域上所有子载波的振幅向量。\n[0016] 2、信号特征求取过程为：\n[0017] 建立无线通信模型：\n[0018]\n[0019] 其中，t是时间，y是接收信号，x是发送信号，h是信道状态信息，z为高斯白噪音，[0020] 计算每个子载波信道状态信息的振幅向量作为信号特征，\n[0021] 在频率域中，OFDM系统在OFDM子载波粒度上信道频率响应为：\n[0022] H＝{H(1),H(2),...,H(N)}\n[0023] 每个子载波的CFR信息为：\n[0024] H(f)＝|H(f)|exp(jsin(∠H(f)))\n[0025] 其中H(f)表示子载波的振幅响应，而∠H(f)表示子载波的相位响应。\n[0026] 3、计算接收机敏感度因子的方法为，\n[0027] 1)将每个子载波的振幅归一化到中心频率上：\n[0028]\n[0029] 其中H(fk)是第k个子载波原始振幅值，Hnorm(fk)第k个子载波归一化后的振幅值，f0是中心频率值，fk是第k个子载波的频率值，\n[0030] 2)求取归一化后的振幅向量的变异系数：\n[0031]\n[0032] 其中，std(Hnorm)和mean(Hnorm)分别是归一化CSI振幅Hnorm的标准差与均值，[0033] 3)求取接收机敏感度因子变量Ks：\n[0034]\n[0035] 其中x为变异系数ν的测量值，μ为变异系数ν的均值，σ为变异系数ν的标准差值。\n[0036] 4、根据接收机敏感度因子估计入侵检测阈值C为：\n[0037] C＝a·Ks+b\n[0038] 其中，a和b为系数。\n[0039] 5、检测信号为入侵检测阈值估计完成后，接收机实时采集到的信号，从检测信号的信道状态信息中提取的信号特征为检测信号特征Hdetect，检测信号特征Hdetect与标准信号特征Hstd相关性为：\n[0040]\n[0041] 有益效果：\n[0042] 本发明提出一个细粒度免校验的室内被动入侵检测的方法，此方法包括四个步骤。通过现场测量接收机敏感度因子大小，利用阈值预测模型预测检测阈值，最后入侵检测系统根据阈值大小实现无线入侵发现。该方法依赖较少的人工劳动，部署快捷方便，大大地降低了入侵检测系统的部署难度。\n[0043] 本发明仅依赖一组接收机与发射机，能够适应不同大小空间下的人体检测需求。\n同时该发明也可以扩展到多组接收机与发射机来适应较大空间下的人体入侵检测，具有灵活的扩展性。\n附图说明\n[0044] 图1是本发明的细粒度免校验的室内被入侵检测的流程图。\n[0045] 图2是本发明的细粒度免校验的室内被动入侵检测的示意图。\n具体实施方式\n[0046] 下面将结合附图对本发明做进一步详细说明。\n[0047] 本发明提供一种细粒度免校验的室内被动入侵检测的方法，利用从普通商业网卡上获取的信道状态信息来评估室内多径信号变化表现的敏感度，进而预测新的入侵检测系统的检测阈值大小，以降低部署开销，增强无线室内入侵检测的应用性。\n[0048] 为实现上诉发明目的，本发明提供一种细粒度免校验的室内被动入侵检测的方法，其包括以下四个步骤：\n[0049] 步骤A：提取信号特征信息；\n[0050] 步骤B：构建标准信号特征文件，计算接收机敏感度因子大小；\n[0051] 步骤C：预测入侵检测阈值大小；\n[0052] 步骤D：计算检测信号特征与标准信号特征相关性，判断是否有入侵人体出现。\n[0053] 步骤A中的信号特征信息是指子载波的振幅信息，其所述提取信号特征信息的步骤包括：\n[0054] 步骤A1：采集信道状态信息；\n[0055] 步骤A2：提取各个子载波的振幅值作为信号特征；\n[0056] 步骤B中的构建标准信号特征文件，计算接收机敏感度因子大小的步骤为：\n[0057] 步骤B1：收集标准信号特征信息，将其存储在本地文件中；\n[0058] 步骤B2：利用标准信号特征文件计算接收机敏感度因子。\n[0059] 步骤C中的预测入侵检测阈值大小是建立在预测入侵检测阈值模型上。预测入侵检测阈值模型是指系统入侵检测阈值与接收机敏感度因子成线性关系。\n[0060] 步骤D中的计算检测信号特征与标准信号特征相关性的步骤为：\n[0061] 步骤D1：采集检测信号特征信息；\n[0062] 步骤D2：计算检测信号特征与标准信号特征相关性；\n[0063] 步骤D3：判断其相关性是否小于阈值，如果小于阈值则判断有入侵人体出现，否则无人出现。\n[0064] 进一步地，接收机和发射机的位置部署可以是任意的，也可以是现场调试的。对于需要较敏感的监测区域，则需要根据敏感度因子来调整接收机的位置。\n[0065] 图1是本发明第一实施例提供的细粒度免校验的室内被动入侵检测的实现流程图，该方法详述如下：\n[0066] 在步骤101中，提取信号特征，从信道状态信息中提取频域上所有子载波的信道状态信息的振幅向量。\n[0067] 在无线传输过程中，无线通信可以被简单建模为：\n[0068]\n[0069] 这里t是时间，y是接收信号，x是发送信号，h是信道相应或者信道状态信息，通常情况下其为一个复数，代表信号振幅与相位信息，z为高斯白噪音。在传统的室内环境中，一个传输信号可以通过多径进行传播，并导致不同的传播长度、路径损失、不同的时延、振幅衰减和相位偏移。而多径环境可以通过时间线性过滤器h(τ)特征化,即信道脉冲响应(Channel Impulse Response，CIR)：\n[0070]\n[0071] 其中，ai、θi和τi分别代表第i个多径的振幅、相位和时延，i＝1,2,…,N，N表示多径数。\n[0072] 在频率域中，OFDM系统在OFDM子载波粒度上提供信道频率响应(Channel Frequency Response，CFR)：\n[0073] H＝{H(1),H(2),...,H(N)}         (3)\n[0074] 其中N为子载波个数。每个子载波的CFR信息是一个复数值，每个子载波又被定义为：\n[0075] H(f)＝|H(f)|exp(jsin(∠H(f)))         (4)\n[0076] 其中H(f)表示子载波的振幅响应，而∠H(f)表示子载波的相位响应。\n[0077] 对于给定的一个带宽，CIR可以通过快速傅里叶变换转换成CFR：\n[0078] H＝FFT(h(τ))        (5)\n[0079] 尽管CIR与CFR在信道相应模型上是等价的。但是在人体检测与定位领域更加偏向使用CFR作为信号特征。目前已可以从商业Intel 5300无线网卡上利用更新后的固件提取具有30个子载波的CFR，并能以信道状态信息的形式提交到用户态进行程序处理。\n[0080] 每一个子载波的信道状态信息是以一个复数的形式来表示的，其包含了每个子载波的振幅与相位信息。通过计算每个子载波信道状态信息的振幅值，程序可以获取一个包含所有子载波振幅的向量，将其作为信号特征。\n[0081] 在步骤102中，当监测区域无人的情况下，系统采集n个数据包并从中随机选取m个数据包，提取每一个数据包中的信号特征 然后构建标准信号特征文件Hstd：\n[0082]\n[0083] 在步骤103中，系统计算接收机敏感度因子。为了刻画接收机对于人体移动的敏感度，我们引入了接收机敏感度因子变量。接收机敏感度因子可以在一定程度上反映出多径信号占接收机总信号比重大小。接收机敏感度因子的计算是利用标准信号特征文件。首先需要将每一个信号特征进行振幅归一化，将每个子载波的振幅归一化到中心频率上：\n[0084]\n[0085] 其中H(fk)和Hnorm(fk)分别是第k个子载波原始振幅值和归一化后的振幅值。f0是中心频率值，fk是第k个子载波的频率值。\n[0086] 为了消除测量尺度影响，获得一个统一的数据量纲，本发明需计算归一化后的振幅向量的变异系数：\n[0087]\n[0088] 其中std(Hnorm)和mean(Hnorm)分别是归一化CSI振幅Hnorm的标准差与均值。因此利用静态标准信号特征文件可以计算得到变异系数集 当敏感度较大时，CV值具有较大值和较宽的分布。当敏感度较小时，CV值则具有较小值和较窄的分布。为了刻画CV的分布与敏感度大小，本发明引入了接收机敏感度因子变量Ks：\n[0089]\n[0090] 其中，x，μ和σ分别代表ν的测量值、均值和标准差值。\n[0091] 在步骤104中，系统根据已计算的接收机敏感度因子对入侵检测阈值进行估计。检测阈值C与敏感度因子Ks有如下关系：\n[0092] C＝a·Ks+b               (10)\n[0093] 利用两者之间成线性关系。预测模型中的参数a和b可通过传统入侵检测阈值测量方法测得。\n[0094] 在步骤105中，接收机收集检测信号的信道状态信息，并从中提取检测信号特征Htest，信号特征提取过程同步骤101。然后系统计算每个检测信号特征Hdetect与标准信号特征文件中信号特征集合Hstd的相关性，其计算方法公式如下：\n[0095]\n[0096] 在步骤106中，系统对比步骤105中获得的Cdetect与检测阈值的大小。如果Cdetect小于系统预测的检测阈值，则说明监测区域内有入侵人体出现，发出入侵警报；否则继续收集检测信号的信道状态信息。\n[0097] 图2是本发明第二实施例提供的细粒度免校验的室内被动入侵检测的方法实现的示意图。信号发射机(TX)采用目前应用最为普及的商用无线路由器(AP)，而信号接收机(RX)为商用电脑，该电脑上配置有Intel NIC 5300网卡，支持802.11n协议，安装Linux操作系统和Linux CSI Tool工具。Linux CSI Tool包含针对Intel NIC 5300网卡的一款驱动程序iwlwifi，它可以获取OFDM系统中30个子载波的信道响应信息，并能够以信道状态信息(CSI)形式提交给相应程序进行处理。信号发射机和信号接收机一般被固定在离地面高度\n1.2m的位置，这样能更准确的对人体的出现进行检测和定位。信号接收机以一定速率(例如设定为每秒钟20个数据包)向信号发射机发送ICMP请求报文，然后信号接收机利用LinuxCSI Tool驱动从信号发射机的应答包中获取CSI信息，并实时将收集的CSI信息通过UDP协议发送给人体检测和定位的装置。系统部署期间，测试人员应首先对系统参数进行设置，包括ICMP包发送速率、应用服务器IP地址参数等相应信息。\n[0098] 当设备部署完成后，系统采集信道状态信息。安装人员开启标准信号特征采集功能，并确保监测区域内无人出现。接收机会以20Hz的频率向发射机发射ICMP报文，接收机通过网卡固件获取信道状态信息，该信道状态信息被上传到用户态数据发送模块。数据发送模块利用UDP协议，将信道状态信息发送给检测系统。检测系统数据接收模块接收机发送来的数据后，利用CSI Tool工具提取信号特征。系统会自动计算天线上信号特征中每个子载波的振幅，并将其存储在标准信号特征文件中。采集过程大概需要5分钟左右。系统将会提取标准信号特征文件中每个数据包的信号特征，计算每个数据包的信号特征变异系数，将其存入内存数组中。当提取所有的数据包的变异系数后，系统按照图1中步骤102的方法计算出检测系统在此场景下的敏感度因子。利用系统预测模型(步骤103)来预测出此场景下的检测阈值。\n[0099] 当完成了环境适应过程后，系统开启在线检测功能。接收机以每秒20Hz的频率向接收机发送ICMP包，并从网卡中获取信道信息，并通过数据发送模块利用UDP协议将其实时地发送给执行检测系统的设备。当检测系统接收到信道信息后，利用CSI Tool工具提取数据包的信号特征，并将其存入到时间窗口缓冲区中。当缓冲区被填满后，系统读取缓冲区内所有的信号特征，求其接收机天线上每个子载波的振幅均值。然后系统读取标准文件中的信号特征，计算接收机天线上的新的信号特征与标准静态信号特征的相关性。如果相关性大于检测阈值则将缓冲区清空，继续下一轮检测；否则意味有入侵人体出现在监测区域。