emd程序对机械振动信号分解

clc;
clear all
fs=15360;
s=load（‘d:/06.txt‘）;

N=length（s）;
n=1:N;
t=n/fs;
figure（1）
subplot（211）;
plot（ts）;
xlabel（‘t/s‘）;
ylabel（‘A/m^2‘）;
nfft1=length（s）;
Y2=abs（fft（snfft1））;
subplot（212）;
plot（（0:nfft1/2-1）/nfft1*fsY2（1:nfft1/2））;
xlabel（‘频率 f/Hz‘）;
ylabel（‘功率谱 P/W‘）;



%EMD分解
Nstd=0;
NE=1;
rslt=eemd（sNstdNE）;

x1=rslt（:1）;
 x2=rslt（:2）;
x3=rslt（:3）;
 x4=rslt（:4）;
 x5=rslt（:5）;
 x6=rslt（:6）;
 x7=rslt（:7）;
 x8=rslt（:8）;
 x9=rslt（:9）;
 figure（2）
 subplot（911）;
 plot（tx1）;
 subplot（912）;
 plot（tx2）;
 subplot（913）;
 plot（tx3）;
 subplot（914）;
 plot（tx4）;
subplot（915）;
plot（tx5）;
 subplot（916）;
 plot（tx6）;
 subplot（917）;
 plot（tx7）;
subplot（918）;
 plot（tx8）;
 subplot（919）;
 plot（tx9）;
 y1=x1;y2=x2;y3=x3;y4=x4