CN108957120A - 一种微弱信号提取方法及系统 - Google Patents

一种微弱信号提取方法及系统 Download PDF

Info

Publication number
CN108957120A
CN108957120A CN201810474461.0A CN201810474461A CN108957120A CN 108957120 A CN108957120 A CN 108957120A CN 201810474461 A CN201810474461 A CN 201810474461A CN 108957120 A CN108957120 A CN 108957120A
Authority
CN
China
Prior art keywords
sequence
signal
value
correlation
transform
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
Application number
CN201810474461.0A
Other languages
English (en)
Other versions
CN108957120B (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.)
North China Institute of Aerospace Engineering
Institute of Geophysical and Geochemical Exploration of CAGS
Original Assignee
North China Institute of Aerospace Engineering
Institute of Geophysical and Geochemical Exploration of CAGS
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 North China Institute of Aerospace Engineering, Institute of Geophysical and Geochemical Exploration of CAGS filed Critical North China Institute of Aerospace Engineering
Priority to CN201810474461.0A priority Critical patent/CN108957120B/zh
Publication of CN108957120A publication Critical patent/CN108957120A/zh
Application granted granted Critical
Publication of CN108957120B publication Critical patent/CN108957120B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis
    • G01R23/165—Spectrum analysis; Fourier analysis using filters
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种微弱信号提取方法及系统。本发明首先对被测的微弱信号进行FFT变换，通过获取FFT变换后的数据计算被测信号的各个频率分量的频率信息，进而进行互相关检测完成对微弱信号的恢复。本发明专利使用FFT变换和互相关结合的方法完成对微弱信号的提取，相比而言可以避免多次迭代，计算量相对自相关互相关结合的方法低。

Description

一种微弱信号提取方法及系统
技术领域
本发明涉及微弱信号提取领域，特别是涉及一种微弱信号提取方法及系统。
背景技术
在对微弱信号测量的过程中传统的互相关检测法需要知道接收信号的重复周期或者频率信息，于接收端产生同接收信号相同频率的参考信号进而完成对微弱信号的提取。但是对于多音信号或者非周期信号，即微弱信号中包含多个频率分量的情况，则需要取得被测信号的自相关函数逐次迭代进而确定被测信号中的频率信息。当被测信号中含有较多的频率分量时这一过程将十分复杂，且计算量很大。
发明内容
针对上述问题，本发明提供了一种微弱信号提取方法及系统。
为实现上述目的，本发明提供了如下方案：
一种微弱信号提取方法，所述方法包括：
按照设定的频率采集外部输入的微弱信号，得到采样序列，所述采样序列包括多个信号点以及多个频率分量；
对所述采样序列进行快速傅里叶变换，得到变换序列；
计算所述变换序列中各段幅度绝对值最高的信号点对应的频率；
根据所述变换序列生成标准正余弦序列；
计算所述标准正余弦序列与所述采样序列的相关值；
根据所述相关值，计算各频率分量的相位值以及幅度值；
通过逐次比较法获取所述各频率分量的幅度值的最大幅度值；
通过幅度阈值以及所述最大幅度值更新所有频率分量的幅度值；
根据所述幅度绝对值最高的信号点对应的频率、所述各频率分量的相位值以及更新后的频率分量的幅度值，对微弱信号进行提取。
可选的，所述计算所述变换序列中各段幅度绝对值最高的信号点对应的频率，具体包括：
将所述变换序列分为多段；
通过逐次比较法获取各段中幅度绝对值最高的信号点以及幅度绝对值最高的信号点的位置；
根据各段所述幅度绝对值最高的信号点的位置，计算各段内所述幅度绝对值最高的信号点对应的频率。
可选的，所述根据所述变换序列生成标准正余弦序列，具体包括：
根据各段内所述幅度绝对值最高的信号点对应的频率，生成标准正余弦序列。
可选的，所述根据所述相关值，计算各频率分量的相位值以及幅度值，具体包括：
根据所述相关值，采用四象限反正切法计算各频率分量的相位值；
根据所述相关值以及所述相位值计算各频率分量的幅度值。
一种微弱信号提取系统，所述系统包括：
前置放大采集电路，用于采集并放大外部输入的微弱信号，得到采样序列；
FPGA，与所述前置放大采集电路连接，用于对所述采样序列进行快速傅里叶变换，得到变换序列；
处理器，与所述FPGA连接，用于根据所述变换序列对所述微弱信号进行提取。
可选的，所述前置放大采集电路包括：
低通滤波器，用于对所述微弱信号进行滤波；
程控增益放大器，与所述的低通滤波器连接，用于对滤波后的信号进行放大；
单端差分变换电路，与所述程控增益放大器连接，用于对放大后的信号进行变换；
高精ADC，与所述单端差分变换电路连接，用于采集变换后的信号，得到采样序列；
基准源，分别与所述单端差分变换电路以及所述高精ADC连接，用于提供基准电压源。
可选的，还包括存储器，与所述处理器连接，用于存储提取后的微弱信号。
可选的，还包括缓存器，与所述处理器连接，用于存储所述处理器生成的标准正余弦序列。
与现有技术相比，本发明具有以下技术效果：本发明首先对被测的微弱信号进行FFT变换，通过获取FFT变换后的数据计算被测信号的各个频率分量的频率信息，进而进行互相关检测完成对微弱信号的恢复。本发明专利使用FFT变换和互相关结合的方法完成对微弱信号的提取，相比而言可以避免多次迭代，计算量相对自相关互相关结合的方法低。本发明专利的方法尤其适用于对多音信号(即包含大量频率分量)的提取。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例微弱信号提取方法的流程图；
图2为本发明实施例微弱信号提取系统的结构框图；
图3为本发明实施例前置放大采集电路的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例微弱信号提取方法的流程图。如图1所示，一种微弱信号提取方法包括：
步骤101：按照设定的频率采集外部输入的微弱信号，得到采样序列，所述采样序列包括多个信号点以及多个频率分量。设采集系统的采样率为fs，对外部输入信号进行采集并获得采样序列：x[n],n∈[1,N]。x[n]为输入信号x(t)的采样序列；
步骤102：对所述采样序列进行快速傅里叶变换，得到变换序列。对采样序列x[n]进行FFT变换得：x(k),k∈[1,N]。x(k)为采样序列x[n]的FFT变换。
步骤103：计算所述变换序列中各段幅度绝对值最高的信号点对应的频率。
具体的，将所述变换序列分为多段；通过逐次比较法获取各段中幅度绝对值最高的信号点以及幅度绝对值最高的信号点的位置；根据各段所述幅度绝对值最高的信号点的位置，计算各段内所述幅度绝对值最高的信号点对应的频率。
将x(k)的数据区间均匀分成N1段，用逐次比较法取得每段数据的最大值x(ki)及相应数据最大值的位置ki，
根据每段数据最大值的位置ki计算该段内主能量的对应频率值主能量指每段区间能量最大的那个点频信号，能量最大即幅度绝对值最高。
N代表采样点的个数，N1为将x[k]数据分成的区间数，fs为采样率，ki每段数据最大值的位置，该段内主能量的对应频率值。
步骤104：根据所述变换序列生成标准正余弦序列。根据各段内所述幅度绝对值最高的信号点对应的频率，生成标准正余弦序列。根据计算得到的主能量频率值于系统存储空间中生成长度为N的标准正余弦序列。
N代表采样点的个数，N1为将x[k]数据分成的区间数。fs为采样率，ki每段数据最大值的位置，该段内主能量的对应频率值。yi[n]生成的标准正弦序列，zi[n]生成的标准余弦序列。
步骤105：计算所述标准正余弦序列与所述采样序列的相关值。
Rxyi[0]为生产的标准正弦序列与采样序列的相关值；
Rxzi[0]为生产的标准余弦序列与采样序列的相关值。
步骤106：根据所述相关值，计算各频率分量的相位值以及幅度值。
根据所述相关值，采用四象限反正切法计算各频率分量的相位值；
为各频率分量的相位值。
根据所述相关值以及所述相位值计算各频率分量的幅度值；
各频率分量的幅度值
步骤107：通过逐次比较法获取所述各频率分量的幅度值的最大幅度值，
Ak_max所述各频率分量的幅度值的最大幅度值。
步骤108：通过幅度阈值以及所述最大幅度值更新所有频率分量的幅度值。设置阈值AMP并对每个频率分量的幅度依照下述规则进行更新，
步骤109：根据所述幅度绝对值最高的信号点对应的频率、所述各频率分量的相位值以及更新后的频率分量的幅度值，对微弱信号进行提取，
根据本发明提供的具体实施例，本发明具有以下技术效果：本发明首先对被测的微弱信号进行FFT变换，通过获取FFT变换后的数据计算被测信号的各个频率分量的频率信息，进而进行互相关检测完成对微弱信号的恢复。本发明专利使用FFT变换和互相关结合的方法完成对微弱信号的提取，相比而言可以避免多次迭代，计算量相对自相关互相关结合的方法低。本发明专利的方法尤其适用于对多音信号(即包含大量频率分量)的提取。
图2为本发明实施例微弱信号提取系统的结构框图；图3为本发明实施例前置放大采集电路的结构框图。
如图2所示，一种微弱信号提取系统，所述系统包括：
前置放大采集电路201，用于采集并放大外部输入的微弱信号，得到采样序列。前置放大采集电路201的实时方案如图3所示。输入被测信号首先经过RC/LC低通滤波器2011进行滤波，滤波后的信号送入程控增益放大器2012(AD8250)中进行放大，放大完成后的信号经过单端差分变换电路2013(AD8475)将单端信号转换为差分信号，转换后的差分信号送入高精ADC2014(24位高精ADC，ADS1271)，进行采集以获得采样序列x[n]。高精度基准源2015(REF5020)为整个前置放大采集电路提供基准电压源。前置放大采集电路201中所选用的器件仅仅是该前置放大采集电路的一种具体实施方案，根据不同的应用场景该前置放大采集电路的实施将不仅仅局限于上文描述的器件。
FPGA202，与所述前置放大采集电路201连接，用于对所述采样序列进行快速傅里叶变换，得到变换序列；
处理器203，与所述FPGA202连接，用于根据所述变换序列对所述微弱信号进行提取。
还包括存储器204，与所述处理器203连接，用于存储提取后的微弱信号。
缓存器205，与所述处理器203连接，用于存储所述处理器203生成的标准正余弦序列。
首先使用前置放大采集电路201对被测微弱信号进行放大并采集以获取采样序列x[n]。使用FPGA202完成步骤102中的FFT运算。步骤103到步骤109于图2中的ARM处理器203中实现，其中步骤104中生成的标准正余弦序列存入系统的SDRAM缓存器205中。最终提取后的微弱信号存入FLASH存储器204中。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种微弱信号提取方法，其特征在于，所述方法包括：
按照设定的频率采集外部输入的微弱信号，得到采样序列，所述采样序列包括多个信号点以及多个频率分量；
对所述采样序列进行快速傅里叶变换，得到变换序列；
计算所述变换序列中各段幅度绝对值最高的信号点对应的频率；
根据所述变换序列生成标准正余弦序列；
计算所述标准正余弦序列与所述采样序列的相关值；
根据所述相关值，计算各频率分量的相位值以及幅度值；
通过逐次比较法获取所述各频率分量的幅度值的最大幅度值；
通过幅度阈值以及所述最大幅度值更新所有频率分量的幅度值；
根据所述幅度绝对值最高的信号点对应的频率、所述各频率分量的相位值以及更新后的频率分量的幅度值，对微弱信号进行提取。
2.根据权利要求1所述的微弱信号提取方法，其特征在于，所述计算所述变换序列中各段幅度绝对值最高的信号点对应的频率，具体包括：
将所述变换序列分为多段；
通过逐次比较法获取各段中幅度绝对值最高的信号点以及幅度绝对值最高的信号点的位置；
根据各段所述幅度绝对值最高的信号点的位置，计算各段内所述幅度绝对值最高的信号点对应的频率。
3.根据权利要求2所述的微弱信号提取方法，其特征在于，所述根据所述变换序列生成标准正余弦序列，具体包括：
根据各段内所述幅度绝对值最高的信号点对应的频率，生成标准正余弦序列。
4.根据权利要求1所述的微弱信号提取方法，其特征在于，所述根据所述相关值，计算各频率分量的相位值以及幅度值，具体包括：
根据所述相关值，采用四象限反正切法计算各频率分量的相位值；
根据所述相关值以及所述相位值计算各频率分量的幅度值。
5.一种微弱信号提取系统，其特征在于，所述系统包括：
前置放大采集电路，用于采集并放大外部输入的微弱信号，得到采样序列；
FPGA，与所述前置放大采集电路连接，用于对所述采样序列进行快速傅里叶变换，得到变换序列；
处理器，与所述FPGA连接，用于根据所述变换序列对所述微弱信号进行提取。
6.根据权利要求5所述的系统，其特征在于，所述前置放大采集电路包括：
低通滤波器，用于对所述微弱信号进行滤波；
程控增益放大器，与所述的低通滤波器连接，用于对滤波后的信号进行放大；
单端差分变换电路，与所述程控增益放大器连接，用于对放大后的信号进行变换；
高精ADC，与所述单端差分变换电路连接，用于采集变换后的信号，得到采样序列；
基准源，分别与所述单端差分变换电路以及所述高精ADC连接，用于提供基准电压源。
7.根据权利要求5所述的系统，其特征在于，还包括存储器，与所述处理器连接，用于存储提取后的微弱信号。
8.根据权利要求5所述的系统，其特征在于，还包括缓存器，与所述处理器连接，用于存储所述处理器生成的标准正余弦序列。
CN201810474461.0A 2018-05-17 2018-05-17 一种微弱信号提取方法及系统 Active CN108957120B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810474461.0A CN108957120B (zh) 2018-05-17 2018-05-17 一种微弱信号提取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810474461.0A CN108957120B (zh) 2018-05-17 2018-05-17 一种微弱信号提取方法及系统

