Underwater acoustic reverberation suppression stationary phase nuclear inversion parameterized phase synchronization transformation method, device and equipment and storage medium

Abstract

The invention discloses an underwater acoustic reverberation suppression stationary phase nuclear inversion parameterized phase synchronization transformation method, device and equipment and a storage medium, and the method comprises the steps: determining the frequency spectrum power information and phase function of a transmitted underwater acoustic signal, carrying out the windowing processing of the received underwater acoustic signal based on the phase alignment information and phase translation information, and obtaining the phase function of the received underwater acoustic signal; performing parameterized phase synchronization frequency-time transformation on each receiving window signal to obtain a standard receiving window signal in phase synchronization with each transmitting window signal; transmitting time-frequency distribution data and receiving time-frequency distribution data are respectively determined, and kernel inversion iteration correction is carried out on each standard receiving window signal; according to the invention, the accuracy of the determination result of the target underwater acoustic signal can be improved, so that the reverberation interference is effectively suppressed, and the detectability of the weak target signal is improved.

Description

Technical Field

[0001] This application belongs to the field of underwater acoustic communication technology, specifically relating to underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation method, apparatus, equipment, and storage medium. Background Technology

[0002] In underwater acoustic communication applications such as active sonar detection, reverberation is an unavoidable strong background interference. Reverberation is mainly formed by the scattering and superposition of emitted sound waves during propagation by the sea surface, seabed, and inhomogeneous bodies within the water. Against a strong reverberation background, the echo of weak target underwater acoustic signals is masked, leading to a significant decrease in the detection performance of target underwater acoustic signals. Therefore, achieving effective signal detection and feature extraction of target underwater acoustic signals has become one of the core challenges in the field of underwater acoustic signal processing.

[0003] In existing technologies, the determination of target underwater acoustic signals mainly involves using short-time Fourier transform to window and segment the received signal to obtain its time-frequency distribution, or designing parameterized time-frequency analysis algorithms with higher time-frequency convergence to determine the target underwater acoustic signal in the received signal. However, existing technologies suffer from low time-frequency resolution and inaccurate target underwater acoustic signal determination results. Summary of the Invention

[0004] The purpose of this application is to provide a method, apparatus, device, and storage medium for underwater acoustic reverberation suppression, stationary phase kernel inversion, and parameterized phase synchronization transformation. This solves the problems of low time-frequency resolution and inaccurate target underwater acoustic signal determination in existing technologies. By performing phase calibration on the signals of each receiving window using phase alignment and phase shift information, a standard receiving window signal synchronized with the phase of each transmitting window signal is obtained. Based on the transmitting and receiving time-frequency distribution data, kernel inversion iterative correction is performed on each standard receiving window signal, and the iterative correction results of each standard receiving window signal are spliced ​​together to obtain the target underwater acoustic signal after reverberation suppression. This achieves the goal of determining the target underwater acoustic signal in the received underwater acoustic signal based on the original characteristics of the transmitted underwater acoustic signal. It effectively overcomes the problems of low time-frequency resolution and large sidelobe interference in conventional time-frequency analysis, significantly suppressing reverberation and other interference signals in the received underwater acoustic signal, improving the time-frequency convergence and detectability of weak target signals, and thus improving the accuracy of the target underwater acoustic signal determination result.

[0005] In a first aspect, embodiments of this application provide a method for underwater acoustic reverberation suppression, stationary phase kernel inversion, and parameterized phase synchronization transformation, the method comprising: The spectral power information of the transmitted underwater acoustic signal is determined based on the preset stationary phase algorithm. The phase function of the transmitted underwater acoustic signal is determined based on the spectral power information. The phase function is then inverted to obtain the phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters. Calculate the phase shift information corresponding to the transmitted underwater acoustic signal based on the preset phase alignment equation and phase alignment information; Based on the preset window center frequency and preset window splitting length, the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed to obtain multiple received window signals. Then, based on the phase alignment information and phase shift information, parameterized phase synchronization frequency-time transformation is performed on each received window signal to obtain a standard received window signal that is phase-synchronized with each transmitted window signal. The standard received window signal is a time-frequency signal. The transmission time-frequency distribution data corresponding to each transmission window signal and the reception time-frequency distribution data corresponding to each standard reception window signal are determined respectively. Based on the transmission time-frequency distribution data and the reception time-frequency distribution data, the kernel inversion iterative correction is performed on each standard reception window signal until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than the preset error threshold. The target underwater acoustic signal after reverberation suppression is obtained by splicing the iterative correction results of the signals from each standard receiving window.

[0006] Furthermore, based on the transmit and receive time-frequency distribution data, the signals of each standard receive window are iteratively corrected using kernel inversion, including: Based on the transmission time-frequency distribution data, the kernel information of the transmitted underwater acoustic signal inversion corresponding to each transmission window signal is determined, and the observation information of the received underwater acoustic signal corresponding to each standard reception window signal is determined based on the reception time-frequency distribution data. Based on the correspondence between each transmission window signal and each standard receiving window signal, the kernel information of the transmitted underwater acoustic signal inversion and the observation information of the received underwater acoustic signal, the kernel inversion iterative correction is performed on each standard receiving window signal to obtain the iterative correction result corresponding to each standard receiving window signal.

[0007] Furthermore, based on the correspondence between each transmission window signal and each standard receiving window signal, the kernel information of the transmitted underwater acoustic signal inversion, and the observation information of the received underwater acoustic signal, the kernel inversion iterative correction of each standard receiving window signal is performed, including: Based on the correspondence between each transmission window signal and each standard receiving window signal, the target transmission underwater acoustic signal inversion kernel information corresponding to the received underwater acoustic signal observation information of each standard receiving window signal is determined in multiple transmission underwater acoustic signal inversion kernel information; Based on the preset iterative algorithm, the received underwater acoustic signal observation information, and the target emitted underwater acoustic signal inversion kernel information, the kernel inversion iterative correction is performed on the signals of each standard receiving window.

[0008] Furthermore, the formula for the preset iterative algorithm is: ; in, It is the first The first of the standard received window signals The result of the second iteration It is the first The i-th iteration result of the standard receive window signal. It is the first The received underwater acoustic signal observation information of the standard receiving window signal. It is the target transmitted underwater acoustic signal inversion kernel information for each of the aforementioned standard receiving window signals. yes The time-frequency distribution pattern, It is a frequency dimension variable. It is a time-dimensional variable.

[0009] Furthermore, based on a preset window center frequency and a preset window splitting length, the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed to obtain multiple received window signals, including: Fourier transforms are performed on the transmitted underwater acoustic signal and the corresponding received underwater acoustic signal to obtain the transmitted spectrum information and the received spectrum information, respectively. The received spectrum information is windowed according to the preset window center frequency and preset window split length to obtain multiple receiving window information. There is an overlap between adjacent receiving windows corresponding to multiple receiving window information. Multiple receiving window signals corresponding to the transmitted underwater acoustic signal are determined based on information from multiple receiving windows.

[0010] Furthermore, by splicing the iterative correction results of the signals from each standard receiving window, the target underwater acoustic signal after reverberation suppression is obtained, including: The iterative correction results of each standard receiving window signal are spliced ​​together according to the window frequency band sequence number corresponding to each standard receiving window signal to obtain the full-band target underwater acoustic signal after underwater acoustic reverberation suppression corresponding to the transmitted underwater acoustic signal. The full-band target underwater acoustic signal is a two-dimensional matrix including frequency and time dimensions.

[0011] Secondly, embodiments of this application provide a device for underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation, the device comprising: The power spectrum information extraction module is used to determine the spectral power information of the transmitted underwater acoustic signal based on a preset stationary phase algorithm; The phase alignment information determination module is used to determine the phase function of the transmitted underwater acoustic signal based on the spectral power information, and to invert the phase function to obtain the phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters; The phase shift information determination module is used to calculate the phase shift information corresponding to the transmitted underwater acoustic signal based on the preset phase alignment equation and phase alignment information. The windowed signal alignment module is used to perform windowed processing on the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal based on the preset window center frequency and preset window splitting length, to obtain multiple received window signals. Based on phase alignment information and phase shift information, parameterized phase synchronization frequency-time transformation is performed on each received window signal to obtain a standard received window signal that is phase-synchronized with each transmitted window signal. The standard received window signal is a time-frequency signal. The underwater acoustic signal correction module is used to determine the transmission time-frequency distribution data corresponding to each transmission window signal and the reception time-frequency distribution data corresponding to each standard reception window signal, and to perform kernel inversion iterative correction on each standard reception window signal based on the transmission time-frequency distribution data and the reception time-frequency distribution data until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than a preset error threshold. The target underwater acoustic signal module is used to stitch together the iterative correction results of the signals from each standard receiving window to obtain the target underwater acoustic signal after underwater reverberation suppression.

[0012] Thirdly, embodiments of this application provide an electronic device including a processor, a memory, and a program or instructions stored in the memory and executable on the processor, wherein the program or instructions, when executed by the processor, implement the steps of the method described in the first aspect.

[0013] Fourthly, embodiments of this application provide a readable storage medium on which a program or instructions are stored, which, when executed by a processor, implement the steps of the method described in the first aspect.

[0014] Fifthly, embodiments of this application also provide a computer program product comprising a computer program stored in a computer-readable storage medium, wherein at least one processor of the device reads from the computer-readable storage medium and executes the computer program, causing the device to perform the method described in the first aspect.

[0015] In this embodiment, the spectral power information of the transmitted underwater acoustic signal is extracted, the phase information of the transmitted underwater acoustic signal is determined based on the spectral power information, and the phase alignment information and phase shift information corresponding to the transmitted underwater acoustic signal are determined based on the phase information. The received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed based on a preset window center frequency and a preset window splitting length to obtain multiple received window signals. The parameterized phase synchronization frequency-time transformation of each received window signal is performed based on the phase alignment information and phase shift information to obtain a standard received window signal that is phase-synchronized with each transmitted window signal. The standard received window signal is a time-frequency signal. The transmitted time-frequency distribution data corresponding to each transmitted window signal and the received time-frequency distribution data corresponding to each standard received window signal are determined respectively. The kernel inversion iterative correction of each standard received window signal is performed based on the transmitted time-frequency distribution data and the received time-frequency distribution data until the mean square error between the correction result of each standard received window signal and the corresponding received time-frequency distribution data is less than a preset error threshold. The iterative correction results of each standard received window signal are spliced ​​together to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression. The above-described method for suppressing underwater acoustic reverberation by stationary phase kernel inversion parameterized phase synchronization transformation solves the problems of low time-frequency resolution and inaccurate target underwater acoustic signal determination in existing technologies. Phase alignment and phase shift information are used to calibrate the phase of each receiving window signal, resulting in a standard receiving window signal that is phase-synchronized with each transmitting window signal. Kernel inversion iterative correction is performed on each standard receiving window signal based on the transmitting and receiving time-frequency distribution data. The iterative correction results of each standard receiving window signal are then combined to obtain the target underwater acoustic signal after reverberation suppression. This method achieves the goal of determining the target underwater acoustic signal in the received underwater acoustic signal based on the original characteristics of the transmitted underwater acoustic signal. It effectively overcomes the problems of low time-frequency resolution and large sidelobe interference in conventional time-frequency analysis, significantly suppresses reverberation and other interference signals in the received underwater acoustic signal, improves the time-frequency aggregation and detectability of weak target signals, and thus improves the accuracy of the target underwater acoustic signal determination. Attached Figure Description

[0016] Figure 1 This is a flowchart of a method for underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation, provided in an embodiment of this application. Figure 2 This is a schematic diagram illustrating the effect of short-time Fourier transform processing of hyperbolic frequency modulated signals provided in the embodiments of this application; Figure 3 This is a schematic diagram illustrating the effect of hyperbolic frequency modulation signal kernel inversion parameterized phase synchronization frequency-time transformation processing provided in the embodiments of this application; Figure 4 This is a comparative schematic diagram of the hyperbolic frequency modulation signal processing results provided in the embodiments of this application; Figure 5This is a flowchart of the iterative correction standard receiving window signal provided in the embodiments of this application; Figure 6 This is a flowchart of another underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method provided in this application embodiment; Figure 7 This is a structural block diagram of a hydroacoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation device provided in an embodiment of this application; Figure 8 This is a structural block diagram of the electronic device provided in the embodiments of this application. Detailed Implementation

[0017] To make the objectives, technical solutions, and advantages of this application clearer, specific embodiments of this application are described in detail below with reference to the accompanying drawings. It is understood that the specific embodiments described herein are merely for explaining this application and not for limiting it. Furthermore, it should be noted that, for ease of description, only the parts relevant to this application are shown in the drawings, not all of them. Before discussing exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although the flowcharts describe operations (or steps) as sequential processes, many of these operations can be performed in parallel, concurrently, or simultaneously. Furthermore, the order of the operations can be rearranged. The process can be terminated when its operation is completed, but may also have additional steps not included in the drawings. The process can correspond to a method, information, procedure, subroutine, subprogram, etc.

[0018] The technical solutions of the embodiments of this application will be clearly described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this application. All other embodiments obtained by those skilled in the art based on the embodiments of this application are within the scope of protection of this application.

[0019] The terms "first," "second," etc., used in the specification and claims of this application are used to distinguish similar objects and not to describe a specific order or sequence. It should be understood that such use of data can be interchanged where appropriate so that embodiments of this application can be implemented in orders other than those illustrated or described herein, and the objects distinguished by "first," "second," etc., are generally of the same class and the number of objects is not limited; for example, a first object can be one or more. Furthermore, in the specification and claims, "and / or" indicates at least one of the connected objects, and the character " / " generally indicates that the preceding and following objects are in an "or" relationship.

[0020] Firstly, this solution can be applied to underwater acoustic signal processing scenarios, particularly those involving the processing of broadband frequency-modulated underwater acoustic signals emitted by sonar. It is suitable for determining target underwater acoustic signals in towed arrays, stereo arrays, area arrays, single hydrophones, and both deep-sea and shallow-sea environments. Phase alignment and phase shift information are used to calibrate the phase of each receiving window signal, resulting in a standard receiving window signal synchronized with each transmitting window signal. Kernel inversion iterative correction is performed on each standard receiving window signal based on the transmitting and receiving time-frequency distribution data. The iterative correction results of each standard receiving window signal are then combined to obtain the target underwater acoustic signal after reverberation suppression. This achieves the goal of determining the target underwater acoustic signal in the received underwater acoustic signal based on the original characteristics of the transmitted underwater acoustic signal. It effectively overcomes the problems of low time-frequency resolution and large sidelobe interference in conventional time-frequency analysis, significantly suppressing reverberation and other interference signals in the received underwater acoustic signal, improving the time-frequency convergence and detectability of weak target signals, and thus enhancing the accuracy of the target underwater acoustic signal determination. Based on the above usage scenarios, it is understood that the execution entity of each step in this solution can be a computer device. The computer device refers to any electronic device with data computing, processing and storage capabilities, such as mobile phones, PCs (Personal Computers), tablet computers and other terminal devices, or it can be a server or other devices. This application embodiment does not limit this.

[0021] The following description, in conjunction with the accompanying drawings, details a method, apparatus, device, and storage medium for underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation, and related embodiments and application scenarios provided in this application.

[0022] Figure 1 This is a flowchart of a method for underwater acoustic reverberation suppression, stationary phase kernel inversion, and parameterized phase synchronization transformation provided in an embodiment of this application. Figure 1 As shown, the specific steps include the following: S101 determines the spectral power information of the transmitted underwater acoustic signal based on a preset stationary phase algorithm, determines the phase function of the transmitted underwater acoustic signal based on the spectral power information, and inverts the phase function to obtain the phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters.

[0023] Among them, the transmitted underwater acoustic signal can be a frequency-modulated signal with a known waveform emitted underwater by a sonar transducer or other device.

[0024] The expression for transmitting underwater acoustic signals is: ; in, It transmits underwater acoustic signals. It is the preset carrier amplitude. It is a preset carrier frequency. It is a modulated signal. It refers to frequency modulation sensitivity.

[0025] Phase alignment information can be information representing the core phase characteristics in the spectral power information of the transmitted underwater acoustic signal. The preset stationary phase algorithm can be a pre-set algorithm used to calculate the spectral power of the transmitted underwater acoustic signal. In this scheme, the spectral power information can be a spectral power function.

[0026] In one embodiment, the spectral power function of the transmitted underwater acoustic signal can be calculated by a preset stationary phase algorithm, and the phase function in the spectral power function can be identified. The phase function is then inverted to obtain the phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters.

[0027] The emitted underwater acoustic signal can be a phase-transmitted frequency-modulated signal, and the phase-transmitted frequency-modulated signal can be approximated using the preset stationary phase method. The spectral power function of the emitted underwater acoustic signal is obtained. .

[0028] Spectral power function The solution formula is: ; in, It is the phase function for transmitting underwater acoustic signals. It is the first derivative of the phase function of the emitted underwater acoustic signal with respect to time. It is the second derivative of the phase function of the emitted underwater acoustic signal with respect to time. Transmitting underwater acoustic signals at the station Amplitude at that point It is a complex exponential phase term, describing the signal at frequency... Phase characteristics at the location, It is a phase correction term used to compensate for phase deviations caused by the positive or negative sign of the second derivative of the phase.

[0029] S102, calculates the phase shift information corresponding to the emitted underwater acoustic signal based on the preset phase alignment equation and phase alignment information.

[0030] The phase shift information can represent the frequency characteristics of the transmitted underwater acoustic signal. In this scheme, the phase alignment information can be a phase alignment function, and the phase shift information can be a phase shift function. Both the phase alignment information and the phase shift information together characterize the phase offset required to align the received signal with the transmitted underwater acoustic signal. The phase alignment information can be represented by a phase alignment function. The phase shift information can be represented by a phase shift function.

[0031] In one embodiment, a frequency-modulated underwater acoustic signal with a known waveform can be transmitted through the transducer of a sonar, and the waveform information of the transmitted frequency-modulated underwater acoustic signal can be recorded. The spectral power information of the transmitted underwater acoustic signal can be extracted based on the waveform information, and the phase information of the transmitted underwater acoustic signal can be determined based on the correlation between the spectral power information and the phase information. Based on the phase information of the transmitted underwater acoustic signal, phase alignment information and phase shift information that need to be adjusted for the phase of the received underwater acoustic signal when the received underwater acoustic signal is aligned with the transmitted underwater acoustic signal can be determined.

[0032] In one embodiment, the spectral power function of the emitted underwater acoustic signal can be taken. The phase portion is obtained by taking the negative value of the acquired phase to obtain the phase alignment function. Substituting the alignment function into the preset phase alignment equation, the phase translation function is obtained. .

[0033] The expression for the preset phase alignment equation is: ; By solving equations Derive the transition function Transition function Substitution The phase shift function is obtained. Phase alignment function and phase translation function Together, they characterize the phase offset between the received and transmitted underwater acoustic signals.

[0034] Figure 2 This is a schematic diagram illustrating the effect of short-time Fourier transform processing of hyperbolic frequency modulated signals provided in the embodiments of this application. For example... Figure 2 As shown in the figure, the time-domain and frequency-domain simulation results of the received underwater acoustic signal under strong reverberation background are displayed. The conventional short-time Fourier transform in the figure has low time-frequency resolution, significant sidelobe interference, and the time-frequency trajectory of the weak target echo is almost completely submerged by the reverberant background, making it difficult to effectively detect and distinguish the target underwater acoustic signal in the received underwater acoustic signal. The specific parameters of the simulation environment for this scheme are: signal sampling frequency: 10000Hz; frequency band: 1000Hz~3000Hz; signal type: hyperbolic frequency modulation (HFM) signal; HFM signal pulse width: 2s; pulse period: 3s; start time: 1s; rectangular window length: 1024 points; number of iterations: 6; signal-to-mixing ratio: Assume there are three echo signals with start times of t=1s, 1.02s, and 1.2s, where the amplitude of the latter two is half that of the first. Simulations are used to process underwater acoustic signals received from multiple targets. Targets with relatively smaller amplitudes are considered weak signals. The background interference energy is set to the standard based on the energy of the first received underwater acoustic signal, with a signal-to-mixing ratio (SMR) of -10dB, to verify the performance of the method under low SMR conditions.

[0035] S103, based on the preset window center frequency and preset window splitting length, the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed to obtain multiple received window signals, and based on the phase alignment information and phase shift information, the received window signals are parametrically phase-synchronized with the frequency-time transformation to obtain a standard received window signal that is phase-synchronized with each transmitted window signal. The standard received window signal is a time-frequency signal.

[0036] The preset window center frequency can be the center frequency of each segmented window when splitting the transmitted and received signals along the frequency axis. The preset window split length can be a preset parameter representing the size of each segmented window. The received underwater acoustic signal can be a signal with the same identifier as the transmitted underwater acoustic signal caused by the transmitted underwater acoustic signal. Windowing processing can be the operation of splitting the spectrum function of the entire received underwater acoustic signal into multiple shorter frequency band signals along the frequency axis. There are signals with overlapping frequency bands between adjacent transmitted window signals among multiple transmitted window signals. There are signals with overlapping frequency bands between adjacent received window signals among multiple received window signals. The window lengths of the transmitted window signal and the received window signal are the same. The standard received window signal can be the signal obtained after phase correction of the received window signal. Parameterized phase synchronization frequency-time transformation can be the operation of aligning the received window signal with the transmitted window signal using phase alignment information and phase shift information.

[0037] In one embodiment, the transmitted underwater acoustic signal and the corresponding received underwater acoustic signal can be converted into frequency domain signals respectively. The frequency domain conversion results of the transmitted underwater acoustic signal and the corresponding received underwater acoustic signal are windowed according to a preset window frequency band length to obtain multiple receiving window signals. The phase parameters of each receiving window signal are phase-calibrated according to phase alignment information and phase shift information to align the phase of each receiving window signal with the phase of the corresponding transmitted window signal, thereby obtaining a standard receiving window signal that is phase-synchronized with each transmitted window signal.

[0038] The frequency domain conversion results of the received underwater acoustic signal can be corrected using phase alignment and phase shift functions. The parameterized phase-synchronized frequency-time transformation result is obtained. The formula for calculating the parameterized phase synchronization frequency-time conversion result is as follows: ; in, It receives the spectrum of underwater acoustic signals. It is a phase alignment function. It is a phase shift function. It is the frequency window obtained by splitting.

[0039] The above transformation achieves phase alignment and translation of the received underwater acoustic signal spectrum, thereby matching the received underwater acoustic signal with the transmitted underwater acoustic signal and obtaining a time-frequency distribution with high time-frequency concentration.

[0040] In one embodiment, the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed based on a preset window center frequency and a preset window splitting length to obtain multiple receiving window signals. This includes: performing Fourier transforms on the transmitted underwater acoustic signal and the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal to obtain transmitted spectrum information and received spectrum information; performing windowing calculations on the received spectrum information according to the preset window center frequency and the preset window splitting length to obtain multiple receiving window information, wherein there is an overlap between adjacent receiving windows corresponding to the multiple receiving window information; and determining multiple receiving window signals corresponding to the transmitted underwater acoustic signal based on the multiple receiving window information.

[0041] In one embodiment, the Fourier transform results of the transmitted and received underwater acoustic signals can be calculated separately to obtain transmitted and received spectrum information. The transmitted and received spectrum information are then split according to a preset window splitting length, with overlaps between adjacent transmitted window information and adjacent received window information. Multiple transmitted window information are used as multiple transmitted window signals corresponding to the transmitted underwater acoustic signals, and multiple received window information are used as multiple received window signals corresponding to the transmitted underwater acoustic signals.

[0042] The frequency spectrum function of the emitted underwater acoustic signal can be obtained by performing a Fourier transform on the emitted underwater acoustic signal. The spectral function is windowed along the frequency axis, and the window function formula is as follows: ; in, It is the result of the window function calculation. It is the center frequency of the frequency window. It is the window length of the rectangular window. It is the frequency axis of the frequency window.

[0043] There is a certain degree of overlap between adjacent transmit window information and adjacent receive window information to ensure the continuity of processing results. However, the overlap rate should not be too large to avoid increasing the number of windows that need to be calculated and significantly reducing the overall calculation speed.

[0044] This scheme obtains the transmitted and received spectrum information by performing Fourier transform on the received underwater acoustic signal, and performs windowing calculation on the received spectrum information according to the preset window center frequency and preset window splitting length to obtain multiple received window signals, which can improve the efficiency of window splitting of transmitted and received underwater acoustic signals.

[0045] S104, determine the transmission time-frequency distribution data corresponding to each transmission window signal and the reception time-frequency distribution data corresponding to each standard reception window signal, and perform kernel inversion iterative correction on each standard reception window signal based on the transmission time-frequency distribution data and the reception time-frequency distribution data until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than the preset error threshold.

[0046] The transmit time-frequency distribution data can be the frequency variation data of the transmit window signal over time. The receive time-frequency distribution data can be the frequency variation data of the receive window signal over time.

[0047] In one embodiment, the transmission time-frequency distribution data corresponding to each transmission window signal can be determined based on the signal frequency corresponding to each timestamp in each transmission window signal, and the reception time-frequency distribution data corresponding to each standard reception window signal can be determined based on the signal frequency corresponding to each timestamp in each reception window signal. Since the time-frequency distribution characteristics of the target underwater acoustic signal are consistent with those of the transmitted underwater acoustic signal, the transmission time-frequency distribution data can be used as the standard for the time-frequency distribution characteristics of the target underwater acoustic signal to perform kernel inversion iterative correction on each standard reception window signal, and the mean square error between each standard reception window signal and the corresponding reception time-frequency distribution data after each round of correction is calculated, until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than a preset error threshold. At this point, the reception time-frequency distribution data is considered to have consistent characteristics with the transmission time-frequency distribution data, and the iterative correction result of each standard reception window signal is obtained. The iterative correction result of each standard reception window signal is used for underwater acoustic reverberation suppression.

[0048] S105: The iterative correction results of the signals from each standard receiving window are spliced ​​together to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression.

[0049] In one embodiment, the iterative correction results of each standard receiving window signal can be spliced ​​together according to the window splitting order to obtain a complete target underwater acoustic signal corresponding to the transmitted underwater acoustic signal.

[0050] In one embodiment, the iterative correction results of each standard receiving window signal are spliced ​​together to obtain the target underwater acoustic signal after underwater reverberation suppression. This includes splicing the iterative correction results of each standard receiving window signal according to the window frequency band number corresponding to each standard receiving window signal to obtain the full-band target underwater acoustic signal after underwater reverberation suppression corresponding to the transmitted underwater acoustic signal. The full-band target underwater acoustic signal is a two-dimensional matrix including frequency and time dimensions.

[0051] The window frequency band number can be a number that indicates the order in which each window frequency band is split.

[0052] In one embodiment, the iterative correction results of each standard receiving window signal can be spliced ​​together according to the window frequency band number corresponding to each standard receiving window signal in the order of window splitting from front to back to obtain the full-band target underwater acoustic signal. At this time, the full-band target underwater acoustic signal is the target underwater acoustic signal after underwater acoustic reverberation suppression.

[0053] This scheme obtains the full-band target underwater acoustic signal by splicing the iterative correction results of each standard receiving window signal according to the window frequency band sequence number corresponding to each standard receiving window signal, which can further improve the accuracy of the target underwater acoustic signal extraction results.

[0054] Figure 3 This is a schematic diagram illustrating the effect of the hyperbolic frequency modulation signal kernel inversion parameterized phase synchronization frequency-time transformation processing provided in the embodiments of this application. For example... Figure 3 As shown, after adopting the method proposed in this application, the time-frequency convergence of the signal is significantly improved, the time-frequency trajectories of the three targets are clearly distinguishable, and the background in the figure is also cleaner. In the simulation, the main factor affecting the background is reverberation interference, which means that this scheme has excellent reverberation suppression performance.

[0055] Figure 4 This is a comparative schematic diagram of the hyperbolic frequency modulation signal processing results provided in the embodiments of this application. For example... Figure 4 As shown in the figure, the results of two signal processing methods, STFT and nuclear inversion parameterized time-frequency transformation, are presented. The signal processing result of this scheme has a narrower main lobe and a more concentrated target energy. That is, this scheme can effectively suppress reverberation interference and background noise interference and has good time-frequency resolution.

[0056] The technical solution provided in this application extracts the spectral power information of the transmitted underwater acoustic signal, determines the phase information of the transmitted underwater acoustic signal based on the spectral power information, and determines the phase alignment information and phase shift information corresponding to the transmitted underwater acoustic signal based on the phase information; performs windowing processing on the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal based on the preset window center frequency and preset window splitting length to obtain multiple received window signals, and performs parameterized phase synchronization frequency-time transformation on each received window signal based on the phase alignment information and phase shift information to obtain standard received window signals that are phase-synchronized with each transmitted window signal. The standard received window signals are time-frequency signals; determines the transmitted time-frequency distribution data corresponding to each transmitted window signal and the received time-frequency distribution data corresponding to each standard received window signal, and performs kernel inversion iterative correction on each standard received window signal based on the transmitted time-frequency distribution data and the received time-frequency distribution data until the mean square error between the correction result of each standard received window signal and the corresponding received time-frequency distribution data is less than a preset error threshold; and splices the iterative correction results of each standard received window signal to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression. The above-described method for suppressing underwater acoustic reverberation by stationary phase kernel inversion parameterized phase synchronization transformation solves the problems of low time-frequency resolution and inaccurate target underwater acoustic signal determination in existing technologies. Phase alignment and phase shift information are used to calibrate the phase of each receiving window signal, resulting in a standard receiving window signal that is phase-synchronized with each transmitting window signal. Kernel inversion iterative correction is performed on each standard receiving window signal based on the transmitting and receiving time-frequency distribution data. The iterative correction results of each standard receiving window signal are then combined to obtain the target underwater acoustic signal after reverberation suppression. This method achieves the goal of determining the target underwater acoustic signal in the received underwater acoustic signal based on the original characteristics of the transmitted underwater acoustic signal. It effectively overcomes the problems of low time-frequency resolution and large sidelobe interference in conventional time-frequency analysis, significantly suppresses reverberation and other interference signals in the received underwater acoustic signal, improves the time-frequency aggregation and detectability of weak target signals, and thus improves the accuracy of the target underwater acoustic signal determination.

[0057] Figure 5 This is a flowchart of the iterative correction standard receiving window signal provided in the embodiments of this application. For example... Figure 5 As shown, the specific steps include the following: S501 determines the inversion kernel information of the transmitted underwater acoustic signal corresponding to each transmission window signal based on the transmission time-frequency distribution data, and determines the observation information of the received underwater acoustic signal corresponding to each standard reception window signal based on the reception time-frequency distribution data.

[0058] The transmitted underwater acoustic signal inversion kernel information can represent the time-frequency characteristics of the transmitted underwater acoustic signal within each frequency window. This kernel information can be represented by the time-frequency distribution magnitude of the transmitted window signal. The received underwater acoustic signal observation information can represent the time-frequency distribution intensity characteristics of the received underwater acoustic signal itself. This observation information can be represented by the time-frequency distribution magnitude of the standard received window signal.

[0059] In one embodiment, the time-frequency distribution modulus of each transmission window signal can be calculated based on the transmission time-frequency distribution data, and the time-frequency distribution modulus can be used as the transmission underwater acoustic signal inversion kernel information corresponding to each transmission window signal. The time-frequency distribution modulus of each standard reception window signal can be calculated based on the reception time-frequency distribution data, and the time-frequency distribution modulus can be used as the reception underwater acoustic signal observation information corresponding to each standard reception window signal.

[0060] The time-frequency distribution modulus of each transmission window signal can be defined as the kernel information for transmitting underwater acoustic signal inversion, and the time-frequency distribution modulus of each standard receiving window signal can be defined as the observation information of the received underwater acoustic signal to be processed. The formula for calculating the kernel information for transmitting underwater acoustic signal inversion is as follows: ; The formula for calculating the received underwater acoustic signal observation information is as follows: ; in, It transmits underwater acoustic signals to retrieve nuclear information. It receives underwater acoustic signal observation information. Because the transmitted underwater acoustic signal is a known signal, parameterized phase synchronization is not required, and conventional time-frequency transformation can be performed directly. .

[0061] S502, based on the correspondence between each transmission window signal and each standard receiving window signal, the kernel information of the transmitted underwater acoustic signal inversion and the observation information of the received underwater acoustic signal, the kernel inversion iterative correction is performed on each standard receiving window signal to obtain the iterative correction result corresponding to each standard receiving window signal.

[0062] In one embodiment, the kernel information of the transmitted underwater acoustic signal inversion can be used as the time-frequency distribution reference standard information of the target underwater acoustic signal. The time-frequency distribution reference standard information corresponding to each standard receiving window signal is determined according to the correspondence between each transmitted window signal and each standard receiving window signal. The kernel inversion iterative correction is performed on each standard receiving window signal according to the time-frequency distribution reference standard information until the time-frequency distribution characteristics of the received underwater acoustic signal observation information match the kernel information of the transmitted underwater acoustic signal inversion, and the iterative correction result corresponding to each standard receiving window signal is obtained.

[0063] In one embodiment, based on the correspondence between each transmission window signal and each standard receiving window signal, the transmission underwater acoustic signal inversion kernel information, and the received underwater acoustic signal observation information, the kernel inversion iterative correction is performed on each standard receiving window signal. This includes: determining the target transmission underwater acoustic signal inversion kernel information corresponding to the received underwater acoustic signal observation information of each standard receiving window signal in multiple transmission underwater acoustic signal inversion kernel information based on the correspondence between each transmission window signal and each standard receiving window signal; and performing kernel inversion iterative correction on each standard receiving window signal based on a preset iterative algorithm, the received underwater acoustic signal observation information, and the target transmission underwater acoustic signal inversion kernel information.

[0064] In one embodiment, based on the correspondence between each transmit window signal and each standard receive window signal, the inversion kernel information of each transmit underwater acoustic signal can be matched with the received underwater acoustic signal observation information of each standard receive window signal to obtain the target transmit underwater acoustic signal inversion kernel information corresponding to the received underwater acoustic signal observation information of each standard receive window signal. According to a preset iterative algorithm, using the target transmit underwater acoustic signal inversion kernel information as a standard, kernel inversion iterative correction is performed on each standard receive window signal until the received underwater acoustic signal observation information and the target transmit underwater acoustic signal inversion kernel information are consistent in characteristics.

[0065] This scheme determines the target transmitted underwater acoustic signal inversion kernel information corresponding to the received underwater acoustic signal observation information of each standard receiving window signal in multiple transmitted underwater acoustic signal inversion kernel information. Based on the preset iterative algorithm, the received underwater acoustic signal observation information and the target transmitted underwater acoustic signal inversion kernel information, the kernel inversion iterative correction is performed on each standard receiving window signal, which can improve the accuracy of the correction of the standard receiving window signal.

[0066] In one embodiment, the formula for the preset iterative algorithm is: ; in, It is the first The first of the standard received window signals The result of the second iteration It is the first The i-th iteration result of the standard receive window signal. It is the first The received underwater acoustic signal observation information of the standard receiving window signal. It is the target transmitted underwater acoustic signal inversion kernel information for each of the aforementioned standard receiving window signals. yes The time-frequency distribution pattern, It is a frequency dimension variable. It is a time-dimensional variable.

[0067] In one embodiment, nuclear information can be retrieved by transmitting underwater acoustic signals to the target. The Lucy-Richardson iterative algorithm is used, with the desired signal obtained from the underwater acoustic signal observation information. The iteration stops when the mean square error between the iteration result and the received underwater acoustic signal observation information is less than a preset mean square error threshold. The iteration result is the kernel inversion result within a single frequency window. .

[0068] The formula for determining when iteration stops is: ; in, This is a preset mean squared error threshold. The preset mean squared error threshold is used to control the convergence of the iteration; when the estimated result... The root mean square value of the inversion kernel is less than When the iteration converges, it is determined that the iteration has converged.

[0069] Nuclear inversion results with different frequency windows The processing results represent different frequency bands, and the kernel inversion results represent different windows. By piecing them together along the frequency axis, we obtain the overall nuclear inversion result. Overall nuclear inversion results The expression is: ; In the formula, It is the window number. It represents the total number of frequency windows.

[0070] The above expression means that the processing results of each frequency window are combined. In this way, when strong interference occurs in a certain frequency window, the frequency window processing can prevent this interference from affecting the processing results of other frequency windows.

[0071] Nuclear inversion results The high-resolution nuclear inversion parameterized phase-synchronized frequency-time transform with the ability to suppress reverberation and background interference is a two-dimensional matrix, where the dimensions of the matrix represent the frequency and time of the received signal, respectively.

[0072] This scheme simplifies the iterative correction steps and improves the efficiency of iterative correction by using a pre-defined iterative algorithm formula to perform kernel inversion iterative correction on the signals of each standard receiving window.

[0073] The technical solution provided in this application determines the inversion kernel information of the transmitted underwater acoustic signal corresponding to each transmission window signal by using transmission time-frequency distribution data, and determines the received underwater acoustic signal observation information corresponding to each standard receiving window signal based on the received time-frequency distribution data. According to the correspondence between each transmission window signal and each standard receiving window signal, the kernel inversion iterative correction is performed on each standard receiving window signal to obtain the iterative correction result corresponding to each standard receiving window signal. This can improve the accuracy of the iterative correction of the standard receiving window signal and is beneficial to the accuracy of target signal determination.

[0074] Figure 6 This is a flowchart of another underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method provided in this application embodiment. For example... Figure 6 As shown, the specific steps include the following: S601, the sonar transducer transmits a frequency-modulated transmission signal of known waveform. It records the waveform information of the transmitted FM signal and acquires the received time-domain signal, including the target echo and reverberation background, through a sonar receiver. For the received time domain signal The received frequency domain signal is obtained by performing a Fourier transform. .

[0075] S602 uses the stationary phase method to approximate the solution of the phase-transmitted frequency modulation signal. The spectral power function of the transmitted signal is obtained. .

[0076] S603, perform a Fourier transform on the transmitted signal to obtain the spectrum function. The spectral function is windowed along the frequency axis, and the window function is... There needs to be some overlap between adjacent frequency window functions to ensure the continuity of the processing results.

[0077] S604, take The phase portion is obtained by taking the negative value of the acquired phase to obtain the phase alignment function. Substituting the alignment function into the phase alignment equation, we obtain the phase translation function. .

[0078] S605 uses a phase alignment function and a phase shift function to correct the Fourier transform result of the received signal. The parameterized phase-synchronized frequency-time transformation result is obtained. .

[0079] S606, the time-frequency distribution modulus of the transmitted signal is defined as the inversion kernel required for the kernel inversion process. The time-frequency distribution magnitude of the received signal is defined as the observation data to be processed. .

[0080] S607, for inversion kernels The Lucy-Richardson iterative algorithm is used, with the desired signal being the observed data. The iteration stops when the mean square error between the estimated result and the observed data is less than a set threshold. The iteration result is the kernel inversion result within a single frequency window. .

[0081] S608, Kernel inversion results with different windows The processing results represent different frequency bands, and the kernel inversion results represent different windows. By splicing the results together, the overall nuclear inversion result is obtained. .

[0082] S609, Nuclear Inversion Results The result of high-resolution time-frequency analysis with suppressed reverberation and background interference is a two-dimensional matrix, where the dimensions of the matrix represent frequency and time, respectively.

[0083] The technical solution provided in this application first extracts the phase features of the transmitted FM signal using the stationary phase method, constructing a phase alignment function and a phase shift function. Then, it performs phase correction and time-frequency focusing on the received signal through parameterized phase synchronization time-time transformation, suppressing reverberation and other interference to form a highly concentrated time-frequency distribution. Next, using the ideal time-frequency modulus of the transmitted signal as the inversion kernel and the time-frequency modulus of the received signal after synchronization processing as the observation data, it uses the Lucy-Richardson iterative algorithm to suppress and invert the reverberant background, ultimately obtaining a high-resolution time-frequency analysis result with a clean background and prominent targets. This method overcomes the time-frequency resolution limitations of fixed window functions, significantly reducing the sidelobes and background levels of the time-frequency distribution, narrowing the main lobe, and thus enhancing the concentration of main lobe energy. This allows for clear differentiation of targets even under low signal-to-mixing conditions, significantly improving the detection capability and feature extraction accuracy of underwater targets in strong reverberant backgrounds.

[0084] Figure 7 This is a structural block diagram of a hydroacoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation device provided in an embodiment of this application. Figure 7 As shown, it specifically includes the following: The power spectrum information extraction module 701 is used to determine the spectral power information of the transmitted underwater acoustic signal based on a preset stationary phase algorithm; The phase alignment information determination module 702 is used to determine the phase function of the transmitted underwater acoustic signal based on the spectral power information, and to invert the phase function to obtain the phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters. The phase shift information determination module 703 is used to calculate the phase shift information corresponding to the emitted underwater acoustic signal based on the preset phase alignment equation and phase alignment information.

[0085] The window signal alignment module 704 is used to perform window processing on the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal based on the preset window center frequency and the preset window splitting length to obtain multiple received window signals. Based on the phase alignment information and phase shift information, the module performs parameterized phase synchronization frequency-time transformation on each received window signal to obtain a standard received window signal that is phase synchronized with each transmitted window signal. The standard received window signal is a time-frequency signal. The underwater acoustic signal correction module 705 is used to determine the transmission time-frequency distribution data corresponding to each transmission window signal and the reception time-frequency distribution data corresponding to each standard reception window signal, and to perform kernel inversion iterative correction on each standard reception window signal based on the transmission time-frequency distribution data and the reception time-frequency distribution data until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than a preset error threshold. The target underwater acoustic signal determination module 706 is used to splice the iterative correction results of the signals from each standard receiving window to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression.

[0086] Furthermore, the underwater acoustic signal correction module 705 is specifically used for: Based on the transmission time-frequency distribution data, the kernel information of the transmitted underwater acoustic signal inversion corresponding to each transmission window signal is determined, and the observation information of the received underwater acoustic signal corresponding to each standard reception window signal is determined based on the reception time-frequency distribution data. Based on the correspondence between each transmission window signal and each standard receiving window signal, the kernel information of the transmitted underwater acoustic signal inversion and the observation information of the received underwater acoustic signal, the kernel inversion iterative correction is performed on each standard receiving window signal to obtain the iterative correction result corresponding to each standard receiving window signal.

