CN107689814A - 一种msk扩频接收机窄带干扰抑制方法 - Google Patents

一种msk扩频接收机窄带干扰抑制方法 Download PDF

Info

Publication number
CN107689814A
CN107689814A CN201610627464.4A CN201610627464A CN107689814A CN 107689814 A CN107689814 A CN 107689814A CN 201610627464 A CN201610627464 A CN 201610627464A CN 107689814 A CN107689814 A CN 107689814A
Authority
CN
China
Prior art keywords
data
module
frequency
time
processing
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.)
Pending
Application number
CN201610627464.4A
Other languages
English (en)
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.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
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 Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201610627464.4A priority Critical patent/CN107689814A/zh
Publication of CN107689814A publication Critical patent/CN107689814A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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/7101—Interference-related aspects the interference being narrowband interference with estimation filters
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种MSK扩频接收机窄带干扰抑制方法，通过A/D数据缓存模块、数据预处理模块、加窗处理模块、时频变换模块、频域滤波模块、时频逆变换模块实现。A/D数据缓存模块对扩频信号的A/D采集数据进行缓存处理；数据预处理模块对缓存数据进行重叠处理；加窗处理模块对数据进行加窗处理；时频变换模块对数据进行时频变换；频域滤波模块对频域数据进行滤波处理；时频逆变换模块对滤波处理后的频域数据进行时频逆变换。本发明解决传统窄带干扰抑制方法运算复杂、实现困难的问题。

Description

一种MSK扩频接收机窄带干扰抑制方法
技术领域
本发明涉及一种窄带干扰抑制方法，特别是一种MSK扩频接收机窄带干扰抑制方法。
背景技术
MSK扩频接收机需要对接收到的MSK直扩信号进行解扩处理，从而恢复出原始信息，当窄带干扰信号的功率远大于扩频信号的功率时，仅依靠扩频系统的处理增益来抑制窄带干扰难以达到很好的性能，因此需要采用其它方法来增强其抗干扰能力。传统的窄带干扰抑制方法多采用自适应滤波的方法，此方法需要预先知道干扰信号统计特性，然后通过调整滤波器的参数进行滤波来进行窄带干扰抑制，此方法原理简单，但运算复杂，实现困难。
发明内容
本发明的目的在于提供一种MSK扩频接收机窄带干扰抑制方法，解决传统窄带干扰抑制方法运算复杂、实现困难的问题。
一种MSK扩频接收机窄带干扰抑制方法，其具体步骤为：
第一步搭建窄带干扰抑制系统
窄带干扰抑制系统，包括：A/D采集数据缓存模块、数据预处理模块、加窗处理模块、时频变换模块、频域滤波模块和时频逆变换模块；所述A/D数据缓存模块的功能为：对MSK直扩信号的A/D采集数据进行缓存处理；数据预处理模块的功能为：对缓存数据进行重叠处理；加窗处理模块的功能为：对数据进行加窗处理；时频变换模块的功能为：对数据进行时频变换；频域滤波模块的功能为：对频域数据进行滤波处理；时频逆变换模块的功能为：对滤波处理后的频域数据进行时频逆变换。在FPGA芯片中运行上述模块。
第二步A/D数据缓存模块对MSK直扩信号的A/D采集数据进行缓存处理
A/D数据缓存模块产生采样时钟，时钟频率为fs，在其控制下A/D芯片对MSK直扩信号进行采集，采样点数为N，A/D数据缓存模块将这N个采样数据x(n),n＝0～N-1进行缓存。
时钟频率fs的选取原则：fs≥4f0，f0为MSK直扩信号的伪码速率。
第三步数据预处理模块对缓存数据进行重叠处理
数据预处理模块对A/D缓存数据分两路分别进行处理，一路将数据x(n),n＝0～N-1不 做处理，直接赋给数据s1(n),n＝0～N-1；另一路将数据x(n),n＝0～N-1进行50％重叠处理，把x(n)的0～N/2-1点作为新数据的后半部分，把x(n)的N/2～N-1点作为新数据的前半部分，新数据为s2(n),n＝0～N-1。
第四步加窗处理模块对数据进行加窗处理
加窗处理模块把数据s1(n)和s2(n)分别与汉明窗函数w(n)作时域乘积运算，得到数据ws1(n)和ws2(n)，其中ws1(n)为原始数据经加窗处理后得到的数据，ws2(n)为重叠数据经加窗处理后得到的数据。
当0≤n≤N-1时，汉明窗函数w(n)＝0.54-0.46cos[2π·n/(N-1)]，其它情况w(n)＝0。
第五步时频变换模块对数据进行时频变换
时频变换模块对数据ws1(n)和ws2(n)分别进行快速傅立叶变换：
把数据ws1(n)和ws2(n)从时域变换为频域数据WS1(k)和WS2(k)，k＝0～N-1，k为频域点数索引值。
第六步频域滤波模块对数据进行滤波处理
频率滤波模块对频域数据WS1(k)和WS2(k)，k＝0～N-1分别进行滤波处理，采用逐点比较方式，如果WS1(k)＞P，令WS1(k)＝0；其它情况WS1(k)保持原值；如果WS2(k)＞P，令WS2(k)＝0，其它情况WS2(k)保持原值。
P为判决门限，是理论计算的固定值，即：在无窄带干扰时，P为扩频信号频域数据中的最大值。
第七步时频逆变换模块对滤波处理的数据进行时频逆变换
时频逆变换模块对滤波处理后的数据进行时频逆变换，把WS1(k)和WS2(k)从频域变换为时域数据is1(n)和is2(n)，n＝0～N-1，n为时域点数索引值；之后对数据is2(n)，n＝0～N-1进行50％叠加处理后，与数据is1(n)相加，得到窄带干扰抑制后的时域数据。
至此，实现了MSK扩频接收机窄带干扰抑制。
本方法解决了传统抑制窄带干扰方法时域运算复杂、实现困难的问题，仿真实验从理论上证明了本方法的有效性和可行性，且本方法已在MSK扩频接收机演示验证板上进行了试验，试验结果证明：本方法能够对窄带干扰进行有效抑制，抑制效果良好，抑制后的数据满足系统要求。
具体实施方式
一种MSK扩频接收机窄带干扰抑制方法，其具体步骤为：
第一步搭建窄带干扰抑制系统
窄带干扰抑制系统，包括：A/D采集数据缓存模块、数据预处理模块、加窗处理模块、时频变换模块、频域滤波模块和时频逆变换模块；所述A/D数据缓存模块的功能为：对MSK直扩信号的A/D采集数据进行缓存处理；数据预处理模块的功能为：对缓存数据进行重叠处理；加窗处理模块的功能为：对数据进行加窗处理；时频变换模块的功能为：对数据进行时频变换；频域滤波模块的功能为：对频域数据进行滤波处理；时频逆变换模块的功能为：对滤波处理后的频域数据进行时频逆变换。在FPGA芯片中运行上述模块。
第二步A/D数据缓存模块对MSK直扩信号的A/D采集数据进行缓存处理
A/D数据缓存模块产生采样时钟，时钟频率为fs，在其控制下A/D芯片对MSK直扩信号进行采集，采样点数为N，A/D数据缓存模块将这N个采样数据x(n),n＝0～N-1进行缓存。
时钟频率fs的选取原则：fs≥4f0，f0为MSK直扩信号的伪码速率。
第三步数据预处理模块对缓存数据进行重叠处理
数据预处理模块对A/D缓存数据分两路分别进行处理，一路将数据x(n),n＝0～N-1不做处理，直接赋给数据s1(n),n＝0～N-1；另一路将数据x(n),n＝0～N-1进行50％重叠处理，把x(n)的0～N/2-1点作为新数据的后半部分，把x(n)的N/2～N-1点作为新数据的前半部分，新数据为s2(n),n＝0～N-1。
第四步加窗处理模块对数据进行加窗处理
加窗处理模块把数据s1(n)和s2(n)分别与汉明窗函数w(n)作时域乘积运算，得到数据ws1(n)和ws2(n)，其中ws1(n)为原始数据经加窗处理后得到的数据，ws2(n)为重叠数据经加窗处理后得到的数据。
当0≤n≤N-1时，汉明窗函数w(n)＝0.54-0.46cos[2π·n/(N-1)]，其它情况w(n)＝0。
第五步时频变换模块对数据进行时频变换
时频变换模块对数据ws1(n)和ws2(n)分别进行快速傅立叶变换：
把数据ws1(n)和ws2(n)从时域变换为频域数据WS1(k)和WS2(k)，k＝0～N-1，k为频域点数索引值。
第六步频域滤波模块对数据进行滤波处理
频率滤波模块对频域数据WS1(k)和WS2(k)，k＝0～N-1分别进行滤波处理，采用逐点比较方式，如果WS1(k)＞P，令WS1(k)＝0；其它情况WS1(k)保持原值；如果WS2(k)＞P，令WS2(k)＝0，其它情况WS2(k)保持原值。
P为判决门限，是理论计算的固定值，即：在无窄带干扰时，P为扩频信号频域数据中的最大值。
第七步时频逆变换模块对滤波处理的数据进行时频逆变换
时频逆变换模块对滤波处理后的数据进行时频逆变换，把WS1(k)和WS2(k)从频域变换为时域数据is1(n)和is2(n)，n＝0～N-1，n为时域点数索引值；之后对数据is2(n)，n＝0～N-1进行50％叠加处理后，与数据is1(n)相加，得到窄带干扰抑制后的时域数据。
至此，实现了MSK扩频接收机窄带干扰抑制。

