Array gain calculation method and device for deep sea near sea surface vertical array buoy

Abstract

The embodiment of the invention provides an array gain calculation method and device for a deep-sea near-sea-surface vertical array buoy, and relates to the technical field of array signal processing. The method is applied to a vertical array buoy and comprises the following steps: acquiring related parameters of the vertical array buoy and noise source parameters; calculating a first contribution value of near-field noise to noise cross-spectrum density and a second contribution value of far-field noise to noise cross-spectrum density based on the vertical array buoy related parameters and the noise source parameters; determining a noise cross-spectral density matrix based on the first contribution value and the second contribution value; determining the noise directivity gain of the vertical array buoy based on the beam scanning pitch angle and the noise source parameter; determining the signal directivity gain of the vertical array buoy based on the beam scanning pitch angle, the received signal and the signal cross-spectral density matrix; based on the noise directivity gain and the signal directivity gain, the array directivity gain of the vertical array buoy is determined, and then the array gain of the vertical array buoy is determined. According to the invention, the accuracy of array gain calculation can be improved.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of array signal processing technology, and more specifically, the embodiments of the present invention relate to a method and apparatus for calculating the array gain of a deep-sea near-shore vertical array buoy. Background Technology

[0002] Array gain is a key parameter in sonar equations that determines detection range; it is defined as the logarithm of the ratio of array signal gain to noise gain. The accurate calculation of this parameter directly affects the accurate evaluation of sonar array performance. With the development of submarine stealth technology, the detection range of deep-sea near-surface vertical array buoys has drastically decreased. The background noise in the detectable area of ​​vertical array buoys mainly comes from overhead noise in the near-vertical direction. Most existing noise gain calculation models assume that the noise is completely uncorrelated, or use normal mode models to calculate the spatial correlation characteristics of the noise field to obtain the array gain. However, normal mode models have large errors in the near-field sound field and cannot accurately calculate the array gain. Summary of the Invention

[0003] In this context, embodiments of the present invention aim to provide a method and apparatus for calculating the array gain of a deep-sea near-surface vertical array buoy, which can improve the accuracy of the array gain calculation.

[0004] In a first aspect of the present invention, a method for calculating the array gain of a deep-sea near-surface vertical array buoy is provided, applied to a vertical array buoy pre-installed in the deep sea, wherein an infinite number of wind-generated noise sources are also installed in the sea area where the vertical array buoy is located, the method comprising: Obtain relevant parameters of the vertical array buoy and noise source parameters; wherein, the relevant parameters of the vertical array buoy include at least the beam scanning elevation angle, the signal cross-spectral density matrix, the received signal, and the buoy coordinates; The first contribution value of near-field noise to the noise cross-spectral density is calculated based on the relevant parameters of the vertical array buoy and the noise source parameters. The second contribution value of far-field noise to the noise cross-spectral density is calculated based on the relevant parameters of the vertical array buoy and the noise source parameters. Based on the first contribution value and the second contribution value, the noise cross-spectral density matrix is ​​determined; Based on the noise cross-spectral density matrix and the beam scanning elevation angle, the noise directivity gain of the vertical array buoy is determined; Based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix, the signal directional gain of the vertical array buoy is determined. Based on the noise directional gain and the signal directional gain, the base array directional gain of the vertical array buoy is determined, and then the array gain of the vertical array buoy is determined.

[0005] In one embodiment of this implementation, the noise source parameters include noise source intensity, horizontal wavenumber, maximum horizontal wavenumber, minimum horizontal wavenumber, and the number of horizontal wavenumber samples; the formula for calculating the first contribution value of the near-field noise to the noise cross-spectral density is: ; ; ; in, This represents the first contribution of near-field noise to the noise cross-spectral density. and This represents the depth of any two elements in the vertical array buoy; The noise source intensity is represented by k; the horizontal wavenumber is represented by k. Indicates the maximum horizontal wavenumber. The minimum horizontal wavenumber is represented by N; the number of horizontal wavenumber samples is represented by a; the radius of the region corresponding to the near-field noise is represented by r; and the calculated horizontal distance is represented by r. Represents the Green's function based on the fast field model; This indicates the conjugate operation. For discrete wavenumbers.

[0006] In one embodiment of this implementation, the formula for calculating the second contribution value of the far-field noise to the noise cross-spectral density is as follows: ; in, This represents the second contribution of far-field noise to the noise cross-spectral density, where m represents the m-th normal mode. The real part of the horizontal wavenumber of the m-th normal mode is given by [the real part of the horizontal wavenumber]. The imaginary part of the horizontal wavenumber of the m-th normal mode is given by [the term in Chinese]. Indicates the depth of the noise source. This represents the amplitude of the m-th normal mode excited by the wind-generated noise source. For depth The amplitude of the m-th normal mode. For the depth of the conjugate operation The amplitude of the m-th normal mode.

[0007] In one embodiment of this implementation, the formula for calculating the noise cross-spectral density matrix is: ; in, This represents the noise cross-spectral density matrix.

[0008] In one embodiment of this implementation, the formula for calculating the noise directivity gain of the vertical array buoy is: ; ; in, This represents the noise directionality gain of the vertical array buoy. Indicates the beam scanning elevation angle, The guide vector represents the beam scanning elevation angle, and H represents the conjugate transpose operation. The noise cross-spectral density matrix represents the noise cross-spectral density matrix. This indicates the operating frequency of the vertical array buoys. This indicates the element spacing of the vertical array buoys. Indicates the number of beam scanning angles. This represents the speed of sound in water. This indicates the number of array elements of the vertical array buoy.

[0009] In one embodiment of this implementation, the formula for calculating the signal directivity gain of the vertical array buoy is: ; ; in, p represents the signal directional gain of the vertical array buoy, and p represents the received signal. This represents the cross-spectral density matrix of the signal.

[0010] In one embodiment of this implementation, the formula for calculating the array directivity gain of the vertical array buoy is: ; in, This represents the array directivity gain of the vertical array buoy. Within the beam scanning range, the array directivity gain... The maximum value is the array gain of the vertical array buoy. .

[0011] In a second aspect of the present invention, a gain calculation device for a deep-sea near-surface vertical array buoy is provided, applied to a vertical array buoy pre-installed in the deep sea, wherein an infinite number of wind-generated noise sources are also provided in the sea area where the vertical array buoy is located, and the device includes: The acquisition unit is used to acquire relevant parameters of the vertical array buoy and noise source parameters; wherein, the relevant parameters of the vertical array buoy include at least the beam scanning elevation angle, the signal cross-spectral density matrix, and the received signal; The first calculation unit is used to calculate the first contribution value of near-field noise to the noise cross-spectral density based on the relevant parameters of the vertical array buoy and the noise source parameters. The second calculation unit is used to calculate the second contribution value of far-field noise to noise cross-spectral density based on the vertical array buoy correlation parameters and the noise source parameters. The first determining unit is used to determine the noise cross-spectral density matrix based on the first contribution value and the second contribution value; The second determining unit is used to determine the noise directivity gain of the vertical array buoy based on the noise cross-spectral density matrix and the beam scanning elevation angle. The third determining unit is used to determine the signal directivity gain of the vertical array buoy based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix. The fourth determining unit is used to determine the array directional gain of the vertical array buoy based on the noise directional gain and the signal directional gain, and then determine the array gain of the vertical array buoy.

[0012] In one embodiment of this implementation, the noise source parameters include noise source intensity, horizontal wavenumber, maximum horizontal wavenumber, minimum horizontal wavenumber, and the number of horizontal wavenumber samples; the formula for calculating the first contribution value of the near-field noise to the noise cross-spectral density is: ; ; ; in, This represents the first contribution of near-field noise to the noise cross-spectral density. and This represents the depth of any two elements in the vertical array buoy; The noise source intensity is represented by k; the horizontal wavenumber is represented by k. Indicates the maximum horizontal wavenumber. The minimum horizontal wavenumber is represented by N; the number of horizontal wavenumber samples is represented by a; the radius of the region corresponding to the near-field noise is represented by r; and the calculated horizontal distance is represented by r. Represents the Green's function based on the fast field model; This indicates the conjugate operation. For discrete wavenumbers.

[0013] In one embodiment of this implementation, the formula for calculating the second contribution value of the far-field noise to the noise cross-spectral density is as follows: ; in, This represents the second contribution of far-field noise to the noise cross-spectral density, where m represents the m-th normal mode. The real part of the horizontal wavenumber of the m-th normal mode is given by [the real part of the horizontal wavenumber]. The imaginary part of the horizontal wavenumber of the m-th normal mode is given by [the term in Chinese]. Indicates the depth of the noise source. This represents the amplitude of the m-th normal mode excited by the wind-generated noise source. For depth The amplitude of the m-th normal mode. For the depth of the conjugate operation The amplitude of the m-th normal mode.

[0014] In a third aspect of the present invention, a computing device is provided, the computing device comprising: at least one processor, a memory, and an input / output unit; wherein the memory is used to store a computer program, and the processor is used to invoke the computer program stored in the memory to execute the method described in any one aspect.

[0015] In a fourth aspect of the present invention, a computer-readable storage medium is provided, comprising instructions which, when executed on a computer, cause the computer to perform the method described in any one of the first aspects.

[0016] In a fifth aspect of the present invention, a computer program product is provided, comprising a computer program that, when executed by a processor, implements the method described in any one of the first aspects.

[0017] The method and apparatus for calculating the array gain of a deep-sea near-surface vertical array buoy according to embodiments of the present invention obtains relevant parameters of the vertical array buoy and noise source parameters, calculates the contributions of near-field noise and far-field noise to the noise cross-spectral density respectively, and then determines the noise cross-spectral density matrix. This comprehensive approach, considering the influence of near-field and far-field noise, makes the noise correlation calculation more accurate. Simultaneously, it combines multiple parameters to determine the noise directional gain and signal directional gain, comprehensively considering factors such as beam scanning elevation angle, received signal, and signal cross-spectral density matrix. Finally, the array gain is determined based on the noise directional gain and signal directional gain. The entire calculation process fully considers various practical factors, improving the accuracy of the calculated array gain. Attached Figure Description

[0018] The above and other objects, features, and advantages of exemplary embodiments of the present invention will become readily apparent from the following detailed description taken in conjunction with the accompanying drawings. Several embodiments of the invention are illustrated in the drawings by way of example and not limitation, wherein: Figure 1 This is a flowchart illustrating a method for calculating the array gain of a deep-sea near-surface vertical array buoy, as provided in an embodiment of the present invention. Figure 2 This is a schematic diagram of the sound velocity profile provided in an embodiment of the present invention; Figure 3A schematic diagram comparing the noise vertical directivity of a method for calculating the array gain of a deep-sea near-surface vertical array buoy according to an embodiment of the present invention with that based on a normal mode model. Figure 4 a is a schematic diagram of the array directional gain of a vertical array buoy, which is a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention. Figure 4 b is a schematic diagram of the array directional gain of a vertical array buoy based on a normal wave model provided in an embodiment of the present invention; Figure 4 c is a schematic diagram showing the difference between the array directional gain and the array directional gain based on the normal mode model in a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention. Figure 5 This is a schematic diagram showing the change of the angle of arrival of the direct sound ray, the sound ray reflected from the sea surface, and the sound ray reflected from the first seabed at the array element 1 of the present invention as a function of horizontal distance. Figure 6 a is a schematic diagram comparing the noise gain of a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention with the noise gain based on a normal mode model. Figure 6 b is a schematic diagram comparing the array gain of a deep-sea near-shore vertical array buoy with the array gain based on the normal mode model, according to an embodiment of the present invention. Figure 7 This is a schematic diagram of the array gain calculation device for a deep-sea near-shore vertical array buoy provided in an embodiment of the present invention; Figure 8 A schematic diagram of the structure of a medium according to an embodiment of the present invention is shown. Figure 9 A schematic diagram of the structure of a computing device according to an embodiment of the present invention is shown.

[0019] In the accompanying drawings, the same or corresponding reference numerals indicate the same or corresponding parts. Detailed Implementation

[0020] The principles and spirit of the invention will now be described with reference to several exemplary embodiments. It should be understood that these embodiments are given merely to enable those skilled in the art to better understand and implement the invention, and are not intended to limit the scope of the invention in any way. Rather, these embodiments are provided to make this disclosure more thorough and complete, and to fully convey the scope of this disclosure to those skilled in the art.

[0021] Those skilled in the art will recognize that embodiments of the present invention can be implemented as a system, apparatus, device, method, or computer program product. Therefore, this disclosure can be specifically implemented in the following forms: entirely hardware, entirely software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0022] According to an embodiment of the present invention, a method and apparatus for calculating the array gain of a deep-sea near-surface vertical array buoy are proposed.

[0023] It should be noted that the number of any elements in the accompanying drawings is for illustrative purposes only and not as a limitation, and any naming is for distinction only and has no limiting meaning.

[0024] The principles and spirit of the present invention will be explained in detail below with reference to several representative embodiments.

[0025] Exemplary methods The following is for reference. Figure 1 , Figure 1 This is a flowchart illustrating a method for calculating the array gain of a deep-sea near-surface vertical array buoy according to an embodiment of the present invention. It should be noted that the embodiments of the present invention can be applied to vertical array buoys pre-installed in the deep sea, where an infinite number of wind-generated noise sources are located.

[0026] Specifically, the wind-generated noise source consists of an infinite number of uncorrelated monopole sound sources, uniformly distributed across an infinitely large plane below the sea surface, with the noise source depth... λ is the wavelength (e.g., uniformly distributed at depths below the sea surface). On the plane, the radius of the noise source is 50km. The calculated sound source frequency is 100Hz, the sea depth is 4km, and the sound velocity profile is as follows. Figure 2 As shown, the speed of sound on the seabed is 2000 m / s, and the density of the seabed is 1.5. The seabed absorption coefficient is 1.2. The uppermost element of the 10-element vertical array buoy is element number 1, with a receiving depth of 30m and an element spacing of 7.5m. Assuming that the marine environment and noise source intensity in this area do not change with distance, the vertical direction... , Cross-spectral density function of sound pressure in noise field between two points for: ; in, Let be the depth between two points in the vertical direction, i.e., the depth between any two vertical array elements, and k be the wave number in the water. For noise source intensity, The horizontal distance of the noise source distribution. For the point source sound pressure Green's function, " indicates conjugate operation.

[0027] Figure 1 The flowchart of the array gain calculation method for deep-sea near-surface vertical array buoys provided in an embodiment of the present invention, shown below, includes: Step S101: Obtain relevant parameters of the vertical array buoy and noise source parameters.

[0028] In this embodiment, the relevant parameters of the vertical array buoy include at least the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix. The noise source parameters include the noise source intensity, the horizontal wavenumber, the maximum horizontal wavenumber, the minimum horizontal wavenumber, and the number of horizontal wavenumber samples.

[0029] In this embodiment of the application, the beam scanning elevation angle is taken as... The beam scanning interval is 1°, for a total of 181 angles.

[0030] Step S102: Calculate the first contribution value of near-field noise to the noise cross-spectral density based on the relevant parameters of the vertical array buoy and the noise source parameters.

[0031] In this embodiment, the location r of the vertical array is taken as the origin of the coordinate system. With the origin as the center and a as the radius, the source plane is divided into two parts: an inner circle and an outer circle. The noise source inside the circle is a near-field noise source with a radius of 3km, and the noise source outside the circle is a far-field noise source.

[0032] In this embodiment of the application, according to the fast field model, the formula for calculating the first contribution value of near-field noise to the noise cross-spectral density is as follows: ; ; ; in, This represents the first contribution of near-field noise to the noise cross-spectral density. and This represents the depth of any two elements in the vertical array buoy; The noise source intensity is represented by k; the horizontal wavenumber is represented by k. Indicates the maximum horizontal wavenumber. The minimum horizontal wavenumber is represented by N; the number of horizontal wavenumber samples is represented by a; the radius of the region corresponding to the near-field noise is represented by r; and the calculated horizontal distance is represented by r. Represents the Green's function based on the fast field model; This indicates the conjugate operation. For discrete wavenumbers.

[0033] For example, if a = 3km, The formula for calculating the first contribution of near-field noise to the noise cross-spectral density can be: ; Step S103: Calculate the second contribution value of far-field noise to noise cross-spectral density based on the vertical array buoy correlation parameters and the noise source parameters.

[0034] In this embodiment of the application, based on normal mode theory, considering seabed emission and absorption, and after ignoring and simplifying different order normal mode interference terms, the formula for calculating the second contribution value of far-field noise to the noise cross-spectral density is as follows: ; in, This represents the second contribution of far-field noise to the noise cross-spectral density, where m represents the m-th normal mode. The real part of the horizontal wavenumber of the m-th normal mode is given by [the real part of the horizontal wavenumber]. The imaginary part of the horizontal wavenumber of the m-th normal mode is given by [the term in Chinese]. Indicates the depth of the noise source. This represents the amplitude of the m-th normal mode excited by the wind-induced noise source. For depth The amplitude of the m-th normal mode. For the depth of the conjugate operation The amplitude of the m-th normal mode.

[0035] For example, if the radius of the noise source is 50km and a=3km, the formula for calculating the second contribution of far-field noise to the noise cross-spectral density can be: Step S104: Determine the noise cross-spectral density matrix based on the first contribution value and the second contribution value.

[0036] In this embodiment of the application, the formula for calculating the noise cross-spectral density matrix is ​​as follows: ; in, This represents the noise cross-spectral density matrix.

[0037] Step S105: Determine the noise directionality gain of the vertical array buoy based on the noise cross-spectral density matrix and the beam scanning elevation angle.

[0038] In this embodiment of the application, the formula for calculating the noise directionality gain of the vertical array buoy is as follows: ; ; in, This represents the noise directionality gain of the vertical array buoy. Indicates the beam scanning elevation angle, The guide vector represents the beam scanning elevation angle, and H represents the conjugate transpose operation. The noise cross-spectral density matrix represents the noise cross-spectral density matrix. This indicates the operating frequency of the vertical array buoys. This indicates the element spacing of the vertical array buoys. Indicates the number of beam scanning angles. This represents the speed of sound in water. This indicates the number of array elements of the vertical array buoy.

[0039] For example, if the sound source frequency is 100Hz, the element spacing is 7.5m, and the vertical array buoy has 10 elements, then the steering vector... It can be expressed as: ; Step S106: Determine the signal directional gain of the vertical array buoy based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix.

[0040] In this embodiment, the signal source depth is 400m. Considering the influence of the ocean channel on sound propagation, the formula for calculating the signal directivity gain of the vertical array buoy is as follows: ; ; in, p represents the signal directional gain of the vertical array buoy, and p represents the received signal. This represents the cross-spectral density matrix of the signal.

[0041] For example, if the number of array elements in a vertical array buoy is 10, then the formula for calculating the signal directivity gain of the vertical array buoy can be: ; Step S107: Based on the noise directivity gain and the signal directivity gain, determine the array directivity gain of the vertical array buoy, and then determine the array gain of the vertical array buoy.

[0042] In this embodiment of the application, the formula for calculating the array gain of the vertical array buoy is: ; in, This represents the array directional gain of the vertical array buoy. Therefore, the following is obtained: Array directional gain within the beam scanning range Beam scanning angle The change in the directional gain of the array is such that the maximum value of the directional gain of the array is the array gain AG of the vertical array.

[0043] Figure 3 This diagram illustrates a comparison between the noise vertical directivity of a deep-sea near-surface vertical array buoy calculation method and the noise vertical directivity based on a normal mode model, according to an embodiment of the present invention. At small angles, the noise primarily originates from the horizontal direction, while at large angles, it mainly originates from overhead noise from the sea surface and seabed. Since the vertical array buoy primarily receives noise from large angles in the near field, the calculation error of the vertical directivity based on the normal mode model is greater at large angles. This is because the normal mode model only considers discrete normal mode modes, neglecting the contribution of the continuous spectrum. The continuous spectrum has a significant impact on noise incident at large grazing angles, and this energy attenuates faster with distance than discrete normal modes. Therefore, when the vertical array buoy beam pointing at large angles, the calculation results of the noise spatial correlation characteristics based on the normal mode model have a large error. The noise spatial correlation characteristic calculation model based on near-field fast field-far-field normal modes used in this invention yields more accurate results.

[0044] Figure 4 a is a schematic diagram of the array directional gain of a vertical array buoy, which is a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention. Figure 4 b is a schematic diagram of the array directional gain of a vertical array buoy based on a normal wave model provided in an embodiment of the present invention; Figure 4 c is a schematic diagram showing the difference between the array directional gain and the array directional gain based on the normal mode model in a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention.

[0045] Specifically, it can be observed that when the beam pointing angle of the vertical array buoy is small, the array directivity gain of the two methods is not significantly different. When the beam pointing angle is greater than 60° or less than -60°, the error of the array directivity gain of the vertical array buoy based on the normal mode model increases basically with the increase of the absolute value of the beam pointing angle, which is consistent with the changing law of the vertical directivity difference between the two methods. The maximum error of the array directivity gain of the vertical array buoy based on the normal mode model reaches 5.8dB.

[0046] Figure 5This is a schematic diagram illustrating the change in the angle of arrival of the direct sound ray, the sound ray reflected from the sea surface, and the sound ray reflected from the first seabed at array element 1 of the present invention as a function of horizontal distance. When the horizontal distance is extremely close, the angles of arrival of the direct sound ray, the sound ray reflected from the sea surface, and the sound ray reflected from the first seabed are all greater than 60° or less than -60°. At this time, the noise received by the array is mainly overhead noise incident from a large angle direction. The noise directionality gain based on the normal mode model has a large error. As the horizontal distance increases, the angle of arrival of the sound ray gradually decreases, and the noise directionality gain error based on the normal mode model also gradually decreases.

[0047] Figure 6 a is a schematic diagram comparing the noise gain of a method for calculating the array gain of a deep-sea near-shore vertical array buoy according to an embodiment of the present invention with the noise gain based on a normal mode model. Figure 6 b is a schematic diagram comparing the array gain of a deep-sea near-shore vertical array buoy with that of an array gain calculation method based on a normal mode model, according to an embodiment of the present invention.

[0048] Specifically, the smaller graphs in the figure represent the noise gain and array gain at horizontal distances of 3.8km-4.8km. When the horizontal distance between the transmitter and receiver is less than 0.4km, the noise gain based on normal mode theory has a significant error compared to the method of this invention. This is because the main contributing sound rays are incident at large angles, and the beam control direction of the vertical array buoy is also at a large angle. Therefore, the calculation of the spatial correlation characteristics of noise based on the normal mode model has an error, leading to an error in the noise gain. As the horizontal distance increases, the angle of arrival of the main contributing sound rays gradually decreases, and the calculation error of the spatial correlation characteristics of noise based on the normal mode model also gradually decreases. Therefore, the noise gain is only slightly greater than that of the method of this invention. When the horizontal distance is greater than 3.8km, at some horizontal distance points, the energy of the first seabed reflected sound ray cannot be ignored, and its angle of arrival is greater than 55°. Therefore, the noise gain based on the normal mode model has an error. In summary, at distances less than 0.4 km and some distances greater than 3.8 km, the noise gain based on the normal mode model has significant errors, which in turn leads to errors in the array gain. The maximum error in the array gain is 5.5 dB. Therefore, compared with the array gain calculation method based on the normal mode model, the method of this invention can achieve accurate prediction of the near-range array gain of deep-sea near-surface vertical array buoys.

[0049] This invention obtains relevant parameters of the vertical array buoy and noise source parameters, calculates the contributions of near-field noise and far-field noise to the noise cross-spectral density, and then determines the noise cross-spectral density matrix. This comprehensive approach, considering the influence of near-field and far-field noise, makes noise correlation calculations more accurate. Simultaneously, it combines multiple parameters to determine the noise directional gain and signal directional gain, comprehensively considering factors such as beam scanning elevation angle, signal cross-spectral density matrix, and received signal. Finally, the array gain is determined based on the noise and signal directional gains. The entire calculation process fully considers various practical factors, improving the accuracy of array gain calculation.

[0050] Exemplary device After introducing the method of exemplary embodiments of the present invention, the following references are made. Figure 7 An exemplary embodiment of the present invention describes a gain calculation device for a deep-sea near-surface vertical array buoy. This device is applied to a vertical array buoy pre-installed in the deep sea, wherein an infinite number of wind-generated noise sources are located within the sea area where the vertical array buoy is situated, and all noise sources are at the same depth. The device includes: The acquisition unit 701 is used to acquire relevant parameters of the vertical array buoy and noise source parameters; wherein, the relevant parameters of the vertical array buoy include at least the beam scanning elevation angle, the signal cross-spectral density matrix, and the received signal; The first calculation unit 702 is used to calculate the first contribution value of near-field noise to noise cross-spectral density based on the relevant parameters of the vertical array buoy and the noise source parameters. The second calculation unit 703 is used to calculate the second contribution value of far-field noise to noise cross-spectral density based on the vertical array buoy correlation parameters and the noise source parameters. The first determining unit 704 is used to determine the noise cross-spectral density matrix based on the first contribution value and the second contribution value; The second determining unit 705 is used to determine the noise directivity gain of the vertical array buoy based on the noise cross-spectral density matrix and the beam scanning elevation angle. The third determining unit 706 is used to determine the signal directivity gain of the vertical array buoy based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix. The fourth determining unit 707 is used to determine the array directional gain of the vertical array buoy based on the noise directional gain and the signal directional gain, and then determine the array gain of the vertical array buoy.

[0051] The noise source parameters include noise source intensity, horizontal wavenumber, maximum horizontal wavenumber, minimum horizontal wavenumber, and the number of horizontal wavenumber samples; the formula for calculating the first contribution value of near-field noise to the noise cross-spectral density is: ; ; ; in, This represents the first contribution of near-field noise to the noise cross-spectral density. and This represents the depth of any two elements in the vertical array buoy; The noise source intensity is represented by k; the horizontal wavenumber is represented by k. Indicates the maximum horizontal wavenumber. The minimum horizontal wavenumber is represented by N; the number of horizontal wavenumber samples is represented by a; the radius of the region corresponding to the near-field noise is represented by r; and the calculated horizontal distance is represented by r. Represents the Green's function based on the fast field model; This indicates the conjugate operation. For discrete wavenumbers.

[0052] In this embodiment, the formula for calculating the second contribution value of far-field noise to the noise cross-spectral density is as follows: ; in, This represents the second contribution of far-field noise to the noise cross-spectral density, where m represents the m-th normal mode. The real part of the horizontal wavenumber of the m-th normal mode is given by [the real part of the horizontal wavenumber]. The imaginary part of the horizontal wavenumber of the m-th normal mode is given by [the term in Chinese]. Indicates the depth of the noise source. This represents the amplitude of the m-th normal mode excited by the wind-induced noise source. For depth The amplitude of the m-th normal mode. For the depth of the conjugate operation The amplitude of the m-th normal mode.

[0053] In this embodiment of the application, the formula for calculating the noise cross-spectral density matrix is ​​as follows: ; in, This represents the noise cross-spectral density matrix.

[0054] In this embodiment of the application, the formula for calculating the noise directionality gain of the vertical array buoy is as follows: ; ; in, This represents the noise directionality gain of the vertical array buoy. Indicates the beam scanning elevation angle, The guide vector represents the beam scanning elevation angle, and H represents the conjugate transpose operation. The noise cross-spectral density matrix represents the noise cross-spectral density matrix. This indicates the operating frequency of the vertical array buoys. This indicates the element spacing of the vertical array buoys. Indicates the number of beam scanning angles. This represents the speed of sound in water. This indicates the number of array elements of the vertical array buoy.

[0055] In this embodiment of the application, the formula for calculating the signal directivity gain of the vertical array buoy is as follows: ; ; in, p represents the signal directional gain of the vertical array buoy, and p represents the received signal. This represents the cross-spectral density matrix of the signal.

[0056] In this embodiment of the application, the formula for calculating the array gain of the vertical array buoy is: ; in, This represents the array directivity gain of the vertical array buoy. Within the beam scanning range, the array directivity gain... The maximum value is the array gain of the vertical array. .

[0057] This invention obtains relevant parameters of the vertical array buoy and noise source parameters, calculates the contributions of near-field noise and far-field noise to the noise cross-spectral density, and then determines the noise cross-spectral density matrix. This comprehensive approach, considering the influence of near-field and far-field noise, makes noise correlation calculations more accurate. Simultaneously, it combines multiple parameters to determine the noise directional gain and signal directional gain, comprehensively considering factors such as beam scanning elevation angle, received signal, and signal cross-spectral density matrix. Finally, the array gain is determined based on the noise directional gain and signal directional gain. The entire calculation process fully considers various practical factors, improving the accuracy of array gain calculation.

[0058] Exemplary media After introducing the methods and apparatus of exemplary embodiments of the present invention, the following references are made. Figure 8 A computer-readable storage medium according to exemplary embodiments of the present invention will be described, please refer to... Figure 8The computer-readable storage medium shown is an optical disc 80, on which a computer program (i.e., a program product) is stored. When the computer program is run by a processor, it performs the steps described in the above method implementation, such as acquiring the relevant parameters of the vertical array buoy and the noise source parameters; wherein the relevant parameters include at least the beam scanning elevation angle, the signal cross-spectral density matrix, and the received signal; calculating the first contribution value of near-field noise to the noise cross-spectral density based on the relevant parameters of the vertical array buoy and the noise source parameters; and calculating the far-field noise based on the relevant parameters of the vertical array buoy and the noise source parameters. The process involves: determining a second contribution value to the noise cross-spectral density; determining a noise cross-spectral density matrix based on the first and second contribution values; determining the noise directivity gain of the vertical array buoy based on the noise cross-spectral density matrix and the beam scanning elevation angle; determining the signal directivity gain of the vertical array buoy based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix; determining the array directivity gain of the vertical array buoy based on the noise directivity gain and the signal directivity gain, and further determining the array gain of the vertical array buoy. The specific implementation methods of each step will not be repeated here. It should be noted that examples of the computer-readable storage medium may also include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory, or other optical and magnetic storage media, which will not be elaborated upon here.

[0059] Exemplary computing device After introducing the methods, apparatus, and media of exemplary embodiments of the present invention, the following references are made. Figure 9 A computing device for calculating the array gain of a deep-sea near-shore vertical array buoy according to an exemplary embodiment of the present invention.

[0060] Figure 9 A block diagram is shown of an exemplary computing device 90 suitable for implementing embodiments of the present invention, which may be a computer system or a server. Figure 9 The computing device 90 shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of the present invention.

[0061] like Figure 9 As shown, the components of the computing device 90 may include, but are not limited to: one or more processors or processing units 901, system memory 902, and bus 903 connecting different system components (including system memory 902 and processing unit 901).

[0062] The computing device 90 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the computing device 90, including volatile and non-volatile media, removable and non-removable media.

[0063] System memory 902 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 9021 and / or cache memory 9022. Computing device 90 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, ROM 9023 may be used to read and write non-removable, non-volatile magnetic media (…). Figure 9 (Not shown in the image, usually referred to as "hard drive"). Although not shown in... Figure 9 The diagram illustrates that disk drives for reading and writing to removable non-volatile disks (e.g., "floppy disks") and optical disc drives for reading and writing to removable non-volatile optical discs (e.g., CD-ROMs, DVD-ROMs, or other optical media) can be provided. In these cases, each drive can be connected to bus 903 via one or more data media interfaces. System memory 902 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of the embodiments of the present invention.

[0064] A program / utility 9025 having a set (at least one) of program modules 9024 may be stored, for example, in system memory 902, and such program modules 9024 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment. Program modules 9024 typically perform the functions and / or methods described in the embodiments of the present invention.

[0065] The computing device 90 can also communicate with one or more external devices 904 (such as a keyboard, pointing device, display, etc.). This communication can be performed through the input / output (I / O) interface 905. Furthermore, the computing device 90 can also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) via a network adapter 906. Figure 9 As shown, network adapter 906 communicates with other modules of computing device 90 (such as processing unit 901) via bus 903. It should be understood that, although... Figure 9 As not shown, it can be used in conjunction with computing device 90 with other hardware and / or software modules.

[0066] The processing unit 901 executes various functional applications and data processing by running programs stored in the system memory 902. For example, it acquires the relevant parameters of the vertical array buoy and the noise source parameters. The relevant parameters of the vertical array buoy include at least the beam scanning elevation angle, the signal cross-spectral density matrix, the received signal, and the buoy coordinates. Based on the relevant parameters of the vertical array buoy and the noise source parameters, it calculates a first contribution value of near-field noise to the noise cross-spectral density. Based on the relevant parameters of the vertical array buoy and the noise source parameters, it calculates a second contribution value of far-field noise to the noise cross-spectral density. Based on the first contribution value and the second contribution value, it determines the noise cross-spectral density matrix. Based on the noise cross-spectral density matrix and the beam scanning elevation angle, it determines the noise directivity gain of the vertical array buoy. Based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix, it determines the signal directivity gain of the vertical array buoy. Based on the noise directivity gain and the signal directivity gain, it determines the array directivity gain of the vertical array buoy, and thus determines the array gain of the vertical array buoy. The specific implementation methods of each step will not be repeated here. It should be noted that although several units / modules or sub-units / sub-modules of the array gain calculation device for deep-sea near-surface vertical array buoys have been mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, according to embodiments of the present invention, the features and functions of two or more units / modules described above can be embodied in one unit / module. Conversely, the features and functions of one unit / module described above can be further divided and embodied by multiple units / modules.

[0067] In the description of this invention, it should be noted that the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0068] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0069] In the several embodiments provided by this invention, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0070] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0071] In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0072] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this invention, essentially, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0073] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

[0074] Furthermore, although the operations of the method of the present invention are described in a specific order in the accompanying drawings, this does not require or imply that these operations must be performed in that specific order, or that all the operations shown must be performed to achieve the desired result. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

[0075] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

Claims

1. A method for calculating the array gain of a deep-sea near-surface vertical array buoy array, characterized in that, The method is applied to a vertical array buoy previously arranged in a deep sea, and infinite wind-generated noise sources are arranged in a sea area where the vertical array buoy is located, and the method comprises the following steps: Obtaining vertical array buoy related parameters and noise source parameters; wherein the vertical array buoy related parameters at least include a beam scanning elevation angle, a signal cross-spectral density matrix, and a received signal; Based on the vertical array buoy related parameters and the noise source parameters, a first contribution value of near-field noise to noise cross-spectral density is calculated; Based on the vertical array buoy related parameters and the noise source parameters, a second contribution value of far-field noise to noise cross-spectral density is calculated; Based on the first contribution value and the second contribution value, a noise cross-spectral density matrix is determined; Based on the noise cross-spectral density matrix and the beam scanning elevation angle, a noise directivity gain of the vertical array buoy is determined; Based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix, a signal directivity gain of the vertical array buoy is determined; Based on the noise directivity gain and the signal directivity gain, a base array directivity gain of the vertical array buoy is determined, and then an array gain of the vertical array buoy is determined.

