CN108921082A - 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法 - Google Patents

基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法 Download PDF

Info

Publication number
CN108921082A
CN108921082A CN201810682880.3A CN201810682880A CN108921082A CN 108921082 A CN108921082 A CN 108921082A CN 201810682880 A CN201810682880 A CN 201810682880A CN 108921082 A CN108921082 A CN 108921082A
Authority
CN
China
Prior art keywords
frequency
combustion
cylinder
hilbert transform
resonant frequency
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
Application number
CN201810682880.3A
Other languages
English (en)
Other versions
CN108921082B (zh
Inventor
杜敬涛
陈曦明
刘龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201810682880.3A priority Critical patent/CN108921082B/zh
Publication of CN108921082A publication Critical patent/CN108921082A/zh
Application granted granted Critical
Publication of CN108921082B publication Critical patent/CN108921082B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising
    • G06F2218/06Denoising by applying a scale-space analysis, e.g. using wavelet analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明提供的是一种基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法。一:对缸压信号进行傅里叶变换获得缸压频谱,利用经验小波算法中的尺度空间算法识别缸压频谱中的特征子模态,并确定子模态的位置,计算子模态的边界频率;二:构造经验小波滤波器,利用构造的小波滤波器对缸压信号进行带通滤波,提取缸压频谱中对应每个燃烧波动的子模态信号,选取振荡能量最大的子模态为燃烧共振信号;三:将步骤二所得到的燃烧共振信号进行希尔伯特变换,并绘制希尔伯特谱,获取共振频率及其对应的共振振幅。本发明具有收敛性好、自适应性强、精确度高的优点,获得的共振频率能够用于柴油机低噪声优化设计和缸内温度与气体质量估计。

Description

基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取 方法
技术领域
本发明涉及的是一种信号处理方法,具体地说是一种柴油机燃烧共振频率的提取方法。
背景技术
柴油机的燃烧共振频率是柴油机低噪声优化与热力分析的重要参数。在燃烧室设计时,可以通过优化燃烧室结构使其固有频率远离燃烧共振频率,从而能够大幅降低柴油机的燃烧噪声。此外,燃烧共振频率还与缸内温度有关,并且可以计算缸内气体质量,因此精确计算共振频率对柴油机的相关研究至关重要。
目前使用的计算柴油机燃烧共振频率的方法主要以信号处理方法为主。BodiscoT在论文《贝叶斯模型确定直喷柴油机共振频率》中利用贝叶斯统计模型对柴油机燃烧共振频率进行预测,但是该方法的局限性在于若统计模型与缸压信号不匹配,则模型计算不收敛。此外,Broatch A在论文《有限元法确定杯状活塞发动机的共振频率—用于估计缸内气体质量》中利用短时傅里叶变换的时频分析方法来获得共振频率。但是短时傅里叶变换由于需要对窗函数的尺寸和类型进行合适地选取,所以共振频率的计算精度也会受此影响。因此如何提供一种即不需要选取经验函数,收敛性又好的燃烧共振频率计算方法,成为本领域研究人员需要考虑的问题。
发明内容
本发明的目的在于提供一种不需要选取经验函数,收敛性好的基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法。
本发明的目的是这样实现的:
步骤一:对缸压信号进行傅里叶变换获得缸压频谱,利用经验小波算法中的尺度空间算法识别缸压频谱中的特征子模态,并确定子模态的位置,计算子模态的边界频率,将整个频谱中对应每个燃烧波动的子模态边界识别出来;
步骤二:根据Meyer小波构造方法来构造经验小波滤波器,利用构造的小波滤波器对缸压信号进行带通滤波,提取缸压频谱中对应每个燃烧波动的子模态信号,在大于4000Hz的子模态信号中选取振荡能量最大的子模态,该信号即为燃烧共振信号;
步骤三:将步骤二所得到的燃烧共振信号进行希尔伯特变换,并绘制希尔伯特谱,获取共振频率及其对应的共振振幅。
本发明还可以包括:
1.步骤一具体包括:
首先对缸压信号f(t)进行傅里叶变换其中t为时间,f为频率。利用尺度空间算法确定缸压频谱中N个有意义的极大值点,并将其作为燃烧波动模态。具体为先定义核函数其中n为尺度参数,则离散尺度空间定义为其中由于傅里叶谱混有噪声,因此尺度空间算法移除特征长度小于的局部细节信息,从而获得频谱中粗略的模态信息,根据尺度空间算法所定义的模态信息,设ω为每个模态的边界,并且ω={ωi}i=1,2,…,N,i为子模态编号,ω0=0,ωN=π,此外,以ωi为中心定义宽度为2τi的频率传递带,其中τi=γωi,并且0<γ<1。
2.步骤二具体包括:
根据Meyer小波构造方法来构造经验小波滤波器,获得尺度函数为和小波函数将构造的经验小波滤波器与缸压信号做卷积,即自适应分割柴油机缸压频谱,在4000Hz以上的子模态信号中选取振荡能量最大的子模态作为燃烧共振信号,振荡能量的评价标准为其中x为子模态信号,j为子模态信号的长度,k为子模态信号节点数,i为4000Hz以上子模态信号的序号。
3.步骤三具体包括:
将步骤二所得到共振模态信号进行希尔伯特变换,希尔伯特变换的解析信号表示为x(t)=m(t)+jl(t),其中m(t)为提取的共振信号,l(t)为共振信号的希尔伯特变换,l(t)表示为*表示卷积,则共振幅值和共振相位可以分别表示为共振频率通过对相位求导获得,即
本发明提出了一种精确计算柴油机燃烧共振频率的方法。对缸压信号进行傅里叶变换获得缸压频谱,利用经验小波变换的“尺度空间算法”对缸压频谱进行准确、自适应地分割并提取燃烧共振子模态,再利用希尔伯特变换处理该子模态得到燃烧共振频率及对应的共振振幅。经验小波变换的自适应性使本发明方法的计算结果不依赖于经验函数,并且对实际缸压信号的处理不存在模型收敛性的问题,拥有较高的精确度。
与现有技术相比,本发明优点与积极效果在于:经验小波变换的“尺度空间算法”能够准确、自适应地分割频谱从而使本方法不需要选择经验函数,并且结合希尔伯特变换计算共振频率的方法适用于任何条件下的缸压信号。相比于现有的共振频率提取技术具有收敛性好、自适应性强、精确度高的优点,获得的共振频率能够用于柴油机低噪声优化设计和缸内温度与气体质量估计。
附图说明
图1是本发明的流程图;
图2是子模态边界识别示意图;
图3是实际缸压信号频谱分割示意图;
图4是计算的燃烧共振频率与振荡幅值示意图。
具体实施方式
下面举例对本发明做更详细的描述。
结合图1,本发明的具体步骤如下:
步骤一:首先对缸压信号f(t)进行傅里叶变换
其中t为时间,f为频率。利用“尺度空间算法”确定缸压频谱中N个有意义的极大值点,并将其作为燃烧波动模态。具体为先定义核函数其中n为尺度参数,则离散尺度空间可以定义为
其中以便确保近似误差可以忽略。由于傅里叶谱经常混有噪声,因此尺度空间算法会移除特征长度小于的局部细节信息,从而获得频谱中粗略的模态信息。
根据“尺度空间算法”所定义的模态信息,设ω为每个模态的边界,并且ω={ωi}i=1,2,…,N,i为子模态编号,ω0=0,ωN=π。此外,以ωi为中心定义宽度为2τi的频率传递带,其中τi=γωi,并且0<γ<1。图2为频谱边界识别简图。
步骤二:根据Meyer小波的构造方法对经验小波滤波器进行构造,本质上依旧是带通滤波器。尺度函数和小波函数分别为
以及
其中β(x)是一任意函数:
因此对每个子模态进行经验小波变换得到的近似信号和细节信号分别为
将构造的经验小波滤波器与缸压信号做卷积,即自适应分割柴油机缸压频谱,分割结果如图3所示。在4000Hz以上的子模态中选取振荡能量最大的子模态作为燃烧共振信号,振荡能量的评价标准为:
其中x为子模态,j为子模态信号的长度,k为子模态信号节点数,i为4000Hz以上子模态信号的序号。
步骤三:由于选取的燃烧共振信号在每个时间点只有一个频率成分,故利用希尔伯特变换处理燃烧共振子信号。希尔伯特变换的解析信号可以表示为:
x(t)=m(t)+jl(t) (9)
其中m(t)为提取的共振信号,l(t)为共振信号的希尔伯特变换,l(t)可以表示为:
*表示卷积。则共振幅值和共振相位可以分别表示为:
然而瞬时频率(共振频率)可以通过对相位求导获得,即
因此,绘制希尔伯特谱如图4,可以直观地展示出计算的共振频率及其对应的振幅。
根据实际应用结果可知本发明方法中,经验小波变换的“尺度空间算法”能够准确且自适应地对缸压频谱进行分割,从而使本方法不需要选择可能会影响计算结果的经验函数,并且结合希尔伯特变换计算共振频率的方法适用于任何条件下的缸压信号,不存在算法收敛性的问题。与现有共振频率计算方法相比,用经验小波变换结合希尔伯特变换的共振频率计算方法,能够有效、精确地计算柴油机燃烧共振频率,从而能够更好地实现柴油机低噪声优化设计和缸内热力过程分析。

Claims (4)

1.一种基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法,其特征是:
步骤一:对缸压信号进行傅里叶变换获得缸压频谱,利用经验小波算法中的尺度空间算法识别缸压频谱中的特征子模态,并确定子模态的位置,计算子模态的边界频率,将整个频谱中对应每个燃烧波动的子模态边界识别出来;
步骤二:根据Meyer小波构造方法来构造经验小波滤波器,利用构造的小波滤波器对缸压信号进行带通滤波,提取缸压频谱中对应每个燃烧波动的子模态信号,在大于4000Hz的子模态信号中选取振荡能量最大的子模态,该信号即为燃烧共振信号;
步骤三:将步骤二所得到的燃烧共振信号进行希尔伯特变换,并绘制希尔伯特谱,获取共振频率及其对应的共振振幅。
2.根据权利要求1所述的基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法,其特征是步骤一具体包括:
首先对缸压信号f(t)进行傅里叶变换其中t为时间,f为频率,利用尺度空间算法确定缸压频谱中N个有意义的极大值点,并将其作为燃烧波动模态,具体为先定义核函数其中n为尺度参数,则离散尺度空间定义为其中由于傅里叶谱混有噪声,因此尺度空间算法移除特征长度小于的局部细节信息,从而获得频谱中粗略的模态信息,根据尺度空间算法所定义的模态信息,设ω为每个模态的边界,并且ω={ωi}i=1,2,…,N,i为子模态编号,ω0=0,ωN=π,此外,以ωi为中心定义宽度为2τi的频率传递带,其中τi=γωi,并且0<γ<1。
3.根据权利要求1所述的基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法,其特征是步骤二具体包括:
根据Meyer小波构造方法来构造经验小波滤波器,获得尺度函数为和小波函数将构造的经验小波滤波器与缸压信号做卷积,即自适应分割柴油机缸压频谱,在4000Hz以上的子模态信号中选取振荡能量最大的子模态作为燃烧共振信号,振荡能量的评价标准为其中x为子模态信号,j为子模态信号的长度,k为子模态信号节点数,i为4000Hz以上子模态信号的序号。
4.根据权利要求1所述的基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法,其特征是步骤三具体包括:
将步骤二所得到共振模态信号进行希尔伯特变换,希尔伯特变换的解析信号表示为x(t)=m(t)+jl(t),其中m(t)为提取的共振信号,l(t)为共振信号的希尔伯特变换,l(t)表示为*表示卷积,则共振幅值和共振相位可以分别表示为共振频率通过对相位求导获得,即
CN201810682880.3A 2018-06-28 2018-06-28 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法 Active CN108921082B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810682880.3A CN108921082B (zh) 2018-06-28 2018-06-28 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810682880.3A CN108921082B (zh) 2018-06-28 2018-06-28 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法

Publications (2)

Publication Number Publication Date
CN108921082A true CN108921082A (zh) 2018-11-30
CN108921082B CN108921082B (zh) 2021-11-12

Family

ID=64422961

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810682880.3A Active CN108921082B (zh) 2018-06-28 2018-06-28 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法

Country Status (1)

Country Link
CN (1) CN108921082B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110389037A (zh) * 2019-07-23 2019-10-29 哈尔滨工程大学 一种柴油机进气量测量系统
CN111666892A (zh) * 2020-06-08 2020-09-15 西南交通大学 一种基于经验小波希尔伯特变换的电力机车空转识别方法
CN113176460A (zh) * 2021-04-28 2021-07-27 河海大学 一种基于改进经验小波变换的电能质量扰动信号检测方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1235242A (zh) * 1998-05-11 1999-11-17 株式会社丰田自动织机制作所 抑制共振的方法和装置
CN101055935A (zh) * 2006-04-13 2007-10-17 株式会社摩拓尼斯 用于运载工具的多波段天线
CN102506995A (zh) * 2011-11-22 2012-06-20 中国建筑材料科学研究总院 一种基于hht变换和相关分析的振动信号处理方法
CN103137225A (zh) * 2013-01-25 2013-06-05 杭州电子科技大学 基于小波变换和希尔伯特变换的核电站松动部件定位方法
CN103767707A (zh) * 2014-02-17 2014-05-07 南通南大四维影通科技有限公司 血样水平依赖性功能磁共振信号震荡频率聚类分析方法
CN106499535A (zh) * 2015-09-08 2017-03-15 现代自动车株式会社 用于控制发动机的装置和方法
CN106568607A (zh) * 2016-11-04 2017-04-19 东南大学 一种基于经验小波变换碰摩声发射故障诊断方法
CN106662740A (zh) * 2014-09-19 2017-05-10 英特尔公司 紧凑型照明系统
CN107316653A (zh) * 2016-04-27 2017-11-03 南京理工大学 一种基于改进的经验小波变换的基频检测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1235242A (zh) * 1998-05-11 1999-11-17 株式会社丰田自动织机制作所 抑制共振的方法和装置
CN101055935A (zh) * 2006-04-13 2007-10-17 株式会社摩拓尼斯 用于运载工具的多波段天线
CN102506995A (zh) * 2011-11-22 2012-06-20 中国建筑材料科学研究总院 一种基于hht变换和相关分析的振动信号处理方法
CN103137225A (zh) * 2013-01-25 2013-06-05 杭州电子科技大学 基于小波变换和希尔伯特变换的核电站松动部件定位方法
CN103767707A (zh) * 2014-02-17 2014-05-07 南通南大四维影通科技有限公司 血样水平依赖性功能磁共振信号震荡频率聚类分析方法
CN106662740A (zh) * 2014-09-19 2017-05-10 英特尔公司 紧凑型照明系统
CN106499535A (zh) * 2015-09-08 2017-03-15 现代自动车株式会社 用于控制发动机的装置和方法
CN107316653A (zh) * 2016-04-27 2017-11-03 南京理工大学 一种基于改进的经验小波变换的基频检测方法
CN106568607A (zh) * 2016-11-04 2017-04-19 东南大学 一种基于经验小波变换碰摩声发射故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郑近德等: "基于改进经验小波变换的时频分析方法及其在滚动轴承故障诊断中的应用", 《电子学报》 *
马腾: "基于CEEMD-HT的汽油机爆震特征识别与强度评价研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110389037A (zh) * 2019-07-23 2019-10-29 哈尔滨工程大学 一种柴油机进气量测量系统
CN111666892A (zh) * 2020-06-08 2020-09-15 西南交通大学 一种基于经验小波希尔伯特变换的电力机车空转识别方法
CN111666892B (zh) * 2020-06-08 2023-04-25 西南交通大学 一种基于经验小波希尔伯特变换的电力机车空转识别方法
CN113176460A (zh) * 2021-04-28 2021-07-27 河海大学 一种基于改进经验小波变换的电能质量扰动信号检测方法

Also Published As

Publication number Publication date
CN108921082B (zh) 2021-11-12

Similar Documents

Publication Publication Date Title
CN108921082A (zh) 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法
CN106443178B (zh) 一种基于IQuinn-Rife综合的正弦信号频率估计方法
CN102519582B (zh) 航空发动机振动信号的盲源分离方法
EP2416293A1 (en) Stereo image processor and stereo image processing method
CN109285561B (zh) 一种基于自适应窗长的船舶螺旋桨空化噪声调制谱特征保真增强方法
CN106096313B (zh) 一种基于奇异谱分解和谱峭度的包络分析方法
CN107607065A (zh) 一种基于变分模态分解的冲击回波信号分析方法
CN107666322A (zh) 一种基于字典学习的自适应微震数据压缩感知方法
CN104297791A (zh) 一种基于地震优势频率的反演方法及系统
CN103956756A (zh) 一种电力系统低频振荡模态辨识方法
Tong et al. Parameter estimation of FH signals based on STFT and music algorithm
CN103743470A (zh) 一种汽车噪声频谱分析方法
CN104515905B (zh) 基于cqt多分辨率的被试的脑电信号自适应频谱分析方法
CN112098102A (zh) 一种基于ewt-scwt的内燃机异响识别与诊断方法
CN111580654A (zh) 一种基于emd的脑电信号的短时特征提取方法
CN101655834A (zh) 一种基于分数阶小波变换的信号分离方法
JP6090000B2 (ja) 周波数解析装置
Verma et al. Signal transforms for feature extraction from vibration signal for air compressor monitoring
CN113445992B (zh) 抽油机运动位移的处理方法和装置
CN105467270B (zh) 基于频谱相似性评价的单端行波测距反射波辨识算法
CN108241171B (zh) 一种采用复值Gauss积分滤波器对地震数据进行滤波和提取三瞬属性的方法
CN112434634A (zh) 一种快速消除土木工程结构健康监测信号尖峰的方法及系统
Li et al. Adaptive feature extraction algorithms and SVM with optimal parameters on fault diagnosis of bearing
CN110389037B (zh) 一种柴油机进气量测量系统
JP4843288B2 (ja) Am/fm混合信号処理装置

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