[0087] Furthermore, the underwater acoustic signal correction module 705 is specifically used for: Based on the correspondence between each transmission window signal and each standard receiving window signal, the target transmission underwater acoustic signal inversion kernel information corresponding to the received underwater acoustic signal observation information of each standard receiving window signal is determined in multiple transmission underwater acoustic signal inversion kernel information; Based on the preset iterative algorithm, the received underwater acoustic signal observation information, and the target emitted underwater acoustic signal inversion kernel information, the kernel inversion iterative correction is performed on the signals of each standard receiving window.

[0088] Furthermore, the formula for the preset iterative algorithm is: ; in, It is the first The first of the standard received window signals The result of the second iteration It is the first The i-th iteration result of the standard receive window signal. It is the first The received underwater acoustic signal observation information of the standard receiving window signal. It is the target transmitted underwater acoustic signal inversion kernel information for each of the aforementioned standard receiving window signals. yes The time-frequency distribution pattern, It is a frequency dimension variable. It is a time-dimensional variable.

[0089] Furthermore, the windowed signal alignment module 704 is specifically used for: Fourier transforms are performed on the transmitted underwater acoustic signal and the corresponding received underwater acoustic signal to obtain the transmitted spectrum information and the received spectrum information, respectively. The received spectrum information is windowed according to the preset window center frequency and preset window split length to obtain multiple receiving window information. There is an overlap between adjacent receiving windows corresponding to multiple receiving window information. Multiple receiving window signals corresponding to the transmitted underwater acoustic signal are determined based on information from multiple receiving windows.

[0090] Furthermore, the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation module 706 is specifically used for: The iterative correction results of each standard receiving window signal are spliced ​​together according to the window frequency band sequence number corresponding to each standard receiving window signal to obtain the full-band target underwater acoustic signal after underwater acoustic reverberation suppression corresponding to the transmitted underwater acoustic signal. The full-band target underwater acoustic signal is a two-dimensional matrix including frequency and time dimensions.

[0091] The technical solution provided in this application includes a power spectrum information extraction module for determining the spectral power information of the transmitted underwater acoustic signal based on a preset stationary phase algorithm; a phase alignment information determination module for determining the phase function of the transmitted underwater acoustic signal based on the spectral power information, and inverting the phase function to obtain phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters; a phase shift information determination module for calculating the phase shift information corresponding to the transmitted underwater acoustic signal based on a preset phase alignment equation and phase alignment information; and a windowed signal alignment module for performing windowed processing on the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal based on a preset window center frequency and a preset window splitting length, obtaining multiple received window signals, and aligning them based on the phase alignment information. The system performs parameterized phase-synchronized frequency-time transformation on each receiving window signal using information and phase shift data to obtain a standard receiving window signal that is phase-synchronized with each transmitting window signal. The standard receiving window signal is a time-frequency signal. The receiving underwater acoustic signal correction module is used to determine the transmitting time-frequency distribution data corresponding to each transmitting window signal and the receiving time-frequency distribution data corresponding to each standard receiving window signal. Based on the transmitting time-frequency distribution data and the receiving time-frequency distribution data, the system performs kernel inversion iterative correction on each standard receiving window signal until the mean square error between the correction result of each standard receiving window signal and the corresponding receiving time-frequency distribution data is less than a preset error threshold. The target underwater acoustic signal determination module is used to stitch together the iterative correction results of each standard receiving window signal to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression. The aforementioned underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation device solves the problems of low time-frequency resolution and inaccurate target underwater acoustic signal determination in existing technologies. It performs phase calibration on each receiving window signal using phase alignment and phase shift information to obtain a standard receiving window signal that is phase-synchronized with each transmitting window signal. Based on the transmitting and receiving time-frequency distribution data, it performs kernel inversion iterative correction on each standard receiving window signal and splices the iterative correction results of each standard receiving window signal to obtain the target underwater acoustic signal after reverberation suppression. This achieves the goal of determining the target underwater acoustic signal in the received underwater acoustic signal based on the original characteristics of the transmitted underwater acoustic signal, effectively overcoming the problems of low time-frequency resolution and large sidelobe interference in conventional time-frequency analysis. It significantly suppresses reverberation and other interference signals in the received underwater acoustic signal, improves the time-frequency convergence and detectability of weak target signals, and thus improves the accuracy of the target underwater acoustic signal determination result.

[0092] The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation system in this application embodiment can be configured in a device, or in a component, integrated circuit, or chip in a terminal. This system can be configured in mobile electronic devices or non-mobile electronic devices. For example, mobile electronic devices can be mobile phones, tablets, laptops, PDAs, in-vehicle electronic devices, wearable devices, ultra-mobile personal computers (UMPCs), netbooks, or personal digital assistants (PDAs), etc., while non-mobile electronic devices can be servers, network-attached storage (NAS), personal computers (PCs), televisions (TVs), ATMs, or self-service machines, etc. This application embodiment does not impose specific limitations.

[0093] The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation system in this application embodiment can be an operating system. This operating system can be Android, iOS, or other possible operating systems; this application embodiment does not specifically limit it.

[0094] The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation system provided in this application embodiment can realize the various processes implemented in the above method embodiments. To avoid repetition, it will not be described again here.

[0095] like Figure 8 As shown, this application embodiment also provides an electronic device 800, including a processor 801, a memory 802, and a program or instructions stored in the memory 802 and executable on the processor 801. When the program or instructions are executed by the processor 801, they implement the various processes of the above-described embodiment of the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method, and can achieve the same technical effect. To avoid repetition, they will not be described again here.

[0096] It should be noted that the electronic devices in the embodiments of this application include the mobile electronic devices and non-mobile electronic devices described above.

[0097] This application also provides a readable storage medium storing a program or instructions. When the program or instructions are executed by a processor, they implement the various processes of the above-described embodiment of the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method, and achieve the same technical effect. To avoid repetition, they will not be described again here.

[0098] The processor is the processor in the electronic device described in the above embodiments. The readable storage medium includes computer-readable storage media, such as computer read-only memory (ROM), random access memory (RAM), magnetic disk, or optical disk.

[0099] This application also provides a program product including program code. When the program product is run on a computer device, the program code causes the computer device to perform the steps of the methods described above according to various exemplary embodiments of this application. For example, the computer device can execute a method for underwater acoustic reverberation suppression, stationary phase kernel inversion, and parameterized phase synchronization transformation as described in an embodiment of this application. The program product can be implemented using any combination of one or more readable media.

[0100] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or system. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or system that includes that element. Furthermore, it should be noted that the scope of the methods and systems in the embodiments of this application is not limited to performing functions in the order shown or discussed, but may also include performing functions substantially simultaneously or in the reverse order, depending on the functions involved. For example, the described methods may be performed in a different order than described, and various steps may be added, omitted, or combined. Additionally, features described with reference to certain examples may be combined in other examples.

[0101] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a computer software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) and includes several instructions to cause a terminal (which may be a mobile phone, computer, server, or network device, etc.) to execute the methods described in the various embodiments of this application.

[0102] The embodiments of this application have been described above with reference to the accompanying drawings. However, this application is not limited to the specific embodiments described above. The specific embodiments described above are merely illustrative and not restrictive. Those skilled in the art can make many other forms under the guidance of this application without departing from the spirit and scope of the claims, and all of these forms are within the protection scope of this application.

[0103] The above description is merely a preferred embodiment and the technical principles employed in this application. This application is not limited to the specific embodiments described herein, and various obvious changes, readjustments, and substitutions that can be made by those skilled in the art will not depart from the scope of protection of this application. Therefore, although this application has been described in detail through the above embodiments, this application is not limited to the above embodiments, and may include more other equivalent embodiments without departing from the concept of this application, the scope of which is determined by the scope of the claims.

Claims

1. A method for underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation, characterized in that, The method includes: The spectral power information of the transmitted underwater acoustic signal is determined based on a preset stationary phase algorithm. The phase function of the transmitted underwater acoustic signal is determined based on the spectral power information. The phase function is then inverted to obtain phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters. Calculate the phase shift information corresponding to the emitted underwater acoustic signal based on the preset phase alignment equation and the phase alignment information; The received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed based on the preset window center frequency and the preset window splitting length to obtain multiple received window signals. Then, based on the phase alignment information and the phase shift information, parameterized phase synchronization frequency-time transformation is performed on each of the received window signals to obtain a standard received window signal that is phase-synchronized with each of the transmitted window signals. The standard received window signal is a time-frequency signal. The transmission time-frequency distribution data corresponding to each of the transmission window signals and the reception time-frequency distribution data corresponding to each of the standard reception window signals are determined respectively. Based on the transmission time-frequency distribution data and the reception time-frequency distribution data, the standard reception window signals are subjected to kernel inversion iterative correction until the mean square error between the correction result of each standard reception window signal and the corresponding reception time-frequency distribution data is less than a preset error threshold. The target underwater acoustic signal after underwater reverberation suppression is obtained by splicing the iterative correction results of each standard receiving window signal.

2. The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method according to claim 1, characterized in that, The step of performing kernel inversion iterative correction on each of the standard receiving window signals based on the transmitted time-frequency distribution data and the received time-frequency distribution data includes: Based on the transmission time-frequency distribution data, determine the transmission underwater acoustic signal inversion kernel information corresponding to each of the transmission window signals, and based on the reception time-frequency distribution data, determine the reception underwater acoustic signal observation information corresponding to each of the standard reception window signals; Based on the correspondence between each transmitted window signal and each standard received window signal, the kernel information of the transmitted underwater acoustic signal inversion, and the observation information of the received underwater acoustic signal, the kernel inversion iterative correction is performed on each standard received window signal to obtain the iterative correction result corresponding to each standard received window signal.

3. The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method according to claim 2, characterized in that, The step of performing kernel inversion iterative correction on each of the standard receiving window signals based on the correspondence between each of the transmitted window signals and each of the standard receiving window signals, the kernel information of the transmitted underwater acoustic signal inversion, and the observation information of the received underwater acoustic signal includes: Based on the correspondence between each of the transmitted window signals and each of the standard received window signals, determine the target transmitted underwater acoustic signal inversion kernel information that corresponds to the received underwater acoustic signal observation information of each of the standard received window signals in the multiple transmitted underwater acoustic signal inversion kernel information; Based on the preset iterative algorithm, the received underwater acoustic signal observation information, and the target emitted underwater acoustic signal inversion kernel information, the kernel inversion iterative correction is performed on each of the standard receiving window signals.

4. The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method according to claim 3, characterized in that, The formula for the preset iterative algorithm is: ； in, It is the first The first of the standard received window signals The result of the second iteration It is the first The i-th iteration result of the standard receive window signal. It is the first The received underwater acoustic signal observation information of the standard receiving window signal. It is the target transmitted underwater acoustic signal inversion kernel information for each of the aforementioned standard receiving window signals. yes The time-frequency distribution pattern, It is a frequency dimension variable. It is a time-dimensional variable.

5. The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method according to claim 1, characterized in that, The received underwater acoustic signal corresponding to the transmitted underwater acoustic signal is windowed based on a preset window center frequency and a preset window splitting length to obtain multiple received window signals, including: Fourier transforms are performed on the transmitted underwater acoustic signal and the corresponding received underwater acoustic signal to obtain the transmitted spectrum information and the received spectrum information, respectively. The received spectrum information is windowed according to the preset window center frequency and preset window split length to obtain multiple receiving window information, and there is an overlap between adjacent receiving windows corresponding to the multiple receiving window information. Based on the information of the multiple receiving windows, multiple receiving window signals corresponding to the transmitted underwater acoustic signal are determined.

6. The underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method according to claim 1, characterized in that, The iterative correction results of splicing the standard receiving window signals are used to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression, including: The iterative correction results of each standard receiving window signal are spliced ​​together according to the window frequency band sequence number corresponding to each standard receiving window signal to obtain the full-band target underwater acoustic signal after underwater acoustic reverberation suppression corresponding to the transmitted underwater acoustic signal. The full-band target underwater acoustic signal is a two-dimensional matrix including frequency and time dimensions.

7. A device for underwater acoustic reverberation suppression, stationary phase kernel inversion, parameterized phase synchronization transformation, characterized in that, The device includes: The power spectrum information extraction module is used to determine the spectral power information of the transmitted underwater acoustic signal based on a preset stationary phase algorithm; The phase alignment information determination module is used to determine the phase function of the transmitted underwater acoustic signal based on the spectral power information, and to invert the phase function to obtain phase alignment information corresponding to the transmitted underwater acoustic signal, which is used to provide basic phase parameters; The phase shift information determination module is used to calculate the phase shift information corresponding to the emitted underwater acoustic signal based on the preset phase alignment equation and the phase alignment information. The windowed signal alignment module is used to perform windowed processing on the received underwater acoustic signal corresponding to the transmitted underwater acoustic signal based on a preset window center frequency and a preset window splitting length to obtain multiple received window signals. Based on the phase alignment information and the phase shift information, the module performs parameterized phase synchronization frequency-time transformation on each of the received window signals to obtain a standard received window signal that is phase-synchronized with each of the transmitted window signals. The standard received window signal is a time-frequency signal. The underwater acoustic signal correction module is used to determine the transmission time-frequency distribution data corresponding to each of the transmission window signals and the reception time-frequency distribution data corresponding to each of the standard reception window signals, and to perform kernel inversion iterative correction on each of the standard reception window signals based on the transmission time-frequency distribution data and the reception time-frequency distribution data until the mean square error between the correction result of each of the standard reception window signals and the corresponding reception time-frequency distribution data is less than a preset error threshold. The target underwater acoustic signal determination module is used to splice the iterative correction results of each of the standard receiving window signals to obtain the target underwater acoustic signal after underwater acoustic reverberation suppression.

8. An electronic device, characterized in that, It includes a processor, a memory, and a program or instructions stored in the memory and running on the processor, wherein when the program or instructions are executed by the processor, they implement the steps of the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method as described in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program or instructions that, when executed by a processor, implement the steps of the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method as described in any one of claims 1-6.

10. A computer program product, characterized in that, The computer program product includes a computer program stored in a computer-readable storage medium. When at least one processor of the electronic device reads and executes the computer program from the computer-readable storage medium, it implements the steps of the underwater acoustic reverberation suppression stationary phase kernel inversion parameterized phase synchronization transformation method as described in any one of claims 1-6.

Images