[go: up one dir, main page]

CN111561991B - A near-field acoustic holography method based on edge filling and Fourier transform - Google Patents

A near-field acoustic holography method based on edge filling and Fourier transform Download PDF

Info

Publication number
CN111561991B
CN111561991B CN202010510380.9A CN202010510380A CN111561991B CN 111561991 B CN111561991 B CN 111561991B CN 202010510380 A CN202010510380 A CN 202010510380A CN 111561991 B CN111561991 B CN 111561991B
Authority
CN
China
Prior art keywords
virtual
holographic surface
holographic
sound
fourier transform
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
Application number
CN202010510380.9A
Other languages
Chinese (zh)
Other versions
CN111561991A (en
Inventor
吕曜辉
汤宁
尚皓楠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Yuyantang Biological Technology Co ltd
Original Assignee
Ocean University of China
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 Ocean University of China filed Critical Ocean University of China
Priority to CN202010510380.9A priority Critical patent/CN111561991B/en
Publication of CN111561991A publication Critical patent/CN111561991A/en
Application granted granted Critical
Publication of CN111561991B publication Critical patent/CN111561991B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Holo Graphy (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于边缘填充和傅里叶变换的近场声全息方法,包括:在声源近场辐射区域内布置全息面H,并测量各测量网格点处的声压;确定声源平面尺寸,在所述全息面H处设置虚拟全息面V,所述虚拟全息面V的面积大于或等于声源平面,且中心点与所述全息面H相重合;对虚拟全息面V填充虚拟测量点数,计算每个虚拟测量点的声压值;确定重建面S坐标(x,y,zS)及波数分量kx、ky、kz,确定传递函数;根据虚拟全息面V的声压数据和传递函数进行声场重建,将空间域卷积转换为波数域乘积形式。本发明在不增加测量工作量与测量成本的前提下,降低卷绕误差和边缘吉布斯效应的影响,显著改善了小全息孔径下的重建结果。

Figure 202010510380

The invention discloses a near-field acoustic holography method based on edge filling and Fourier transform. Source plane size, a virtual holographic surface V is set at the holographic surface H, the area of the virtual holographic surface V is greater than or equal to the sound source plane, and the center point coincides with the holographic surface H; fill the virtual holographic surface V The number of virtual measurement points, calculate the sound pressure value of each virtual measurement point; determine the reconstruction surface S coordinates (x, y, z S ) and wavenumber components k x , ky , k z , determine the transfer function; According to the virtual holographic surface V The sound field is reconstructed from the sound pressure data and transfer function, and the spatial domain convolution is converted into the wavenumber domain product form. On the premise of not increasing the measurement workload and measurement cost, the invention reduces the influence of the winding error and the edge Gibbs effect, and significantly improves the reconstruction result under the small holographic aperture.

Figure 202010510380

Description

Near-field acoustic holography method based on edge filling and Fourier transform
Technical Field
The invention relates to an acoustic processing technology, in particular to a near-field acoustic holography method based on Fourier transform for a small-aperture holographic surface.
Background
Near-field Acoustic Holography (NAH) is an important sound source identification, positioning and sound field visualization technology, and not only can identify and position a noise source, but also can predict the radiation characteristics of the sound source in a sound field.
Among the near-field acoustic holography technologies, the near-field acoustic holography technology based on Fourier transform has the advantages of high calculation speed, simple theory, easy realization and the like, so the near-field acoustic holography technology is most widely applied to practical engineering. However, due to the influence of a winding error and an edge gibbs effect, the reconstruction effect of the near-field acoustic holography technology based on fourier transform is poor under a small holographic aperture, and in order to obtain higher reconstruction accuracy, the size of the holographic aperture is often required to be at least twice of the area of a sound source, which causes many problems in practical engineering application. On one hand, the measuring surface is large, the number of measuring points is very large, and the measuring workload and the measuring cost are greatly increased; on the other hand, the technology is difficult to be applied to large-size structures such as automobiles, airplanes, submarine shells and the like.
Therefore, the research on how to improve the reconstruction precision of the Fourier transform-based near-field acoustic holography under the small-aperture holographic surface has great significance for practical engineering application.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, provides a near-field acoustic holography method based on edge filling and Fourier transform, and improves the reconstruction accuracy under a small holographic aperture.
In order to achieve the purpose, the invention adopts the following technical scheme:
a near-field acoustic holography method based on edge filling and fourier transform, comprising the steps of:
(1) arranging a holographic surface H in a near-field radiation area of a sound source, wherein the holographic surface H is a square matrix in a sound field generated by the sound source, measuring grid points are distributed on the holographic surface H, and the sound pressure at the measuring grid points is measured;
(2) determining the size of a sound source plane, and setting a virtual holographic surface V at the holographic surface H, wherein the area of the virtual holographic surface V is larger than or equal to that of the sound source plane, and the central point of the virtual holographic surface V is coincident with the holographic surface H;
(3) filling virtual measurement points into the virtual holographic surface V, and calculating the sound pressure value of each virtual measurement point;
(4) determining reconstruction surface S coordinates (x, y, z)S) And the wave number component kx、ky、kzDetermining a transfer function according to the reconstruction surface coordinates and the wave number components;
(5) and (4) reconstructing a sound field according to the sound pressure data of the virtual holographic surface V and the transfer function, and converting the spatial domain convolution into a wave number domain product form.
Further, in the step (3), the number of virtual measurement points is filled in the virtual holographic surface V in a manner that the scanning lines scan line by line, and the specific method is as follows: setting a scan parallel to the x-axisA line for scanning the virtual holographic surface by moving a distance d in the positive direction of the y-axis, and when there are four intersections between the scanning line and the boundaries of the virtual holographic surface V and the holographic surface H, x is respectively set from left to right0,x1,x2,x3At x0And x1Supplementing a virtual measuring point at x every distance d2And x3Supplementing a virtual measuring point every distance d; when the scanning line only has an intersection point with the boundary of the virtual holographic surface V, supplementing a virtual measuring point every distance d between the two intersection points; when the boundaries of the scanning line, the holographic surface H and the virtual holographic surface V have no intersection point, no treatment is carried out; where d is the spacing of the holographic surface measurement grid points.
Further, the step (3) obtains the sound pressure value of each virtual measurement point by using the sound pressure value on the holographic surface H through a berger recurrence method.
Further, the step (3) further comprises: and filtering the sound pressure data of the virtual holographic surface V by using Turkey filtering.
Further, the transfer function in the step (4) is a two-dimensional space Fourier transform G of a Green functionD(kx,kyZ) of the form:
Figure BDA0002528028400000021
i is the unit of imaginary number and z is the coordinate of the spatial point on the z-axis.
Further, the reconstruction formula in the step (5) is specifically:
Figure BDA0002528028400000022
in the formula, p (x, y, z)V) Sound pressure value, F, of virtual hologram plane Vx、FyFourier transform of the x, y axes, respectively, zV、zsRespectively, the z-axis coordinates of the virtual holographic surface and the reconstruction surface.
Compared with the prior art, the invention has the following beneficial effects:
the invention provides a near-field acoustic holography method based on edge filling and Fourier transform, which reduces the influence of a winding error and an edge Gibbs effect on the premise of not increasing the measurement workload and the measurement cost, and obviously improves the reconstruction result under a small holographic aperture. Compared with the traditional Fourier transform near-field acoustic holography algorithm, the method provided by the invention greatly reduces the requirements on the holographic aperture, and has higher reconstruction precision and application range.
Drawings
FIG. 1 is a flow chart of a near-field acoustic holography method based on edge filling and Fourier transform according to the present invention;
FIG. 2 is a schematic diagram of the process before and after the virtual measurement point is supplemented;
FIG. 3 is a schematic view of the completion of the supplement of the virtual hologram surface according to the present invention;
FIG. 4 is a time domain and frequency domain diagram of the Turkey window of the present invention;
FIG. 5 is a sound pressure reconstruction contrast diagram of a conventional Fourier transform-based near-field acoustic holography algorithm under a small aperture condition;
FIG. 6 is a sound pressure reconstruction contrast diagram under a small holographic aperture condition according to the method of the present invention.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
The invention provides a near-field acoustic holography method based on edge filling and Fourier transform, which has good reconstruction precision for a small-aperture holographic surface. Referring to fig. 1, the method generally comprises an edge filling algorithm and a near-field acoustic holography algorithm based on fourier transform, wherein the edge filling algorithm fills virtual measurement points in an original holographic measurement surface, the algorithm adopts a berger recurrence method to obtain sound pressure values of the virtual measurement points, and edges are smoothed through Turkey filtering, so that a new virtual holographic surface is obtained; and (3) performing convolution conversion on the virtual holographic surface sound pressure data and the space domain of the Green function into a wave number domain product form based on a Fourier transform near-field acoustic holography algorithm to reconstruct the sound field. By utilizing the method provided by the invention, the holographic aperture does not need to be increased to reduce errors, the influence of the winding error and the edge Gibbs effect is reduced on the premise of not increasing the measurement workload and the measurement cost, and the reconstruction precision of the Fourier transform-based near-field acoustic holographic algorithm under the small-aperture holographic surface is obviously improved.
A detailed description of the steps of the near-field acoustic holography method based on edge filling and fourier transform is given below in conjunction with a simulation example.
In step 001, the present embodiment performs numerical simulation by taking 5 point sound sources as an example. The simulated point sound source has a radius of 0.01m, a frequency of 1500Hz, a vibration speed of 2.5m/s, and coordinates (-0.5, -0.5,0) (-0.5,0.5,0) (0,0,0) (0.5,0.5,0) (0.5, -0.5,0) in a spatial rectangular coordinate system. The distances between the hologram surface H and the sound source surface are 0.05m, respectively. The numerical aperture of the holographic surface, that is, the holographic aperture, is 0.5m × 0.5m, and the measurement point pitch is 0.05 m.
And step 002, for each individual point sound source, obtaining a sound pressure simulation theoretical value at the holographic surface H and the reconstruction surface S through the following point sound source radiation formula. And then overlapping the sound pressure to obtain a final theoretical value.
Figure BDA0002528028400000031
In the formula (I), the compound is shown in the specification,
Figure BDA0002528028400000032
v0、r0respectively the vibration speed and radius of the point sound source, omega and k respectively the frequency and wave number of the point sound source, i is an imaginary unit, rho0R is the distance from the point source to each point on the holographic or reconstruction surface, which is the density of the air medium. It should be understood that the above steps 001 and 002 are the calculation of the sound pressure in the holographic plane with the simulation example. In the practical application process, a holographic surface H can be arranged in a near-field radiation area of a sound source, the holographic surface H is a square matrix in a sound field generated by the sound source, measuring grid points are distributed on the holographic surface H, and a sensor is used for measuring sound pressure at the measuring grid points.
And step 003, placing a virtual holographic surface V at the position of the holographic surface H, wherein the numerical aperture of the virtual holographic surface is 1.5m multiplied by 1.5 m. In order to obtain a reconstruction value with higher precision, the numerical aperture and the number of measurement points of the reconstruction surface are the same as those of the virtual holographic surface, and the distance between the virtual holographic surface and the sound source surface is 0.03 m. In the example, the area of the sound source surface is 1m multiplied by 1m, the area of the holographic surface is 0.5m multiplied by 0.5m, the area of the sound source surface is obviously larger than that of the holographic surface, and the area of the virtual holographic surface is ensured to be larger than or equal to that of the sound source surface to be detected, so that the effectiveness of the method is verified.
Step 004, supplementing virtual measuring points on the virtual holographic surface V.
The method comprises the following specific steps:
(1) fig. 2 is a schematic diagram of a process of supplementing virtual measurement points, which is a schematic diagram before supplementation and a schematic diagram after one supplementation from left to right, wherein black dots represent original measurement points on the holographic surface H, hollow circles represent supplemented virtual measurement points, and black squares represent intersections of scanning lines with boundaries of the virtual holographic surface V and the holographic surface H. The supplementary method is as follows: a scanning line is arranged parallel to the x-axis, and the virtual holographic surface is scanned by moving a distance d (d is the distance between the measuring points of the holographic surface, and is 0.05m in the embodiment) in the positive direction of the y-axis, which is shown by a dotted line. When the scanning line has four intersections with the boundaries of the virtual holographic surface V and the holographic surface H, x is respectively set from left to right0,x1,x2,x3At x0And x1Supplementing a virtual measurement point for each distance d of the interval. For x2And x3The same process is carried out between the two, and one virtual measuring point is supplemented at every distance d. When the scanning line has an intersection only with the boundary of the virtual hologram plane V, go to (2). When the scanning line has no intersection with the boundaries of the holographic surface H and the virtual holographic surface V, no processing is performed.
(2) When the scanning line has two intersections with the virtual hologram surface V, as shown by x in FIG. 24,x5As shown, two intersection points supplement one virtual measurement point per distance d apart. And finally, replacing the point where the scanning line intersects with the virtual holographic surface with a virtual measuring point. After all the virtual measurement points are set, as shown in fig. 3, the open circles represent the supplementary virtual measurement points, and the black dots represent the original measurement points on the hologram surface H.
And 005, determining the coordinates of the virtual measurement points, and calculating the sound pressure value of each virtual measurement point.
In the method, the sound pressure value of each point on the virtual holographic surface V is reconstructed by utilizing the sound pressure value on the holographic surface H through the principle of the Berger recursion method. To obtain the point measurement values of the entire virtual hologram surface, the virtual measurement values of all rows on the virtual hologram surface are first derived in the positive and negative directions of the x-axis, and then the virtual measurement values of all columns on the virtual hologram surface are derived in the positive and negative directions of the y-axis, the virtual measurement values of the columns including the boundary filling portions that did not exist before. And setting L as the measuring point number of the original holographic surface in the x direction and S as the measuring point number of the virtual holographic surface in the x direction. Assuming that the sound pressure value of a certain line of the original holographic surface is PHrow=[p(rHrow1),p(rHrow2),···,p(rHrowL)]T. The holographic surface is expanded rightwards at first, and the right side sound pressure value is
Figure BDA0002528028400000051
The sound pressure value on the left side is
Figure BDA0002528028400000052
The sound pressure value of a certain row on the right side of the holographic surface is as follows:
Figure BDA0002528028400000053
in the formula, P is the number of the supplementary measuring points in a certain line on the right side of the original holographic surface and is less than or equal to L.
The sound pressure value of the same row on the left side of the holographic surface is as follows:
Figure BDA0002528028400000054
wherein P is the number of measurement points supplemented to a certain line on the left side of the original holographic surface and is less than or equal to L, aiThe method is a P-order autoregressive parameter obtained by fitting of a Berger algorithm, and the solving process is as follows:
is provided with
Figure BDA0002528028400000055
And
Figure BDA0002528028400000056
p-order forward and backward prediction errors, respectively, where,
Figure BDA0002528028400000057
the reflection coefficient formula is:
Figure BDA0002528028400000058
order to
Figure BDA0002528028400000059
Figure BDA00025280284000000510
ap(p)=kp
ρp=(1-|kp|2p-1
The coefficient a can be obtained1(p)=k1,ρ1=(1|k1|2)rp(0)。
By
Figure BDA00025280284000000511
Can obtain the product
Figure BDA00025280284000000512
And
Figure BDA00025280284000000513
and bring it into the reflection coefficientK can be obtained from the formula2And by analogy, all the order autoregressive parameters can be obtained until P is equal to P.
The sound pressure value calculated as described above is used as a virtual point on the virtual hologram surface. The virtual measuring point and the original measuring point on the holographic surface jointly form sound pressure data of the virtual holographic surface.
Step 006, filtering the sound pressure data of the virtual holographic surface V by using the Tukey window, thereby improving the problem that the signal-to-noise ratio of the sound pressure signal at the edge of the holographic surface is low.
As shown in fig. 4, the Tukey window is a superposition of a rectangular window and two cosine windows, and is characterized in that the main side lobe ratio of the sampled signal is high, and the side lobes converge quickly. This is used here because the constant part of the window is placed right above the original aperture, while the cosine part covers the area that fills the border. In this way, the spatial window does not affect the acoustic information and leakage is reduced.
Step 007 of determining spatial wavenumber component kx,kyAnd kzThe method comprises the following specific steps:
(1) k in the simulation example of the present invention according to the Nyquist sampling theoremx,kyThe value ranges are defined as follows:
Figure BDA0002528028400000061
Figure BDA0002528028400000062
Δ x and Δ y are sampling intervals in the x and y directions, respectively, and are determined by dividing the hologram surface size by the number of measurement points.
(2) Wave number component k in z directionzAccording to kx,kyAnd the sound source frequency k.
Figure BDA0002528028400000063
Step 008, respectively reconstructing a sound field through a Fourier transform near-field acoustic holography algorithm according to the sound pressure data of the holographic surface H and the sound pressure data of the virtual holographic surface V, and comparing the reconstruction effects of the sound field and the sound field.
Sound field reconstruction is performed according to the sound pressure data of the holographic surface H by a conventional fourier transform-based near-field acoustic holography algorithm, and the reconstruction result is shown in fig. 5, wherein the theoretical sound pressure of the reconstruction surface is calculated in step 002. The reconstruction formula based on the Fourier transform near-field acoustic holography algorithm is specifically as follows:
Figure BDA0002528028400000064
in the formula, p (x, y, z)H) Sound pressure data, transfer function G, obtained for holographic surface HD(kx,kyZ) is a two-dimensional spatial Fourier transform of the Green function,
Figure BDA0002528028400000065
Fx、Fyfourier transform of the x, y axes respectively. z is a radical ofV、zsRespectively, the z-axis coordinates of the virtual holographic surface and the reconstruction surface.
The reconstruction is carried out by using the method of the invention, the reconstruction result is shown in figure 6, the concrete steps are consistent with the reconstruction according to the holographic surface H, and only the input holographic surface sound pressure value p (x, y, z)H) Replacing the sound pressure value p (x, y, z) of the virtual holographic surface V after Tukey window filteringV) And the corresponding parameters such as the size of the reconstruction surface are changed. The concrete reconstruction formula is as follows:
Figure BDA0002528028400000071
as can be seen from fig. 5 and 6, the sound pressure reconstructed by the method of the present invention is substantially consistent with the theoretical sound pressure of the reconstruction surface, and compared with the conventional near-field acoustic holography technique of fourier transform, the reconstruction accuracy is greatly improved.
The components not specified in this embodiment can be implemented by the prior art.
The invention provides a near-field acoustic holography method based on edge filling and Fourier transform, which can be used in the fields of noise level evaluation, mechanical noise control, equipment fault diagnosis, tone quality design, loudspeaker systems, indoor acoustics and the like. For example, the position of the noise source is determined by reconstructing the sound field distribution of the whole three-dimensional space, which has great significance for effectively carrying out noise source control and noise source sound radiation characteristic research. The method greatly reduces the requirement of the traditional Fourier transform-based near-field acoustic holography algorithm on the holographic aperture, and has higher reconstruction precision and application range.
The above description is only a preferred embodiment of the present invention and there are many ways and ways to implement this solution. It should be noted that, for those skilled in the art, without departing from the principle of the present invention, several improvements and modifications can be made, and these improvements and modifications should also be construed as the protection scope of the present invention.

Claims (4)

1.一种基于边缘填充和傅里叶变换的近场声全息方法,其特征在于,包括以下步骤:1. a near-field acoustic holography method based on edge filling and Fourier transform, is characterized in that, comprises the following steps: (1)在声源近场辐射区域内布置全息面H,在由声源产生的声场中,所述全息面H为方阵,在所述全息面H上分布有各测量网格点,并测量各测量网格点处的声压;(1) Arrange the holographic surface H in the near-field radiation area of the sound source. In the sound field generated by the sound source, the holographic surface H is a square matrix, and each measurement grid point is distributed on the holographic surface H. Measure the sound pressure at each measurement grid point; (2)确定声源平面尺寸,在所述全息面H处设置虚拟全息面V,所述虚拟全息面V的面积大于或等于声源平面,且中心点与所述全息面H相重合;(2) Determine the size of the sound source plane, set a virtual holographic surface V at the holographic surface H, the area of the virtual holographic surface V is greater than or equal to the sound source plane, and the center point coincides with the holographic surface H; (3)通过扫描线逐行扫描的方式对虚拟全息面V填充虚拟测量点数,计算每个虚拟测量点的声压值,其中所述扫描线逐行扫描的方式对虚拟全息面V填充虚拟测量点数包括:平行于x轴设置一条扫描线,沿y轴的正方向每次移动距离d对虚拟全息面进行扫描,当扫描线与虚拟全息面V、全息面H边界有四个交点时,设从左至右分别为x0,x1,x2,x3,在x0与x1之间每间隔d距离补充一个虚拟测量点,在x2和x3之间每间隔d距离补充一个虚拟测量点;当扫描线只与虚拟全息面V边界有交点时,在两个交点之间每间隔d距离补充一个虚拟测量点;当扫描线与全息面H、虚拟全息面V的边界均没有交点时,不做任何处理;其中d为全息面测量网格点的间距;所述每个虚拟测量点的声压值利用全息面H上的声压值通过伯格递推法获得;(3) Filling the virtual holographic surface V with the number of virtual measurement points by scanning line by line, and calculating the sound pressure value of each virtual measuring point, wherein the scanning line scanning line by line fills the virtual holographic surface V with virtual measurement points The number of points includes: setting a scanning line parallel to the x-axis, and scanning the virtual holographic surface along the positive direction of the y-axis by moving distance d each time, when the scanning line and the virtual holographic surface V and the holographic surface H have four intersection points, set From left to right are x 0 , x 1 , x 2 , x 3 , one virtual measurement point is added every d distance between x 0 and x 1 , and one virtual measuring point is added every d distance between x 2 and x 3 Virtual measurement point; when the scan line only intersects with the boundary of the virtual holographic surface V, a virtual measurement point is added every d distance between the two intersection points; when the scan line has no boundary between the holographic surface H and the virtual holographic surface V When the intersection point, do not do any processing; wherein d is the distance between the holographic surface measurement grid points; the sound pressure value of each virtual measurement point is obtained by using the sound pressure value on the holographic surface H by the Berger recurrence method; (4)确定重建面S坐标(x,y,zS)及波数分量kx、ky、kz,根据所述重建面坐标及波数分量确定传递函数;(4) Determine the reconstructed surface S coordinates (x, y, z S ) and the wavenumber components k x , ky , and k z , and determine the transfer function according to the reconstructed surface coordinates and the wavenumber components; (5)根据虚拟全息面V的声压数据和传递函数进行声场重建,将空间域卷积转换为波数域乘积形式。(5) According to the sound pressure data and transfer function of the virtual holographic surface V, the sound field is reconstructed, and the spatial domain convolution is converted into the wavenumber domain product form. 2.根据权利要求1所述的基于边缘填充和傅里叶变换的近场声全息方法,其特征在于,所述步骤(3)中还包括:利用Turkey滤波对所述虚拟全息面V的声压数据进行滤波。2. the near-field acoustic holography method based on edge filling and Fourier transform according to claim 1, is characterized in that, in described step (3), also comprises: utilize Turkey filtering to the sound of described virtual holographic surface V pressure data for filtering. 3.根据权利要求1所述的基于边缘填充和傅里叶变换的近场声全息方法,其特征在于,所述步骤(4)中传递函数为格林函数的二维空间傅里叶变换GD(kx,ky,z),形式如下:3. the near-field acoustic holography method based on edge filling and Fourier transform according to claim 1, is characterized in that, in described step (4), transfer function is the two-dimensional space Fourier transform G D of Green's function (k x , ky , z), in the form:
Figure FDA0002899063350000011
Figure FDA0002899063350000011
i为虚数单位,z为空间点在z轴的坐标。i is the imaginary unit, and z is the coordinate of the space point on the z-axis.
4.根据权利要求3所述的基于边缘填充和傅里叶变换的近场声全息方法,其特征在于,所述步骤(5)中重建公式具体为:4. the near-field acoustic holography method based on edge filling and Fourier transform according to claim 3, is characterized in that, in described step (5), the reconstruction formula is specifically:
Figure FDA0002899063350000021
Figure FDA0002899063350000021
式中,p(x,y,zV)为虚拟全息面V的声压值,Fx、Fy分别为x,y轴的傅里叶变换,zV、zs分别是虚拟全息面及重建面所在z轴坐标。In the formula, p(x, y, z V ) is the sound pressure value of the virtual holographic surface V, F x and F y are the Fourier transform of the x and y axes, respectively, z V and z s are the virtual holographic surface and The z-coordinate of the reconstructed surface.
CN202010510380.9A 2020-06-08 2020-06-08 A near-field acoustic holography method based on edge filling and Fourier transform Active CN111561991B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010510380.9A CN111561991B (en) 2020-06-08 2020-06-08 A near-field acoustic holography method based on edge filling and Fourier transform

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010510380.9A CN111561991B (en) 2020-06-08 2020-06-08 A near-field acoustic holography method based on edge filling and Fourier transform

Publications (2)

Publication Number Publication Date
CN111561991A CN111561991A (en) 2020-08-21
CN111561991B true CN111561991B (en) 2021-03-30

Family

ID=72072584

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010510380.9A Active CN111561991B (en) 2020-06-08 2020-06-08 A near-field acoustic holography method based on edge filling and Fourier transform

Country Status (1)

Country Link
CN (1) CN111561991B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114894297B (en) * 2022-05-07 2024-11-12 湖南科技大学 A sound field separation method based on virtual double holographic surfaces

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1385265A (en) * 1971-05-28 1975-02-26 Thomson Csf Acoustical holography system for acoustic image conversion
US3964052A (en) * 1975-03-28 1976-06-15 Holosonics, Inc. Acoustical imaging system
CN102494755A (en) * 2011-12-06 2012-06-13 合肥工业大学 Unsteady sound field separation method
CN103712684A (en) * 2013-12-25 2014-04-09 广西科技大学 Sound field rebuilding method
CN103822698A (en) * 2013-12-25 2014-05-28 广西科技大学 Sound field reconstruction method based on multi-ball domain wave superposition method
CN103929706A (en) * 2013-01-11 2014-07-16 克里佩尔有限公司 Arrangement and method for measuring the direct sound radiated by acoustical sources
CN110765673A (en) * 2019-12-30 2020-02-07 上海索辰信息科技有限公司 Noise sound field reconstruction method based on equivalent source method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101125391A (en) * 2007-09-01 2008-02-20 南昌大学 Composite scanning filling method for rapid prototyping machine
CN101539455B (en) * 2009-04-22 2010-09-01 合肥工业大学 Method for re-establishing moving sound source by adopting moving equivalent source method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1385265A (en) * 1971-05-28 1975-02-26 Thomson Csf Acoustical holography system for acoustic image conversion
US3964052A (en) * 1975-03-28 1976-06-15 Holosonics, Inc. Acoustical imaging system
CN102494755A (en) * 2011-12-06 2012-06-13 合肥工业大学 Unsteady sound field separation method
CN103929706A (en) * 2013-01-11 2014-07-16 克里佩尔有限公司 Arrangement and method for measuring the direct sound radiated by acoustical sources
CN103712684A (en) * 2013-12-25 2014-04-09 广西科技大学 Sound field rebuilding method
CN103822698A (en) * 2013-12-25 2014-05-28 广西科技大学 Sound field reconstruction method based on multi-ball domain wave superposition method
CN110765673A (en) * 2019-12-30 2020-02-07 上海索辰信息科技有限公司 Noise sound field reconstruction method based on equivalent source method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
陈允峰等.《基于新型传感器阵列的声全息测试分析方法研究》.《传感器与微系统》.2017,第36卷(第8期), *
陶文俊.《近场声全息重建精度提高方法研究》.《中国优秀硕士学位论文全文数据库 工程科技II辑》.2019,(第2期), *

Also Published As

Publication number Publication date
CN111561991A (en) 2020-08-21

Similar Documents

Publication Publication Date Title
JP5496192B2 (en) Method for reconstructing an acoustic field
CN102680071B (en) Noise source identification method adopting vibration speed measurement and partial near-field acoustical holography method
CN112116616B (en) Phase information extraction method, storage medium and device based on convolutional neural network
CN105181121A (en) High-precision near-field acoustic holography algorithm adopting weighted iteration equivalent source method
CN107478325B (en) Free field reduction method of unsteady planar sound source
CN102494755B (en) Unsteady sound field separation method
CN106569191A (en) Method of acquiring target RCS by using high resolution imaging
CN113063490A (en) A sound field separation method based on double-sided measurement of sound pressure and particle vibration velocity
CN104678359B (en) A kind of porous sound holographic method of sound field identification
CN112926231B (en) A near-field acoustic holography measurement method in a limited space based on the equivalent source method
CN113758559A (en) Method for measuring shallow sea channel based on scattering sound field separation algorithm
CN112577592B (en) Finite space plane near-field acoustic holography measurement method based on space Fourier transform
CN110379011A (en) A kind of underwater topography dot cloud hole method for repairing and mending based on improvement B-spline Curve
CN111561991B (en) A near-field acoustic holography method based on edge filling and Fourier transform
CN107566969A (en) A kind of enclosed environment internal low-frequency Reconstruction of Sound Field method
CN114485917A (en) A Sound Field Reconstruction Method Based on Planar Array Scanning
CN110208764A (en) Dynamic sea large scene echo simulation method based on electromagnetism Scattering Calculation
CN112415497B (en) Method and system for separating scattering sound field in shallow sea waveguide by finite element combined wave superposition
CN116540220A (en) Target three-dimensional passive positioning method and device
CN110045333A (en) A kind of sound source three-dimensional positioning method based on Kalman filtering
Yang et al. Development and calibration of acoustic video camera system for moving vehicles
CN101515033A (en) Multilayer stereoscopic grid array for recognizing noise source in a beam shaping method
CN117233266A (en) Full waveform inversion guided wave tomography method based on cyclic neural network
CN115782480B (en) Rolling tire structure noise time domain prediction method
CN113126029B (en) A multi-sensor pulsed sound source localization method suitable for deep-sea reliable sound path environment

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
TR01 Transfer of patent right

Effective date of registration: 20240812

Address after: Room 401, Building 18, No. 333 Yingshanhong Road, Huangdao District, Qingdao City, Shandong Province, 266000

Patentee after: Qingdao Yuyantang Biological Technology Co.,Ltd.

Country or region after: China

Address before: 266100 Shandong Province, Qingdao city Laoshan District Songling Road No. 238

Patentee before: OCEAN University OF CHINA

Country or region before: China

TR01 Transfer of patent right