CN106959432A - 一种基于小波分解低频系数的海上作业平台人员定位方法 - Google Patents
一种基于小波分解低频系数的海上作业平台人员定位方法 Download PDFInfo
- Publication number
- CN106959432A CN106959432A CN201710177034.1A CN201710177034A CN106959432A CN 106959432 A CN106959432 A CN 106959432A CN 201710177034 A CN201710177034 A CN 201710177034A CN 106959432 A CN106959432 A CN 106959432A
- Authority
- CN
- China
- Prior art keywords
- energy block
- uwb signal
- low frequency
- wavelet decomposition
- kurtosis
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于海上作业平台人员定位领域,具体是一种基于小波分解低频系数的海上作业平台人员定位方法。主要分为六个步骤:A.安装在油区工作人员身上的定位标签发送UWB信号,定位基站接收该信号;B.定位基站对接收到的UWB信号进行小波分解;C.取小波分解后的低频系数构造能量块;D.建立阈值与能量块峭度的映射关系;E.实际测量时,根据收到的UWB信号,进行实际测距;F.服务器再根据定位标签与不同定位基站之间的距离差,利用双曲线算法,实现人员定位。本发明首次将UWB信号、小波分解、小波分解低频系数构造能量块等应用于海上作业平台人员定位中,在多径、非直视、多金属等海上油区特殊作业条件下可以提高人员测距定位精度。
Description
技术领域
本发明属于无线定位、测距技术领域,特别涉及超宽带无线定位技术,具体是一种基于小波分解低频系数的海上作业平台人员定位方法。
背景技术
随着我国国民经济的快速发展,能源问题日益突出,同时由于我国陆上石油资源的日益枯竭,海上石油开发已成为能源发展的战略重点。然而海上生产环境受天气、海浪等外界因素影响,易遇到突发性、灾难性气候,如遇人员落水的情况,需投入大量人力、物力进行拉网式搜索,并且搜救过程及效率受天气、海浪等外界因素影响,导致搜救周期长、搜救效果不理想。随着海上作业日益频繁,施工人员在作业平台上以及落水后的精确定位成为海上作业平台安全生产需要重点解决的问题。
目前常用的无线定位技术包括:①GNSS定位(Global NavigationSatelliteSystem,全球卫星导航系统),包括GPS(Global Positioning System,全球定位系统)和BDS(BeiDou Navigation Satellite System,北斗卫星导航系统),这些定位技术在室外理想情况下虽然可以实现小于10米的定位,但其信号一旦受到遮挡或者受到多径衰减等影响,将无法实现定位,所以在海上作业平台上基本无法使用。②短距离定位技术:目前常用的短距离无线定位技术主要包括UWB(Ultra-WideBand,超宽带)、Wi-Fi、红外线、超声波等。
UWB是一种新兴的高速、低成本、低功耗、高带宽的无线通信技术,与目前传统的无线定位技术相比具有以下5个主要的特点,所以更适合在海上作业平台使用。
1)定位精度高:由于UWB采用纳秒级窄脉冲,其持续时间远小于多径传播时延,所以在接收端多径信号能有效分离,从而具有较高的多径分辨率,因此可以实现厘米级别的高精度的测距和定位。
2)穿透能力强:由于UWB信号具有丰富的低频分量,所以具有很强的穿透树叶和障碍物的能力,从而可以实现室内和地下的精确定位。
3)发射功率低:由于UWB信号功率谱密度很低,所以所需平均功率很小。
4)安全性高:由于UWB信号功率谱密度很低,所以很难被检测到,所以安全性好,保密性高。
5)对其他系统干扰少:相对于其他通信系统,UWB信号所产生的干扰相当于宽带白噪声,这样可以使UWB系统与其他现有窄带通信系统共存,提高无线频谱资源的利用率,所以可以缓解海上作业平台日益紧张的无线频谱资源。
基于UWB的IEEE802.15.4a标准是IEEE于2007年最终通过的第一个精确测距定位的无线物理层标准,正是由于UWB独特的通信机制使其具有高精度、低成本、抗多径干扰、穿透能力强等特点,成为海上作业平台短距离、高精度无线定位的首选。
本发明专利正是在IEEE802.15.4a标准基础上提出的基于超宽带信号小波分解低频系数的提高海上作业平台人员定位精度的方法。该方法可以实现多径、多金属环境、有遮挡的海上作业平台环境下的高精度测距定位,从而为海上作业平台人员安全提供保障。
发明内容
鉴于目前海上作业平台人员定位系统多采用GPS,而GPS又有其不可克服的局限与不足,本发明提出了一种基于超宽带信号小波分解低频系数的海上作业平台人员定位方法,该方法可以大大提高传统定位方法中的多径、多金属环境、有遮挡环境下的定位精度,为海上作业平台作业人员提供安全保障。
一种基于超宽带信号小波分解低频系数的海上作业平台人员定位方法,包括以下步骤:
A.UWB信号发送:安装在作业平台工作人员身上的定位标签发送UWB信号,定位基站接收该信号。如果要实现三维定位则基站数量不能少于4个,如果要实现二维定位则基站数量不能少于3个。
B.UWB信号小波分解:定位基站对接收到的UWB信号进行小波分解,计算方法为:
其中,c0表示接收到的UWB信号,j∈{0,1,2,…}表示小波分解的层数,cj表示UWB信号经小波分解到j层上得到的低频部分的系数,dj表示UWB信号经小波分解到j层上得到的高频部分的系数。h(n)表示低通滤波器,g(n)表示高通滤波器,其通过选取的小波类型决定。
C.取低频系数构造能量块:每隔一定的区间对低频系数进行积分计算,构造一系列能量块,方法为:
其中,Tf表示UWB信号一个帧的长度,Tb表示积分周期即每个能量块的长度,r(t)表示接收到的UWB信号,z[n]表示能量块。
D.建立阈值与能量块峭度的映射关系:由低频能量块的峭度来计算阈值,具体方法如下:
1)生成指纹数据:在不同的信噪比下进行多次仿真,分别计算每一次仿真能量块的峭度并进行取整,同时归一化阈值取{0.1,0.2,…,1.0}之间不同值时,计算测距误差。峭度的计算公式如下:
其中,Nb表示能量块的个数,表示能量块的均值,δ表示能量块的标准差。
2)选取最优阈值:当能量块峭度取不同值时,分别计算取不同归一化阈值时测距误差的平均值,取平均误差最小时的归一化阈值作为此峭度下选取的归一化阈值。
3)曲线拟合:根据能量块峭度和归一化阈值的对应关系,利用曲线拟合建立起峭度和归一化阈值的映射关系。
E.测距计算:收到UWB信号后,在进行实际测距时,执行如下步骤:
1)构造信号低频系数能量块:对接收到的UWB信号进行小波变换,取低频系数构造能量块;
2)计算阈值:计算能量块峭度,根据建立的映射关系计算出阈值;
3)计算UWB信号到达时间:将能量块与阈值进行比较,取第一个超过阈值的能量块作为信号到达的时间TOA;
4)计算距离,计算方法如下:
其中,c表示电磁波在空气中的传播速度,表示估计的信号到达时间TOA。
F.定位计算:多个定位基站将上述测距结果通过网络传输到服务器,服务器再根据定位标签与不同定位基站之间的距离差,利用传统的双曲线无线定位算法,实现海上平台工作人员的定位。
本发明对比已有技术具有以下显著优点:
1.通过UWB信号发送接收、UWB信号小波分解、构造基于低频系数的能量块、建立阈值与能量块峭度的映射关系、测距计算、定位计算六个步骤,可以实现海上作业平台人员的高精度定位。
2.在多径、多金属、有遮挡等环境下,该方法测距精度优于其他常用的方法。
附图说明
图1是本发明方法的主要流程图。
图2是CM1信道接收信号和小波低频系数平均绝对误差随信噪比的变化。
图3是CM2信道接收信号和小波低频系数平均绝对误差随信噪比的变化。
具体实施方式
下面结合附图,说明本发明的实施方式。实施例中选用的信道是IEEE802.15.4a标准的CM1信道和CM2信道,主要实施步骤如图1所示,具体包括:
A.UWB信号发送:安装在作业平台工作人员身上的定位标签发送UWB信号,定位基站接收该信号;UWB信号采用的波形是二阶高斯脉冲、调制方式是PPM-TH-UWB、采样频率是50GHz。
B.UWB信号小波分解:定位基站对接收到的UWB信号进行小波分解,计算方法为:
其中,c0表示接收到的UWB信号,j∈{0,1}表示小波分解的层数,在本例中对接收到的UWB信号进行两层小波分解,cj表示UWB信号经小波分解到j层上得到的低频部分的系数,dj表示UWB信号经小波分解到j层上得到的高频部分的系数。h和g分别是由“db6”小波决定的低通和高通滤波器,h和g分别为[-0.0011,0.0048,0.0006,-0.0316,0.0275,0.0975,-0.1298,-0.2263,0.3153,0.7511,0.4946,0.1115],和[-0.1115,0.4946,-0.7511,0.3153,0.2263,-0.1298,-0.0975,0.0275,0.0316,0.0006,-0.0048,-0.0011]。
C.取低频系数构造能量块:每隔一定的区间对低频系数进行积分计算,构造一系列能量块,方法为:
其中,Tf表示UWB信号一个帧的长度,Tb表示积分周期即每个能量块的长度,r(t)表示接收到的UWB信号,z[n]表示能量块。
D.建立阈值与能量块峭度的映射关系:由低频能量块的峭度来计算阈值,具体方法如下:
1)生成指纹数据:在信噪比分别为{5dB,6dB,…,25dB}时进行多次仿真,分别计算每一次仿真能量块的峭度并进行取整,同时归一化阈值取{0.1,0.2,…,1.0}之间不同值时,计算测距误差。峭度的计算公式如下:
其中,Nb表示能量块的个数,表示能量块的均值,δ表示能量块的标准差。
2)选取最优阈值:当能量块峭度取不同值时,分别计算取不同归一化阈值时测距误差的平均值,取平均误差最小时的归一化阈值作为此峭度下选取的归一化阈值。
3)曲线拟合:将能量块峭度作为横坐标,将最佳归一化阈值作为纵坐标,采用最小二乘法,建立能量块峭度k和最佳归一化阈值ξ的映射关系,如下所示:
CM1信道:ξ=-3.1079*10-8k3+2.5191*10-5k 2-7.0462*10-3*k+0.82585
CM2信道:ξ=-2.6988*10-8k3+2.2681*10-5k 2-6.5214*10-3*k+0.7816
E.测距计算:收到UWB信号后,在进行实际测距时,执行如下步骤:
1)构造能量块:对接收的UWB信号进行小波变换,取低频系数构造能量块;
2)计算阈值:计算能量块峭度,根据建立的映射关系计算出阈值;
3)计算UWB信号到达时间:将能量块与阈值进行比较,取第一个超过阈值的能量块作为信号到达的时间TOA,与传统的能量接收方法的对比分别如图1、图2所示。
4)计算距离,计算方法如下:
其中,c表示电磁波在空气中的传播速度,表示估计的信号到达时间TOA。
F.定位计算:多个定位基站将上述测距结果通过网络传输到服务器,服务器再根据定位标签与不同定位基站之间的距离差,利用传统的双曲线无线定位算法,实现海上平台工作人员的定位。
Claims (1)
1.一种基于小波分解低频系数的海上作业平台人员定位方法其特征在于包括以下步骤:
A.UWB信号发送:安装在作业平台工作人员身上的定位标签发送UWB信号,定位基站接收该信号,如果要实现三维定位则基站数量不能少于4个,如果要实现二维定位则基站数量不能少于3个;
B.UWB信号小波分解:定位基站对接收到的UWB信号进行小波分解,计算方法为:
其中,c0表示接收到的UWB信号,j∈{0,1,2,…}表示小波分解的层数,cj表示UWB信号经小波分解到j层上得到的低频部分的系数,dj表示UWB信号经小波分解到j层上得到的高频部分的系数,h(n)表示低通滤波器,g(n)表示高通滤波器,其通过选取的小波类型决定;
C.取低频系数构造能量块:每隔一定的区间对低频系数进行积分计算,构造一系列能量块,方法为:
其中,Tf表示UWB信号一个帧的长度,Tb表示积分周期即每个能量块的长度,r(t)表示接收到的UWB信号,z[n]表示能量块;
D.建立阈值与能量块峭度的映射关系:由低频能量块的峭度来计算阈值,具体方法如下:
1)生成指纹数据:在不同的信噪比下进行多次仿真,分别计算每一次仿真能量块的峭度并进行取整,同时归一化阈值取{0.1,0.2,…,1.0}之间不同值时,计算测距误差,峭度的计算公式如下:
其中,Nb表示能量块的个数,表示能量块的均值,δ表示能量块的标准差;
2)选取最优阈值:当能量块峭度取不同值时,分别计算取不同归一化阈值时测距误差的平均值,取平均误差最小时的归一化阈值作为此峭度下选取的归一化阈值;
3)曲线拟合:根据能量块峭度和归一化阈值的对应关系,利用曲线拟合建立起峭度和归一化阈值的映射关系;
E.测距计算:收到UWB信号后,在进行实际测距时,执行如下步骤:
1)构造能量块:对接收到的UWB信号进行小波变换,取低频系数构造能量块;
2)计算阈值:计算能量块峭度,根据建立的映射关系计算出阈值;
3)计算UWB信号到达时间:将能量块与阈值进行比较,取第一个超过阈值的能量块作为信号到达的时间TOA;
4)计算距离,计算方法如下:
其中,c表示电磁波在空气中的传播速度,表示估计的信号到达时间TOA;
F.定位计算:多个定位基站将上述测距结果通过网络传输到服务器,服务器再根据定位标签与不同定位基站之间的距离差,利用传统的双曲线无线定位算法,实现海上作业平台工作人员的定位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710177034.1A CN106959432B (zh) | 2017-03-23 | 2017-03-23 | 一种基于小波分解低频系数的海上作业平台人员定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710177034.1A CN106959432B (zh) | 2017-03-23 | 2017-03-23 | 一种基于小波分解低频系数的海上作业平台人员定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106959432A true CN106959432A (zh) | 2017-07-18 |
CN106959432B CN106959432B (zh) | 2019-06-18 |
Family
ID=59471015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710177034.1A Expired - Fee Related CN106959432B (zh) | 2017-03-23 | 2017-03-23 | 一种基于小波分解低频系数的海上作业平台人员定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106959432B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108092736A (zh) * | 2017-12-26 | 2018-05-29 | 广州海格通信集团股份有限公司 | 一种基于小波变换的低检测概率波形方法 |
CN108169707A (zh) * | 2017-11-30 | 2018-06-15 | 中原智慧城市设计研究院有限公司 | 基于超宽带能量谱特征向量的被动式定位算法 |
CN111890126A (zh) * | 2020-07-03 | 2020-11-06 | 太原理工大学 | 一种基于声压能量峭度指标的早期车削颤振预警监测方法 |
CN113179549A (zh) * | 2021-04-25 | 2021-07-27 | 深圳大漠大智控技术有限公司 | 一种低信噪比下获取基站与标签距离的方法及其相关组件 |
CN114222241A (zh) * | 2021-10-29 | 2022-03-22 | 中国石油大学(华东) | 一种基于小波变换模极大值的无线定位方法 |
CN114697860A (zh) * | 2022-03-28 | 2022-07-01 | 广东电网有限责任公司广州供电局 | 一种基于神经网络的多功能定位装置 |
CN115293219A (zh) * | 2022-09-29 | 2022-11-04 | 中国电建集团华东勘测设计研究院有限公司 | 一种融合小波和峭度的脉冲信号去噪方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011257158A (ja) * | 2010-06-04 | 2011-12-22 | Toyota Motor Corp | レーダ装置 |
CN104219762A (zh) * | 2014-10-09 | 2014-12-17 | 中国石油大学(华东) | 一种基于均方差的超宽带无线定位方法 |
CN105589063A (zh) * | 2016-01-12 | 2016-05-18 | 中国海洋大学 | 基于偏度的脉冲无线电60GHz测距方法 |
WO2018025794A1 (ja) * | 2016-08-04 | 2018-02-08 | シャープ株式会社 | 基地局装置、ロケーションサーバーおよび通信方法 |
-
2017
- 2017-03-23 CN CN201710177034.1A patent/CN106959432B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011257158A (ja) * | 2010-06-04 | 2011-12-22 | Toyota Motor Corp | レーダ装置 |
CN104219762A (zh) * | 2014-10-09 | 2014-12-17 | 中国石油大学(华东) | 一种基于均方差的超宽带无线定位方法 |
CN105589063A (zh) * | 2016-01-12 | 2016-05-18 | 中国海洋大学 | 基于偏度的脉冲无线电60GHz测距方法 |
WO2018025794A1 (ja) * | 2016-08-04 | 2018-02-08 | シャープ株式会社 | 基地局装置、ロケーションサーバーおよび通信方法 |
Non-Patent Citations (2)
Title |
---|
李娟等: ""UWB接收信号积分能量块的统计特性分析"", 《研究与设计》 * |
翟双等: ""IR-UWB系统中基于归一化门限的TOA估计"", 《北京邮电大学学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169707A (zh) * | 2017-11-30 | 2018-06-15 | 中原智慧城市设计研究院有限公司 | 基于超宽带能量谱特征向量的被动式定位算法 |
CN108092736A (zh) * | 2017-12-26 | 2018-05-29 | 广州海格通信集团股份有限公司 | 一种基于小波变换的低检测概率波形方法 |
CN111890126A (zh) * | 2020-07-03 | 2020-11-06 | 太原理工大学 | 一种基于声压能量峭度指标的早期车削颤振预警监测方法 |
CN111890126B (zh) * | 2020-07-03 | 2022-03-11 | 太原理工大学 | 一种基于声压能量峭度指标的早期车削颤振预警监测方法 |
CN113179549A (zh) * | 2021-04-25 | 2021-07-27 | 深圳大漠大智控技术有限公司 | 一种低信噪比下获取基站与标签距离的方法及其相关组件 |
CN114222241A (zh) * | 2021-10-29 | 2022-03-22 | 中国石油大学(华东) | 一种基于小波变换模极大值的无线定位方法 |
CN114697860A (zh) * | 2022-03-28 | 2022-07-01 | 广东电网有限责任公司广州供电局 | 一种基于神经网络的多功能定位装置 |
CN115293219A (zh) * | 2022-09-29 | 2022-11-04 | 中国电建集团华东勘测设计研究院有限公司 | 一种融合小波和峭度的脉冲信号去噪方法 |
CN115293219B (zh) * | 2022-09-29 | 2023-04-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种融合小波和峭度的脉冲信号去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106959432B (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106959432A (zh) | 一种基于小波分解低频系数的海上作业平台人员定位方法 | |
CN102638888B (zh) | 基于信号统计的室内定位指纹分组方法 | |
CN104076349B (zh) | 一种基于多普勒频移的被动式移动目标定位方法 | |
CN103795479B (zh) | 一种基于特征值的协作频谱感知方法 | |
CN102064895B (zh) | Rssi与模式匹配相结合的无源定位方法 | |
CN106231549B (zh) | 一种基于重构算法的60GHz脉冲室内定位方法 | |
CN105960017B (zh) | 基于超宽带节点网络的免携带设备定位方法 | |
CN104020444B (zh) | 一种双轮差动式机器人室内组合定位方法 | |
CN105866750A (zh) | 一种多个gps卫星微弱回波信号检测的方法 | |
Jiang et al. | Ultra-wide band applications in industry: a critical review | |
CN103874020A (zh) | 非直达径环境下单接收器的超宽带定位方法 | |
CN104270190A (zh) | 基于电离层资料的同步自适应短波通信选频方法 | |
CN106255059A (zh) | 一种基于几何方式的无设备目标定位方法 | |
Zhu et al. | EKF localization based on TDOA/RSS in underground mines using UWB ranging | |
CN104602191B (zh) | 一种海上搜救无线传感网节点定位方法 | |
Liu et al. | Research on wavelet threshold denoising method for UWB tunnel personnel motion location | |
Li et al. | Shadowing correlation and a novel statistical model for inland river radio channel | |
CN102832966B (zh) | 基于非视距抑制的矿井超宽带定位方法及系统 | |
CN105891800B (zh) | 单发三收天线雷达抗多径散射检测方法 | |
CN107202973A (zh) | 基于超宽带能量谱特征向量定位人员的方法 | |
Li et al. | Measurement‐based wireless channel analysis and modelling for shipping environments | |
CN113203985B (zh) | 一种短波同频信号直接定位方法 | |
CN101984562B (zh) | 一种窄带信道增益估计方法 | |
Lee et al. | Stand-off distances for non-line-of-sight maritime mobile applications in 5 GHz band | |
Lu et al. | Enhanced LoRaWAN RSSI indoor localization based on BP neural network |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190618 Termination date: 20200323 |