CN104243064A - 一种非视距密集多径场景下的扩频toa测距方法 - Google Patents
一种非视距密集多径场景下的扩频toa测距方法 Download PDFInfo
- Publication number
- CN104243064A CN104243064A CN201410437925.2A CN201410437925A CN104243064A CN 104243064 A CN104243064 A CN 104243064A CN 201410437925 A CN201410437925 A CN 201410437925A CN 104243064 A CN104243064 A CN 104243064A
- Authority
- CN
- China
- Prior art keywords
- correlation function
- peak
- slip
- slip correlation
- toa
- 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.)
- Granted
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明提出了一种室内非视距密集多径场景下基于到达时间(Time of Arrival,TOA)的测距方法。将接收端扩频信号与原发射端扩频信号进行滑动相关,得到滑动相关函数,在滑动相关函数中进行首径搜索确定首径时延,进而由TOA测距求出待测距离。为了提高测距精度,对首径搜索中的判决门限进行了优化,找到了最优的判决门限,使平均测距误差最小。测量结果表明,在典型室内非视距场景下,本方法的首径搜索平均误差小于24.67ns,平均测距误差小于7.4m。
Description
技术领域
本发明属于信号检测与估计技术领域,特别涉及一种到达时间估计的测距方法。
背景技术
室内定位技术能够在住宅楼、办公楼、图书馆、地下车库、货品仓库等复杂环境下实现对人员以及物品的迅速定位。室内定位在人们的生活工作及科学研究中起到了非常重要的作用。实现非视距(Non Line of Sight,NLOS)室内场景下有效对抗多径效应的测距方法是实现室内定位技术的基本条件。对于视距(Line of Sight,LOS)路径,接收端可以通过寻找最大径功率信号来求取首径信号,进而求取信号到达时间(Time of Arrival,TOA),得到目标距离。然而室内NLOS环境存在密集多径,大量多径信号叠加得到的最大功率相关峰不一定是首径信号所在的相关峰。这种现象导致首径信号难以提取,TOA结果误差较大。基于此背景,本发明提出了一种非视距密集多径场景扩频TOA测距方法。首先,应用滑动相关技术将接收端扩频信号与发射端扩频信号进行滑动相关,得到滑动相关函数。其次,设计首径搜索算法,对滑动相关函数进行首径搜索得到首径所在相关峰的TOA。最终,采用TOA估计待测距离。在此基础上,为了提高测距精度,通过实际测量对首径搜索中使用的判决门限进行优化。
基于TOA的扩频测距是利用扩频序列的相关特性解算出扩频码在空间的传输时延,进而求出目标体的距离。扩频测距具有测距精度高、抗干扰性能好、适应性强等优点。基于TOA的扩频测距方法目前主要应用于导航定位、卫星测控***中。在测距***中,对距离的测定就是对到达扩频序列和本地扩频序列时间差的估计和跟踪。扩频测距***的具体工作原理:在发送端,信源经过扩频和调制,由发射机发射出去。接收端,信号经过相关等处理,由TOA估计算法得到发射信号和接收信号的时延差,进而求得到达时间TOA,从而计算出待测目标和发送端之间的距离,完成测距任务。
TOA技术由于其算法相对简单且精度较高得到广范应用。GPS(GlobalPosition System)等TOA测距***主要应用于LOS场景,这类场景下存在大功率视距路径,基于滑动相关方法获取的相关函数中最大相关峰的时延与首径时延基本相等,测距***具有较高精度。而在室内NLOS场景下,由于存在密集多径,在滑动相关的过程中,很多密集多径在同一个相关峰内大量叠加,导致通过滑动相关测得的最大相关峰时延未必是首径时延,甚至最大相关峰未必是首径所在的相关峰。现有文献都采用波峰能量检测将能量最大的波峰作为首径信号所在的相关峰得到TOA。其缺陷是,在室内NLOS下由于存在密集的多径,在作相关时很多多径在同一个波峰内大量叠加,使得最大峰值不一定是首径信号,测距存在较大误差。在室内NLOS测距领域,迄今为止,尚没有能够有效抵抗密集多径影响的TOA测距方法。
发明内容
为了克服现有技术中的缺陷,本发明提出了一种非视距NLOS密集多径场景下的扩频TOA测距方法,该方法能够有效抵抗密集多径对TOA估计的影响。
本发明通过以下技术方案实现:
一种非视距密集多径场景下的扩频到达时间测距方法,包括以下步骤:
a.扩频测距***发射扩频码信号,经过无线信道传播后接收机获得接收信号,将接收端信号与发射信号进行滑动相关,获取n个滑动相关函数;这里的n个滑动相关函数是由n个收发天线间的接收信号与其发射信号进行滑动相关得到的;
b.在信道测量时,同时测量收发天线间的直线距离di,1≤i≤n,进而得到精确时延其中c为光速;
c.选取判决门限p;
d.对于第i个滑动相关函数,将所述滑动相关函数波形的所有峰值找到,记最大峰值为am;
e.首径搜索:由判决门限p得到搜索区域D=(pam,am),即所述滑动相关函数中由上下两个边界am、pam确定的一块条形区域;找到D内的所有波峰,时延最小的峰值记为ap,ap即在判决门限p下首径搜索得到的首径信号所在的相关峰的峰值,记ap的时延为τi;。
f.首径搜索误差ei=τi-ξi,对应于不同的判决门限p,平均误差为
g.判断是否为最小值,若不是最小值,则回到步骤c,重新选取判决门限p;若是最小值,则此时对应的p为最优判决门限,τi为首径时延;
h.根据τi获得测距结果。
附图说明
图1是本发明的扩频TOA测距方法的滑动相关过程示意图;
图2是典型滑动相关函数图,其中首径是最大峰;
图3是典型滑动相关函数图,其中首径不是最大峰;
图4是本发明的扩频TOA测距方法的首径搜索流程图;
图5是本发明的首径搜索示意图,其中首径在最大峰;
图6是本发明的首径搜索示意图,其中首径不在最大峰;
图7是采用本发明的首径搜索方法得到的不同判决门限下首径搜索的平均误差图。
具体实施方式
下面结合附图说明及具体实施方式对本发明进一步说明。
由于存在密集多径,为了获得TOA,必须求取首径信号的时延。扩频测距***中,发射端的扩频信号一般是由扩频码调制的载波,扩频码(如m序列)具有良好的自相关性,其经过无线信道传播后,可以根据其相关波形确定传输时延。基于此原理,本发明将扩频测距***中接收端信号与原发射信号进行滑动相关,通过滑动相关函数来提取首径。滑动相关的原理是,将接收到的信号延时不同的时间,分别与原发射信号进行匹配滤波,得到不同延时的匹配结果,即要得到的滑动相关函数。具体的滑动相关过程如附图1所示。
通过以上过程得到的典型滑动相关函数波形如附图2、3所示。由附图2、3可知,典型的滑动相关函数是由若干个相关峰构成。由于存在多径效应,滑动相关函数的波峰实际上是由多个到达径信号的叠加而成。由于室内NLOS场景下,信道存在密集多径,多径的功率和时延都是随机的。因而,最大功率的相关峰未必是首径所在的相关峰。必须采用合理的方法从滑动相关函数中准确提取首径信号所在的相关峰。
本发明中采用左右差分法寻找滑动相关峰,具体为:
对于滑动相关函数曲线的某点,记其坐标为(k,f(k)),左移一个单位后该点对应的滑动相关函数值为f(k-1),右移一个单位后该时刻对应的滑动相关函数值为f(k+1),若有
成立,则该点为滑动相关波峰;若该式不成立,则该点不是滑动相关波峰。
传统方法认为最大相关峰就是首径所在的相关峰(如附图2),然而由于信道存在密集多径分量,多径分量相互叠加导致首径信号所在相关峰的功率有可能小于其他多径所在相关峰的功率,此时首径不在滑动相关函数的最大相关峰内,利用最大相关峰的时延进行TOA估计将会引入很大误差(如附图3)。综合考虑以上情况,本发明提出了一种NLOS密集多径场景扩频TOA测距方法的首径搜索方法,其流程图如附图4所示,包括如下步骤:
1)将滑动相关函数波形的所有峰值找到,记最大峰值为am,如附图2中的最大波峰为第一个波峰,且此时am=10-5.075。
2)确定判决门限p(dB),由p得到搜索区域D=(pam,am),即滑动相关函数中由上下两个边界确定的一块条形区域。
3)找到D内的所有波峰,时延最小的峰值记为ap,ap即在判决门限p下首径搜索得到的首径信号所在的相关峰的峰值,记ap的时延为τi。
4)τi即为首径时延。
该方法实际上是在滑动相关函数波形中确定一个搜索区域。即附图5、6中两条虚线之间区域。对于不在搜索区域的相关峰,噪声干扰严重,信号功率足够小,判决时舍弃这些相关峰。对于在搜索区域内的相关峰,寻找其中时延最小的相关峰,即为首径所在的相关峰。
附图5中,当判决门限p=-36.65dB时,搜索区域范围为(10-6.909,10-5.075)。图中有A1,B1,C1,D1,E1,F1,G1共7个相关峰,而在搜索区域内只有A1,B1,G1三个相关峰,时延最小的峰是A1,时延为100ns,进而首径搜索获得的首径信号时延为100ns。附图5的滑动相关函数首径信号时延就是峰值最大的相关峰的时延。
附图6中,当判决门限p=-25.3dB时,搜索区域范围(10-6.843,10-5.579)。图中有A2,B2,C2,D2,E2,F2共6个相关峰,位于搜索区域内有A2,B2,C2,D2,F2共5个相关峰,时延最小的峰是A2峰,时延为40ns,则首径搜索获得的首径信号时延为40ns。附图6的滑动相关函数首径信号时延并不是峰值最大的相关峰的时延,与附图5有所区别。
使用首径搜索方法获得TOA后,由信号传播速度,进而得到待测距离。TOA的精度直接决定测距精度。
上述方法寻找的首径具有较高的精度。为了提高首径搜索精度,有必要对判决门限p进行优化,寻找最优的判决门限,使平均搜索误差降到最小,大幅提高测距精度。
本发明采用的首径搜索方法,运用到了判决门限p,对于不同的p值,首径搜索的效果不同。为了提高首径搜索的精度,必须对所有判决门限值,分别使用上述首径搜索方法,从而确定一个最优的判决门限。
为了获得精确的首径时延,在信道测量时,同时测量收发天线间的直线距离di(1≤i≤n),进而得到精确时延其中c为光速。在本次信道测量时,n取1440。
首径搜索误差ei即为图2中滑动相关函数得到的首径τi与精确时延ξi之间的差值ei=τi-ξi,对应于不同的判决门限p,求取对应的平均误差
取使公式(2)达到最小的判决门限p,此时的门限p即为本发明得到的最优判决门限。
为了确定最优判决门限,需要获取足够数量的滑动相关函数。为此,进行了典型室内NLOS信道测量,得到了1440个信道测量数据,对应不同的信道数据,分别将发射端扩频信号通过信道,仿真出接收信号,将接收信号与发射信号作滑动相关,得到1440个滑动相关函数后,应用本发明的首径搜索方法,得到不同判决门限下首径搜索的平均误差,如附图7所示。
由图7可知,存在一个判决门限,使得首径搜索平均误差达到最小,该门限值为-18.7dB,此时的首径搜索平均误差为测距平均误差为7.4m,实现了判决门限的优化。
本发明提出了一种非视距密集多径场景下的扩频TOA测距方法。将接收端扩频信号与原发射端扩频信号进行滑动相关,得到滑动相关函数,在滑动相关函数中进行首径搜索确定首径时延,进而由TOA测距求出待测距离。为了提高测距精度,对首径搜索中的判决门限进行了优化,找到了最优的判决门限,使平均测距误差最小。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (3)
1.一种非视距密集多径场景下的扩频到达时间测距方法,其特征在于:所述方法包括以下步骤:
a.将接收端信号与发射信号进行滑动相关,获取n个滑动相关函数;
b.在信道测量时,同时测量收发天线间的直线距离di,1≤i≤n,进而得到精确时延其中c为光速;
c.选取判决门限p;
d.对于第i个滑动相关函数,将所述滑动相关函数波形的所有峰值找到,记最大峰值为am;
e.首径搜索:由判决门限p得到搜索区域D=(pam,am),即所述滑动相关函数中由上下两个边界am、pam确定的一块条形区域;找到D内的所有波峰,时延最小的峰值记为ap,ap即在判决门限p下首径搜索得到的首径信号所在的相关峰的峰值,记ap的时延为τi;。
f.首径搜索误差ei=τi-ξi,对应于不同的判决门限p,平均误差为
g.判断是否为最小值,若不是最小值,则回到步骤c,重新选取判决门限p;若是最小值,则此时对应的p为最优判决门限,τi为首径时延;
h.根据τi获得测距结果。
2.根据权利要求1所述的方法,其特征在于:所述步骤a具体为:扩频测距***发射扩频码信号,经过无线信道传播后接收机获得接收信号,将接收到的信号延时不同的时间,分别与原发射信号进行匹配滤波,得到不同延时的匹配结果,即为所获取滑动相关函数。
3.根据权利要求1所述的方法,其特征在于:所述步骤b中采用左右差分法寻找滑动相关峰,具体为:
对于所述滑动相关函数曲线的某点,记其坐标为(k,f(k)),左移一个单位后该点对应的滑动相关函数值为f(k-1),右移一个单位后该时刻对应的滑动相关函数值为f(k+1),若有
成立,则该点为滑动相关波峰;若该式不成立,则该点不是滑动相关波峰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410437925.2A CN104243064B (zh) | 2014-08-29 | 2014-08-29 | 一种非视距密集多径场景下的扩频toa测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410437925.2A CN104243064B (zh) | 2014-08-29 | 2014-08-29 | 一种非视距密集多径场景下的扩频toa测距方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104243064A true CN104243064A (zh) | 2014-12-24 |
CN104243064B CN104243064B (zh) | 2016-07-06 |
Family
ID=52230499
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410437925.2A Active CN104243064B (zh) | 2014-08-29 | 2014-08-29 | 一种非视距密集多径场景下的扩频toa测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104243064B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107483395A (zh) * | 2016-06-08 | 2017-12-15 | 恩智浦有限公司 | 用于通信装置的处理模块和其方法 |
CN107959513A (zh) * | 2016-10-13 | 2018-04-24 | 大唐移动通信设备有限公司 | 一种采用时延参数进行测距的方法和装置 |
CN111381225A (zh) * | 2020-04-09 | 2020-07-07 | 中国电子科技集团公司第二十八研究所 | 一种基于vdb***的民航进场着陆高精度测距***及测距方法 |
CN114826331A (zh) * | 2022-04-22 | 2022-07-29 | 青岛鼎信通讯股份有限公司 | 一种低压线路载波通信信号多径时延特性测试方法 |
CN116299166A (zh) * | 2023-05-24 | 2023-06-23 | 四川思凌科微电子有限公司 | 一种用于chirp信号的低复杂度融合测距方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001281330A (ja) * | 2000-03-30 | 2001-10-10 | Clarion Co Ltd | スペクトル拡散通信測距装置 |
WO2007056662A2 (en) * | 2005-11-02 | 2007-05-18 | Sirf Technology, Inc. | System and method for detecting multipath effects in a gps receiver |
CN101394233A (zh) * | 2007-09-21 | 2009-03-25 | 哈尔滨工业大学深圳研究生院 | 室内视距环境下脉冲超宽带多径信号建模的方法及*** |
US20120063491A1 (en) * | 2010-09-13 | 2012-03-15 | Korea Advanced Institute Of Science And Technology | Spread spectrum signal receiver, method for multipath super-resolution thereof, and recording medium thereof |
CN103118429A (zh) * | 2013-03-12 | 2013-05-22 | 南京烽火星空通信发展有限公司 | Wifi定位方法 |
CN103220781A (zh) * | 2013-04-24 | 2013-07-24 | 哈尔滨工业大学深圳研究生院 | 无线传感器网络静止目标定位方法及*** |
-
2014
- 2014-08-29 CN CN201410437925.2A patent/CN104243064B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001281330A (ja) * | 2000-03-30 | 2001-10-10 | Clarion Co Ltd | スペクトル拡散通信測距装置 |
WO2007056662A2 (en) * | 2005-11-02 | 2007-05-18 | Sirf Technology, Inc. | System and method for detecting multipath effects in a gps receiver |
CN101394233A (zh) * | 2007-09-21 | 2009-03-25 | 哈尔滨工业大学深圳研究生院 | 室内视距环境下脉冲超宽带多径信号建模的方法及*** |
US20120063491A1 (en) * | 2010-09-13 | 2012-03-15 | Korea Advanced Institute Of Science And Technology | Spread spectrum signal receiver, method for multipath super-resolution thereof, and recording medium thereof |
CN103118429A (zh) * | 2013-03-12 | 2013-05-22 | 南京烽火星空通信发展有限公司 | Wifi定位方法 |
CN103220781A (zh) * | 2013-04-24 | 2013-07-24 | 哈尔滨工业大学深圳研究生院 | 无线传感器网络静止目标定位方法及*** |
Non-Patent Citations (4)
Title |
---|
史庭祥等: "时空调制扩频通信定位***抗多径测角测距方法", 《应用科学学报》 * |
汪洋等: "UWB脉冲信号室内空旷环境路径损耗特性研究", 《哈尔滨工业大学学报》 * |
潘仲明等: "超声波扩频测距及其信道自适应均衡技术", 《国防科技大学学报》 * |
王忠等: "扩频测距***中多径跟踪误差研究", 《***工程与电子技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107483395A (zh) * | 2016-06-08 | 2017-12-15 | 恩智浦有限公司 | 用于通信装置的处理模块和其方法 |
US10715355B2 (en) | 2016-06-08 | 2020-07-14 | Nxp B.V. | Processing module for a communication device and method therefor |
US10771288B2 (en) | 2016-06-08 | 2020-09-08 | Nxp B.V. | Processing module for a communication device and method therefor |
CN107483395B (zh) * | 2016-06-08 | 2020-11-06 | 恩智浦有限公司 | 用于通信装置的处理模块和其方法 |
CN107959513A (zh) * | 2016-10-13 | 2018-04-24 | 大唐移动通信设备有限公司 | 一种采用时延参数进行测距的方法和装置 |
CN107959513B (zh) * | 2016-10-13 | 2019-11-22 | 大唐移动通信设备有限公司 | 一种采用时延参数进行测距的方法和装置 |
CN111381225A (zh) * | 2020-04-09 | 2020-07-07 | 中国电子科技集团公司第二十八研究所 | 一种基于vdb***的民航进场着陆高精度测距***及测距方法 |
CN111381225B (zh) * | 2020-04-09 | 2022-03-22 | 中国电子科技集团公司第二十八研究所 | 一种基于vdb***的民航进场着陆高精度测距***及测距方法 |
CN114826331A (zh) * | 2022-04-22 | 2022-07-29 | 青岛鼎信通讯股份有限公司 | 一种低压线路载波通信信号多径时延特性测试方法 |
CN116299166A (zh) * | 2023-05-24 | 2023-06-23 | 四川思凌科微电子有限公司 | 一种用于chirp信号的低复杂度融合测距方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104243064B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104243064B (zh) | 一种非视距密集多径场景下的扩频toa测距方法 | |
CN101606076B (zh) | 对距离估计或涉及距离估计的改进 | |
US8587478B1 (en) | Localization method of multiple jammers based on TDOA method | |
US20050124368A1 (en) | Method distinguishing line of sight (los) from non-line of sight (nlos) in cdma mobile communication system | |
CN104181503A (zh) | 一种基于rss的led定位***及方法 | |
CN101883426A (zh) | 高精度的超宽带无线定位方法 | |
CN103308883B (zh) | 一种基于单天线的到达角估计方法 | |
Zhang et al. | A low complexity NLOS error mitigation method in UWB localization | |
CN106879068A (zh) | 一种强多径环境下信号的到达时间估计方法 | |
CN104459675A (zh) | 一种基于测距的物体定位和跟踪方法及使用该方法的定位设备 | |
CN104219761A (zh) | 一种基于最大斜率的超宽带无线定位方法 | |
CN103592624B (zh) | 一种基于接收信号强度的测距方法 | |
Hao et al. | A method for improving UWB indoor positioning | |
Zhu et al. | EKF localization based on TDOA/RSS in underground mines using UWB ranging | |
Svecova et al. | Object localization using round trip propagation time measurements | |
Kabakchiev et al. | Moving Target FSR Shadow Detection using GPS signals | |
Zhang et al. | An improved Taylor series based location algorithm for IEEE 802.15. 4a channels | |
Li et al. | Performance of TOA estimation techniques in indoor multipath channels | |
CN1474196A (zh) | 无线定位中高精度、抗多径的窄带相关测距*** | |
KR100740702B1 (ko) | 시변조 uwb 신호의 직선 경로 신호 검출방법 및 그 장치 | |
Mary et al. | Test measurements of improved UWB localization technique for precision automobile parking | |
Awarkeh et al. | Improvement of the angle measurement accuracy for indoor UWB localization | |
CN102662168A (zh) | 基于uwb技术的楼宇定位方法 | |
He et al. | Modelling ranging error of indoor TOA–based position system based on receive signal strength | |
Elbahhar et al. | Indoor positioning system based on the ultra wide band for transport applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |