CN105301607B - 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法 - Google Patents

一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法 Download PDF

Info

Publication number
CN105301607B
CN105301607B CN201510810632.9A CN201510810632A CN105301607B CN 105301607 B CN105301607 B CN 105301607B CN 201510810632 A CN201510810632 A CN 201510810632A CN 105301607 B CN105301607 B CN 105301607B
Authority
CN
China
Prior art keywords
module
gnss signal
signal
way
gnss
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
Application number
CN201510810632.9A
Other languages
English (en)
Other versions
CN105301607A (zh
Inventor
肖永田
马百良
王志国
李廷
王庆海
莫钧
屠恩源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
WUHAN MENGXIN TECHNOLOGY Co Ltd
Original Assignee
WUHAN MENGXIN TECHNOLOGY Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by WUHAN MENGXIN TECHNOLOGY Co Ltd filed Critical WUHAN MENGXIN TECHNOLOGY Co Ltd
Priority to CN201510810632.9A priority Critical patent/CN105301607B/zh
Publication of CN105301607A publication Critical patent/CN105301607A/zh
Application granted granted Critical
Publication of CN105301607B publication Critical patent/CN105301607B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • 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
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69—Spread spectrum techniques
    • H04B1/707—Spread spectrum techniques using direct sequence modulation
    • H04B1/7097—Interference-related aspects
    • H04B1/71—Interference-related aspects the interference being narrowband interference
    • H04B1/7102—Interference-related aspects the interference being narrowband interference with transform to frequency domain

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法，其中，多频GNSS信号窄带干扰抑制系统包括至少一个预处理模块、至少一个缓存处理模块、至少一个缓存控制模块、至少一个延迟模块、复用模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块、解复用模块、至少一个叠加处理模块以及分时复用控制模块。本发明能够对单频率以及多频率的GNSS信号进行窄带干扰抑制，不仅具有良好的抑制效果，还能够降低硬件资源成本以及系统的功耗，实用性较强。

Description

一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法
技术领域
本发明涉及卫星导航数字信号处理技术领域，具体涉及一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法。
背景技术
全球导航卫星系统(Global Navigation Satellite System，GNSS)信号在全球范围内提供全天候、实时、连续的高精度位置信息、速度和时间信息，在人们生活中发挥着日益重要的作用。
GNSS信号采用扩频通信技术，本身具有强抗干扰能力，其抗干扰能力由扩频处理增益及带宽决定。但在实际通信中，由于信号发射功率、设备实现成本及受增大带宽代价的限制，扩频通信的抗干扰能力是有限的。GNSS信号到达接收天线的信号十分微弱，通常淹没在噪声之中，容易受到外界干扰，特别是人为的干扰。当干扰信号功率超过扩频通信的抗干扰容限时，通信性能就得不到保证，必须使用抗干扰技术来保证通信质量。按照干扰带宽相对于GNSS信号带宽的大小，将其分为宽带干扰和窄带干扰两种。其中，对GNSS接收机影响最大、最常见的干扰类型是窄带干扰。扩频通信中的窄带干扰抑制(Narrow BandInterference，NBI)技术通常在解扩处理之前进行，削弱窄带干扰部分，提高接收信号的载噪比。
现有技术中的窄带干扰抑制装置主要分为时域的基于自适应滤波的窄带干扰抑制装置和频域的基于FFT(Fast Fourier Transform，快速傅立叶变换)/IFFT(InverseFast Fourier Transform，快速傅立叶逆变换)的窄带干扰抑制装置。时域滤波具有硬件实现复杂度高、收敛时间缓慢、群时延变化等缺点；频域的基于FFT/IFFT的窄带干扰抑制装置是采用FFT计算接收到的信号的频谱特性，在频谱中查找干扰出现的频点，并对该频点的信号进行相应的滤波处理，从而降低干扰在接收信号中的影响。频域算法由于可以利用FFT/IFFT快速算法减少计算量，且能保证干扰抑制系统的线性相位，不存在收敛性问题，得到了广泛应用。
目前常见的针对单频率GNSS信号的频域滤波主要有以下两种方式：
第一种方式，简单的基于FFT/IFFT的干扰抑制技术，将GNSS信号通过加窗处理、FFT计算、NBI滤波处理以及IFFT计算后输出，该方法虽然可以去除窄带干扰，但加窗处理在降低FFT时的频谱泄漏，避免干扰扩散至其他频点的同时，存在引起输入信号畸变，导致输出信号的信噪比受到损失的缺点。
第二种方式，基于OFFT(Overlapped Fast Fourier Transform，叠加的快速傅立叶变换)的干扰抑制技术，在第一种方式的基础上，采用将GNSS信号分为两路，分别经1/2FFT块长度延迟和不延迟，分别独立进行加窗、FFT、干扰频率剔除、IFFT后，然后按照一定的算法进行合路，该方法具有设计容易、可同时滤除多个窄带干扰等优点，但需要4个FFT(2个FFT模块和2个IFFT模块)引擎，硬件实现复杂度高，且频率分辨率较低，对信号衰减较大。
随着GNSS技术的快速发展和应用需求的不断提高，利用多系统信号进行联合定位的多模多频GNSS系统表现出了极大的优势。多模多频接收机可同时接收多个频率的卫星信号，利用多系统联合定位技术，使定位精度和可靠性得到更大的提高，迅速成为导航领域的研究和应用热点。但目前在多模多频接收机中使用OFFT技术进行NBI滤波时，若同时接收N个频率的卫星信号，需要复制N个单频率信号的OFFT处理模块，此时硬件资源成本将成倍增加，同时还增加了系统功耗。
另外，在FFT/IFFT硬件实现时，通常采用典型的DIF(Decimation in Frequency，频域抽取)进行FFT处理。即FFT计算数据按顺序输入，计算结果逆序输出，采用这种结构进行IFFT时，输入数据之前要对数据进行缓存以调整顺序的处理，这不仅占用大量存储空间，且加大了流水处理延时，既提高了系统对存储资源的需求，又降低了系统实时性。
为解决上述问题，提出一种逆序输入顺序输出的FFT结构设计方法，解决了IFFT带来的额外存储需求和流水迟滞问题，但是，GNSS输入中频信号按时域顺序输入，若按此方法，必须在加窗处理后FFT计算前对输入数据顺序按“倒位序”处理，以给FFT输入的数据按FFT蝶形运算顺序，同时保证了延迟M/2处理模块和加窗处理模块按时域顺序处理，但是，“倒位序”处理需要增加存储器进行存储，耗费较大的硬件资源。
发明内容
本发明所要解决的技术问题是提供一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法，能够克服上述现有技术的不足，实现对单频率GNSS信号以及多频率GNSS信号的窄带干扰抑制。
本发明解决上述技术问题的技术方案如下：
依据本发明的一个方面，提供了一种单频GNSS信号窄带干扰抑制装置，所述装置包括延迟模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及叠加处理模块；
所述延迟模块，用于接收一路GNSS信号,并将其延迟M/2个采样点输出给所述加窗处理模块，其中M为FFT块长度，具体为2的指数幂；
所述加窗处理模块，用于接收一路GNSS信号，以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，并对两路时域GNSS信号分别进行加窗处理；
所述FFT计算模块，包括实部输入与虚部输入，用于对经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号；
所述NBI滤波处理模块，用于对经过FFT计算的两路频域GNSS信号进行滤波处理；
所述IFFT计算模块，包括实部输入与虚部输入，用于对经过滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号，并将其输出；
所述叠加处理模块，用于对接收的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，并将其输出。
依据本发明的另一个方面，提供了一种单频GNSS信号窄带干扰抑制方法，所述方法包括：
S301、将单频率的时域GNSS信号分为两路，其中一路信号延迟M/2个采样点，另一路信号不延迟，并将延迟后的以及未延迟的两路GNSS信号分别进行加窗处理；
S302、将经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号；
S303、将经过FFT计算的两路频域GNSS信号进行滤波处理，并将经滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号；
S304、将经过IFFT计算的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，以进行后续的解扩处理。
依据本发明的再一个方面，提供了一种多频GNSS信号窄带干扰抑制系统，所述系统包括单频GNSS信号窄带干扰抑制装置以及至少一个预处理模块、至少一个缓存处理模块、至少一个缓存控制模块、复用模块、解复用模块以及分时复用控制模块；
所述预处理模块，包括格式转换单元以及信号检测单元，用于在RF时钟域下对接收到的GNSS信号进行格式转换，并对格式转换后的GNSS信号进行干扰信号检测，以确定所述GNSS信号是否需要进行窄带干扰抑制处理，且根据确定结果生成通道干扰抑制使能控制信号并发送给所述缓存控制模块以及将检测之后的GNSS信号输出给所述缓存处理模块；还用于将所述缓存处理模块返回的GNSS信号输出，以进行后续解扩处理；
所述缓存处理模块，用于在RF时钟域下写入干扰信号检测之后的GNSS信号，然后在所述缓存控制模块的控制下读出所述GNSS信号，并根据GNSS信号是否需要进行窄带干扰抑制处理的确定结果，将读出的GNSS信号返回给对应连接的所述预处理模块输出或者发送给对应连接的所述延迟模块以及所述复用模块；
所述缓存控制模块，用于在系统时钟域下向所述分时复用控制模块发送所述通道干扰抑制使能控制信号以及缓存数据就绪指示信号，并接收所述分时复用控制模块反馈的通道缓存读使能信号，进而控制所述缓存处理模块工作；
所述复用模块，用于根据所述分时复用控制模块发送的多通道复用控制信号，将依次接收的多个所述缓存处理模块发送的未延迟的多个频率的GNSS信号，以及多个所述延迟模块发送的延迟的多个频率的GNSS信号整合到一个通道中，并按照频率的不同，分时段每依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块；
所述解复用模块，用于根据所述分时复用控制模块发送的多通道解复用控制信号，对分时段接收的处于一个通道中的多个包含两路信号的不同频率的GNSS信号进行解复用，使多个两路GNSS信号按照频率的不同处在多个通道中，并分别发送给多个所述叠加处理模块。
依据本发明的另一个方面，提供了一种多频GNSS信号窄带干扰抑制方法，其特征在于，所述方法包括：
S401、N个所述预处理模块在RF时钟域下，对接收的N个GNSS信号分别进行预处理，并将预处理之后的GNSS信号输出给对应连接的N个所述缓存处理模块，同时还通过N个所述缓存控制模块分别将通道干扰抑制使能控制信号发送给所述分时复用控制模块；
S402、在RF时钟域下N个所述缓存处理模块在写入地址依次相差M/4的情况下分别写入N个GNSS信号，在依次写满之后，其对应的N个所述缓存控制模块依次发送缓存数据就绪指示信号给所述分时复用控制模块；所述分时复用控制模块根据所述缓存数据就绪指示信号以及之前接收的所述通道干扰抑制使能控制信号，生成通道缓存读使能信号并依次反馈给N个所述缓存控制模块；
S403、N个所述缓存控制模块根据所述缓存读使能信号分别控制N个所述缓存处理模块读出N个GNSS信号，若不需要对GNSS信号进行窄带干扰抑制，将GNSS信号返回给对应连接的所述预处理模块输出以进行后续解扩处理，结束当前GNSS信号的抑制流程；若需要对GNSS信号进行窄带干扰抑制，将GNSS信号发送给对应连接的所述延迟模块以及所述复用模块；
S404、N个所述延迟模块在系统时钟域下，依次分别将接收的GNSS信号延迟M/2个采样点，并将延迟后的GNSS信号发送给所述复用模块，所述复用模块根据所述分时复用控制模块发送的多通道复用控制信号，依次将接收的延迟的N个频率的GNSS信号以及未延迟的N个频率的GNSS信号整合到一个通道中，并按照频率的不同分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块，直到N个不同频率的包含两路信号的GNSS信号全部被发送至所述加窗处理模块；
S405、所述加窗处理模块以及FFT计算模块依次对接收的延迟以及未延迟的两路GNSS信号进行加窗处理以及FFT计算处理，并将得到的计算结果发送给所述NBI滤波处理模块，所述NBI滤波处理模块根据所述分时复用控制模块发送的通道NBI滤波配置信号对滤波参数进行配置，进而对两路GNSS信号进行滤波处理，并将滤波结果发送给所述IFFT计算模块，所述IFFT计算模块对两路GNSS信号进行IFFT计算，并将计算之后的两路GNSS信号发送给所述解复用模块，同时向所述分时复用控制模块发送IFFT完成指示信号，所述分时复用控制模块根据所述完成指示信号向所述解复用控制模块发送多通道解复用控制信号；
S406、所述解复用模块根据所述多通道解复用控制信号，使分时段接收的N个两路GNSS信号按照频率的不同处在N个通道中，并分别发送给对应连接的N个所述叠加处理模块,N个所述叠加处理模块分别将接收的两路时域GNSS信号进行叠加处理，并将得到的一路时域GNSS信号分别返回给对应连接的N个所述缓存处理模块以及预处理模块输出，以进行后续解扩处理。
本发明提供的一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法，实现了对单频率GNSS信号以及多频率的GNSS信号的窄带干扰抑制。其中，所述FFT计算模块与IFFT计算模块均设有实部输入与虚部输入，能够同时对两路GNSS信号进行计算，实现了FFT与IFFT计算模块的共享复用，使FFT/IFFT计算模块的硬件成本降低一半；所述加窗处理模块能够对两路GNSS信号进行加窗处理，共享一个窗函数系数ROM表，并根据窗函数的对称性只存储M/2个加窗函数系数，如此将存储窗函数系数的ROM资源节省为原来的1/4；所述NBI滤波处理模块同时对两路GNSS信号进行滤波处理，并能根据不同频率的GNSS信号配置不同的滤波参数，保持滤波的独立性，有效抑制窄带干扰；在对N个频率的GNSS信号进行窄带干扰抑制时，所述加窗处理模块、FFT计算模块、NBI滤波处理模块以及IFFT计算模块的采用分时复用的机制，将硬件资源降低为原来的1/N；所述缓存处理模块将GNSS信号按照FFT蝶形运算顺序读出，并且所述延迟模块以及加窗处理模块均按照FFT蝶形运算顺序处理GNSS信号，将FFT计算模块前的“倒位序”去掉，节省M个数据存储单元；所述分时复用控制单元能够根据输入的GNSS信号的频率的数量自动调整系统时钟频率，降低系统功耗；另外所述预处理模块能够对GNSS信号进行干扰信号检测，实现对GNSS信号的自适应窄带干扰抑制处理，并且在所有的GNSS信号都不需要进行窄带干扰抑制的情况下，关闭所述延迟模块、叠加处理模块、复用模块、加窗模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块，降低系统功耗。本发明能够对单频率以及多频率的GNSS信号进行窄带干扰抑制，不仅具有良好的抑制效果，还能够降低硬件资源成本以及系统的功耗，并具有较强的实用性。
附图说明
图1为本发明实施例一的一种单频GNSS信号窄带干扰抑制装置示意图；
图2为本发明实施例二的一种单频GNSS信号窄带干扰抑制方法流程图；
图3为本发明实施例三的一种多频GNSS信号窄带干扰抑制系统示意图；
图4为本发明实施例四的一种多频GNSS信号窄带干扰抑制方法流程图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
实施例一、一种单频GNSS信号窄带干扰抑制装置。下面结合图1对本实施例提供的装置进行详细说明。
参见图1，本实施例提供的装置包括延迟模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及叠加处理模块。
所述延迟模块，用于接收一路GNSS信号，并将其延迟M/2个采样点输出给所述加窗处理模块，其中M为FFT块长度，具体为2的指数幂。
所述加窗处理模块，用于接收一路GNSS信号，以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，并对两路时域GNSS信号分别进行加窗处理。
具体的，所述加窗处理模块用于接收一路GNSS信号，以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，其中M为FFT块长度，具体为2的指数幂，并将两路时域GNSS信号分别乘以在加窗函数的系数ROM(Read Only Memory，只读存储器)表中查找得到的同一加窗系数，进行加窗处理，实现加窗处理模块关于两路GNSS信号的共享复用。其中，所述加窗函数具体包括矩形窗、三角窗、汉宁窗、海明窗以及布莱克曼窗，所述加窗函数均为对称函数，若FFT的块长度为M，对应于加窗函数的M个系数分别为F1,F2,…,FM，则有Fi＝FM+1-i(i＝1,2,…,M/2-1)，加窗函数系数ROM表中只需存储M/2个加窗系数，加窗函数系数ROM表大小只需M/2，可将存储窗函数系数的ROM表资源节省为原来所需的1/4。
所述FFT计算模块，包括实部输入与虚部输入，用于对经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号。
具体的，所述FFT计算模块将经过加窗处理的两路GNSS信号中的一路GNSS信号实部输入进行FFT计算，另一路GNSS信号虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号，将两路GNSS信号从时域数据转换到频域数据，实现FFT计算模块关于两路GNSS信号的共享复用，使FFT计算模块资源降为原来的一半。
所述NBI滤波处理模块，用于对经过FFT计算的两路频域GNSS信号进行滤波处理。
具体的，所述NBI滤波处理模块可采用各种滤波抑制算法，例如门限检测法、K谱线算法以及中值滤波法等，其中，输入的GNSS信号为单频率时，所述NBI滤波处理模块的滤波参数为提前设置好的，保持滤波的独立性。
所述IFFT计算模块，包括实部输入与虚部输入，用于对经过滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号，并将其输出。
具体的，所述IFFT计算模块将经过滤波处理的两路频域GNSS信号中的一路GNSS信号实部输入进行IFFT计算，另一路GNSS信号虚部输入进行IFFT计算，并将计算得到的复数IFFT信号分离计算得到两路实数IFFT信号，将两路GNSS信号从频域数据转换到时域数据，实现IFFT计算模块关于两路GNSS信号的共享复用，使IFFT计算模块资源降为原来的一半。
所述叠加处理模块，用于对接收的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，并将其输出。
实施例二、一种单频GNSS信号窄带干扰抑制方法。下面结合图2对本实施例提供的方法进行详细说明。
参见图1及图2，S301、将单频率的时域GNSS信号分为两路，其中一路信号延迟M/2个采样点，另一路信号不延迟，并将延迟后的以及未延迟的两路GNSS信号分别进行加窗处理。
具体的，将单频率的时域GNSS信号分为两路，延迟模块接收其中一路信号并将其延迟M/2个采样点，并将延迟后的GNSS信号发送给所述加窗处理模块，所述加窗处理模块接收另一路GNSS信号以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，并将其分别乘以在加窗函数系数ROM表中查找得到的同一加窗系数，进行加窗处理，其中，所述加窗函数具体包括矩形窗、三角窗、汉宁窗、海明窗以及布莱克曼窗，所述加窗函数均为对称函数，若FFT的块长度为M，对应于加窗函数的M个系数分别为F1,F2,…,FM，则有Fi＝FM+1-i(i＝1,2,…,M/2-1)，加窗函数的系数ROM表中只需存储M/2个加窗系数，加窗函数系数ROM表大小只需M/2；然后所述加窗处理模块将处理过后的两路GNSS信号发送给所述FFT计算模块。
S302、将经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号，将两路GNSS信号从时域数据转换到频域数据。
S303、将经过FFT计算的两路频域GNSS信号进行滤波处理，并将经滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号。
具体的，所述NBI滤波处理模块对经过FFT计算的两路频域GNSS信号进行滤波处理，所述NBI滤波处理模块可采用各种滤波抑制算法，例如门限检测法、K谱线算法以及中值滤波法等，其中，输入的GNSS信号为单频率时，所述NBI滤波处理模块的滤波参数为提前设置好的，所述NBI滤波处理模块将滤波处理的结果发送给所述IFFT计算模块。
所述IFFT计算模块将经过滤波处理的两路频域GNSS信号中的一路GNSS信号实部输入进行IFFT计算，另一路GNSS信号虚部输入进行IFFT计算，并将计算得到的复数IFFT信号分离计算得到两路实数IFFT信号，将两路GNSS信号从频域数据转换到时域数据，然后将所述IFFT计算模块将计算结果发送给所述叠加处理模块。
S304、将经过IFFT计算的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，以进行后续的解扩处理。
实施例三、一种多频GNSS信号窄带干扰抑制系统。下面结合图1以及图3对本实施例提供的系统进行详细说明。
参见图1及图3，本实施例提供的系统包括单频GNSS信号窄带干扰抑制装置以及至少一个预处理模块、至少一个缓存处理模块、至少一个缓存控制模块、复用模块、解复用模块以及分时复用控制模块。
所述预处理模块，包括格式转换单元以及信号检测单元，用于接收未经过处理的GNSS信号，所述格式转换单元在RF时钟域下将不同位宽以及不同编码格式的GNSS信号转换为相同的位宽以及编码格式，所述信号检测单元对格式转换后的GNSS信号进行干扰信号检测，以确定所述GNSS信号是否需要进行窄带干扰抑制处理，且根据确定结果生成通道干扰抑制使能控制信号并发送给所述缓存控制模块，以及将检测之后的GNSS信号输出给所述缓存处理模块；还用于将所述缓存处理模块返回的GNSS信号输出，以进行后续解扩处理。
所述缓存处理模块，用于在RF时钟域下写入干扰信号检测之后的GNSS信号，然后在所述缓存控制模块的控制下读出所述GNSS信号，并根据GNSS信号是否需要进行窄带干扰抑制处理的确定结果，将读出的GNSS信号返回给对应连接的所述预处理模块输出或者发送给对应连接的所述延迟模块以及所述复用模块。
所述缓存控制模块，用于在系统时钟域下向所述分时复用控制模块发送所述通道干扰抑制使能控制信号以及缓存数据就绪指示信号，并接收所述分时复用控制模块反馈的通道缓存读使能信号，进而控制所述缓存处理模块工作。
具体的，所述缓存处理模块采用大小为2个FFT块长度的双端口SRAM(StaticRandom Access Memory，静态随机存储器)以进行乒乓操作；所述缓存处理模块在RF时钟域下写入预处理之后的GNSS信号，所述缓存控制模块根据接收的所述缓存读使能信号控制所述缓存处理模块读出GNSS信号：若所述缓存读使能信号表明不需要对GNSS信号进行窄带干扰抑制，所述缓存处理模块在RF时钟域下顺序读出当前GNSS信号，并返回给对应连接的所述预处理模块输出；若所述缓存读使能信号表明需要对GNSS信号进行窄带干扰抑制，在系统时钟域下，按照FFT蝶形运算顺序读出当前GNSS信号，并发送给对应连接的所述延迟模块以及所述复用模块；所述缓存处理模块在读出当前GNSS信号的同时，在RF时钟域下继续写入下一个经预处理之后的GNSS信号，实现乒乓操作；另外，所述缓存处理模块在所述当前GNSS信号写入之前，提前一个RF时钟将当前GNSS信号之前的经过所述叠加处理模块处理并返回的GNSS信号读出，并返回给对应连接的所述预处理模块输出。
所述延迟模块与所述单频GNSS信号窄带干扰抑制装置中的延迟模块完全相同，用于将接收的GNSS信号延迟M/2个采样点,并将延迟后的GNSS信号输出。
所述复用模块，用于根据所述分时复用控制模块发送的多通道复用控制信号，将依次接收的多个所述缓存处理模块发送的未延迟的多个频率的GNSS信号，以及多个所述延迟模块发送的延迟的多个频率的GNSS信号整合到一个通道中，并按照频率的不同，分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块。
所述加窗处理模块与所述单频GNSS信号窄带干扰抑制装置中的加窗处理模块完全相同，用于接收延迟M/2个采样点的GNSS信号和未延迟的GNSS信号，并对两路时域GNSS信号分别进行加窗处理；所述加窗处理模块分时复用对多个频率的两路GNSS信号进行加窗处理。
其中，由于所述缓存处理模块按照FFT蝶形运算顺序读出当前GNSS信号，因此所述延迟模块以及加窗处理模块均按照FFT蝶形运算顺序对所述GNSS信号进行处理，若所述缓存处理模块在RF时钟域下写入M个数据采样点的顺序分别为D0、D1、D2、…、DM-1，那么在系统时钟域下读出的顺序为D0、DM/2、D1、DM/2+1、D2、DM/2+2、…、DM/2-1、DM-1，进而所述加窗处理模块从加窗函数ROM表中读出对应的M个加窗系数的顺序为F0、FM/2、F1、FM/2+1、F2、FM/2+2、…、FM/2-1、FM-1,加窗函数系数ROM表中只需存储M/2个加窗系数，加窗函数系数ROM表大小只需M/2。
所述缓存处理模块按FFT蝶形运算顺序读出GNSS信号，并且所述延迟模块以及加窗处理模块都按照FFT蝶形运算顺序进行处理，将FFT计算模块前的“倒位序”处理去掉，共节省M个数据存储单元(每一路GNSS信号节省M/2个数据存储单元)。
所述FFT计算模块与所述单频GNSS信号窄带干扰抑制装置中的FFT计算模块完全相同，用于对经过加窗处理的两路GNSS信号一路实部输入进行FFT计算，一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离得到两路实数FFT信号，将两路GNSS信号从时域数据转换到频域数据；所述FFT计算模块分时复用对多个频率的两路信号进行FFT计算。
所述NBI滤波处理模块与所述单频GNSS信号窄带干扰抑制装置中的NBI滤波处理模块完全相同，用于对经过FFT计算的两路频域GNSS信号进行滤波处理；所述NBI滤波处理模块分时复用对多个频率的两路信号进行滤波处理，并接收所述分时复用控制模块发送的针对不同频率GNSS信号的通道NBI滤波配置信号对滤波参数进行配置，进而保证不同频率GNSS信号滤波的独立性。
所述IFFT计算模块与所述单频GNSS信号窄带干扰抑制装置中的IFFT计算模块完全相同，用于对经过滤波处理的两路频域GNSS信号一路实部输入进行IFFT计算，一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离得到两路实数IFFT信号，将两路GNSS信号从频域数据转换到时域数据，然后将其输出；所述IFFT计算模块分时复用对多个频率的两路信号进行IFFT计算。
所述叠加处理模块与所述单频GNSS信号窄带干扰抑制装置中的叠加处理模块完全相同，用于对接收的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，并将其输出给所述缓存处理模块。
所述解复用模块，用于根据所述分时复用控制模块发送的多通道解复用控制信号，对分时段接收到的处于一个通道中的多个包含两路信号的不同频率的GNSS信号进行解复用，使多个两路GNSS信号按照频率的不同处在多个通道中，并分别发送给多个所述叠加处理模块。
所述分时复用控制模块，用于接收所述缓存控制模块发送的所述通道干扰抑制使能控制信号以及所述缓存数据就绪指示信号，然后向所述缓存控制模块反馈所述通道缓存读使能信号；用于向所述复用模块发送所述多通道复用控制信号控制所述复用模块工作以及向所述NBI滤波处理模块发送通道NBI滤波配置信号；用于接收IFFT计算模块发送的IFFT完成指示信号，并根据所述完成指示信号向所述解复用模块发送多通道解复用控制信号控制所述解复用模块工作。
另外，所述分时复用控制模块，还用于在所有的GNSS信号都不需要进行窄带干扰抑制时，控制所有的延迟模块、所有的叠加处理模块、复用模块、加窗模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块全部关闭,以降低系统的功耗。
另外，在实际接入的GNSS信号的频率个数为K时，所述分时复用控制模块调整系统时钟频率为RF时钟频率的K+1倍，这样若RF时钟域输入M个采样点数据的时间为T，在系统时钟域内读出M个采样点数据时间只需T/(K+1)。
N个不同频率的GNSS信号依次在所述延迟模块、加窗处理模块、FFT计算模块、NBI滤波处理模块以及IFFT计算模块中处理所需总时间是固定的，例如，设GNSS信号在一个RF时钟域下所需时间为T，所述延迟模块在系统时钟域下将N个不同频率的GNSS信号延迟M/2个采样点所需时间为T/2(N+1)，N个不同频率的GNSS信号以及延迟M/2个采样点的GNSS信号在系统时钟域下从加窗处理模块到NBI滤波处理模块依次处理所需总时间为N*T/(N+1)，所以总的处理时间约为(2N+1)*T/2(N+1)。
所述多频GNSS信号窄带干扰抑制系统具体适用于2到5个频率的GNSS信号的窄带干扰抑制，图3中示出的多频GNSS信号窄带干扰抑制系统为N取值为3时的系统。
实施例四、一种多频GNSS信号窄带干扰抑制方法。下面结合图1至图4对本实施例提供的方法进行详细说明。
参见图1至图4，S401、N个所述预处理模块在RF时钟域下，对接收的N个GNSS信号分别进行预处理，并将预处理之后的GNSS信号输出给对应连接的N个所述缓存处理模块，同时还通过N个所述缓存控制模块分别将通道干扰抑制使能控制信号发送给所述分时复用控制模块。
具体的，在RF时钟域下，N个所述预处理模块中的N个格式转换单元分别对接收的N个未经过处理的GNSS信号进行格式转换，将不同位宽以及不同编码格式的GNSS信号转换为相同的位宽以及编码格式的GNSS信号，N个信号检测单元对转换后的GNSS信号进行干扰信号检测，确定是否需要对所述GNSS信号进行窄带干扰抑制，然后根据确定结果生成通道干扰抑制使能控制信号，并分别通过N个所述缓存控制模块将通道干扰抑制使能控制信号发送给所述分时复用控制模块，同时还将检测之后的GNSS信号输出给对应连接的所述缓存处理模块，其中，N的取值范围为2到5。
S402、在RF时钟域下N个所述缓存处理模块在写入地址依次相差M/4的情况下分别写入N个GNSS信号，在依次写满之后，其对应的N个所述缓存控制模块依次发送缓存数据就绪指示信号给所述分时复用控制模块；所述分时复用控制模块根据所述缓存数据就绪指示信号以及之前接收的所述通道干扰抑制使能控制信号，生成通道缓存读使能信号并依次反馈给N个所述缓存控制模块。
具体的，N个所述缓存处理模块在RF时钟域下写入预处理之后的GNSS信号，写入地址依次相差M/4，第一个所述缓存处理模块从写入地址0开始写入，当写到M/4个采样点时，第二个所述缓存处理模块才从写入地址0开始写入；当第一个所述缓存处理模块写到M/2个采样点，第二个所述缓存处理模块写到M/4个采样点时，第三个所述缓存处理模块才从写入地址0开始写入，以此类推至第N个所述缓存处理模块；后续不停写入，各个所述缓存处理模块的写入地址不断加1，在写入地址写满后均会清零重新开始计数；这样保证N个所述缓存处理模块在写满之后生成的缓存数据指示信号刚好错开M/4个RF时钟域时间，依次写满之后，其对应的N个所述缓存控制模块依次发送所述缓存数据就绪指示信号给所述分时复用控制模块。
所述分时复用控制模块根据所述缓存数据就绪指示信号以及步骤S401中接收的对应的所述通道干扰抑制使能控制信号，生成是否对GNSS信号进行窄带干扰抑制的通道缓存读使能信号，并依次发送给N个所述缓存控制模块。
S403、N个所述缓存控制模块根据所述缓存读使能信号分别控制N个所述缓存处理模块读出N个GNSS信号，若不需要对GNSS信号进行窄带干扰抑制，将GNSS信号返回给对应连接的所述预处理模块输出以进行后续解扩处理，结束当前GNSS信号的抑制流程；若需要对GNSS信号进行窄带干扰抑制，将GNSS信号发送给对应连接的所述延迟模块以及所述复用模块。
具体的，N个所述缓存控制模块依次根据所述缓存读使能信号控制N个所述缓存处理模块读出N个GNSS信号：若所述缓存读使能信号表明不需要对GNSS信号进行窄带干扰抑制，所述缓存处理模块在RF时钟域下顺序读出当前GNSS信号，并发送给对应连接的所述预处理模块，所述预处理模块将其输出进行后续解扩处理，结束当前GNSS信号的抑制流程；若所述缓存读使能信号表明需要对GNSS信号进行窄带干扰抑制，所述缓存处理模块在系统时钟域下，按照FFT蝶形运算顺序读出当前GNSS信号，并发送给对应连接的所述延迟模块以及所述复用模块；所述缓存处理模块在读出当前GNSS信号的同时，在RF时钟域下继续写入下一个经预处理之后的GNSS信号，实现乒乓操作。
S404、N个所述延迟模块在系统时钟域下，依次分别将接收的GNSS信号延迟M/2个采样点，并将延迟后的GNSS信号发送给所述复用模块，所述复用模块根据所述分时复用控制模块发送的多通道复用控制信号，依次将接收的延迟的N个频率的GNSS信号以及未延迟的N个频率的GNSS信号整合到一个通道中，并按照频率的不同分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块，直到N个不同频率的包含两路信号的GNSS信号全部被发送至所述加窗处理模块。
具体的，N个所述延迟模块在系统时钟域下依次分别将接收的GNSS信号延迟M/2个采样点,并将延迟后的GNSS信号输出给所述复用模块，所述复用模块根据所述分时复用控制模块发送的多通道复用控制信号，依次将接收的延迟的N个频率的GNSS信号以及未延迟的N个频率的GNSS信号整合到一个通道中，其中，每个频率的GNSS信号均包含延迟与未延迟的两路信号，并按照频率的不同，分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块，直到N个不同频率的包含两路信号的GNSS信号全部被发送至所述加窗处理模块。
S405、所述加窗处理模块以及FFT计算模块依次对接收的延迟以及未延迟的两路GNSS信号进行加窗处理以及FFT计算处理，并将得到的计算结果发送给所述NBI滤波处理模块，所述NBI滤波处理模块根据所述分时复用控制模块发送的通道NBI滤波配置信号对滤波参数进行配置，进而对两路GNSS信号进行滤波处理，并将滤波结果发送给所述IFFT计算模块，所述IFFT计算模块对两路GNSS信号进行IFFT计算，并将计算之后的两路GNSS信号发送给所述解复用模块，同时向所述分时复用控制模块发送IFFT完成指示信号，所述分时复用控制模块根据所述完成指示信号向所述解复用控制模块发送多通道解复用控制信号。
具体的，所述加窗处理模块分时段重复使用，对N个频率的GNSS信号分时段分别进行加窗处理，在处理时，对依次接收的延迟以及未延迟的两路GNSS信号分别乘以在加窗函数系数ROM表中查找得到的同一加窗系数，进行加窗处理，其中，所述加窗函数具体包括矩形窗、三角窗、汉宁窗、海明窗以及布莱克曼窗，所述加窗函数均为对称函数。
其中，由于所述缓存处理模块按照FFT蝶形运算顺序读出当前GNSS信号，所以所述延迟模块以及加窗处理模块均按照FFT蝶形运算顺序进行处理，若所述缓存处理模块在RF时钟域下写入M个数据采样点的顺序分别为D0、D1、D2、…、DM-1，那么在系统时钟域下读出的顺序为D0、DM/2、D1、DM/2+1、D2、DM/2+2、…、DM/2-1、DM-1，进而所述加窗处理模块从加窗函数ROM表中读出对应的M个加窗系数的顺序为F0、FM/2、F1、FM/2+1、F2、FM/2+2、…、FM/2-1、FM-1,加窗函数系数ROM表中只需存储M/2个加窗系数，加窗函数系数ROM表大小只需M/2。
所述缓存处理模块按FFT蝶形运算顺序读出GNSS信号，并且所述延迟模块以及加窗处理模块都按照FFT蝶形运算顺序进行处理，将FFT计算模块前的“倒位序”处理去掉，共节省M个数据存储单元(每一路GNSS信号节省M/2个数据存储单元)。
所述加窗处理模块将加窗处理的结果发送给所述FFT计算模块，所述FFT计算模块分时段重复使用，对N个频率的GNSS信号分时段分别进行FFT计算处理，具体的，所述FFT计算模块将经过加窗处理的两路GNSS信号中的一路GNSS信号实部输入进行FFT计算，另一路GNSS信号虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离得到两路实数FFT信号，将两路GNSS信号从时域数据转换到频域数据，实现FFT计算模块关于两路GNSS信号的共享复用，并将得到的计算结果发送给所述NBI滤波处理模块。
所述NBI滤波处理模块分时段重复使用，根据所述分时复用控制模块发送的针对N个不同频率GNSS信号的通道NBI滤波配置信号对滤波参数进行配置，对经过FFT计算的N个不同频率的两路频域GNSS信号分时段进行滤波处理，所述NBI滤波处理模块可采用各种滤波抑制算法，例如门限检测法、K谱线算法以及中值滤波法等，所述NBI滤波处理模块将滤波处理的结果发送给所述IFFT计算模块。
所述IFFT计算模块分时段重复使用,对N个频率的GNSS信号分时段分别进行IFFT计算处理，具体的，所述IFFT计算模块将经过滤波处理的两路GNSS信号中的一路GNSS信号实部输入进行IFFT计算，另一路GNSS信号虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离得到两路实数IFFT信号，将两路GNSS信号从频域数据转换到时域数据，然后将得到的计算结果发送给所述解复用模块，同时还向所述分时复用控制模块发送IFFT完成指示信号，所述分时复用控制模块根据所述完成指示信号向所述解复用控制模块发送多通道解复用控制信号。
S406、所述解复用模块根据所述多通道解复用控制信号，使分时段接收的N个两路GNSS信号按照频率的不同处在N个通道中，并分别发送给对应连接的N个所述叠加处理模块,N个所述叠加处理模块分别将接收的两路时域GNSS信号进行叠加处理，并将得到的一路时域GNSS信号分别返回给对应连接的N个所述缓存处理模块以及预处理模块输出，以进行后续解扩处理。
具体的，所述解复用模块根据所述分时复用控制模块发送的所述多通道解复用控制信号，对分时段接收到的处于一个通道中的N个频率的包含两路信号的GNSS信号进行解复用，使N个两路GNSS信号按照频率的不同处在N个通道中，并分别发送给对应连接的N个所述叠加处理模块，N个所述叠加处理模块分别将接收的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，并分别返回给对应连接的N个所述缓存处理模块，N个所述缓存处理模块分别将N个频率的GNSS信号读出返回给对应的所述预处理模块，其中，所述缓存处理模块是在所述缓存处理模块写入下一个GNSS信号之前，提前一个RF时钟，读出经叠加处理之后的GNSS信号；N个所述预处理模块分别将N个频率的GNSS信号输出，以进行后续解扩处理。
另外，在所有频率的GNSS信号都不需要进行窄带干扰抑制时，所述分时复用控制模块控制所有的延迟模块、所有的叠加处理模块、复用模块、加窗模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块全部关闭。
在实际接入的GNSS信号的频率个数为K时，所述分时复用控制模块调整系统时钟频率为RF时钟频率的K+1倍，这样若RF时钟域输入M个采样点数据的时间为T，在系统时钟域内读出M个采样点数据时间只需T/(K+1)。
N个不同频率的GNSS信号依次在所述延迟模块、加窗处理模块、FFT计算模块、NBI滤波处理模块以及IFFT计算模块中处理所需总时间是固定的，例如，设GNSS信号在一个RF时钟域下所需时间为T，所述延迟模块在系统时钟域下将N个不同频率的GNSS信号延迟M/2个采样点所需时间为T/2(N+1)，N个不同频率的GNSS信号以及延迟M/2个采样点的GNSS信号在系统时钟域下从加窗处理模块到NBI滤波处理模块依次处理所需总时间为N*T/(N+1)，所以总的处理时间约为(2N+1)*T/2(N+1)。
本发明提供的一种单频及多频GNSS信号窄带干扰抑制装置、系统及其方法，实现了对单频率GNSS信号以及多频率的GNSS信号的窄带干扰抑制。其中，所述FFT计算模块与IFFT计算模块均设有实部输入与虚部输入，能够同时对两路GNSS信号进行计算，实现了FFT与IFFT计算模块的共享复用，使FFT/IFFT计算模块的硬件成本降低一半；所述加窗处理模块能够对两路GNSS信号进行加窗处理，共享一个窗函数系数ROM表，并根据窗函数的对称性只存储M/2个加窗函数系数，如此将存储窗函数系数的ROM资源节省为原来的1/4；所述NBI滤波处理模块同时对两路GNSS信号进行滤波处理，并能根据不同频率的GNSS信号配置不同的滤波参数，保持滤波的独立性，有效抑制窄带干扰；在对N个频率的GNSS信号进行窄带干扰抑制时，所述加窗处理模块、FFT计算模块、NBI滤波处理模块以及IFFT计算模块的采用分时复用的机制，将硬件资源降低为原来的1/N；所述缓存处理模块将GNSS信号按照FFT蝶形运算顺序读出，并且所述延迟模块以及加窗处理模块均按照FFT蝶形运算顺序处理GNSS信号，将FFT计算模块前的“倒位序”去掉，节省M个数据存储单元；所述分时复用控制单元能够根据输入的GNSS信号的频率的数量自动调整系统时钟频率，降低系统功耗；另外所述预处理模块能够对GNSS信号进行干扰信号检测，实现对GNSS信号的自适应窄带干扰抑制处理，并且在所有的GNSS信号都不需要进行窄带干扰抑制的情况下，关闭所述延迟模块、叠加处理模块、复用模块、加窗模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块，降低系统功耗。本发明能够对单频率以及多频率的GNSS信号进行窄带干扰抑制，不仅具有良好的抑制效果，还能够降低硬件资源成本以及系统的功耗，并具有较强的实用性。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种单频GNSS信号窄带干扰抑制装置，其特征在于，所述装置包括延迟模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及叠加处理模块；
所述延迟模块，用于接收一路GNSS信号，并将其延迟M/2个采样点输出给所述加窗处理模块，其中M为FFT块长度，具体为2的指数幂；
所述加窗处理模块，用于接收一路GNSS信号，以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，并对两路时域GNSS信号分别进行加窗处理；
所述FFT计算模块，包括实部输入与虚部输入，用于对经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号；
所述NBI滤波处理模块，用于对经过FFT计算的两路频域GNSS信号进行滤波处理；
所述IFFT计算模块，包括实部输入与虚部输入，用于对经过滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号，并将其输出；
所述叠加处理模块，用于对接收的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，并将其输出。
2.如权利要求1所述的一种单频GNSS信号窄带干扰抑制装置，其特征在于，所述加窗处理模块具体用于：
接收一路GNSS信号，以及所述延迟模块发送的延迟M/2个采样点的GNSS信号，并将两路时域GNSS信号分别乘以在加窗函数的系数ROM表中查找得到的同一加窗系数，进行加窗处理；
所述加窗函数具体包括矩形窗、三角窗、汉宁窗、海明窗以及布莱克曼窗，所述加窗函数均为对称函数，若FFT的块长度为M，对应于加窗函数的M个系数分别为F1,F2,…,FM，则有Fi＝FM+1-i(i＝1,2,…,M/2-1)，加窗函数系数ROM表中只需存储M/2个加窗系数。
3.一种多频GNSS信号窄带干扰抑制系统，其特征在于，所述系统包括权利要求1或2所述的单频GNSS信号窄带干扰抑制装置以及至少一个预处理模块、至少一个缓存处理模块、至少一个缓存控制模块、复用模块、解复用模块以及分时复用控制模块；
所述预处理模块，包括格式转换单元以及信号检测单元，用于在RF时钟域下对接收到的GNSS信号进行格式转换，并对格式转换后的GNSS信号进行干扰信号检测，以确定所述GNSS信号是否需要进行窄带干扰抑制处理，且根据确定结果生成通道干扰抑制使能控制信号并发送给所述缓存控制模块以及将检测之后的GNSS信号输出给所述缓存处理模块；还用于将所述缓存处理模块返回的GNSS信号输出，以进行后续解扩处理；
所述缓存处理模块，用于在RF时钟域下写入干扰信号检测之后的GNSS信号，然后在所述缓存控制模块的控制下读出所述GNSS信号，并根据GNSS信号是否需要进行窄带干扰抑制处理的确定结果，将读出的GNSS信号返回给对应连接的所述预处理模块输出或者发送给对应连接的所述延迟模块以及所述复用模块；
所述缓存控制模块，用于在系统时钟域下向所述分时复用控制模块发送所述通道干扰抑制使能控制信号以及缓存数据就绪指示信号，并接收所述分时复用控制模块反馈的通道缓存读使能信号，进而控制所述缓存处理模块工作；
所述复用模块，用于根据所述分时复用控制模块发送的多通道复用控制信号，将依次接收的多个所述缓存处理模块发送的未延迟的多个频率的GNSS信号，以及多个所述延迟模块发送的延迟的多个频率的GNSS信号整合到一个通道中，并按照频率的不同，分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给所述加窗处理模块；
所述解复用模块，用于根据所述分时复用控制模块发送的多通道解复用控制信号，对分时段接收的处于一个通道中的多个包含两路信号的不同频率的GNSS信号进行解复用，使多个两路GNSS信号按照频率的不同处在多个通道中，并分别发送给多个所述叠加处理模块。
4.如权利要求3所述的一种多频GNSS信号窄带干扰抑制系统，其特征在于，所述缓存处理模块具体用于：
在RF时钟域下写入干扰信号检测之后的GNSS信号，并根据所述缓存控制模块的控制信号读出当前GNSS信号：在不需要对GNSS信号进行窄带干扰抑制时，在RF时钟域下顺序读出GNSS信号并返回给对应连接的所述预处理模块输出；在需要对GNSS信号进行窄带干扰抑制时，在系统时钟域下按照FFT蝶形运算顺序读出GNSS信号并发送给对应连接的所述延迟模块以及所述复用模块；在读出当前GNSS信号的同时，在RF时钟域下继续写入下一个经干扰信号检测之后的GNSS信号，实现乒乓操作。
5.如权利要求3所述的一种多频GNSS信号窄带干扰抑制系统，其特征在于，所述分时复用控制模块，用于在所有的GNSS信号都不需要进行窄带干扰抑制时，控制所有的延迟模块、所有的叠加处理模块、复用模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块全部关闭；以及用于在实际接入的GNSS信号的频率个数为K时，调整系统时钟频率为RF时钟频率的K+1倍，其中K为正整数。
6.一种单频GNSS信号窄带干扰抑制方法，其特征在于，所述方法包括：
S301、将单频率的时域GNSS信号分为两路，其中一路信号延迟M/2个采样点，另一路信号不延迟，并将延迟后的以及未延迟的两路GNSS信号分别进行加窗处理；
S302、将经过加窗处理的两路GNSS信号中的一路实部输入进行FFT计算，另一路虚部输入进行FFT计算，并将计算得到的复数FFT信号进行分离计算得到两路实数FFT信号；
S303、将经过FFT计算的两路频域GNSS信号进行滤波处理，并将经滤波处理的两路频域GNSS信号中的一路实部输入进行IFFT计算，另一路虚部输入进行IFFT计算，并将计算得到的复数IFFT信号进行分离计算得到两路实数IFFT信号；
S304、将经过IFFT计算的两路时域GNSS信号进行叠加处理，得到一路时域GNSS信号，以进行后续的解扩处理。
7.如权利要求6所述的一种单频GNSS信号窄带干扰抑制方法，其特征在于，所述S301中所述将两路GNSS信号进行加窗处理具体包括：
将延迟M/2个采样点的GNSS信号以及未延迟的GNSS信号，分别乘以在加窗函数系数ROM表中查找得到的同一加窗系数，进行加窗处理；
所述加窗函数具体包括矩形窗、三角窗、汉宁窗、海明窗以及布莱克曼窗，所述加窗函数均为对称函数，若FFT的块长度为M，对应于加窗函数的M个系数分别为F1,F2,…,FM，则有Fi＝FM+1-i(i＝1,2,…,M/2-1)，加窗函数系数ROM表中只需存储M/2个加窗系数。
8.一种多频GNSS信号窄带干扰抑制方法，其特征在于，所述方法包括：
S401、N个预处理模块在RF时钟域下，对接收的N个GNSS信号分别进行预处理，并将预处理之后的GNSS信号输出给对应连接的N个缓存处理模块，同时还通过N个缓存控制模块分别将通道干扰抑制使能控制信号发送给分时复用控制模块；
S402、在RF时钟域下N个所述缓存处理模块在写入地址依次相差M/4的情况下分别写入N个GNSS信号，在依次写满之后，其对应的N个所述缓存控制模块依次发送缓存数据就绪指示信号给所述分时复用控制模块；所述分时复用控制模块根据所述缓存数据就绪指示信号以及之前接收的所述通道干扰抑制使能控制信号，生成通道缓存读使能信号并依次反馈给N个所述缓存控制模块；
S403、N个所述缓存控制模块根据所述缓存读使能信号分别控制N个所述缓存处理模块读出N个GNSS信号，若不需要对GNSS信号进行窄带干扰抑制，将GNSS信号返回给对应连接的所述预处理模块输出以进行后续解扩处理，结束当前GNSS信号的抑制流程；若需要对GNSS信号进行窄带干扰抑制，将GNSS信号发送给对应连接的延迟模块以及复用模块；
S404、N个所述延迟模块在系统时钟域下，依次分别将接收的GNSS信号延迟M/2个采样点，并将延迟后的GNSS信号发送给所述复用模块，所述复用模块根据所述分时复用控制模块发送的多通道复用控制信号，依次将接收的延迟的N个频率的GNSS信号以及未延迟的N个频率的GNSS信号整合到一个通道中，并按照频率的不同分时段依次发送频率相同的延迟以及未延迟的两路GNSS信号给加窗处理模块，直到N个不同频率的包含两路信号的GNSS信号全部被发送至所述加窗处理模块；
S405、所述加窗处理模块以及FFT计算模块依次对接收的延迟以及未延迟的两路GNSS信号进行加窗处理以及FFT计算处理，并将得到的计算结果发送给NBI滤波处理模块，所述NBI滤波处理模块根据所述分时复用控制模块发送的通道NBI滤波配置信号对滤波参数进行配置，进而对两路GNSS信号进行滤波处理，并将滤波结果发送给IFFT计算模块，所述IFFT计算模块对两路GNSS信号进行IFFT计算，并将计算之后的两路GNSS信号发送给解复用模块，同时向所述分时复用控制模块发送IFFT完成指示信号，所述分时复用控制模块根据所述完成指示信号向所述解复用控制模块发送多通道解复用控制信号；
S406、所述解复用模块根据所述多通道解复用控制信号，使分时段接收的N个两路GNSS信号按照频率的不同处在N个通道中，并分别发送给对应连接的N个叠加处理模块,N个所述叠加处理模块分别将接收的两路时域GNSS信号进行叠加处理，并将得到的一路时域GNSS信号分别返回给对应连接的N个所述缓存处理模块以及预处理模块输出，以进行后续解扩处理。
9.如权利要求8所述的一种多频GNSS信号窄带干扰抑制方法，其特征在于，所述步骤S403具体包括：
N个所述缓存控制模块根据所述缓存读使能信号分别控制N个所述缓存处理模块读出N个当前GNSS信号：在不需要对GNSS信号进行窄带干扰抑制时，在RF时钟域下顺序读出GNSS信号并返回给对应连接的所述预处理模块输出；在需要对GNSS信号进行窄带干扰抑制时，在系统时钟域下按照FFT蝶形运算顺序读出GNSS信号并发送给对应连接的所述延迟模块以及所述复用模块；N个所述缓存控制模块在读出当前GNSS信号的同时，依次分别在RF时钟域下继续写入下一个经预处理之后的GNSS信号，实现乒乓操作；
所述延迟处理模块以及加窗处理模块均按照FFT蝶形运算顺序对接收的GNSS信号进行处理。
10.如权利要求8所述的一种多频GNSS信号窄带干扰抑制方法，其特征在于，所述方法还包括：
在所有的GNSS信号都不需要进行窄带干扰抑制时，所述分时复用控制模块控制所有的延迟模块、所有的叠加处理模块、复用模块、加窗处理模块、FFT计算模块、NBI滤波处理模块、IFFT计算模块以及解复用模块全部关闭；
在实际接入的GNSS信号的频率个数为K时，所述分时复用控制模块调整系统时钟频率为RF时钟频率的K+1倍，其中，K为正整数。
CN201510810632.9A 2015-11-20 2015-11-20 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法 Active CN105301607B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510810632.9A CN105301607B (zh) 2015-11-20 2015-11-20 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510810632.9A CN105301607B (zh) 2015-11-20 2015-11-20 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法

Publications (2)

Publication Number Publication Date
CN105301607A CN105301607A (zh) 2016-02-03
CN105301607B true CN105301607B (zh) 2017-08-08

Family

ID=55199085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510810632.9A Active CN105301607B (zh) 2015-11-20 2015-11-20 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法

Country Status (1)

Country Link
CN (1) CN105301607B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111399001B (zh) * 2020-03-31 2022-05-06 和芯星通科技(北京)有限公司 一种处理宽带干扰方法和装置
CN113253315B (zh) * 2021-07-07 2021-10-22 湖南跨线桥航天科技有限公司 一种基于时分复用技术的导航信号相关峰监测装置及方法
CN113595587B (zh) * 2021-08-16 2022-07-08 中国电子科技集团公司第五十四研究所 一种用于变换域干扰抑制的自适应幅度恢复系统
CN116049612B (zh) * 2023-04-03 2023-07-28 南京楚航科技有限公司 一种基于zynq-fpga的fft复用方法、系统及装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201910793U (zh) * 2010-12-29 2011-07-27 成都国星通信有限公司 扩频接收机前端频域滤波抗干扰装置
CN102073621B (zh) * 2011-01-19 2012-08-22 东南大学 一种应用于fft/ifft的基4蝶形单元电路及其处理方法
WO2013090568A1 (en) * 2011-12-15 2013-06-20 Northrop Grumman Guidance And Electronics Company, Inc. System and method for detection of rf signal spoofing
CN102944884B (zh) * 2012-11-02 2015-08-05 锐迪科科技有限公司 Gnss接收机检测并消除窄带干扰的方法
CN203894413U (zh) * 2013-12-16 2014-10-22 重庆九洲星熠导航设备有限公司 一种基于fpga的卫星导航信号抗窄带干扰处理机
CN104749591A (zh) * 2015-04-09 2015-07-01 中国人民解放军国防科学技术大学 面向全球卫星导航系统的多模并行射频接收机
CN104898132A (zh) * 2015-05-14 2015-09-09 电子科技大学 一种阈值处理和空频自适应算法结合的导航抗干扰算法

Also Published As

Publication number Publication date
CN105301607A (zh) 2016-02-03

Similar Documents

Publication Publication Date Title
CN109212559A (zh) 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法
CN105301607B (zh) 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法
CN107461611B (zh) 基于小波和emd重构相结合的泄漏检测方法及泄漏检测装置
CN102944884B (zh) Gnss接收机检测并消除窄带干扰的方法
CN203894413U (zh) 一种基于fpga的卫星导航信号抗窄带干扰处理机
CN109164483A (zh) 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置
CN104360355B (zh) 抗干扰方法和装置
CN103913778B (zh) 多个近端参考线圈的核磁共振信号实时噪声抵消装置
CN104898132A (zh) 一种阈值处理和空频自适应算法结合的导航抗干扰算法
CN102053276A (zh) 一种地震数字信号的复数道集二维滤波方法
CN104076369A (zh) 基于自适应门限判决的频域抗干扰方法及装置
CN205157788U (zh) 单频gnss信号窄带干扰抑制装置及多频gnss信号窄带干扰抑制系统
CN203465417U (zh) 基于自适应噪声抵消器的声纳阵列信号处理装置
Priputin et al. Implementation of a Discrete Frequency Demodulator Using Filter Bank Technologies on an ARM Processor
CN101345548A (zh) 基于可调流水线频率变换tpft的抗窄带干扰的方法
CN107809256A (zh) 一种在窄带干扰下的短波抑制方法
CN109474356A (zh) 宽带多通道信号能量检测系统及方法
US9705608B2 (en) Method and system for interference cancellation of data channel
Carter Time delay estimation
CN112489669A (zh) 一种音频信号处理方法、装置、设备和介质
CN102571031B (zh) 一种用于实现小波变换的宽带低插损声表面波滤波器组
CN104483679A (zh) 一种gnss数字中频信号窄带干扰抑制方法及相应的电路
Yao et al. Iterative threshold algorithm for narrow-band interference suppression
CN111812682A (zh) 一种抗窄带干扰电路
Raj et al. An overview of MDCT for Time Domain Aliasing Cancellation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant