CN110764120B - 一种高灵敏卫星导航信号捕获方法 - Google Patents
一种高灵敏卫星导航信号捕获方法 Download PDFInfo
- Publication number
- CN110764120B CN110764120B CN201911082548.4A CN201911082548A CN110764120B CN 110764120 B CN110764120 B CN 110764120B CN 201911082548 A CN201911082548 A CN 201911082548A CN 110764120 B CN110764120 B CN 110764120B
- Authority
- CN
- China
- Prior art keywords
- signal
- node
- fft
- local
- priority queue
- 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
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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/246—Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radio Relay Systems (AREA)
Abstract
本发明涉及一种高灵敏卫星导航信号捕获方法,对接收信号进行滤波、下变频、数字采样后,再与中频载波混频,然后抽取滤波处理、混频、分组后构成二叉树,对二叉树进行寻优,对最优累加结果进行FFT后和各卫星PRN码的FFT值的共轭值进行相乘再进行IFFT处理,选出最大值即为多普勒、码相位。本发明通过对二叉树进行寻优,利用比特信息,可以实现不同场景下微弱卫星信号的快速捕获。
Description
技术领域
本发明涉及卫星导航技术领域,具体涉及一种基于接收机的弱信号处理方法。
背景技术
在城市峡谷,森林、室内等复杂环境下,GNSS信号由于受到遮挡,会产生严重衰减,且存在多径和干扰,信号功率将低于正常GPS接收机的工作范围,其可用性及定位精度都会大大降低。高灵敏度GNSS接收技术已成为导航技术研究热点。
高灵敏度技术主要是通过相干积分和非相干积分提高信号的处理增益,将GNSS信号从背景噪声中提取出来。相干积分时间长短的选取受到数据比特翻转和多普勒频差的限制,如果可以得到数据位值和边沿位置及实时多普勒频差进行码多普勒补偿,则可以进行长时间的相干积分,提高信号增益。
目前对弱信号捕获,主要包括相干检测,非相干检测,差分检测。相干检测捕获算法在高斯白噪声的条件下可以提供最大的增益。但是,最大相干积分时间会受到晶振稳定性,比特翻转,多普勒等因素影响;非相干累积对导航数据和多普勒频率引起相位变化不敏感,但也存在“平方损失”,随着C/N0的降低,平方损失越严重。
Psiaki提出了半比特交替捕获技术。将收到的20ms数据块分成两段,保证其中一段的数据完全是在同一比特,分别对两段数据进行相干积分,然后,相干积分的结果进一步进行非相干累积。Ziedan给出一种基于导航数据预测的弱信号捕获算法,但是数据预测的准确性和算法的复杂度仍然无法有效地解决。
除了导航数据翻转外,长相干积分时间导致频率槽搜索增多,导致频率域搜索次数显著增加,限制了相干积分时间。
发明内容
要解决的技术问题
为了克服现有卫星弱信号捕获技术的缺点与不足,本发明提供了一种高灵敏卫星导航信号捕获方法,该方法适用于美国的GPS及我国的“北斗”系统,可以实现不同场景下微弱卫星信号的快速捕获。
技术方案
一种高灵敏卫星导航信号捕获方法,其特征在于步骤如下:
步骤1:将接收到的射频信号经过滤波、下变频、数字采样后,得到基带信号signal_rec,signal_rec和本地中频载波local_if混频后得到复信号signal_if;
步骤2:对signal_if进行Q倍CIC抽取滤波处理,将信号的采样率由中频采样率fs降为fs′,得到signal_dif;Q倍大于2倍带宽;
步骤3:signal_dif与local_freq_bin进行混频得到data,所述的local_freq_bin为多普勒补偿;
步骤4:对data进行分组,每fs′/1000个点为一组,共N组,构成深度为floor(N/20)+1的二叉树,内部各节点度为2,代表两个不同的状态路径,累积度量值为acc(k);
步骤5:对二叉树进行寻优,首先进行初始化,优先队列包含开始节点n0,初始化累积度量值acc(n0)=0;
步骤6:寻优器不断的将优先队列中具有最大累积度量值的节点nk输出,这个节点在传统Viterbi网格图中具有最佳的路径;
步骤7:将寻优器输出节点nk插入优先队列,插入的内容包括nk所具有的时刻、状态以及前向指针;如果该节点不是终端节点,还要将nk的后继节点nk+1插入到优先队列结构中,插入节点nk+1的累积度量为acc(nk+1)=acc(nk)+d(nk,k+1)ak,将节点插入优先队列后,返回步骤6;如果网格图已经扩展到了终端节点,则停止寻优;
步骤8:利用优先队列对数据进行累加,得到数据b1,k,b2,k,…bn,k;累加方式如下所示:
设置b1,0,b2,0,…bn,k初值为0,
其中k=1,2....m,dk∈{-1,1};
步骤9:对最优累加结果进行FFT得到rec_fft;
步骤10:利用本地码表提前生成各卫星PRN码的FFT值,并对其进行共轭处理,得到信号local_fft;
步骤11:对步骤9和10得到的值进行相乘得到rslt_fft;
步骤12:对rslt_fft进行IFFT处理,得到结果,将其存入acqu_mat,返回步骤3,直至遍历所有local_freq_bin,得到acqu_mat;
步骤13:比较acqu_mat,得到最大值,其对应值为多普勒、码相位。
有益效果
本发明提出的一种高灵敏卫星导航信号捕获方法,通过对二叉树进行寻优,利用比特信息,可以实现不同场景下微弱卫星信号的快速捕获。
附图说明
图1为本发明捕获信号流程图;
图2为本发明数据分组图;
图3位本发明二叉树结构图。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
下面结合附图并通过具体实施方式以GPS信号为例来进一步说明本发明的技术方案。
(1)将接收到99ms的射频信号经过滤波、下变频、数字采样后,得到基带信号signal_rec,signal_rec和本地中频载波local_if混频后得到复信号signal_if;
(2)对signal_if进行CIC抽取滤波处理,将信号的采样率由中频采样率16.384M降为2.046M,得到signal_dif;
(3)signal_dif与local_freq_bin进行混频,local_freq_bin为1KHz范围内的精细多普勒补偿,得到data;
(4)对data进行分组,每2046个点为一组,共99组,起始点为k,k=1,2...19每20组为一块,构成深度为5的二叉树,内部各节点度为2,代表两个不同的状态路径[1,-1],路径累积度量值为acc(k);
(5)对二叉树进行寻优,首先进行初始化,优先队列包含开始节点n0,初始化累积度量值acc(n0)=0;
(6)寻优器不断的将优先队列中具有最大累积度量值的节点nk输出,这个节点在传统Viterbi网格图中具有最佳的路径。
(7)将寻优器输出节点nk插入优先队列,插入的内容包括nk所具有的时刻、状态以及前向指针。如果该节点不是终端节点,还要将nk的后继节点nk+1插入到优先队列结构中,插入节点nk+1的累积度量为acc(nk+1)=acc(nk)+d(nk,k+1)ak,将节点插入优先队列后,返回步骤(6);如果网格图已经扩展到了最后底层节点,则停止寻优。
(8)利用优先队列对数据进行累加,得到数据b1,k,b2,k,…bn,k。累加方式如下所示:
设置b1,0,b2,0,…bn,k初值为0,
其中k=1,2....m,dk∈{-1,1};
得到累加值rec_data;
(9)对最优累加结果rec_data进行FFT得到rec_fft。
(10)利用本地码表提前生成各卫星PRN码的FFT值,并对其进行共轭处理,得到信号local_fft。
(11)对步骤(9),(10)得到的值进行相乘得到rslt_fft。
(12)对rslt_fft进行IFFT处理,得到结果,将其存入acqu_mat(:,k),返回步骤(3),直至遍历所有local_freq_bin,得到acqu_mat,返回步骤(4),k=19时,转至步骤(13)。
(13)比较acqu_mat,得到最大值,其对应值为捕获多普勒,码相位。
Claims (1)
1.一种高灵敏卫星导航信号捕获方法,其特征在于步骤如下:
步骤1:将接收到的射频信号经过滤波、下变频、数字采样后,得到基带信号signal_rec,signal_rec和本地中频载波local_if混频后得到复信号signal_if;
步骤2:对signal_if进行Q倍CIC抽取滤波处理,将信号的采样率由中频采样率fs降为fs′,得到signal_dif;Q倍大于2倍带宽;
步骤3:signal_dif与local_freq_bin进行混频得到data,所述的local_freq_bin为多普勒补偿;
步骤4:对data进行分组,每fs′/1000个点为一组,共N组,构成深度为floor(N/20)+1的二叉树,内部各节点度为2,代表两个不同的状态路径,累积度量值为acc(k);
步骤5:对二叉树进行寻优,首先进行初始化,优先队列包含开始节点n0,初始化累积度量值acc(n0)=0;
步骤6:寻优器不断的将优先队列中具有最大累积度量值的节点nk输出,这个节点在传统Viterbi网格图中具有最佳的路径;
步骤7:将寻优器输出节点nk插入优先队列,插入的内容包括nk所具有的时刻、状态以及前向指针;如果该节点不是终端节点,还要将nk的后继节点nk+1插入到优先队列结构中,插入节点nk+1的累积度量为acc(nk+1)=acc(nk)+d(nk,k+1)ak,将节点插入优先队列后,返回步骤6;如果网格图已经扩展到了终端节点,则停止寻优;
步骤8:利用优先队列对数据进行累加,得到数据b1,k,b2,k,…bn,k;累加方式如下所示:
设置b1,0,b2,0,…bn,k初值为0,
其中k=1,2....m,dk∈{-1,1};
步骤9:对最优累加结果进行FFT得到rec_fft;
步骤10:利用本地码表提前生成各卫星PRN码的FFT值,并对其进行共轭处理,得到信号local_fft;
步骤11:对步骤9和10得到的值进行相乘得到rslt_fft;
步骤12:对rslt_fft进行IFFT处理,得到结果,将其存入acqu_mat,返回步骤3,直至遍历所有local_freq_bin,得到acqu_mat;
步骤13:比较acqu_mat,得到最大值,其对应值为多普勒、码相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911082548.4A CN110764120B (zh) | 2019-11-07 | 2019-11-07 | 一种高灵敏卫星导航信号捕获方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911082548.4A CN110764120B (zh) | 2019-11-07 | 2019-11-07 | 一种高灵敏卫星导航信号捕获方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110764120A CN110764120A (zh) | 2020-02-07 |
CN110764120B true CN110764120B (zh) | 2023-03-28 |
Family
ID=69336476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911082548.4A Active CN110764120B (zh) | 2019-11-07 | 2019-11-07 | 一种高灵敏卫星导航信号捕获方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110764120B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276374A (zh) * | 2007-03-30 | 2008-10-01 | 索尼株式会社 | 内容管理装置、图像显示装置、图像拾取装置及处理方法 |
CN102608626A (zh) * | 2012-03-09 | 2012-07-25 | 暨南大学 | 一种高灵敏度卫星导航信号捕获方法及系统 |
CN103278829A (zh) * | 2013-05-06 | 2013-09-04 | 东南大学 | 一种基于gpu的并行导航卫星信号跟踪方法及其系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090303251A1 (en) * | 2008-06-10 | 2009-12-10 | Andras Balogh | Displaying, processing and storing geo-located information |
US8442173B2 (en) * | 2010-02-09 | 2013-05-14 | Analog Devices, Inc. | Apparatus and method for clock and data recovery |
US20170367684A1 (en) * | 2015-01-05 | 2017-12-28 | Innomind Technology Corporation | Systems and methods for super-resolution compact ultrasound imaging |
-
2019
- 2019-11-07 CN CN201911082548.4A patent/CN110764120B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276374A (zh) * | 2007-03-30 | 2008-10-01 | 索尼株式会社 | 内容管理装置、图像显示装置、图像拾取装置及处理方法 |
CN102608626A (zh) * | 2012-03-09 | 2012-07-25 | 暨南大学 | 一种高灵敏度卫星导航信号捕获方法及系统 |
CN103278829A (zh) * | 2013-05-06 | 2013-09-04 | 东南大学 | 一种基于gpu的并行导航卫星信号跟踪方法及其系统 |
Non-Patent Citations (3)
Title |
---|
Data fusion of high-resolution satellite imagery and LiDAR data for automatic building extraction;Sohn G 等;《ISPRS Journal of Photogrammetry and Remote Sensing》;20071231;第62卷(第01期);第43-63页 * |
GPS导航信号的质量评估研究;付娟;《中国优秀硕士学位论文全文数据库信息科技辑》;20190215(第02期);第I136-1506页 * |
卫星导航信号捕获的一种高效相关算法;戴志春 等;《全球定位系统》;20181231;第43卷(第06期);第33-36页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110764120A (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102162852B (zh) | 大规模频率偏移下微弱gnss信号捕获方法及其系统 | |
CN101029923B (zh) | 全球定位系统室内定位与导航技术的软件接收机 | |
CN1122380C (zh) | 用于确定频率同步信号的位置的方法、设备与系统 | |
CN104536016B (zh) | 一种gnss新体制信号捕获装置及方法 | |
CN102928853B (zh) | 一种捕获北斗d1卫星导航系统弱信号的方法 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
CN102486539B (zh) | 一种改进的导航卫星信号跟踪方法 | |
CN101520505B (zh) | 全球卫星导航系统弱信号下的自适应相干累积联合捕获方法 | |
CN104765052B (zh) | 一种geo导航卫星高灵敏度载波跟踪方法 | |
CN113009522B (zh) | 多普勒频率残差校正的长时间相干积分捕获算法模块 | |
CN108169772B (zh) | 一种加窗fft的卫星信号捕获方法 | |
CN104898136A (zh) | 一种北斗二代卫星b1信号的捕获方法及其系统 | |
CN106506037B (zh) | 一种非相干扩频通信系统的解调方法 | |
CN103941269A (zh) | 用于卫星导航系统的pn码捕获方法 | |
CN107037457A (zh) | 一种基于Inmarsat系统的星基增强接收机 | |
CN103760578A (zh) | 一种gnss卫星导航信号的无模糊跟踪方法 | |
CN110764120B (zh) | 一种高灵敏卫星导航信号捕获方法 | |
CN108627863B (zh) | 基于dft的gnss导航电文数据解调方法、装置及gnss接收终端 | |
CN106019335B (zh) | 一种卫星导航信号的斜率位同步法 | |
Li et al. | Highly sensitive weak signal acquisition method for GPS/compass | |
CN109547129B (zh) | 一种步进加权Tong检测方法 | |
CN107843905B (zh) | 一种快速高动态gnss频率搜索方法 | |
CN108833320B (zh) | 一种基于Keystone变换的宽带DS/FH信号参数估计方法及系统 | |
CN108226968B (zh) | 一种导航信号快速捕获方法 | |
WO2013140910A1 (ja) | 信号サーチ方法、信号サーチプログラム、信号サーチ装置、gnss信号受信装置、および情報機器端末 |
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 |