【管道泄漏】LMD算法和ELMD算法管道泄漏信号处理（含SNR、MAE）对比【含Matlab源码 3450期】

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、局部分解理论研究

局部均值分解算法(Local Mean Decomposition， LMD) 作为处理非平稳随机信号的一种手段，得到了广泛应用，并成熟地应用于机械故障诊断、信号特征提取与分析等方面。LMD算法最大的特点就在于其对信号的自适应分解能力， 这种自适应能力主要体现在该方法能够通过数据自身特点，通过特定手段，将原始信号分为不同模态函数， 针对不同模态函数进一步处理。与此同时， 局部均值分解算法(LMD) 相较于模态分解的创始算法经验模态分解算法(Empirical Mode Decomposition， EMD) 而言，其具备端点效应小、迭代次数少等优势。本章将详述LMD算法基本原理及分解流程，并针对LMD存在的基本问题， 进行改进研究。

1 局部均值分解算法理论研究
2005年， JonathanS.Smith提出的一种新的非线性和非平稳信号分析方法一局部均值分解算法(LMD) ， 并应用于脑电信号的分析中， 取得不错的效果。此外， 局部均值分解在机械故障诊断中也得到良好的应用。局部均值分解可以依据信号本身的特征进行自适应分解的，产生具有真实物理意义的PF分量，并由此得到能够清晰准确反映出信号能量在空间各尺度上分布规律的时频分布，有利于更加细致的对信号特征进行分析。
1.1局部均值算法(LMD) 分解流程
利用LMD分解， 可以将原始信号分解并产生若干有效PF分量， 将所有PF分量相加即可重构出原始信号。其中每个PF分量都是一个纯调频信号和包络信号的乘积，且每个PF分量的瞬时频率具有实际物理意义。LMD分解的具体流程如下：(1)原始信号x(t)，找出x(t)上的所有局部极值点n，由相邻的两个极值点n，n.计算出的一个均值m，即

随后将得到的所有平均值m，用折线连接，连接过程采用滑动平均方法进行平滑处理，进而得到局部均值函数m，(t)。同时，利用相邻极值点计算包络估计值a，，即





通过公式(2-12)可以发现，信号重构过程中，不会出现信息的丢失，保证了信号的完整性。具体流程图如下图1.1所示

(1)总体局部均值分解算法研究
模态混叠现象的发生一方面来自于算法本身，另一方面则会受到原始信号频率特征的影响。当待处理信号确定之后，信号特征频率变不会发生改变，因此，为抑制模态混叠现象的出现，从算法本身出发，引入噪声辅助信号处理方法，即在信号中加入白噪声来平滑脉冲干扰。因为在LMD分解过程中， 需要对极值点进行处理， 得到局部均值函数及包络估计函数，极值点的分布就会影响到包络的拟合情况。如果信号的极值点分布不均匀， 则极易产生模态混叠。因此， 在LMD中， 借助高斯白噪声辅助法， 对模态混叠进行抑制， 得到总体局部均值分解(Ensemble Local Mean Decomposition ELMD) 算法。ELMD分解， 是指在LMD分解前， 将不同有限幅值的白噪声信号加入待分解信号，利用白噪声均值为零，频谱能量分布均匀的特性，使得白噪声可以均匀的分布在整个时频空间中，并且不同时间尺度的信号会自动分布到与背景噪声相关的适当尺度上去。对于单次试验，由于噪声的添加，使得每次结果都会受到噪声的影响产生偏差，这是由于在分解过程中，信号包括原始信号及附加的白噪声。但随着试验的次数增加，由于零均值的特性，噪声将会相互抵消，进而得到消除，唯一持久稳固的部分便是信号本身， 因此可以认定集成均值的结果就是最终分解结果， 即EL MD算法。
总体局部均值EL MD分解算法流程图如图2.13所示， 即对原始信号分别加入n组不同高斯白噪声， 分别进行LMD分解， 随后将得到的n组PF分量进行平均处理，得到最终的LMD分解结果。

⛄二、部分源代码

clear;clc;close all;
%x1=xlsread(‘signal.xlsx’);
%VarName2=x1(:,2);
%x=VarName2’;
%%t=linspace(0,1,10110);
%N = length(x);
%fs=1000;
%Ts=1/fs;
%t=0:Ts:N*Ts;
%t=t(1:N);

T = 5000 ;
t = (1:T)/T;
fs=1/T
freqs = 2pi(t-0.5-1/T)/(fs);
os1=cos(30pit)
os2=cos(45pit)+0.8
os=os1+os2
x=awgn(os,20,‘measured’,‘dB’);
T=length(x);

%绘制仿真信号和其频谱图
figure(1)
subplot(211)
plot(t,x)
subplot(212)
y2=x;
L=length(y2);
NFFT = 2^nextpow2(L);
Y = fft(y2,NFFT)/L;
f = fs/2linspace(0,1,NFFT/2);
plot(f,2abs(Y(1:NFFT/2)))
% ELMD分析
Nstd=0.4;
NR=8;
modos = mlmd(x,Nstd,NR);
line=size(modos,1);
NN = length(modos(1,:));
n = linspace(0,1,NN);
figure(‘Color’,[1 1 1]);
subplot(line+1,1,1),plot(t,x),ylabel(‘X(t)’);
subplot(line+1,1,2),plot(t,modos(1,:)),ylabel(‘PF_1(t)’);
subplot(line+1,1,3),plot(t,modos(2,:)),ylabel(‘PF_2(t)’);
subplot(line+1,1,4),plot(t,modos(3,:)),ylabel(‘PF_3(t)’);
subplot(line+1,1,5),plot(t,modos(4,:)),ylabel(‘PF_4(t)’);
subplot(line+1,1,6),plot(t,modos(5,:)),ylabel(‘PF_5(t)’);
subplot(line+1,1,7),plot(t,modos(6,:)),ylabel(‘PF_6(t)’);
xlabel(‘Time / s’);

for i = 1:1:line
cc(i)=min(min(corrcoef(modos(i,:), x)));
end

figure
plot(cc,‘-g<’,‘LineWidth’,1.5,‘MarkerEdgeColor’,‘b’,‘MarkerFaceColor’,‘b’,‘MarkerSize’,5);
set(gca,‘XGrid’, ‘on’, ‘YGrid’, ‘on’);
legend(‘CC’);
xlabel(‘IMF’);
ylabel(‘相关系数’);

%重构信号
cg_ev=modos(4,:)+modos(5,:)+modos(3,:);
figure;
subplot(311);plot(t,os);xlabel(‘t/s’);ylabel(‘幅值’);legend(‘纯净信号’);
hold on
subplot(312);plot(t,x);xlabel(‘t/s’);ylabel(‘幅值’);legend(‘含噪信号’);
hold on
subplot(313);
plot(t,cg_ev);xlabel(‘t/s’);ylabel(‘幅值’);legend(‘重构信号’);

%误差信号
err=os-cg_ev;
figure
plot(t,err);xlabel(‘t/s’);ylabel(‘幅值’);legend(‘LMD误差信号’);

% ELMD分解
modes = mlmd(x,Nstd,NR);
[a, b]=size(modes);

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a或2019b

2 参考文献
[1]李博健.改进LMD算法在管道泄漏中的应用研究[D].东北石油大学

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合