Claims (1)

1.一种MSK扩频接收机窄带干扰抑制方法，其特征在于具体步骤为：
第一步搭建窄带干扰抑制系统
窄带干扰抑制系统，包括：A/D采集数据缓存模块、数据预处理模块、加窗处理模块、时频变换模块、频域滤波模块和时频逆变换模块；所述A/D数据缓存模块的功能为：对MSK直扩信号的A/D采集数据进行缓存处理；数据预处理模块的功能为：对缓存数据进行重叠处理；加窗处理模块的功能为：对数据进行加窗处理；时频变换模块的功能为：对数据进行时频变换；频域滤波模块的功能为：对频域数据进行滤波处理；时频逆变换模块的功能为：对滤波处理后的频域数据进行时频逆变换；在FPGA芯片中运行上述模块；
第二步A/D数据缓存模块对MSK直扩信号的A/D采集数据进行缓存处理
A/D数据缓存模块产生采样时钟，时钟频率为fs，在其控制下A/D芯片对MSK直扩信号进行采集，采样点数为N，A/D数据缓存模块将这N个采样数据x(n),n＝0～N-1进行缓存；
时钟频率fs的选取原则：fs≥4f0，f0为MSK直扩信号的伪码速率；
第三步数据预处理模块对缓存数据进行重叠处理
数据预处理模块对A/D缓存数据分两路分别进行处理，一路将数据x(n),n＝0～N-1不做处理，直接赋给数据s1(n),n＝0～N-1；另一路将数据x(n),n＝0～N-1进行50％重叠处理，把x(n)的0～N/2-1点作为新数据的后半部分，把x(n)的N/2～N-1点作为新数据的前半部分，新数据为s2(n),n＝0～N-1；
第四步加窗处理模块对数据进行加窗处理
加窗处理模块把数据s1(n)和s2(n)分别与汉明窗函数w(n)作时域乘积运算，得到数据ws1(n)和ws2(n)，其中ws1(n)为原始数据经加窗处理后得到的数据，ws2(n)为重叠数据经加窗处理后得到的数据；
当0≤n≤N-1时，汉明窗函数w(n)＝0.54-0.46cos[2π·n/(N-1)]，其它情况w(n)＝0；
第五步时频变换模块对数据进行时频变换
时频变换模块对数据ws1(n)和ws2(n)分别进行快速傅立叶变换：
把数据ws1(n)和ws2(n)从时域变换为频域数据WS1(k)和WS2(k)，k＝0～N-1，k为频域点数索引值；
第六步频域滤波模块对数据进行滤波处理
频率滤波模块对频域数据WS1(k)和WS2(k)，k＝0～N-1分别进行滤波处理，采用逐点比较方式，如果WS1(k)＞P，令WS1(k)＝0；其它情况WS1(k)保持原值；如果WS2(k)＞P，令WS2(k)＝0，其它情况WS2(k)保持原值；
P为判决门限，是理论计算的固定值，即：在无窄带干扰时，P为扩频信号频域数据中的最大值；
第七步时频逆变换模块对滤波处理的数据进行时频逆变换
时频逆变换模块对滤波处理后的数据进行时频逆变换，把WS1(k)和WS2(k)从频域变换为时域数据is1(n)和is2(n)，n＝0～N-1，n为时域点数索引值；之后对数据is2(n)，n＝0～N-1进行50％叠加处理后，与数据is1(n)相加，得到窄带干扰抑制后的时域数据；
至此，实现了MSK扩频接收机窄带干扰抑制。
CN201610627464.4A 2016-08-03 2016-08-03 一种msk扩频接收机窄带干扰抑制方法 Pending CN107689814A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610627464.4A CN107689814A (zh) 2016-08-03 2016-08-03 一种msk扩频接收机窄带干扰抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610627464.4A CN107689814A (zh) 2016-08-03 2016-08-03 一种msk扩频接收机窄带干扰抑制方法

Publications (1)

Publication Number Publication Date
CN107689814A true CN107689814A (zh) 2018-02-13

Family

ID=61150836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610627464.4A Pending CN107689814A (zh) 2016-08-03 2016-08-03 一种msk扩频接收机窄带干扰抑制方法

Country Status (1)

Country Link
CN (1) CN107689814A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108540168A (zh) * 2018-03-09 2018-09-14 西安电子科技大学 一种抗窄带干扰长码扩频同步头实时捕获的装置及方法
CN109100687A (zh) * 2018-10-29 2018-12-28 北京遥感设备研究所 一种雷达设备lfm脉冲信号pslr参数确定方法
CN110380756A (zh) * 2019-07-09 2019-10-25 北京遥感设备研究所 一种高动态扩频通信设备pn码序列被动相关检测平台及方法
CN111641421A (zh) * 2020-05-08 2020-09-08 北京遥感设备研究所 一种高动态数字接收机窄带干扰抑制系统及方法
CN115913279A (zh) * 2022-11-22 2023-04-04 中国电子科技集团公司第十研究所 一种复杂干扰环境的多载波扩频信号捕获方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0857379A1 (en) * 1995-10-24 1998-08-12 Motorola, Inc. Method and apparatus for interference suppression in spread spectrum signals
CN101841349A (zh) * 2010-04-12 2010-09-22 北京理工大学 一种用于直扩系统的mpsk窄带干扰抑制方法
CN102082748A (zh) * 2010-12-29 2011-06-01 成都国星通信有限公司 扩频接收机前端频域滤波抗干扰装置及其实现方法
CN103763229A (zh) * 2014-01-15 2014-04-30 西安理工大学 基于频域滤波的cdma通信抗干扰系统及方法
CN104007453A (zh) * 2014-05-07 2014-08-27 西安理工大学 概率搜索辅助的频域空域联合抗干扰方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0857379A1 (en) * 1995-10-24 1998-08-12 Motorola, Inc. Method and apparatus for interference suppression in spread spectrum signals
CN101841349A (zh) * 2010-04-12 2010-09-22 北京理工大学 一种用于直扩系统的mpsk窄带干扰抑制方法
CN102082748A (zh) * 2010-12-29 2011-06-01 成都国星通信有限公司 扩频接收机前端频域滤波抗干扰装置及其实现方法
CN103763229A (zh) * 2014-01-15 2014-04-30 西安理工大学 基于频域滤波的cdma通信抗干扰系统及方法
CN104007453A (zh) * 2014-05-07 2014-08-27 西安理工大学 概率搜索辅助的频域空域联合抗干扰方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108540168A (zh) * 2018-03-09 2018-09-14 西安电子科技大学 一种抗窄带干扰长码扩频同步头实时捕获的装置及方法
CN109100687A (zh) * 2018-10-29 2018-12-28 北京遥感设备研究所 一种雷达设备lfm脉冲信号pslr参数确定方法
CN110380756A (zh) * 2019-07-09 2019-10-25 北京遥感设备研究所 一种高动态扩频通信设备pn码序列被动相关检测平台及方法
CN110380756B (zh) * 2019-07-09 2021-02-26 北京遥感设备研究所 一种高动态扩频通信设备pn码序列被动相关检测平台及方法
CN111641421A (zh) * 2020-05-08 2020-09-08 北京遥感设备研究所 一种高动态数字接收机窄带干扰抑制系统及方法
CN115913279A (zh) * 2022-11-22 2023-04-04 中国电子科技集团公司第十研究所 一种复杂干扰环境的多载波扩频信号捕获方法及系统

Similar Documents

Publication Publication Date Title
CN107689814A (zh) 一种msk扩频接收机窄带干扰抑制方法
CN106814396A (zh) 一种基于vmd的矿山微震信号的降噪滤波方法
WO2020107269A1 (zh) 自适应语音增强方法和电子设备
CN104007453B (zh) 概率搜索辅助的频域空域联合抗干扰方法
Sun et al. Cross-correlation analysis and time delay estimation of a homologous micro-seismic signal based on the Hilbert–Huang transform
CN102664657A (zh) 一种自适应中值门限频域抗干扰装置及方法
CN109241823B (zh) 一种基于变分模态分解和支持向量回归的信号预测方法
CN103427870B (zh) 基于稀疏傅里叶变换的深度扩频捕获与干扰抑制方法
CN111881858A (zh) 一种微震信号多尺度去噪方法、装置及可读存储介质
CN102982808A (zh) 基于小波变换的语音去噪装置及方法
CN111680737A (zh) 差异性信噪比条件下的雷达辐射源个体识别方法
CN101527698A (zh) 基于希尔伯特黄变换和自适应陷波的非平稳干扰抑制方法
CN105573104A (zh) 基于改进emd的手表检测降噪方法
CN105871413A (zh) 低信噪比直接序列扩频信号检测方法
CN108802175A (zh) 一种基于emd的脉冲涡流检测信号降噪方法
CN103915102B (zh) 一种lfm水声多途信号的噪声抑制方法
CN103438983B (zh) 一种信号随机平均谱的数据处理方法
CN111641421A (zh) 一种高动态数字接收机窄带干扰抑制系统及方法
Zhao et al. Pipeline leak fault feature extraction based on wavelet packet analysis and application
CN114242103A (zh) 语音降噪方法、装置、计算机设备和存储介质
Vishwakarma et al. De-noising of Audio Signal using Heavy Tailed Distribution and comparison of wavelets and thresholding techniques
CN108205127B (zh) 一种基于稀疏表示的水声信号处理方法
Ma et al. Research of the optimal wavelet selection on entropy function
Jiang et al. Adaptive Speech Enhancement Algorithm Based on Hilbert-Huang Transform.
CN104361892A (zh) 一种干扰信号与语音调制信号混叠的窄带干扰识别方法

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180213