CN112304632A - Transient modulation evaluation method for describing human ear perception - Google Patents
Transient modulation evaluation method for describing human ear perception Download PDFInfo
- Publication number
- CN112304632A CN112304632A CN202010654276.7A CN202010654276A CN112304632A CN 112304632 A CN112304632 A CN 112304632A CN 202010654276 A CN202010654276 A CN 202010654276A CN 112304632 A CN112304632 A CN 112304632A
- Authority
- CN
- China
- Prior art keywords
- modulation
- matrix
- sound signal
- frequency
- human ear
- 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.)
- Pending
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 27
- 230000008447 perception Effects 0.000 title claims abstract description 15
- 238000011156 evaluation Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 230000005236 sound signal Effects 0.000 claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 15
- 238000000034 method Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 3
- 210000005069 ears Anatomy 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; 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
- G10L25/51—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention relates to a transient modulation evaluation method for describing human ear perception, which comprises the following steps: synchronously acquiring an engine rotating speed signal and an in-vehicle sound signal, and dividing an original sound signal matrix S according to a self-defined rotating speed increment; a weighting filtering is carried out on the sound signals to obtain a filtered sound signal matrix SA(ii) a For the filtered sound signal matrix SAPerforming Hilbert transform and calculating to obtain corresponding envelope matrix EA(ii) a For envelope matrix EAFourier transformation is carried out, and a modulation depth matrix D under different modulation frequencies is obtained through calculationA(ii) a And determining the modulation depth of each modulation order under different rotating speeds based on a peak value holding principle according to the corresponding relation between the modulation frequency and the modulation order. The invention can obtain the modulation depth of the modulation order under different transient working conditionsAnd (3) describing transient modulation components perceived by human ears, and providing a parameterized evaluation method for the whole vehicle dynamic sound quality benchmarks of different styles.
Description
Technical Field
The invention relates to the field of sound quality analysis and evaluation, in particular to a transient modulation evaluation method for describing human ear perception.
Background
With the development of automobile NVH control technology, the characteristic of dynamic sound quality is created, and the transmission of brand DNA becomes a development hotspot of each large automobile enterprise. Relevant studies have shown that the quality of sporty sound should be sound with the right half-order ratio and forming certain "sense of modulation" attributes.
Currently, methods for objectively evaluating a sound modulation phenomenon are less researched. Chinese patent publication No. CN101598596A proposes a method for analyzing noise sources, particularly vehicle noise, and realizes objective description of impulse type noise phenomenon by calculation of modulation characteristic parameters. The method has the disadvantages that the band-pass filtering needs to be adopted for multiple times to carry out frequency band decomposition on input signals respectively, and then modulation information is solved, when the transient working condition is analyzed by the processing mode, because a sound signal matrix needs to be established according to the rotating speed, the analysis efficiency is obviously reduced by the multiple times of band-pass filtering, and the current engineering application is only limited to vector type data analysis such as steady-state working conditions.
Disclosure of Invention
It is an object of the present invention to propose a method for evaluating transient modulations describing the perception of the human ear, in order to alleviate or eliminate the above-mentioned technical problems.
The invention relates to a transient modulation evaluation method for describing human ear perception, which comprises the following steps:
step 1, synchronously acquiring an engine rotating speed signal and an in-vehicle sound signal, and dividing an original sound signal matrix S according to a self-defined rotating speed increment;
Step 3, a filtered sound signal matrix SAPerforming Hilbert transform and calculating to obtain correspondingEnvelope matrix EA;
Step 4, aiming at envelope line matrix EAFourier transformation is carried out, and a modulation depth matrix D under different modulation frequencies is obtained through calculationA;
And 5, determining the modulation depth of each modulation order under different rotating speeds based on a peak value holding principle according to the corresponding relation between the modulation frequency and the modulation order.
Further, the step 1 specifically comprises:
collecting engine speed signal and sound signal with the same time sampling rate fs according to initial speed R0And the rotation speed increment delta R determines the central point N (R) of each data blockn)midThen, the start and stop points of each data block are determined according to the frequency resolution df:
N(Rn)1=N(Rn)mid-fs/2df;N(Rn)end=N(Rn)mid+fs/2df;
wherein n is 1,2, …, nR,nRThe total number of data blocks of the rotating speed sequence;
the finally divided sound signal matrix is S.
Further, step 2 specifically comprises: first, define the coefficient w of A weighting filterA(f):
Wherein f represents frequency;
then, the sound signal is FFT-transformed: f (f) is FFT [ S ], where f (f) is the corresponding spectrum matrix;
and then, weighting A to the frequency spectrum matrix, and obtaining a filtered sound signal matrix S through FFT inverse transformationA:
Further, step 3 specifically comprises:
to SAHilbert transformation is performed in sections, and an absolute value is taken to obtain an envelope matrix EA:
EA=|Hilbert[SA]|。
Further, step 4 specifically comprises:
first, for EAAnd carrying out FFT (fast Fourier transform) on the segments, and taking an absolute value to obtain an amplitude spectrum:
wherein, FAIs a corresponding spectrum matrix, A0Is an amplitude matrix corresponding to 0Hz, AiAn amplitude matrix corresponding to the non-zero frequency;
then, according to A0And AiCalculating a modulation depth matrix DA:
Further, step 5 specifically comprises:
firstly, according to the modulation frequency fmDetermining the modulation order O from the speed Rm:
Wherein f ismIs a frequency spectrum matrix FAA corresponding modulation frequency;
then, according to the customized order width OwDetermining the upper limit and the lower limit of the modulation frequency:
and then determining the modulation depth of each modulation order under different rotating speeds based on the peak value holding principle:
according to the invention, the original noise signal is processed in a segmented manner according to the rotating speed, A weighting filtering, Hilbert transform and fast Fourier transform are carried out successively to obtain a transient modulation cloud picture, and the modulation depth of the modulation order under different transient working conditions is extracted and obtained based on the peak value holding principle by further combining the corresponding relation of the rotating speed with the same order, so that the transient modulation component sensed by human ears is described, and a parameterized evaluation method is provided for the whole vehicle dynamic sound quality benchmarks of different styles.
Drawings
Fig. 1 is an algorithm flow chart of transient modulation analysis.
Fig. 2 is a transient modulation cloud picture of noise signals in a certain sports car.
Fig. 3 is a transient modulation cloud chart of noise signals in a comfortable car.
Fig. 4 is a comparison of the 1 st order modulation curves of the in-vehicle noise signal.
Detailed Description
The invention will be further explained with reference to the drawings.
A method for evaluating transient modulation describing human ear perception as shown in fig. 1, comprising the steps of:
step 1, synchronously acquiring an engine rotating speed signal and an in-vehicle sound signal, and dividing an original sound signal matrix S according to a self-defined rotating speed increment; the step 1 specifically comprises the following steps:
collecting engine speed signal and sound signal with the same time sampling rate fs according to initial speed R0And the rotation speed increment delta R determines the central point N (R) of each data blockn)midThen, the start and stop points of each data block are determined according to the frequency resolution df:
N(Rn)1=N(Rn)mid-fs/2df;N(Rn)end=N(Rn)mid+fs/2df;
wherein n is 1,2, …, nR,nRThe total number of data blocks of the rotating speed sequence;
the finally divided sound signal matrix is S.
first, define the coefficient w of A weighting filterA(f):
Wherein f represents frequency;
then, the sound signal is FFT-transformed: f (f) is FFT [ S ], where f (f) is the corresponding spectrum matrix;
and then, weighting A to the frequency spectrum matrix, and obtaining a filtered sound signal matrix S through FFT inverse transformationA:
Step 3, a filtered sound signal matrix SAPerforming Hilbert transform and calculating to obtain corresponding envelope matrix EA(ii) a The step 3 specifically comprises the following steps:
to SAHilbert transformation is performed in sections, and an absolute value is taken to obtain an envelope matrix EA:
EA=|Hilbert[SA]|。
Step 4, aiming at envelope line matrix EAFourier transformation is carried out, and a modulation depth matrix D under different modulation frequencies is obtained through calculationA(ii) a The step 4 specifically comprises the following steps:
first, for EAAnd carrying out FFT (fast Fourier transform) on the segments, and taking an absolute value to obtain an amplitude spectrum:
wherein, FAIs a corresponding spectrum matrix, A0Is an amplitude matrix (i.e. DC component matrix) corresponding to 0Hz, AiAn amplitude matrix (i.e., an AC component matrix) corresponding to a non-zero frequency;
then, according to A0And AiCalculating a modulation depth matrix DA:
Step 5, determining the modulation depth of each modulation order under different rotating speeds based on a peak value holding principle according to the corresponding relation between the modulation frequency and the modulation order, wherein the step 5 specifically comprises the following steps:
firstly, according to the modulation frequency fmAnd the rotational speed R (in rpm) determines the modulation order Om:
Wherein f ismIs a frequency spectrum matrix FAA corresponding modulation frequency;
then, according to the customized order width OwDetermining the upper limit and the lower limit of the modulation frequency:
and then determining the modulation depth of each modulation order under different rotating speeds based on the peak value holding principle:
the transient modulation evaluation method for describing human ear perception is applied to a certain sports car and a certain comfortable car for verification.
Fig. 2 is a transient modulation cloud chart (three-dimensional graph of rotation speed, modulation frequency and modulation depth, which shows that the order modulation characteristics are significant, and low-order modulations of 0.5 order, 1 order, 1.5 order, 2 order and the like are most significant) of noise signals in a certain sports car.
Fig. 3 is a transient modulation cloud chart (three-dimensional graph of rotation speed, modulation frequency and modulation depth) of a noise signal in a comfortable car, and it can be found that the order modulation characteristic is not obvious relative to a sports car.
Fig. 4 is a comparison of 1-order modulation curves of noise signals in a car, and it can be found that the modulation depth of a comfortable car in a full rotation speed range changes smoothly, while the modulation depth of a sports car in a rotation speed range above 2krpm is obviously increased, and the subjective feeling of the human ears on the particles can be effectively evaluated.
Therefore, the transient modulation component perceived by the human ear can be quantized according to the transient modulation evaluation method for describing human ear perception, and a parameterized evaluation method is provided for the whole vehicle dynamic sound quality benchmarks of different styles.
Claims (6)
1. A transient modulation evaluation method for describing human ear perception is characterized by comprising the following steps:
step 1, synchronously acquiring an engine rotating speed signal and an in-vehicle sound signal, and dividing an original sound signal matrix S according to a self-defined rotating speed increment;
step 2, carrying out A weighting filtering on the sound signals to obtain a filtered sound signal matrix SA;
Step 3, a filtered sound signal matrix SAPerforming Hilbert transform and calculating to obtain corresponding envelope matrix EA;
Step 4, aiming at envelope line matrix EAFourier transformation is carried out, and a modulation depth matrix D under different modulation frequencies is obtained through calculationA;
And 5, determining the modulation depth of each modulation order under different rotating speeds based on a peak value holding principle according to the corresponding relation between the modulation frequency and the modulation order.
2. The method for evaluating transient modulation describing human ear perception according to claim 1, wherein the step 1 specifically comprises:
collecting engine speed signal and sound signal with the same time sampling rate fs according to initial speed R0And the rotation speed increment delta R determines the central point N (R) of each data blockn)midThen, the start and stop points of each data block are determined according to the frequency resolution df:
N(Rn)1=N(Rn)mid-fs/2df;N(Rn)end=N(Rn)mid+fs/2df;
wherein n is 1,2, …, nR,nRThe total number of data blocks of the rotating speed sequence;
the finally divided sound signal matrix is S.
3. The method for evaluating transient modulation describing human ear perception according to claim 1, wherein the step 2 specifically comprises: first, define the coefficient w of A weighting filterA(f):
Wherein f represents frequency;
then, the sound signal is FFT-transformed: f (f) is FFT [ S ], where f (f) is the corresponding spectrum matrix;
and then, weighting A to the frequency spectrum matrix, and obtaining a filtered sound signal matrix S through FFT inverse transformationA:
4. The method for evaluating transient modulation describing human ear perception according to claim 1, wherein step 3 specifically comprises:
to SAHilbert transformation is performed in sections, and an absolute value is taken to obtain an envelope matrix EA:
EA=|Hilbert[SA]|。
5. The method for evaluating transient modulation describing human ear perception according to claim 1, wherein step 4 specifically comprises:
first, for EAAnd carrying out FFT (fast Fourier transform) on the segments, and taking an absolute value to obtain an amplitude spectrum:
wherein, FAIs a corresponding spectrum matrix, A0Is an amplitude matrix corresponding to 0Hz, AiAn amplitude matrix corresponding to the non-zero frequency;
then, according to A0And AiCalculating a modulation depth matrix DA:
6. The method for evaluating transient modulation describing human ear perception according to claim 1, wherein step 5 specifically comprises:
firstly, according to the modulation frequency fmDetermining the modulation order O from the speed Rm:
Wherein f ismIs a frequency spectrum matrix FAA corresponding modulation frequency;
then, according to the customized order width OwDetermining the upper limit and the lower limit of the modulation frequency:
and then determining the modulation depth of each modulation order under different rotating speeds based on the peak value holding principle:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010654276.7A CN112304632A (en) | 2020-07-08 | 2020-07-08 | Transient modulation evaluation method for describing human ear perception |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010654276.7A CN112304632A (en) | 2020-07-08 | 2020-07-08 | Transient modulation evaluation method for describing human ear perception |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112304632A true CN112304632A (en) | 2021-02-02 |
Family
ID=74483495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010654276.7A Pending CN112304632A (en) | 2020-07-08 | 2020-07-08 | Transient modulation evaluation method for describing human ear perception |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112304632A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113758713A (en) * | 2021-08-09 | 2021-12-07 | 重庆长安汽车股份有限公司 | Adaptive rough acoustic frequency band identification method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598596A (en) * | 2008-06-06 | 2009-12-09 | 福特环球技术公司 | Analyze the particularly method and apparatus of the noise of vehicle of noise source |
CN107642426A (en) * | 2017-08-31 | 2018-01-30 | 清华大学苏州汽车研究院(相城) | A kind of automobile engine noise initiative control method and system |
CN109194306A (en) * | 2018-08-28 | 2019-01-11 | 重庆长安汽车股份有限公司 | A kind of method and device quantifying automobile noise modulation problems |
CN109211570A (en) * | 2018-11-12 | 2019-01-15 | 吉林大学 | A kind of electric drive power assembly system noise source separation recognition methods |
CN109747575A (en) * | 2018-12-05 | 2019-05-14 | 江苏大学 | A kind of multi-mode in-vehicle sound quality optimization system based on order optimization |
CN110285982A (en) * | 2019-07-31 | 2019-09-27 | 重庆长安汽车股份有限公司 | Evaluation method, device and the controller of the sound quality relative variation of automobile or engine |
CN110427709A (en) * | 2019-08-06 | 2019-11-08 | 华研慧声(苏州)电子科技有限公司 | Sound quality order balance optimizing and active sounding design method based on genetic algorithm |
-
2020
- 2020-07-08 CN CN202010654276.7A patent/CN112304632A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598596A (en) * | 2008-06-06 | 2009-12-09 | 福特环球技术公司 | Analyze the particularly method and apparatus of the noise of vehicle of noise source |
CN107642426A (en) * | 2017-08-31 | 2018-01-30 | 清华大学苏州汽车研究院(相城) | A kind of automobile engine noise initiative control method and system |
CN109194306A (en) * | 2018-08-28 | 2019-01-11 | 重庆长安汽车股份有限公司 | A kind of method and device quantifying automobile noise modulation problems |
CN109211570A (en) * | 2018-11-12 | 2019-01-15 | 吉林大学 | A kind of electric drive power assembly system noise source separation recognition methods |
CN109747575A (en) * | 2018-12-05 | 2019-05-14 | 江苏大学 | A kind of multi-mode in-vehicle sound quality optimization system based on order optimization |
CN110285982A (en) * | 2019-07-31 | 2019-09-27 | 重庆长安汽车股份有限公司 | Evaluation method, device and the controller of the sound quality relative variation of automobile or engine |
CN110427709A (en) * | 2019-08-06 | 2019-11-08 | 华研慧声(苏州)电子科技有限公司 | Sound quality order balance optimizing and active sounding design method based on genetic algorithm |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113758713A (en) * | 2021-08-09 | 2021-12-07 | 重庆长安汽车股份有限公司 | Adaptive rough acoustic frequency band identification method |
CN113758713B (en) * | 2021-08-09 | 2023-06-23 | 重庆长安汽车股份有限公司 | Adaptive recognition method for rough audio frequency band |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1865494B1 (en) | Engine sound processing device | |
CN113758713B (en) | Adaptive recognition method for rough audio frequency band | |
CN108713228B (en) | Vehicle engine sound management | |
CN110803102B (en) | In-vehicle engine order sound analysis method and electric vehicle active sound production system | |
CN106168538B (en) | A kind of ITD of rolling bearing, spectrum kurtosis and smooth iteration envelope Analysis Method | |
CN112304632A (en) | Transient modulation evaluation method for describing human ear perception | |
CN108429995A (en) | Acoustic processing device, sound processing method and storage medium | |
CN112216300A (en) | Noise reduction method and device for sound in driving cab of mixer truck and mixer truck | |
CN106383028A (en) | Gear case fault diagnosis method | |
CN108827452A (en) | A kind of noise rating method of permanent-magnet synchronous hub motor | |
JP2002160512A (en) | Air pressure of tire estimating device | |
CN110767249A (en) | Loudness saliency-based automobile transmission system transient impact noise evaluation method | |
CN110285982B (en) | Method and device for evaluating relative change of sound quality of automobile or engine and controller | |
CN110427709B (en) | Acoustic quality order equalization optimization and active sounding design method based on genetic algorithm | |
CN109194306B (en) | Method and device for quantifying automobile noise modulation problem | |
CN113752943A (en) | Method for controlling tone of electric vehicle based on motor vibration | |
Vold et al. | Vold-Kalman Order Tracking: New Methods for Vehicle Sound Quality and Drive-Train NVH Applications | |
CN106553084A (en) | A kind of lathe flutter on-line monitoring method based on wavelet package transforms and approximate entropy feature | |
US9959852B2 (en) | Vehicle engine sound extraction | |
US6112167A (en) | Tire tread noise treatment | |
CN115662470A (en) | Audio particle extraction method, sound wave synthesis method, device, equipment and medium | |
Liang et al. | A rattle signal denoising and enhancing method based on wavelet packet decomposition and mathematical morphology filter for vehicle | |
CN112380485A (en) | Method for evaluating NVH subjective and objective consistency of new energy automobile electric drive assembly | |
WO2022061709A1 (en) | Electric vehicle sound effect simulation method, device, and system, and automobile | |
CN114446267A (en) | Active sound synthesis method for vehicle |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210202 |
|
RJ01 | Rejection of invention patent application after publication |