CN109884585B - 一种基于波形边缘检测的室内信号到达时间测量方法 - Google Patents
一种基于波形边缘检测的室内信号到达时间测量方法 Download PDFInfo
- Publication number
- CN109884585B CN109884585B CN201910122291.4A CN201910122291A CN109884585B CN 109884585 B CN109884585 B CN 109884585B CN 201910122291 A CN201910122291 A CN 201910122291A CN 109884585 B CN109884585 B CN 109884585B
- Authority
- CN
- China
- Prior art keywords
- signal
- state information
- channel state
- time
- max
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
- Radio Transmission System (AREA)
Abstract
一种基于波形边缘检测的室内信号达到时间测量方法,该方法基于对波形的采样,采用空‑频‑时三维搜索算法,计算某一子载波间隔一个采样周期的相邻区间分别对应的信道状态信息的相位差,通过相位差判断估计信道状态信息值的信号处理区间是否进入波形边缘,得到信号到达时间的估计值。该方法克服了传统子空间方法的缺陷,为信号到达时间估计提供了一条新的思路。
Description
技术领域
本发明涉及一种室内定位技术,具体涉及室内定位技术的室内信号到达时间测量方法。
背景技术
位置信息,是移动互联网的基本要素,是物联网中重要的感知信息,和人们的生产生活密切相关。室外定位和基于位置的服务已经发展成熟,基于卫星导航和地图的位置服务被广泛应用,并成为各种移动设备使用最多的应用之一。而人类80%以上的时间是在室内环境中活动,如果能够解决室内定位问题,就能实现室内和室外的无缝衔接,是导航领域的重大突破。
近十五年来,各种各样的室内定位技术层出不穷,包括谷歌、微软、苹果、博通等在内的一些科技巨头,以及世界著名的大学都在研究室内定位技术,其未来主要的应用领域包括:室内精准导航、大数据分析、个性化营销、社交网络等。随着智能手机的普及,以及移动互联网的发展,室内定位市场必将蓬勃发展、需求不断攀升,其商业化也将带来一波创新高潮。
室内定位技术主要有:红外线、超声波、射频识别、蓝牙、超宽带、以及无线局域网(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的室内定位技术具有广阔的前景和重要的意义。
Wi-Fi室内定位技术主要分为两类,一类是根据信号的指纹参数特征进行指纹匹配;另一类是根据信号的到达时间(TOA)进行几何解算。采用几何解算的现有技术方案中,最常用的是子空间算法。
假设室内多径信道的等效低通信道冲击响应可以表示为
其中,Lp是多径的数量。αk和τk分别是复衰落系数和传播时延。τk(0≤k≤Lp-1)按照时延长短从小到大排列,τ0代表第一径(如果直达径存在,则代表直达径)的传播时延,即TOA。对上式进行傅立叶变换,信道的频域响应可以表示为
由于Wi-Fi信号采用正交频分复用(OFDM)体制进行传输,假设OFDM信号一共有L个子载波,测量得到第l个子载波的信道状态信息(信道状态信息)可以表示为
x=H+w=Va+w
其中
x=[x(0) x(1) … x(L-1)]T
H=[H(f0) H(f1) … H(fL-1)]T
w=[ω(0) ω(1) … ω(L-1)]T
v(τk)=[1 exp(-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、时间搜索:确定子载波和接收天线后,判断判决准则的相位从某一固定值到发生改变的采样时刻,所述判断准则为:其相位定义为:其中是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式,定义为:TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);发射端信号在t时刻发射,该信号在t+t0时刻到达接收端,Δt是符号同步定时误差,Ts是采样间隔,T是符号时长,di是第i个子载波的已知符号;
S4、计算并输出:遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值,输出信号的到达时间。
进一步地,
所述的步骤S1包括:
步骤S1空域搜索范围从1到N,N为天线数量;
Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通过正交频分复用(OFDM)波形传输,假设OFDM波形的子载波数量为L,发射端信号在t时刻发射,在t+t0时刻到达接收端;传输信号的等效复信号为:
其中,αk和τk分别是第k条路径的衰落系数和时延,δ(t)是冲击函数,接收端等效复信号可以表示为:
其中,LPF表示低通滤波;发射端信号在t=0时刻发射,则信号在t=t0时刻到达接收端;信道状态信息通过非高吞吐量模式下的长训练符号L-LTF(Non-HT Long Trainingfield,)的两个TFFT部分进行估计,TFFT是完成快速傅里叶变换(FFT)运算的时域区间,推导出
其中,di是第i个子载波的已知符号,Δt是符号同步定时误差,Ts是采样间隔,TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);
所述的步骤S2包括:
SCImax为子载波的编号的最大值,频域搜索范围为从-SCImax到SCImax,0除外;
所述的步骤S3包括:
相位差:Δ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值;
所述的步骤S4包括:
本发明的优点在于新方法克服了子空间方法的缺陷,为TOA估计提供了一条新的思路。
附图说明
图1物理层汇聚协议数据单元的信号格式
图2室内多径信道传输时间示意图
具体实施方式
本专利基于信道状态信息,因此需要描述信道状态信息的获取。Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过OFDM波形传输(OFDM:Orthogonal FrequencyDivision Multiplexing即正交频分复用技术),汇聚协议数据单元的信号格式如图1所示,假设OFDM波形的子载波数量为L,载波频率为fc,传输信号可以表示为Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通常是通过正交频分复用(OFDM)波形传输,假设OFDM波形的子载波数量为L,发射端信号在t时刻发射,在t+t0时刻到达接收端;传输信号的等效复信号为:
其中,αk和τk分别是第k条路径的衰落系数和时延,δ(t)是冲击函数,接收端等效复信号可以表示为:
其中,表示低通滤波;发射端信号在t=0时刻发射,则信号在t=t0时刻到达接收端;信道状态信息通过L-LTF(Non-HT Long Training field,非高吞吐量模式下的长训练符号)的两个TFFT(完成快速傅里叶变换(FFT)运算的时域区间)部分进行估计,推导出
其中,是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式,di是第i个子载波的已知符号,Δt是符号同步定时误差,Ts是采样间隔,TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…)。
用于信道状态信息估计的L-LTF部分采用窗函数进行成形,其数学表达式为
其中,TFFT=3.2μs,T=2TFFT=6.4μs,G=1.6μs,TTR=100ns。
信道状态信息信号处理区间的终止位置为t0-Δt+TFFT+mTs(m=0,1,2,…),如图2所示,假设室内多径信道有M条径,若终止位置为t0-Δt+TFFT+δTs,信号处理区间恰好没有跨入窗函数右边缘,即第一径没有跨入窗函数右边缘;若终止位置为t0-Δt+TFFT+(δ+1)Ts,信号处理区间恰好跨入窗函数右边缘,即第一径刚好跨入窗函数右边缘。此时,第一径跨入窗函数右边缘的时长为ξ1;第二径跨入窗函数右边缘的时长为ξ2;第U径跨入窗函数右边缘的时长为ξU。假设判决准则为我们采用以下准则
从(6)式可以得出判断第一径是否刚好跨入窗函数右边缘方法:取某一子载波相邻信号处理区间(即利用[t0-Δt+(m-1)Ts,t0-Δt+TFFT+(m-1)Ts]和[t0-Δt+mTs,t0-Δt+TFFT+mTs])计算出信道状态信息值后取其共轭的相位,判断相位从某一固定值到发生变化的m=δ,输出t0-Δt+δTs作为第一径达到接收机的参考值,即TOA的估计值。
具体的算法可以描述为
以上仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内所做的任何修改、等同替换和改进等,均应包含在本发明创造的保护范围之内。
Claims (2)
1.一种基于波形边缘检测的室内信号到达时间测量方法,其特征在于:所述方法包括以下步骤:
S1、空域搜索:对接收天线从1到N进行编号,从1到N进行空域搜索,确定信道状态信息幅度值;
S2、频域搜索:子载波数量为L,对子载波从-SCImax到+SCImax进行编号,并进行频域搜索,对于某一个子载波,确定信道状态信息幅度值最大的接收天线的编号,其中L=2SCImax;
S3、时间搜索:确定子载波和接收天线后,判断判决准则的相位从某一固定值到发生改变的采样时刻,所述判断准则为:其相位定义为:其中是完成符号同步后从第m个采样点开始计算第i个子载波信道状态信息的表达式,定义为:TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);发射端信号在t时刻发射,该信号在t+t0时刻到达接收端,Δt是符号同步定时误差,Ts是采样间隔,T是符号时长,di是第i个子载波的已知符号;
S4、计算并输出:遍历不同子载波并将出现频率最高的那个采样时刻作为采样时刻的估计值,输出信号的到达时间。
2.根据权利要求1所述的一种基于波形边缘检测的室内信号到达时间测量方法,其特征在于,
所述的步骤S1包括:
步骤S1空域搜索范围从1到N,N为天线数量;
Wi-Fi系统的物理层汇聚协议数据单元(PPDU)通过正交频分复用(OFDM)波形传输,假设OFDM波形的子载波数量为L,发射端信号在t时刻发射,在t+t0时刻到达接收端;传输信号的等效复信号为:
其中,αk和τk分别是第k条路径的衰落系数和时延,δ(t)是冲击函数,接收端等效复信号可以表示为:
其中,LPF表示低通滤波;发射端信号在t=0时刻发射,则信号在t=t0时刻到达接收端;信道状态信息通过非高吞吐量模式下的长训练符号L-LTF(Non-HT Long Training field)的两个TFFT部分进行估计,TFFT是完成快速傅里叶变换(FFT)运算的时域区间,推导出
其中,di是第i个子载波的已知符号,Δt是符号同步定时误差,Ts是采样间隔,TFFT是估计信道状态信息需要覆盖的信号处理区间,区间的起始位置为t=t0-Δt+mTs(m=0,1,2,…);
所述的步骤S2包括:
SCImax为子载波的编号的最大值,频域搜索范围为从-SCImax到SCImax,0除外;
所述的步骤S3包括:
arg [ · ] 表示取相位值;其中δ是arg[Ψ(i,m)]的计算结果由某一恒定值发生改变时对应的区间起始位置t=t0-Δt+mTs中的m值;
相位差:Δ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值;
所述的步骤S4包括:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910122291.4A CN109884585B (zh) | 2019-02-19 | 2019-02-19 | 一种基于波形边缘检测的室内信号到达时间测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910122291.4A CN109884585B (zh) | 2019-02-19 | 2019-02-19 | 一种基于波形边缘检测的室内信号到达时间测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109884585A CN109884585A (zh) | 2019-06-14 |
CN109884585B true CN109884585B (zh) | 2020-03-17 |
Family
ID=66928438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910122291.4A Active CN109884585B (zh) | 2019-02-19 | 2019-02-19 | 一种基于波形边缘检测的室内信号到达时间测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109884585B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104812061B (zh) * | 2015-03-24 | 2019-03-22 | 成都希盟泰克科技发展有限公司 | 一种基于mimo-ofdm信道状态信息的室内测距及定位方法 |
CN106851554A (zh) * | 2016-10-09 | 2017-06-13 | 滁州学院 | 基于OFDM信道状态信息相位变化特征的高精度室内WiFi定位系统 |
CN107347210B (zh) * | 2017-07-04 | 2019-10-29 | 江苏先安科技有限公司 | 一种基于信道状态信息的精确目标定位方法 |
-
2019
- 2019-02-19 CN CN201910122291.4A patent/CN109884585B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109884585A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100820618B1 (ko) | 스펙트럼 추정을 이용한 통신 시스템에 대한 채널 추정 | |
US7336647B2 (en) | System and method for ranging | |
CN110351655B (zh) | 一种基于信号多径传播测量的室内定位方法及系统 | |
US10349213B2 (en) | Method and apparatus for determining an angle of arrival of a wireless signal | |
CN106879068B (zh) | 一种强多径环境下信号的到达时间估计方法 | |
CN103095638B (zh) | 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 | |
EP3125483A1 (en) | Communication system determining time of arrival using matching pursuit | |
WO2010139234A1 (zh) | 最大多普勒频偏的估计方法和装置 | |
EP3554138A1 (en) | Method for determining candidate line-of-sight path, and wireless local area network device | |
JP6626246B2 (ja) | 多帯域ワイヤレスローカルエリアネットワーク内の早期フレーム帯域評価 | |
US9374245B2 (en) | Observation of the true channel from band-limited frequency domain observations | |
TW200847705A (en) | Method of synchronizing multi-carrier systems and multi-carrier system | |
CN112954791B (zh) | 一种基于子载波筛选的信道状态信息csi定位方法 | |
CN111539323B (zh) | 一种循环前缀线性调频信号的频率估计方法与装置 | |
Reddy et al. | An improved time-of-arrival estimation for WLAN-based local positioning | |
CN109884585B (zh) | 一种基于波形边缘检测的室内信号到达时间测量方法 | |
Han et al. | A new high precise indoor localization approach using single access point | |
KR20050020622A (ko) | 수신 장치 및 수신 방법 | |
CN110290084A (zh) | 一种基于数据频率能量峰值的短波信道盲符号同步方法 | |
CN113422748B (zh) | 窄带物联网时频偏的估计方法、装置及存储介质 | |
CN110611629B (zh) | 一种频率偏差的估计方法、装置及通信设备 | |
Zhang et al. | Pilot-based delay and Doppler estimation in 6G integrated communication and sensing networks | |
WO2013040772A1 (zh) | 一种基于信道估计的定位方法及装置 | |
CN109782223B (zh) | 一种基于信号指纹匹配的室内定位方法 | |
Haghparast et al. | OFDM ranging in multipath channels using time reversal method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |