CN113108893B - Sound field reconstruction system and method based on sound pressure and particle vibration velocity - Google Patents
Sound field reconstruction system and method based on sound pressure and particle vibration velocity Download PDFInfo
- Publication number
- CN113108893B CN113108893B CN202110302791.3A CN202110302791A CN113108893B CN 113108893 B CN113108893 B CN 113108893B CN 202110302791 A CN202110302791 A CN 202110302791A CN 113108893 B CN113108893 B CN 113108893B
- Authority
- CN
- China
- Prior art keywords
- sound
- vibration velocity
- particle vibration
- sound pressure
- sound field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a sound field reconstruction system and method based on sound pressure and particle vibration velocity, wherein the system comprises: the multi-sensor scanning module is used for scanning by using an array formed by the acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of the measuring surface; the multi-channel real-time data acquisition module is used for carrying out multi-channel real-time acquisition on the sound pressure and particle vibration velocity data based on labview; the filtering module is used for filtering data points with amplitude values smaller than a preset value on a measuring surface by using the amplitude filtering coefficient of the sound field measuring surface; and the sound field reconstruction module is used for performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound field reconstruction of a sound intensity field. According to the method, the particle vibration velocity is introduced to reconstruct the sound field on the basis of the original sound pressure-based equivalent source method near-field acoustic holography, so that the precision of sound field reconstruction and the positioning precision of a sound source are effectively improved.
Description
Technical Field
The invention relates to the technical field of sound field reconstruction, in particular to a sound field reconstruction system and method based on sound pressure and particle vibration velocity.
Background
The existing noise monitoring is usually carried out by means of measuring instruments such as a sound level meter, a frequency analyzer and the like, an obtained result only contains one aspect of sound information, and a large amount of noise information cannot be effectively obtained. The sound field distribution condition in a certain space can be obtained through the acoustic imaging technology, the noise source can be accurately and visually positioned, and the running conditions of corresponding equipment and systems can be monitored.
Existing acoustic imaging techniques are mainly implemented by near field acoustic holography (NAH). The NAH technology requires that the distance between the holographic measuring surface and the sound source surface is a fraction of the wavelength, data are collected through a near field, acoustic physical quantities such as sound pressure, vibration velocity, sound intensity and sound power in the sound source and the sound field are reconstructed through a space sound field transformation algorithm, and the sound field can be predicted by comparing the position far away from the holographic surface. Compared with the traditional acoustic holography technology, the NAH technology not only records propagation wave components, but also records evanescent wave components, so that the NAH can reconstruct a sound field on the premise of not being limited by Rayleigh resolution criterion. The reconstruction of various acoustic quantities of a sound source and any point of the whole sound field can be realized only by measuring the complex sound pressure or the vibration velocity of the holographic surface by the NAH, and the NAH technology adopts non-contact measurement, so that the method is suitable for environments which cannot directly measure vibration such as high temperature, high heat and the like. The NAH technology is a powerful tool for sound source identification and sound field visualization, has wide application prospects in mechanical fault diagnosis and mechanical equipment noise measurement and control engineering, and is a powerful tool for engineering technicians to solve the problems of environment and product noise. By virtue of the advantages, the NAH is valued by numerous acoustic researchers and organizations at home and abroad since its birth, and through the development of decades, the NAH has made important progress in the aspects of theory, reconstruction algorithm and application.
Currently, the mainstream near-field acoustic holography carries out sound field reconstruction and prediction by measuring sound pressure, and as is well known, particle vibration velocity is also worth paying attention as one of acoustic quantities. Especially, with the improvement of the method and the requirement of practical engineering, the sound field reconstruction method based on the sound pressure alone cannot meet the requirement of people on the reconstruction accuracy.
Disclosure of Invention
The invention aims to provide a sound field reconstruction system and method based on sound pressure and particle vibration velocity, which introduces particle vibration velocity quantity to reconstruct a sound field on the basis of the original sound pressure-based equivalent source method near-field acoustic holography, constructs a set of sound field reconstruction system including data acquisition, storage and processing, greatly improves the precision of sound field reconstruction and the positioning precision of a sound source, and proves the beneficial effect of the particle vibration velocity on the aspect of sound field reconstruction.
To solve the above technical problem, an embodiment of the present invention provides the following solutions:
in one aspect, a sound field reconstruction system based on sound pressure and particle vibration velocity is provided, including:
the multi-sensor scanning module is used for scanning by using an array formed by the acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of the measuring surface;
the multichannel real-time data acquisition module is used for carrying out multichannel real-time acquisition on the sound pressure and particle vibration velocity data based on Iabview;
the filtering module is used for filtering data points with amplitude values smaller than a preset value on a measuring surface by using the amplitude filtering coefficient of the sound field measuring surface;
and the sound field reconstruction module is used for performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound field reconstruction of a sound intensity field.
Preferably, the multi-sensor scanning module comprises a plurality of acoustic vector sensors, a sensor support, a reference microphone, an XY linear module and a scanning stepping controller;
the acoustic vector sensors are used for simultaneously measuring sound pressure and particle vibration velocity data, the sensor support is used for accommodating and fixing a plurality of acoustic vector sensors, the sensor support is connected with the XY linear module, the scanning stepper controller controls the sensor support to move up and down for scanning, and the scanning stepper controller controls the sensor array to move back and forth to adjust the distance between the sensor array and equipment.
Preferably, the amplitude filter coefficient C p And C u Is defined as:
wherein p is max And u max Respectively the maximum sound pressure and particle vibration speed on the measuring surface, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtration, wherein C p And C u Take 15 dB.
Preferably, the sound field reconstruction module is specifically configured to:
aiming at different characteristics of sound pressure and particle vibration velocity, a Tikhonov regularization method is used for solving the problem of uncertainty generated in the matrix inversion process, a Bayesian method is used for solving regularization parameters in sound field reconstruction of the sound pressure, broadband sound holography based on sound pressure differential filtering is used for solving sound field reconstruction of the particle vibration velocity, and a sound intensity field is obtained by calculating the reconstructed sound pressure and the particle vibration velocity.
Preferably, the Bayesian method obtains the optimal regularization parameter by minimizing the following function:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let theta be lambda 2 The derivative of J (theta) with respect to theta is obtained:
combining the above formulas, and solving for J (lambda) by gradient descent method 2 ) Finally, the regularization parameter lambda is obtained.
Preferably, the sound field reconstruction module is further specifically configured to:
reconstructing sound pressure P near an equivalent source surface using Bayesian regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, and P is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing broadband acoustic holography, and the equivalent source intensity q is equivalent by utilizing an equivalent source intensity filter gamma for each iteration preset step number k And (3) carrying out filtering treatment:
q k+1 ＝q k ·γ
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
On one hand, the sound field reconstruction method based on the sound pressure and the particle vibration velocity is provided, and comprises the following steps:
scanning by using an array formed by acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of a measuring surface;
multi-channel real-time acquisition is carried out on sound pressure and particle vibration velocity data based on labview;
filtering data points with amplitude values smaller than a preset value on a measuring surface by using an amplitude value filtering coefficient of the sound field measuring surface;
and performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound field reconstruction of a sound intensity field.
Preferably, the amplitude filtering coefficient C p And C u Is defined as:
wherein p is max And u max Respectively the maximum sound pressure and particle vibration speed on the measuring surface, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtration, wherein C p And C u Take 15 dB.
Preferably, aiming at different characteristics of sound pressure and particle vibration velocity, a Tikhonov regularization method is used for solving the problem of uncertainty generated in the matrix inversion process, a Bayesian method is used for solving regularization parameters in sound field reconstruction of the sound pressure, broadband sound holography based on sound pressure differential filtering is used for solving sound field reconstruction of the particle vibration velocity, and a sound intensity field is obtained by calculation of the reconstructed sound pressure and the particle vibration velocity;
the Bayesian method obtains optimal regularization parameters by minimizing the following functions:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let theta be lambda 2 The derivative of J (θ two) with respect to θ is obtained:
combining the above formulas, and solving for J (lambda) by gradient descent method 2 ) Finally, the regularization parameter lambda is obtained.
Preferably, the sound field reconstruction step comprises:
reconstructing sound pressure P near an equivalent source surface using Bayesian regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, rho is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing broadband acoustic holography, and the equivalent source intensity q is equivalent by utilizing an equivalent source intensity filter gamma for each iteration preset step number k And (3) carrying out filtering treatment:
q k+1 ＝q k ·γ
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
the invention provides a sound field reconstruction system and a sound field reconstruction method for near-field acoustic imaging, which are used for obtaining complete measurement surface data by scanning with fewer sensors; the amplitude filtering coefficient of the measuring surface is provided, and adverse effects caused by measuring errors and sound source instability are eliminated; the particle vibration velocity is added into the calculation of sound field reconstruction, so that the reconstruction precision is improved; the near-field acoustic holography based on the equivalent source method is improved and applied to the calculation of the mass point vibration velocity, the problem of uncertainty in the matrix inversion process is solved by using a standard Tikhonov regularization method, parameters are solved by using a Bayes method in the regularization process of the sound pressure quantity, the equivalent source strength of the mass point vibration velocity quantity is solved by using broadband acoustic holography based on sound pressure differential filtering, and the positioning accuracy of a sound source is improved; finally, the sound pressure quantity and the mass point vibration velocity quantity are combined to obtain a sound intensity field, so that the reconstruction effect is better.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic structural diagram of a sound field reconstruction system based on sound pressure and particle vibration velocity according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a multi-sensor scanning module provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a calculation process of a sound field reconstruction module according to an embodiment of the present invention;
FIG. 4 is a flowchart of a sound field reconstruction method based on sound pressure and particle vibration velocity according to an embodiment of the present invention;
fig. 5a to 5f are schematic diagrams of measurement results and reconstruction results of sound pressure, particle vibration velocity and sound intensity of a dual sound source of an anechoic chamber according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
An embodiment of the present invention first provides a sound field reconstruction system based on sound pressure and particle velocity, as shown in fig. 1, the system includes:
the multi-sensor scanning module is used for scanning by using an array formed by the acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of the measuring surface;
the multi-channel real-time data acquisition module is used for carrying out multi-channel real-time acquisition on the sound pressure and particle vibration velocity data based on labview;
the filtering module is used for filtering data points with amplitude values smaller than a preset value on a measuring surface by using the amplitude filtering coefficient of the sound field measuring surface;
and the sound field reconstruction module is used for performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound field reconstruction of a sound intensity field.
According to the method, the particle vibration velocity is introduced to carry out sound field reconstruction on the basis of the original near-field acoustic holography based on the equivalent source method of sound pressure, a set of sound field reconstruction system including data acquisition, storage and processing is built, the precision of sound field reconstruction and the positioning precision of a sound source are greatly improved, and the beneficial effect of the particle vibration velocity on the aspect of sound field reconstruction is verified.
Furthermore, the multi-sensor scanning module comprises a plurality of acoustic vector sensors, a sensor bracket, a reference microphone, an XY linear module and a scanning stepping controller;
the acoustic vector sensors are used for simultaneously measuring sound pressure and particle vibration velocity data, the sensor support is used for accommodating and fixing the acoustic vector sensors, the sensor support is connected with the XY linear module, the scanning stepper controller controls the sensor support to move up and down for scanning, and the scanning stepper controller controls the sensor array to move back and forth to adjust the distance between the sensor array and equipment.
FIG. 2 presents a schematic view of a multi-sensor scanning module in accordance with an embodiment of the present invention. Preferably, the transducer holder is about 1.7m long and can accommodate up to 29 acoustic vector transducers (not shown); the sensor support is connected with the XY linear module, and the scanning precision can reach 0.1 mm.
Further, the amplitude filtering coefficient C p And C u Is defined as:
wherein p is max And u max Respectively the maximum sound pressure and particle vibration speed on the measuring surface, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtration, wherein C p And C u Take 15 dB.
Further, the sound field reconstruction module is specifically configured to:
aiming at different characteristics of sound pressure and particle vibration velocity, a Tikhonov regularization method is used for solving the problem of uncertainty generated in the matrix inversion process, a Bayesian method is used for solving regularization parameters in sound field reconstruction of the sound pressure, broadband sound holography based on sound pressure differential filtering is used for solving sound field reconstruction of the particle vibration velocity, and a sound intensity field is obtained by calculating the reconstructed sound pressure and the particle vibration velocity. The specific calculation process of the sound field reconstruction module is shown in fig. 3.
The Bayesian method obtains an optimal regularization parameter by minimizing the following function:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let theta be lambda 2 The derivative of J (theta) with respect to theta is obtained:
the gradient descent method can be used to solve J (. lamda.) in the formula (3) in the combination formula (4) 2 ) Finally, the regularization parameter lambda is obtained.
Further, the sound pressure P near the equivalent source surface is reconstructed by Bayes regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, rho is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing broadband acoustic holography, and the equivalent source intensity q is equivalent by utilizing an equivalent source intensity filter gamma for each iteration preset step number k And (3) carrying out filtering treatment:
q k+1 ＝q k ·γ (7)
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
Accordingly, an embodiment of the present invention further provides a sound field reconstruction method based on sound pressure and particle velocity, as shown in fig. 4, the method includes the following steps:
scanning by using an array formed by acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of a measuring surface;
multi-channel real-time acquisition is carried out on sound pressure and particle vibration velocity data based on labview;
filtering data points with amplitude values smaller than a preset value on a measuring surface by using an amplitude value filtering coefficient of the sound field measuring surface;
and performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound field reconstruction of a sound intensity field.
Further, the amplitude filtering coefficient C p And C u Is defined as:
wherein p is max And u max Respectively the maximum sound pressure and particle vibration speed on the measuring surface, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtration, wherein C p And C u Take 15 dB.
Further, aiming at different characteristics of sound pressure and particle vibration velocity, a Tikhonov regularization method is used for solving the problem of uncertainty generated in the matrix inversion process, a Bayes method is adopted for regularization parameter calculation in sound field reconstruction of the sound pressure, broadband sound holography based on sound pressure differential filtering is used for solving in sound field reconstruction of the particle vibration velocity, and a sound intensity field is obtained by calculation of the reconstructed sound pressure and the particle vibration velocity;
the Bayesian method obtains the optimal regularization parameter by minimizing the following function:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let θ equal λ 2 The derivative of J (theta) with respect to theta is obtained:
the gradient descent method can be used to solve J (. lamda.) in the formula (3) in the combination formula (4) 2 ) Finally, the regularization parameter lambda is obtained.
Further, the sound field reconstruction step includes:
reconstructing sound pressure P near an equivalent source surface using Bayesian regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, rho is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing broadband acoustic holography, and the equivalent source intensity q is equivalent by utilizing an equivalent source intensity filter gamma for each iteration preset step number k And (3) carrying out filtering treatment:
q k+1 ＝q k ·γ (7)
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
Taking a sound-deadening chamber dual-sound source as an example, sound field reconstruction is performed by using the method of the invention, and fig. 5 a-5 f are schematic diagrams of measurement results and reconstruction results of sound pressure, particle vibration velocity and sound intensity of the sound-deadening chamber dual-sound source. The method can effectively improve the accuracy of sound field reconstruction and the positioning accuracy of the sound source, combines the sound pressure quantity and the particle vibration velocity quantity to obtain the sound intensity field, and has better reconstruction effect.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (2)
1. A sound field reconstruction system based on sound pressure and particle vibration velocity, comprising:
the multi-sensor scanning module is used for scanning by using an array formed by the acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of the measuring surface;
the multi-sensor scanning module comprises a plurality of acoustic vector sensors, a sensor bracket, a reference microphone, an XY linear module and a scanning stepping controller;
the acoustic vector sensors are used for simultaneously measuring sound pressure and particle vibration velocity data, the sensor support is used for accommodating and fixing a plurality of acoustic vector sensors, the sensor support is connected with the XY linear module, the scanning stepping controller controls the sensor support to move up and down for scanning, and the scanning stepping controller controls the sensor array to move back and forth to adjust the distance between the sensor array and equipment;
the multi-channel real-time data acquisition module is used for carrying out multi-channel real-time acquisition on the sound pressure and particle vibration velocity data based on labview;
the filtering module is used for filtering data points with amplitude values smaller than a preset value on a measuring surface by using the amplitude filtering coefficient of the sound field measuring surface;
the amplitude filter coefficient C p And C u Is defined as:
wherein p is max And u max Respectively, maximum on the measuring planeSound pressure and particle vibration velocity, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtration, wherein C p And C u Taking 15 dB;
the sound field reconstruction module is used for performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound intensity field sound field reconstruction;
the sound field reconstruction module is specifically configured to:
aiming at different characteristics of sound pressure and particle vibration velocity, solving the problem of uncertainty generated in the matrix inversion process by using a Tikhonov regularization method, solving regularization parameters by using a Bayesian method for sound field reconstruction of the sound pressure, solving the sound field reconstruction of the particle vibration velocity by using broadband sound holography based on sound pressure differential filtering, and calculating a sound intensity field by using the reconstructed sound pressure and the particle vibration velocity;
the Bayesian method obtains optimal regularization parameters by minimizing the following functions:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let theta be lambda 2 The derivative of J (theta) with respect to theta is obtained:
combining the above formulas, and solving for J (lambda) by gradient descent method 2 ) Finally obtaining a regularization parameter lambda;
the sound field reconstruction module is further specifically configured to:
reconstructing sound pressure P near an equivalent source surface using Bayesian regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, and P is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing the broadband acoustic holography, and each iteration of a preset step number carries out filtering processing on the equivalent source intensity qk by utilizing an equivalent source intensity filter gamma:
q k+1 ＝q k ·γ
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
2. A sound field reconstruction method based on sound pressure and particle vibration velocity is characterized by comprising the following steps:
scanning by using an array formed by acoustic vector sensors to obtain the sound pressure and particle vibration velocity data of a measuring surface;
the method comprises the following steps that a multi-sensor scanning module is based on, wherein the multi-sensor scanning module comprises a plurality of acoustic vector sensors, a sensor bracket, a reference microphone, an XY linear module and a scanning stepping controller;
the acoustic vector sensors are used for simultaneously measuring sound pressure and particle vibration velocity data, the sensor support is used for accommodating and fixing a plurality of acoustic vector sensors, the sensor support is connected with the XY linear module, the scanning stepping controller controls the sensor support to move up and down for scanning, and the scanning stepping controller controls the sensor array to move back and forth to adjust the distance between the sensor array and equipment;
multi-channel real-time acquisition is carried out on sound pressure and particle vibration velocity data based on labview;
filtering data points with amplitude values smaller than a preset value on a measuring surface by using an amplitude value filtering coefficient of the sound field measuring surface;
the amplitude filter coefficient C p And C u Is defined as:
wherein p is max And u max Respectively the maximum sound pressure and particle vibration speed on the measuring surface, p cutoff And u cutoff Cutting off sound pressure and particle vibration velocity during filtering, wherein C p And C u Taking 15 dB;
performing sound field reconstruction by using a near-field acoustic holography method based on an equivalent source method, wherein the sound field reconstruction comprises sound pressure sound field reconstruction, particle vibration velocity sound field reconstruction and sound intensity field sound field reconstruction;
aiming at different characteristics of sound pressure and particle vibration velocity, solving the problem of uncertainty generated in the matrix inversion process by using a Tikhonov regularization method, solving regularization parameters by using a Bayesian method for sound field reconstruction of the sound pressure, solving the sound field reconstruction of the particle vibration velocity by using broadband sound holography based on sound pressure differential filtering, and calculating a sound intensity field by using the reconstructed sound pressure and the particle vibration velocity;
the Bayesian method obtains optimal regularization parameters by minimizing the following functions:
wherein σ k Is the k-th singular value, y, of the transfer function matrix G k Is the kth element of the vector y, y ═ U H P and U are left singular matrixes obtained by SVD decomposition of the transfer function matrix, and P is sound pressure on the measuring surface;
let theta be lambda 2 The derivative of J (theta) with respect to theta is obtained:
combining the above formulas, and solving for J (lambda) by gradient descent method 2 ) Finally obtaining a regularization parameter lambda;
the sound field reconstruction step comprises:
reconstructing sound pressure P near an equivalent source surface using Bayesian regularization 1 And P 2 ，P 1 And P 2 The sound pressure on two planes separated by delta is the vibration velocity V of mass points near an equivalent source according to a wave equation 1 From P 1 And P 2 The difference approximation of (d) yields:
wherein, omega is wave number, rho is air density; when Δ is sufficiently small, V 1 Equivalent to the real particle vibration speed; then, for V 1 Normalization processing is performed, and the normalized values are used as a filter gamma of particle vibration velocity equivalent source intensity:
on the basis, iterative solution is carried out on the particle vibration velocity equivalent source intensity by utilizing the broadband acoustic holography, and the equivalent source intensity q is obtained by utilizing an equivalent source intensity filter gamma for each iteration preset step number k And (3) carrying out filtering treatment:
q k+1 ＝q k ·γ
and finally, realizing the reconstruction of the particle vibration velocity by using the particle vibration velocity equivalent source strength q meeting the convergence condition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110302791.3A CN113108893B (en) | 2021-03-22 | 2021-03-22 | Sound field reconstruction system and method based on sound pressure and particle vibration velocity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110302791.3A CN113108893B (en) | 2021-03-22 | 2021-03-22 | Sound field reconstruction system and method based on sound pressure and particle vibration velocity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113108893A CN113108893A (en) | 2021-07-13 |
| CN113108893B true CN113108893B (en) | 2022-08-16 |
Family
ID=76710329
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110302791.3A Active CN113108893B (en) | 2021-03-22 | 2021-03-22 | Sound field reconstruction system and method based on sound pressure and particle vibration velocity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113108893B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114136648B (en) * | 2021-10-20 | 2023-06-13 | 中国航发四川燃气涡轮研究院 | Pneumatic excitation identification method for aeroengine fan movable blade based on acoustic array |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4905202A (en) * | 1987-12-28 | 1990-02-27 | Richard L. Scully | Detection panel and method for acoustic holography |
| CN1370310A (en) * | 1999-08-18 | 2002-09-18 | 哈茉尼彩色音技术开发公司 | Method and appats. for electing harmonic color using harmonics, and method and appts. for converting sound to color or color to sound |
| CN101566496A (en) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | Method for sound field separation by double plane vibration speed measurement and equivalent source method |
| CN102089633A (en) * | 2008-07-08 | 2011-06-08 | 布鲁尔及凯尔声音及振动测量公司 | Method for reconstructing an acoustic field |
| CN104008287A (en) * | 2014-05-23 | 2014-08-27 | 北京科技大学 | Sound field reconstructing and ghost image restraining method based on PSO-MVDR |
| JP2019030623A (en) * | 2016-11-15 | 2019-02-28 | 炭 親良 | Beamforming method, measurement imaging device, and communication device |
| CN112437392A (en) * | 2020-12-10 | 2021-03-02 | 科大讯飞(苏州)科技有限公司 | Sound field reconstruction method and device, electronic equipment and storage medium |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6615143B2 (en) * | 2001-03-01 | 2003-09-02 | Wayne State University | Method and apparatus for reconstructing and acoustic field |
| EP1682856B1 (en) * | 2003-11-10 | 2014-01-08 | Brüel & Kjaer Sound & Vibration Measurement A/S | A method of determining the sound pressure resulting from a surface element of a sound emitting surface |
| CN102680071B (en) * | 2012-05-21 | 2013-11-20 | 哈尔滨工程大学 | Noise source identification method adopting vibration speed measurement and partial near-field acoustical holography method |
| CN102901559B (en) * | 2012-09-27 | 2014-08-06 | 哈尔滨工程大学 | Sound field separating method based on single-surface measurement and local acoustical holography method |
| CN107389789B (en) * | 2017-03-03 | 2020-07-07 | 南京航空航天大学 | Ultrasonic tomography method based on Bayesian regularization |
| CN111664932A (en) * | 2020-05-22 | 2020-09-15 | 重庆大学 | Sound source identification method based on Bayesian compressed sensing |
-
2021
- 2021-03-22 CN CN202110302791.3A patent/CN113108893B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4905202A (en) * | 1987-12-28 | 1990-02-27 | Richard L. Scully | Detection panel and method for acoustic holography |
| CN1370310A (en) * | 1999-08-18 | 2002-09-18 | 哈茉尼彩色音技术开发公司 | Method and appats. for electing harmonic color using harmonics, and method and appts. for converting sound to color or color to sound |
| CN102089633A (en) * | 2008-07-08 | 2011-06-08 | 布鲁尔及凯尔声音及振动测量公司 | Method for reconstructing an acoustic field |
| CN101566496A (en) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | Method for sound field separation by double plane vibration speed measurement and equivalent source method |
| CN104008287A (en) * | 2014-05-23 | 2014-08-27 | 北京科技大学 | Sound field reconstructing and ghost image restraining method based on PSO-MVDR |
| JP2019030623A (en) * | 2016-11-15 | 2019-02-28 | 炭 親良 | Beamforming method, measurement imaging device, and communication device |
| CN112437392A (en) * | 2020-12-10 | 2021-03-02 | 科大讯飞(苏州)科技有限公司 | Sound field reconstruction method and device, electronic equipment and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| Experimental reconstruction of non-stationary sound and vibration sources by means of Transient Planar Near-field Acoustic Holography;van Dalen, PL;《PROCEEDINGS OF INTERNATIONAL CONFERENCE ON NOISE AND VIBRATION ENGINEERING (ISMA2012) / INTERNATIONAL CONFERENCE ON UNCERTAINTY IN STRUCTURAL DYNAMICS (USD2012)》;20120919;3631-3643 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113108893A (en) | 2021-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109885945B (en) | Boundary element method near-field acoustic holographic transformation method in half-space environment | |
| CN102089633B (en) | Method for reconstructing an acoustic field | |
| CN103743469B (en) | Piezoelectric underwater sound transducer acoustic radiation mode measurement method and system | |
| CN109444861B (en) | Plane sonar array impedance characteristic near-field acoustic holography calibration and measurement method | |
| CN113063490B (en) | Sound field separation method based on sound pressure and particle vibration velocity double-sided measurement | |
| CN113848252B (en) | Corrosion imaging method based on ultrasonic guided wave and convolutional neural network | |
| CN110487393A (en) | The unstable state free field restoring method measured using single side acoustic pressure and particle velocity | |
| CN113687307B (en) | Self-adaptive wave beam forming method under low signal-to-noise ratio and reverberation environment | |
| CN113108893B (en) | Sound field reconstruction system and method based on sound pressure and particle vibration velocity | |
| CN108318123A (en) | A kind of near field sound chromatography test method of underwater acoustic array mutual radiation impedance | |
| CN113671041A (en) | Sound insulation measuring system based on near-field acoustic holography and measuring method thereof | |
| CN109375198B (en) | Low-frequency sonar array impedance characteristic calibration method | |
| CN109798975B (en) | Free field real-time reduction method of unsteady-state planar sound source by adopting sound pressure and particle acceleration measurement | |
| CN114543979B (en) | Prediction method for sound source direct radiation far-field acoustic quantity based on near-field acoustic holography in bounded space | |
| Damani et al. | The Characteristics of Sub-Convective Wall Pressure Fluctuations in a Turbulent Boundary Layer | |
| CN114630260B (en) | Method for reconstructing sound field in semi-space uniform flow | |
| CN117849683A (en) | Magnetic resonance spectrum reconstruction method based on hypercomplex neural network | |
| CN111257833A (en) | Sound source identification method based on Laplace norm for fast iterative shrinkage threshold | |
| Mosland | Reciprocity calibration method for ultrasonic piezoelectric transducers in air | |
| CN110161459B (en) | Rapid positioning method for amplitude modulation sound source | |
| Feng et al. | An Equivalent Source Method Based on Sound Pressure Differential Filter for Particle Velocity Reconstruction | |
| Okubo et al. | Three-Dimensional Sound Wave Propagation Reproduction by CE-FDTD Simulation Applying Actual Radiation Characteristics | |
| CN109270492B (en) | Regularization parameter selection method for large holographic distance | |
| CN115782480B (en) | Rolling tire structure noise time domain prediction method | |
| Huang et al. | Theory and methodology of electromagnetic ultrasonic guided wave imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |