%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))