本发明专利技术公开了一种利用信道状态信息对民用小型无人机进行实时入侵检测的方法，包括信道状态信息获取、信道状态信息特征提取、基于这些特征的状态分类等步骤；通过IEEE 802.11n标准的设备获得信道状态信息，对信道状态信息进行特征提取，基于提取得到的特征通过分类算法进行状态分类；当环境周围存在无人机时进行报警。本发明专利技术采用符合IEEE 802.11n标准的设备，通过检测信道状态信息的方法，避免传统的基于视觉和音频识别的无人机检测方法中成本高与干扰因素大的问题。在检测过程中，通过算法计算与设备功耗无关的相关特征，通过支持向量机(Support Vector Machine)与长短期记忆网络(Long Short‑Term Memory)进行分类，采用阈值算法，进行无人机的入侵检测。
【技术实现步骤摘要】
一种利用信道状态信息对民用小型无人机进行实时入侵检测的方法
本专利技术属于信号处理与特征提取领域，涉及民用小型无人机入侵时对信道状态信息的影响问题，尤其是民用小型无人机入侵一般住宅时对信道的影响。
技术介绍
当今，民用小型无人直升机越来越普及，被广泛地用于航拍等。然而，由于民用无人机监管力度不足以及监管困难的问题，民用无人机被用于“黑飞”、“滥飞”的事件时常发生。考虑到民用无人机不仅对机场这类重要场有设备安全、人生安全的威胁，一些重要机构甚至是一般住宅都需要实时检测无人机的入侵，因为民用小型无人机会可能携带爆炸物或放射性物质等有害物品，并且可能对隐私甚至是重要机密进行窃听或者窥视。现有检测无人机主要通过雷达检测、机器视觉和音频识别的方法。基于雷达的方法设备昂贵，并且属于主动发射电磁波的检测方式，不仅辐射大而且功耗大，不适合用于一般楼宇或者住宅的无人机检测。另外，如今一些无人机对雷达是隐形的。基于机器视觉的方法，受天气和遮挡物的影响巨大，并且需要多个摄像头的部署，同样不适合用于一般楼宇或者住宅的无人机检测。基于音频识别的方法，其全向麦克风或者麦克风阵列设备同样昂贵。针对现有检测无人机的方法中存在的问题，设计一种低成本、不受天气影响、低功耗的适合一般楼宇或者住宅部署的无人机检测方法显得尤为重要。
技术实现思路
本专利技术的目的在于针对现有检测无人机的方法中存在的问题，设计一种低成本、不受天气影响、低功耗的适合一般楼宇或者住宅部署的无人机检测方法，该方法通过检测无人机对信道状态信息的影响，通过机器学习算法，实现无人机的检测。本专利技术解决其技术问题是通过以下技术方案来实现的：一种利用信道状态信息(ChannelStateInformation，简称CSI)对民用小型无人机进行实时入侵检测的方法，包括以下步骤：(1)信道状态信息获取，该步骤包括以下子步骤：(1.1)部署无线发射机与接收机；(1.2)提取接收机中的信道状态信息，读出一组包含N个子载波的信道频率响应信息，即CSI矩阵，该矩阵的规模为Package×Ntx×Nrx×N，其中Package为CSI数据包的数量，Ntx为发射机使用天线数，Nrx为接收机使用天线数，N取30或56；(2)信道状态信息特征提取，该步骤包括以下子步骤：(2.1)在步骤(1.2)获得的CSI矩阵中，将发射天线tx与接收天线rx之间的信道频率响应记为H(tx,rx)，则H(tx,rx)＝[H(f1),H(f2),…,H(fN)]，其中0≤tx＜Ntx,0≤rx＜Nrx；H(fk)为Package行的复数向量，可以表示为1≤k≤N；(2.2)设时间窗口为K，每次分析K个CSI数据包的信息；在一次分析中，每对天线之间的归一化CSI矩阵可以表示为其中Hi(1≤i≤K)为在一个时间窗口内，序号为i的数据包中包含的H(tx,rx)；(2.3)使用线性变换方法降低相位偏移；若N＝30，具体包括以下子步骤：(a)将测量得到的相位表示为：其中φi为真实相位，δ为接收机相对于发射机的时钟偏移，其产生的对应相位偏移为β是未知的相位偏置，Z是测量噪声，ki表示第i个子载波的子载波索引，Nc表示FFT的大小；(b)定义算子其中i∈{-28,-24,-20,-16,-12,-8,-4,-1,1,5,9,13,17,21,25,28},j∈{-26,-22,-18,-14,-10,-6,-2,3,7,11,15,19,23,27}忽略测量噪声Z，并计算：根据算子的定义，b可以化简为(c)对测量所得的相位进行以下线性运算：其中为能够反映真实相位的相位指标；可以化简为：将中测量得到的相位替换为若N＝56，具体包括以下子步骤：(A)将测量得到的相位表示为：其中φi为真实相位，δ为接收机相对于发射机的时钟偏移，其产生的对应相位偏移为β是未知的相位偏置，Z是测量噪声，ki表示第i个子载波的子载波索引，Nc表示FFT的大小；(B)忽略测量噪声Z，计算：因为N＝56时，ki是对称的，所以有b可以化简为(C)对测量所得的相位进行以下线性运算：其中为能够反映真实相位的相位指标；可以化简为：将中测量得到的相位替换为(2.4)去除中的异常值；(2.5)记和分别为归一化后的CSI幅度序列向量和相位序列向量，计算的协方差矩阵和的协方差矩阵(2.6)计算的最大特征值m1与第二大特征值m2：计算的最大特征值n1与第二大特征值n2；(2.7)在Ntx×Nrx对天线的数据得到的m1,m2,n1,n2中，分别去掉最大与最小值，取剩下数据的平均值，得到：(3)将步骤(2.7)中得到的作为支持向量机分类算法或者长短期记忆网络算法的输入进行学习，将得到的模型用于判断环境周围有无无人机入侵。进一步地，所述步骤(1.2)中，若N＝30，代表CSI包含正交频分复用技术中的30个子载波，其子载波索引为{-28,-26,-24,-22,-20,-18,-16,-14,-12,-10,-8,-6,-4,-2,-1,1,3,5,7,9,11,13,15,17,19,21,23,25,27,28}；若N＝56，代表CSI包含正交频分复用技术中的56个子载波，其子载波索引为{-28,-27,…,-2,-1,1,2,…,27,28}。进一步地，所述步骤(2.4)具体为：利用Hampel辨识器，将中落在区间[μ-γσ,μ+γσ]外的测量值判定为异常值，并滤去，其中μ为测量值序列的中位数，σ为绝对中位差，γ为滤波器参数。进一步地，所述步骤(2.5)中，协方差矩阵计算公式如下：其中，表示向量与的协方差，表示归一化的N向量。本专利技术与现有技术相比具有以下显著的优点：1)发射机可以采用无线路由器，接收机可以采用无线网卡，都是方便获得的低成本设备。2)设备部署方便，不需要部署设备阵列。3)设备功耗低。4)设备受天气影响小。5)改进了相位噪声去除算法，使得在子载波索引不对称的情况下仍能完全消去未知相位偏置与由时钟偏移带来的相位偏移。附图说明图1是本专利技术的无人机检测方法流程图。图2是本专利技术中发射机与接收机的部署示例；具体实施方式下面结合附图和具体实施例对本专利技术作进一步详细说明。如图1所示，本专利技术提供的一种利用信道状态信息对民用小型无人机进行实时入侵检测的方法，包括如下步骤：(1)信道状态信息获得，该步骤包括以下子步骤：(1.1)按图2所示部署无线发射机与接收机。尽量增大无线信道在室外的路径，可以增大无人机对信道状态的影响，从而使无人机入侵时的特征更明显。其中无线发射机可以使用无线路由器，接收机可以采用无线网卡。只需要这样较低成本的设备，甚至是家家户户本身就配备有的设备即可作为信号采集装置；(1.2)提取无线网卡中的信道状态信息，读出一组包含N个子载波的信道频率响应信息，即CSI矩阵，该矩阵的规模为Package×Ntx×Nr××N，其中Package为CSI数据包的数量，Ntx为发射机使用天线数，Nrx为接收机使用天线数，N取30或56；若N＝30，代表CSI包含正交频分复用技术中的30个子载波，其子载波索引为{-28,-26,-24,-22,-20,-18,-16,-14,-12,-10,-8,-6,-4,-2,-1,1,3,5,7,9,11,13,15,17,19,21,23,25,本文档来自技高网...
【技术保护点】
1.一种利用信道状态信息对民用小型无人机进行实时入侵检测的方法，其特征在于，包括以下步骤：(1)信道状态信息获取，该步骤包括以下子步骤：(1.1)部署无线发射机与接收机；(1.2)提取接收机中的信道状态信息，读出一组包含N个子载波的信道频率响应信息，即CSI矩阵，该矩阵的规模为Package×Ntx×Nrx×N，其中Package为CSI数据包的数量，Ntx为发射机使用天线数，Nrx为接收机使用天线数，N取30或56；(2)信道状态信息特征提取，该步骤包括以下子步骤：(2.1)在步骤(1.2)获得的CSI矩阵中，将发射天线tx与接收天线rx之间的信道频率响
【技术特征摘要】
1.一种利用信道状态信息对民用小型无人机进行实时入侵检测的方法，其特征在于，包括以下步骤：(1)信道状态信息获取，该步骤包括以下子步骤：(1.1)部署无线发射机与接收机；(1.2)提取接收机中的信道状态信息，读出一组包含N个子载波的信道频率响应信息，即CSI矩阵，该矩阵的规模为Package×Ntx×Nrx×N，其中Package为CSI数据包的数量，Ntx为发射机使用天线数，Nrx为接收机使用天线数，N取30或56；(2)信道状态信息特征提取，该步骤包括以下子步骤：(2.1)在步骤(1.2)获得的CSI矩阵中，将发射天线tx与接收天线rx之间的信道频率响应记为H(tx,rx)，则H(tx,rx)＝[H(f1),H(f2),…,H(fN)]，其中0≤tx＜Ntx,0≤rx＜Nrx；H(fk)为Package行的复数向量，可以表示为(2.2)设时间窗口为K，每次分析K个CSI数据包的信息；在一次分析中，每对天线之间的归一化CSI矩阵可以表示为其中Hi(1≤i≤K)为在一个时间窗口内，序号为i的数据包中包含的H(tx,rx)；(2.3)使用线性变换方法降低相位偏移；若N＝30，具体包括以下子步骤：(a)将测量得到的相位表示为：其中φi为真实相位，δ为接收机相对于发射机的时钟偏移，其产生的对应相位偏移为β是未知的相位偏置，Z是测量噪声，ki表示第i个子载波的子载波索引，Nc表示FFT的大小；(b)定义算子其中i∈{-28,-24,-20,-16,-12,-8,-4,-1,1,5,9,13,17,21,25,28},j∈{-26,-22,-18,-14,-10,-6,-2,3,7,11,15,19,23,27}忽略测量噪声Z，并计算：根据算子的定义，b可以化简为(c)对测量所得的相位进行以下线性运算：其中为能够反映真实相位的相位指标；可以化简为：将中测量得到的相位替换为若N＝56，具体包括以下子步骤：(A)将测量得到的相位表示为：其中φi为真实相位，δ为接收机相对于发射机的时钟偏移，其产生的对应相位偏移为β是未知的相位偏置，Z是测量噪声，ki表示第i个子载波的子载波索引，N...
【专利技术属性】
技术研发人员：邓江毅，郑润禾，杨晗，徐文渊，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。