局部均值分解(local meanvdecomposition，LMD)是 Smith 等人于2005年提出，该算法可以自适应将复杂的非平稳信号分解成一系列的乘积函数(product function，PF)，每个 PF函数有明确的物理意义(由包络信号和纯调频信号相乘得到)，包络信号即瞬时幅值函数，纯调频信号可以求得瞬时频率。

一、LMD算法

1、基本原理

LMD算法的具体步骤如下：
1）找出输入信号 x(t)x(t)x(t) 的所有局部极值点 nin_{i}ni，计算相邻局部极值点的均值 mim_{i}mi 和包络估计值 aia_{i}ai，计算公式为：
mi=(ni+ni+1)/2；ai=∣ni−ni+1∣/2m_{i}=(n_{i}+n_{i+1})/2； a_{i}=|n_{i}-n_{i+1}|/2 mi=(ni+ni+1)/2；ai=∣ni−ni+1∣/2
2）将 mim_{i}mi 各点用直线连接并采用移动平均方法进行平滑处理，得到局部均值函数 m11(t)m_{11}(t)m11(t)，同理，将 aia_{i}ai 各点用直线连接并采用移动平均方法进行平滑处理，得到包络估计函数 a11(t)a_{11}(t)a11(t)。

3）将局部均值函数 m11(t)m_{11}(t)m11(t)从原始信号 x(t)x(t)x(t) 中分离出得到 h11(t)h_{11}(t)h11(t)，并用包络估计函数 a11(t)a_{11}(t)a11(t) 对其进行解调，可得S11(t)S_{11}(t)S11(t)。
h11(t)=x(t)−m11(t)；S11(t)=h11(t)/a11(t)h_{11}(t)=x(t) - m_{11}(t)； S_{11}(t)=h_{11}(t)/a_{11}(t) h11(t)=x(t)−m11(t)；S11(t)=h11(t)/a11(t)
4）将 S11(t)S_{11}(t)S11(t) 视为输入信号 x(t)x(t)x(t)，重复上述步骤得到 a12(t)a_{12}(t)a12(t)。若 a12(t)≠1a_{12}(t)≠1a12(t)=1，则说明 S11(t)S_{11}(t)S11(t) 不是纯调频函数，重复上述迭代步骤 nnn 次，直到 S1n(t)S_{1n}(t)S1n(t) 为纯调频函数，即 limn→∞a1n(t)=1\underset{n→∞}{lim}a_{1n}(t)=1n→∞lima1n(t)=1。
5）将迭代过程中的所有包络估计函数相乘得到包络信号 a1(t)a_{1}(t)a1(t) ，即瞬时幅值函数。
a1(t)=a11(t)a12(t)...a1n(t)=∏k=1na1k(t)a_{1}(t)=a_{11}(t)a_{12}(t)...a_{1n}(t)=\prod ^n_{k=1}a_{1k}(t) a1(t)=a11(t)a12(t)...a1n(t)=k=1∏na1k(t)
6）将包络信号 a1(t)a_{1}(t)a1(t) 和纯调频信号 S1n(t)S_{1n}(t)S1n(t) 相乘得到原始信号 x(t)x(t)x(t) 的第一个分量PF1(t)PF_{1}(t)PF1(t)。
PF1(t)=a1(t)S1n(t)PF_{1}(t)=a_{1}(t)S_{1n}(t) PF1(t)=a1(t)S1n(t)
PF1(t)PF_{1}(t)PF1(t) 分量的瞬时幅值就是包络信号 a1(t)a_{1}(t)a1(t) ，瞬时频率为：
f1(t)=12πd[arccos(S1n(t))]dtf_{1}(t)=\frac{1}{2π}\frac{d[arccos(S_{1n}(t))]}{dt} f1(t)=2π1dtd[arccos(S1n(t))]
7）将PF1(t)PF_{1}(t)PF1(t) 分量从原始信号 x(t)x(t)x(t) 中分离出来：
u1(t)=x(t)−PF1(t)u_{1}(t)=x(t)-PF_{1}(t) u1(t)=x(t)−PF1(t)
8）将 u1(t)u_{1}(t)u1(t) 作为原始信号 x(t)x(t)x(t) 重复以上步骤，循环 kkk 次，
直到 uk(t)u_{k}(t)uk(t) 为单调函数。

利用LMD算法对信号进行分解后，各分量与原始信号的关系如下：
x(t)=∑i=1kPFi(t)+uk(t)x(t)=\sum_{i=1}^k{PF_{i}(t)+u_{k}(t)} x(t)=i=1∑kPFi(t)+uk(t)

2、算法的优缺点

LMD 方法是一种基于 EMD 方法发展而来的新的时频分解方法，它的出发点在于从原信号中更加合理地自适应分解出符合信号固有时频属性的稳态子波成分，同样属于一种经验分解方法，相对于 EMD 在端点效应、过包络、欠包络等方面有一定的优势，减少了迭代次数，提高了运算效率。但是在端点效应和模态混叠问题上仍有较大的改进空间。

二、LMD应用+MATLAB代码

利用LMD算法对信号进行分解后，通过筛选最优的IMF分量进行重构，可以实现信号的去噪；也可以将LMD与近似熵、样本熵结合，作为一种特征提取方法。

1、去噪

首先使用LMD分解，获得一系列由高频到低频分布的乘积函数（product functions，PF）；通过计算原始信号与各个PF分量之间的相关系数，确定含噪信号与有效信号之间的分界位置，将分界分量之前的分量剔除，实现初步降噪。然后，针对LMD分解结果中的残留噪声，使用SVD法，以加权能量贡献率（percent of contribution to total energy，PCTE）作为奇异值阶数的确定方法，对分界PF分量进行降噪处理，实现二次滤波。

2、特征提取

1）对机械振动信号进行 LMD 分解，得到若干具有物理意义的乘积函数分量，再通过相关性分析，筛选出与原始信号相关性最大的3个分量作为数据源，求取其近似熵作为特征向量。
2）利用 LMD 算法将经预处理的运动想象脑电信号分解为多个 PF 分量，根据左右手运动想象脑电信号所处的频段选出包含 α 节律和 β 节律的 PF 分量。采用模糊熵算法提取已选择的 PF 分量的特征，使用共空间模式算法对得到的特征进行投影映射。最后，对比投影前后的特征向量，并选取其中最具区分度的特征输入到支持向量机中用于分类。
3）提出一种基于 LMD 分解的特征提取与识别算法。对原始信号进行 LMD 分解，得到若干 PF 分量，然后通过自相关原理重构信号，并提取重构信号的样本熵特征和能量特征。

3、改进VMD分解

1）快速自适应局部均值分解（Fast and Adaptive Local Mean Decomposition，FALMD）
采用顺序统计滤波器求取信号上下包络线的均值来获得局部均值函数及包络估计函数，然后将信号分解为若干PF分量及一个残余分量。该算法一方面改变了 LMD 严格的终止条件，提高了运算速率，另一方面减少了对极值点的依赖，在一定程度上抑制了端点效应。

2）改进LMD算法在微电网电能质量扰动信号检测中的应用
从边界条件、包络估计和筛选停止准则等方面对 LMD 算法进行改进：提出四点波形曲率延拓，寻找最优匹配波形用以改善“端点效应”。采用三次样条函数插值提高计算速度，使得筛选过程更快，间接减小了“端点效应”和“模态混叠”的影响。进一步提出了自适应筛选停止准则，通过内外层循环判据确定筛选停止条件，从而抑制“模态混叠”。

3）Robust Local Mean Decomposition (RLMD)
RLMD是一种改进的局部均值分解，由一组优化策略支持。优化策略可以处理LMD中的边界条件、包络估计和筛选停止准则。它同时从混合信号中提取一组单分量信号（称为乘积函数）及其相关解调信号（即AM信号和FM信号），这是与其他自适应信号处理方法（如EMD）相比最吸引人的特征。RLMD可用于时频分析。
MATLAB代码以及对应的案例在下面的链接里：
信号分解方法之鲁棒型局部均值分解(RLMD)算法

参考文献

[1] Smith J S. The local mean decomposition and its application to EEG perception data[J]. Journal of the Royal Society Interface, 2005, 2(5): 443-454.
[2] 林江刚, 胡正新, 李晶, 等. 低转速下基于AE信号与LMD的滚动轴承故障诊断[J]. 动力工程学报, 2019, 292(04): 293-298.
[3] 董林鹭, 蒋若辰, 徐奴文, 等. 基于LMD–SVD的微震信号降噪方法研究[J]. 工程科学与技术, 2019, 51(05): 126-136.
[4] 徐艳春, 高永康, 李振兴, 等. 改进LMD算法在微电网电能质量扰动信号检测中的应用[J]. 电网技术, 2019, 422(01): 332-341.
[5] 张坤, 马朝永, 胥永刚, 等. 快速自适应局部均值分解及轴承故障诊断应用[J]. 振动工程学报, 2020.
[6] Zhiliang Liu, Yaqiang Jin, Ming J. Zuo, et al. Time-frequency representation based on robust local mean decomposition for multi-component AM-FM signal analysis[J]. Mechanical Systems and Signal Processing, 2017, 95: 468-487.

【信号分解】LMD算法及其应用相关推荐

【信号分解】基于LMD算法和ELMD算法实现管道泄漏信号处理附matlab代码
1 内容介绍在科技水平相当发达的今天,互联网+.大数据慢慢渗透进人们的生活当中,但科技的进步不仅仅要体现在生活质量水平的提高.经济的快速发展,更应该体现在对社会资源的合理利用.自建国以来,我国管道 ...
类EMD的“信号分解方法”及MATLAB实现（第六篇）——LMD
继续完善"类EMD"方法系列,本篇是继EEMD.CEEMD.CEEMDAN.VMD.ICEEMDAN后的第6篇,想要看前几种方法的点击链接可以跳转. LMD(local mean ...
二、信号分解—＞局部均值分解（LMD）学习笔记
LMD学习笔记一.总述局部均值分解算法(LMD), LMD算法最大的特点就在依据信号本身的特征对信号的自适应 ...
lmd matlab 信号处理程序,LMD经验模态分解matlab程序.doc
LMD经验模态分解matlab程序 LMD经验模态分解matlab程序--原味的曾经也用滑动平均写过LMD,其实滑动平均的EMD才是原汁原味的居于均值分解. 分享给有需要的人,程序写的不好,只是希望 ...
语音跟踪：信号分解、锁相、鸡尾酒会效应、基于PR的信号分离
NLP中关于语音的部分,其中重要的一点是语音信号从背景噪音中分离.比如在一个办公室场景中,有白天的底噪-类似于白噪音的噪音.空调的声音.键盘的啪啪声.左手边45度7米元的地方同事讨论的声音.右手边1. ...
自适应_自适应信号分解综述
真实世界中的复杂系统是由相互影响的内在机制控制着,这些内在机制在多重时间和空间尺度上运行,表现出复杂的多组成成分.多层次结构和突现性等特点,这使得理解和刻画复杂系统变得困难.一个重要的方法是分析复杂 ...
经验模式分解EMD算法原理
目录简介假设条件与原理假设条件基本原理 EMD的优缺点存在的问题简介经验模态分解(Empirical Mode Decomposition, EMD)算法是由 NE. Huang 等人提 ...
类EMD的“信号分解方法”及MATLAB实现（第七篇）——EWT
这是"类EMD"方法系列的第7篇,前几篇分别是EEMD.CEEMD.CEEMDAN.VMD.ICEEMDAN.LMD,想要看前几种方法的点击链接可以跳转. 经验小波变换(empir ...
时间序列信号处理（三）——局部均值分解LMD的python实现
1.介绍局部均值分解(LMD)作为近年来出现的一种新的自适应时频分析方法,能够依据信号的自身特点将复杂的多分量调幅调频信号分解为有限个的单分量调幅调频信号之和,进而求取瞬时频率和瞬时幅值并进行组合, ...

【信号分解】LMD算法及其应用

目录