2. The method of claim 1, wherein, The noise source parameters include noise source intensity, horizontal wave number, maximum horizontal wave number, minimum horizontal wave number, and the number of horizontal wave number samples; the calculation formula of the first contribution value of near-field noise to noise cross-spectral density is: ； ； ； wherein, represents a first contribution value of the near-field noise to the noise cross-spectral density, and represents the depth of any two elements in the vertical array buoy; represents the noise source strength; k represents the horizontal wave number; represents the maximum horizontal wave number, represents the minimum horizontal wave number; N represents the sampling number of the horizontal wave number, a represents the radius of the region corresponding to the near-field noise; r represents the calculated horizontal distance; represents the Green function based on the fast-field model; represents the conjugate operation; is a discrete wave number.

3. The method of claim 2, wherein, The calculation formula of the second contribution value of far-field noise to noise cross-spectral density is: ； wherein, represents a second contribution value of the far-field noise to the noise cross-spectral density, m represents the mth order normal mode, represents a real part of the horizontal wave number of the mth order normal mode, represents an imaginary part of the horizontal wave number of the mth order normal mode, represents a depth of the noise source, represents an amplitude of the mth order normal mode excited by the wind-generated noise source, is an amplitude of the mth order normal mode at a depth is an amplitude of the mth order normal mode at a depth which has been subjected to a conjugate operation.​ 4. The method of claim 3, wherein, The calculation formula of the noise cross-spectral density matrix is: ； wherein denotes the noise cross-spectral density matrix.

5. The method of claim 4, wherein, The calculation formula of the noise directivity gain of the vertical array buoy is: ； ； wherein, represents the noise directivity gain of the vertical array buoy, represents the beam scanning elevation angle, represents the steering vector of the beam scanning elevation angle, H represents the conjugate transpose operation, represents the noise cross-spectral density matrix, represents the operating frequency of the vertical array buoy, represents the element spacing of the vertical array buoy, represents the number of beam scanning angles, represents the sound speed in water, represents the number of elements of the vertical array buoy.

6. The method of claim 5, wherein, The calculation formula of the signal directivity gain of the vertical array buoy is: ； ； wherein represents the signal directivity gain of the vertical array of buoys, p represents the received signal, represents the signal cross-spectral density matrix.

7. The method of claim 6, wherein, The calculation formula of the base array directivity gain of the vertical array buoy is: ； wherein, represents the base array directivity gain of the vertical array buoy, and the base array directivity gain of the vertical array buoy has a maximum value of the base array gain of the vertical array buoy .

8. A device for calculating the array gain of a deep-sea near-surface vertical array buoy, characterized in that, The device is applied to a vertical array buoy previously arranged in a deep sea, and wind-generated noise sources are arranged in a sea area where the vertical array buoy is located, and the device comprises: An obtaining unit is configured to obtain vertical array buoy related parameters and noise source parameters; wherein the vertical array buoy related parameters at least include a beam scanning elevation angle, a signal cross-spectral density matrix, a received signal, and a buoy coordinate; A first calculation unit is configured to calculate a first contribution value of near-field noise to noise cross-spectral density based on the vertical array buoy related parameters and the noise source parameters; A second calculation unit is configured to calculate a second contribution value of far-field noise to noise cross-spectral density based on the vertical array buoy related parameters and the noise source parameters; A first determination unit is configured to determine a noise cross-spectral density matrix based on the first contribution value and the second contribution value; A second determination unit is configured to determine a noise directivity gain of the vertical array buoy based on the noise cross-spectral density matrix and the beam scanning elevation angle; A third determination unit is configured to determine a signal directivity gain of the vertical array buoy based on the beam scanning elevation angle, the received signal, and the signal cross-spectral density matrix; A fourth determination unit is configured to determine a base array directivity gain of the vertical array buoy based on the noise directivity gain and the signal directivity gain, and then determine an array gain of the vertical array buoy.

9. The method of claim 8, wherein, The noise source parameters include noise source intensity, horizontal wave number, maximum horizontal wave number, minimum horizontal wave number and sampling number of horizontal wave number; and the calculation formula of the first contribution value of the near-field noise to the noise cross-spectral density is: ； ； ； wherein, represents a first contribution value of the near-field noise to the noise cross-spectral density, and represents a depth of any two array elements in the vertical array buoy; represents a noise source strength; k represents a horizontal wave number; represents a maximum horizontal wave number, represents a minimum horizontal wave number; N represents a sampling number of horizontal wave numbers, a represents a radius of a region corresponding to the near-field noise; and r represents a horizontal distance of calculation; represents a Green function based on a fast-field model; represents a conjugate operation, is a discrete wave number.

10. The method of claim 9, wherein, The calculation formula of the second contribution value of the far-field noise to the noise cross-spectral density is: ； wherein, represents a second contribution value of the far-field noise to the noise cross-spectral density, m represents the mth order normal mode, represents a real part of the horizontal wave number of the mth order normal mode, represents an imaginary part of the horizontal wave number of the mth order normal mode, represents a depth of the noise source, represents an amplitude of the mth order normal mode excited by the wind-generated noise source, is an amplitude of the mth order normal mode at a depth of the mth order normal mode at a depth is an amplitude of the mth order normal mode at a depth of the mth order normal mode at a depth

Images