WO2014206212A1 - Digital image generation method through 3d spatial distribution of sound quality objective parameters - Google Patents
Digital image generation method through 3d spatial distribution of sound quality objective parameters Download PDFInfo
- Publication number
- WO2014206212A1 WO2014206212A1 PCT/CN2014/079954 CN2014079954W WO2014206212A1 WO 2014206212 A1 WO2014206212 A1 WO 2014206212A1 CN 2014079954 W CN2014079954 W CN 2014079954W WO 2014206212 A1 WO2014206212 A1 WO 2014206212A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- sound
- field
- sound pressure
- dimensional
- loudness
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 238000012800 visualization Methods 0.000 claims abstract description 7
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 20
- 241000282414 Homo sapiens Species 0.000 claims description 10
- 238000001093 holography Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims 1
- 238000010168 coupling process Methods 0.000 claims 1
- 238000005859 coupling reaction Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 230000035807 sensation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 210000003477 cochlea Anatomy 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H3/00—Holographic processes or apparatus using ultrasonic, sonic or infrasonic waves for obtaining holograms; Processes or apparatus for obtaining an optical image from them
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R5/00—Stereophonic arrangements
- H04R5/027—Spatial or constructional arrangements of microphones, e.g. in dummy heads
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/002—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means for representing acoustic field distribution
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/40—Details of arrangements for obtaining desired directional characteristic by combining a number of identical transducers covered by H04R1/40 but not provided for in any of its subgroups
- H04R2201/401—2D or 3D arrays of transducers
Definitions
- the invention relates to near-field acoustic holography technology, sound source recognition and positioning technology, sound field visualization technology, measurement and evaluation technology of objective parameters of sound quality and noise control technology.
- the traditional noise control is mainly aimed at reducing the sound pressure level of the sound field response.
- the research results show that the reduction of sound pressure level can not fully improve the subjective hearing experience of human voice.
- the effect of noise on people is not only related to the sound pressure level, but also to the frequency composition of the sound, the physical characteristics of the human auditory system, and the psychological characteristics.
- the two sounds of the same sound pressure level will cause a huge difference in the psychoacoustic sensitivity of people due to the different frequency composition. Therefore, the evaluation of noise needs to introduce a quality index that reflects the subjective and objective feelings of the person.
- the invention is applied to the visualization of the objective parameter distribution of the sound quality in the three-dimensional sound field and the spatial localization of the sound source which is most closely related to human hearing.
- the existing sound quality analysis and measurement method evaluates the sound quality of the sound in the specified three-dimensional space, only the objective parameter information of the sound quality of a certain measurement position specified in the three-dimensional space can be obtained, and the distribution of the sound quality information in the entire three-dimensional space cannot be obtained.
- visualizing the image it is not only impossible to evaluate the sound quality of the sound field in the entire three-dimensional space, but also to give the spatial position of the sound source related to the sound quality.
- the present invention overcomes the following shortcomings of the prior art: 1. Overcoming the existing optical quality objective parameter measurement, calculation analysis and evaluation methods and techniques alone, can not give the sound quality distribution and the visualization image in the entire three-dimensional space by one measurement, Defects in sound source position information related to objective parameters of sound quality cannot be given. 2. Overcoming the use of near-field acoustic holography to analyze the three-dimensional spatial sound field distribution without considering the subjective auditory sensation of human beings, can not give the defect of the sound source position in the sound field which is closely related to human subjective hearing sensation.
- the invention provides a method for generating a three-dimensional digital image of an acoustic quality objective parameter combined with a near-field acoustic holography method and a sound quality objective parameter analysis method.
- the sound quality objective parameter three-dimensional digital image generating method proposed by the invention has the sound field distribution, and can also give the sound source spatial position information which affects the subjective auditory feeling of the human being, and provides the sound in the whole three-dimensional space.
- the quality of the objective parameter distribution of the image thus providing a visual, direct guidance for noise reduction and sound quality improvement.
- the microphone array can be a rigid surface spherical microphone array, a hollow sphere array, a microphone array conforming to the sound source structure or a planar array.
- a spherical microphone array is arranged in the closed sound field, and sound field sound pressure information is measured and recorded.
- Holographic sound pressure data is obtained from an array of planar microphones and other arrays of conformal microphones of any shape in open and semi-open sound fields.
- Three-dimensional sound field acoustic volume reconstruction According to the measured holographic sound pressure data, the distribution information (sound pressure, normal particle velocity and normal sound intensity, etc.) of the acoustic field of the three-dimensional space is obtained by the near-field acoustic holography method. The form of the three-dimensional image is given.
- n 0
- the 4TM in the formula is determined by the following formula:
- ⁇ is the holographic sound pressure data collected on the microphone array; A (r, A ⁇ iy) is the reconstructed sound pressure at the specified position (r, A) in the three-dimensional space; Y n m ⁇ 0, ( ⁇ % spherical harmonic
- the function, ( ) is the ball Bessel function, is the ball Hank function; represents the conjugate, "," represents the derivative, and N is the number of spherical harmonic expansion items.
- the single point sound quality loudness calculation model in the sound field is as follows: C GE ⁇ + 2E THRQ ) _ ( 2E THRQ ) E THR Q ⁇ E T ⁇ 10
- N/ is the characteristic loudness of the zth filter, which is the audible valve energy level, which is the energy level of the signal, and the unit is dB
- C is the fixed value of 0.046871, when _500, the fixed value is 2.3067
- cochlear The low-frequency gain G is 1, a is 0.2, and when ⁇ 500Hz, E THR Q, « can be obtained according to the discrete data interpolation calculation provided by ANSI.
- G is the cochlear low frequency gain. Therefore the total loudness formula is:
- the coupled three-dimensional matrix mapping model of the sound pressure field and the loudness field in space is established.
- (r, ⁇ ) is the spherical coordinate of the third point in the three-dimensional space, which is the angular frequency.
- M are the number of frequencies corresponding to different sound sources.
- W is an auditory filter matrix consisting of 372 filters W, representing the response of the human ear to all frequencies in the audible band
- N' is the characteristic loudness vector.
- the sum of the points in the characteristic loudness vector N' can be obtained by the formula (4) to obtain the loudness of the specified point of the sound field, and the three-dimensional distribution result of the sound field loudness can be obtained by repeating the calculation process for the three-dimensional spatial node of the sound field.
- Other sound quality objective quantity three-dimensional distribution results can also be obtained using an analysis procedure similar to the loudness calculation, see Figure 1.
- the objective parameter values of the sound quality at various positions in the measured space can be given, and the spatial distribution can be given in the form of a three-dimensional image, thereby identifying the position of the sound source having the greatest influence on the subjective hearing of the human.
- Figure 1 Flow chart for calculating the three-dimensional distribution of objective parameters of sound quality according to the present invention
- Figure 2 Schematic diagram of a spherical microphone array
- FIG. 1 Schematic diagram of spherical microphone array and sound source distribution of dual source
- FIG. 4(d) Comparison of the sharpness and sound pressure calculation results of the sound field of the two-point sound source of 3.5 kHz (70 dB) and 7 kHz (76 dB).
- the present invention will be further described below by way of specific examples. Referring to the drawings:
- the image of the objective parameters of the sound quality in the three-dimensional space provides visual and direct guidance for noise reduction and sound quality improvement.
- the microphone array can be a rigid surface spherical microphone array, a hollow sphere array, a microphone array conforming to the sound source structure or a planar array.
- a spherical microphone array is arranged in the closed sound field, and sound field sound pressure information is measured and recorded.
- Holographic sound pressure data is obtained from an array of planar microphones and other arrays of conformal microphones of any shape in open and semi-open sound fields.
- Three-dimensional sound field acoustic volume reconstruction According to the measured holographic sound pressure data, the distribution information (sound pressure, normal particle velocity and normal sound intensity, etc.) of the acoustic field of the three-dimensional space is obtained by the near-field acoustic holography method. The form of the three-dimensional image is given.
- the holographic sound pressure measured by the spherical microphone array is used to reconstruct the acoustic quantity distribution in the three-dimensional space.
- the sound field transformation formula is as follows:
- the 4TM in the formula is determined by the following formula: p t ( ⁇ , ⁇ , ⁇ ) Y n m ( ⁇ , ⁇ )* sin ⁇
- ⁇ is the holographic sound pressure data collected on the microphone array; A (r, A ⁇ iy) is the reconstructed sound pressure at the specified position (r, A) in the three-dimensional space;
- m ( ) is the spherical harmonic function, ( )
- ( ) is the ball Hank function; represents the conjugate, "," represents the derivative, and N is the number of spherical harmonic expansion items.
- N/ is the characteristic loudness of the Zth filter, which is the audible valve energy level, which is the energy level of the signal, and the unit is dB
- C is the fixed value of 0.046871, when _500, the fixed value is 2.3067
- the cochlea The low-frequency gain G is 1, a is 0.2, and when ⁇ 500Hz, E THR Q, « can be obtained according to the discrete data interpolation calculation provided by ANSI.
- G is the cochlear low frequency gain.
- the total loudness formula is: According to the three-dimensional spatial distribution results of sound pressure in the sound field calculated in 2, combined with the calculation model of the objective parameters of single point sound quality in space, the coupled three-dimensional matrix mapping model of sound pressure field and loudness field in space is established.
- (r, ⁇ ) is the spherical coordinate of the third point in the three-dimensional space, which is the angular frequency
- (o m 2 f m , is the frequency
- m OM M is the number of frequencies corresponding to different sound sources.
- the sound pressure at any frequency at the first point is the vector composed of the sound pressure reconstruction values at the frequencies at the third point in the sound field
- W is the auditory filter matrix composed of 372 filters W, representing the person The response of the ear to all frequencies in the audible frequency band, N' is the characteristic loudness vector.
- the sum of the parts in the characteristic loudness vector N' can be obtained by the formula (4) to obtain the loudness of the specified point of the sound field, and the three-dimensional space node of the sound field Repeating this calculation process can obtain the three-dimensional distribution of sound field loudness.
- Other sound quality objective quantity three-dimensional distribution results can also be obtained by an analysis process similar to loudness calculation, see Figure 1.
- spherical arrays are used for measurement. Array, as shown in Figure 2, on the sphere 1.
- the parameters of sound source 1 are set to lkHz, 75dB, placed 0.3m on the positive X-axis of the space rectangular coordinate system (0.3m, 0, 0);
- the parameter of the sound source 2 is set to 3.5 kHz, 69 dB, and is placed at 0.3 m (-0.3 m, 0, 0) on the X negative half axis of the space rectangular coordinate system, that is, the angle between the two sound sources is 180.
- the radius of the spherical microphone array ( Figure 2) is "0.1m.
- the three-dimensional spatial distribution of sound pressure, loudness and sharpness at a radius of 0.2 m is reconstructed using the calculation method given by the present invention.
- FIG. 4(a) is a comparison diagram of the calculation results of the three-dimensional spatial distribution of sound pressure and loudness when the sound source 1 of 1 kHz, 75 dB, and the sound source 2 of 3.5 kHz and 69 dB are simultaneously present.
- Fig. 4(b) is a comparison diagram of the calculation results of the three-dimensional spatial distribution of the sound pressure and the sharpness when the sound source 1 of the 1 kHz, 75 dB, and the sound source 2 of the 3.5 kHz, 69 dB sound source are simultaneously present.
- the same two-sound source sound field model is used, but the parameter of the sound source 1 is set to 7 kHz, 70 dB, and placed at 0.3 m (0.3 m, 0, 0) on the x positive half-axis of the space rectangular coordinate system;
- the parameter of source 2 is set to 3.5 kHz, 76 dB, and is placed at 0.3 m (-0.3 m, 0, 0) on the x negative half-axis of the space rectangular coordinate system, and the angle between the two sources is still 180°.
- Figure 4 (c) and Figure 4 (d) are 7 kHz, 70 B sound source 1 and 3.5 kHz, 76 dB sound source 2 two sound sources using this method to reconstruct the spherical sound pressure, loudness and A comparison of the results of the three-dimensional spatial distribution of sharpness.
- the method provided by the present invention can obtain the three-dimensional spatial distribution information of the objective parameters of the sound quality of the sound field, and recognize the subjective subject matter.
- the sound source position information closely related to hearing Figure 4 shows the three-dimensional spatial distribution map of the objective parameters of sound quality such as loudness and sharpness, realizing the sound source localization based on human subjective auditory feeling, comparing Fig. 4(a) and Fig. 4(b), figure 4(c) and Fig.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
A digital image generation method through 3D spatial distribution of sound quality objective parameters is performed by the following steps: Step 1. recording holograpahic acoustic pressure data on a hologrpahic measurement surface by using a microphone array; Step 2. reconstructing acoustic quantities in a 3D acoustic field; Step 3. calculating the distribution of a sound quality objective parameter (such as loudness, acuteness, and roughness) at each point in the 3D space according to acoustic quantity information inside the 3D space obtained in Step 2 and based on an acoustic pressure-sound quality objective parameter mapping model, and giving the distribution in the form of a 3D digital image, thereby implementing 3D spatial visualization of sound quality objective parameters.
Description
维空间分布数字图像生成方法 技术领域: Dimensional spatial distribution digital image generation method
本发明涉及近场声全息技术、声源识别定位技术、声场可视化技 术、 声品质客观参量的测量评价技术以及噪声控制技术。 The invention relates to near-field acoustic holography technology, sound source recognition and positioning technology, sound field visualization technology, measurement and evaluation technology of objective parameters of sound quality and noise control technology.
背景技术: Background technique:
传统的噪声控制均以降低声场响应的声压级为主要目标,而研究 结果表明声压级的降低并不能充分改善人对声音的主观听觉感受。噪 声对于人的影响不仅与声压级有关, 还与声音的频率组成、人类听觉 系统的物理特性以及心理特性有关。相同声压级的两种声音由于各自 频率组成的不同, 会导致人在心理感觉上响度的巨大差异, 因此对于 噪声的评价需要引入能够反映人对声音主观和客观感受的量化指标 一声品质, 作为对噪声评价的参考。 The traditional noise control is mainly aimed at reducing the sound pressure level of the sound field response. The research results show that the reduction of sound pressure level can not fully improve the subjective hearing experience of human voice. The effect of noise on people is not only related to the sound pressure level, but also to the frequency composition of the sound, the physical characteristics of the human auditory system, and the psychological characteristics. The two sounds of the same sound pressure level will cause a huge difference in the psychoacoustic sensitivity of people due to the different frequency composition. Therefore, the evaluation of noise needs to introduce a quality index that reflects the subjective and objective feelings of the person. A reference to noise evaluation.
本发明应用于三维声场中声品质客观参量分布的可视化以及与 人听觉感受最密切的声源的空间定位。现有声品质分析和测量方法对 指定三维空间内的声音进行声品质评价时,只能得到三维空间中指定 的某一测量位置的声品质客观参量信息,不能得到整个三维空间内声 品质信息的分布以及可视化图像,因而不仅不能对整个三维空间声场 的声品质好坏进行评价, 更不能给出与声品质相关的声源的空间位 置。 The invention is applied to the visualization of the objective parameter distribution of the sound quality in the three-dimensional sound field and the spatial localization of the sound source which is most closely related to human hearing. When the existing sound quality analysis and measurement method evaluates the sound quality of the sound in the specified three-dimensional space, only the objective parameter information of the sound quality of a certain measurement position specified in the three-dimensional space can be obtained, and the distribution of the sound quality information in the entire three-dimensional space cannot be obtained. As well as visualizing the image, it is not only impossible to evaluate the sound quality of the sound field in the entire three-dimensional space, but also to give the spatial position of the sound source related to the sound quality.
发明内容:
本发明要克服现有技术的下列缺点: 1. 克服单独使用现有声品 质客观参量测量、计算分析和评价方法及技术, 通过一次测量不能给 出整个三维空间内的声品质分布以及可视化图像,也不能给出与声品 质客观参量相关的声源位置信息的缺陷。 2. 克服单独采用近场声全 息方法分析三维空间声场分布时没有考虑人的主观听觉感受因素,不 能给出声场中与人主观听觉感受密切相关的声源位置的缺陷。本发明 提供近场声全息方法与声品质客观参量分析方法相结合的声品质客 观参量三维数字图像生成方法。 本发明提出的声品质客观参量三维数字图像生成方法,其计算结 果在给出声场分布的同时,也能给出影响人主观听觉感受最大的声源 空间位置信息, 并提供整个三维空间内声品质客观参量分布的图像, 从而为降噪和声品质改善提供可视化的、 最直接的指导。 Summary of the invention: The present invention overcomes the following shortcomings of the prior art: 1. Overcoming the existing optical quality objective parameter measurement, calculation analysis and evaluation methods and techniques alone, can not give the sound quality distribution and the visualization image in the entire three-dimensional space by one measurement, Defects in sound source position information related to objective parameters of sound quality cannot be given. 2. Overcoming the use of near-field acoustic holography to analyze the three-dimensional spatial sound field distribution without considering the subjective auditory sensation of human beings, can not give the defect of the sound source position in the sound field which is closely related to human subjective hearing sensation. The invention provides a method for generating a three-dimensional digital image of an acoustic quality objective parameter combined with a near-field acoustic holography method and a sound quality objective parameter analysis method. The sound quality objective parameter three-dimensional digital image generating method proposed by the invention has the sound field distribution, and can also give the sound source spatial position information which affects the subjective auditory feeling of the human being, and provides the sound in the whole three-dimensional space. The quality of the objective parameter distribution of the image, thus providing a visual, direct guidance for noise reduction and sound quality improvement.
本发明按如下歩骤进行: The invention proceeds as follows:
1. 利用传声器阵列 (传声器阵列可以是刚性表面的球形传声器 阵列、 空心球阵列、 与声源结构共形的传声器阵列或平面阵列)记录 全息测量面上的全息声压数据。 在封闭声场布置球形传声器阵列, 测量并记录声场声压信息。在 开放和半开放声场可用平面传声器阵列及其它任意形状共形传声器 阵列获得全息声压数据。 1. Record the holographic sound pressure data on the holographic measurement surface using a microphone array (the microphone array can be a rigid surface spherical microphone array, a hollow sphere array, a microphone array conforming to the sound source structure or a planar array). A spherical microphone array is arranged in the closed sound field, and sound field sound pressure information is measured and recorded. Holographic sound pressure data is obtained from an array of planar microphones and other arrays of conformal microphones of any shape in open and semi-open sound fields.
2. 三维声场声学量重构 根据测量到的全息声压数据,通过近场声全息方法得到三维空间 声场声学量的分布信息 (声压、 法向粒子速度和法向声强等), 并以
三维图像的形式给出。 2. Three-dimensional sound field acoustic volume reconstruction According to the measured holographic sound pressure data, the distribution information (sound pressure, normal particle velocity and normal sound intensity, etc.) of the acoustic field of the three-dimensional space is obtained by the near-field acoustic holography method. The form of the three-dimensional image is given.
采用球形传声器阵列测量的全息声压来重构三维空间的声学: 分布, 声场变换公式如下: n = 0 The holographic sound pressure measured by a spherical microphone array is used to reconstruct the acoustics in three dimensions: distribution, the sound field transformation formula is as follows: n = 0
上式中: (r,^ )为声场中三维空间任意一点的球坐标; α为球形传声 器阵列的半径, k为赚, k = , «为角频率, c为声速, ω 2π;[ ,In the above formula: (r, ^) is the spherical coordinate of any point in the three-dimensional space in the sound field; α is the radius of the spherical microphone array, k is earned, k = , « is the angular frequency, c is the speed of sound, ω 2π;
/为频率。 Α (" 为传声器阵列上采集的全息声压数据; A(r,A^iy)为三维空间指定位置(r,A )处的重构声压; Yn m{0,(^ % 球谐函数, ( )为球贝赛尔函数, 为球汉克函数; 表示 共轭, ","表示导数, N为球谐函数扩展项数。 当指定整个三维空间 中声场重构点的位置, 则获得整个声场的声压分布。 / for the frequency. Α (" is the holographic sound pressure data collected on the microphone array; A (r, A^iy) is the reconstructed sound pressure at the specified position (r, A) in the three-dimensional space; Y n m {0, (^ % spherical harmonic The function, ( ) is the ball Bessel function, is the ball Hank function; represents the conjugate, "," represents the derivative, and N is the number of spherical harmonic expansion items. When specifying the position of the sound field reconstruction point in the entire three-dimensional space, then The sound pressure distribution of the entire sound field is obtained.
3. 根据 2 中得到的三维空间内的声场信息, 计算三维空间各点 声品质客观参量(如响度、 尖锐度、 粗糙度等) 的分布, 并以三维图 像的形式给出, 实现声品质客观参量的三维可视化。 3. According to the sound field information in the three-dimensional space obtained in 2, calculate the distribution of the objective parameters (such as loudness, sharpness, roughness, etc.) of the sound quality of each point in the three-dimensional space, and give them in the form of three-dimensional images to realize the objective quality of sound quality. 3D visualization of parameters.
声场中单点声品质响度计算模型如下:
C GE^ + 2ETHRQ ) _ ( 2ETHRQ ) ETHRQ < ET < 10 The single point sound quality loudness calculation model in the sound field is as follows: C GE^ + 2E THRQ ) _ ( 2E THRQ ) E THR Q < E T < 10
上式中 N/为第 z个滤波器的特征响度, 为可听阀能量级, 为 信号的能量级, 单位均为 dB, C为定值 0.046871, 当_ 500 时, 为定值 2.3067, 耳蜗低频增益 G为 1, a为 0.2, 而当 <500Hz 时, ETHRQ, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
In the above formula, N/ is the characteristic loudness of the zth filter, which is the audible valve energy level, which is the energy level of the signal, and the unit is dB, C is the fixed value of 0.046871, when _500, the fixed value is 2.3067, cochlear The low-frequency gain G is 1, a is 0.2, and when <500Hz, E THR Q, « can be obtained according to the discrete data interpolation calculation provided by ANSI. G is the cochlear low frequency gain. Therefore the total loudness formula is:
根据 2中计算得到的声场中声压的三维空间分布结果,并结合空 间中单点声品质客观参量的计算模型,建立空间中声压场与响度场的 耦合三维矩阵映射模型, According to the three-dimensional spatial distribution results of the sound pressure in the sound field calculated in 2, combined with the calculation model of the single point sound quality objective parameters in the space, the coupled three-dimensional matrix mapping model of the sound pressure field and the loudness field in space is established.
Wl " ^372 W l " ^372
Wl "^372 W l "^372
= \Ν[{τ,θ,φΧ(τ,θ, )入 1x372
或简写为
式中: (r,^ )为三维空间第 点的球坐标, 为角频率,
为频率, mO.M , M为不同声源对应的频率组成个数。 为声 场中第 ζ点处任一频率下的声压, Ρ,为声场中第 ζ点处各频率下的声压
重构值所组成的矢量, W为由 372个滤波器 W组成的听觉滤波器矩阵, 表示人耳对可听频带内所有频率的响应, N'为特征响度矢量。 由公式 (4)对特征响度矢量 N'中的各项求和便可获得声场指定点的响度, 而 对声场三维空间节点重复这一计算过程便可获得声场响度三维分布 结果。其它声品质客观量三维分布结果也可采用与响度计算类似的分 析流程获得, 见附图 1。 = \Ν[{τ,θ,φΧ(τ,θ, ) into 1x372 Or abbreviated as Where: (r,^ ) is the spherical coordinate of the third point in the three-dimensional space, which is the angular frequency. For frequency, mO.M, M are the number of frequencies corresponding to different sound sources. Is the sound pressure at any frequency at the third point in the sound field, Ρ, the sound pressure at each frequency at the third point in the sound field A vector consisting of reconstructed values, W is an auditory filter matrix consisting of 372 filters W, representing the response of the human ear to all frequencies in the audible band, and N' is the characteristic loudness vector. The sum of the points in the characteristic loudness vector N' can be obtained by the formula (4) to obtain the loudness of the specified point of the sound field, and the three-dimensional distribution result of the sound field loudness can be obtained by repeating the calculation process for the three-dimensional spatial node of the sound field. Other sound quality objective quantity three-dimensional distribution results can also be obtained using an analysis procedure similar to the loudness calculation, see Figure 1.
通过上述方法能够给出被测空间内各个位置的声品质客观参量 值, 并以三维图像的形式给出其空间分布, 进而可识别出对人主观听 觉影响最大的声源位置。 附图说明: Through the above method, the objective parameter values of the sound quality at various positions in the measured space can be given, and the spatial distribution can be given in the form of a three-dimensional image, thereby identifying the position of the sound source having the greatest influence on the subjective hearing of the human. BRIEF DESCRIPTION OF THE DRAWINGS:
图 1. 本发明所述的声品质客观参量三维分布的计算流程图 Figure 1. Flow chart for calculating the three-dimensional distribution of objective parameters of sound quality according to the present invention
图 2. 球形传声器阵列示意图 Figure 2. Schematic diagram of a spherical microphone array
图 3. 球形传声器阵列及双声源声场分布示意图 Figure 3. Schematic diagram of spherical microphone array and sound source distribution of dual source
图 4(a). 3.5kHz(69dB)和 lkHz(75dB)两点声源声场的响度和声压计 算结果对比 Figure 4(a). Comparison of the loudness and sound pressure of the two-point source sound field of 3.5kHz (69dB) and lkHz (75dB)
图 4(b). 3.5kHz(69dB)和 lkHz(75dB)两点声源声场的尖锐度和声压计 算结果对比 Figure 4(b). Comparison of the sharpness and sound pressure of the sound field of the two-point sound source of 3.5 kHz (69 dB) and l kHz (75 dB)
图 4(c). 3.5kHz(70dB)和 7kHz(76dB)两点声源声场的响度和声压计算 结果对比 Figure 4(c). Comparison of loudness and sound pressure calculations for sound fields of 3.5 kHz (70 dB) and 7 kHz (76 dB) sound sources
图 4(d). 3.5kHz(70dB)和 7kHz(76dB)两点声源声场的尖锐度和声压计 算结果对比 具体实施方案:
下面通过具体实施例子对本发明作进一歩的描述。 参照附图: 本发明提出的声品质客观参量三维数字图像生成方法,其计算结 果在给出声场分布的同时,也能给出影响人主观听觉感受最大的声源 空间位置信息, 并提供整个三维空间内声品质客观参量分布的图像, 从而为降噪和声品质改善提供可视化的、 最直接的指导。 Figure 4(d). Comparison of the sharpness and sound pressure calculation results of the sound field of the two-point sound source of 3.5 kHz (70 dB) and 7 kHz (76 dB). The present invention will be further described below by way of specific examples. Referring to the drawings: The method for generating a three-dimensional digital image of an objective parameter of sound quality proposed by the present invention, the calculation result of which gives the distribution of the sound field, and also gives the spatial position information of the sound source which affects the subjective auditory feeling of the person, and provides the whole The image of the objective parameters of the sound quality in the three-dimensional space provides visual and direct guidance for noise reduction and sound quality improvement.
本发明按如下歩骤进行: The invention proceeds as follows:
1. 用传声器阵列 (传声器阵列可以是刚性表面的球形传声器 阵列、 空心球阵列、 与声源结构共形的传声器阵列或平面阵列)记录 全息测量面上的全息声压数据。 在封闭声场布置球形传声器阵列, 测量并记录声场声压信息。在 开放和半开放声场可用平面传声器阵列及其它任意形状共形传声器 阵列获得全息声压数据。 1. Record the holographic sound pressure data on the holographic measurement surface using a microphone array (the microphone array can be a rigid surface spherical microphone array, a hollow sphere array, a microphone array conforming to the sound source structure or a planar array). A spherical microphone array is arranged in the closed sound field, and sound field sound pressure information is measured and recorded. Holographic sound pressure data is obtained from an array of planar microphones and other arrays of conformal microphones of any shape in open and semi-open sound fields.
2. 三维声场声学量重构 根据测量到的全息声压数据,通过近场声全息方法得到三维空间 声场声学量的分布信息 (声压、 法向粒子速度和法向声强等), 并以 三维图像的形式给出。 2. Three-dimensional sound field acoustic volume reconstruction According to the measured holographic sound pressure data, the distribution information (sound pressure, normal particle velocity and normal sound intensity, etc.) of the acoustic field of the three-dimensional space is obtained by the near-field acoustic holography method. The form of the three-dimensional image is given.
采用球形传声器阵列测量的全息声压来重构三维空间的声学量 分布, 声场变换公式如下: The holographic sound pressure measured by the spherical microphone array is used to reconstruct the acoustic quantity distribution in the three-dimensional space. The sound field transformation formula is as follows:
Pt (r ^, ∑ P t (r ^, ∑
式中的 4™由有下式确定:
pt (α, θ, φ) Yn m (θ, φ)* sin θάθάφ The 4TM in the formula is determined by the following formula: p t (α, θ, φ) Y n m (θ, φ)* sin θάθάφ
0 JO ^ 0 JO ^
上式中: (r,^ )为声场中三维空间任意一点的球坐标; a为球形传声 器阵列的半径, k为赚, k c, «为角频率, c为声速, ω二 2π;[ ,In the above formula: (r, ^) is the spherical coordinate of any point in the three-dimensional space in the sound field; a is the radius of the spherical microphone array, k is earned, k c, « is the angular frequency, c is the speed of sound, ω 2 2π;
/为频率。 Α (" 为传声器阵列上采集的全息声压数据; A(r,A^iy)为三维空间指定位置(r,A )处的重构声压; „m( )为 球谐函数, ( )为球贝赛尔函数, ( )为球汉克函数; 表示 共轭, ","表示导数, N为球谐函数扩展项数。 当指定整个三维空间 中声场重构点的位置, 则获得整个声场的声压分布。 / for the frequency. Α (" is the holographic sound pressure data collected on the microphone array; A (r, A^iy) is the reconstructed sound pressure at the specified position (r, A) in the three-dimensional space; „ m ( ) is the spherical harmonic function, ( ) For the ball Bessel function, ( ) is the ball Hank function; represents the conjugate, "," represents the derivative, and N is the number of spherical harmonic expansion items. When the position of the sound field reconstruction point in the entire three-dimensional space is specified, the whole is obtained. Sound pressure distribution of the sound field.
3. 根据 2 中得到的三维空间内的声场信息, 计算三维空间各点 声品质客观参量(如响度、 尖锐度、 粗糙度等) 的分布, 并以三维图 像的形式给出, 实现声品质客观参量的三维可视化。 3. According to the sound field information in the three-dimensional space obtained in 2, calculate the distribution of the objective parameters (such as loudness, sharpness, roughness, etc.) of the sound quality of each point in the three-dimensional space, and give them in the form of three-dimensional images to realize the objective quality of sound quality. 3D visualization of parameters.
声场中单点声品质响度计算模型如下: The single point sound quality loudness calculation model in the sound field is as follows:
C GF -\2F E <EI < 10C GF -\2F E <E I < 10
上式中 N/为第 Z个滤波器的特征响度, 为可听阀能量级, 为 信号的能量级, 单位均为 dB, C为定值 0.046871, 当_ 500 时, 为定值 2.3067, 耳蜗低频增益 G为 1, a为 0.2, 而当 <500Hz 时, ETHRQ, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
根据 2中计算得到的声场中声压的三维空间分布结果,并结合空 间中单点声品质客观参量的计算模型,建立空间中声压场与响度场的 耦合三维矩阵映射模型, In the above formula, N/ is the characteristic loudness of the Zth filter, which is the audible valve energy level, which is the energy level of the signal, and the unit is dB, C is the fixed value of 0.046871, when _500, the fixed value is 2.3067, the cochlea The low-frequency gain G is 1, a is 0.2, and when <500Hz, E THR Q, « can be obtained according to the discrete data interpolation calculation provided by ANSI. G is the cochlear low frequency gain. Therefore the total loudness formula is: According to the three-dimensional spatial distribution results of sound pressure in the sound field calculated in 2, combined with the calculation model of the objective parameters of single point sound quality in space, the coupled three-dimensional matrix mapping model of sound pressure field and loudness field in space is established.
或简写为: Or abbreviated as:
式中: (r,^ )为三维空间第 点的球坐标, 为角频率, (om = 2 fm , 为频率, m O.M , M为不同声源对应的频率组成个数。 为声 场中第 ζ点处任一频率下的声压, 为声场中第 ζ点处各频率下的声压 重构值所组成的矢量, W为由 372个滤波器 W组成的听觉滤波器矩阵, 表示人耳对可听频带内所有频率的响应, N'为特征响度矢量。 由公式 (4)对特征响度矢量 N'中的各项求和便可获得声场指定点的响度, 而 对声场三维空间节点重复这一计算过程便可获得声场响度三维分布 结果。其它声品质客观量三维分布结果也可采用与响度计算类似的分 析流程获得, 见附图 1。 本实施例子中, 均以球形阵列作为测量阵, 如图 2所示, 球面上
1. 如图 3所示, 在空间中布置两个脉动球源: 声源 1 的参数设 置为 lkHz、 75dB, 放置在空间直角坐标系的 X 正半轴上 0.3m 处 (0.3m,0,0); 声源 2的参数设置为 3.5kHz、 69dB , 放置在空间直角坐 标系的 X负半轴上 0.3 m处 (-0.3m,0,0), 即两声源间的夹角 为 180° , 球形传声器阵列 (如图 2)的半径《为 0.1m。 采用本发明给出的计算方 法重构半径为 0.2m处的声压、 响度和尖锐度的三维空间分布图。 图 4(a)是 lkHz、 75dB的声源 1和 3.5kHz、 69dB的声源 2两声源同时存 在时的声压和响度的三维空间分布计算结果的比较图。 图 4(b)是 lkHz、 75dB的声源 1和 3.5kHz、 69dB的声源 2两声源同时存在时 的声压和尖锐度的三维空间分布计算结果的比较图。 Where: (r,^) is the spherical coordinate of the third point in the three-dimensional space, which is the angular frequency, (o m = 2 f m , is the frequency, m OM , M is the number of frequencies corresponding to different sound sources. The sound pressure at any frequency at the first point is the vector composed of the sound pressure reconstruction values at the frequencies at the third point in the sound field, and W is the auditory filter matrix composed of 372 filters W, representing the person The response of the ear to all frequencies in the audible frequency band, N' is the characteristic loudness vector. The sum of the parts in the characteristic loudness vector N' can be obtained by the formula (4) to obtain the loudness of the specified point of the sound field, and the three-dimensional space node of the sound field Repeating this calculation process can obtain the three-dimensional distribution of sound field loudness. Other sound quality objective quantity three-dimensional distribution results can also be obtained by an analysis process similar to loudness calculation, see Figure 1. In this example, spherical arrays are used for measurement. Array, as shown in Figure 2, on the sphere 1. As shown in Figure 3, arrange two pulsating ball sources in space: The parameters of sound source 1 are set to lkHz, 75dB, placed 0.3m on the positive X-axis of the space rectangular coordinate system (0.3m, 0, 0); The parameter of the sound source 2 is set to 3.5 kHz, 69 dB, and is placed at 0.3 m (-0.3 m, 0, 0) on the X negative half axis of the space rectangular coordinate system, that is, the angle between the two sound sources is 180. °, the radius of the spherical microphone array (Figure 2) is "0.1m. The three-dimensional spatial distribution of sound pressure, loudness and sharpness at a radius of 0.2 m is reconstructed using the calculation method given by the present invention. Fig. 4(a) is a comparison diagram of the calculation results of the three-dimensional spatial distribution of sound pressure and loudness when the sound source 1 of 1 kHz, 75 dB, and the sound source 2 of 3.5 kHz and 69 dB are simultaneously present. Fig. 4(b) is a comparison diagram of the calculation results of the three-dimensional spatial distribution of the sound pressure and the sharpness when the sound source 1 of the 1 kHz, 75 dB, and the sound source 2 of the 3.5 kHz, 69 dB sound source are simultaneously present.
2. 同样采用 1所述的双声源声场模型, 但声源 1 的参数设置为 7kHz、70dB,放置在空间直角坐标系的 x正半轴上 0.3m处 (0.3m,0,0); 声源 2的参数设置为 3.5kHz、 76dB , 放置在空间直角坐标系的 x负 半轴上 0.3m处 (-0.3m,0,0), 两声源间的夹角 仍为 180°。 图 4(c)和图 4(d)分别是 7kHz, 70B的声源 1和 3.5kHz, 76dB的声源 2两声源采 用本方法重构半径为 0.2m处的球面的声压、 响度和尖锐度的三维空 间分布计算结果的比较图。 2. The same two-sound source sound field model is used, but the parameter of the sound source 1 is set to 7 kHz, 70 dB, and placed at 0.3 m (0.3 m, 0, 0) on the x positive half-axis of the space rectangular coordinate system; The parameter of source 2 is set to 3.5 kHz, 76 dB, and is placed at 0.3 m (-0.3 m, 0, 0) on the x negative half-axis of the space rectangular coordinate system, and the angle between the two sources is still 180°. Figure 4 (c) and Figure 4 (d) are 7 kHz, 70 B sound source 1 and 3.5 kHz, 76 dB sound source 2 two sound sources using this method to reconstruct the spherical sound pressure, loudness and A comparison of the results of the three-dimensional spatial distribution of sharpness.
对比分析上述图给出的结果,与传统根据声压来识别定位声源的 声全息方法不同,本发明提供的方法能够获得声场的声品质客观参量 的三维空间分布信息,并识别出与人主观听觉密切相关的声源位置信 息, 图 4给出了响度和尖锐度等声品质客观参量的三维空间分布图, 实现了依据人的主观听觉感受的声源定位, 对比图 4(a)和图 4(b)、 图
4(c)和图 4(d), 可以发现声压最大的空间位置与响度最大的空间位置 并不相同,而响度最大的空间位置与尖锐度最大的空间位置也并不相 同。 因此, 要根据人主观听觉感受的定位关键声源并采取相应的措施 才能真正实现声场降噪和声品质改善的目的。
Comparing and analyzing the results given in the above figure, different from the traditional acoustic holography method for identifying the localized sound source according to the sound pressure, the method provided by the present invention can obtain the three-dimensional spatial distribution information of the objective parameters of the sound quality of the sound field, and recognize the subjective subject matter. The sound source position information closely related to hearing, Figure 4 shows the three-dimensional spatial distribution map of the objective parameters of sound quality such as loudness and sharpness, realizing the sound source localization based on human subjective auditory feeling, comparing Fig. 4(a) and Fig. 4(b), figure 4(c) and Fig. 4(d), it can be found that the spatial position with the highest sound pressure is not the same as the spatial position with the loudest loudness, and the spatial position with the loudest loudness is not the same as the spatial position with the sharpest sharpness. Therefore, it is necessary to locate the key sound source according to the subjective auditory feeling and take corresponding measures to achieve the purpose of sound field noise reduction and sound quality improvement.
Claims
权 利 要 求 书Claims
歩骤 1. 利用传声器阵列,传声器阵列可以是刚性表面的球形传声 器阵列、 空心球阵列、 与声源结构共形的传声器阵列或平面阵列, 记 录全息测量面上的全息声压数据。 Step 1. Using a microphone array, the microphone array can be a rigid surface spherical microphone array, a hollow sphere array, a microphone array conforming to the sound source structure or a planar array, recording holographic sound pressure data on the holographic measurement surface.
在封闭内声场布置球形传声器阵列, 测量并记录声场全息声压等 声学量信息。 在开放和半开放声场可用平面传声器阵列及其它任意形 状共形传声器阵列获得全息声压数据。 A spherical microphone array is arranged in the closed inner sound field, and acoustic information such as sound field holographic sound pressure is measured and recorded. Holographic sound pressure data is obtained from an array of planar microphones and other arrays of conformal microphones in open and semi-open sound fields.
歩骤 2. 三维声场声学量重构 Step 2. Three-dimensional sound field acoustic volume reconstruction
根据测量到的全息声压数据, 通过近场声全息方法得到三维空间 声场声学量的分布信息 (声压、 法向粒子速度和法向声强等), 并以 三维图像的形式给出; According to the measured holographic sound pressure data, the distribution information (sound pressure, normal particle velocity and normal sound intensity) of the acoustic field of the three-dimensional space is obtained by the near-field acoustic holography method, and is given in the form of a three-dimensional image;
如采用球形传声器阵列测量的全息声压来重构三维空间的声学量 分 The acoustic volume of a three-dimensional space is reconstructed using a holographic sound pressure measured by a spherical microphone array.
式中的 4™由有下式确定: The 4TM in the formula is determined by the following formula:
(2) 上式中: (r, ^ )为声场中三维空间任意一点的球坐标; α为球形传声 器阵列的半径, 为波数, k = 0)lc , «为角频率, c为声速, ω二 2π;[ , (2) In the above formula: (r, ^ ) is the spherical coordinate of any point in the three-dimensional space in the sound field; α is the radius of the spherical microphone array, which is the wave number, k = 0) lc , « is the angular frequency, c is the speed of sound, ω Two 2π;[ ,
/为频率。 Α (" »为传声器阵列上采集的全息声压数据; A(r,A^ iy)为三维空间指定位置 (r,^ )处的重构声压; „m ( )为球 谐函数, ( )为球贝赛尔函数, ( )为球汉克函数; 表示共 轭, ","表示导数, N为球谐函数扩展项数。 当指定整个三维空间中声 场重构点的位置, 则获得整个声场的声压分布; / for the frequency. Α (" » is the holographic sound pressure data collected on the microphone array; A (r, A^ iy) is the reconstructed sound pressure at the specified position (r, ^) in the three-dimensional space; „ m ( ) is the spherical harmonic function, ( ) is the ball Bessel function, ( ) is the ball Hank function; represents the conjugate, "," represents the derivative, and N is the number of spherical harmonic extensions. When the position of the sound field reconstruction point in the entire three-dimensional space is specified, Sound pressure distribution of the entire sound field;
歩骤 Step
3. 根据歩骤 2中得到的三维空间内的声压分布,计算三维空
间各点声品质客观参量的分布, 并以三维图像的形式给出, 实现声品 质客观参量的三维可视化; 3. Calculate the three-dimensional space according to the sound pressure distribution in the three-dimensional space obtained in step 2. The distribution of the objective parameters of the sound quality between the points is given in the form of a three-dimensional image to realize the three-dimensional visualization of the objective parameters of the sound quality;
声场 Sound field
上式中 N,'为第 z个滤波器的特征响度, 为可听阀能量级, 为信 号的能量级, 单位均为 dB, C为定值 0.046871, 当 _ 500Hz时, ETHRQ 为定值 2.3067, 耳蜗低频增益 G为 1, 《为 0.2, 而当 <500HZ 时, E丽 Q, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
In the above formula, N, 'is the characteristic loudness of the z-th filter, which is the audible valve energy level, which is the energy level of the signal, the unit is dB, C is the fixed value 0.046871, when _ 500Hz, E THRQ is the fixed value 2.3067, the cochlear low-frequency gain G is 1, "0.2", and when <500H Z , E-Q, « can be obtained according to the discrete data interpolation calculation provided by ANSI. G is the cochlear low frequency gain. Therefore the total loudness formula is:
根据歩骤 2中计算得到的声场中声压的三维空间分布结果,并 结合空间中单点声品质客观参量的计算模型,建立空间中声压场与 响度场的耦合, 矩阵映射模型为: According to the three-dimensional spatial distribution results of the sound pressure in the sound field calculated in the second step, combined with the calculation model of the objective parameters of the single point sound quality in the space, the coupling between the sound pressure field and the loudness field in the space is established. The matrix mapping model is:
Wl - W372 W l - W 372
Wl - W372 W l - W 372
[ pt {τ,θ,φ,ωχ) -- pt (r, θ, φ. ■ [Ν[(Γ,Θ, )- 2(Γ,Θ, )] 1x372
或简写为: [ p t {τ,θ,φ,ω χ ) -- p t (r, θ, φ. ■ [Ν[(Γ,Θ, )- 2 (Γ,Θ, )] 1x372 Or abbreviated as:
P,W = N' (6) 式中: (r,^ )为三维空间第 点的球坐标, 为角频率, om =2 fm , 为频率, m二 l ..M, M为不同声源对应的频率组成个数。 为声 场中第 z点处指定频率下的声压, 为声场中第 z点处各频率下的声压 重构值所组成的矢量, W为由 372个滤波器 W组成的听觉滤波器矩阵,
表示人耳对可听频带内所有频率的响应, 为特征响度矢量; 由公式P,W = N' (6) where: (r,^) is the spherical coordinate of the third point in space, is the angular frequency, o m =2 f m , is the frequency, m two l..M, M is different The number of frequencies corresponding to the sound source is composed. The sound pressure at the specified frequency at the z-th point in the sound field is the vector composed of the sound pressure reconstruction values at the z-th points in the sound field, and W is the auditory filter matrix composed of 372 filters W. Representing the response of the human ear to all frequencies in the audible band, which is the characteristic loudness vector;
(4)对特征响度矢量 N'中的各项求和便可获得声场指定点的响度,而对 声场三维空间节点重复这一计算过程便可获得声场响度三维分布结 果。
(4) The sum of the elements in the characteristic loudness vector N' can obtain the loudness of the specified point of the sound field, and the calculation process of the three-dimensional space node of the sound field can obtain the three-dimensional distribution result of the sound field loudness.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310261258.2 | 2013-06-26 | ||
CN201310261258.2A CN103389155B (en) | 2013-06-26 | 2013-06-26 | Digital image generation method of three-dimensional spatial distribution of sound quality objective parameters |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014206212A1 true WO2014206212A1 (en) | 2014-12-31 |
Family
ID=49533496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2014/079954 WO2014206212A1 (en) | 2013-06-26 | 2014-06-16 | Digital image generation method through 3d spatial distribution of sound quality objective parameters |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103389155B (en) |
WO (1) | WO2014206212A1 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680376A (en) * | 2017-01-23 | 2017-05-17 | 华南理工大学 | Sound insulation measurement system and measurement method based on three-dimensional sound intensity array |
CN108344503A (en) * | 2018-03-09 | 2018-07-31 | 农业部南京农业机械化研究所 | A kind of high-speed transplanter noise qualities evaluation system |
CN109100685A (en) * | 2018-07-12 | 2018-12-28 | 南京信息工程大学 | A kind of passive acoustic direction blending algorithm of two-sided quaternary cross battle array |
CN109782230A (en) * | 2019-01-21 | 2019-05-21 | 柳州市展虹科技有限公司 | A kind of small-sized acoustical holography measurement of free found field and inverting device |
CN109917338A (en) * | 2019-01-21 | 2019-06-21 | 柳州市展虹科技有限公司 | A kind of small-sized acoustical holography measurement of free found field and inverting device intelligence control system |
CN111812587A (en) * | 2020-07-06 | 2020-10-23 | 上海交通大学 | Sound field test analysis method and system based on machine vision and holographic method |
CN113239573A (en) * | 2021-06-05 | 2021-08-10 | 西北工业大学 | Closed space sound field reconstruction method based on grid fluctuation-free modeling |
CN114386296A (en) * | 2021-11-29 | 2022-04-22 | 哈尔滨工程大学 | Numerical calculation method for three-dimensional sound field in reverberation pool |
CN114383855A (en) * | 2021-11-29 | 2022-04-22 | 江铃汽车股份有限公司 | Method and device for detecting sound quality of electric seat, storage medium and electronic equipment |
CN114543979A (en) * | 2022-02-17 | 2022-05-27 | 浙江工业大学 | Method for predicting far-field acoustic quantity directly radiated by sound source based on near-field acoustic holography in bounded space |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103389155B (en) * | 2013-06-26 | 2015-05-27 | 浙江工业大学 | Digital image generation method of three-dimensional spatial distribution of sound quality objective parameters |
CN103616071B (en) * | 2013-12-09 | 2016-03-09 | 浙江工业大学 | Patch near field acoustic holography-sound quality objective parameter distributed in three dimensions method for visualizing |
CN105785320A (en) * | 2016-04-29 | 2016-07-20 | 重庆大学 | Function type delay summation method for identifying solid sphere array three-dimensional sound source |
CN106124044B (en) * | 2016-06-24 | 2019-05-07 | 重庆大学 | Medicine ball identification of sound source low sidelobe ultrahigh resolution acoustic picture fast acquiring method |
CN106568501B (en) * | 2016-10-25 | 2020-06-23 | 浙江工业大学 | Near-field detection method for sound quality objective parameters of low-noise product |
CN110440906B (en) * | 2018-05-04 | 2020-12-18 | 重庆海扶医疗科技股份有限公司 | Sound field intensity distribution detection method and device of ultrasonic transducer |
JP7205192B2 (en) | 2018-11-22 | 2023-01-17 | 日本電信電話株式会社 | sound pickup device |
CN109668626A (en) * | 2018-12-25 | 2019-04-23 | 东莞材料基因高等理工研究院 | A kind of sound quality evaluation method based on human-computer interaction interface |
CN113515048B (en) * | 2021-08-13 | 2023-04-07 | 华中科技大学 | Method for establishing fuzzy self-adaptive PSO-ELM sound quality prediction model |
CN114577330A (en) * | 2022-01-29 | 2022-06-03 | 中国船舶重工集团公司第七一五研究所 | Underwater far field sound field measurement system based on spherical transducer surface vibration measurement |
CN117407650B (en) * | 2023-12-13 | 2024-04-09 | 中汽研新能源汽车检验中心(天津)有限公司 | Noise quality level evaluation method for driving motor system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1453359A2 (en) * | 2003-02-27 | 2004-09-01 | Modulation Sciences, Inc. | Apparatus for generating and displaying images for determining the quality of audio reproduction |
CN101556187A (en) * | 2009-05-07 | 2009-10-14 | 广东美的电器股份有限公司 | Statistically optimal near-field acoustical holography used for visual recognition of air-conditioner noise sources and operation method thereof |
CN103167373A (en) * | 2011-12-09 | 2013-06-19 | 现代自动车株式会社 | Technique for localizing sound source |
CN103389155A (en) * | 2013-06-26 | 2013-11-13 | 浙江工业大学 | Digital image generation method of three-dimensional spatial distribution of sound quality objective parameters |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183298B (en) * | 2011-03-02 | 2012-12-12 | 浙江工业大学 | Method for separating non-free sound field on irregular single holographic sound pressure measurement plane |
CN102901950B (en) * | 2012-09-20 | 2014-08-20 | 浙江工业大学 | Method for recognizing three-dimensional coordinates of sound sources via planar arrays |
-
2013
- 2013-06-26 CN CN201310261258.2A patent/CN103389155B/en active Active
-
2014
- 2014-06-16 WO PCT/CN2014/079954 patent/WO2014206212A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1453359A2 (en) * | 2003-02-27 | 2004-09-01 | Modulation Sciences, Inc. | Apparatus for generating and displaying images for determining the quality of audio reproduction |
CN101556187A (en) * | 2009-05-07 | 2009-10-14 | 广东美的电器股份有限公司 | Statistically optimal near-field acoustical holography used for visual recognition of air-conditioner noise sources and operation method thereof |
CN103167373A (en) * | 2011-12-09 | 2013-06-19 | 现代自动车株式会社 | Technique for localizing sound source |
CN103389155A (en) * | 2013-06-26 | 2013-11-13 | 浙江工业大学 | Digital image generation method of three-dimensional spatial distribution of sound quality objective parameters |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680376A (en) * | 2017-01-23 | 2017-05-17 | 华南理工大学 | Sound insulation measurement system and measurement method based on three-dimensional sound intensity array |
CN106680376B (en) * | 2017-01-23 | 2023-10-27 | 华南理工大学 | Sound insulation measurement system and method based on three-dimensional sound intensity array |
CN108344503A (en) * | 2018-03-09 | 2018-07-31 | 农业部南京农业机械化研究所 | A kind of high-speed transplanter noise qualities evaluation system |
CN109100685A (en) * | 2018-07-12 | 2018-12-28 | 南京信息工程大学 | A kind of passive acoustic direction blending algorithm of two-sided quaternary cross battle array |
CN109782230B (en) * | 2019-01-21 | 2023-03-31 | 柳州市展虹科技有限公司 | Free sound field small-sized acoustic holographic measurement and inversion device |
CN109782230A (en) * | 2019-01-21 | 2019-05-21 | 柳州市展虹科技有限公司 | A kind of small-sized acoustical holography measurement of free found field and inverting device |
CN109917338A (en) * | 2019-01-21 | 2019-06-21 | 柳州市展虹科技有限公司 | A kind of small-sized acoustical holography measurement of free found field and inverting device intelligence control system |
CN111812587A (en) * | 2020-07-06 | 2020-10-23 | 上海交通大学 | Sound field test analysis method and system based on machine vision and holographic method |
CN113239573A (en) * | 2021-06-05 | 2021-08-10 | 西北工业大学 | Closed space sound field reconstruction method based on grid fluctuation-free modeling |
CN113239573B (en) * | 2021-06-05 | 2024-05-07 | 西北工业大学 | Closed space sound field reconstruction method based on gridless fluctuation modeling |
CN114383855A (en) * | 2021-11-29 | 2022-04-22 | 江铃汽车股份有限公司 | Method and device for detecting sound quality of electric seat, storage medium and electronic equipment |
CN114386296A (en) * | 2021-11-29 | 2022-04-22 | 哈尔滨工程大学 | Numerical calculation method for three-dimensional sound field in reverberation pool |
CN114543979A (en) * | 2022-02-17 | 2022-05-27 | 浙江工业大学 | Method for predicting far-field acoustic quantity directly radiated by sound source based on near-field acoustic holography in bounded space |
CN114543979B (en) * | 2022-02-17 | 2024-05-03 | 浙江工业大学 | Prediction method for sound source direct radiation far-field acoustic quantity based on near-field acoustic holography in bounded space |
Also Published As
Publication number | Publication date |
---|---|
CN103389155B (en) | 2015-05-27 |
CN103389155A (en) | 2013-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2014206212A1 (en) | Digital image generation method through 3d spatial distribution of sound quality objective parameters | |
Ahrens et al. | An analytical approach to sound field reproduction using circular and spherical loudspeaker distributions | |
Pollow et al. | Calculation of head-related transfer functions for arbitrary field points using spherical harmonics decomposition | |
JP2021502015A (en) | Playback system using acoustic holographic recording and metamaterial layers | |
Mehra et al. | Source and listener directivity for interactive wave-based sound propagation | |
JP5024792B2 (en) | Omnidirectional frequency directional acoustic device | |
WO2016065719A1 (en) | Error model-based multi-area sound reproduction method and device | |
Sakamoto et al. | Sound-space recording and binaural presentation system based on a 252-channel microphone array | |
JP2009512364A (en) | Virtual audio simulation | |
CN103616071B (en) | Patch near field acoustic holography-sound quality objective parameter distributed in three dimensions method for visualizing | |
Debertolis et al. | Archaeoacoustics in ancient sites | |
JP7206027B2 (en) | Head-related transfer function learning device and head-related transfer function reasoning device | |
CN107566969A (en) | A kind of enclosed environment internal low-frequency Reconstruction of Sound Field method | |
Bader | Microphone array | |
Nakanishi et al. | Two-dimensional sound field recording with multiple circular microphone arrays considering multiple scattering | |
CN107301153B (en) | Head-related transfer function modeling method based on self-adaptive Fourier decomposition | |
Hiipakka | Estimating pressure at the eardrum for binaural reproduction | |
Rönkkö | Measuring acoustic intensity field in upscaled physical model of ear | |
Wang et al. | Image source method based on the directional impulse responses | |
Sheaffer et al. | A spherical array approach for simulation of binaural impulse responses using the finite difference time domain method | |
Farahikia | Scattering and Fluid-Structure Interaction Analyses Using Computational Acoustics | |
JP7529623B2 (en) | Acoustic material property estimation program, device and method, and acoustic simulation program | |
Douros et al. | Comparison between 2d and 3d models for speech production: a study of french vowels | |
Seo et al. | VR environment-based evaluation of impact factors on the urban soundscape recognition | |
Al-Sheikh et al. | Head related transfer function interpolation based on finite impulse response models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14816728 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 21/06/2016) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 14816728 Country of ref document: EP Kind code of ref document: A1 |