CN110286392B - 一种基于降采样的频域抗窄带干扰实现方法 - Google Patents
一种基于降采样的频域抗窄带干扰实现方法 Download PDFInfo
- Publication number
- CN110286392B CN110286392B CN201910701608.XA CN201910701608A CN110286392B CN 110286392 B CN110286392 B CN 110286392B CN 201910701608 A CN201910701608 A CN 201910701608A CN 110286392 B CN110286392 B CN 110286392B
- Authority
- CN
- China
- Prior art keywords
- signals
- frequency domain
- signal
- sampling
- band interference
- 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/015—Arrangements for jamming, spoofing or other methods of denial of service of such systems
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Noise Elimination (AREA)
Abstract
本发明公开了一种基于降采样的频域抗窄带干扰实现方法,属于卫星导航通信领域。本发明将利用变频方法实现信号的降采样,然后利用降采样后的多路信号对抗窄带干扰模块的复用,同时通过降采样的方案实现在同样资源下的频域抗窄带干扰方法的频率分辨率。该方法克服了现有频域抗窄带干扰方法资源占用率高,频率分辨率差,后端信号恶化严重的问题。本发明中的降采样处理和分时复用方法,极大节省了硬件实现资源。本发明能够在传统处理方法同样资源上处理多路信号,并且频率分辨率至少提高2倍,性能远远优于传统方法。
Description
技术领域
本发明涉及一种基于降采样的频域抗窄带干扰实现方法,属于卫星导航通信中阵列信号处理领域。
背景技术
卫星导航系统作为全球的空间信息基础设施,为全球各类用户提供时间空间基准和所有与位置相关的实时动态信息,是国家重大空间和信息化基础设施。同时导航信号也同样存在功率低,易受到各种各样电磁干扰的影响,使得导航服务能力极大降低的缺陷。窄带干扰信号通常是指带宽为信号带宽10%以下的干扰信号,有着易生成,发射功率强的特点,是当前影响导航性能下降的一个重要因素,所以如何有效抑制窄带干扰对导航信号的影响,是当前导航领域的一个关键课题。
传统频域抗干扰算法是将AD采样信号加窗处理后直接通过傅立叶变换由时域转化到频域进行处理,通过设置合适阈值,将检测到的窄带干扰信号滤除,然后在变换回时域输出。基于降采样的频域抗窄带干扰方法是对传统频域抗干扰算法的一种优化,通过降采样的方法对信号进行二次处理,然后采用先存储后处理的方案,通过这种处理将处理信号带宽下降到理想范围内,提高变换后的频域分辨率,提升频域抗干扰能力。同时通过不同信号之间分时协调处理,同时在同一模块中处理多路导航信号,极大的提高资源利用率,降低硬件使用资源。相比传统频域抗干扰算法,基于降采样的优化算法具有抗干扰性能好,资源利用率高的优点。
发明内容
本发明的目的是为了克服传统频域抗干扰能力差,资源占用率高的问题,提出基于降采样的频域抗干扰实现算法,普遍适用于当前卫星导航各类接收机终端。
本发明技术方案的基本思路是:在传统频域抗窄带干扰的基础上,在信号输入前端进行变频和降采样处理,然后将信号进行缓存后进行频域抗窄带处理,达到提高抗干扰性能和资源复用的目的。
本发明采用的技术方案为:
一种基于降采样的频域抗窄带干扰实现方法,包括如下步骤:
步骤1,将每一路AD中频采样信号分别进行二次变频处理,得到正交数字信号;
步骤2,将每一路正交数字信号分别进行降采样处理,将降采样处理后的信号进行信号加窗,然后将加窗后的信号由时域变换为频域,在频域中计算出阈值,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,将剥离窄带干扰后的有用信号变换到时域中;
步骤3,将每一路无窄带干扰的时域信号分别输出给接收机终端进行捕获跟踪解调。
其中,步骤2具体包括如下步骤:
(202)将每一路降采样后的信号分别采用乒乓处理的方式进行缓存;
(203)控制相应时序,将缓存里的每路信号及延迟N/2点的信号分别进行加窗处理,加窗后的信号表示为:
X3(nt)=X2(nt)*windows(nt)
X'3(nt)=X2(nt-N/2)*windows(nt)
(204)将每一组加窗后的信号分别进行傅立叶变化,得到一组傅立叶变化后的频域信号;
(205)每一组傅立叶变化后的频域信号中每个频域信号分别计算一个阈值,计算阈值公式为当频域信号幅值大于阈值时信号正常通过,频域信号幅值小于阈值时信号幅值钳制为零,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,公式表示为F(n)表示傅里叶变化后的频域信号;
(206)将步骤(205)处理后的多组无窄带干扰的信号分别通过反傅立叶变换变换到时域中;
(207)通过时序控制,将多组反傅立叶变换后的时域信号分别缓存到不同的缓存器中。
本发明与现有技术相比所取得的有益效果为:
本发明提出的一种基于降采样的频域抗窄带干扰实现方法,算法采用降采样和时分复用技术,克服了传统频域抗窄带干扰频率分辨率低,资源占用率高和使用率低的缺点,有效减少了抗干扰对后端接收机终端的影响,极大的节省了硬件资源。
附图说明
图1为本发明的基于降采样的频域抗窄带干扰流程图;
图2为本发明的优化后的频域抗窄带干扰流程图。
具体实施方式
下面结合附图对本发明作进一步解释说明。
一种基于降采样的频域抗窄带干扰实现方法,如图1所示,包括如下步骤:
步骤1,对接收信号进行变频处理
将每一路AD中频采样信号进行二次变频处理,得到路正交数字信号接收信号。
步骤2,对变频信号进行频域抗干扰处理
将每一路正交数字信号分别进行降采样处理,将降采样处理后的信号进行信号加窗,然后将加窗后的信号由时域变换为频域,在频域中计算出阈值,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,将剥离窄带干扰后的有用信号变换到时域中;具体为:
第(1)步,进行降采样处理
第(2)步,降采样数据缓存
将每一路降采样后的信号分别采用乒乓处理的方式进行缓存;缓存器输入时钟为降采样后的工作时钟,输出时钟可以是AD信号采样时钟Fs或者更高工作时钟。
第(3)步,数据加窗处理,如图2;
控制相应时序,将缓存里的每路信号及延迟N/2点的信号分别进行加窗处理,加窗的目的是为了抑制后续时域/频域变换时的频谱展宽问题,使得变换后的频域更接近真实频谱,加窗后的信号表示为:
X3(nt)=X2(nt)*windows(nt)
X'3(nt)=X2(nt-N/2)*windows(nt)
第(4)步,时域频域变换
将每一组加窗后的信号分别输出给傅立叶变化模块,输入信号为复数数据,复数数据的实部和虚部分别为正交两路信号。得到一组傅立叶变化后的频域信号;
第(5)步,阈值计算
每一组傅立叶变化后的频域信号中每个频域信号分别计算一个阈值,阈值当信号幅值大于阈值时信号正常通过,信号幅值小于阈值时信号幅值钳制为零,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,公式表示为F(n)表示傅里叶变化后的频域信号;
第(6)步,频域时域变换
将处理后的多组无窄带干扰的信号分别通过反傅立叶变换变换到时域中;可以选择进行去窗处理,也可以省去窗步骤。
第(7)步,时域数据缓存
通过时序控制,将多路处理信号分别缓存到不同的缓存器中,缓存器输入时钟为前端缓存器输出时钟,缓存器输出时钟为降采样的采样时钟。通过输入输出时钟不一致性,保障了输出数据的连续性和处理模块的复用性。
步骤3,将抗干扰处理后的数据输出给接收机终端。采用时分复用技术,可以处理多路信号。
将经过频域处理后无窄带干扰的信号输出,此时各路时域信号为连续信号,并且输出频率不尽相同,每路信号的工作频率与前端降采样后采样频率相同。将处理后的信号输出给接收机终端捕获跟踪解调。
虽然结合附图描述了本发明的实施方式,但本发明的保护范围并不局限于此,对于本领域技术人员来说,在不脱离本发明原理的前提下,还可以做出若干变形和改进,这些也应视为属于本发明的保护范围。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (1)
1.一种基于降采样的频域抗窄带干扰实现方法,其特征在于,包括如下步骤:
步骤1,将每一路AD中频采样信号分别进行二次变频处理,得到正交数字信号;
步骤2,将每一路正交数字信号分别进行降采样处理,将降采样处理后的信号进行信号加窗,然后将加窗后的信号由时域变换为频域,在频域中计算出阈值,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,将剥离窄带干扰后的有用信号变换到时域中;
步骤3,将每一路无窄带干扰的时域信号分别输出给接收机终端进行捕获跟踪解调;
其中,步骤2包括如下步骤:
(202)将每一路降采样后的信号分别采用乒乓处理的方式进行缓存;缓存器输入时钟为降采样后的工作时钟,输出时钟是AD信号采样时钟或者更高工作时钟;
(203)控制相应时序,将缓存里的每路信号及延迟N/2点的信号分别进行加窗处理,加窗后的信号表示为:
X3(nt)=X2(nt)*windows(nt)
X′3(nt)=X2(nt-N/2)*windows(nt)
(204)将每一组加窗后的信号分别进行傅立叶变化,得到一组傅立叶变化后的频域信号;
(205)每一组傅立叶变化后的频域信号中每个频域信号分别计算一个阈值,计算阈值公式为当频域信号幅值大于阈值时信号正常通过,频域信号幅值小于阈值时信号幅值钳制为零,通过将频域信号与阈值比较,将窄带干扰信号从有用信号频谱中剥离,公式表示为F(n)表示傅里叶变化后的频域信号;
(206)将步骤(205)处理后的多组无窄带干扰的信号分别通过反傅立叶变换变换到时域中;
(207)通过时序控制,将多组反傅立叶变换后的时域信号分别缓存到不同的缓存器中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910701608.XA CN110286392B (zh) | 2019-07-31 | 2019-07-31 | 一种基于降采样的频域抗窄带干扰实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910701608.XA CN110286392B (zh) | 2019-07-31 | 2019-07-31 | 一种基于降采样的频域抗窄带干扰实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110286392A CN110286392A (zh) | 2019-09-27 |
CN110286392B true CN110286392B (zh) | 2021-06-22 |
Family
ID=68024488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910701608.XA Active CN110286392B (zh) | 2019-07-31 | 2019-07-31 | 一种基于降采样的频域抗窄带干扰实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110286392B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111812682A (zh) * | 2020-07-24 | 2020-10-23 | 华力智芯(成都)集成电路有限公司 | 一种抗窄带干扰电路 |
CN115604066B (zh) * | 2022-10-08 | 2024-05-03 | 杭州万高科技股份有限公司 | 一种基于离散傅立叶变换的低成本窄带干扰频率估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483679A (zh) * | 2014-12-02 | 2015-04-01 | 泰斗微电子科技有限公司 | 一种gnss数字中频信号窄带干扰抑制方法及相应的电路 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101345548B (zh) * | 2008-08-22 | 2011-12-28 | 清华大学 | 基于可调流水线频率变换tpft的抗窄带干扰的方法 |
US9585621B2 (en) * | 2013-03-25 | 2017-03-07 | Iliya Mitov | Technique for real-time removal of power line interference in ECG |
US9214972B2 (en) * | 2013-06-13 | 2015-12-15 | Qualcomm Technologies International, Ltd. | Method and apparatus for on-demand interference rejection in multi-band GNSS receivers |
FR3012620A1 (fr) * | 2013-10-31 | 2015-05-01 | Thales Sa | Procede de detection d'interferences dans un signal de radio-navigation par satellite par detection d'une deformation de la fonction de correlation |
CN109212559A (zh) * | 2018-09-26 | 2019-01-15 | 永州市金蚂蚁新能源机械有限公司 | 一种单频及多频gnss信号窄带干扰抑制装置、系统及其方法 |
-
2019
- 2019-07-31 CN CN201910701608.XA patent/CN110286392B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483679A (zh) * | 2014-12-02 | 2015-04-01 | 泰斗微电子科技有限公司 | 一种gnss数字中频信号窄带干扰抑制方法及相应的电路 |
Non-Patent Citations (1)
Title |
---|
数字下变频中抽取滤波器的设计及FPGA实现;周云等;《微电子技术》;20151231;第45-50页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110286392A (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107612865B (zh) | 一种应用于低压电力线载波通信的信号降噪方法 | |
CN110286392B (zh) | 一种基于降采样的频域抗窄带干扰实现方法 | |
CN102437865B (zh) | 提高直扩跳频混合扩频系统频域抗干扰性能的方法 | |
CN114362837B (zh) | 一种扩频卫星信号窄带干扰自适应消除装置 | |
CN116299577B (zh) | 一种应用于北斗三号基带芯片的窄带抗干扰装置及方法 | |
CN113447893B (zh) | 一种雷达脉冲信号频谱自动检测方法、系统及介质 | |
CN108923877B (zh) | 一种pcma定时捕获和跟踪方法 | |
CN107171693A (zh) | 基于波形重构的宽带干扰抑制方法 | |
CN113253313A (zh) | 一种改进pmf-fft系统多普勒频率估计精度的方法 | |
CN114070702B (zh) | 一种基于叠加与快速频域相关联合的pss检测方法 | |
CN112003632A (zh) | 一种直接序列扩频通信多址干扰抑制办法 | |
CN103560799B (zh) | 一种有效抑制中等强度测距仪脉冲干扰的正交频分复用接收机装置 | |
CN112968850A (zh) | 微弱bpsk信号载波捕获方法 | |
CN116527081A (zh) | 一种适用于中压载波系统的脉冲噪声抑制方法 | |
CN100553243C (zh) | 波束空间正交频分复用调制系统自适应波束形成方法 | |
CN108401582B (zh) | 一种适应tdma/fdma体制的geo卫星移动通信系统初始频率同步方法 | |
CN114615113B (zh) | 一种利用前导序列对基系数降噪的多符号bem信道估计方法 | |
CN1845538B (zh) | 一种时域信道同步跟踪器 | |
CN102801663B (zh) | 一种用于深空探测器的多档码速率估计方法 | |
CN112600784B (zh) | 一种基于二次差分相关的大频偏位同步方法 | |
CN205157788U (zh) | 单频gnss信号窄带干扰抑制装置及多频gnss信号窄带干扰抑制系统 | |
CN108540411B (zh) | 一种利用快速傅里叶变换的低复杂度bem信道估计方法 | |
CN113406672A (zh) | 一种基于自适应加窗的分段线性调频干扰消除方法 | |
CN103346827B (zh) | 稳健gsc波束形成器及其波束优化方法 | |
Ji et al. | FPGA implementation of SC-FDE timing synchronization algorithm |
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 |