Código MTALAB para calcular indicadores de características en el dominio del tiempo
1. Documentos de referencia
Documentos de referencia en el sitio web oficial de MATLAB
2. código MATLAB
function valfactor=time_statistical_compute(x)
%对时域信号进行统计量分析
%val返回有量纲指标,factor返回无量纲指标
N=length(x);
p1=mean(x); %均值
xx=x-p1;
p11=sqrt(sum(xx.^2)/(N-1));%标准差
p2=rms(x); %均方根值 p2=sqrt(sum(x.^2)/N);
p3=(sum(sqrt(abs(x)))/N).^2; %方根幅值
p4=sum(abs(x))/N; %绝对平均值
p5=sum(xx.^3)/((N-1)*p11^3); %偏斜度
p6=kurtosis(x); %峭度 p6=sum(x.^4)/N;
p7=sum((xx).^2)/N; %方差
p8=max(x);%最大值
p9=min(x);%最小值
p10=p8-p9;%峰峰值
%%以上都是有量纲统计量,以下是无量纲统计量
f1=p2/p4; %波形指标
% f1=p2/(sum(x)/N);
f2=p8/p2; %峰值指标 E[MAX(X)]=P8
f3=p8/p4; %脉冲指标
f4=p8/p3; %裕度指标
f5=p5/((N-1)*(p11)^3); %偏斜度指标
% f5=p5/(p2^3); %偏斜度指标
f6=p6/((N-1)*(p11)^4); %峭度指标
% f6=p6/(p2^4); %峭度指标
%valfactor=[val;factor];
valfactor=[p1;p2;p3;p5;p6;p8;p11;f1;f2;f3;f4];