在之前的文章(Mr.括号:信号时域分析方法的理解(峰值因子、脉冲因子、裕度因子、峭度因子、波形因子和偏度等))里对时域指标做过一些分析。最近由于新建立一个公众号(括号的城堡)将会对以前讲过的内容进行系统的梳理。内容将在公众号首发,欢迎大家关注。

时域特征值是衡量信号特征的重要指标,时域特征值通常分为有量纲参数与无量纲参数。

所谓“量纲”,简单地理解就是“单位”。有量纲的参数就是有单位的,比如平均值,一段温度信号(单位℃)的平均值依旧是℃;无量纲的参数没有单位,无量纲量常写作两个有量纲量之积或比,但其最终的纲量互相消除后会得出无量纲量,比如,应变是量度形变的量,定义为长度差与原先长度之比。

有量纲的特征值往往具有直观的物理含义,是最为常用的特征指标。有量纲特征值主要包括:最大值、最小值、峰峰值、均值、方差、标准差、均方值、均方根值(RMS)、均方误差(MSE)、均方根误差(RMSE)、方根幅值等。

1.均值

均值、方差、均方值、均方根值之间有内在的联系。

均值是信号的平均,是一阶矩,可以表示为:

2.均方值

均方值是信号的平方的平均(信号→平方→平均值),代表了信号的能量,是二阶矩,可以表示为:

3.方差

方差是每个样本值与全体样本值的平均数之差的平方值的平均数,代表了信号能量的动态分量(均值的平方是静态分量),反应数据间的离散程度,是二阶中心距,可以表示为:

方差的不同表达方式,可以看出方差的几种理解方式:

(1)

式中可以看出:方差描述的是信号的离散程度,也就是变量离其期望值的距离。

(2)

式中可以看出:方差即平方的期望(均方值)减掉期望的平方。

(3)从物理含义上讲,均方值代表信号的能量,期望的平方代表信号的直流分量,而方差代表信号的交流分量。

4.标准差

标准差又叫均方差,是方差的算数平方根。标准差反应的是数据的离散程度

问题来了,方差和标准差都表示数据的离散程度,那么既然有了方差,为什么还要有标准差呢?

为了和原始信号统一量纲。

举个例子,假设北京一年的平均气温是20℃,气温标准差是10℃;乌鲁木齐一年的平均气温是20℃,气温标准差是15℃。这样会对气温的离散程度有一个直观理解,但如果说北京的气温方差是100,乌鲁木齐是225,就很不方便理解了。

5.均方根

均方根(RMS)又叫有效值。将所有值平方求和,求其均值,再开平方,就得到均方根值。或者说均方根值等于均方值的算数平方根。

其物理含义可以这样理解:让交流电与直流电分别通过同一电阻,若两者在相同的时间内所消耗的电能相等(或产生的焦耳热相同),那么该直流电的数值就叫做交流电的有效值。

6.均方误差

均方误差(MSE)是某种意义上的方差,均方误差是指参数估计值与参数真值之差平方的数学期望值。如果我们把随机变量的数学期望E认为是参数估计值(未来的),把随机变量本身作为参数真值,那么均方误差就是普通方差。

均方误差MSE可以评价数据的变化(偏离)程度,MSE的值越小(相互之间的比较,而不是跟参数真值的比较),说明预测模型描述实验数据具有更好的精确度。

均方误差在机器学习中常作为一种误差量度。

7.均方根误差

均方根误差(RMSE)就是均方误差的算术平方根:

均方误差与均方根误差,正如方差与标准差一样,是量纲上的区别,应用不同场合。

7. 特征提取的MATLAB代码实现

笔者最近在做这样一件事情。由于搞科研写论文时,对特征提取的需要往往是集中性的、多种类的、需求各异的,所以我把常用的特征值指标:

有量纲特征值8个——最大值、最小值、峰峰值、均值、方差、标准差、均方值、均方根值(RMS)

无量纲特征值6个——峭度、偏度、波形因子、峰值因子、脉冲因子、裕度因子

频域特征值5个——重心频率、均方频率、均方根频率、频率方差、频率标准差 (熵相关特征若干——这个后续补充。还会有其他常用、不常用的特征指标)

以上19种全都集中到一个封装函数里,实现一行代码完成特征提取。

比如提取数据“波形因子”就可以像这样写:

fea = genFeature(data,{'waveformF'}) %对data数据求波形因子

如果提取数据“最大值、最小值、峰峰值、均值、方差、标准差、均方值...”这全部19种特征,就可以这样写:

fea =genFeature(data,{'max','min','mean','peak','arv','var','std','kurtosis','skewness', 'rms','waveformF','peakF','impulseF','clearanceF'})%对data数据求波形因子

也就是说需要提取哪个特征,在函数中直接指定就可以了。输出的fea变量里就会得到相应的这些特征值,顺序也是与输入的排序保持一致的。

这个函数的介绍如下:

function fea = genFeature(data,fs,featureNamesCell)
% 时域信号特征提取
% 输入:
% data:待特征提取的时域信号
% fs:采样频率,如果不提取频域特征,fs值可以设置为1
% featureNamesCell:拟进行特征提取的特征名称,该变量为cell类型,其中包含的特征名称为字符串,特征名称需要在下边列表中:
% 目前支持的特征(2021.10.29):
% max :最大值
% min :最小值
% mean :平均值
% peak :峰峰值
% arv  :整流平均值
% var  :方差
% std  :标准差
% kurtosis  :峭度
% skewness  :偏度
% rms       :均方根
% waveformF :波形因子
% peakF     :峰值因子
% impulseF  :脉冲因子
% clearanceF:裕度因子
% FC:重心频率
% MSF:均方频率
% RMSF:均方根频率
% VF:频率方差
% RVF:频率标准差
%
% 输出:
% fea:数据data的特征值数组,其特征值顺序与featureNamesCell一一对应

需要上边这个函数文件以及测试代码的同学,可以在下述链接中获取:

时域、频域特征提取MATLAB代码

这个函数会持续更新,有哪些想要加进去的特征指标,同学们可以在评论区留言,笔者会考虑纳入到这个“特征提取指标全家桶”中。

总结

欢迎批评指正!

时域分析——有量纲特征值含义一网打尽相关推荐

  1. 时域特征偏度_时域分析——有量纲特征值含义一网打尽

    在之前的文章(Mr.括号:信号时域分析方法的理解(峰值因子.脉冲因子.裕度因子.峭度因子.波形因子和偏度等))里对时域指标做过一些分析.最近由于新建立一个公众号(括号的城堡)将会对以前讲过的内容进行系 ...

  2. 时域特征提取_时域分析——无量纲特征值含义一网打尽

    上一篇文章中我们讲到了有量纲的特征值含义.(Mr.括号:时域分析--有量纲特征值含义一网打尽)然而在信号表征时,有量纲指标虽然对信号特征比较敏感,但也会因工作条件(如负载)的变化而变化,并极易受环境干 ...

  3. 时域分析——无量纲特征值含义一网打尽

    上一篇文章中我们讲到了有量纲的特征值含义.(Mr.括号:时域分析--有量纲特征值含义一网打尽)然而在信号表征时,有量纲指标虽然对信号特征比较敏感,但也会因工作条件(如负载)的变化而变化,并极易受环境干 ...

  4. 概率论中的“矩”是什么意思

    在前两篇文章(时域分析--有量纲特征值含义一网打尽.时域分析--无量纲特征值含义一网打尽)中提到了"矩"这个概念.例如期望是一阶矩,方差是二阶矩等等.要怎么理解"矩&qu ...

  5. 三岁看大,七岁看老——基于退化模型进行剩余有效寿命预测的案例讲解

    本篇是寿命预测系列的第二个案例,在看这篇之前可以先看一下这里. 使用退化模型进行机电产品剩余寿命预测在工程应用中比较常见,退化模型法又有着几种不同的具体形式,博主见到比较多的包括基于随机过程的退化模型 ...

  6. 衡量数据“像不像”——协方差与相关系数

    协方差与相关系数也是在时域分析时常见的两个概念,他们都是用来描述数据"像不像"的. 1.协方差 对于信号X和Y,其协方差定义为: 可以描述为:(信号X减X期望)乘以(信号Y减Y期望 ...

  7. 信号时域分析方法的理解(峰值因子、脉冲因子、裕度因子、峭度因子、波形因子和偏度等)

    最近要做关于信号处理的一些工作,涉及时域.频域分析,滤波.神经网络以及贝叶斯等内容.在此做一些整理,作为梳理思路的笔记,也希望能对读者有所帮助. 文章如要转载请私信与我联系,并注明来源知乎专栏与信号处 ...

  8. 时域特征值提取的MATLAB代码实现(均方根、峰值因子、脉冲因子、裕度因子、峭度因子、波形因子和偏度等)

    时域计算代码相对简单,下列代码主要来自求助matlab时域分析程序-峰值.均方根.裕度.峭度_百度知道.不过原文中求"峭度因子"的方法不对,即"Kr = sum(y.^4 ...

  9. 特征值和特征向量的理解 浅显易懂 肯定有收获

    特征值和特征向量的理解 前言 一.矩阵是什么? 二.举个例子 1.计算特征值与特征向量 2.用特征向量表示任意向量 三.理解其他结论 1.对角化分解 2.矩阵的特征值分别是原矩阵特征值的倒数 3.特征 ...

  10. 如何理解特征值和特征向量

    学完线性代数的同学,可能会对线性代数的很多概念有所疑惑. 这个东西有什么用?那个玩意定义出来有什么意义? 本文将探讨线性代数中及其重要的两个概念:特征值与特征向量. (PS:下文中的矩阵AA均认为是方 ...

最新文章

  1. NYOJ 366 D的小L
  2. 独家 | 手把手教你推导决策树算法
  3. MySQL5.6多实例部署
  4. golang中的redigo
  5. 14.4 线程通讯-生产者与消费者
  6. POJ 2142——扩展欧几里得
  7. 面试题3:二维数组中的查找
  8. java web/Maven项目结构分析
  9. 8个前端技术学习网站
  10. Matlab 嵌套传递函数简化_MATLAB的数据处理方法及图形绘制详解
  11. emWin6.12模拟器发布,更新内容较多,增加环形控件,WIFI二维码(2020-04-09)
  12. 球缺体积和球冠表面积的计算公式及应用
  13. 老男孩教育运维班100台规模集群全网数据备份项目上机实战
  14. 99.扩展183.信息系统运维维护服务
  15. PDF怎么合并?教你2个免费一键合并PDF的方法
  16. 什么是挂载,Linux挂载如何实现详解
  17. 那些长期喝咖啡的人,现在都怎么样了?
  18. 虚拟环境下安装pytorch成功但import不成功
  19. Spring boot+Shiro+ spring MVC+swagger UI +Mybatis+mysql+Vue +Element UI 之二 vue 环境演示
  20. 取模运算性质_取模运算的性质

热门文章

  1. D3.js 绘制地图
  2. openlayers 地图上加图标_Openlayers绘制地图标注
  3. css固定定位如何基于父级来定为而不是视口
  4. 巅峰战舰服务器维护,《巅峰战舰》8月5日更新公告
  5. imitate wechat - 0
  6. EXSi虚拟机缺少vmdk文件报错问题
  7. raw格式镜像文件转vmdk
  8. JP1081B/9700_USB网卡驱动
  9. python 邮件抄送是什么意思_python使用SMTP发邮件时使用Cc(抄送)和Bcc(密送)...
  10. Ubuntu18.04-ROS-Python-通过周立功Can发送16进制整数