Method for obtaining band attenuation scattering sound waves through three-dimensional S-domain transformation

Abstract

The invention discloses a method for obtaining band attenuation scattering sound waves through three-dimensional S-domain transformation. The method comprises the steps that a binary matrix corresponding to a target scatterer is obtained; performing first transformation processing on the binary matrix to obtain a target form of the binary matrix; the first transformation processing is three-dimensional Laplace transformation; obtaining a target scattering sound pressure function based on the target form of the binary matrix, the first propagation distance and the scattering sound pressure function; and determining a target sound pressure field according to the target scattering sound pressure function and the updated first propagation distance. According to the method disclosed by the invention, by introducing the sound absorption attenuation factor, the complexity of processing the sound wave signal is improved, and the robustness to noise and interference is enhanced, so that the propagation characteristic and the attenuation characteristic of the scattering sound wave are described more accurately.

Description

Technical Field

[0001] The present invention belongs to the technical field of acoustic wave detection, and more specifically, relates to a method for obtaining attenuated scattered acoustic waves by using three-dimensional S-domain transformation. Background Art

[0002] The propagation of sound waves in a medium will generally attenuate to a certain extent with the propagation distance. Common forms include sound absorption and sound scattering. The phenomenon that sound waves encounter obstacles during propagation, causing the sound waves to deviate from the original propagation direction and diverge in different directions is called sound scattering. Sound scattering is used in geological exploration, medical imaging or industrial fields, and is one of the important means of non-destructive testing. In the prior art, sound scattering is mostly solved by analytical methods, numerical methods, etc., but it is not considered that sound absorption is also one of the causes of attenuation during acoustic propagation. Therefore, the solution to sound scattering in the prior art is still inaccurate. Summary of the invention

[0003] In order to solve the deficiencies in the prior art, the purpose of the present invention is to solve the above-mentioned defects and further propose a method for obtaining attenuated scattered sound waves by using three-dimensional S-domain transformation.

[0004] The present invention adopts the following technical solution:

[0005] The first aspect of the present invention discloses a method for obtaining attenuated scattered sound waves by using three-dimensional S-domain transformation, the method comprising:

[0006] Get the binary matrix corresponding to the target scatterer;

[0007] Performing a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transform;

[0008] Obtaining a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function;

[0009] The target sound pressure field is determined according to the target scattered sound pressure function and the updated first propagation distance.

[0010] In some possible embodiments, obtaining a binary matrix corresponding to the target scatterer includes:

[0011] Acquire propagation parameters corresponding to the target scatterer, wherein the propagation parameters include the wave number and sound velocity of the target scatterer in the propagation medium, and the wave number of the target scatterer in the obstacle;

[0012] Based on a preset binary matrix calculation formula, the binary matrix corresponding to the target scatterer is determined.

[0013] In some possible embodiments, before obtaining the binary matrix corresponding to the target scatterer, the method further includes:

[0014] Establishing a standard three-dimensional coordinate system; the standard three-dimensional coordinate system is determined based on the position of the obstacle;

[0015] Based on the standard three-dimensional coordinate system, the unit point source position and the target measurement point are determined.

[0016] In some possible embodiments, after determining the unit point source position and the target measurement point, the method further includes:

[0017] Determine a first propagation function and a second propagation function; the first propagation function corresponds to the propagation of sound waves in a propagation medium, and the second propagation function corresponds to the propagation of sound waves in an obstacle;

[0018] Determine the target sound field function corresponding to the unit point source position at the target measurement point;

[0019] Based on the Green function and the Dirac function, the target sound field function, the first propagation function and the second propagation function are subjected to function deformation processing to obtain a scattered sound pressure function.

[0020] In some possible embodiments, obtaining the target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function includes:

[0021] Substituting the target form of the binary matrix and the first propagation distance into the scattered sound pressure function;

[0022] The diffuse sound pressure expression after substitution is subjected to a second transformation process to obtain a target diffuse sound pressure function; the second transformation process is a two-dimensional inverse Laplace transform.

[0023] In some possible embodiments, determining the target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance includes:

[0024] Iterating the first propagation distance to obtain a plurality of updated first propagation distances;

[0025] Substituting the multiple updated first propagation distances into the target scattered sound pressure function to obtain multiple sound pressure fields corresponding to the multiple updated first propagation distances;

[0026] The multiple sound pressure fields are superimposed to obtain the target sound pressure field.

[0027] The second aspect of the present invention discloses a device for obtaining attenuated scattered sound waves by using three-dimensional S-domain transformation, the device comprising:

[0028] A matrix acquisition module is used to obtain the binary matrix corresponding to the target scatterer;

[0029] A first transformation module, configured to perform a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transformation;

[0030] A function determination module, configured to obtain a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function;

[0031] The sound pressure calculation module is used to determine the target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance.

[0032] The third aspect of the present invention discloses a device, which includes a processor and a memory, wherein the memory stores at least one instruction and at least one program, and the at least one instruction and the at least one program are loaded and executed by the processor to implement the method of obtaining attenuated scattered sound waves using three-dimensional S-domain transformation as described above.

[0033] The fourth aspect of the present invention discloses a computer storage medium, in which at least one instruction and at least one program are stored. The at least one instruction and the at least one program are loaded and executed by a processor to implement the method of obtaining attenuated scattered sound waves using three-dimensional S-domain transformation as described above.

[0034] Compared with the prior art, the present invention has the following advantages:

[0035] In the present invention, a binary matrix corresponding to a target scatterer is first obtained, and then a three-dimensional Laplace transform is performed on the binary matrix to obtain a target form of the binary matrix, and the target form of the binary matrix and the sound propagation distance are substituted into the scattered sound pressure function obtained in advance, and then a two-dimensional Laplace inverse transform is performed on the scattered sound pressure function after the substitution of the above two to obtain a target scattered sound pressure function; wherein, the determination of the scattered sound pressure function introduces a sound absorption attenuation factor. Compared with the prior art, the present invention calculates the sound pressure field based on the target scattered sound pressure function with the introduction of the sound absorption attenuation factor, which can more accurately describe the propagation and attenuation characteristics of scattered sound waves in a three-dimensional medium, and improves the complexity of the processable sound wave signal, and enhances the robustness to noise and interference; at the same time, in its practical application, for example, it helps to improve the resolution and contrast of the image in ultrasonic imaging, or detect defects and damage inside the material in the industrial field, thereby improving the accuracy of detection and imaging. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 It is a flowchart corresponding to the method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation provided by an embodiment of the present invention;

[0037] Figure 2 It is a schematic diagram of a corresponding process of determining a binary matrix provided by an embodiment of the present invention;

[0038] Figure 3 is a schematic diagram of a process flow corresponding to position determination provided by an embodiment of the present invention;

[0039] Figure 4 is a schematic diagram corresponding to a standard three-dimensional coordinate system provided by a specific embodiment of the present invention;

[0040] Figure 5 is a schematic diagram of a corresponding process for determining a scattered sound pressure function provided by an embodiment of the present invention;

[0041] Figure 6 is a schematic diagram of a corresponding flow chart of determining a target scattered sound pressure function provided by an embodiment of the present invention;

[0042] Figure 7 is a schematic diagram of a corresponding process for determining a target sound pressure field provided by an embodiment of the present invention;

[0043] Figure 8 It is a schematic diagram of the structure corresponding to the device for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0044] The following will be combined with the drawings in the embodiments of this specification to clearly and completely describe the technical solutions in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0045] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product, or server that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products, or devices.

[0046] Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference numerals in the accompanying drawings represent elements with the same or similar functions. Although various aspects of the embodiments are shown in the accompanying drawings, the drawings are not necessarily drawn to scale unless otherwise specified.

[0047] The word “exemplary” is used exclusively herein to mean “serving as an example, example, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments.

[0048] The term "and / or" herein is only a description of the association relationship of the associated objects, indicating that there may be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" herein represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C can represent including any one or more elements selected from the set consisting of A, B, and C.

[0049] In addition, in order to better illustrate the present disclosure, numerous specific details are given in the following specific embodiments. It should be understood by those skilled in the art that the present disclosure can also be implemented without certain specific details. In some examples, methods, means, components and circuits well known to those skilled in the art are not described in detail in order to highlight the subject matter of the present disclosure.

