CN107748375A - 一种基于ceemd‑ht算法的gnss实时多径抑制算法 - Google Patents
一种基于ceemd‑ht算法的gnss实时多径抑制算法 Download PDFInfo
- Publication number
- CN107748375A CN107748375A CN201710967811.2A CN201710967811A CN107748375A CN 107748375 A CN107748375 A CN 107748375A CN 201710967811 A CN201710967811 A CN 201710967811A CN 107748375 A CN107748375 A CN 107748375A
- Authority
- CN
- China
- Prior art keywords
- signal
- mode
- multipath
- frequency
- mode decomposition
- 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.)
- Pending
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/22—Multipath-related issues
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及滤除GNSS信号中多径信号的方法,具体为一种基于CEEMD‑HT算法的GNSS实时多径抑制算法。本发明在GNSS接收机跟踪环路中应用了模态分解辅助载波跟踪的方法。首先,对接收到的直达和多径混叠的信号进行模态分解,模态分解可将杂乱的混合信号分解为规律的不同频率和幅度的多个信号;其次,根据鉴相器的结果得出包含多普勒频移的信号频率;然后将模态分解得到的频率与鉴相器结果比较,选择最接近的频率进入跟踪环路。本发明在载波层面上对接收信号进行处理,可对短多径起到抑制作用。对于动态环境,直达信号相对多径信号是有规律性的,模态分解可将随机变化的多径信号进行有效的滤除。因此,本技术方案可达到对动态多径进行实时抑制的目的。
Description
技术领域
本发明涉及滤除GNSS信号中多径信号的方法,具体为一种基于CEEMD-HT算法的GNSS实时多径抑制算法,来提高GNSS接收机的抗多径性能,可适用于高动态环境。
背景技术
多径效应是全球卫星导航系统(GNSS,Global Navigation Satellite System)的主要误差源,其引起的测量误差严重影响接收机的定位精度。目前,对多径效应的抑制主要从对天线的改进和对接收信号的抗多径处理两大方面进行。国内外对多径抑制技术的研究也取得了很多成果。Van Dierendonck提出的窄相关技术可有效抑制长多径,但是对短多径无法进行抑制。Townsend B提出的最大似然估计算法可对多径信号进行估计,一定程度上减少了多径干扰,但该算法计算复杂度高。研究者们还提出了基于天线阵列、能量算子、粒子滤波、自适应滤波、盲源分离、地图匹配等方法。虽然这些方法可在一定程度上抑制多径效应,从而提高GNSS接收机的定位精度,但普遍存在复杂度高以及无法对高动态多径进行实时抑制的问题。
发明内容
本发明在GNSS接收机跟踪环路中应用了模态分解辅助载波跟踪的方法。首先,对接收到的直达和多径混叠的信号进行模态分解,模态分解可将杂乱的混合信号分解为规律的不同频率和幅度的多个信号;其次,根据鉴相器的结果得出包含多普勒频移的载波频率;然后将模态分解得到的频率与鉴相器结果比较,选择频率最接近的模态作为近似直达信号进入跟踪环路。本发明在载波层面上对接收信号进行处理,可对短多径起到抑制作用。对于动态环境,直达信号相对多径信号是有规律性的,模态分解可将随机变化的多径信号进行有效的滤除,筛选出近似直达信号进入码环相关器,即可得出高精度的定位结果。因此,本技术方案可达到对动态多径进行实时抑制的目的。
本发明是采用如下的技术方案实现的:一种基于CEEMD-HT算法的GNSS实时多径抑制算法,包括以下步骤:
(1)接收机天线接收到的信号首先经射频前端变成中频信号r(t);
(2)中频信号r(t)采用互补集总经验模态分解方法分解:
<1>中频信号r(t)添加大小相同方向相反的白噪声w(t)得到r+(t)和
<2>r+(t)和r-(t)分别采用经验模态分解方法进行模态分解,分别得到各阶模态Ci+(t)和Ci-(t),
<3>加入N对白噪声,对分解出的各阶模态进行统计平均,分别得到各模态的统计平均值各模态的统计平均值再统计平均,得到中频信号r(t)的估计信号式中Cr(t)为估计信号中的各阶模态,rm(t)为第m阶残余;
(3)对估计信号中各模态进行希尔伯特变换,得到估计信号瞬时值其中 是瞬时幅值,是瞬时相位,是瞬时角频率,P为柯西主值;
(4)将估计信号瞬时值Ar(t)与跟踪环路码环中得到的多普勒载波频率比较,筛选出与多普勒载波频率最接近的模态信号作为近似直达信号,近似直达信号进入载波环进行定位运算。
本发明针对多径干扰对GNSS接收机的影响,利用直达信号的规律性和多径信号的随机性,使用一种基于CEEMD-HT算法的多径实时抑制方法,可以有效地减少多径信号的影响,提高接收机的定位精度。
附图说明
图1为本发明算法框图。
图2为CEEMD-HT算法流程图。
图3为跟踪环路示意图。
图4为频率选择流程图。
具体实施方式
一种基于CEEMD-HT算法的GNSS实时多径抑制算法,包括以下步骤:
(1)接收机天线接收到的信号首先经射频前端变成中频信号r(t)。
(2)中频信号r(t)采用互补集总经验模态分解方法分解:
<1>中频信号r(t)添加大小相同方向相反的白噪声w(t)得到r+(t)和
<2>r+(t)和r-(t)分别采用经验模态分解方法进行模态分解,分别得到各阶模态Ci+(t)和Ci-(t),
对于信号r+(t)采用经验模态分解EMD的步骤如下:
(1)用三次样条法绘制出由信号极大值决定的上包络和极小值决定的下包络;
(2)对上下包络求均值得到m(t);
(3)提取信号h(t),h(t)=r+(t)-m(t);
(4)重复以上步骤直至m(t)接近零,此时的h(t)被指定为一个模态C(t);
(5)计算残余
(6)对每个ri(t)进行以上步骤得到每个模态和残余,最终得到
其中Ci+(t)是各阶模态,rn(t)是第n阶残余;
对于信号r-(t)采用经验模态分解EMD的步骤和信号r+(t)采用经验模态分解EMD的步骤相同。
<3>加入N对白噪声,对分解出的各阶模态进行统计平均,分别得到各模态的统计平均值各模态的统计平均值再统计平均,得到中频信号r(t)的估计信号式中Cr(t)为估计信号中的各阶模态,rm(t)为第m阶残余。
(3)对估计信号中各模态进行希尔伯特变换(HT),得到估计信号瞬时值其中 是瞬时幅值,是瞬时相位,是瞬时角频率,P为柯西主值;
(4)将估计信号瞬时值Ar(t)与跟踪环路码环中得到的多普勒载波频率fd比较,筛选出与多普勒载波频率最接近的模态信号作为近似直达信号,近似直达信号进入载波环进行定位运算。
经验模态分解方法EMD(Empirical Mode Decomposition)是美国宇航局在20世纪末提出的,被认为是对傅里叶变换的重大突破,但是该方法分解出的模态中存在模态混叠现象。于是,Yeh等人在EMD的基础上进行了改进提出了互补集总经验模态分解方法CEEMD(Complementary Ensemble Empirical Mode Decomposition)。该方法通过添加成对的大小相同方向相反的高斯白噪声可消除模态混叠现象。对CEEMD分解出的模态进行希尔伯特变换HT(Hilbert Transform)之后可得到信号瞬时频率、瞬时幅度和瞬时相位,方便对信号进行处理。CEEMD-HT算法流程如图2所示。
Claims (1)
1.一种基于CEEMD-HT算法的GNSS实时多径抑制算法,其特征在于包括以下步骤:
(1)接收机天线接收到的信号首先经射频前端变成中频信号r(t);
(2)中频信号r(t)采用互补集总经验模态分解方法分解:
<1>中频信号r(t)添加大小相同方向相反的白噪声w(t)得到r+(t)和
<2>r+(t)和r-(t)分别采用经验模态分解方法进行模态分解,分别得到各阶模态Ci+(t)和
<3>加入N对白噪声,对分解出的各阶模态进行统计平均,分别得到各模态的统计平均值各模态的统计平均值再统计平均,得到中频信号r(t)的估计信号式中Cr(t)为估计信号中的各阶模态,rm(t)为第m阶残余;
(3)对估计信号中各模态进行希尔伯特变换,得到估计信号瞬时值其中 是瞬时幅值,是瞬时相位,是瞬时角频率,P为柯西主值;
(4)将估计信号瞬时值Ar(t)与跟踪环路码环中得到的多普勒载波频率比较,筛选出与多普勒载波频率最接近的模态信号作为近似直达信号,近似直达信号进入载波环进行定位运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710967811.2A CN107748375A (zh) | 2017-10-18 | 2017-10-18 | 一种基于ceemd‑ht算法的gnss实时多径抑制算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710967811.2A CN107748375A (zh) | 2017-10-18 | 2017-10-18 | 一种基于ceemd‑ht算法的gnss实时多径抑制算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107748375A true CN107748375A (zh) | 2018-03-02 |
Family
ID=61252614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710967811.2A Pending CN107748375A (zh) | 2017-10-18 | 2017-10-18 | 一种基于ceemd‑ht算法的gnss实时多径抑制算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107748375A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594177A (zh) * | 2018-03-16 | 2018-09-28 | 西安电子科技大学 | 基于改进hht的雷达信号调制方式分析方法、信号处理系统 |
CN108776326A (zh) * | 2018-04-13 | 2018-11-09 | 北京邮电大学 | 一种多径效应抑制方法、装置及设备 |
CN109001768A (zh) * | 2018-07-31 | 2018-12-14 | 太原理工大学 | 一种应用于天线中的改进双极化顺序ml多径抑制方法 |
CN112180408A (zh) * | 2020-09-29 | 2021-01-05 | 中山大学 | 一种基于智能终端的多径误差提取方法和相关装置 |
CN113189624A (zh) * | 2021-04-30 | 2021-07-30 | 中山大学 | 一种自适应分类的多径误差提取方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1657567A2 (en) * | 2004-11-12 | 2006-05-17 | Northrop Grumman Corporation | Multipath height finding method |
CN105891815A (zh) * | 2016-05-03 | 2016-08-24 | 西南大学 | 一种基于广播信号无源定位的联合估计算法 |
CN106556647A (zh) * | 2016-11-18 | 2017-04-05 | 湘潭大学 | 一种冲击回波数据处理方法 |
CN106650218A (zh) * | 2016-10-19 | 2017-05-10 | 上海电机学院 | 基于ceemd算法和希尔伯特变换的谐波分析方法 |
-
2017
- 2017-10-18 CN CN201710967811.2A patent/CN107748375A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1657567A2 (en) * | 2004-11-12 | 2006-05-17 | Northrop Grumman Corporation | Multipath height finding method |
CN105891815A (zh) * | 2016-05-03 | 2016-08-24 | 西南大学 | 一种基于广播信号无源定位的联合估计算法 |
CN106650218A (zh) * | 2016-10-19 | 2017-05-10 | 上海电机学院 | 基于ceemd算法和希尔伯特变换的谐波分析方法 |
CN106556647A (zh) * | 2016-11-18 | 2017-04-05 | 湘潭大学 | 一种冲击回波数据处理方法 |
Non-Patent Citations (2)
Title |
---|
LI JING 等: "GPR Signal Denoising and Target Extraction", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
SU YONGQING 等: "The Application of Empirical Mode Decomposition in Data Detrending for Ionosphere Amplitude Scintillation Study", 《INTERNATIONAL CONFERENCE ON CONTROL ENGINEERING AND AUTOMATION (ICCEA 2014)》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594177A (zh) * | 2018-03-16 | 2018-09-28 | 西安电子科技大学 | 基于改进hht的雷达信号调制方式分析方法、信号处理系统 |
CN108776326A (zh) * | 2018-04-13 | 2018-11-09 | 北京邮电大学 | 一种多径效应抑制方法、装置及设备 |
CN108776326B (zh) * | 2018-04-13 | 2022-04-05 | 北京邮电大学 | 一种多径效应抑制方法、装置及设备 |
CN109001768A (zh) * | 2018-07-31 | 2018-12-14 | 太原理工大学 | 一种应用于天线中的改进双极化顺序ml多径抑制方法 |
CN109001768B (zh) * | 2018-07-31 | 2022-09-13 | 太原理工大学 | 一种应用于天线中的改进双极化顺序ml多径抑制方法 |
CN112180408A (zh) * | 2020-09-29 | 2021-01-05 | 中山大学 | 一种基于智能终端的多径误差提取方法和相关装置 |
CN112180408B (zh) * | 2020-09-29 | 2023-06-23 | 中山大学 | 一种基于智能终端的多径误差提取方法和相关装置 |
CN113189624A (zh) * | 2021-04-30 | 2021-07-30 | 中山大学 | 一种自适应分类的多径误差提取方法及装置 |
CN113189624B (zh) * | 2021-04-30 | 2023-10-03 | 中山大学 | 一种自适应分类的多径误差提取方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107748375A (zh) | 一种基于ceemd‑ht算法的gnss实时多径抑制算法 | |
CN102162852B (zh) | 大规模频率偏移下微弱gnss信号捕获方法及其系统 | |
CN105866750B (zh) | 一种多个gps卫星微弱回波信号检测的方法 | |
CN106646541A (zh) | 一种基于差分相关积分的北斗弱信号捕获方法 | |
WO2015172622A1 (zh) | 一种高频地波雷达射频干扰抑制的方法 | |
CN103728634B (zh) | 双天线a-gnss接收机系统及方法 | |
US8417207B2 (en) | High-performance cellular telephone receiver | |
CN105158773A (zh) | 一种用于卫星导航接收系统的多波束指向抗干扰处理方法 | |
CN102590829A (zh) | 用于卫星导航系统的完全自适应陷波器及其陷波方法 | |
CN107991695A (zh) | 基于补零算法和差分相干算法的北斗弱信号捕获方法 | |
CN105629266B (zh) | 卫星导航欺骗式和压制式干扰盲自适应的联合抑制方法 | |
CN103281275B (zh) | 一种msk/gmsk直接序列扩频信号接收机 | |
CN105704081B (zh) | 多个偏移载波调幅信号的优选和解调方法 | |
CA2776191A1 (en) | Method and system for locating interferences affecting a satellite-based radionavigation signal | |
CN104360355A (zh) | 抗干扰方法和装置 | |
CN105388499A (zh) | 伪卫星地基增强系统中利用双天线克服远近效应的方法 | |
CN105372678B (zh) | 一种正弦boc调制信号的无模糊跟踪方法 | |
CN104280749B (zh) | 基于soc的北斗导航系统接收机抗攻击系统及其抗攻击方法 | |
CN105242286A (zh) | 一种基于小波域滤波码捕获方法 | |
US9124356B2 (en) | Reconfigurable correlator for a navigation system | |
CN106990393B (zh) | 一种基于盲波束形成的导航接收机抗干扰系统 | |
CN114217329A (zh) | 一种基于串行搜索的短码捕获方法 | |
CN114578387A (zh) | 多径干扰消除方法及装置、估计方法及装置、接收机 | |
CN102508266B (zh) | 一种抑制多径干扰及噪声的滤波器方法 | |
CN116338740A (zh) | 一种基于pcf重构的北斗b1c信号无模糊捕获计算方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180302 |