Publications (2)

Publication Number Publication Date
CN108957120A true CN108957120A (zh) 2018-12-07
CN108957120B CN108957120B (zh) 2020-08-11

Family

ID=64499259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810474461.0A Active CN108957120B (zh) 2018-05-17 2018-05-17 一种微弱信号提取方法及系统

Country Status (1)

Country Link
CN (1) CN108957120B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111122974A (zh) * 2019-12-31 2020-05-08 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种用于未知信号频率分析或已知信号频率校准的方法
CN111490789A (zh) * 2020-04-15 2020-08-04 中国科学技术大学 基于伪中位数累加的周期性弱信号检测方法及装置
CN115824394A (zh) * 2023-01-06 2023-03-21 山东省科学院海洋仪器仪表研究所 一种噪声测量中的微弱信号相位对齐方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101251390A (zh) * 2008-04-10 2008-08-27 中国科学院长春光学精密机械与物理研究所 一种基于时频变换的微弱信号检测装置
CN102213603A (zh) * 2011-04-13 2011-10-12 北京航空航天大学 一种基于单片机的低频微弱信号检测仪
CN103256044B (zh) * 2012-12-19 2016-01-06 电子科技大学 一种随钻声波信号处理装置
CN105553443A (zh) * 2015-12-10 2016-05-04 贵州省计量测试院 一种强噪声条件下微弱信号提取及数字处理系统
US20170115333A1 (en) * 2015-10-21 2017-04-27 Stojan Radic Sub-noise detection of a fast random event

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101251390A (zh) * 2008-04-10 2008-08-27 中国科学院长春光学精密机械与物理研究所 一种基于时频变换的微弱信号检测装置
CN102213603A (zh) * 2011-04-13 2011-10-12 北京航空航天大学 一种基于单片机的低频微弱信号检测仪
CN103256044B (zh) * 2012-12-19 2016-01-06 电子科技大学 一种随钻声波信号处理装置
US20170115333A1 (en) * 2015-10-21 2017-04-27 Stojan Radic Sub-noise detection of a fast random event
CN105553443A (zh) * 2015-12-10 2016-05-04 贵州省计量测试院 一种强噪声条件下微弱信号提取及数字处理系统

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111122974A (zh) * 2019-12-31 2020-05-08 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种用于未知信号频率分析或已知信号频率校准的方法
CN111122974B (zh) * 2019-12-31 2021-09-17 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种用于未知信号频率分析或已知信号频率校准的方法
CN111490789A (zh) * 2020-04-15 2020-08-04 中国科学技术大学 基于伪中位数累加的周期性弱信号检测方法及装置
CN111490789B (zh) * 2020-04-15 2023-03-10 中国科学技术大学 基于伪中位数累加的周期性弱信号检测方法及装置
CN115824394A (zh) * 2023-01-06 2023-03-21 山东省科学院海洋仪器仪表研究所 一种噪声测量中的微弱信号相位对齐方法
CN115824394B (zh) * 2023-01-06 2023-05-05 山东省科学院海洋仪器仪表研究所 一种噪声测量中的微弱信号相位对齐方法