[0050] Figure 1 is a flow chart corresponding to the method for obtaining scattered sound waves with attenuation by using three-dimensional S-domain transformation provided by an embodiment of the present invention. The execution subject may be a system corresponding to the method for obtaining scattered sound waves with attenuation by using three-dimensional S-domain transformation, such as a system including a sound wave generator, a sound wave measuring device and other devices or structures that can implement the method for obtaining scattered sound waves with attenuation by using three-dimensional S-domain transformation; please refer to Figure 1 In one embodiment, a method for obtaining attenuated scattered sound waves using three-dimensional S-domain transform comprises the following steps:

[0051] Step S101: obtaining a binary matrix corresponding to a target scatterer;

[0052] In a specific embodiment, a scatterer refers to an obstacle or medium non-uniformity encountered by sound waves during propagation, which causes the sound waves to deviate from the original propagation direction. Scatterers can be objects of various shapes and sizes, such as underground rocks, tissue structures in medical imaging, or material defects and damage in the industrial field. The target scatterer can be a scatterer corresponding to the deviation of the target sound wave in the process of solving the scattered sound pressure; wherein the sound wave can be emitted by a sound wave generator, or by any other structure or device that can realize the propagation of sound waves; a three-dimensional scanner can be used to collect data on the target scatterer, and the data here is equivalent to a four-dimensional array. When scanning the scatterer, the three-dimensional coordinates of the scanning point can be obtained, and the constant obtained by judging whether the scanning point is in the medium or in the target is, therefore, it is equivalent to obtaining a four-dimensional array.

[0053] The propagation process of sound waves can be divided into the propagation of sound waves in obstacles and the propagation of sound waves in propagation media. Therefore, when solving the binary matrix corresponding to the target scatterer, different distinctions can be made based on the propagation medium and obstacles. Specifically, Figure 2 A schematic diagram of a flow chart of determining a binary matrix corresponding to an embodiment of the present invention; and a flow chart of determining a binary matrix corresponding to a target scatterer, such as Figure 2 As shown, step S101 may include:

[0054] Step S201: Acquire propagation parameters corresponding to the target scatterer, wherein the propagation parameters include the wave number and sound velocity of the target scatterer in the propagation medium, and the wave number of the target scatterer in the obstacle;

[0055] In a specific embodiment, a binary matrix of a target scatterer, i.e., a target sound wave, can be calculated based on multiple propagation parameters corresponding to the target scatterer, wherein the propagation parameters may include at least the wave number and sound speed of the sound wave propagating in the medium, and the wave number and sound speed of the sound wave propagating in the obstacle, etc.; wherein the wave number is the ratio of the frequency of the sound wave propagating in the propagation medium and / or obstacle to the sound speed, and the size of the wave number represents the propagation characteristics of the sound wave in the propagation medium and / or obstacle, and the larger the wave number, the higher the frequency, the shorter the wavelength of the sound wave, and the stronger the propagation ability.

[0056] Step S202: Determine the binary matrix corresponding to the target scatterer based on a preset binary matrix calculation formula.

[0057] In a specific embodiment, the preset binary matrix calculation formula may be:

[0058]

[0059] Where k0 is the wave number corresponding to the propagation of the sound wave in the obstacle, c0 is the sound speed of the sound wave in the propagation medium, and c1 is the sound speed of the sound wave in the obstacle. Based on formula (1), the binary matrix corresponding to the target scatterer can be determined.

[0060] Figure 3 is a schematic diagram of a flow chart corresponding to position determination provided by an embodiment of the present invention; Figure 3 As shown, before obtaining the binary matrix corresponding to the target scatterer, the method further includes:

[0061] Step S301: establishing a standard three-dimensional coordinate system; the standard three-dimensional coordinate system is determined based on the position of the obstacle;

[0062] In a specific embodiment, Figure 4 is a schematic diagram corresponding to a standard three-dimensional coordinate system provided by a specific embodiment of the present invention; Figure 4 As shown, the obstacle can be a solid sphere, a hollow sphere or other three-dimensional obstacles of any shape and density. The present invention does not limit the specific shape of the obstacle. Figure 4 The establishment of the standard three-dimensional coordinate system is illustrated by taking the obstacle as a solid sphere.

[0063] The standard three-dimensional coordinate system can be a Cartesian coordinate system; the obstacle can be set inside the propagation medium. When the obstacle is a solid sphere, the center of the solid sphere can be regarded as the origin of the standard three-dimensional coordinate system, and its x-axis direction is the sound wave propagation direction. The distance on the x-axis from the origin is the sound wave propagation distance.

[0064] Step S302: Determine the unit point source position and the target measurement point based on the standard three-dimensional coordinate system.

[0065] In a specific embodiment, the unit point source position may be the position where the sound pressure is generated during the propagation of the sound wave; the position corresponding to the target measurement point may be the position where the sound pressure is measured; specifically, in the established standard three-dimensional coordinate system, the unit point source position may be the position where the sound wave is generated, and the target measurement point position may be the target position where the sound pressure of the sound wave is measured.

[0066] Establish as above and as Figure 4 After the standard three-dimensional coordinate system is shown, it can be considered that the scattered sound pressure P corresponding to the propagation of the sound wave in this space s(r) conforms to the Helmholtz equation. Based on this, the position of the unit point source can be set as: r′=ix′+iy′+kz′, and the position of the target measurement point that generates the sound pressure is: r=ix+jy+kz; where i, j, k can be the basis vectors of the standard three-dimensional coordinate system, representing the directions of the x, y, and z axes respectively; x, y, and z can be used to represent the coordinates in the three-dimensional standard coordinate system.

[0067] Further, Figure 5 is a schematic diagram of a corresponding flow chart of determining a scattered sound pressure function provided by an embodiment of the present invention; Figure 5 As shown, after determining the unit point source position and the target measurement point, the method further includes:

[0068] Step S501: determining a first propagation function and a second propagation function; the first propagation function corresponds to the propagation of sound waves in a propagation medium, and the second propagation function corresponds to the propagation of sound waves in an obstacle;

[0069] In a specific embodiment, the first propagation function and the second propagation function may be further solved based on the unit point source position and the target measurement point position obtained in step S302; wherein the first propagation function may correspond to the mathematical relationship of the propagation of sound waves in the propagation medium, and the second propagation function may correspond to the mathematical relationship of the propagation of sound waves in the obstacle;

[0070] The first propagation function can be:

[0071]

[0072] And, the second propagation function can be:

[0073]

[0074] Wherein, in formula (2) and formula (3), k can represent the wave number corresponding to the propagation of the sound wave in the propagation medium, k0 can represent the wave number corresponding to the propagation of the sound wave in the obstacle, c0 is the sound speed of the sound wave in the propagation medium, c1 is the sound speed of the sound wave in the obstacle, p(r) can be the total sound pressure, p(r)=p i (r)+p s (r), where p i (r) is the incident sound pressure.

[0075] Step S502: determining a target sound field function corresponding to the unit point source position at the target measurement point;

[0076] In a specific embodiment, the target sound field function corresponding to the unit point source position at the target measurement point position is determined. The target sound field function may be:

[0077]

[0078] Among them, g(r,r′) is the Green function and δ(rr′) is the Dirac function.

[0079] Step S503: Based on the Green function and the Dirac function, function deformation processing is performed on the target sound field function, the first propagation function and the second propagation function to obtain a scattered sound pressure function.

[0080] In a specific embodiment, the scattered sound pressure is the sound pressure generated in different directions when the sound wave is scattered after encountering an obstacle during the propagation process; during the propagation process, the scattered sound wave will attenuate due to sound absorption, sound scattering, etc., so the scattered sound pressure can be used as an important parameter to describe the propagation and attenuation characteristics of the scattered sound wave.

[0081] Multiply the above formula (4) by the scattered sound pressure p s (r), implement function deformation processing, and obtain formula (5), which is:

[0082]

[0083] Then, the formula (2) corresponding to the first propagation function is multiplied by the Green function g(r,r′) to achieve function deformation processing and obtain formula (6), which is:

[0084]

[0085] And, multiplying the formula (3) corresponding to the second propagation function by the Green function g(r, r′) to implement function deformation processing, and obtain formula (7), which is:

[0086]

[0087] Furthermore, using formula (6) and formula (7), we make a difference with formula (5) respectively, and integrate them three times to eliminate the Dirac function, and we get formula (8):

[0088]

[0089] And, formula (9) is obtained as:

[0090]

[0091] When the sound wave propagation condition belongs to weak scattering, that is, when the density and sound speed of the obstacle and the medium are similar, or it can be considered that the scattered sound pressure is much smaller than the incident sound pressure, the weak scattering condition is satisfied; then formula (8) can be written as:

[0092] p s (r)=0 (10)

[0093] And formula (9) can be written as:

[0094]

[0095] Therefore, the scattered sound pressure function is obtained by combining the preset binary matrix calculation formula of formula (1). The scattered sound pressure function can be expressed as:

[0096] p s (r)=∫∫∫q(r′)g(r, r′)p i (r′)dr′ (12)

[0097] Step S102: performing a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transform;

[0098] In a specific embodiment, the first transformation processing can be a three-dimensional Laplace transform, which can improve the accuracy and efficiency of the sound wave measurement process and the sound wave analysis process; since the attenuation of sound waves includes both sound scattering and sound absorption, therefore, considering the attenuation of the amplitude of the incident wave caused by sound absorption with the propagation distance during the propagation process, a sound absorption attenuation factor can be introduced to improve the accuracy of obtaining attenuated scattered sound waves based on the target scattered sound pressure function, as well as to improve the complexity of the processable sound wave signals and enhance the robustness to noise and interference.

[0099] Specifically, the incident plane wave can be expressed as:

[0100] P i =A x e -jωt (13)

[0101] Among them, P i is the incident plane wave, A x =A0e -μx .

[0102] Furthermore, the Green's function g(r, r') can be expressed in the form of plane wave decomposition in a three-dimensional sound field, and formula (14) can be obtained, which is:

[0103]

[0104] The integral variable satisfies For the sound pressure on the forward channel, when the first propagation distance x corresponding to the target measurement point is x=l0>0, the incident wave is normalized, and the normalized incident wave and the plane wave decomposition form of the Green's function are substituted into the integral equation to obtain formula (15):

[0105]

[0106] At the same time, since s1 = μ + j (α - k0), s2 = jβ, s3 = jγ, the target form Q of the binary matrix can be solved as:

[0107]

[0108] Step S103: obtaining a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function;

[0109] After solving the target form Q(s1, s2, s3) of the binary matrix according to the above steps, the target scattered sound pressure function is obtained according to the target form of the binary matrix, the first propagation distance and the scattered sound pressure function.

[0110] In a specific embodiment, Figure 6 is a schematic diagram of a corresponding flow chart of determining a target scattered sound pressure function provided by an embodiment of the present invention; Figure 6 As shown, the target scattered sound pressure function can be obtained by the following steps:

[0111] Step S601: Substituting the target form of the binary matrix and the first propagation distance into the scattered sound pressure function;

[0112] In a specific embodiment, formula (16) and the first propagation distance are substituted into formula (15), where the first propagation distance may be the distance from the origin to the target position, and more specifically may be the propagation distance x=10. Therefore, the expression of the two-dimensional Laplace transform of the scattered sound pressure after substituting the above two is obtained, that is, the two-dimensional S-domain transform of the scattered sound pressure with attenuated scattered sound waves can be expressed as formula (17):

[0113]

[0114] Step S602: performing a second transformation process on the S-domain expression of the scattered sound pressure after substitution to obtain a target scattered sound pressure function; the second transformation process is a two-dimensional inverse Laplace transform.

[0115] In a specific embodiment, further, the second transformation processing may be a two-dimensional Laplace inverse transform; by performing Laplace inverse transform on β and γ on both sides of formula (17), the scattered sound pressure with attenuated scattered sound waves can be obtained, that is, the target scattered sound pressure function formula (18), which is:

[0116]

[0117] Step S104: determining a target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance.

[0118] In a specific embodiment, the first propagation distance can be re-determined, that is, the data corresponding to the first propagation distance is updated, and the updated first propagation distance is substituted into the target scattered sound pressure function to obtain the first sound pressure field corresponding to the updated first propagation distance; further, the data corresponding to the updated first propagation distance is updated again, and the first propagation distance after the update is substituted into the target scattered sound pressure function to obtain the second sound pressure field corresponding to the first propagation distance after the update; the iteration of the first propagation distance is repeatedly executed, and the first propagation distance after multiple iterations is substituted into the target scattered sound pressure function respectively to obtain multiple sound pressure fields corresponding to the first propagation distance after multiple iterations; when the number of iterations reaches the iteration threshold, the updating of the first propagation distance is stopped, and the multiple sound pressure fields are superimposed to obtain the target sound pressure field.

[0119] Specifically, Figure 7 is a schematic diagram of a corresponding flow chart of determining a target sound pressure field according to an embodiment of the present invention; Figure 7 As shown, step S104 may include:

[0120] Step S701: iterating the first propagation distance to obtain a plurality of updated first propagation distances;

[0121] In a specific embodiment, after obtaining the target scattered sound pressure function, the first propagation distance is updated to obtain the second propagation distance, the second propagation distance is updated to obtain the third propagation distance, the third propagation distance is updated to obtain the fourth propagation distance, and the first propagation distance is iterated n times in sequence until the number of iterations n reaches the iteration threshold, and multiple updated first propagation distances are obtained, and the multiple updated first propagation distances may include the second propagation distance, the third propagation distance... and the nth propagation distance. Among them, the iteration threshold can be preset, and can also be adjusted according to different obstacles and propagation media to ensure the accuracy of the target scattered sound pressure function solution process.

[0122] Step S702: Substituting the multiple updated first propagation distances into the target scattered sound pressure function to obtain multiple sound pressure fields corresponding to the multiple updated first propagation distances;

[0123] In a specific embodiment, each propagation distance is substituted into the target scattered sound pressure function to solve the sound pressure field, and a second sound pressure field corresponding to the second propagation distance, a third sound pressure field corresponding to the third propagation distance... and an nth sound pressure field corresponding to the nth propagation distance are obtained.

[0124] Step S703: superimposing the multiple sound pressure fields to obtain the target sound pressure field.

[0125] In a specific embodiment, multiple sound pressure fields are superimposed to obtain a target sound pressure field, thereby completing the acquisition of attenuated scattered sound waves.

[0126] In summary, the method for obtaining scattered sound waves with attenuation by using three-dimensional S-domain transformation disclosed in the present invention first obtains a binary matrix corresponding to a target scatterer, and then performs a three-dimensional Laplace transform on the binary matrix to obtain a target form of the binary matrix, and substitutes the target form of the binary matrix and the sound propagation distance into the scattered sound pressure function obtained by pre-solution, and then performs a two-dimensional Laplace inverse transform on the scattered sound pressure function after substituting the above two to obtain a target scattered sound pressure function; wherein, the determination of the scattered sound pressure function introduces an acoustic absorption attenuation factor. Compared with the prior art, the present invention calculates the sound pressure field based on the target scattered sound pressure function with the acoustic absorption attenuation factor, which can more accurately describe the propagation and attenuation characteristics of scattered sound waves in a three-dimensional medium, and improves the complexity of processable sound wave signals, and enhances the robustness to noise and interference; at the same time, in its practical application, for example, it helps to improve the resolution and contrast of images in ultrasonic imaging, or detect defects and damage inside materials in the industrial field, thereby improving the accuracy of detection and imaging.

[0127] The following is a description of the device for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation provided by the present invention. The system control device described below and the method for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation described above can be referenced to each other.

[0128] like Figure 8 As shown, a device for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation comprises:

[0129] A matrix acquisition module 810 is used to acquire a binary matrix corresponding to a target scatterer;

[0130] A first transformation module 820 is used to perform a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transform;

[0131] A function determination module 830, configured to obtain a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function;

[0132] The sound pressure calculation module 840 is used to determine the target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance.

[0133] In some other embodiments, the matrix acquisition module 810 includes:

[0134] A propagation parameter determination module, used to obtain propagation parameters corresponding to the target scatterer, wherein the propagation parameters include the wave number and sound speed of the target scatterer in the propagation medium, and the wave number of the target scatterer in the obstacle;

[0135] The matrix calculation module is used to determine the binary matrix corresponding to the target scatterer based on a preset binary matrix calculation formula.

[0136] In some other embodiments, the device further comprises:

[0137] A coordinate system establishment module, used to establish a standard three-dimensional coordinate system; the standard three-dimensional coordinate system is determined based on the position of the obstacle;

[0138] The position determination module is used to determine the unit point source position and the target measurement point based on the standard three-dimensional coordinate system.

[0139] In some other embodiments, the device further comprises:

[0140] A first determination module, used to determine a first propagation function and a second propagation function; the first propagation function corresponds to the propagation of sound waves in a propagation medium, and the second propagation function corresponds to the propagation of sound waves in an obstacle;

[0141] A second determination module is used to determine a target sound field function corresponding to the unit point source position at the target measurement point;

[0142] The first function processing module is used to perform function deformation processing on the target sound field function, the first propagation function and the second propagation function based on the Green function and the Dirac function to obtain a scattered sound pressure function.

[0143] In some other embodiments, the function determination module 830 includes:

[0144] A parameter substitution module, used for substituting the target form of the binary matrix and the first propagation distance into the scattered sound pressure function;

[0145] The second function processing module is used to perform a second transformation on the scattered sound pressure expression after substitution to obtain a target scattered sound pressure function; the second transformation is a two-dimensional inverse Laplace transform.

[0146] In some other embodiments, the sound pressure calculation module 840 includes:

[0147] a parameter updating module, configured to iterate the first propagation distance to obtain a plurality of updated first propagation distances;

[0148] A third determination module is used to substitute the multiple updated first propagation distances into the target scattered sound pressure function to obtain multiple sound pressure fields corresponding to the multiple updated first propagation distances;

[0149] The sound pressure superposition module is used to superimpose the multiple sound pressure fields to obtain the target sound pressure field.

[0150] The device and method embodiments in the device embodiment are based on the same inventive concept and are used to implement the above-mentioned method of obtaining attenuated scattered sound waves using three-dimensional S-domain transformation.

[0151] An embodiment of the present invention also provides an electronic device, comprising: a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement a method for obtaining attenuated scattered sound waves using three-dimensional S-domain transform as described in any of the method embodiments.

[0152] The embodiment of the present invention also provides a storage medium, and the computer-readable storage medium can be a tangible device that keeps and stores the instructions used by the instruction execution device. The computer-readable storage medium can be an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the above. More specific examples (non-exhaustive list) of computer-readable storage media include: portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), static random access memories (SRAM), portable compact disk read-only memories (CD-ROM), digital versatile disks (DVD), memory sticks, floppy disks, mechanical encoding devices, such as punch cards or convex structures in grooves on which instructions are stored, and any suitable combination of the above. The computer-readable storage medium used here is not interpreted as a transient signal itself, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagated by waveguides or other transmission media (e.g., light pulses through optical fiber cables), or electrical signals transmitted by wires.

[0153] The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, optical fiber transmission, wireless transmission, routers, firewalls, switches, gateway computers, and / or edge servers. The network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing / processing device.

[0154] The computer program instructions for performing the operation of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Smalltalk, C++, etc., and conventional procedural programming languages ​​such as "C" language or similar programming languages. Computer-readable program instructions may be executed entirely on a user's computer, partially on a user's computer, as an independent software package, partially on a user's computer, partially on a remote computer, or completely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., using an Internet service provider to connect via the Internet). In some embodiments, an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), may be personalized by utilizing the state information of the computer-readable program instructions, and the electronic circuit may execute the computer-readable program instructions, thereby realizing various aspects of the present disclosure.

[0155] Various aspects of the present disclosure are described herein with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer-readable program instructions.

[0156] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device that implements the functions / actions specified in one or more boxes in the flowchart and / or block diagram is generated. These computer-readable program instructions can also be stored in a computer-readable storage medium, and these instructions cause the computer, programmable data processing device, and / or other equipment to work in a specific manner, so that the computer-readable medium storing the instructions includes a manufactured product, which includes instructions for implementing various aspects of the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0157] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operating steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0158] The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system, method and computer program product according to multiple embodiments of the present disclosure. In this regard, each square frame in the flow chart or block diagram can represent a part of a module, program segment or instruction, and a part of the module, program segment or instruction includes one or more executable instructions for realizing the specified logical function. In some alternative implementations, the functions marked in the square frame can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square frames can actually be executed substantially in parallel, and they can sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each square frame in the block diagram and / or flow chart, and the combination of the square frames in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of special-purpose hardware and computer instructions.

[0159] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that the specific implementation modes of the present invention can still be modified or replaced by equivalents, and any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A method for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation, characterized in that: The method comprises: Get the binary matrix corresponding to the target scatterer; Performing a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transform; Obtaining a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function; The target sound pressure field is determined according to the target scattered sound pressure function and the updated first propagation distance.

2. The method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation according to claim 1, characterized in that: The step of obtaining a binary matrix corresponding to the target scatterer includes: Acquire propagation parameters corresponding to the target scatterer, wherein the propagation parameters include the wave number and sound velocity of the target scatterer in the propagation medium, and the wave number of the target scatterer in the obstacle; Based on a preset binary matrix calculation formula, the binary matrix corresponding to the target scatterer is determined.

3. The method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation according to claim 1, characterized in that: Before obtaining the binary matrix corresponding to the target scatterer, the method further includes: Establishing a standard three-dimensional coordinate system; the standard three-dimensional coordinate system is determined based on the position of the obstacle; Based on the standard three-dimensional coordinate system, the unit point source position and the target measurement point are determined.

4. The method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation according to claim 3, characterized in that: After determining the unit point source position and the target measurement point, the method further includes: Determine a first propagation function and a second propagation function; the first propagation function corresponds to the propagation of sound waves in a propagation medium, and the second propagation function corresponds to the propagation of sound waves in an obstacle; Determine the target sound field function corresponding to the unit point source position at the target measurement point; Based on the Green function and the Dirac function, the target sound field function, the first propagation function and the second propagation function are subjected to function deformation processing to obtain a scattered sound pressure function.

5. The method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation according to claim 1, characterized in that: The step of obtaining a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function includes: Substituting the target form of the binary matrix and the first propagation distance into the scattered sound pressure function; The diffuse sound pressure expression after substitution is subjected to a second transformation process to obtain a target diffuse sound pressure function; the second transformation process is a two-dimensional inverse Laplace transform.

6. The method for obtaining attenuated scattered sound waves by three-dimensional S-domain transformation according to claim 1, characterized in that: The determining the target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance includes: Iterating the first propagation distance to obtain a plurality of updated first propagation distances; Substituting the multiple updated first propagation distances into the target scattered sound pressure function to obtain multiple sound pressure fields corresponding to the multiple updated first propagation distances; The multiple sound pressure fields are superimposed to obtain the target sound pressure field.

7. A device for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation, characterized in that: The device comprises: A matrix acquisition module is used to obtain the binary matrix corresponding to the target scatterer; A first transformation module, configured to perform a first transformation process on the binary matrix to obtain a target form of the binary matrix; the first transformation process is a three-dimensional Laplace transformation; A function determination module, configured to obtain a target scattered sound pressure function based on the target form of the binary matrix, the first propagation distance, and the scattered sound pressure function; The sound pressure calculation module is used to determine the target sound pressure field according to the target scattered sound pressure function and the updated first propagation distance.

8. An electronic device, comprising a processor and a memory, wherein the memory stores at least one instruction and at least one program, and the at least one instruction and the at least one program are loaded and executed by the processor to implement the method for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation as described in any one of claims 1 to 6.

9. A computer storage medium, wherein at least one instruction and at least one program are stored in the computer storage medium, wherein the at least one instruction and the at least one program are loaded and executed by a processor to implement a method for obtaining attenuated scattered sound waves using three-dimensional S-domain transformation as described in any one of claims 1 to 6.

Images