技术领域
[0001] 本
发明涉及一种室内
定位技术,具体涉及室内定位技术的室内信号到达时间测量方法。
背景技术
[0002]
位置信息,是移动互联网的基本要素,是
物联网中重要的
感知信息,和人们的生产生活密切相关。室外定位和基于位置的服务已经发展成熟,基于卫星导航和地图的位置服务被广泛应用,并成为各种移动设备使用最多的应用之一。而人类80%以上的时间是在室内环境中活动,如果能够解决室内定位问题,就能实现室内和室外的无缝衔接,是导航领域的重大突破。
[0003] 近十五年来,各种各样的室内定位技术层出不穷,包括谷歌、微软、苹果、博通等在内的一些科技巨头,以及世界著名的大学都在研究室内定位技术,其未来主要的应用领域包括:室内精准导航、
大数据分析、个性化营销、社交网络等。随着智能手机的普及,以及移动互联网的发展,室内定位市场必将蓬勃发展、需求不断攀升,其商业化也将带来一波创新高潮。
[0004] 室内定位技术主要有:红外线、
超声波、
射频识别、蓝牙、超宽带、以及无线局域网(WLAN,即Wi-Fi,基于IEEE 802.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的室内定位技术具有广阔的前景和重要的意义。
[0005] Wi-Fi室内定位技术主要分为两类,一类是根据信号的指纹参数特征进行指纹匹配;另一类是根据信号的到达时间(TOA)进行几何解算。采用几何解算的
现有技术方案中,最常用的是子空间
算法。
[0006] 假设室内多径信道的等效低通信道冲击响应可以表示为
[0007]
[0008] 其中,Lp是多径的数量。αk和τk分别是复衰落系数和传播时延。τk(0≤k≤Lp-1)按照时延长短从小到大排列,τ0代表第一径(如果直达径存在,则代表直达径)的传播时延,即TOA。对上式进行
傅立叶变换,信道的频域响应可以表示为
[0009]
[0010] 由于Wi-Fi信号采用
正交频分复用(OFDM)体制进行传输,假设OFDM信号一共有L个
子载波,测量得到第l个子载波的信道状态信息(信道状态信息)可以表示为[0011]
[0012] 其中,l=0,1,…,L-1,f0是信号下边缘载波
频率,ω(l)是噪声,均值为零,方差为用矩阵和矢量表示为:
[0013] x=H+w=Va+w
[0014] 其中
[0015] x=[x(0) x(1) … x(L-1)]T
[0016] H=[H(f0) H(f1) … H(fL-1)]T
[0017] w=[ω(0) ω(1) … ω(L-1)]T
[0018]
[0019] v(τk)=[1 exp(-j2πΔfτk) … exp[-j2π(L-1)Δfτk]]
[0020]
[0021] α'k=αkexp(-j2πf0τk)
[0022] 若采用MUSIC超
分辨率算法,首先求解x的互相关矩阵
[0023]
[0024] 其中,A=E{aaH}。将Rxx进行特征值分解后按照特征值从大到小排列,第i个特征值对应的特征矢量为qi,假设L>Lp,噪声空间可以表示为
[0025]
[0026] MUSIC谱可以表示为:
[0027]
[0028] 式的第一个峰值对应的τ值就是TOA的估计值τ0。
[0029] 采用子空间的方法估计TOA存在三个方面的问题:一、需要满足子载
波数量大于多径数量的条件,在实际中这个条件很难满足。室内多径数量通常趋向于无穷多,即使只考虑
能量较大的多径,此条件也很难满足。二、由于信号带宽有限,导致
时间分辨率不足。例如通常的Wi-Fi信号带宽为20MHz,时间分辨率为1/20MHz=5ns,空间分辨率为1.5米,使得实现亚米级定位精度存在一定困难。三、估计x的互相关矩阵Rxx需要一定数量的
采样快拍,而实际信号体制能提供的快拍数量可能无法满足要求。
发明内容
[0030] 本发明设计了一种新的TOA估计方法,新方法基于对波形的采样,计算某一子载波间隔一个采样周期的相邻区间分别对应的信道状态信息的
相位差,通过
相位差判断估计信道状态信息值的
信号处理区间是否进入波形边缘,得到TOA的估计值。
[0031] 一种基于波形边缘检测的室内信号到达时间测量方法,包括以下步骤:
[0032] S1、
空域搜索:对接收天线从1到N进行编号,从1到N进行空域搜索,确定信道状态信息幅度值;
[0033] S2、频域搜索:子载波数量为L,对子载波从-SCImax到+SCImax进行编号,并进行频域搜索,对于某一个子载波,确定信道状态信息幅度值最大的接收天线的编号,其中L=2SCImax;
[0034] S3、时间搜索:确
定子载波和接收天线后,判断判决准则的相位从某一固定值到发生改变的采样时刻 ,所述判断准则为 : 其 相位 定义为:其中 是完成符号同步后从第m个采样点开始计算第i
个子载波信道状态信息的表达式,定义为:
TFFT是估计信道状态信息需
要
覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);发射端信号在t时刻发射,该信号在t+t0时刻到达接收端,Δt是符号同步定时误差,Ts是采样间隔,T是符号时长,di是第i个子载波的已知符号;
[0035] S4、计算并输出:遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值,
输出信号的到达时间。
[0036] 进一步地,
[0037] 所述的步骤S1包括:
[0038] 步骤S1空域搜索范围从1到N,N为天线数量;
[0039] 是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式, 定义为公式(5),推导过程如下:
[0040] Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通过正交频分复用(OFDM)波形传输,假设OFDM波形的子载波数量为L,发射端信号在t时刻发射,在t+t0时刻到达接收端;传输信号的等效复信号为:
[0041]
[0042] 其中, 是传输的调制信号,Ps是发射的信号功率,
载波频率为fc, 是相偏,T是符号时长,q指第q个子载波,j为虚数单位,Δfc是频偏;
[0043] 多径信道由M条路径组成, 是t时刻多径信道的信道冲击响应函数,表示为:
[0044]
[0045] 其中,αk和τk分别是第k条路径的衰落系数和时延,δ(t)是冲击函数,接收端等效复信号可以表示为:
[0046]
[0047] 其中, 是接收端噪声,*是卷积操作;经过下变频和低通滤波,接收端信号可以表示为:
[0048]
[0049] 其中,LPF表示低通滤波;发射端信号在t=0时刻发射,则信号在t=t0时刻到达接收端;信道状态信息通过非高吞吐量模式下的长训练符号L-LTF(Non-HT Long Training field,)的两个TFFT部分进行估计,TFFT是完成快速傅里叶变换(FFT)运算的时域区间,推导出
[0050]
[0051] 其中,di是第i个子载波的已知符号,Δt是符号同步定时误差,Ts是采样间隔,TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);
[0052] 空域搜索遍历从1到N所有天线,得到对应信道状态信息 的幅度其中, 是从t=ts+mTs时刻开始第k个天线接收到的第i个子载波计算
出的信道状态信息;
[0053] 所述的步骤S2包括:
[0054] SCImax为子载波的编号的最大值,频域搜索范围为从-SCImax到SCImax,0除外;
[0055] 对于m号子载波,-SCImax≤m≤SCImax,如果 号天线对应的幅度值最大,则输出其编号 作为信道状态信息幅度值最大的接收天线;
[0056]
[0057] 所述的步骤S3包括:
[0058] 时间搜索范围从2到 是上取整符;
[0059] 计算信道状态信息共轭值的相位
[0060] 其中, 为判决准则,定义为公式(6)
[0061]
[0062] 表示取相位值;其中δ是arg[Ψ(i,m)]的计算结果由某一恒定值发生改变时对应的区间起始位置t=t0-Δt+mTs中的m值;
[0063] 相位差:Δk(i,m)=arg[Ψk(i,m)]-arg[Ψk(i,m-1)]|,计算间隔一个采样周期相位值的差; 表示第i个子载波满足Δk(i,m)/Δk(i,m-1)>α或者Δk(i,m)/Δk(i,m-1)<β的区间起始位置t=t0-Δt+mTs中的m值;
[0064] 如果Δk(i,m)/Δk(i,m-1)>α或者Δk(i,m)/Δk(i,m-1)<β,则 其中α和β为
门限值,α>1,0<β<1;当
信噪比区间为[0dB,20dB]时,α=1.25,β=0.8;
[0065] 所述的步骤S4包括:
[0066] 从-SCImax到SCImax遍历子载波,统计 的次数,其中,n为任意正整数;如果出现的次数最多,则输出信号到达时间的估计值为: Ts是采样间隔,ts=t0-Δt是L-LTF部分符号同步的起始时间。
[0067] 本发明的优点在于新方法克服了子空间方法的
缺陷,为TOA估计提供了一条新的思路。
附图说明
[0068] 图1物理层汇聚协议数据单元的信号格式
[0069] 图2室内多径信道传输时间示意图
具体实施方式
[0070] 本
专利基于信道状态信息,因此需要描述信道状态信息的获取。Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过OFDM波形传输(OFDM:Orthogonal Frequency Division Multiplexing即正交频分复用技术),汇聚协议数据单元的信号格式如图1所示,假设OFDM波形的子载波数量为L,载波频率为fc,传输信号可以表示为Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过正交频分复用(OFDM)波形传输,假设OFDM波形的子载波数量为L,发射端信号在t时刻发射,在t+t0时刻到达接收端;传输信号的等效复信号为:
[0071]
[0072] 其中, 是传输的调制信号,Ps是发射的信号功率,载波频率为fc, 是相偏,T是符号时长,q指第q个子载波,j为虚数单位,Δfc是频偏;
[0073] 多径信道由M条路径组成, 是t时刻多径信道的信道冲击响应函数,表示为:
[0074]
[0075] 其中,αk和τk分别是第k条路径的衰落系数和时延,δ(t)是冲击函数,接收端等效复信号可以表示为:
[0076]
[0077] 其中, 是接收端噪声,*是卷积操作;经过下变频和低通滤波,接收端信号可以表示为:
[0078]
[0079] 其中, 表示低通滤波;发射端信号在t=0时刻发射,则信号在t=t0时刻到达接收端;信道状态信息通过L-LTF(Non-HT Long Training field,非高吞吐量模式下的长训练符号)的两个TFFT(完成快速傅里叶变换(FFT)运算的时域区间)部分进行估计,推导出
[0080]
[0081] 其中, 是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式,di是第i个子载波的已知符号,Δt是符号同步定时误差,Ts是采样间隔,TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…)。
[0082] 用于信道状态信息估计的L-LTF部分采用窗函数进行成形,其数学表达式为[0083]
[0084] 其中,TFFT=3.2μs,T=2TFFT=6.4μs,G=1.6μs,TTR=100ns。
[0085] 信道状态信息信号处理区间的终止位置为t0-Δt+TFFT+mTs(m=0,1,2,…),如图2所示,假设室内多径信道有M条径,若终止位置为t0-Δt+TFFT+δTs,信号处理区间恰好没有跨入窗函数右边缘,即第一径没有跨入窗函数右边缘;若终止位置为t0-Δt+TFFT+(δ+1)Ts,信号处理区间恰好跨入窗函数右边缘,即第一径刚好跨入窗函数右边缘。此时,第一径跨入窗函数右边缘的时长为ξ1;第二径跨入窗函数右边缘的时长为ξ2;第U径跨入窗函数右边缘的时长为ξU。假设判决准则为 我们采用以下准则
[0086]
[0087] 来判断第一径是否刚好跨入窗函数右边缘。其中, 表示取相位值。
[0088] 从(6)式可以得出判断第一径是否刚好跨入窗函数右边缘方法:取某一子载波相邻信号处理区间(即利用[t0-Δt+(m-1)Ts,t0-Δt+TFFT+(m-1)Ts]和[t0-Δt+mTs,t0-Δt+TFFT+mTs])计算出信道状态信息值后取其共轭的相位,判断相位从某一固定值到发生变化的m=δ,输出t0-Δt+δTs作为第一径达到接收机的参考值,即TOA的估计值。
[0089] 具体的算法可以描述为
[0090]
[0091]
[0092]
[0093]
[0094] 以上仅为本发明创造的较佳
实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内所做的任何
修改、等同替换和改进等,均应包含在本发明创造的保护范围之内。