CN111863028B - 一种发动机声音合成方法及系统 - Google Patents

一种发动机声音合成方法及系统 Download PDF

Info

Publication number
CN111863028B
CN111863028B CN202010696931.5A CN202010696931A CN111863028B CN 111863028 B CN111863028 B CN 111863028B CN 202010696931 A CN202010696931 A CN 202010696931A CN 111863028 B CN111863028 B CN 111863028B
Authority
CN
China
Prior art keywords
engine sound
frequency
signals
signal
splicing
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.)
Active
Application number
CN202010696931.5A
Other languages
English (en)
Other versions
CN111863028A (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.)
Dragon Totem Technology Hefei Co ltd
Original Assignee
Jiangmen Polytechnic
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 Jiangmen Polytechnic filed Critical Jiangmen Polytechnic
Priority to CN202010696931.5A priority Critical patent/CN111863028B/zh
Publication of CN111863028A publication Critical patent/CN111863028A/zh
Application granted granted Critical
Publication of CN111863028B publication Critical patent/CN111863028B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60Q—ARRANGEMENT OF SIGNALLING OR LIGHTING DEVICES, THE MOUNTING OR SUPPORTING THEREOF OR CIRCUITS THEREFOR, FOR VEHICLES IN GENERAL
    • B60Q5/00—Arrangement or adaptation of acoustic signal devices
    • B60Q5/005—Arrangement or adaptation of acoustic signal devices automatically actuated
    • B60Q5/006—Arrangement or adaptation of acoustic signal devices automatically actuated indicating risk of collision between vehicles or with pedestrians
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60Q—ARRANGEMENT OF SIGNALLING OR LIGHTING DEVICES, THE MOUNTING OR SUPPORTING THEREOF OR CIRCUITS THEREFOR, FOR VEHICLES IN GENERAL
    • B60Q5/00—Arrangement or adaptation of acoustic signal devices
    • B60Q5/005—Arrangement or adaptation of acoustic signal devices automatically actuated
    • B60Q5/008—Arrangement or adaptation of acoustic signal devices automatically actuated for signaling silent vehicles, e.g. for warning that a hybrid or electric vehicle is approaching
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Mechanical Engineering (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Multimedia (AREA)
  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)

Abstract

本发明涉及一种发动机声音合成方法及系统。该方法包括：对每个发动机声音段进行低通滤波，得到低频谐波信号；对每个发动机声音段进行高通窄带滤波，得到高频窄带合成信号；采用希尔伯特变换方法，对每个发动机声音段的高频窄带信号进行分解，得到对应的包络信号和载波信号；将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号；将所有发动机声音段的载波信号进行拼接，得到载波拼接信号；将所有发动机声音段的包络信号进行拼接，得到包络拼接信号；将包络拼接信号和载波拼接信号相乘，得到高频窄带合成信号；将低频谐波拼接信号和高频窄带合成信号相加，合成发动机声音。本发明可以提高发动机声音合成的真实度。

Description

一种发动机声音合成方法及系统
技术领域
本发明涉及声音合成领域，特别是涉及一种发动机声音合成方法及系统。
背景技术
人们已经习惯于各种汽车声音的存在，而且能够分辨出各种发动机声音的品质。电动汽车行驶在道路上，安静的汽车在降低了道路噪声的同时，也为行人带来了潜在危险。因此，发动机声音的合成在多个领域有应用价值，例如，在电动汽车领域，合成的发动机声音可以作为警示音提醒行人；在虚拟驾驶领域，播放合成的发动机声音可以增加驾驶的真实感。
现有技术中，有的文献中公开发动机声音合成技术是将声音视为谐波与随机噪声组成。在频谱处理和恢复过程中，一些细节必然会丢失，而这些细节对于发动机声音的自然度往往非常重要。在频谱建模中，为降低计算复杂度，人们不得不使用简化模型。一般方法将发动机声音视为基于发动机转速的一系列低频谐波和高频随机噪声组成。这类方法能够有效的模拟发动机的低频声音，但对于高频声音来说，采用随机噪声处理，与真实的发动机声差异明显，合成的效果较差。
还有相关文献中公开发动机声音合成技术是基于样本的合成。该方法能够最大程度的保留原始噪声特征，根据需要选择合适的声音片段进行播放，不同声音片段之间的衔接是关键。为此，有人提出了基于PSOLA(基音同步叠加)的发动机声音合成方法。该方法假设各噪声由最小片段组成，片段宽度由基频决定。但由于PSOLA是在时域范围内根据基频确定的长度进行拼接，无法确保使每个频率的波形都同时对齐，在高频段相位不匹配，导致谱线衔接处出现跳变或产生虚假频率，此外由于发动机声音片段很短，当发动机转速缓慢变化或者稳定时，循环播放相同的声音会不自然，拼接段的相位不匹配的影响会更加明显。
发明内容
本发明的目的是提供一种发动机声音合成方法及系统，以提高发动机声音合成的真实度。
为实现上述目的，本发明提供了如下方案：
一种发动机声音合成方法，包括：
获取发动机声音样本；所述发动机声音样本包括多个发动机声音段；
对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号；
对每个所述发动机声音段进行高通窄带滤波，得到所述发动机声音段的多组高频窄带信号；
采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号；
将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号；
将所有发动机声音段的载波信号进行拼接，得到载波拼接信号；
将所有发动机声音段的包络信号进行拼接，得到包络拼接信号；
将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号；
将所述低频谐波拼接信号和所述高频窄带合成信号相加，合成发动机声音。
可选的，所述对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号，具体包括：
根据每个所述发动机声音段的基频，确定每个所述发动机声音段的三次谐波频率；
基于所述三次谐波频率，根据1/3倍频程的频带划分，确定每个所述发动机声音段对应的低通滤波器的频率上限；
采用所述低通滤波器对对应的所述发动机声音段进行低通滤波，得到每个所述发动机声音段对应的低频谐波信号；每个所述发动机声音段对应的低频谐波信号包括发动机基频的
Figure BDA0002591579780000021
阶次谐波信号、1阶次谐波信号、2阶次谐波信号和3阶次谐波信号。
可选的，所述对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号，具体包括：
对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，得到每个所述发动机声音段对应的多组高频窄带信号；
可选的，所述采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号，具体包括：
采用希尔伯特变换方法，对每个所述发动机声音段的每组高频窄带信号进行分解，得到每个所述发动机声音段中每组高频窄带信号对应的包络信号和载波信号。每组高频窄带信号对应的的包络信号均包括发动机基频的1/4阶次信号、1/2阶次信号、1阶次信号、2阶次信号和3阶次信号。
可选的，所述将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号，具体包括：
获取当前发动机声音段的低频谐波信号相位；
通过基音同步技术确定下一发动机声音段的低频谐波信号初始相位；
采用过零点对齐的方法使所述下一发动机声音段的低频谐波信号初始相位与所述当前发动机声音段的低频谐波信号相位对齐，在低频谐波信号拼接段对当前发动机声音段的低频谐波信号和下一发动机声音段的低频谐波信号进行拼接；低频谐波信号段拼接段的长度为基频波长的1/2；
依次对所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号。
可选的，所述将所有发动机声音段的载波信号进行拼接，得到载波拼接信号，具体包括：
对于第i组高频窄带信号，获取当前发动机声音段中第i组高频窄带信号的末端相位；i＝1、2、3、4、5…；
根据所述当前发动机声音段中第i组高频载波信号的末端相位，确定下一发动机声音段中第i组高频载波信号的起始点；所述下一发动机声音段中第i组高频载波信号的起始点的相位与所述当前发动机声音段中第i组高频载波信号的末端相位相同；
将所述当前发动机声音段中第i组高频载波信号的末端和所述下一发动机声音段中第i组高频载波信号的起始点对齐，在第i组高频载波信号载波拼接段，对当前发动机声音段中第i组载波信号和下一发动机声音段中第i组载波信号进行拼接；当第i组高频载波信号的频率小于2kHz时，第i组载波信号拼接段的长度为1/3倍频程中心频率波长的1/2；当第i组高频窄带信号的频率大于或等于2kHz时，第i组载波信号拼接段的长度为1/3倍频程中心频率波长；
依次对所有发动机声音段中第i组高频窄带信号的载波信号进行拼接，得到第i各组载波拼接信号；
依次得到每一组载波拼接信号。
可选的，所述将所有发动机声音段的包络信号进行拼接，得到包络拼接信号，具体包括：
对于第i组高频窄带信号，获取当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值；
根据当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值，确定下一发动机声音段中第i组高频窄带信号的峰值点；所述下一发动机声音段中第i组高频窄带信号的峰值点的峰值与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值相等；
将下一发动机声音段中第i组高频窄带信号的峰值点与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值位置对齐，在第i组包络信号拼接段，对当前发动机声音段中第i组高频窄带信号的包络信号和下一发动机声音段中第i组高频窄带信号的包络信号进行拼接；所述第i组包络信号拼接段为基频波长的1.5倍，且所述下一发动机声音段中第i组高频窄带信号的峰值点位于所述第i组包络信号拼接段的1/2位置；
依次对所有发动机声音段中第i组高频窄带信号的包络信号进行拼接，得到第i组包络拼接信号；
依次得到每一组包络拼接信号。
可选的，所述将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号，具体包括：
对于第i组高频窄带信号，将第i组包络拼接信号和第i组载波拼接信号相乘，得到第i组窄带信号；
将所有组窄带信号相加，得到所述高频窄带合成信号。
本发明还提供一种发动机声音合成系统，包括：
发动机声音样本获取模块，用于获取发动机声音样本；所述发动机声音样本包括多个发动机声音段；
低通滤波模块，用于对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号；
高通窄带滤波模块，用于对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号；
希尔伯特变换模块，用于采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号；
低频谐波信号拼接模块，用于将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号；
载波信号拼接模块，用于将所有发动机声音段的载波信号进行拼接，得到载波拼接信号；
包络信号拼接模块，用于将所有发动机声音段的包络信号进行拼接，得到包络拼接信号；
高频窄带信号计算模块，用于将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号；
发动机声音合成模块，用于将所述低频谐波拼接信号和所述高频窄带合成信号相加，合成发动机声音。
可选的，所述高通窄带滤波模块具体包括：
高通窄带滤波单元，用于对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，得到每个所述发动机声音段对应的多组高频窄带信号；
所述希尔伯特变换模块，具体包括：
分解单元，用于采用希尔伯特变换方法，对每个所述发动机声音段的每组高频窄带信号进行分解，得到每个所述发动机声音段中每组高频窄带信号对应的包络信号和载波信号。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明考虑了声音信号的低频谐波信号和高频窄带信号，将声音样本分解为低频谐波信号和高频窄带信号，然后对分解后的低频谐波信号和高频窄带信号分别进行拼接，最后进行发动机声音合成，信号失真度低，从而使合成的发动机信号更接近与真实声音，提高合成发动机声音的真实度。
此外，在对低频谐波信号和高频窄带信号拼接时，考虑了相位匹配问题，进而在相同相位处对低频谐波信号和高频窄带信号分别进行拼接，避免谱线衔接处出现跳变或产生虚假频率，提高发动机声音的真实度和自然度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明发动机声音合成方法的流程示意图；
图2为本发明低频谐波信号拼接的示意图；
图3为本发明载波信号拼接的示意图；
图4为本发明包络信号拼接的示意图；
图5为本发明发动机声音合成系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明发动机声音合成方法的流程示意图。如图1所示，本发明发动机声音合成方法包括以下步骤：
步骤100：获取发动机声音样本。所述发动机声音样本包括多个发动机声音段。例如，可以采集汽车在真实道路上行驶时，从怠速到2400转/分，以50转为间隔的15秒长发动机声音片段作为发动机声音段，进而得到多个15秒的发动机声音段，即发动机声音样本。具体的，根据实际情况和实际车辆进行相应的发动机声音样本采集。
步骤200：对每个发动机声音段进行低通滤波，得到每个发动机声音段的低频谐波信号。具体过程如下：
根据每个所述发动机声音段的基频，确定每个所述发动机声音段的三次谐波频率。例如，转速为2400转/分的发动机声音段对应的基频为80Hz，则该发动机声音段的三次谐波为240Hz。
基于所述三次谐波频率，根据1/3倍频程的频带划分，确定每个所述发动机声音段对应的低通滤波器的频率上限。例如，转速为2400转/分的发动机声音段对应低频滤波器的频率上限为282Hz。
采用所述低通滤波器对对应的所述发动机声音段进行低通滤波，得到每个所述发动机声音段对应的低频谐波信号。低于该低通滤波器频率上限的信号为低频谐波信号，每个所述发动机声音段对应的低频谐波信号包括发动机基频的
Figure BDA0002591579780000071
阶次谐波信号、1阶次谐波信号、2阶次谐波信号和3阶次谐波信号。
步骤300：对每个发动机声音段进行高通窄带滤波，得到每个发动机声音段的高频窄带信号。具体的，对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，可以得到每个所述发动机声音段对应的多组高频窄带信号。
步骤400：采用希尔伯特变换方法，对每个发动机声音段的高频窄带信号进行分解，得到每个发动机声音段高频窄带信号对应的包络信号和载波信号。本发明中对每一组高频窄带信号均进行分解得到包络信号和载波信号，即对每一组高频窄带信号均采用希尔伯特变换方法进行分解，进而可以得到每个发动机声音段中每组高频窄带信号对应的包络信号和载波信号。
步骤500：将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号。经过低通滤波后，低频谐波信号以基频为主。如图2所示，当与下一个发动机声音段的低频谐波信号进行拼接时，首先获取当前发动机声音段的低频谐波信号的相位，进而通过基音同步技术，确定下一发动机声音段的低频谐波信号的初始相位。然后，在低频谐波信号拼接段对相邻两个发动机声音段的低频谐波信号进行拼接，在低频谐波信号拼接段的长度为基频波长的1/2。考虑到两段低频谐波信号虽以基频为主，但也包含其它频率成分，两段新会的幅值相位会有一些差异，因此采用过零点对齐使相邻两个发动机声音段的低频谐波信号的相位对齐，即使得下一发动机声音段的低频谐波信号初始相位与当前发动机声音段的低频谐波信号相位对齐，然后再进行拼接。依次对所有相邻两个发动机声音段的低频谐波信号进行拼接后，得到低频谐波拼接信号。
步骤600：将所有发动机声音段的各组载波信号分别进行拼接，得到载波拼接信号。如图3所示，本步骤对载波信号进行拼接时，是将相邻发动机声音段的同一组高频窄带信号的载波信号进行拼接。由于所有组高频窄带信号对应的载波信号
Figure BDA0002591579780000081
的幅值均为1，因此，对载波信号进行拼接时，只需要对齐前后载波信号的相位即可。具体过程如下：
对于第i组高频窄带信号，获取当前发动机声音段中第i组载波信号的末端相位；i＝1、2、3、4、5…。第i组载波信号为第i组高频窄带信号对应的载波信号。
根据所述当前发动机声音段中第i组载波信号的末端相位，确定下一发动机声音段中第i组载波信号的起始点。所述下一发动机声音段中第i组载波信号的起始点的相位与所述当前发动机声音段中第i组载波信号的末端相位相同。
将所述当前发动机声音段中第i组载波信号的末端和所述下一发动机声音段中第i组载波信号的起始点对齐，在第i组载波信号拼接段，对当前发动机声音段中第i组载波信号和下一发动机声音段中第i组载波信号进行拼接。当第i组载波信号的频率小于2kHz时，第i组载波信号拼接段的长度为1/3倍频程中心频率波长的1/2；当第i组载波信号的频率大于或等于2kHz时，由于采样点较少，第i组载波信号拼接段的长度为1/3倍频程中心频率波长。
依次对所有相邻的发动机声音段中第i组载波信号进行拼接，得到第i组载波拼接信号，第i组载波拼接信号是将所有发动机声音段中第i组载波信号全部拼接的结果。
按照上述方式，依次得到每一组载波拼接信号，即依次得到第1组载波拼接信号、第2组载波拼接信号、第2组载波拼接信号、第3组载波拼接信号和第4组载波拼接信号，以此类推。
步骤700：将所有发动机声音段的包络信号进行拼接，得到包络拼接信号。如图4所示，与载波信号拼接过程类似，包络信号拼接时，也是将相邻发动机声音段的同一组高频窄带信号的包络信号进行拼接。高频窄带信号中的包络信号值均为正值，包络信号频率以发动机基频为主，包含1/4、1/2阶次频率及谐波。考虑到包络信号中有基频的1/4和1/2阶次频率存在，且均为正值，幅值变化较大，因此增加包络拼接段的长度，将包络拼接段设定为基频波长的1.5倍。具体的，包络信号拼接过程如下：
对于第i组高频窄带信号，获取当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值，通常一个基频周期内约包括2个峰值点。
根据当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值，确定下一发动机声音段中第i组高频窄带信号的峰值点；所述下一发动机声音段中第i组高频窄带信号的峰值点的峰值与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值相等。
将下一发动机声音段中第i组高频窄带信号的峰值点与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值位置对齐，在第i组包络信号拼接段，对当前发动机声音段中第i组高频窄带信号的包络信号和下一发动机声音段中第i组高频窄带信号的包络信号进行拼接；所述第i组包络信号拼接段为基频波长的1.5倍，且所述下一发动机声音段中第i组高频窄带信号的峰值点位于所述第i组包络信号拼接段的1/2位置，即在峰值点左右分别延展1/4基频波长为包络拼接段。
依次对所有发动机声音段中第i组高频窄带信号的包络信号进行拼接，得到第i组包络拼接信号。
按照上述步骤依次得到每一组包络拼接信号，即依次得到第1组包络拼接信号、第2组包络拼接信号、第3组包络拼接信号、第4组包络拼接信号和第5组包络拼接信号，以此类推。
步骤800：将各组包络拼接信号和载波拼接信号分别相乘，然后相加，得到高频窄带合成信号。具体的，将相同组的包络拼接信号和载波拼接信号相乘，可以得到对应组的高频窄带信号。例如，对于第i组高频窄带信号，将第i组包络拼接信号和第i组载波拼接信号相乘，得到第i组窄带信号。然后将所有组窄带信号相加，可以得到高频窄带合成信号。
步骤900：将低频谐波拼接信号和高频窄带合成信号相加，合成发动机声音。合成的发动机声音模型如下：
Figure BDA0002591579780000101
式中，
Figure BDA0002591579780000102
为低频谐波信号，n为谐波阶次，An为第n阶次谐波幅值，ω0为基频，
Figure BDA0002591579780000103
为第n阶次谐波初始相位；
Figure BDA0002591579780000104
为高频窄带信号，
Figure BDA0002591579780000105
为第i组载波信号
Figure BDA0002591579780000106
为第i组包络信号，Bi为第i组高频窄带信号的载波信号幅值，ωi为第i组高频窄带信号的载波中心频率，
Figure BDA0002591579780000107
为第i组高频窄带信号的载波初始相位，Mi为第i组高频窄带信号的调幅因子，i为高频窄带信号的组号，kij为发动机阶次，
Figure BDA0002591579780000108
为第i组高频窄带信号中kij阶次信号的相位，N为高频窄带信号分组的个数。
本发明还提供一种发动机声音合成系统，图5为本发明发动机声音合成系统的结构示意图。如图5所示，本发明发动机声音合成系统包括以下结构：
发动机声音样本获取模块501，用于获取发动机声音样本；所述发动机声音样本包括多个发动机声音段。
低通滤波模块502，用于对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号。
高通窄带滤波模块503，用于对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号。
希尔伯特变换模块504，用于采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号。
低频谐波信号拼接模块505，用于将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号。
载波信号拼接模块506，用于将所有发动机声音段的载波信号进行拼接，得到载波拼接信号。
包络信号拼接模块507，用于将所有发动机声音段的包络信号进行拼接，得到包络拼接信号。
高频窄带合成信号计算模块508，用于将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号。
发动机声音合成模块509，用于将所述低频谐波拼接信号和所述高频窄带合成信号相加，合成发动机声音。
作为另一实施例，本发明发动机声音合成系统中，所述高通窄带滤波模块503具体包括：高通窄带滤波单元，用于对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，得到每个所述发动机声音段对应的多组高频窄带信号。
所述希尔伯特变换模块504，具体包括：分解单元，用于采用希尔伯特变换方法，对每个所述发动机声音段的每组高频窄带信号进行分解，得到每个所述发动机声音段中每组高频窄带信号对应的包络信号和载波信号。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种发动机声音合成方法，其特征在于，包括：
获取发动机声音样本；所述发动机声音样本包括多个发动机声音段；
对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号；
对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号；
采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号；
将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号；
将所有发动机声音段的载波信号进行拼接，得到载波拼接信号；
将所有发动机声音段的包络信号进行拼接，得到包络拼接信号；
将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号；
将所述低频谐波拼接信号和所述高频窄带合成信号相加，合成发动机声音。
2.根据权利要求1所述的发动机声音合成方法，其特征在于，所述对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号，具体包括：
根据每个所述发动机声音段的基频，确定每个所述发动机声音段的三次谐波频率；
基于所述三次谐波频率，根据1/3倍频程的频带划分，确定每个所述发动机声音段对应的低通滤波器的频率上限；
采用所述低通滤波器对对应的所述发动机声音段进行低通滤波，得到每个所述发动机声音段对应的低频谐波信号；每个所述发动机声音段对应的低频谐波信号包括发动机基频的12阶次谐波信号、1阶次谐波信号、2阶次谐波信号和3阶次谐波信号。
3.根据权利要求1所述的发动机声音合成方法，其特征在于，所述对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号，具体包括：
对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，得到每个所述发动机声音段对应的多组高频窄带信号。
4.根据权利要求3所述的发动机声音合成方法，其特征在于，所述采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号，具体包括：
采用希尔伯特变换方法，对每个所述发动机声音段的每组高频窄带信号进行分解，得到每个所述发动机声音段中每组高频窄带信号对应的包络信号和载波信号；每组高频窄带信号对应的包络信号均包括发动机基频的1/4阶次信号、1/2阶次信号、1阶次信号、2阶次信号和3阶次信号。
5.根据权利要求1所述的发动机声音合成方法，其特征在于，所述将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号，具体包括：
获取当前发动机声音段的低频谐波信号相位；
通过基音同步技术确定下一发动机声音段的低频谐波信号初始相位；
采用过零点对齐的方法使所述下一发动机声音段的低频谐波信号初始相位与所述当前发动机声音段的低频谐波信号相位对齐，在低频谐波信号拼接段对当前发动机声音段的低频谐波信号和下一发动机声音段的低频谐波信号进行拼接；低频谐波信号段拼接段的长度为基频波长的1/2；
依次对所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号。
6.根据权利要求4所述的发动机声音合成方法，其特征在于，所述将所有发动机声音段的载波信号进行拼接，得到载波拼接信号，具体包括：
对于第i组高频窄带信号，获取当前发动机声音段中第i组载波信号的末端相位；第i组载波信号为第i组高频窄带信号对应的载波信号；
根据所述当前发动机声音段中第i组载波信号的末端相位，确定下一发动机声音段中第i组载波信号的起始点；所述下一发动机声音段中第i组载波信号的起始点的相位与所述当前发动机声音段中第i组载波信号的末端相位相同；
将所述当前发动机声音段中第i组载波信号的末端和所述下一发动机声音段中第i组载波信号的起始点对齐，在第i组载波信号拼接段，对当前发动机声音段中第i组载波信号和下一发动机声音段中第i组载波信号进行拼接；当第i组载波信号的频率小于2kHz时，第i组载波信号拼接段的长度为1/3倍频程中心频率波长的1/2；当第i组载波信号的频率大于或等于2kHz时，第i组载波信号拼接段的长度为1/3倍频程中心频率波长；
依次对所有发动机声音段中第i组载波信号进行拼接，得到第i组载波拼接信号；
依次得到每一组载波拼接信号。
7.根据权利要求6所述的发动机声音合成方法，其特征在于，所述将所有发动机声音段的包络信号进行拼接，得到包络拼接信号，具体包括：
对于第i组高频窄带信号，获取当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值；
根据当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值，确定下一发动机声音段中第i组高频窄带信号的峰值点；所述下一发动机声音段中第i组高频窄带信号的峰值点的峰值与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值相等；
将下一发动机声音段中第i组高频窄带信号的峰值点与当前发动机声音段中第i组高频窄带信号中一个基频周期内的最大峰值位置对齐，在第i组包络信号拼接段，对当前发动机声音段中第i组高频窄带信号的包络信号和下一发动机声音段中第i组高频窄带信号的包络信号进行拼接；所述第i组包络信号拼接段为基频波长的1.5倍，且所述下一发动机声音段中第i组高频窄带信号的峰值点位于所述第i组包络信号拼接段的1/2位置；
依次对所有发动机声音段中第i组高频窄带信号的包络信号进行拼接，得到第i组包络拼接信号；
依次得到每一组包络拼接信号。
8.根据权利要求7所述的发动机声音合成方法，其特征在于，所述将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号，具体包括：
对于第i组高频窄带信号，将第i组包络拼接信号和第i组载波拼接信号相乘，得到第i组窄带信号；
将所有组窄带信号相加，得到所述高频窄带合成信号。
9.一种发动机声音合成系统，其特征在于，包括：
发动机声音样本获取模块，用于获取发动机声音样本；所述发动机声音样本包括多个发动机声音段；
低通滤波模块，用于对每个所述发动机声音段进行低通滤波，得到每个所述发动机声音段的低频谐波信号；
高通窄带滤波模块，用于对每个所述发动机声音段进行高通窄带滤波，得到每个所述发动机声音段的高频窄带信号；
希尔伯特变换模块，用于采用希尔伯特变换方法，对每个所述发动机声音段的高频窄带信号进行分解，得到每个所述发动机声音段的对应的包络信号和载波信号；
低频谐波信号拼接模块，用于将所有发动机声音段的低频谐波信号进行拼接，得到低频谐波拼接信号；
载波信号拼接模块，用于将所有发动机声音段的载波信号进行拼接，得到载波拼接信号；
包络信号拼接模块，用于将所有发动机声音段的包络信号进行拼接，得到包络拼接信号；
高频窄带信号计算模块，用于将所述包络拼接信号和所述载波拼接信号相乘，得到高频窄带合成信号；
发动机声音合成模块，用于将所述低频谐波拼接信号和所述高频窄带合成信号相加，合成发动机声音。
10.根据权利要求9所述的发动机声音合成系统，其特征在于，所述高通窄带滤波模块具体包括：
高通窄带滤波单元，用于对每个所述发动机声音段，按1/3倍频程进行高通窄带滤波，得到每个所述发动机声音段对应的多组高频窄带信号；
所述希尔伯特变换模块，具体包括：
分解单元，用于采用希尔伯特变换方法，对每个所述发动机声音段的每组高频窄带信号进行分解，得到每个所述发动机声音段中每组高频窄带信号对应的包络信号和载波信号。
CN202010696931.5A 2020-07-20 2020-07-20 一种发动机声音合成方法及系统 Active CN111863028B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010696931.5A CN111863028B (zh) 2020-07-20 2020-07-20 一种发动机声音合成方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010696931.5A CN111863028B (zh) 2020-07-20 2020-07-20 一种发动机声音合成方法及系统

Publications (2)

Publication Number Publication Date
CN111863028A CN111863028A (zh) 2020-10-30
CN111863028B true CN111863028B (zh) 2023-05-09

Family

ID=73001011

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010696931.5A Active CN111863028B (zh) 2020-07-20 2020-07-20 一种发动机声音合成方法及系统

Country Status (1)

Country Link
CN (1) CN111863028B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000149167A (ja) * 1998-11-10 2000-05-30 Nissan Motor Co Ltd 警報装置
JP2005010621A (ja) * 2003-06-20 2005-01-13 Matsushita Electric Ind Co Ltd 音声帯域拡張装置及び帯域拡張方法
JP2005275420A (ja) * 2005-04-28 2005-10-06 Yamaha Corp 音声分析及び合成装置、方法、プログラム
JP2009063968A (ja) * 2007-09-10 2009-03-26 Yamaha Corp エンジン音合成装置
CN201721384U (zh) * 2010-04-19 2011-01-26 小山工业株式会社 电动车辆用模拟发动机声音产生装置
CN102113048A (zh) * 2008-07-11 2011-06-29 株式会社音响设计 合成声音
CN103544949A (zh) * 2012-07-12 2014-01-29 哈曼贝克自动系统股份有限公司 发动机声音合成
CN205652027U (zh) * 2015-01-13 2016-10-19 Ask工业股份公司 车辆的发动机声音增强系统
CN110108467A (zh) * 2019-05-20 2019-08-09 长沙理工大学 基于便携式移动设备的主动发声测速方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014021318A1 (ja) * 2012-08-01 2014-02-06 独立行政法人産業技術総合研究所 音声分析合成のためのスペクトル包絡及び群遅延の推定システム及び音声信号の合成システム

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000149167A (ja) * 1998-11-10 2000-05-30 Nissan Motor Co Ltd 警報装置
JP2005010621A (ja) * 2003-06-20 2005-01-13 Matsushita Electric Ind Co Ltd 音声帯域拡張装置及び帯域拡張方法
JP2005275420A (ja) * 2005-04-28 2005-10-06 Yamaha Corp 音声分析及び合成装置、方法、プログラム
JP2009063968A (ja) * 2007-09-10 2009-03-26 Yamaha Corp エンジン音合成装置
CN102113048A (zh) * 2008-07-11 2011-06-29 株式会社音响设计 合成声音
CN201721384U (zh) * 2010-04-19 2011-01-26 小山工业株式会社 电动车辆用模拟发动机声音产生装置
CN103544949A (zh) * 2012-07-12 2014-01-29 哈曼贝克自动系统股份有限公司 发动机声音合成
CN205652027U (zh) * 2015-01-13 2016-10-19 Ask工业股份公司 车辆的发动机声音增强系统
CN110108467A (zh) * 2019-05-20 2019-08-09 长沙理工大学 基于便携式移动设备的主动发声测速方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Ming C.Lin.Physically Based Sound Synthesis for Large-Scale Virtual Environments.《IEEE computer Graphics and Applications》.2007,全文. *
贾炎秋.汽车引擎声音模拟系统设计.《中国优秀硕士学位论文全文数据库》.2020,(第第4期期),全文. *

