CN103630916B - 基于双fft频域滤波码捕获方法 - Google Patents
基于双fft频域滤波码捕获方法 Download PDFInfo
- Publication number
- CN103630916B CN103630916B CN201310627674.XA CN201310627674A CN103630916B CN 103630916 B CN103630916 B CN 103630916B CN 201310627674 A CN201310627674 A CN 201310627674A CN 103630916 B CN103630916 B CN 103630916B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- code
- signals
- frequency domain
- 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
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
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)
Abstract
本发明公开了一种基于双FFT频域滤波码捕获方法,包括以下步骤:1)对接受到的中频信号进行下变频;2)将基带信号分为Qms数据,分别与本地码频域值的复共轭相乘,再分别取傅里叶逆变换;3)分段处理;4)快速傅里叶变换;5)低通滤波;6)滤波后的频域信号进行快速傅里叶反变换;7)对得到的信号进行非相干积累,构造判决统计量;8)门限判决,大于门限即为检测信号。本发明的方法将经过PMF的信号先经过一个低通滤波器,滤除有用信号带外的频率,再对滤波后的信号进行重构,有效地滤除了噪声,这样信号再进行非相干积累,可以更大程度地提高信号的信噪比,进而提高检测概率。
Description
技术领域:
本发明属于导航信号检测技术领域,可用于GPS、北斗弱信号捕获。
背景技术:
全球卫星定位***(Global Navigation Satellite System,GNSS)信号是经过直接序列扩频调制的扩频信号,对经过直接扩频调制的卫星导航信号的捕获是GNSS***需要解决的首要问题。串行捕获是一种捕获技术,但由于此方法采用时域计算接受码和本地码相关,一次计算只能得到一个码相位的相关值且计算量大,所以出现了并行捕获算法,将时域计算接受码和本地码相关转化为频域方法计算且频域有相应的快速傅里叶变换(FastFourier Transformation,FFT)算法节省了大量的捕获时间,但由于需要捕获接收GNSS弱信号,需要大数据量的相关积累,所以出现了部分匹配滤波(partial matched filter,PMF)算法来计算接受码和本地码相关,其中具有代表性的方法是基于非相干积累的部分匹配滤波(partial matched filter based on non-coherent accumulation,PMF_NCA),就是将进行下变频后的基带信号进行部分匹配滤波,然后进行非相干积累,最后将积累的判决量进行阈值判定,采用非相干积累有两点好处:第一提高了信噪比进而提高了检测概率,第二克服了导航数据调制所带来的比特翻转问题。但是非相干积累采用将信号平方进行积累,所以将噪声幅度也给平方了,所以非相干积累提高检测概率的能力是有限的。
由于非相干积累将噪声幅度也增大了,所以本发明方法就是在非相干积累之前根据信号特点加以自适应滤波器,将滤波后的信号再进行非相干积累,提高了信噪比,进而提高检测概率,由于本方法提高了信号的信噪比,所以对GNSS弱信号捕获具有重要意义。
发明内容:
本发明的目的在于针对现有技术的不足,提出一种提高检测概率的基于双FFT频域滤波码捕获方法。
为了实现上述目的,本发明的方法包括如下步骤:
(1)在下变频单元中,对接受到的中频信号r(n)进行下变频处理:
其中sr(n)为基带信号,P(n)为中频信号r(n)的噪声且为复高斯噪声,实部、虚部都服从均值为0方差为δ2高斯分布,f1为信号频率,Δf为搜索步长,i为搜索次数,c(n)为GPS的C/A码或北斗信号的测距码。
(2)在频域相关单元,将sr(n)信号分为Q个数据段,每段为相关信号的码周期,表示为sr1(n),…,srq(n),…,srQ(n),分别将它们与本地码频域值的复共轭相乘,再分别取傅里叶逆变换,可以表示为:
其中FFT()表示取信号快速傅里叶变换,IFFT()表示求信号的快速傅里叶逆变换,conj()表示求信号的复共轭。
(3)在分段处理单元,上式Rq(n)中设0≤n≤p,当本地码与接收码同步时,n=p,即码相位取p时同步,则
s=[R1(p)...Rq(p)...RQ(p)]
其中Rq(p)为第q段同步时的相关值,即
其中Co为每段的采样点数,f0为基带信号的频率,W(q)为基带信号噪声且为复高斯噪声,实部、虚部都服从均值为0方差为Co*δ2高斯分布。
(4)在FFT单元中,对s向量进行快速傅里叶变换:
(5)在低通滤波器单元,对信号sr(k)进行低通滤波:
kmax为低通滤波器的截止频率,sf(k)为频率滤波后的频域信号。
(6)在IFFT单元,对sf(k)进行快速傅里叶反变换:
(7)在非相干积累单元,对信号sf(n)进行非相干积累,并构造判决统计量Z:
(8)门限判决,此时统计量若没有信号则服从卡方分布,则总的虚警概率为:
可以通过设定的虚警概率求得门限,跟统计量进行比较,大于门限即为检测到信号。
本发明的有益效果:与PMF_NCA相比,本发明的方法具有如下有点:
1)由于采用频域相关算法使本方法能一次相关计算所***相位对应的相关值;
2)由于本方法采用频域滤波算法,使滤波后的有用信号能量没有减少且使滤波后噪声信号的能量和方差,提高了检测概率。
附图说明:
图1是F_F检测框图;
图2是F_F检测详细检测流图;
图3为本发明方法与现有PMF_NCA方法在高斯白噪声信道下对GPS信号处理后的检测概率曲线;
图4为本发明方法与现有PMF_NCA方法在衰落信道下对GPS信号处理后的检测概率曲线;
图5为本发明方法与现有PMF_NCA方法在高斯白噪声信道下对北斗信号处理后的检测概率曲线;
图6为本发明方法与现有PMF_NCA方法在衰落信道下对北斗信号处理后的检测概率曲线。
具体实施方式:
为使本发明的优点和技术方案更加清楚明白,下面结合以下仿真数据和图像对本发明的技术方案作进一步说明。
1.仿真条件
以北斗、GPS信号为例进行了计算机仿真,在同一虚警概率下分别对F_F方法、PMF_NCA方法的蒙特卡洛仿真的检测概率结果进行对比,具体仿真参数:无偏码***,调制方式为BPSK,对于GPS信号:采用1023阶C/A码扩频,码片速率为Rc=1.023Mc/s,用于捕获总数据长度为8ms,***采样频率为1.023MHZ,部分匹配滤波器的数据长度为1ms,即Q=8,kmax=Q/2,虚警概率为Pfa=0.02;对于北斗信号:码片速率为Rc=2.046Mc/s,用于捕获总数据长度为8ms,***采样频率为2.046MHZ,部分匹配滤波器的数据长度为2ms,即Q=4,kmax=Q/2。信道采用高斯白噪声和一条可分辨路径的瑞利衰落信道,信道的多普勒频率为10HZ,
2.仿真内容
让不同频偏的信号通过高斯信道、衰落信道在经过本发明的方法和PMF_NCA方法进行处理得出的检测概率进行对比,从而通过检测概率的高低区分出两种方法的好坏。
仿真1:分别用本发明方法和现有的PMF_NCA方法对通过高斯信道的GPS信号进行处理所得检测概率曲线,如图3所示,其中:
图3a为频偏f为0HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图3b为频偏f为51.15HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图3c为频偏f为102.3HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
表1表示在一定信噪比下分别用本发明方法和PMF_NCA方法处理后的检测概率值的大小对比。
表1本发明方法与现有PMF_NCA方法在高斯白噪声信道下检测概率
仿真2:分别用本发明方法和现有的PMF_NCA方法对通过衰落信道的GPS信号进行处理所得检测概率曲线,如图4,其中:
图4a为频偏f为0HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图4b为频偏f为51.15HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图4c为频偏f为102.3HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
表2表示在一定信噪比下分别用本发明方法和PMF_NCA方法处理后的检测概率值的大小对比。
表2本发明方法与现有PMF_NCA方法在衰落信道下检测概率
仿真3:分别用本发明方法和现有的PMF_NCA方法对通过高斯信道的北斗信号进行处理所得检测概率曲线,如图5,其中:
图5a为频偏f为0HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图5b为频偏f为51.15HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图5c为频偏f为102.3HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
仿真4:分别用本发明方法和现有的PMF_NCA方法对通过衰落信道的北斗信号进行处理所得检测概率曲线,如图6,其中:
图6a为频偏f为0HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图6b为频偏f为51.15HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
图6c为频偏f为102.3HZ时分别用本发明方法和现有的PMF_NCA方法处理后检测概率曲线对比图。
3.仿真结果分析
通过实验结果可以发现,在同一信噪比下,无论是北斗信号、还是GPS信号,无论通过高斯信道还是衰落信道,用本发明的方法都会比PMF_NCA方法检测概率高,可以说明本发明方法的有效性和普遍适用性。
上述实施例仅用以说明本发明的技术方案,而非对本发明的构思和保护范围进行限定,本领域的普通技术人员对本发明的技术方案进行修改或者等同替换,而不脱离技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围中。
Claims (1)
1.基于双FFT频域滤波码捕获方法,其特征在于:包括以下步骤:
(1)对接受到的中频信号r(n)进行下变频处理:
其中sr(n)为基带信号,P(n)为中频信号r(n)的噪声且为复高斯噪声,实部、虚部都服从均值为0方差为δ2高斯分布,f1为信号频率,Δf为搜索步长,i为搜索次数;
(2)将sr(n)信号分为Qms数据,表示为sr1(n),…,srq(n),…,srQ(n),分别将它们的频域值与本地码频域值的复共轭相乘,再分别取傅里叶逆变换,表示为:
其中FFT()表示取信号快速傅里叶变换,IFFT()表示求信号的快速傅里叶逆变换,conj()表示求信号的复共轭,c(n)表示本地C/A码;
(3)上式Rq(n)中设0≤n≤p,当本地码与接收码同步时,n=p,即码相位取p时同步,则
s=[R1(p)...Rq(p)...RQ(p)]
其中Rq(p)为第qms同步时的相关值,即
其中Co为1ms的点数,f0为基带信号的频率,W(q)为基带信号噪声且为复高斯噪声,实部、虚部都服从均值为0方差为Co*δ2高斯分布;
(4)对s向量进行快速傅里叶变换,得到sr(k):
(5)对信号sr(k)进行低通滤波:
kmax为低通滤波器的截止频率,sf(k)为频率滤波后的频域信号;
(6)对频率滤波后的频域信号sf(k)进行快速傅里叶反变换:
(7)对信号sf(n)进行非相干积累,并构造判决统计量Z:
(8)门限判决,此时若没有信号则统计量服从卡方分布,则总的虚警概率为:
通过设定的虚警概率求得门限,跟统计量进行比较,大于门限即为检测到信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310627674.XA CN103630916B (zh) | 2013-11-29 | 2013-11-29 | 基于双fft频域滤波码捕获方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310627674.XA CN103630916B (zh) | 2013-11-29 | 2013-11-29 | 基于双fft频域滤波码捕获方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630916A CN103630916A (zh) | 2014-03-12 |
CN103630916B true CN103630916B (zh) | 2017-04-26 |
Family
ID=50212166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310627674.XA Active CN103630916B (zh) | 2013-11-29 | 2013-11-29 | 基于双fft频域滤波码捕获方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630916B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360360A (zh) * | 2014-11-07 | 2015-02-18 | 中国电子科技集团公司第二十研究所 | 一种gps卫星捕获判决的高效算法 |
CN105242286A (zh) * | 2015-05-29 | 2016-01-13 | 西安电子科技大学宁波信息技术研究院 | 一种基于小波域滤波码捕获方法 |
US20160377711A1 (en) * | 2015-06-26 | 2016-12-29 | Delphi Technologies, Inc. | Radar signal processing for automated vehicles |
CN106093981B (zh) * | 2016-06-06 | 2018-09-25 | 电子科技大学 | 基于优化并行码相位搜索的gps捕获电路 |
CN106656303B (zh) * | 2016-09-30 | 2019-08-02 | 西安电子科技大学 | 一种基于星载天线捕获跟踪指向***的信号检测方法 |
CN107172554A (zh) * | 2017-05-02 | 2017-09-15 | 贝尔顺(中国)科技股份有限公司 | 一种听力检测的控制***及控制方法 |
CN109752740A (zh) * | 2018-12-28 | 2019-05-14 | 航天信息股份有限公司 | 一种利用多通道分段相关和fft进行信号捕获的方法及*** |
CN112910498B (zh) * | 2021-01-25 | 2022-02-15 | 上海航天电子通讯设备研究所 | Pmf-fft测控信号捕获装置和方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1152255A2 (en) * | 2000-05-03 | 2001-11-07 | Lucent Technologies Inc. | Fast GPS signal acquisition in GPS systems |
CN101005293A (zh) * | 2006-12-31 | 2007-07-25 | 浙江大学 | 高动态卫星导航接收机中的信号快速捕获方法及装置 |
CN203086464U (zh) * | 2013-01-21 | 2013-07-24 | 西安空间无线电技术研究所 | 一种基于频域的pn码快速捕获装置 |
-
2013
- 2013-11-29 CN CN201310627674.XA patent/CN103630916B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1152255A2 (en) * | 2000-05-03 | 2001-11-07 | Lucent Technologies Inc. | Fast GPS signal acquisition in GPS systems |
CN101005293A (zh) * | 2006-12-31 | 2007-07-25 | 浙江大学 | 高动态卫星导航接收机中的信号快速捕获方法及装置 |
CN203086464U (zh) * | 2013-01-21 | 2013-07-24 | 西安空间无线电技术研究所 | 一种基于频域的pn码快速捕获装置 |
Non-Patent Citations (1)
Title |
---|
基于码域频域FFT的伪码捕获的改进与实现;王鹏宇等;《***工程与电子技术》;20080831;第30卷(第08期);1407-1411 * |
Also Published As
Publication number | Publication date |
---|---|
CN103630916A (zh) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103630916B (zh) | 基于双fft频域滤波码捕获方法 | |
CN104280750B (zh) | 基于部分匹配滤波fft算法的长码捕获***及方法 | |
CN103746722B (zh) | 一种跳频信号跳周期和起跳时间估计方法 | |
CN103645483B (zh) | 一种弱信号环境下北斗信号捕获方法 | |
CN103412287B (zh) | 基于lvd的线性调频信号参数估计方法 | |
CN102914782B (zh) | 一种适用于gps弱信号的快速捕获方法 | |
CN105577229B (zh) | Cpu辅助gpu实现扩频信号快速捕获方法 | |
CN102323601B (zh) | 一种gnss-boc调制信号的捕获方法 | |
CN104181556B (zh) | 一种基于重叠差分循环相干积分的boc调制信号捕获方法 | |
CN103675758B (zh) | 一种双曲调频信号周期斜率和起始频率估计方法 | |
CN103634026A (zh) | 一种基于fpga的数字零中频自适应陷波方法 | |
CN107085222A (zh) | 一种 boc 信号捕获方法和卫星导航接收机 | |
CN105785324A (zh) | 基于mgcstft的线性调频信号参数估计方法 | |
CN105553506B (zh) | 一种长码扩频信号的快速捕获方法及装置 | |
CN104901909B (zh) | 一种α非高斯噪声下chirp信号的参数估计方法 | |
CN105242286A (zh) | 一种基于小波域滤波码捕获方法 | |
CN101527698A (zh) | 基于希尔伯特黄变换和自适应陷波的非平稳干扰抑制方法 | |
CN102158248A (zh) | 基于循环谱理论的ds/fh扩频信号参数估计方法 | |
CN102546499A (zh) | 一种实线性调频信号的分数阶信道化接收方法 | |
CN106879068A (zh) | 一种强多径环境下信号的到达时间估计方法 | |
CN109379310A (zh) | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 | |
CN103427870A (zh) | 基于稀疏傅里叶变换的深度扩频捕获与干扰抑制方法 | |
CN103941269A (zh) | 用于卫星导航***的pn码捕获方法 | |
CN105445767B (zh) | 基于平均模糊函数的boc信号参数盲估计方法 | |
CN102655419B (zh) | 一种捕获扩频信号的判决函数计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |