Código MTALAB para calcular indicadores de características en el dominio del tiempo

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];

3. Representación del indicador

inserte la descripción de la imagen aquí

Supongo que te gusta

Origin blog.csdn.net/weixin_39458727/article/details/126191927
Recomendado
Clasificación