Also Published As

Publication number Publication date
CN108957120B (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
CN108957120A (zh) 一种微弱信号提取方法及系统
CN103837740A (zh) 一种高精度数字瞬时测频方法及装置
CN203054093U (zh) 利用同步采样在混合信号中提取直流信号的装置
CN103185837B (zh) 一种电力系统频率测量的方法
CN106597215B (zh) 一种基于短时窗内暂态行波的故障测距方法及系统
CN110313903B (zh) 一种脉搏波频域特征参数提取方法及装置
CN108710027B (zh) 通道间相位差、幅度差的高精度测量方法
US20170296081A1 (en) Frame based spike detection module
CN103969508B (zh) 一种实时高精密的电力谐波分析方法及装置
CN109031261A (zh) 一种时差估计方法及装置
JP2012233824A (ja) パッシブレーダ装置、誘導装置及び電波検出方法
CN101968963A (zh) 音频信号压缩采样系统
CN103549947A (zh) 智能手机平台心电图波形实时准确绘制方法
CN113884996A (zh) 一种雷达专用测试设备多时宽脉冲信号校正方法和系统
CN112964931A (zh) 基于双通道欠采样的非理想多阻尼谐波信号参数测量方法
CN109633243B (zh) 一种基于多相位采样的束流信号峰值幅度精确提取方法
Huang et al. Sub-Nyquist sampling of multiple exponentially damped sinusoids with feedback structure
CN103605904A (zh) 基于误差估算的自补偿电力系统幅值算法
CN103513275A (zh) 具有采样率可调及自检功能的地震信号采集节点
CN108287328B (zh) 基于脉冲序列和波形重构的阶跃响应上升时间测试方法
CN106597049B (zh) 基于链表数组的多幅波形包络提取方法
KR20130110269A (ko) 샘플링 신호들 간의 차이를 이용하여 엔벨로프를 검출하는 방법 및 장치
CN111551785A (zh) 基于无迹卡尔曼滤波的频率与谐波检测方法
CN105322974B (zh) 一种干扰信号消除方法及用户终端
CN113050131A (zh) 一种基于预处理fft和栅栏效应校正的捕获方法

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