CN105911565A - 一种抑制窄带干扰的方法及装置 - Google Patents
一种抑制窄带干扰的方法及装置 Download PDFInfo
- Publication number
- CN105911565A CN105911565A CN201610215072.7A CN201610215072A CN105911565A CN 105911565 A CN105911565 A CN 105911565A CN 201610215072 A CN201610215072 A CN 201610215072A CN 105911565 A CN105911565 A CN 105911565A
- Authority
- CN
- China
- Prior art keywords
- signal
- initial
- panel
- delay
- carried out
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 261
- 230000008569 process Effects 0.000 claims abstract description 196
- 230000001629 suppression Effects 0.000 claims abstract description 45
- 230000001934 delay Effects 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 35
- 230000003595 spectral effect Effects 0.000 claims description 20
- 230000009466 transformation Effects 0.000 claims description 11
- 238000012935 Averaging Methods 0.000 claims description 5
- 239000004744 fabric Substances 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 230000003044 adaptive effect Effects 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 7
- 230000005764 inhibitory process Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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
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)
- Noise Elimination (AREA)
Abstract
本申请公开了一种抑制窄带干扰的方法及装置，其中，所述方法包括对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；将所述各路处理信号进行叠加处理，获得输出信号。由于在所述干扰抑制处理的过程中采用了最大旁瓣低于布莱克曼窗函数最大旁瓣的预设窗函数，因此所述方法对窄带干扰的抑制效果较好；并且由于所述方法采用至少四路信号同时处理的措施避免了所述输出信号的信噪比损失较大而导致卫星导航定位系统性能急剧下降的情况出现。
Description
技术领域
本发明涉及抗干扰技术领域，更具体地说，涉及一种抑制窄带干扰的方法及装置。
背景技术
卫星导航定位系统广泛应用于生活中的各个领域中，卫星导航定位系统采用扩频通信技术，其本身能够抑制环境中一定强度的干扰信号，但是当卫星导航定位系统受到人为的较强的窄带干扰时，整个系统的性能会急剧下降。
现有技术中主流的抑制窄带干扰的技术是变换域抗干扰技术，其主要思想是利用有效信号和窄带干扰信号在不同变换域中的差异极大的特性，利用这种特征将窄带干扰信号检测出，从而对含有干扰的子频谱空间进行频谱幅度的衰减或置零处理，再通过正交反变换恢复到时域，达到抑制窄带干扰的目的。
在变换域干扰技术中，应用最为广泛的是一种基于两路重叠加窗的抑制干扰的方法，即对两路数据进行处理，并对其中一路进行50％时钟延迟处理，窗函数采用布莱克曼窗函数或汉明窗函数。但该技术的主要缺点是应用的窗函数最大旁瓣较高，即使具有较低最大旁瓣的布莱克曼窗函数也只有-57dB，对干扰信号的抑制效果较差，已经无法满足越来越高的抗干扰要求，而受限于该方法的处理流程，如果采用较低最大旁瓣的窗函数又会带来较大的信噪比损失，从而导致卫星导航定位系统的性能急剧下降。
发明内容
为解决上述技术问题，本发明提供了一种抑制窄带干扰的方法及装置，以实现对窄带干扰较好的抑制效果，且不会带来较大信噪比损失的目的。
为实现上述技术目的，本发明实施例提供了如下技术方案：
一种抑制窄带干扰的方法，应用于卫星导航定位系统，所述方法包括：
对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；
分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；
将所述各路处理信号进行叠加处理，获得输出信号。
优选的，对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号为对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号；
对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号包括：
将所述卫星导航定位系统中的中频信号作为第一初始信号；
对所述中频信号进行25％时钟延迟处理获得第二初始信号；
对所述中频信号进行50％时钟延迟处理获得第三初始信号；
对所述中频信号进行75％时钟延迟处理获得第四初始信号。
优选的，分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号包括：
对所述第一初始信号进行干扰抑制处理获得第一处理信号；
对所述第二初始信号进行干扰抑制处理获得第二处理信号；
对所述第三初始信号进行干扰抑制处理获得第三处理信号；
对所述第四初始信号进行干扰抑制处理获得第四处理信号。
优选的，将所述各路处理信号进行叠加处理，获得输出信号包括：
将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；
将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；
将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；
将所述第四处理信号和所述第一信号、第二信号及第三信号叠加，获得所述输出信号。
优选的，所述干扰抑制处理包括：
取出多个采样点数据，并对所述多个采样点数据以预设窗函数进行加窗处理；
对加窗处理后的多个采样点数据进行傅里叶变换，将所述加窗处理后的多个采样点数据由时域变换为频域；
确定干扰门限，将在频域中高于干扰门限的采样点数据判定为干扰谱线，并对所述干扰谱线进行抑制处理；
对频域中的多个采样点数据进行傅里叶逆变换。
优选的，取出多个采样点数据包括：
取出2N个采样点数据，N为正整数。
优选的，N的取值范围为256-4096，包括端点值。
优选的，对所述干扰谱线进行抑制处理包括：
将所述干扰谱线的幅值置为零。
优选的，所述预设窗函数为布莱克曼哈里斯Blackman-Harris窗函数或凯撒Kaiser窗函数。
一种抑制窄带干扰的装置，应用于卫星导航定位系统，所述装置包括：
延迟处理单元，用于对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；
干扰抑制单元，用于分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；
叠加单元，用于将所述各路处理信号进行叠加处理，获得输出信号。
优选的，所述延迟处理单元包括第一延迟单元、第二延迟单元、第三延迟单元和第四延迟单元；所述干扰抑制单元包括第一抑制单元、第二抑制单元、第三抑制单元和第四抑制单元；其中，
所述第一延迟单元，用于将所述卫星导航定位系统中的中频信号作为第一初始信号；
所述第二延迟单元，用于对所述中频信号进行25％时钟延迟处理获得第二初始信号；
所述第三延迟单元，用于对所述中频信号进行50％时钟延迟处理获得第三初始信号；
所述第四延迟单元，用于对所述中频信号进行75％时钟延迟处理获得第四初始信号；
所述第一抑制单元，用于对所述第一初始信号进行干扰抑制处理获得第一处理信号；
所述第二抑制单元，用于对所述第二初始信号进行干扰抑制处理获得第二处理信号；
所述第三抑制单元，用于对所述第三初始信号进行干扰抑制处理获得第三处理信号；
所述第四抑制单元，用于对所述第四初始信号进行干扰抑制处理获得第四处理信号。
从上述技术方案可以看出，本发明实施例提供了一种抑制窄带干扰的方法及装置，其中，所述方法首先对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；然后对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号；最后将所述各路处理信号进行叠加处理，获得最终的输出信号。由于在所述干扰抑制处理的过程中采用了最大旁瓣低于布莱克曼窗函数最大旁瓣的预设窗函数，因此所述方法对窄带干扰的抑制效果较好；并且由于所述方法采用同时对至少四路不同百分比时钟延迟的初始信号进行干扰抑制处理后叠加，获得最终输出信号的措施避免了所述输出信号的信噪比损失较大而导致卫星导航定位系统性能急剧下降的情况出现。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本申请的一个实施例提供的一种抑制窄带干扰的方法的流程示意图；
图2为本申请的一个实施例提供的一种干扰抑制处理的流程示意图；
图3为本申请的一个优选实施例提供的一种抑制窄带干扰的方法的流程示意图；
图4和图5分别为本申请的一个实施例提供的加入强窄带干扰的GPS信号经过本申请实施例提供的抑制窄带干扰的方法处理前后的信号频域图；
图6为本申请的一个实施例提供的一种抑制窄带干扰的装置的结构示意图；
图7为本申请的一个优选实施例提供的一种抑制窄带干扰的装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本申请实施例提供了一种抑制窄带干扰的方法，应用于卫星导航定位系统，如图1所示，所述方法包括：
S101：对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号。
需要说明的是，在本申请的一个实施例中，对所述卫星导航定位系统中的中频信号进行延迟处理后，获得了五路不同百分比时钟延迟的初始信号，各路时钟延迟的百分比分别为0％、20％、40％、60％和80％。
在本申请的另一个实施例中，对所述中频信号进行延迟处理，获得八路不同百分比时钟延迟的初始信号，各路时钟延迟的百分比分别为0％、12.5％、25％、37.5％、50％、62.5％、75％和87.5％。本申请对此并不做限定，具体视实际情况而定。
还需要说明的是，对所述中频信号进行延迟处理，获得不同百分比时钟延迟的初始信号的路数为X时，各路时钟延迟的百分比分布根据X的不同而不同，具体分布方式已为本领域技术人员所熟知，本申请在此不做赘述。
但在本申请的一个优选实施例中，对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号包括为对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号；
对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号包括：
将所述卫星导航定位系统中的中频信号作为第一初始信号；
对所述中频信号进行25％时钟延迟处理获得第二初始信号；
对所述中频信号进行50％时钟延迟处理获得第三初始信号；
对所述中频信号进行75％时钟延迟处理获得第四初始信号。
S102：分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号,所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣。
在上述实施例的基础上，在本申请的一个实施例中，分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号包括：
对所述第一初始信号进行干扰抑制处理获得第一处理信号；
对所述第二初始信号进行干扰抑制处理获得第二处理信号；
对所述第三初始信号进行干扰抑制处理获得第三处理信号；
对所述第四初始信号进行干扰抑制处理获得第四处理信号。
在上述实施例的基础上，在本申请的另一个实施例中，将所述各路处理信号进行叠加处理，获得输出信号包括：
将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；
将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；
将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；
将所述第四处理信号和所述第一信号、第二信号及第三信号叠加，获得所述输出信号。
在上述实施例的基础上，在本申请的一个具体实施例中，如图2所示，所述干扰抑制处理包括：
A101：取出多个采样点数据，并对所述多个采样点数据以预设窗函数进行加窗处理；
A102：对加窗处理后的多个采样点数据进行傅里叶变换，将所述加窗处理后的多个采样点数据由时域变换为频域；
A103：确定干扰门限，将在频域中高于干扰门限的采样点数据判定为干扰谱线，并对所述干扰谱线进行抑制处理；
A104：对频域中的多个采样点数据进行傅里叶逆变换。
需要说明的是，在实际应用中，由于计算资源与存储空间有限，在傅里叶变换运算中，单位时间内只能对由2N个点组成的离散信号进行处理，N为正整数，因此必须对进行干扰抑制处理的信号进行2N个点的截断操作，即取出多个采样点数据即为取出2N个采样点数据。在本申请的一个实施例中，N的取值范围为256-4096，包括端点值。这种对于信号在时域的截断相当于对信号加了矩形窗，会造成频域的频谱泄漏现象。为了减小这种频谱泄漏的现象，需要对截断后的信号进行加窗处理，加窗处理使得信号的频谱更加集中在有限长度内，使得信号在截断处的频谱平滑过渡，减少了在截断处边缘的不连续性，而并不影响信号的波形。因此N的取值越高，经过加窗处理后的信号的频谱泄漏程度越小。但由于硬件限制，N的取值不可能无限的高，因此在本申请的一个优选实施例中，N的取值为1024。在对硬件要求较低的情况下实现所述方法的抑制窄带干扰的要求。
同时加窗处理也会对信号带来一定的信噪比损耗。选择不同的窗函数，带来的信噪比损耗也不同，各窗函数信噪比损耗大小是实际应用中选择窗函数的一个重要因素。选择窗函数的原则是选用主瓣尽量宽且最大旁瓣尽量小。
表1 常用窗函数的基本参数
常用的窗函数基本参数如表1所示。从表1可以看出，布莱克曼哈里斯窗函数的最大旁瓣最低，可以对更强的窄带干扰进行抑制。因此在上述实施例的基础上，在本申请的一个优选实施例中，所述预设窗函数为布莱克曼哈里斯窗函数。但在本申请的其他实施例中，所述预设窗函数还可以为凯撒窗函数。但是布莱克曼哈里斯窗函数为余弦窗，更便于硬件实现，而凯撒窗函数为非余弦窗，较难由硬件实现。因此在实际应用中，更广泛采用布莱克曼哈里斯窗函数。
还需要说明的是，采用最大旁瓣较低的窗函数进行加窗处理虽然可以对更强的窄带干扰进行抑制，但同时带来的信噪比损失也较大，本申请通过将信号同时分为四路以上的信号同时进行干扰抑制处理，并且将经过干扰抑制处理的信号进行叠加的方法补回信噪比的损失。
在上述实施例的基础上，在本申请的一个优选实施例中，对加窗处理后的多个采样点数据进行傅里叶变换的过程中采用蝶形算法实现快速运算。由于采用蝶形算法进行傅里叶变换的方法和具体流程已为本领域技术人员所熟知，本申请在此不做赘述。
在上述实施例的基础上，本申请的一个实施例提供了一种具体的确定干扰门限的方法：在频域中采用门限检测技术确定干扰门限。
需要说明的是，所述门限检测技术包括固定门限方法和自适应门限方法。其中，固定门限方法是人为设置一个固定值作为干扰门限。而自适应门限法而是根据采样点数据的谱线分布，将谱线的幅值的平方与一个特定系数的乘积作为干扰门限。在本申请的一个优选实施例中，通过自适应门限法确定干扰门限。通过自适应门限法确定干扰门限的方法适用范围更为广泛。
在上述实施例的基础上，在本申请的又一个实施例中，对所述干扰谱线进行抑制处理包括：
将所述干扰谱线的幅值置为零。
S103：将所述各路处理信号进行叠加处理，获得输出信号。
需要说明的是，本申请通过将所述各路处理信号进行叠加处理获得输出信号的方法弥补采用最大旁瓣低于布莱克曼窗函数的预设窗函数而带来的信噪比损失。
在上述实施例的基础上，在本申请的再一个实施例中，将所述各路处理信号进行叠加处理，获得输出信号包括：
将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；
将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；
将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；
将所述第四处理信号和所述第一处理信号、第二处理信号及第三处理信号叠加，获得所述输出信号。
在上述实施例的基础上，本申请的一个具体优选实施例提供了一种抑制窄带干扰的方法，如图3所示，包括：
S201：将所述卫星导航定位系统中的中频信号作为第一初始信号；对所述中频信号进行25％时钟延迟处理获得第二初始信号；对所述中频信号进行50％时钟延迟处理获得第三初始信号；对所述中频信号进行75％时钟延迟处理获得第四初始信号；
S202：对所述第一初始信号进行干扰抑制处理获得第一处理信号；对所述第二初始信号进行干扰抑制处理获得第二处理信号；对所述第三初始信号进行干扰抑制处理获得第三处理信号；对所述第四初始信号进行干扰抑制处理获得第四处理信号，所述干扰抑制处理中采用的窗函数为最大旁瓣低于布莱克曼窗函数的预设窗函数；
S203：将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；将所述第四处理信号和所述第一信号、第二信号及第三信号叠加，获得所述输出信号。
需要说明的是，经发明人研究发现，将所述中频信号分为四路进行干扰抑制处理后再叠加的优点在于，在利用最少的硬件资源的情况下即可实现对窄带干扰较好的抑制效果。而将所述中频信号分为五路或以上时，并不会提升对窄带干扰的抑制效果，反而会造成硬件资源的浪费。
还需要说明的是，应用本实施例提供的所述方法，对加入强窄带干扰的GPS信号进行了测试；测试用GPS信号的干信比为120dB，其信号中心频率为1572.42MHz，数字基带中频为4.092MHz，采样率为16.368MHz。GPS信号测试前的信号频域图如图4所示，经过本实施例提供的抑制窄带干扰方法处理之后的信号频域图如图5所示。附图4和附图5的横坐标为采样点，纵坐标为幅值。对比图4和图5可以看出，本实施例提供的所述方法具有较好的窄带干扰抑制效果。
相应的，本申请还提供了一种抑制窄带干扰的装置，应用于卫星导航定位系统，如图6所示，所述装置包括：
延迟处理单元A100，用于对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；
干扰抑制单元A200，用于分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；
叠加单元A300，用于将所述各路处理信号进行叠加处理，获得输出信号。
在上述实施例的基础上，在本申请的一个实施例中，如图7所示，所述延迟处理单元A100包括第一延迟单元、第二延迟单元、第三延迟单元和第四延迟单元；所述干扰抑制单元A200包括第一抑制单元、第二抑制单元、第三抑制单元和第四抑制单元；其中，
所述第一延迟单元A110，用于将所述卫星导航定位系统中的中频信号作为第一初始信号；
所述第二延迟单元A120，用于对所述中频信号进行25％时钟延迟处理获得第二初始信号；
所述第三延迟单元A130，用于对所述中频信号进行50％时钟延迟处理获得第三初始信号；
所述第四延迟单元A140，用于对所述中频信号进行75％时钟延迟处理获得第四初始信号；
所述第一抑制单元A210，用于对所述第一初始信号进行干扰抑制处理获得第一处理信号；
所述第二抑制单元A220，用于对所述第二初始信号进行干扰抑制处理获得第二处理信号；
所述第三抑制单元A230，用于对所述第三初始信号进行干扰抑制处理获得第三处理信号；
所述第四抑制单元A240，用于对所述第四初始信号进行干扰抑制处理获得第四处理信号。
在上述实施例的基础上，在本申请的另一个实施例中，所述叠加单元用于将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；将所述第四处理信号和所述第一信号、第二信号及第三信号叠加，获得所述输出信号。
在上述实施例的基础上，在本申请的一个具体实施例中，所述干扰抑制处理包括：
A101：取出多个采样点数据，并对所述多个采样点数据以预设窗函数进行加窗处理；
A102：对加窗处理后的多个采样点数据进行傅里叶变换，将所述加窗处理后的多个采样点数据由时域变换为频域；
A103：确定干扰门限，将在频域中高于干扰门限的采样点数据判定为干扰谱线，并对所述干扰谱线进行抑制处理；
A104：对频域中的多个采样点数据进行傅里叶逆变换。
需要说明的是，在实际应用中，由于计算资源与存储空间有限，在傅里叶变换运算中，单位时间内只能对由2N个点组成的离散信号进行处理，因此必须对进行干扰抑制处理的信号进行2N个点的截断操作，即取出多个采样点数据即为取出2N个采样点数据。在本申请的一个实施例中，N的取值范围为256-4096，包括端点值。在本申请的一个优选实施例中，N的取值为1024。
在上述实施例的基础上，在本申请的一个优选实施例中，对加窗处理后的多个采样点数据进行傅里叶变换的过程中采用蝶形算法实现快速运算。由于采用蝶形算法进行傅里叶变换的方法和具体流程已为本领域技术人员所熟知，本申请在此不做赘述。
在上述实施例的基础上，本申请的一个实施例提供了一种具体的确定干扰门限的方法：在频域中采用门限检测技术确定干扰门限。
需要说明的是，所述门限检测技术包括固定门限方法和自适应门限方法。其中，固定门限方法是人为设置一个固定值作为干扰门限。而自适应门限法而是根据采样点数据的谱线分布，将谱线的幅值的平方与一个特定系数的乘积作为干扰门限。在本申请的一个优选实施例中，通过自适应门限法确定干扰门限。通过自适应门限法确定干扰门限的方法适用范围更为广泛。
在上述实施例的基础上，在本申请的又一个实施例中，对所述干扰谱线进行抑制处理包括：
综上所述，本发明实施例提供了一种抑制窄带干扰的方法及装置，其中，所述方法首先对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；然后对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号；最后将所述各路处理信号进行叠加处理，获得最终的输出信号。由于在所述干扰抑制处理的过程中采用了最大旁瓣低于布莱克曼窗函数最大旁瓣的预设窗函数，因此所述方法对窄带干扰的抑制效果较好；并且由于所述方法采用同时对至少四路不同百分比时钟延迟的初始信号进行干扰抑制处理后叠加，获得最终输出信号的措施避免了所述输出信号的信噪比损失较大而导致卫星导航定位系统性能急剧下降的情况出现。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (11)
1.一种抑制窄带干扰的方法，应用于卫星导航定位系统，其特征在于，所述方法包括：
对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；
分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；
将所述各路处理信号进行叠加处理，获得输出信号。
2.根据权利要求1所述的方法，其特征在于，对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号为对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号；
对所述卫星导航定位系统中的中频信号进行延迟处理，获得四路不同百分比时钟延迟的初始信号包括：
将所述卫星导航定位系统中的中频信号作为第一初始信号；
对所述中频信号进行25％时钟延迟处理获得第二初始信号；
对所述中频信号进行50％时钟延迟处理获得第三初始信号；
对所述中频信号进行75％时钟延迟处理获得第四初始信号。
3.根据权利要求2所述的方法，其特征在于，分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号包括：
对所述第一初始信号进行干扰抑制处理获得第一处理信号；
对所述第二初始信号进行干扰抑制处理获得第二处理信号；
对所述第三初始信号进行干扰抑制处理获得第三处理信号；
对所述第四初始信号进行干扰抑制处理获得第四处理信号。
4.根据权利要求3所述的方法，其特征在于，将所述各路处理信号进行叠加处理，获得输出信号包括：
将所述第一处理信号进行75％的时钟延迟处理，获得第一信号；
将所述第二处理信号进行50％的时钟延迟处理，获得第二信号；
将所述第三处理信号进行25％的时钟延迟处理，获得第三信号；
将所述第四处理信号和所述第一信号、第二信号及第三信号叠加，获得所述输出信号。
5.根据权利要求1所述的方法，其特征在于，所述干扰抑制处理包括：
取出多个采样点数据，并对所述多个采样点数据以预设窗函数进行加窗处理；
对加窗处理后的多个采样点数据进行傅里叶变换，将所述加窗处理后的多个采样点数据由时域变换为频域；
确定干扰门限，将在频域中高于干扰门限的采样点数据判定为干扰谱线，并对所述干扰谱线进行抑制处理；
对频域中的多个采样点数据进行傅里叶逆变换。
6.根据权利要求5所述的方法，其特征在于，取出多个采样点数据包括：
取出2N个采样点数据，N为正整数。
7.根据权利要求6所述的方法，其特征在于，N的取值范围为256-4096，包括端点值。
8.根据权利要求5所述的方法，其特征在于，对所述干扰谱线进行抑制处理包括：
将所述干扰谱线的幅值置为零。
9.根据权利要求1-8任一项所述的方法，其特征在于，所述预设窗函数为布莱克曼哈里斯Blackman-Harris窗函数或凯撒Kaiser窗函数。
10.一种抑制窄带干扰的装置，应用于卫星导航定位系统，其特征在于，所述装置包括：
延迟处理单元，用于对所述卫星导航定位系统中的中频信号进行延迟处理，获得至少四路不同百分比时钟延迟的初始信号；
干扰抑制单元，用于分别对所述各路初始信号进行干扰抑制处理获得与各路初始信号一一对应的处理信号，所述干扰抑制处理中采用的预设窗函数的最大旁瓣低于布莱克曼窗函数的最大旁瓣；
叠加单元，用于将所述各路处理信号进行叠加处理，获得输出信号。
11.根据权利要求10所述的装置，其特征在于，所述延迟处理单元包括第一延迟单元、第二延迟单元、第三延迟单元和第四延迟单元；所述干扰抑制单元包括第一抑制单元、第二抑制单元、第三抑制单元和第四抑制单元；其中，
所述第一延迟单元，用于将所述卫星导航定位系统中的中频信号作为第一初始信号；
所述第二延迟单元，用于对所述中频信号进行25％时钟延迟处理获得第二初始信号；
所述第三延迟单元，用于对所述中频信号进行50％时钟延迟处理获得第三初始信号；
所述第四延迟单元，用于对所述中频信号进行75％时钟延迟处理获得第四初始信号；
所述第一抑制单元，用于对所述第一初始信号进行干扰抑制处理获得第一处理信号；
所述第二抑制单元，用于对所述第二初始信号进行干扰抑制处理获得第二处理信号；
所述第三抑制单元，用于对所述第三初始信号进行干扰抑制处理获得第三处理信号；
所述第四抑制单元，用于对所述第四初始信号进行干扰抑制处理获得第四处理信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610215072.7A CN105911565B (zh) | 2016-04-08 | 2016-04-08 | 一种抑制窄带干扰的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610215072.7A CN105911565B (zh) | 2016-04-08 | 2016-04-08 | 一种抑制窄带干扰的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105911565A true CN105911565A (zh) | 2016-08-31 |
| CN105911565B CN105911565B (zh) | 2018-12-21 |
Family
ID=56744803
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610215072.7A Active CN105911565B (zh) | 2016-04-08 | 2016-04-08 | 一种抑制窄带干扰的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105911565B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526620A (zh) * | 2016-11-04 | 2017-03-22 | 中国电子科技集团公司第五十四研究所 | 基于变换域滤波的窄带干扰抑制方法 |
| CN106850043A (zh) * | 2017-01-18 | 2017-06-13 | 成都科脉通信技术有限公司 | 基于超短波通信的窄带卫星通信系统 |
| CN107789008A (zh) * | 2017-11-29 | 2018-03-13 | 声泰特（成都）科技有限公司 | 一种基于通道数据的自适应超声波束合成方法和系统 |
| CN108551351A (zh) * | 2018-04-13 | 2018-09-18 | 中国科学院微电子研究所 | 抑制窄带干扰的方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020118133A1 (en) * | 2001-02-23 | 2002-08-29 | Desjardins Gerard A. | System and method for computing navigation information in the presence of interference |
| EP1862816B1 (en) * | 2006-06-04 | 2013-03-20 | SISVEL Technology Srl | Time-frequency interfering detector |
| CA2877790A1 (en) * | 2012-06-27 | 2014-01-03 | Airbus Defence And Space Limited | Authentication of satellite navigation signals |
| CN203894413U (zh) * | 2013-12-16 | 2014-10-22 | 重庆九洲星熠导航设备有限公司 | 一种基于fpga的卫星导航信号抗窄带干扰处理机 |
-
2016
- 2016-04-08 CN CN201610215072.7A patent/CN105911565B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020118133A1 (en) * | 2001-02-23 | 2002-08-29 | Desjardins Gerard A. | System and method for computing navigation information in the presence of interference |
| EP1862816B1 (en) * | 2006-06-04 | 2013-03-20 | SISVEL Technology Srl | Time-frequency interfering detector |
| CA2877790A1 (en) * | 2012-06-27 | 2014-01-03 | Airbus Defence And Space Limited | Authentication of satellite navigation signals |
| CN203894413U (zh) * | 2013-12-16 | 2014-10-22 | 重庆九洲星熠导航设备有限公司 | 一种基于fpga的卫星导航信号抗窄带干扰处理机 |
Non-Patent Citations (1)
| Title |
|---|
| 孟东等: "反加窗算法及其在扩频系统窄带干扰抑制中的应用", 《电子与信息学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526620A (zh) * | 2016-11-04 | 2017-03-22 | 中国电子科技集团公司第五十四研究所 | 基于变换域滤波的窄带干扰抑制方法 |
| CN106526620B (zh) * | 2016-11-04 | 2018-10-19 | 中国电子科技集团公司第五十四研究所 | 基于变换域滤波的窄带干扰抑制方法 |
| CN106850043A (zh) * | 2017-01-18 | 2017-06-13 | 成都科脉通信技术有限公司 | 基于超短波通信的窄带卫星通信系统 |
| CN107789008A (zh) * | 2017-11-29 | 2018-03-13 | 声泰特（成都）科技有限公司 | 一种基于通道数据的自适应超声波束合成方法和系统 |
| CN108551351A (zh) * | 2018-04-13 | 2018-09-18 | 中国科学院微电子研究所 | 抑制窄带干扰的方法及装置 |
| CN108551351B (zh) * | 2018-04-13 | 2020-08-04 | 中国科学院微电子研究所 | 抑制窄带干扰的方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105911565B (zh) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105911565A (zh) | 一种抑制窄带干扰的方法及装置 | |
| CN101388688B (zh) | 一种用于直接序列扩频通信系统的扫频干扰抑制方法 | |
| CN105634543B (zh) | 一种抗窄带干扰方法及系统 | |
| CN203894413U (zh) | 一种基于fpga的卫星导航信号抗窄带干扰处理机 | |
| CN104898132A (zh) | 一种阈值处理和空频自适应算法结合的导航抗干扰算法 | |
| CN107733541A (zh) | 频谱感知的方法、装置、设备以及计算机可读存储介质 | |
| CN109429249A (zh) | 一种基于mr定位数据的天线方位角优化方法及设备 | |
| CN110837003B (zh) | 一种基于三角窗的双窗全相位dft同步相量测量方法及系统 | |
| CN104931942B (zh) | 一种消除雷达多径干扰的方法 | |
| CN102263601A (zh) | 一种宽带多信号检测方法 | |
| CN107870316A (zh) | 一种基于时差计算的tdoa定位方法、装置及系统 | |
| CN106872934B (zh) | L型电磁矢量传感器阵列解相干esprit参数估计方法 | |
| CN109814064A (zh) | 一种基于三阵元l型直角阵干涉仪测向方法和装置 | |
| CN103605870B (zh) | 一种相控阵雷达任务综合优先级计算方法 | |
| CN103105614A (zh) | 基于惯导辅助的空时域联合抗干扰方法 | |
| CN104215939A (zh) | 一种融合广义对称结构信息的知识辅助空时自适应处理方法 | |
| CN108551351B (zh) | 抑制窄带干扰的方法及装置 | |
| CN102820904A (zh) | 一种用于广播网定位信号的互相关干扰减轻方法及装置 | |
| CN103701736B (zh) | Lte系统中的干扰抑制方法及装置 | |
| Singh et al. | Propagation modes of low-and very-low-latitude whistlers | |
| CN102708707B (zh) | 管制空域中单元窗口的自适应调整方法及装置 | |
| CN112817014B (zh) | 一种适用于时频域自适应的窄带抗干扰方法 | |
| CN102497340A (zh) | 基于迭加窗抑制算法的时频门处理方法 | |
| Klinngam et al. | The occurrence of equatorial spread-F at conjugate stations in Southeast Asia | |
| Wang et al. | DOA estimation of smart antenna signal based on MUSIC algorithm |
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 |