一种基于波形边缘检测的室内信号到达时间测量方法技术

一种基于波形边缘检测的室内信号达到时间测量方法，该方法基于对波形的采样，采用空‑频‑时三维搜索算法，计算某一子载波间隔一个采样周期的相邻区间分别对应的信道状态信息的相位差，通过相位差判断估计信道状态信息值的信号处理区间是否进入波形边缘，得到信号到达时间的估计值。该方法克服了传统子空间方法的缺陷，为信号到达时间估计提供了一条新的思路。

【技术实现步骤摘要】
一种基于波形边缘检测的室内信号达到时间测量方法
本专利技术涉及一种室内定位技术，具体涉及室内定位技术的室内信号达到时间测量方法。
技术介绍
位置信息，是移动互联网的基本要素，是物联网中重要的感知信息，和人们的生产生活密切相关。室外定位和基于位置的服务已经发展成熟，基于卫星导航和地图的位置服务被广泛应用，并成为各种移动设备使用最多的应用之一。而人类80％以上的时间是在室内环境中活动，如果能够解决室内定位问题，就能实现室内和室外的无缝衔接，是导航领域的重大突破。近十五年来，各种各样的室内定位技术层出不穷，包括谷歌、微软、苹果、博通等在内的一些科技巨头，以及世界著名的大学都在研究室内定位技术，其未来主要的应用领域包括：室内精准导航、大数据分析、个性化营销、社交网络等。随着智能手机的普及，以及移动互联网的发展，室内定位市场必将蓬勃发展、需求不断攀升，其商业化也将带来一波创新高潮。室内定位技术主要有：红外线、超声波、射频识别、蓝牙、超宽带、以及无线局域网(WLAN，即Wi-Fi，基于IEEE802.11协议)等。与Wi-Fi相比，其它室内定位系统均需要部署专用的硬件设备，导致成本高昂；Wi-Fi技术自诞生以来，拥有庞大的客户群体且已经广泛部署。截止到2017年末，Wi-Fi设备累计出货量接近200亿台。Wi-Fi不仅可以连接笔记本和智能手机等上网设备，还可以连接电视、音响、灯具等家用电器，甚至将渗透进入智能联网汽车等领域。未来，Wi-Fi将是物联网和5G网络的重要支撑技术。人们除了通过Wi-Fi设备进行信息传输外，还可以利用Wi-Fi设备进行精确定位。尤其是室内，卫星导航等技术由于信号遮挡，导航接收机接收到的信号受到衰落和散射，无法正确解析当前位置，导致位置信息无法获得。因此，利用广泛部署的Wi-Fi设备接续定位功能具有重要的现实意义。目前基于Wi-Fi的室内定位精度为3米左右，还没有实现高精度的定位。因此，研究基于Wi-Fi的室内定位技术具有广阔的前景和重要的意义。Wi-Fi室内定位技术主要分为两类，一类是根据信号的指纹参数特征进行指纹匹配；另一类是根据信号的到达时间(TOA)进行几何解算。采用几何解算的现有技术方案中，最常用的是子空间算法。假设室内多径信道的等效低通信道冲击响应可以表示为其中，Lp是多径的数量。αk和τk分别是复衰落系数和传播时延。τk(0≤k≤Lp-1)按照时延长短从小到大排列，τ0代表第一径(如果直达径存在，则代表直达径)的传播时延，即TOA。对(1)式进行傅立叶变换，信道的频域响应可以表示为由于Wi-Fi信号采用正交频分复用(OFDM)体制进行传输，假设OFDM信号一共有L个子载波，测量得到第l个子载波的信道状态信息(信道状态信息)可以表示为其中，l＝0,1,…,L-1，f0是信号下边缘载波频率，ω(l)是噪声，均值为零，方差为(3)式可以用矩阵和矢量表示为x＝H+w＝Va+w其中x＝[x(0)x(1)...x(L-1)]TH＝[H(f0)H(f1)...H(fL-1)]Tw＝[ω(0)ω(1)...ω(L-1)]Tv(τk)＝[1exp(-j2π△fτk)…exp[-j2π(L-1)△fτk]]α'k＝αkexp(-j2πf0τk)若采用MUSIC超分辨率算法，首先求解x的互相关矩阵其中，A＝E{aaH}。将Rxx进行特征值分解后按照特征值从大到小排列，第i个特征值对应的特征矢量为qi，假设L>Lp，噪声空间可以表示为MUSIC谱可以表示为：式的第一个峰值对应的τ值就是TOA的估计值τ0。采用子空间的方法估计TOA存在三个方面的问题：一、需要满足子载波数量大于多径数量的条件，在实际中这个条件很难满足。室内多径数量通常趋向于无穷多，即使只考虑能量较大的多径，此条件也很难满足。二、由于信号带宽有限，导致时间分辨率不足。例如通常的Wi-Fi信号带宽为20MHz，时间分辨率为1/20MHz＝5ns，空间分辨率为1.5米，使得实现亚米级定位精度存在一定困难。三、估计x的互相关矩阵Rxx需要一定数量的采样快拍，而实际信号体制能提供的快拍数量可能无法满足要求。
技术实现思路
本专利技术设计了一种新的TOA估计方法，新方法基于对波形的采样，计算某一子载波间隔一个采样周期的相邻区间分别对应的信道状态信息的相位差，通过相位差判断估计信道状态信息值的信号处理区间是否进入波形边缘，得到TOA的估计值。一种基于波形边缘检测的室内信号达到时间测量方法，其特征在于：所述方法包括以下步骤:S1、空域搜索：对接收天线从1到N进行编号，从1到N进行空域搜索，确定信道状态信息幅度值；S2、频域搜索：子载波数量为L，对子载波从-SCImax到+SCImax进行编号，并进行频域搜索，对于某一个子载波，确定信道状态信息幅度值最大的接收天线的编号，其中L＝2SCImax；S3、时间搜索：确定子载波和接收天线后，判断准则相位从某一固定值到发生改变的采样时刻；S4、计算并输出：遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值，输出信号的到达时间。所述的步骤S1、S2、S3、S4包括以下步骤：其中，S1进行空域搜索包括以下步骤：步骤S1空域搜索范围从1到N，N为天线数量；是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式，定义为公式(5)，推导过程如下：Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过正交频分复用(OFDM)波形传输，假设OFDM波形的子载波数量为L，发射端信号在t时刻发射，在t+t0时刻到达接收端；传输信号的等效复信号为：其中，是传输的调制信号，Ps是发射的信号功率，载波频率为fc，是相偏，T是符号时长，q指第q个子载波，j为虚数单位，△fc是频偏；多径信道由M条路径组成，是t时刻多径信道的信道冲击响应函数，表示为：其中，αk和τk分别是第k条路径的衰落系数和时延，δ(t)是冲击函数，接收端等效复信号可以表示为：其中，是接收端噪声，*是卷积操作；经过下变频和低通滤波，接收端信号可以表示为：其中，LPF[·]表示低通滤波；发射端信号在t＝0时刻发射，则信号在t＝t0时刻到达接收端；信道状态信息通过L-LTF(Non-HTLongTrainingfield，非高吞吐量模式下的长训练符号)的两个TFFT(完成快速傅里叶变换(FFT)运算的时域区间)部分进行估计，推导出其中，di是第i个子载波的已知符号，△t是符号同步定时误差，Ts是采样间隔，TFFT是估计信道状态信息需要覆盖的信号处理区间，区间的起始位置为t＝t0-△t+mTs(m＝0,1,2,…)；空域搜索遍历从1到N所有天线，得到对应信道状态信息幅度S2进行频域搜索包括以下步骤：SCImax为子载波的编号的最大值频域搜索范围为从-SCImax到SCImax，0除外；是从t＝ts+mTs时刻开始第k个天线接收到的第i个子载波计算出的信道状态信息，号天线对应的幅度值最大，则输出其编号作为信道状态信息幅度值最大的接收天线；S3中时间搜索包括以下步骤：时间搜索范围从2到是上取整符；计算信道状态信息共轭值的相位其中，为判决准则，定义为公式(6)arg[·]表示取相位值相位差：△k(i,m)＝|arg[Ψk(i,m)]-ar本文档来自技高网...

【技术保护点】
1.一种基于波形边缘检测的室内信号达到时间测量方法，其特征在于：所述方法包括以下步骤：S1、空域搜索：对接收天线从1到N进行编号，从1到N进行空域搜索，确定信道状态信息幅度值；S2、频域搜索：子载波数量为L，对子载波从‑SCImax到+SCImax进行编号，并进行频域搜索，对于某一个子载波，确定信道状态信息幅度值最大的接收天线的编号，其中L＝2SCImax；S3、时间搜索：确定子载波和接收天线后，判断准则相位从某一固定值到发生改变的采样时刻；S4、计算并输出：遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值，输出信号的到达时间。

【技术特征摘要】
1.一种基于波形边缘检测的室内信号达到时间测量方法，其特征在于：所述方法包括以下步骤：S1、空域搜索：对接收天线从1到N进行编号，从1到N进行空域搜索，确定信道状态信息幅度值；S2、频域搜索：子载波数量为L，对子载波从-SCImax到+SCImax进行编号，并进行频域搜索，对于某一个子载波，确定信道状态信息幅度值最大的接收天线的编号，其中L＝2SCImax；S3、时间搜索：确定子载波和接收天线后，判断准则相位从某一固定值到发生改变的采样时刻；S4、计算并输出：遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值，输出信号的到达时间。2.根据权利要求1所述的一种基于波形边缘检测的室内信号达到时间测量方法，其特征在于，所述的步骤S1、S2、S3、S4包括以下步骤：其中，权利要求1中所述的S1进行空域搜索包括以下步骤：步骤S1空域搜索范围从1到N，N为天线数量；是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式，定义为公式(5)，推导过程如下：Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过正交频分复用(OFDM)波形传输，假设OFDM波形的子载波数量为L，发射端信号在t时刻发射，在t+t0时刻到达接收端；传输信号的等效复信号为：其中，是传输的调制信号，Ps是发射的信号功率，载波频率为fc，是相偏，T是符号时长，q指第q个子载波，j为虚数单位，△fc是频偏；多径信道由M条路径组成，是t时刻多径信道的信道冲击响应函数，表示为：其中，αk和τk分别是第k条路径的衰落系数和时延，δ(t)是冲击函数，接收端等效复信号可以表示为：其中，是接收端噪声，*是卷积操作；经过下变频和低通滤波，接收端信号可以表示为：其中，LPF[·]表示低通滤波；发射端信号在t＝0时刻发射，则信号在t＝t0时刻到达接收端；信道状态信息通过L-LTF(Non-HTLongTrainingfield，非高吞吐...

【专利技术属性】
技术研发人员：李敏，王琰，于全，柯明星，张驭龙，
申请(专利权)人：军事科学院系统工程研究院网络信息研究所，
类型：发明
国别省市：北京,11