CN109871733A - 一种自适应海杂波信号去噪方法 - Google Patents
一种自适应海杂波信号去噪方法 Download PDFInfo
- Publication number
- CN109871733A CN109871733A CN201811130695.XA CN201811130695A CN109871733A CN 109871733 A CN109871733 A CN 109871733A CN 201811130695 A CN201811130695 A CN 201811130695A CN 109871733 A CN109871733 A CN 109871733A
- Authority
- CN
- China
- Prior art keywords
- imf
- component
- signal
- sea clutter
- rank
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种自适应海杂波信号去噪方法,属于雷达信号处理领域。该方法包括如下步骤:1)对含噪海杂波信号x(n)进行自适应噪声完备互补集成经验模态分解,得到m+1个IMF分量;2)计算x(n)与各IMF分量之间的互相关系数,并确定互相关系数中第一个局部极小值所对应的IMF分量imfk(n);3)利用ICA对识别出前k个含噪分量进行去噪,重构去噪后的各IMF分量,得到消噪后的信号xd(n);4)对消噪后的信号xd(n)再次进行CEEMDAN分解,并根据步骤2)确定出互相关系数中一个局部极小值所在的层数n;5)判断k是否等于n,如果不相等,则将n赋值给k,并跳转到步骤3)继续向下执行;否则重构各IMF模态分量,输出最终消噪后的海杂波信号。本方法能够在保留海杂波原有混沌特性的基础上有效地抑制其中的噪声。
Description
技术领域
本发明涉及一种自适应海杂波信号去噪方法,属于雷达信号处理领域。
背景技术
海杂波是指雷达照射海面的后向散射回波,其产生主要受到海风、海浪、潮 汐等不同因素的影响。研究表明,海杂波具有混沌特性,是一种典型的非平稳信 号,利用海杂波的混沌特性可以监视海面情况,实现海面低小目标的检测,因此 海杂波背景下的微弱信号检测一直是信号处理领域的研究热点之一。而实际的海 杂波数据通常会受到诸如雷达测量噪声和海面动态噪声的影响,这些噪声的引入 无疑增加了目标信号的检测难度,传统的频谱分析方法已不再适用于海杂波的去 噪,所以研究一种有效的去噪方法对于分析海杂波物理特性和提高信号检测精度 具有重要的意义。
目前与本发明相关的报道有:1“Uncovering nonlinear dynamics-the casestudy of sea clutter”(Proceedings of the IEEE,2002,90(5):860-881)在分 析海杂波混沌特性的基础上提出了一种保障海杂波物理特性不受影响的去噪方 法。2“EMD-BasedSignal Filtering”(IEEE transactions on instrumentation and measurement,2007,56(6):2196-2202)介绍了利用经验模态分解 (Empirical Mode Decomposition,EMD)法对海杂波信号进行分解,并采用不同 的阈值对各固有模态函数(Intrinsic Mode Function,IMF)进行去噪,较好地抑制 了噪声的干扰。3“基于EMD算法的海杂波信号去噪”(《雷达科学与技术》2010 年第8卷第2期)针对实测海杂波信号非线性非平稳的特点,利用小波分解实现了海杂波的去噪,但小波分解方式的不唯一影响了海杂波的去噪效果,不具 普适性。4“一种基于集成经验模态分解的海杂波去噪方法”(CN201510209153.1) 该方法将集成经验模态分解算法用于海杂波去噪处理中,采用自相关函数分选 有用信号与噪声分量,选用SG(Savitzky-Golay)滤波对噪声模态分量进行削 噪。
目前,基于传统的基于经验模态分解去噪的方法大多是将前几层高频IMF (固有模态函数)分量当做不必要的噪声而滤除,但随着EMD(经验模态分解) 理论体系不断发展,已有研究证实(“Development of EMD-Based Denoising Methods Inspired by WaveletThresholding”(IEEE Transactions on Signal Processing, 2009,57(4):1351-1362.))高频IMF分量中仍含有较多的细节成分,如果简单的 对这些分量进行舍弃会造成细节信息丢失较多,且噪声也不能得到彻底的去除。 因此如何有效地滤除高频IMF分量中噪声,最大程度地提取细节信息是提升去 噪效果的关键。由于各高频IMF分量中的噪声方差、系数分布模型等都是未知 信息,使得难以对其中的噪声分布进行准确估计。独立成分分析(independent component analysis,ICA)的出现则给IMF分量的去噪提供了新思路,该方法是一 种全新的信号处理和数据分析方法,可以在源信号和混合参数均未知的前提下,仅根据信号间的高阶统计特性恢复出源信号,非常适合于IMF分量的消噪。
发明内容
本发明提出了一种自适应海杂波信号去噪方法,克服原有模型的局限,进一 步在保留海杂波内在特性的基础上提高去噪效果。
本发明为解决其技术问题采用如下技术方案:
一种自适应海杂波信号去噪方法,包括如下步骤:
1)对含噪海杂波信号x(n)进行CEEMDAN,得到m+1个IMF分量 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)},imf1(n)为x(n)经CEEMDAN分解得到的第一阶 IMF分量,imf2(n)为x(n)经CEEMDAN分解得到的第二阶IMF分量,imfm(n)为x(n)经 CEEMDAN分解得到的第m阶IMF分量,imfm+1(n)为x(n)经CEEMDAN分解得到的第 m+1阶IMF分量;
2)计算原信号x(n)与分解所得各IMF分量 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)}之间的互相关系数,并在这些互相关系数中 寻找第一个局部极小值所对应的IMF分量imfk(n),即为有用信号与噪声的分界, 即临界模态分量,k<m+1;
3)利用ICA对识别出前k个含噪分量{imf1(n),imf2(n),...,imfk(n)}进行去噪,则第i个含噪分量去噪后的信号记为imfi d(n),i=1,2,....,k,重构去噪后的IMF 分量和剩余分量,得到消噪后的信号xd(n),则有:
其中:imfi(n)为x(n)经CEEMDAN分解后所得的第i个IMF分量;
4)对消噪重构后的信号xd(n)再次进行CEEMDAN分解,并根据步骤2)确 定出互相关系数中一个局部极小值所在的层数n;
5)判断k是否等于n,如果不相等,则将n赋值给k,并跳转到步骤3)继 续向下执行;否则重构各IMF模态分量,输出最终消噪后的海杂波信号。
所述步骤1)的具体过程如下:
Step 1产生N个与原海杂波信号x(n)等长的高斯白噪声序列,并对这N个 白噪声进行EMD分解,取每个噪声分解后的第一阶IMF分量作为新噪声加入到 x(n)中,则有:
xi(n)=x(n)+β0E1[wi(n)] (2)
其中,xi(n)为第i次加入噪声后的信号,wi(n)为产生的第i个白噪声序列, E1[wi(n)]为对wi(n)进行EMD分解后得到的第一阶IMF模态分量,β0为常数, i=1,2,3......N;
Step 2找出xi(n)的所有极值点,利用三次样条插值函数分别拟合由极值点构 成的上下两条包络线,并计算上下包络线的平均值,记为M[xi(n)];对这N组包 络线均值做平均处理,得到第一阶残差分量r1(n),则有:
Step 3计算第一阶IMF分量c1(n),则有:
c1(n)=x(n)-r1(n) (4)
Step 4将r1(n)作为原始信号,分别向其中加入N个白噪声分解后的第二阶 IMF分量,重复Step 2、3求得原信号的第二阶IMF分量c2(n),则有:
其中:β1为常数,E2[wi(n)]为对wi(n)进行EMD分解后得到的第二阶IMF模 态分量,M[]为包络线平均值算子;
Step 5重复Step 2-4分别求得第3、4......m阶的残差分量rm(n),则有:
其中:rm-1(n)为第m-1阶残差分量,βm-1为为常数,Em[wi(n)]为对wi(n)进行 EMD分解后得到的第m阶IMF模态分量;
Step 6根据残差分量rm(n)进一步求出第m阶IMF分量cm(n),则有:
cm(n)=rm-1(n)-rm(n) (7)
Step 7最终得到m+1个IMF分量{c1(n),c2(n),......,cm(n),rm(n)},记为 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)}。
所述步骤2)的具体过程如下:
按公式(8)计算原信号x(n)与各IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)} 之间的互相关系数:
式中,R(x(n),imfi(n))为x(n)与imfi(n)之间的互相关系数,imfi(n)为第i个IMF分 量,N为原始信号长度,且
寻找R(x(n),imfi(n))中第一个局部极小值所对应的IMF分量imfk(n),即为临 界模态分量,k<m+1。
所述步骤3)的具体过程如下:
a)原含噪海杂波信号CEEMDAN分解得到m+1个IMF分量 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)},利用互相关系数确定出前k个含噪分量 {imf1(n),imf2(n),...,imfk(n)},并令:
其中:xp(n)为第二阶IMF分量到第m+1累加和,imfi(n)为第i个IMF分量;
b)对imf1(n)分量进行P次循环移位,构造出P个噪声信号 {imf1 1(n),imf1 2(n),...,imf1 P(n)},其中:imf1 1(n)为imf1(n)经一次循环移位得到的分 量,imf1 2(n)为imf1(n)经两次循环移位得到的分量,imf1 P(n)为imf1(n)经P次循环 移位得到的分量;
c)将imf1 α(n)与xp(n)分别累加,重构出P个新的含噪信号:
xα(n)=xp(n)+imf1 α(n),α=1,2,....,P (10)
其中:imf1 α(n)为构造出P个新的含噪信号,xα(n)为重构出的第α个新含噪信号;
d)分别对xα(n)进行m+1层CEEMDAN分解,得到分解后的模态分量: {imf′1,α(n),imf′2,α(n),...,imf′m+1,α(n)},其中:imf′1,α(n)为xα(n)经CEEMDAN分解得到的 第一阶IMF分量,imf′2,α(n)为xα(n)经CEEMDAN分解得到的第二阶IMF分量, imf′m+1,α(n)为xα(n)经CEEMDAN分解得到的第m+1阶IMF分量;
e)利用ICA对前k个含噪IMF分量进行去噪,其中,第i层含噪分量imfi(n)去噪时的输入样 本为{imf′i,1(n),imf′i,2(n),...,imf′i,P(n)},其中:imf′i,P(n)为xP(n)经CEEMDAN分解得到 的第i阶IMF分量,xP(n)就是xα(n)中的第P个,因为α=1,2,....,P,去噪后的信号记为 imfi d(n),i=1,2,....,k;
f)重构去噪后的IMF分量和剩余分量,得到消噪后的信号xd(n)。
所述消噪后的信号xd(n):
其中:imfi d(n)为k个含噪IMF去噪后的信号,imfi(n)为x(n)经CEEMDAN 分解得到的第i阶IMF分量,i=1,...,m+1,k为确定出前k个含噪分量。
本发明的有益效果如下:
为了进一步提高海杂波的去噪效果,本申请结合CEEMDAN(自适应噪声完 备互补集成经验模态分解)分解理论和独立成分分析的去噪特点,提出了基于各 IMF分量互相关系数分布的自适应去噪方法。对含有噪声的混沌海杂波信号进行 CEEMDAN分解,根据分解所得各IMF分量与原信号互相关系数的第一个局部 极小值所在层数别出信号和噪声的分解,并对含噪分量自适应地进行ICA消噪。 实验结果表明,该方法能够在保留海杂波原有混沌特性的基础上有效地抑制其中 的噪声,与其它去噪方法相比,在均方根误差(Root MeanSquared Error,RMSE) 指标上的提升效果明显。
附图说明
图1为一种自适应海杂波信号去噪方法的流程图。
图2(a)为17#海杂波信号去噪前的真实值和预测值;图2(b)为17#海杂波 信号去噪前的预测误差。
图3(a)为17#海杂波信号经EMD模态阈值法去噪后的真实值和预测值;图 3(b)为17#海杂波信号经EMD模态阈值法去噪后的预测误差。
图4(a)为17#海杂波信号经小波阈值法去噪后的真实值和预测值,图4(b) 为17#海杂波信号经小波阈值法去噪后的预测误差。
图5(a)为17#海杂波信号经EEMD-SG法去噪后的真实值和预测值,图5(b) 为17#海杂波信号经EEMD-SG法去噪后的预测误差。
图6(a)为17#海杂波信号经一种自适应海杂波信号去噪方法去噪后的真实 值和预测值,图6(b)为17#海杂波信号经一种自适应海杂波信号去噪方法去噪 后的预测误差。
具体实施方式
下面结合附图对本发明做进一步详细说明。
混合去噪算法流程如图1所示,具体步骤如下:
步骤1:对含噪海杂波信号x(n)进行CEEMDAN分解,得到一系列IMF 模态分量。
步骤2:寻找各IMF分量与原始含噪海杂波信号互相关系数中第一个局 部极小值所在层数k,并基于各分量的互相关特性识别出含噪IMF分量。
步骤3:利用快速独立成分分析(fast independent component analysis,FastICA)对前k个含噪分量进行去噪,并重构。
步骤4:对重构后的信号再次进行CEEMDAN分解。
步骤5:求得互相关系数中的第一个局部极小值所在层数。
步骤6:判断是否k等于n,如条件不成立则k=n,重复上述步骤3、4、 5直至条件满足,否则跳转到步骤6。
步骤7:重构各IMF模态分量,输出消噪后的海杂波信号。
所述步骤1)的具体实施为:
Step 1产生N个与原海杂波信号x(n)等长的高斯白噪声序列,并对这N个 白噪声进行EMD分解,取每个噪声分解后的第一阶IMF分量作为新噪声加入到 x(n)中,则有:
xi(n)=x(n)+β0E1[wi(n)] (2)
其中,xi(n)为第i次加入噪声后的信号,wi(n)为产生的第i个白噪声序列, E1[wi(n)]为对wi(n)进行EMD分解后得到的第一阶IMF模态分量,β0为常数, i=1,2,3......N。
Step 2找出xi(n)的所有极值点,利用三次样条插值函数分别拟合由极值点构 成的上下两条包络线,并计算上下包络线的平均值,记为M[xi(n)]。对这N组包 络线均值做平均处理,得到第一阶残差分量r1(n),则有:
Step 3计算第一阶IMF分量c1(n),则有:
c1(n)=x(n)-r1(n) (4)
Step 4将r1(n)作为原始信号,分别向其中加入N个白噪声分解后的第二阶 IMF分量,重复Step 2、3求得原信号的第二阶IMF分量c2(n),则有:
其中:β1为常数,E2[wi(n)]为对wi(n)进行EMD分解后得到的第二阶IMF模 态分量,M[]为包络线平均值算子。
Step 5重复Step 2-4分别求得第3、4......m阶的残差分量rm(n),则有:
其中:rm-1(n)为第m-1阶残差分量,βm-1为为常数,Em[wi(n)]为对wi(n)进行 EMD分解后得到的第m阶IMF模态分量,
Step 6根据残差分量rm(n)进一步求出第m阶IMF分量cm(n),则有:
cm(n)=rm-1(n)-rm(n) (7)
Step 7最终得到m+1个IMF分量{c1(n),c2(n),......,cm(n),rm(n)},记为 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)}。
所述步骤2)的具体实施为:
按公式(8)计算原信号x(n)与各IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)} 之间的互相关系数:
式中,imfi(n)为第i个IMF分量,N为原始信号长度,且
寻找R(x(n),imfi(n))中第一个局部极小值所对应的IMF分量imfk(n),即为临 界模态分量,k<m+1。
所述步骤3)的具体实施为:
Step 1原含噪海杂波信号CEEMDAN分解得到m+1个IMF分量 {imf1(n),imf2(n),...,imfm(n),imfm+1(n)},利用互相关系数确定出前k个含噪分量 {imf1(n),imf2(n),...,imfk(n)},并令:
Step 2对imf1(n)分量进行P次循环移位,构造出P个噪声信号 {imf1 1(n),imf1 2(n),...,imf1 P(n)}。
Step 3将imf1 α(n)与xp(n)分别累加,重构出P个新的含噪信号:
xα(n)=xp(n)+imf1 α(n),α=1,2,....,P (10)
其中:imf1 α(n)为构造出P个新的含噪信号。
Step 4分别对xα(n)进行m+1层CEEMDAN分解,得到分解后的模态分量: {imf′1,α(n),imf′2,α(n),...,imf′m+1,α(n)}。
Step 5利用ICA对前k个含噪IMF分量进行去噪,其中,第i层含噪分量 imfi(n)去噪时的输入样本为{imf′i,1(n),imf′i,2(n),...,imf′i,P(n)},含噪IMF分量去噪 后的信号记为imfi d(n),i=1,2,....,k。
Step 6重构去噪后的IMF分量和剩余分量,得到消噪后的信号xd(n)。
此外,本申请去噪过程中涉及的ICA算法采用的是FastICA算法。
所述步骤4)的具体实施为:
对消噪重构后的信号xd(n)再次进行CEEMDAN分解,并根据步骤2)确定出 互相关系数中一个局部极小值所所在的层数n。
所述步骤5)的具体实施为:
判断k是否等于n,如果不相等,则将n赋值给k,并跳转到步骤3)继续向 下执行;否则重构各IMF模态分量,输出最终消噪后的海杂波信号。
为了说明本申请所提方法的有效性,本实例选用加拿大McMaster大学的 IPIX雷达17#海杂波目标距离单元(目标位置在300-400和550-560处)的2900个 数据点(VV极化)进行验证。将前1900个数据点设为训练样本,后1000个数据 点设为测试样本,利用RBF(径向基函数)神经网络分别对去噪前的原始海杂 波数据和经EMD模态阈值法、小波阈值法、EEMD-SG法(专利, CN201510209153.1)以及本文方法去噪后的海杂波数据进行单步预测,通过预测 误差判定微弱目标存在情况,检测效果如图2(a)、图2(b),图3(a)、图3(b),图 4(a)、图4(b),图5(a)、图5(b),图6(a)、图6(b)所示。
从图2(a)、图2(b),图3(a)、图3(b),图4(a)、图4(b),图5(a)、图5(b),图6(a)、 图6(b)可以看出,在未经任何消噪处理的原始海杂波数据中,微弱信号已完全湮 没在噪声背景下,预测误差不存在明显突起,无法检测到微弱信号的存在。利用 RBF神经网络对经EMD模态阈值法消噪后的海杂波数据进行预测,预测误差仅 在n=300-400处出现了明显的尖峰,但鉴于该方法对噪声的抑制并不彻底, n=550-560处的微弱信号依旧未检测出来。而经小波阈值法,EEMD-SG(EEMD 分解和SG滤波联合去噪)法和本申请方法去噪后的预测误差在n=300-400处和 n=550-560处均出现了明显尖峰,可以初步判定海杂波背景中微弱信号的存在。
表1 17#海杂波去噪性能对比
由表1可知,本申请所提方法去噪后的RMSE(6.3558*10-4)较未经任何去噪处 理的RMSE(0.0058)提高了近一个数量级,与其它去噪方法相比,在RMSE上也 提高了47.04%-69.73%,进一步表明本申请所提去噪算法能够在保障海杂波内在 物理特征的基础上有效地去除噪声干扰,对于提高海杂波背景下微弱信号的检测 精度,具有很高的实用价值。
Claims (5)
1.一种自适应海杂波信号去噪方法,其特征在于,包括如下步骤:
1)对含噪海杂波信号x(n)进行CEEMDAN,得到m+1个IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)},imf1(n)为x(n)经CEEMDAN分解得到的第一阶IMF分量,imf2(n)为x(n)经CEEMDAN分解得到的第二阶IMF分量,imfm(n)为x(n)经CEEMDAN分解得到的第m阶IMF分量,imfm+1(n)为x(n)经CEEMDAN分解得到的第m+1阶IMF分量;
2)计算原信号x(n)与分解所得各IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)}之间的互相关系数,并在这些互相关系数中寻找第一个局部极小值所对应的IMF分量imfk(n),即为有用信号与噪声的分界,即临界模态分量,k<m+1;
3)利用ICA对识别出前k个含噪分量{imf1(n),imf2(n),...,imfk(n)}进行去噪,则第i个含噪分量去噪后的信号记为imfi d(n),i=1,2,....,k,重构去噪后的IMF分量和剩余分量,得到消噪后的信号xd(n),则有:
其中:imfi(n)为x(n)经CEEMDAN分解后所得的第i个IMF分量;
4)对消噪重构后的信号xd(n)再次进行CEEMDAN分解,并根据步骤2)确定出互相关系数中一个局部极小值所在的层数n;
5)判断k是否等于n,如果不相等,则将n赋值给k,并跳转到步骤3)继续向下执行;否则重构各IMF模态分量,输出最终消噪后的海杂波信号。
2.根据权利要求1所述的一种自适应海杂波信号去噪方法,其特征在于,所述步骤1)的具体过程如下:
Step1产生N个与原海杂波信号x(n)等长的高斯白噪声序列,并对这N个白噪声进行EMD分解,取每个噪声分解后的第一阶IMF分量作为新噪声加入到x(n)中,则有:
xi(n)=x(n)+β0E1[wi(n)] (2)
其中,xi(n)为第i次加入噪声后的信号,wi(n)为产生的第i个白噪声序列,E1[wi(n)]为对wi(n)进行EMD分解后得到的第一阶IMF模态分量,β0为常数,i=1,2,3......N;
Step2找出xi(n)的所有极值点,利用三次样条插值函数分别拟合由极值点构成的上下两条包络线,并计算上下包络线的平均值,记为M[xi(n)];对这N组包络线均值做平均处理,得到第一阶残差分量r1(n),则有:
Step3计算第一阶IMF分量c1(n),则有:
c1(n)=x(n)-r1(n) (4)
Step4将r1(n)作为原始信号,分别向其中加入N个白噪声分解后的第二阶IMF分量,重复Step2、3求得原信号的第二阶IMF分量c2(n),则有:
其中:β1为常数,E2[wi(n)]为对wi(n)进行EMD分解后得到的第二阶IMF模态分量,M[]为包络线平均值算子;
Step5重复Step2-4分别求得第3、4......m阶的残差分量rm(n),则有:
其中:rm-1(n)为第m-1阶残差分量,βm-1为为常数,Em[wi(n)]为对wi(n)进行EMD分解后得到的第m阶IMF模态分量,
Step6根据残差分量rm(n)进一步求出第m阶IMF分量cm(n),则有:
cm(n)=rm-1(n)-rm(n) (7)
Step7最终得到m+1个IMF分量{c1(n),c2(n),......,cm(n),rm(n)},记为{imf1(n),imf2(n),...,imfm(n),imfm+1(n)}。
3.根据权利要求1所述的一种自适应海杂波信号去噪方法,其特征在于,所述步骤2)的具体过程如下:
按公式(8)计算原信号x(n)与各IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)}之间的互相关系数:
式中,R(x(n),imfi(n))为x(n)与imfi(n)之间的互相关系数,imfi(n)为第i个IMF分量,N为原始信号长度,且i=1,...,m+1;
寻找R(x(n),imfi(n))中第一个局部极小值所对应的IMF分量imfk(n),即为临界模态分量,k<m+1。
4.根据权利要求1所述的一种自适应海杂波信号去噪方法,其特征在于,所述步骤3)的具体过程如下:
a)原含噪海杂波信号CEEMDAN分解得到m+1个IMF分量{imf1(n),imf2(n),...,imfm(n),imfm+1(n)},利用互相关系数确定出前k个含噪分量{imf1(n),imf2(n),...,imfk(n)},并令:
其中:xp(n)为第二阶IMF分量到第m+1累加和,imfi(n)为第i个IMF分量;
b)对imf1(n)分量进行P次循环移位,构造出P个噪声信号{imf1 1(n),imf1 2(n),...,imf1 P(n)},其中:imf1 1(n)为imf1(n)经一次循环移位得到的分量,imf1 2(n)为imf1(n)经两次循环移位得到的分量,imf1 P(n)为imf1(n)经P次循环移位得到的分量;
c)将imf1 α(n)与xp(n)分别累加,重构出P个新的含噪信号:
xα(n)=xp(n)+imf1 α(n),α=1,2,....,P (10)
其中:imf1 α(n)为构造出P个新的含噪信号,xα(n)为重构出的第α个新含噪信号;
d)分别对xα(n)进行m+1层CEEMDAN分解,得到分解后的模态分量:{imf'1,α(n),imf'2,α(n),...,imf'm+1,α(n)},其中:imf'1,α(n)为xα(n)经CEEMDAN分解得到的第一阶IMF分量,imf'2,α(n)为xα(n)经CEEMDAN分解得到的第二阶IMF分量,imf'm+1,α(n)为xα(n)经CEEMDAN分解得到的第m+1阶IMF分量;
e)利用ICA对前k个含噪IMF分量进行去噪,其中,第i层含噪分量imfi(n)去噪时的输入样本为{imf'i,1(n),imf'i,2(n),...,imf'i,P(n)},其中:imf'i,P(n)为xP(n)经CEEMDAN分解得到的第i阶IMF分量,xP(n)就是xα(n)中的第P个,因为α=1,2,....,P,去噪后的信号记为imfi d(n),i=1,2,....,k;
f)重构去噪后的IMF分量和剩余分量,得到消噪后的信号xd(n)。
5.根据权利要求4所述的一种自适应海杂波信号去噪方法,其特征在于,所述消噪后的信号xd(n):
其中:imfi d(n)为k个含噪IMF分量去噪后的信号,imfi(n)为x(n)经CEEMDAN分解得到的第i阶IMF分量,i=1,...,m+1,k为确定出前k个含噪分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811130695.XA CN109871733B (zh) | 2018-09-27 | 2018-09-27 | 一种自适应海杂波信号去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811130695.XA CN109871733B (zh) | 2018-09-27 | 2018-09-27 | 一种自适应海杂波信号去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109871733A true CN109871733A (zh) | 2019-06-11 |
CN109871733B CN109871733B (zh) | 2023-05-19 |
Family
ID=66916915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811130695.XA Active CN109871733B (zh) | 2018-09-27 | 2018-09-27 | 一种自适应海杂波信号去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109871733B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN112098516A (zh) * | 2020-09-29 | 2020-12-18 | 国家电网有限公司 | 一种超声检测用传感器及其信号处理方法 |
CN112766127A (zh) * | 2021-01-12 | 2021-05-07 | 南京信息工程大学 | 基于互补集合模态分解和sg滤波的雷云点电荷定位方法 |
CN112800831A (zh) * | 2020-12-24 | 2021-05-14 | 长安大学 | 一种用于时变重力场的emd滤波方法及系统 |
CN113158769A (zh) * | 2021-03-03 | 2021-07-23 | 安徽大学 | 基于CEEMDAN和FastICA的机电设备轴承振动信号除噪方法 |
CN113378485A (zh) * | 2021-08-12 | 2021-09-10 | 南京天朗防务科技有限公司 | 一种基于meemd与sg滤波的海杂波去噪方法及装置 |
CN114264724A (zh) * | 2021-12-23 | 2022-04-01 | 大连海洋大学 | 基于ceemdan联合小波包阈值的超声信号处理方法 |
CN114936571A (zh) * | 2022-04-01 | 2022-08-23 | 西南交通大学 | 噪声抑制方法、装置、设备及存储介质 |
CN117648557A (zh) * | 2024-01-30 | 2024-03-05 | 山东科技大学 | 一种基于soh联合降噪的soh预测方法及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765979A (zh) * | 2015-04-28 | 2015-07-08 | 南京信息工程大学 | 一种基于集成经验模态分解的海杂波去噪方法 |
CN106228529A (zh) * | 2016-09-05 | 2016-12-14 | 上海理工大学 | 一种激光散斑图像处理分析方法 |
CN107272066A (zh) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | 一种含噪地震信号初至走时拾取方法及装置 |
-
2018
- 2018-09-27 CN CN201811130695.XA patent/CN109871733B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765979A (zh) * | 2015-04-28 | 2015-07-08 | 南京信息工程大学 | 一种基于集成经验模态分解的海杂波去噪方法 |
CN106228529A (zh) * | 2016-09-05 | 2016-12-14 | 上海理工大学 | 一种激光散斑图像处理分析方法 |
CN107272066A (zh) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | 一种含噪地震信号初至走时拾取方法及装置 |
Non-Patent Citations (1)
Title |
---|
罗志增等: "基于CEEMDAN-ICA的单通道脑电信号眼电伪迹滤除方法", 《HTTP://KNS.CNKI.NET/KCMS/DETAIL/32.1322.TN.20180829.1831.030.HTML》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989005B (zh) * | 2019-12-23 | 2021-12-28 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN112098516B (zh) * | 2020-09-29 | 2022-10-04 | 国家电网有限公司 | 一种超声检测用传感器及其信号处理方法 |
CN112098516A (zh) * | 2020-09-29 | 2020-12-18 | 国家电网有限公司 | 一种超声检测用传感器及其信号处理方法 |
CN112800831A (zh) * | 2020-12-24 | 2021-05-14 | 长安大学 | 一种用于时变重力场的emd滤波方法及系统 |
CN112766127A (zh) * | 2021-01-12 | 2021-05-07 | 南京信息工程大学 | 基于互补集合模态分解和sg滤波的雷云点电荷定位方法 |
CN112766127B (zh) * | 2021-01-12 | 2023-08-22 | 南京信息工程大学 | 基于互补集合模态分解和sg滤波的雷云点电荷定位方法 |
CN113158769A (zh) * | 2021-03-03 | 2021-07-23 | 安徽大学 | 基于CEEMDAN和FastICA的机电设备轴承振动信号除噪方法 |
CN113378485A (zh) * | 2021-08-12 | 2021-09-10 | 南京天朗防务科技有限公司 | 一种基于meemd与sg滤波的海杂波去噪方法及装置 |
CN114264724A (zh) * | 2021-12-23 | 2022-04-01 | 大连海洋大学 | 基于ceemdan联合小波包阈值的超声信号处理方法 |
CN114936571A (zh) * | 2022-04-01 | 2022-08-23 | 西南交通大学 | 噪声抑制方法、装置、设备及存储介质 |
CN114936571B (zh) * | 2022-04-01 | 2023-05-05 | 西南交通大学 | 噪声抑制方法、装置、设备及存储介质 |
CN117648557A (zh) * | 2024-01-30 | 2024-03-05 | 山东科技大学 | 一种基于soh联合降噪的soh预测方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN109871733B (zh) | 2023-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109871733A (zh) | 一种自适应海杂波信号去噪方法 | |
CN105913393B (zh) | 一种自适应小波阈值图像去噪方法及装置 | |
CN107688120A (zh) | 基于模糊熵的含噪信号处理方法及迭代奇异谱软阈值去噪方法 | |
CN102697495B (zh) | 基于总体平均经验模式分解的二代小波肌电信号消噪方法 | |
CN104765979A (zh) | 一种基于集成经验模态分解的海杂波去噪方法 | |
CN104636609A (zh) | 一种基于经验模态分解与小波分析的信号联合去噪方法 | |
CN107144829A (zh) | 一种高效的激光雷达回波信号去噪方法 | |
CN113378661A (zh) | 一种基于改进小波阈值和相关检测的直流电能信号去噪方法 | |
CN103995950A (zh) | 基于空域相关修正阈值的变小波系数局部放电信号消噪方法 | |
CN109709585B (zh) | 去除gps坐标时间序列中有色噪声的方法 | |
CN113537112B (zh) | 基于改进小波阈值及vmd的激光超声表面缺陷成像降噪方法 | |
CN105676205A (zh) | 一种机载LiDAR波形数据高斯分解方法 | |
CN104635223A (zh) | 基于经验模式分解与分数阶傅里叶变换的激光回波去噪方法 | |
CN113887398A (zh) | 一种基于变分模态分解和奇异谱分析的gpr信号去噪方法 | |
CN101685158A (zh) | 基于隐马尔科夫树模型的sar图像去噪方法 | |
CN108802175A (zh) | 一种基于emd的脉冲涡流检测信号降噪方法 | |
CN103576131B (zh) | 距离扩展目标内禀模态特征能量智能融合检测方法 | |
CN105809649A (zh) | 基于变分多尺度分解的sar图像与可见光图像融合方法 | |
CN113589253A (zh) | 一种基于伪时域的小波变换算法对微弱回波信号进行检测的方法 | |
Li et al. | Magnetotelluric signal-noise separation method based on SVM–CEEMDWT | |
CN112595782A (zh) | 一种基于eemd算法的超声波横波起跳点识别方法及系统 | |
CN103605120A (zh) | 距离扩展目标内禀模态局部奇异值有序融合检测方法 | |
CN101930605B (zh) | 基于二维混合变换的sar图像目标提取方法及系统 | |
CN106815434A (zh) | 基于小波和emd的噪声量时频分析去噪效果评价方法 | |
CN116383605A (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 | ||
CB02 | Change of applicant information |
Address after: 210032 No. 219 Ning six road, Jiangbei new district, Nanjing, Jiangsu Applicant after: Nanjing University of Information Science and Technology Address before: 211500 Yuting Square, 59 Wangqiao Road, Liuhe District, Nanjing City, Jiangsu Province Applicant before: Nanjing University of Information Science and Technology |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |