matlab求解时域特征参数,matlab 实现数据时域特征提取

matlab 实现数据时域特征提取

2020-07-31 02:15

阅读数 64

%数据时域特征提取 clc; close ; clear all; %数据导入 data=xlsread('E:\排球\原始数据 z方向.xlsx',1)';

%初始参数设置 [m,n]=size(data); D=[]; DA=[]; %求解特征 for i=1:1:m d=data(i,:)

d=d(~isnan(d));%去除NAN值 ave=mean(d);%均值 u=std(d);%标准差 time=length(d);%时间

theta=var(d);%方差 area=sum(abs(d));%面积 maxv=max(d);%最大值 minv=min(d);%最小值

[dd,minp,maxp]=premnmx(d); %求解信息熵 entropy=yyshang(dd,9);%信息熵

D=[D;ave;maxv;minv;u;area;time;theta;entropy]; DA=[DA,D]; D=[]; end DA=DA';

其中yyshang是求解熵的函数，代码如下

function Hx=yyshang(y,duan) %不以原信号为参考的时间域的信号熵 %输入：maxf:原信号的能量谱中能量最大的点

%y:待求信息熵的序列 %duan:待求信息熵的序列要被分块的块数 %Hx:y的信息熵

%duan=10;%将序列按duan数等分，如果duan=10,就将序列分为10等份 x_min=min(y); x_max=max(y);

maxf(1)=abs(x_max-x_min); maxf(2)=x_min; duan_t=1.0/duan;

jiange=maxf(1)*duan_t; % for i=1:10 %

pnum(i)=length(find((y_p>=(i-1)*jiange)&(y_p

pnum(1)=length(find(y

pnum(i)=length(find((y>=maxf(2)+(i-1)*jiange)&(y

pnum(duan)=length(find(y>=maxf(2)+(duan-1)*jiange)); %sum(pnum)

ppnum=pnum/sum(pnum);%每段出现的概率 %sum(ppnum) Hx=0; for i=1:duan if ppnum(i)==0

Hi=0; else Hi=-ppnum(i)*log2(ppnum(i)); end Hx=Hx+Hi; end end