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