MTALAB计算时域特征指标的代码_matlab时域指标

1、参考文档

MATLAB官网的参考文档

2、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];
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

3、指标效果图