专注于工业智能预警系统研发, 通过机理算法和数据驱动算法分析振动信号、音频、DCS、PLC信号、SCADA信号等设备运行状态数据对机器设备进行看病预诊,为机器设备健康运行保驾护航。 网站正在不断建设和完善过程中,欢迎大家给予建议和参与社区建设
52phm,专注于预测性维护知识学习和交流,欢迎广大从事预测性维护行业人员投稿,投稿请联系管理员(wx: www52phmcn),投稿内容可以是:
振动信号频域分析首先需要把信号的时域波形借助离散傅里叶变换转化为频谱信息,公式如下:
式中:
x(kΔt)为振动信号的采样值;N为采样点数;Δt为采样间隔;k为时域离散值的序号。
求得频谱信息后,可根据频域统计指标公式计算相应的值,公式如下:
def get_fre_domain_features(f,y):
fre_line_num = len(y)
p1 = y.mean()
p2 = math.sqrt(sum((y-p1)**2)/fre_line_num)
p3 = sum((y-p1)**3)/(fre_line_num*p2**3)
p4 = sum((y-p1)**4)/(fre_line_num*p2**4)
p5 = sum(f*y)/sum(y)
p6 = math.sqrt(sum((f-p5)**2*y)/fre_line_num)
p7 = math.sqrt(sum(f**2*y)/sum(y))
p8 = math.sqrt(sum(f**4*y)/sum(f**2*y))
p9 = sum(f**2*y)/math.sqrt(sum(y)*sum(f**4*y))
p10 = p6/p5
p11 = sum((f-p5)**3*y)/(p6**3*fre_line_num)
p12 = sum((f-p5)**4*y)/(p6**4*fre_line_num)
p13 = sum(abs(f-p5)*y)/(math.sqrt(p6)*fre_line_num)
p = [p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13]
return p
补充:---------------------------------------------------------------------
程序中的f是傅里叶变换之后的频率轴,y是傅里叶变换之后的幅值,也就是平时画的频谱图的x轴和y轴。
def nextpow2(x):
if x == 0:
return 0
else:
return int(np.ceil(np.log2(x)))
def Do_fft(sig,Fs):#输入信号和采样频率
xlen = len(sig)
sig = sig - sig.mean()
NFFT = 2**nextpow2(xlen)
yf = np.fft.fft(sig,NFFT)/xlen*2
yf = abs(yf[0:int(NFFT/21)])
f = Fs/2*np.linspace(0,1,int(NFFT/21))
f = f[:]
return f,yf
#频域离散值的序号
f,y = Do_fft(sig, Fs)
p = get_fre_domain_features(f,y)
版权声明:遵循 CC 4.0 BY-SA 版权协议原文链接:https://blog.csdn.net/baidu_38963740/article/details/111824612
2021-12-07 23:41:14
互联网
950
分类:算法开发
专栏:振动信号预处理
2021-12-10 17:18:58
互联网
902
分类:算法开发
专栏:故障诊断
2021-12-14 18:49:46
互联网
1376
分类:算法开发
专栏:振动信号预处理
关注公众号进群
让志同道合读者学习交流
评价一个故障诊断系统的性能指标有: 1)故障检测的及时性:是指系统在发生故障后,故障诊断系统在最短时间内检测到故障的能力。故障发生到被检测出的时间越短说明故障检测的及时性越好。 2)早期检测的灵敏度:是指故障诊断系统对微小故障信号的检测能力。故障诊断系统能检测到的故障信号越小说明其早期检测的灵敏度越高。 3)故障的误报率和漏报率:误报指系统没有出去故障却被错误检测出发生故障
2021-12-13 13:51:37
互联网
1961
分类:算法开发
专栏:故障诊断
1. 故障诊断概念故障诊断主要研究如何对系统中出现的故障进行检测、分离和辨识 , 即判断故障是否发生 , 定位故障发生的部位和种类 , 以及确定故障的大小和发生的时间等 。2. 故障诊断方法故障诊断防范可分为定性分析和定量分析两大类 , 如图 1 所示。 其中 , 定量分析方法又分为基于解析模型的方法和数据驱动的方法 , 后者又进一步包括机器学习类方法、多元统计分析类方法、信号处理类方...
2021-12-13 14:48:01
互联网
2029
分类:算法开发
专栏:故障诊断
DSP期末总复习基于教材:《数字信号处理——原理、实现及应用(第3版)》 高西全等 编著一、基本概念1.1 信号1.2 系统1.3 卷积1.4 差分方程二、傅里叶变换2.1 连续信号:CTFT和CTFS2.2 时域离散信号:DTFT和DTFS2.3 周期信号的傅里叶变换2.4 离散傅里叶变换(DFT)2.5 快速傅里叶变换(FFT)三、Z变换3.1 Z变换3.2 Z变换的性质3.3 逆Z变换3.4 Z变换的应用四、模拟信号数字处理五、滤波器5.1 一些散落的实例5.
2021-12-14 10:44:43
互联网
795
分类:算法开发
专栏:数字信号处理
学习——信号调制识别 (一) 看了《通信信号调制识别技术及其发展》这一论文后,将学习到的知识记录在这篇博客里。1、通信信号调制识别技术 调制识别问题从本质上来说是一种典型的模式识别问题。其基本框架如图所示: 调制识别由三部分组成:信号预处理、提取特征参数和分类识别。信号预处理部分包括载波同步、频率下变频、噪声抑制以及对信噪比、符号周期、载波频率等参数的估计。特征提取部分...
2021-12-14 22:00:52
互联网
872
分类:算法开发
专栏:数字信号处理
光谱原始数据通常存在伪影,噪声等问题,而预处理步骤可能对准确率造成影响,同时对不同数据的适应性差,神经网络能提供端到端的方法(或者省去某个预处理步骤),因而具有一定前景。数据驱动的深度学习技术可以发现大数据集中复杂的结构,并从数据中提取关键特征。一.综述文章:Deep learning for vibrational spectral analysis: Recent progress and a practical guide(2019)自2017年提出,将深度学习应用于光谱处理这一领域,迄
2021-12-14 22:12:02
互联网
740
分类:论文速递
专栏:未分组
机械振动(按频率分类):低频振动(小于10HZ)、中频振动(10HZ-1000HZ)、高频振动(大于1000)HZ 软件积分误差中趋势项产生的原因:所谓趋势项是在测试信号中存在线性项或缓变的非线性项成分。工程实际测量的信号大部分是复杂周期信号与随机信号 的 混 合,而 且 周 期 信 号 往 往 是 研 究对象。趋势项的存在会使数值积分的结果产生很大的误差,严重地背离
2021-12-17 09:40:29
互联网
718
分类:算法开发
专栏:振动信号预处理
共振频率与固有频率是不是同一个?从数值上来说,它们是相等的。但是两个概念是不同的。当一个装置成型时,他本身发生的振动的频率是固定的,这一频率就是固有频率。比如一个单摆做好后,他的振动频率等于2*Pi*(l/g)^(1/2),l是单摆的长度,g是重力加速度,所以这个单摆的振幅无论多大,加在下面的东西多重,只要是没有外界的干扰,都以一个频率振动(固有频率)。而当我们用一个周期的力推这个单摆
2021-12-17 09:41:37
互联网
837
分类:算法开发
专栏:振动信号预处理
概述语音信号是一种非平稳的时变信号,它携带着大量信息。在语音编码、语音合成、语音识别和语音增强等语音处理中,都需要提取语音中包含的各种信息语音处理的目的对语音信号进行分析,提取特征参数,用于后续处理加工语音信息,如语音增强和语音合成中的应用根据所分析的参数类型,语音信号可以分成:时域分析最简单、最直观直接对语音信号的时域波形进行分析特征参数:语音的短时能力、平均幅度、短时...
2022-01-13 16:53:53
互联网
790
分类:算法开发
专栏:语音信号预处理
数据集和代码均已上传到Github中,欢迎大家下载使用。Github地址:https://github.com/JasonZhang156/Sound-Recognition-Tutorial如果这个教程对您有所帮助,请不吝贡献您的小星星Q^Q.数据分析本节针对ESC-10数据集进行基本的数据分析,包括数据样本数,数据类别数,每类声音样本数等信息。并且对每类样本的声音波形,功率谱进...
2022-03-23 22:57:45
互联网
647
分类:算法开发
专栏:语音信号预处理
故障诊断之基于振动信号的阶比谱分析
2022-05-31 11:08:40
互联网
1569
分类:算法开发
专栏:振动信号预处理
从事设备故障预测与健康管理行业多年的PHM算法工程师(机器医生)、国际振动分析师, 实践、研发和交付的项目涉及“化工、工业机器人、风电机组、钢铁、核电、机床、机器视觉”等领域。专注于工业智能预警系统研发, 通过机理算法和数据驱动算法分析振动信号、音频、DCS、PLC信号、SCADA信号等设备运行状态数据对机器设备进行看病预诊,为机器设备健康运行保驾护航。