Also Published As

Publication number Publication date
CN111863028A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
EP2685448B1 (en) Engine sound synthesis
CN103943113A (zh) 一种歌曲去伴奏的方法和装置
WO1990013887A1 (en) Musical signal analyzer and synthesizer
CN111798831B (zh) 一种声音粒子合成方法及装置
CN102750955B (zh) 基于残差信号频谱重构的声码器
CN104783780A (zh) 心电信号除噪方法及装置
CN103198836B (zh) 基于Hilbert变换的欠定声音信号分离方法及装置
CN111863028B (zh) 一种发动机声音合成方法及系统
CN111739544B (zh) 语音处理方法、装置、电子设备及存储介质
WO2019062197A1 (zh) 一种新能源车电机噪声信号提取方法及系统
CN116206624B (zh) 一种车辆声浪合成方法、装置、存储介质及设备
Disch et al. Modulation and delay line based digital audio effects
CN107017005A (zh) 一种基于dft的双通道语声分离方法
Parker Efficient dispersion generation structures for spring reverb emulation
US10825443B2 (en) Method and system for implementing a modal processor
Pendharkar Auralization of road vehicles using spectral modeling synthesis
JP2002175099A (ja) 雑音抑制方法および雑音抑制装置
Hanna et al. Time scale modification of noises using a spectral and statistical model
US11488574B2 (en) Method and system for implementing a modal processor
CN113593592B (zh) 基于多域扩张的音频增强方法及装置
CN117351927A (zh) 一种汽车加速声的合成方法、系统及设备
CN115395981B (zh) 一种高速跳频信号的数学建模与生成方法
CN110287948B (zh) 一种基于能量分离的魏格纳-维利时频分解方法
AU2020200340B2 (en) Cross Product Enhanced Subband Block Based Harmonic Transposition
CN117524183A (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
TR01 Transfer of patent right

Effective date of registration: 20240117

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Address before: 529000 No. 6, ChaoLian Avenue, ChaoLian Town, Pengjiang district, Jiangmen City, Guangdong Province

Patentee before: JIANGMEN POLYTECHNIC

TR01 Transfer of patent right