时域特征提取

频域特征提取：https://blog.csdn.net/qq_37240982/article/details/107425375
时频域特征提取：https://blog.csdn.net/qq_37240982/article/details/107425452

时域、频域、时频域特征提取matlab程序的zip：
https://www.lanzoui.com/b01bp72xa

时域特征

时域信息是以时间为变量，描绘出信号的波形[22]。时域信号包括量纲特征参数以及无量纲特征参数。根据工作状况的差异，有量纲特征值的大小相应发生改变，而且工作环境对有量纲特征值有很大影响，具有表现不够稳定的缺陷，给工程应用带来一定困难，而无量纲指标对负载及转动速度的改变不敏感，能够更直观地表现出旋转设备运行中的状态信息。因此有量纲特征指标时常与无量纲指标一起使用。有量纲特征值容易理解、计算简单，故常被研究者使用。其中信号的平均值并不能反映样本数据的动态变化，通常用于检测，但其它参数的计算通常会用到平均值。工程应用中随着设备健康状态改变直至完全失效，均方根值、绝对均值、方根幅值都会有相应的增加。常用的有量纲时域特征参数集合，其表达式如下表 2.1 所示。

无量纲特征值中峰值指标反映了波形中峰值的尖峰程度。脉冲指标是信号峰值除以信号绝对值平均值；与峰值指标类似，脉冲指标也可以表征振动信号中含有瞬时尖峰与否。峭度指标能够很好地描述变量的分布。裕度指标通常可以表征部件的疲劳磨损程度。实际工程应用中当样本数据倍增，无量纲特征保持不变，但有量纲特征变大；在设备前期健康状态改变中，出现较多小幅度的尖峰，可以检测出均方根值变化不大，但由于裕度指标、峭度指标和脉冲指标对脉冲对健康状态改变的灵敏度高，以上三个指标会先上升再下降。结果表明：均方根值稳定好但对前期健康状态改变灵敏度低，裕度指
标、峭度指标和脉冲指标稳定性较均方根值差但更能表征设备前期健康状态改变的信息；表 2.2 为无量纲特征对故障的敏感性及稳定性的比较，常用的无量纲时域特征参数如表 2.3 所示。

时域特征提取matlab程序如下：

function [ timestruct ] = timeDomainFeatures( src)
% 计算时域统计特征
% src  N*1维矩阵 是要计算特征的源信号

if nargin>2
     error(message('参数有误，只有只能是一个或2个输入参数'));
end

%////////////////////////////////////////////////////////////////////////\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
timestruct.max=max(src);%1.最大值
timestruct.min=min(src);%2.最小值
timestruct.peak=max(abs(src)); %3.峰值
timestruct.p2p=max(src)-min(src);%4.峰峰值
timestruct.mean=mean(src);%5.均值
timestruct.averageAmplitude=mean(abs(src));%6.绝对平均值（平均幅值）
timestruct.rootAmplitude=mean(sqrt(abs(src)))^2;%7.方根幅值
timestruct.var=var(src,1);%8.方差  有偏
timestruct.std=std(src,1);%9.标准差
timestruct.rms=sqrt(sum(src.^2)/length(src));%10.有效值（均方根）
timestruct.kurtosis=kurtosis(src,1);%11.峭度
timestruct.skewness=skewness(src,1);%12.偏度
timestruct.shapeFactor=timestruct.rms/timestruct.averageAmplitude;%13.波形因子
timestruct.peakingFactor=timestruct.peak/timestruct.rms;%14.峰值因子（波峰因子）
timestruct.pulseFactor=timestruct.peak/timestruct.averageAmplitude;%15.脉冲因子
timestruct.marginFactor=timestruct.peak/timestruct.rootAmplitude;%16.裕度因子
timestruct.clearanceFactor=timestruct.peak/timestruct.rms^2;%17.余隙因子


%//////////////////////////////////////////////////////////////////////////\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
% if nargin==2
%      switch(option)
%          case 'max'
%              disp('max value=');
%              disp(timestruct.max);
%          case 'min'
%              disp('min value=');
%              disp(timestruct.min);           
%      end 
% end 

end

时域、频域、时频域特征提取matlab程序的zip：
https://www.lanzoui.com/b01bp72xa