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 PDF

Info

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
Application number
PCT/CN2014/079954
Other languages
French (fr)
Chinese (zh)
Inventor
卢奂采
金江明
胡伟杰
Original Assignee
浙江工业大学
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 浙江工业大学 filed Critical 浙江工业大学
Publication of WO2014206212A1 publication Critical patent/WO2014206212A1/en

Links

Classifications

    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H3/00Holographic processes or apparatus using ultrasonic, sonic or infrasonic waves for obtaining holograms; Processes or apparatus for obtaining an optical image from them
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/0443Digital holography, i.e. recording holograms with digital recording means
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/08Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R5/00Stereophonic arrangements
    • H04R5/027Spatial or constructional arrangements of microphones, e.g. in dummy heads
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/002Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means for representing acoustic field distribution
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R2201/00Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
    • H04R2201/40Details 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/4012D 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
Figure imgf000005_0001
Figure imgf000005_0001
式中的 4™由有下式确定:
Figure imgf000005_0002
The 4TM in the formula is determined by the following formula:
Figure imgf000005_0002
上式中: (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'=
Figure imgf000006_0001
N'=
Figure imgf000006_0001
上式中 N/为第 z个滤波器的特征响度, 为可听阀能量级, 为 信号的能量级, 单位均为 dB, C为定值 0.046871, 当_ 500 时, 为定值 2.3067, 耳蜗低频增益 G为 1, a为 0.2, 而当 <500Hz 时, ETHRQ, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
Figure imgf000006_0002
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:
Figure imgf000006_0002
根据 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
Figure imgf000006_0003
或简写为
Figure imgf000006_0004
式中: (r,^ )为三维空间第 点的球坐标, 为角频率,
Figure imgf000006_0005
为频率, mO.M , M为不同声源对应的频率组成个数。 为声 场中第 ζ点处任一频率下的声压, Ρ,为声场中第 ζ点处各频率下的声压 重构值所组成的矢量, W为由 372个滤波器 W组成的听觉滤波器矩阵, 表示人耳对可听频带内所有频率的响应, N'为特征响度矢量。 由公式 (4)对特征响度矢量 N'中的各项求和便可获得声场指定点的响度, 而 对声场三维空间节点重复这一计算过程便可获得声场响度三维分布 结果。其它声品质客观量三维分布结果也可采用与响度计算类似的分 析流程获得, 见附图 1。
= \Ν[{τ,θ,φΧ(τ,θ, ) into 1x372
Figure imgf000006_0003
Or abbreviated as
Figure imgf000006_0004
Where: (r,^ ) is the spherical coordinate of the third point in the three-dimensional space, which is the angular frequency.
Figure imgf000006_0005
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 ^, ∑
n =
Figure imgf000008_0001
n =
Figure imgf000008_0001
式中的 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'=
Figure imgf000009_0001
N'=
Figure imgf000009_0001
上式中 N/为第 Z个滤波器的特征响度, 为可听阀能量级, 为 信号的能量级, 单位均为 dB, C为定值 0.046871, 当_ 500 时, 为定值 2.3067, 耳蜗低频增益 G为 1, a为 0.2, 而当 <500Hz 时, ETHRQ, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
Figure imgf000010_0001
根据 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:
Figure imgf000010_0001
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.
Figure imgf000010_0002
Figure imgf000010_0002
或简写为: Or abbreviated as:
Figure imgf000010_0003
Figure imgf000010_0003
式中: (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
Figure imgf000013_0001
间分布数字图像生成方法,按如下; 进行:
Figure imgf000013_0001
An inter-distributed digital image generation method, as follows;
歩骤 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.
Pt
Figure imgf000013_0002
P t
Figure imgf000013_0002
式中的 4™由有下式确定:  The 4TM in the formula is determined by the following formula:
Figure imgf000013_0003
Figure imgf000013_0003
(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:
Figure imgf000014_0001
N:
Figure imgf000014_0001
上式中 N,'为第 z个滤波器的特征响度, 为可听阀能量级, 为信 号的能量级, 单位均为 dB, C为定值 0.046871, 当 _ 500Hz时, ETHRQ 为定值 2.3067, 耳蜗低频增益 G为 1, 《为 0.2, 而当 <500HZ 时, E丽 Q, «均可根据 ANSI提供的离散数据插值计算获得。 G为耳 蜗低频增益。 因此总响度公式为:
Figure imgf000014_0002
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:
Figure imgf000014_0002
根据歩骤 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
Figure imgf000014_0003
或简写为:
[ p t {τ,θ,φ,ω χ ) -- p t (r, θ, φ. ■ [Ν[(Γ,Θ, )- 2 (Γ,Θ, )] 1x372
Figure imgf000014_0003
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.
PCT/CN2014/079954 2013-06-26 2014-06-16 Digital image generation method through 3d spatial distribution of sound quality objective parameters WO2014206212A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201310261258.2A CN103389155B (en) 2013-06-26 2013-06-26 Digital image generation method of three-dimensional spatial distribution of sound quality objective parameters
CN201310261258.2 2013-06-26

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)

* Cited by examiner, † Cited by third party
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
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

Families Citing this family (11)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
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
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
CN103389155A (en) 2013-11-13
CN103389155B (en) 2015-05-27

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
JP5024792B2 (en) Omnidirectional frequency directional acoustic 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
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
Mehra et al. Wave-based sound propagation for VR applications
Farahikia Scattering and Fluid-Structure Interaction Analyses Using Computational Acoustics
JP7529623B2 (en) Acoustic material property estimation program, device and method, and acoustic simulation program
Jin Spectral analysis and resolving spatial ambiguities in human sound localization
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
Hu Spatial Domain Modeling of Head-Related Transfer Functions

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