CN110108467A - 基于便携式移动设备的主动发声测速方法 - Google Patents
基于便携式移动设备的主动发声测速方法 Download PDFInfo
- Publication number
- CN110108467A CN110108467A CN201910416881.8A CN201910416881A CN110108467A CN 110108467 A CN110108467 A CN 110108467A CN 201910416881 A CN201910416881 A CN 201910416881A CN 110108467 A CN110108467 A CN 110108467A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- envelope
- harmonic
- spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
基于便携式移动设备的主动发声测速方法,将移动设备的发声单元对准被测试轴后,通过发声单元发出设定频率的声音,同时麦克风记录接收到的声音信号。使用有限脉冲响应滤波器对原始信号进行带通滤波,然后通过希尔伯特变换得到信号的包络,再利用快速傅里叶变换得到包络信号的频谱。最后通过谐波乘积谱法在频域累积包络信号中各阶谐波的能量,搜索谐波乘积谱中的最大幅值对应频率,准确确定转速对应的基频分量。本发明方法利用便携式移动设备的发声单元、麦克风硬件以及计算显示能力,无须额外的专用测试硬件,通过包络解调和谐波乘积谱的信号处理方法,对轴的瞬时转速进行准确估计,具有良好的鲁棒性和抗干扰能力。
Description
技术领域
本发明涉及机械旋转轴的转速测量技术和信号处理技术领域,特别涉及基于便携式移动设备的主动发声测速方法。
背景技术
机械系统中轴的瞬时转速是反映机械设备运行状态的重要参数,同时也是变转速工况下进行机械故障诊断的必要信息,机械系统中轴瞬时转速的准确获取对机械设备的健康状况监测与故障诊断具有重要意义。
轴的转速测量一般通过专门的测速硬件实现,如:编码器、光电测速传感器,把这些硬件安装在靠近旋转部件或与部件直接耦合的位置,再把所测信号转换成角位移或者转速显示,但采用硬件测速方法的测试成本高,还存在无法安装测速装置的现场情况,给测速技术的普及应用带来了困难。
专利申请号为201810414579.4的专利提出从安装在机械设备表面的振动信号还原轴的瞬时转速信息,进而用于阶次跟踪和故障特征提取。但由于机械振动信号结构十分复杂,从振动信号中获取瞬时转速信息涉及冗繁的信号处理过程,容易受到干扰振动分量的影响,存在鲁棒性不高的问题。此外,振动信号的获取也依赖专门的测试软硬件。
发明内容
本发明的目的在于使技术人员可以使用便携式移动设备(智能手机、平板电脑、笔记本电脑),通过提出的发明方法对机械设备中轴的瞬时转速进行精确测量,而无须额外的专用测试硬件。利用便携式移动设备的发声单元、麦克风硬件以及计算显示能力,通过包络解调和谐波乘积谱的信号处理方法,对轴的瞬时转速进行准确估计,具有良好的鲁棒性和抗干扰能力。
为了达到上述目的,本发明采取的技术方案为:
基于便携式移动设备的主动发声测速方法,包括以下步骤:
步骤一:将移动设备的发声单元对准被测试轴,通过发声单元发出设定频率的声音,同时麦克风记录接收到的声音信号,记为x(n)。
步骤二:构建有限脉冲响应(finite impulse response,FIR)带通滤波器,对信号x(n)进行带通滤波,得到滤波信号y(n)。
FIR带通滤波器中心频率设置为发声频率,频带宽度设为被测轴极限转频的8倍。
步骤三:由滤波信号y(n)通过希尔伯特(Hilbert)变换得到信号y(n)的包络z0(n),再作去均值处理和去趋势项处理得到z(n)。
利用希尔伯特变换得到信号的包络的原理是让原信号产生一个90°的相移动,从而与原信号构成一个解析信号,将解析信号的幅值作为原信号的包络。
去趋势项的作用是减少滤波信号中由于移动设备移动和外部测试环境引起的低频分量的干扰。
步骤四:对包络信号z(n)作快速傅里叶变换(FFT),得到测速信号不同频率下对应的幅值,记为Zn(ejw)。
步骤五:通过谐波乘积谱法(Harmonic Product Spectrum,HPS)将Zn(ejw)中多阶谐波的能量进行累积,搜索谐波乘积谱中的最大幅值对应频率,即为被测轴转频f。
信号z(n)的谐波乘积谱定义为:其中,Zn(ejrw)是频谱 Zn(ejw)在频域压缩r倍的结果,使Zn(ejw)的r阶谐波峰值位置总与Zn(ejw)的基频频率对齐,从而累积各次谐波的能量。
本发明相比于现有技术,具有以下优势:
a)基于便携式移动设备的主动发声测速方法,相比传统硬件测速方法成本更低,具有便捷易用的优势。
b)本发明采用的主动发声方法相比基于振动信号的测速方法,检波信号信噪比高,信号处理过程无需先验知识,有利于实现瞬时转速的实时测量与显示。
c)本发明提出的方法能够克服声音信号中的倍频分量干扰与无关频率分量干扰,准确确定转速对应的基频分量,具有良好鲁棒性。
附图说明
图1为本发明实施例试验示意图。
图2为本发明方法的流程图。
图3为本发明实施例的原始声音信号x(n)。
图4为本发明实施例的原始信号的频谱图。
图5为本发明方法构建的FIR带通滤波器的频率响应。
图6为本发明实施例的带通滤波后的信号y(n)。
图7为本发明实施例的带通滤波后信号的频谱图。
图8为本发明实施例的滤波信号的包络信号z(n)。
图9为本发明实施例的滤波信号的包络谱。
图10为本发明实施例的谐波乘积谱图Pn(ejw)。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
以德帆(DEFAN)YS7124型号三相交流电动机与魅族15组成实施例试验台,电动机转轴与移动设备的相对位置、移动设备发出的载波信号频率和载波幅值在测试过程中保持不变,试验示意图如图1所示。
具体参数如下:1)电机转速为900RPM(15Hz);2)移动设备为智能手机(魅族 15);3)载波频率为5kHz;4)测速移动设备与转轴轴心距离为5cm;5)测试装备对声音信号进行高频采样和数据处理,采样过程的频率为44100Hz,采样时间为1.47s。
应用本发明对原始数据分析确定轴转速,如图2所示,基于便携式移动设备的主动发声测速方法,包括以下步骤:
步骤一:设置载波频率为5kHz,移动设备的发声单元对准被测试轴,距离5cm,麦克风记录接收到的声音信号记为x(n)。x(n)通过快速傅里叶(FFT)变换得到频谱,如图4 所示,在载波频率附近轴转速信息被掩盖,且整个频带上有许多杂波干扰,难以确定轴转速。
步骤二:构建有限脉冲响应(finite impulse response,FIR)带通滤波器。FIR带通滤波器的中心频率设置为载频频率,带宽设为轴极限转频的8倍。本发明实施例中的FIR滤波器的频率响应特性如图5所示。信号x(n)滤波后的信号y(n)如图6所示。
步骤三:通过希尔伯特(Hilbert)变换得到信号y(n)的包络z0(n),并作去均值处理和去趋势项处理得到z(n)。作为示例,由步骤二中的滤波信号y(n)得到的包络信号z(n)如图8所示。
利用希尔伯特变换得到信号的包络的原理是让原信号产生一个90°的相移动,从而与原信号构成一个解析信号,将解析信号的幅值作为原信号的包络。
去趋势项的作用是减少滤波信号中由于移动设备移动和外部测试环境引起的低频分量的干扰。
步骤四:对包络信号z(n)作快速傅里叶变换(FFT),得到测速信号不同频率下对应的幅值,记为Zn(ejw)。实施例中,包络谱如图9所示,如果直接通过包络谱幅值极值对应频率来估计转频,将受到高阶谐波(30Hz和45Hz分量)和无关杂波(0.9Hz和48Hz分量) 的干扰,而产生错误的转频估计结果。
步骤五:利用谐波乘积谱法(Harmonic Product Spectrum,HPS)将Zn(ejw)中多个谐波的能量累积起来,再搜索谐波乘积谱中的最大幅值对应频率,即为被测轴转频f。实施例中,谐波乘积谱的频域压缩倍数r设为3,结果如图10所示,谱图中幅值最大位置准确对应了设电机实际转频(15Hz),有效克服了倍频和杂波干扰的影响。
本发明提出的基于便携式移动设备的主动发声测速方法,克服了信号中转频的高阶谐波和其他杂波的干扰,可以对轴的转速进行准确估计,具有良好的鲁棒性。
Claims (1)
1.基于便携式移动设备的主动发声测速方法,其特征在于,包括以下步骤:步骤一:将移动设备的发声单元对准被测试轴,通过发声单元发出设定频率的声音,同时麦克风记录接收到的声音信号,记为x(n)。步骤二:构建有限脉冲响应(finite impulse response,FIR)带通滤波器,对信号x(n)进行带通滤波,得到滤波信号y(n)。FIR带通滤波器中心频率设置为发声频率,频带宽度设为被测轴极限转频的8倍。步骤三:由滤波信号y(n)通过希尔伯特(Hilbert)变换得到信号y(n)的包络z0(n),再作去均值处理和去趋势项处理得到z(n)。利用希尔伯特变换得到信号的包络的原理是让原信号产生一个90°的相移动,从而与原信号构成一个解析信号,将解析信号的幅值作为原信号的包络。去趋势项的作用是减少滤波信号中由于移动设备移动和外部测试环境引起的低频分量的干扰。步骤四:对包络信号z(n)作快速傅里叶变换(FFT),得到测速信号不同频率下对应的幅值,记为Zn(ejw)。步骤五:通过谐波乘积谱法(Harmonic Product Spectrum,HPS)将Zn(ejw)中多阶谐波的能量进行累积,搜索谐波乘积谱中的最大幅值对应频率,即为被测轴转频f。信号z(n)的谐波乘积谱定义为:其中,Zn(ejrw)是频谱Zn(ejw)在频域压缩r倍的结果,使Zn(ejw)的r阶谐波峰值位置总与Zn(ejw)的基频频率对齐,从而累积各次谐波的能量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910416881.8A CN110108467B (zh) | 2019-05-20 | 2019-05-20 | 基于便携式移动设备的主动发声测速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910416881.8A CN110108467B (zh) | 2019-05-20 | 2019-05-20 | 基于便携式移动设备的主动发声测速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110108467A true CN110108467A (zh) | 2019-08-09 |
CN110108467B CN110108467B (zh) | 2020-10-20 |
Family
ID=67490961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910416881.8A Active CN110108467B (zh) | 2019-05-20 | 2019-05-20 | 基于便携式移动设备的主动发声测速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110108467B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110514865A (zh) * | 2019-09-12 | 2019-11-29 | 核工业理化工程研究院 | 一种多路同步测量旋转机械设备转速及振幅的便携式测试仪 |
CN111863028A (zh) * | 2020-07-20 | 2020-10-30 | 江门职业技术学院 | 一种发动机声音合成方法及系统 |
CN113702046A (zh) * | 2021-09-13 | 2021-11-26 | 长沙理工大学 | 基于移动设备的变转速工况下轴承故障诊断方法 |
CN115586345A (zh) * | 2022-12-09 | 2023-01-10 | 深圳市双合电气股份有限公司 | 一种基于定子电流的电机自适应转速估计方法及系统 |
DE102022201317A1 (de) | 2022-02-09 | 2023-08-10 | Robert Bosch Gesellschaft mit beschränkter Haftung | Verfahren und Vorrichtung zum Erkennen einer Drehgeschwindigkeit einer Maschine |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1339166A (zh) * | 1999-01-22 | 2002-03-06 | 魁北克水电公司 | 高压机电开关系统中的振动声学特征的处理方法 |
CN1926413A (zh) * | 2004-10-18 | 2007-03-07 | 日本精工株式会社 | 机械设备的异常诊断系统 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN104464722A (zh) * | 2014-11-13 | 2015-03-25 | 北京云知声信息技术有限公司 | 基于时域和频域的语音活性检测方法和设备 |
CN104819766A (zh) * | 2015-05-13 | 2015-08-05 | 西安交通大学 | 基于谐噪比的包络解调频带确定方法 |
CN204575672U (zh) * | 2015-03-05 | 2015-08-19 | 中国人民解放军63956部队 | 一种非接触式发动机测速仪 |
CN105203914A (zh) * | 2015-07-22 | 2015-12-30 | 广东电网有限责任公司电力科学研究院 | 一种突发短路下的变压器绕组状态诊断方法 |
CN106053871A (zh) * | 2016-07-25 | 2016-10-26 | 昆明理工大学 | 一种通过滚动轴承滚道剥落故障对应冲击特征提取转速的方法 |
CN108871742A (zh) * | 2018-05-03 | 2018-11-23 | 西安交通大学 | 一种改进的无键相故障特征阶次提取方法 |
EP3239667B1 (en) * | 2012-01-30 | 2018-11-28 | S.P.M. Instrument AB | An analysis system |
-
2019
- 2019-05-20 CN CN201910416881.8A patent/CN110108467B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1339166A (zh) * | 1999-01-22 | 2002-03-06 | 魁北克水电公司 | 高压机电开关系统中的振动声学特征的处理方法 |
CN1926413A (zh) * | 2004-10-18 | 2007-03-07 | 日本精工株式会社 | 机械设备的异常诊断系统 |
EP3239667B1 (en) * | 2012-01-30 | 2018-11-28 | S.P.M. Instrument AB | An analysis system |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN104464722A (zh) * | 2014-11-13 | 2015-03-25 | 北京云知声信息技术有限公司 | 基于时域和频域的语音活性检测方法和设备 |
CN204575672U (zh) * | 2015-03-05 | 2015-08-19 | 中国人民解放军63956部队 | 一种非接触式发动机测速仪 |
CN104819766A (zh) * | 2015-05-13 | 2015-08-05 | 西安交通大学 | 基于谐噪比的包络解调频带确定方法 |
CN105203914A (zh) * | 2015-07-22 | 2015-12-30 | 广东电网有限责任公司电力科学研究院 | 一种突发短路下的变压器绕组状态诊断方法 |
CN106053871A (zh) * | 2016-07-25 | 2016-10-26 | 昆明理工大学 | 一种通过滚动轴承滚道剥落故障对应冲击特征提取转速的方法 |
CN108871742A (zh) * | 2018-05-03 | 2018-11-23 | 西安交通大学 | 一种改进的无键相故障特征阶次提取方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110514865A (zh) * | 2019-09-12 | 2019-11-29 | 核工业理化工程研究院 | 一种多路同步测量旋转机械设备转速及振幅的便携式测试仪 |
CN111863028A (zh) * | 2020-07-20 | 2020-10-30 | 江门职业技术学院 | 一种发动机声音合成方法及系统 |
CN111863028B (zh) * | 2020-07-20 | 2023-05-09 | 江门职业技术学院 | 一种发动机声音合成方法及系统 |
CN113702046A (zh) * | 2021-09-13 | 2021-11-26 | 长沙理工大学 | 基于移动设备的变转速工况下轴承故障诊断方法 |
CN113702046B (zh) * | 2021-09-13 | 2024-06-11 | 长沙理工大学 | 基于移动设备的变转速工况下轴承故障诊断方法 |
DE102022201317A1 (de) | 2022-02-09 | 2023-08-10 | Robert Bosch Gesellschaft mit beschränkter Haftung | Verfahren und Vorrichtung zum Erkennen einer Drehgeschwindigkeit einer Maschine |
CN115586345A (zh) * | 2022-12-09 | 2023-01-10 | 深圳市双合电气股份有限公司 | 一种基于定子电流的电机自适应转速估计方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110108467B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110108467A (zh) | 基于便携式移动设备的主动发声测速方法 | |
McCormick et al. | Cyclostationarity in rotating machine vibrations | |
Wang et al. | Rolling element bearing fault diagnosis via fault characteristic order (FCO) analysis | |
CN112101245B (zh) | 基于频域窗函数的短时傅里叶变换机械冲击特征提取方法 | |
CN110617964A (zh) | 用于滚动轴承故障诊断的同步压缩变换阶比分析法 | |
CN111665051A (zh) | 基于能量权重法的强噪声变转速条件下轴承故障诊断方法 | |
CN110907162B (zh) | 一种变转速下无转速计的旋转机械故障特征提取方法 | |
CN109101768A (zh) | 基于压缩感知的叶端定时传感器布局优化设计方法 | |
EP2237002A1 (en) | Timing analysis | |
CN105466550B (zh) | 非均匀欠采样叶端定时振动信号重构方法及其装置 | |
CN108871742B (zh) | 一种改进的无键相故障特征阶次提取方法 | |
CN102353500B (zh) | 一种用于动平衡测量的不平衡信号提取方法 | |
CN107843740B (zh) | 一种融合振动与声音信号频谱特征的转速测量方法 | |
CN108761117B (zh) | 便携式电流检测转速测试仪 | |
CN109101769A (zh) | 一种基于压缩感知的叶端定时传感器数目确定方法 | |
CN109520738A (zh) | 基于阶次谱和包络谱的旋转机械滚动轴承故障诊断方法 | |
CN103116032A (zh) | 用于获取风力发电机组的转速的方法和设备 | |
CN106872017A (zh) | 一种磁电式低频振动传感器的幅相补偿方法 | |
CN114061678A (zh) | 一种科氏流量计数字驱动方法 | |
Shi et al. | A dual-guided adaptive decomposition method of fault information and fault sensitivity for multi-component fault diagnosis under varying speeds | |
CN101368870B (zh) | 一种用于机械转子单截面轴振分析的幅值频谱的制作方法 | |
CN112782421B (zh) | 一种基于音频的转速识别方法 | |
CN112465068A (zh) | 一种基于多传感器数据融合的旋转设备故障特征提取方法 | |
CN107247155A (zh) | 一种汽轮发电机组转子瞬时转速测试装置及其测试方法 | |
CN103913271B (zh) | 非平稳转速时转子动不平衡信号的提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |