CN113705335A - 一种低频稀疏线谱信号时频超分辨率提取方法 - Google Patents
一种低频稀疏线谱信号时频超分辨率提取方法 Download PDFInfo
- Publication number
- CN113705335A CN113705335A CN202110796227.1A CN202110796227A CN113705335A CN 113705335 A CN113705335 A CN 113705335A CN 202110796227 A CN202110796227 A CN 202110796227A CN 113705335 A CN113705335 A CN 113705335A
- Authority
- CN
- China
- Prior art keywords
- line spectrum
- spectrum signal
- frequency
- iteration
- estimated
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 140
- 238000000605 extraction Methods 0.000 title claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 21
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 4
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种低频稀疏线谱信号时频超分辨率提取方法,能够实现时频超分辨率的低频稀疏线谱提取。该方法包括如下步骤:获得测量数据,初始化待估计参数,包括每一维线谱信号分量的幅值、频率和相位。执行期望最大算法EM流程;在每次迭代的E步骤,注入测量数据,结合上一次迭代的待估计参数的估计值,对线谱信号各分量进行估计。在每次迭代的M步骤,首先通过预先设置的频率分辨率,对线谱信号各分量的频率进行枚举,枚举的线谱信号分量的频率代入线谱信号的稀疏线谱先验模型,得到线谱信号各分量的幅值和相位,然后利用线谱频率规律知识优化更新线谱信号各分量的频率。迭代结束后获得待估计参数的取值作为线谱信号的参数估计结果。
Description
技术领域
本发明涉及信号处理技术领域,具体涉及一种低频稀疏线谱信号时频超分辨率提取方法。
背景技术
低频稀疏线谱信号的测量在故障检测、海洋声学探测、海洋声学通信等应用中有着广阔的应用前景与研究价值。然而复杂时变的测量环境往往会给测量带来非平稳、低信噪比等挑战,这给低频稀疏线谱信号的测量带来了极端的困难,如何在极短时间内实现低频稀疏线谱信号的高精度测量是急需解决的问题。
通常对线谱信号的检测有以下几类方法
1、基于基函数展开的方法
傅里叶变换、分数阶傅里叶变换、小波变换等方法,直接将测量数据以相应基函数展开到不同的基函数空间进行处理。然而该方法需要在时间分辨率与频率分辨率上进行权衡,无法同时保证高精度的时间分辨率和频率分辨率。同时,由于窗函数的影响,会给测量带来旁瓣谱的干扰,进而影响测量结果。
2、基于相关性的方法
随机共振、自相关、匹配滤波等方法,基于构造数据的模型信息或通过短时延迟相关。往往会受到较大的环境影响,且可能会造成测量信号的相位模糊。
3、基于压缩感知的相关方法
压缩感知方法主要利用了线谱的稀疏特性,通过构建字典矩阵对信号进行重构,核心是通过远低于奈奎斯特采样定律的采样频率对信号进行恢复,降低采样数据存储。但是该方法严重依赖信号的高信噪比,在信噪比较低的情况下会有不可用的情况。
4、先验与后验相结合的方法
卡尔曼滤波线性最小均方误差估计方法利用了模型的先验信息与测量数据的后验信息。但传统用法是用在了机动目标跟踪上,该应用场景下往往无法获取大量的测量数据。
因此目前缺少一种能够实现时频超分辨率的低频稀疏线谱提取方法。
发明内容
有鉴于此,本发明提供了一种低频稀疏线谱信号时频超分辨率提取方法,能够基于低频线谱的稀疏特性和频率规律约束,将短时间窗内极高采样率获得的测量数据同先验知识进行结合,进而实现时频超分辨率的低频稀疏线谱提取。
为达到上述目的,本发明的技术方案包括如下步骤:
S1、预先设置线谱信号的频带范围为(fL,fH),fL,fH分别为最低频率和最高频率;设定数据采样频率fs>>5fH对线谱信号进行数据采样,并设定时间窗为Tw,则采样获得的数据长度为N=Tw/fs,获取的测量数据为y(n),n为第n个采样点n=0,1,...N-1。
S2、根据测量数据y(n)长度N与预设线谱信号的维数M;初始化待估计参数,包括每一维线谱信号分量的幅值、频率和相位。
S3、执行期望最大算法EM流程,EM流程包括期望E步骤和最大化M步骤,设定迭代结束条件。
在每次迭代的E步骤,注入测量数据y(n),结合上一次迭代的待估计参数的估计值,对线谱信号各分量进行估计。
在每次迭代的M步骤,首先通过预先设置的频率分辨率,对线谱信号各分量的频率进行枚举,枚举的线谱信号分量的频率代入线谱信号的稀疏线谱先验模型,得到线谱信号各分量的幅值和相位,然后利用线谱频率规律知识优化更新线谱信号各分量的频率。
迭代结束后获得待估计参数的取值作为线谱信号的参数估计结果。
进一步地,初始化待估计参数,具体为:
进一步地,在每次迭代的E步骤,注入测量数据y(n),结合上一次迭代的待估计参数的估计值,对线谱信号各分量进行估计,具体为:
根据测量数据y(n)与上一轮迭代的M步骤得到的即第k次迭代的得到的第m维线谱信号分量的待估计参数其中包含第k次迭代得到的第m维线谱信号分量的幅值估计值第k次迭代得到的第m维线谱信号分量的频率估计值第k次迭代得到的第m维线谱信号分量的相位估计值
进一步地,M步骤,具体包括如下流程:
其中χ为不同维度线谱信号的频率知识。
优化第k+1次迭代得到的待估计参数的估计值其中包含第k+1次迭代得到的第m维线谱信号分量的幅值估计值第k+1次迭代得到的第m维线谱信号分量的频率估计值第k+1次迭代得到的第m维线谱信号分量的相位估计值
进一步地,迭代结束条件,具体为:
有益效果:
1、本发明提供的一种低频稀疏线谱信号时频超分辨率提取方法,主要基于贝叶斯统计学习思想,将测量数据(约束1)、线谱稀疏先验模型(约束2)和频率规律先验模型(约束3)作为约束条件,利用期望最大(Expectation Maximum,EM)算法,估计稀疏线谱信号每一维信号分量的参数。由于时间窗大小与频率搜索分辨率可控,所以本算法可实现时频超分辨率。本算法额外需要线谱信号稀疏与线谱频率规律的先验知识,同时需要在短时间窗内以极高采样率增加数据量,以便算法能够从数据中更好的学习到更多知识。
2、本发明提供的一种低频稀疏线谱信号时频超分辨率提取方法,针对低频线谱信号的恢复场景,提出了使用短时间窗提高时间分辨率,同时使用极高的采样率,通过采集更多的数据降低信号测量的不确定度。
3、本发明提供的一种低频稀疏线谱信号时频超分辨率提取方法,结合低频线谱的稀疏性先验,使用了最大似然估计(MLE)方法确定了各线谱分量中频率、幅度、相位存在的关系,通过控制频率的搜索间隔提高频率分辨率。
4、本发明提供的一种低频稀疏线谱信号时频超分辨率提取方法,结合低频线谱可能存在频率规律约束,提出了引入频率规律的先验模型,以提高算法的准确度与收敛速度。
附图说明
图1为本发明提供的一种低频稀疏线谱信号时频超分辨率提取方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种低频稀疏线谱信号时频超分辨率提取方法,其流程如图1所示,包括如下步骤:
S1、数据读取流程:
预先设置线谱信号的频带范围为(fL,fH),fL,fH分别为最低频率和最高频率;设定数据采样频率fs>>5fH对线谱信号进行数据采样,并设定时间窗为Tw,则采样获得的数据长度为N=Tw/fs,获取的测量数据为y(n),n为第n个采样点n=0,1,...N-1。
S2、初始化流程:
根据测量数据y(n)长度N与预设线谱信号的维数M;初始化待估计参数,包括每一维线谱信号分量的幅值、频率和相位;初始化待估计参数与第m维线谱信号的初始测量值 为第m维线谱信号的初始化幅值;为第m维线谱信号的初始化频率;为第m维线谱信号的初始化相位。
S3、执行期望最大算法EM流程:
在期望最大算法EM流程,期望最大算法(Expectation-Maximization,EM)是一种以迭代方式求解最大似然估计的算法,EM算法有两个步骤,期望步骤(E-step)和最大步骤(M-step),当完成k轮迭代后得到估计值。
EM流程包括期望E步骤和最大化M步骤,设定迭代结束条件;
在每次迭代的E步骤,注入测量数据y(n),结合上一次迭代的待估计参数的估计值,对线谱信号各分量进行估计。E步骤具体为:
根据测量数据y(n)与上一轮迭代的M步骤得到的即第k次迭代的得到的第m维线谱信号分量的待估计参数其中包含第k次迭代得到的第m维线谱信号分量的幅值估计值第k次迭代得到的第m维线谱信号分量的频率估计值第k次迭代得到的第m维线谱信号分量的相位估计值
在每次迭代的M步骤,首先通过预先设置的频率分辨率,对线谱信号各分量的频率进行枚举,枚举的线谱信号分量的频率代入线谱信号的稀疏线谱先验模型,得到线谱信号各分量的幅值和相位,然后利用线谱频率规律知识优化更新线谱信号各分量的频率。
M步骤,具体包括如下流程:
其中χ为不同维度线谱信号的频率知识,例如船螺旋桨噪声频谱存在调制谐波关系、钢琴存在着基频泛频关系、仪器设备异常故障时存在频率规律的知识。
优化第k+1次迭代得到的待估计参数的估计值其中包含第k+1次迭代得到的第m维线谱信号分量的幅值估计值第k+1次迭代得到的第m维线谱信号分量的频率估计值第k+1次迭代得到的第m维线谱信号分量的相位估计值
输出估计结果流程:迭代结束后获得待估计参数的取值作为线谱信号的参数估计结果。
迭代结束条件,具体为:
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种低频稀疏线谱信号时频超分辨率提取方法,其特征在于,包括如下步骤:
S1、预先设置线谱信号的频带范围为(fL,fH),fL,fH分别为最低频率和最高频率;设定数据采样频率fs>>5fH对所述线谱信号进行数据采样,并设定时间窗为Tw,则采样获得的数据长度为N=Tw/fs,获取的测量数据为y(n),n为第n个采样点n=0,1,...N-1;
S2、根据所述测量数据y(n)长度N与预设线谱信号的维数M;初始化待估计参数,包括每一维线谱信号分量的幅值、频率和相位;
S3、执行期望最大算法EM流程,所述EM流程包括期望E步骤和最大化M步骤,设定迭代结束条件;
在每次迭代的E步骤,注入所述测量数据y(n),结合上一次迭代的信号参数估计值,对线谱信号各分量进行估计;
在每次迭代的M步骤,首先通过预先设置的频率分辨率,对线谱信号各分量的频率进行枚举,枚举的线谱信号分量的频率代入线谱信号的稀疏线谱先验模型,得到线谱信号各分量的幅值和相位,然后利用线谱频率规律知识优化更新线谱信号各分量的频率;
迭代结束后获得待估计参数的取值作为线谱信号的参数估计结果。
4.如权利要求3所述的方法,其特征在于,所述M步骤,具体包括如下流程:
其中χ为不同维度线谱信号的频率知识;
优化第k+1次迭代得到的待估计参数的估计值其中包含第k+1次迭代得到的第m维线谱信号分量的幅值估计值第k+1次迭代得到的第m维线谱信号分量的频率估计值第k+1次迭代得到的第m维线谱信号分量的相位估计值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110796227.1A CN113705335B (zh) | 2021-07-14 | 2021-07-14 | 一种低频稀疏线谱信号时频超分辨率提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110796227.1A CN113705335B (zh) | 2021-07-14 | 2021-07-14 | 一种低频稀疏线谱信号时频超分辨率提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113705335A true CN113705335A (zh) | 2021-11-26 |
CN113705335B CN113705335B (zh) | 2022-11-08 |
Family
ID=78648595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110796227.1A Active CN113705335B (zh) | 2021-07-14 | 2021-07-14 | 一种低频稀疏线谱信号时频超分辨率提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113705335B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115730198A (zh) * | 2022-11-10 | 2023-03-03 | 浙江衡玖医疗器械有限责任公司 | 一种超声信号tof自动提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919240A (zh) * | 2018-04-23 | 2018-11-30 | 东南大学 | 一种基于组稀疏结构的水声目标辐射噪声调制谱重构方法 |
CN109946518A (zh) * | 2019-03-27 | 2019-06-28 | 河南天通电力有限公司 | 基于贝叶斯方法的电力谐波信号分析方法与分析设备 |
-
2021
- 2021-07-14 CN CN202110796227.1A patent/CN113705335B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919240A (zh) * | 2018-04-23 | 2018-11-30 | 东南大学 | 一种基于组稀疏结构的水声目标辐射噪声调制谱重构方法 |
CN109946518A (zh) * | 2019-03-27 | 2019-06-28 | 河南天通电力有限公司 | 基于贝叶斯方法的电力谐波信号分析方法与分析设备 |
Non-Patent Citations (1)
Title |
---|
王逸林等: "基于稀疏重构的色噪声背景下未知线谱信号估计", 《电子与信息学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115730198A (zh) * | 2022-11-10 | 2023-03-03 | 浙江衡玖医疗器械有限责任公司 | 一种超声信号tof自动提取方法 |
CN115730198B (zh) * | 2022-11-10 | 2024-03-19 | 浙江衡玖医疗器械有限责任公司 | 一种超声信号tof自动提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113705335B (zh) | 2022-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Bayesian compressive sensing for approximately sparse signals and application to structural health monitoring signals for data loss recovery | |
Luisier et al. | A new SURE approach to image denoising: Interscale orthonormal wavelet thresholding | |
Blumensath et al. | Gradient pursuits | |
Hu et al. | Compressed sensing of complex sinusoids: An approach based on dictionary refinement | |
Wu et al. | Adaptive sparsity matching pursuit algorithm for sparse reconstruction | |
CN109946518B (zh) | 基于贝叶斯方法的电力谐波信号分析方法与分析设备 | |
Carrillo et al. | Iterative hard thresholding for compressed sensing with partially known support | |
CN105807241B (zh) | 一种利用先验信息的指数信号去噪方法 | |
CN113705335B (zh) | 一种低频稀疏线谱信号时频超分辨率提取方法 | |
CN113917490A (zh) | 激光测风雷达信号去噪方法及装置 | |
Maskan et al. | Demixing sines and spikes using multiple measurement vectors | |
CN109917330B (zh) | 一种存在相位误差时基于稀疏正交匹配追踪理论的到达角估计方法 | |
CN112087235A (zh) | 基于伪逆感知字典的稀疏度自适应doa估计方法及系统 | |
Eappen et al. | Semiparametric algorithm for processing MST radar data | |
CN113392700B (zh) | 基于Fisher判别字典学习的主动声呐目标分类方法 | |
Juhlin et al. | Fast gridless estimation of damped modes | |
CN112327260B (zh) | 一种sar回波数据中脉冲式干扰信号的抑制方法和装置 | |
CN115097391A (zh) | 一种基于自动编码器的合成孔径雷达干扰抑制方法 | |
CN105303009A (zh) | 基于压缩感知与正则mfocuss的超分辨谱估计方法 | |
CN110927678B (zh) | 自适应稀疏分数阶模糊函数杂波抑制及动目标检测方法 | |
Qu et al. | Wavelet reconstruction of nonuniformly sampled signals | |
CN116680503B (zh) | 双稀疏优化阵列天线的卫星信号稳健捕获方法及相关设备 | |
Al Irkhis | Fast CS algorithm for chirp signals parameters estimation using multi-resolution DCFT atom transformation | |
Mourad et al. | Automatic threshold estimation for iterative shrinkage algorithms used with compressed sensing | |
CN114662045B (zh) | 基于框架集的p阶张量深度学习的多维地震数据去噪方法 |
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 |