clc
clear?all
%?z=load（‘C:\Users\wangkun\Desktop\軸承診斷\test2.mat‘）;
%?x=z.X105_DE_time（1:1024）;
%?load?（‘IR007_0.mat‘）;?
load（‘F:\軸承數據\CaseWesternReserveUniversityData-master\12k_Drive_End_OR007@3_2_146.mat‘）
a=X146_DE_time;

x=a（10001:15000）;
N=5000;
fs=12000;
n=0:N-1;
f=n*fs/N;
lag=N;
n=0:N-1;
t=n/fs;??????
imf=emd（x）;
[mn]=size（imf）;?%imf為一m*n階矩陣，m是imf分量，n為數據點
emd_visu（x1:length（x）imfm）;?%實信號的信號重構及emd結果顯示函數
for?i=1:m
??a1（i）=kurtosis（imf（i:））;%峭度
??b1（i）=mean（imf（i:））;?%均值;
??c1（i）=var（imf（i:））;??%方差;
??d1（i）=std（imf（i:））;??%均方值
??e1（i）=std（imf（i:））.^0.5;?%均方根值
??f1（i）=s