testkdesa.m

能量算子的matlab实现程序及其详细算法

%function  testPDESA
%测试泵数据的gSE计算
clear
clc
load KsWave.txt
load Khilbert.txt
x=KsWave;
fs=5120;
Ns=length(x);
t=[0:Ns-1]/fs;

y=DESA(x);
subplot(2,1,1),plot(t,x,'g-'),grid on,title('原始数据');
subplot(2,1,2),plot(t,y,'b -'),grid on,title('能量算子包络结果/Hibert包络结果')
hold on , plot(t,Khilbert,'k -')
legend('能量算子包络','Hibert变换包络');

RMS=sqrt(sum(y.^2)/length(y))