Method for reconstructing an acoustic field

• US 8,731,851 B2
• Filed: 06/26/2009
• Issued: 05/20/2014
• Est. Priority Date: 07/08/2008
• Status: Active Grant

• Current Assignee: Brüel & Kjær Sound & Vibration Measurement A/S (Spectris PLC)
• Original Assignee: Brüel & Kjær Sound & Vibration Measurement A/S (Spectris PLC)
• Inventors: Gomes, Jesper, Hald, Jrgen
• Primary Examiner(s): Breene, John
• Assistant Examiner(s): Zhang, Ruihua
• Application Number:
  US13/002,918
• Publication Number:
  US 20110172936A1
• Time in Patent Office: 1,789 Days
• Field of Search: 702/56, 367/8, 736/46
• US Class Current: 702/56
• CPC Class Codes: G01H 3/125 for representing acoustic field distribution using optical means G01H9/002; sonar systems for imaging G01S7/56, G01S15/89; acoustic holography G03H3/00

Abstract

Disclosed herein is a method of reconstructing a sound field. The method comprises receiving measured values indicative of a first acoustic quantity measured at a set of measurement locations; defining a set of virtual source locations; and computing a second acoustic quantity for at least one target location from one or more wave functions each representative of a respective sound field originating from a respective one of the defined set of virtual source locations; wherein the one or more wave functions are weighted by respective one or more weighting factors, and wherein computing comprises determining the one or more weighting factors from a least-norm fit of the one or more wave functions to the received measured values.

Claims

1. A method of reconstructing a sound field at at least one target location, the sound field being generated by at least one sound source, the method comprising:

receiving, at an interface circuit, one or more measured values of a first acoustic quantity measured at a set of measurement locations;

defining, by a computer, a set of virtual source locations and a set of wave functions each representative of a respective sound field originating from a respective one of the defined set of virtual source locations;

wherein the defining comprises (a) defining a respective scaling surface for each of the defined virtual source locations, the scaling surface being located outside a first source-free region, and (b) scaling each of the set of wave functions to each have a predetermined amplitude on the scaling surface associated with the virtual source location;

computing, by the computer, a second acoustic quantity for at least one target location from a superposition of the set of scaled wave functions multiplied by respective expansion coefficients, the measurement locations and the at least one target location being located in the first source-free region, and the virtual source locations being located outside the first source-free region; and

storing the second acoustic quantity in a memory.

2. A method according to claim 1, wherein computing comprises determining the expansion coefficients from a least-squares fit of the superposition of the set of wave functions to the received measured values.

3. A method according to claim 1, wherein the first source-free region is defined by a first surface;

wherein the at least one sound source is comprised in an object having an object surface; and

wherein at least a part of the object surface defines at least a part of the first surface.

4. A method according to claim 1, wherein the computing the second acoustic quantity comprises

computing each of a set of correlation functions, each correlation function being indicative of a correlation of the set of wave functions at a first one of said measurement locations with the set of wave functions at a second location, the second location being chosen from the target location and the measurement locations;

and

computing the second acoustic quantity from at least the computed correlation functions and the one or more measured values.

5. A method according to claim 4, wherein a first one of the set of correlation functions is a square correlation matrix having a dimensionality equal to the number of measurement locations, and wherein each element of the square correlation matrix is indicative of a correlation of the set of wave functions evaluated at respective measurement locations.

6. A method according to claim 1, wherein each of the one or more wave functions is a solution of the Helmholtz equation.

7. A method according to claim 6, wherein each of the set of wave functions is chosen from a spherical wave function, a multipole wave function, a monopole wave function, a dipole wave function.

8. A method according to claim 3, wherein the object surface defines a first spatial domain including the measurement locations and the target location, and a second spatial domain, different from the first spatial domain, the second spatial domain including the virtual source locations.

9. A method according to claim 8, wherein all virtual source locations have the same distance from the object surface.

10. A method according to claim 1, wherein the measurement locations are arranged in a regular grid.

11. A method according to claim 1, wherein the first acoustic quantity is chosen from a sound pressure, a particle velocity, a sound pressure gradient, and a combination thereof.

12. A method according to claim 1, wherein the second acoustic quantity is chosen from a sound pressure, a particle velocity, and a sound intensity.

16. The method according to claim 1, wherein the total number of wave functions is larger than the number of measurement locations, and wherein the computing comprises determining the expansion coefficients from a least-norm fit of the superposition of the set of wave functions to the received measured values.

17. The method according to claim 1, wherein the first acoustic quantity is the same as the second acoustic quantity.

18. The method according to claim 1, wherein each of the wave functions is scaled to have the predetermined amplitude on the scaling surface at a predetermined distance from the corresponding virtual source location.

13. A processing apparatus for reconstructing a sound field at at least one target location, the sound field being generated by at least one sound source, the processing apparatus comprising:

an interface for receiving measured values of a first acoustic quantity measured at a set of measurement locations; and

a processing unit configured to;

define a set of virtual source locations and a set of wave functions each representative of a respective sound field originating from a respective one of the defined set of virtual source locations;

for each of the defined virtual source locations, define a respective scaling surface located outside a first source-free region;

scale each of the set of wave functions to each have a predetermined amplitude on the scaling surface associated with the virtual source location; and

compute a second acoustic quantity for at least one target location from a superposition of the set of scaled wave functions multiplied by respective expansion coefficients, the measurement locations and the at least one target location being located in the first source-free region, and the virtual source locations being located outside the first source-free region.

14. The processing unit of claim 13 in combination with a set of transducers for measuring a first acoustic quantity at a set of measurement locations and connectable in communication connection to the apparatus so as to forward the measured acoustic quantity to the processing apparatus.

15. The processing apparatus of claim 13 used for reconstructing sound fields generated by noise sources of a vehicle.