A method for predicting ship-wide acoustic compatibility based on time-frequency space energy

Abstract

The present invention discloses a method for predicting the acoustic compatibility of a whole ship based on time-frequency space energy, belonging to the field of acoustic compatibility technology. The method comprises steps S1, establishing an acoustic compatibility technical parameter library based on multidimensional characteristics; S2, calculating the interference level at the interfered acoustic device; S3, calculating the interference intensity entering the front-end receiving channel of the interfered acoustic device, and judging whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device based on the interference intensity; S4, judging whether the interfered acoustic device generates in-band interference based on the multidimensional characteristics of acoustic compatibility; and S5, based on the calculation and judgment results of S3 and S4, quantitatively evaluating the comprehensive impact of the interference, and generating a prediction result of the acoustic compatibility of the whole ship based on time-frequency space energy. The method more accurately predicts the acoustic compatibility status between shipboard acoustic devices, supporting more reasonable acoustic compatibility design of the ship as a whole. The method effectively supports the acoustic compatibility design of the whole and equipment in the early design stages, thereby optimizing the acoustic compatibility of the whole ship.

Description

Technical Field

[0001] The present invention belongs to the technical field of acoustic compatibility, and in particular relates to a method for predicting the acoustic compatibility of a whole ship based on time-frequency space energy. Background Art

[0002] With the continuous development and application of acoustic equipment throughout the ship, acoustic compatibility has become increasingly prominent, becoming a key factor affecting ship performance and operational safety. Traditional acoustic compatibility analysis methods primarily rely on single-dimensional parameters, such as frequency, time, or space, making it difficult to comprehensively and accurately assess and predict acoustic interference between multiple devices. For example, some methods focus solely on the frequency domain characteristics of the interference source, ignoring the impact of time and space domain characteristics on acoustic compatibility, resulting in insufficient prediction accuracy in complex environments. Furthermore, previous ship-wide acoustic compatibility design primarily relied on interference magnitude, employing far-field theoretical calculations and simulations to estimate inter-device acoustic compatibility. However, with the increasing number of underwater functions within a ship, the number of onboard acoustic devices has increased, operating frequency bands have expanded, operating modes have become more diverse, and frequency band overlap has become a prominent issue. Traditional analysis methods based on interference magnitude are no longer able to address complex acoustic compatibility scenarios, making it difficult to comprehensively and accurately predict the acoustic compatibility status between onboard acoustic devices. This makes it difficult to propose correct and reasonable anti-interference measures and technical requirements during the overall ship and equipment design phases, and thus, cannot support ship-wide acoustic compatibility design.

[0003] Therefore, it is urgent to propose a full-ship acoustic compatibility prediction method based on time-frequency space energy to solve the above problems. Summary of the Invention

[0004] In view of this, the purpose of the present invention is to provide a method for predicting the acoustic compatibility of the entire ship based on time-frequency space energy, which can be used to predict the acoustic interference between ship-borne acoustic equipment, support the proposal of quantitative acoustic compatibility technical requirements and preliminary anti-interference measures, and support the acoustic compatibility design of the ship as a whole and acoustic equipment, which is conducive to improving the acoustic compatibility of the entire ship, solving the problem of no basis for the acoustic compatibility design of the ship as a whole, system and acoustic equipment, and improving underwater performance.

[0005] In order to achieve the above object, the present invention provides the following technical solutions:

[0006] A method for predicting the acoustic compatibility of a whole ship based on time-frequency space energy comprises the following steps:

[0007] S1. Establish an acoustic compatibility technical parameter library based on multi-dimensional characteristics, including the acoustic compatibility technical parameters of the interference source and the acoustic compatibility technical parameters of the interfered acoustic equipment;

[0008] S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device based on the energy domain characteristics, spatial domain characteristics, and interference propagation attenuation of the interference source acoustic device;

[0009] S3. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device based on the interference level at the interfered acoustic device and relevant acoustic compatibility technical parameters, and determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device based on the interference intensity;

[0010] S4. Determine whether the interfered acoustic device generates in-band interference based on the multi-dimensional characteristics of acoustic compatibility;

[0011] S5. Based on the calculation and judgment results of S3 and S4, the comprehensive impact of the interference is quantitatively evaluated to generate the full-ship acoustic compatibility prediction results of time-frequency space energy.

[0012] As a further preference of the present invention, the acoustic compatibility multi-dimensional technical parameters of the interference source described in S1 include the energy domain of the interference source, the time domain of the interference source, the frequency domain of the interference source and the spatial domain of the interference source; the acoustic compatibility technical parameters of the interfered acoustic device include the acoustic compatibility technical parameters of the receiving sensor, the acoustic compatibility technical parameters of the front-end receiving channel and the acoustic compatibility technical parameters of the digital signal processing; wherein the acoustic compatibility technical parameters of the receiving sensor include the spatial domain and frequency domain of the receiving sensor; the acoustic compatibility technical parameters of the front-end receiving channel include the frequency domain and energy domain of the front-end receiving channel; the acoustic compatibility technical parameters of the digital signal processing include the frequency domain, spatial domain and time domain of the digital signal processing.

[0013] As a further preferred embodiment of the present invention, S2 specifically includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, spatial domain characteristics and interference propagation attenuation of the interference source acoustic device, and the formula is as follows:

[0014]

[0015] in, is the interference level at the interfered acoustic device, SL is the emission sound source level of the interference source; DIi is the emission directivity of the interference source. When the interfered acoustic device is in the main lobe azimuth of the interference source, DIi=0. When it is in the side lobe azimuth, it is the main lobe to side lobe ratio; TL is the interference propagation attenuation, including distance propagation loss and hull structure shielding. The formula is as follows:

[0016]

[0017] Among them, g represents the acceleration of gravity, R represents the distance between the two devices, represents the seawater absorption coefficient, Indicates the attenuation of interference caused by hull structure obstruction.

[0018] As a further preferred embodiment of the present invention,

[0019] Said S3 specifically includes:

[0020] S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device:

[0021] The intensity of interference picked up by the receiving transducer of the interfered acoustic device The formula is as follows:

[0022]

[0023] Where P represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; N represents the attenuation of the interference source direction compared to the detection direction of the interfered acoustic device;

[0024] S32. Determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device:

[0025] According to the frequency response characteristics of amplification and filtering in the front-end receiving channel of the interfered acoustic device, the output of each level of the front-end receiving channel of the interfered acoustic device is calculated and compared with the corresponding limit voltage. If the limit voltage is exceeded, it indicates saturation; if not, it indicates non-saturation;

[0026] Use the following formula to determine whether the frequency response characteristics of the preamplifier are saturated:

[0027]

[0028] in, is the frequency response characteristic of the preamplifier, Indicates the limiting voltage of the preamplifier; when ≥0, it means saturation occurs. When <0, no saturation occurs;

[0029] Use the following formula to determine whether the pre-filter is saturated:

[0030]

[0031] in, is the frequency response characteristic of the pre-filter, Indicates the limiting voltage of the pre-filter; when ≥0, it means saturation occurs. No saturation occurs.

[0032] As a further preferred embodiment of the present invention, S4 specifically includes:

[0033] Based on the multi-dimensional characteristics of the interference source and the interfered acoustic device in the time domain, frequency domain, and spatial domain, the following formula is used to determine whether the interfered acoustic device generates in-band interference:

[0034]

[0035] Where: M represents the background noise level when the interfered acoustic equipment is demonstrated, Y represents the emission bandwidth of the interference source, g represents the acceleration of gravity, Indicates the attenuation of the interference source's transmission band compared to the interference frequency. Indicates the frequency domain interference suppression capability of the interfered acoustic equipment, Indicates the spatial interference suppression capability of the interfered acoustic equipment. Indicates the time domain interference suppression capability of the interfered acoustic equipment;

[0036] Calculate the frequency domain interference suppression capability based on the frequency domain processing method of the interfered acoustic device , when using line spectrum processing, =0; when broadband processing is used, The calculation method is as follows:

[0037]

[0038] in, represents the processing bandwidth of the interfered acoustic device, Indicates the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic equipment;

[0039] Calculation of spatial interference suppression based on the receiving directivity of the interfered acoustic device , when the interference source is in the detection main lobe direction of the interfered acoustic equipment, =0; when in the side lobe, The main-side lobe ratio of this azimuth is greater than the main lobe;

[0040] Calculates time-domain interference suppression based on the overlap between the interference source's transmission and the victim's acoustic device's processing period, as well as the victim's acoustic device's processing integration time. , The calculation method is shown in the following formula:

[0041]

[0042] in, Indicates the overlap duration, Indicates the integration duration of the disturbed acoustic device.

[0043] As a further preferred embodiment of the present invention, the S5 specifically includes: S51, constructing an interference comprehensive impact index: as shown in the following formula:

[0044]

[0045] in, and is the weight coefficient, =0.6, =0.4, is the calculation result of S3, and is the corresponding limiting voltage threshold, K is the calculation result of S4, and NL is the background noise level; and is the nonlinear correction factor, =1.2, =1.5.

[0046] As a further preference of the present invention, the S5 also includes S52, dividing the interference level based on the CII value. If the CII value is less than 0.3, it is a slight interference, and only the interference event needs to be recorded without active suppression; if 0.3≤CII<0.7, it is moderate interference, and frequency domain notching or spatial domain nulling is enabled; if CII≥0.7, it is severe interference, and the device working frequency band is forced to switch or time domain asynchronous sampling is started.

[0047] As a further preferred embodiment of the present invention, the S5 further includes S53, confidence assessment, which introduces a Monte Carlo simulation to verify the reliability of the result, as shown below:

[0048]

[0049] in, represents the mean value of CII, which is obtained by randomly perturbing the input parameters multiple times. is the standard deviation of CII;

[0050] If the confidence level is less than 0.8, the process returns to S51 to recalibrate the parameters.

[0051] The beneficial effects of the present invention are:

[0052] This paper proposes a ship-wide acoustic compatibility prediction method based on time-frequency-space energy. This method includes determining whether saturation and in-band interference occur in the front-end receiving channels of interfered acoustic devices. This method considers the multi-dimensional characteristics of the interference source and the transmitted and received signals of the interfered acoustic devices, including time, frequency, and space, as well as factors such as interference propagation attenuation. To address the current complex acoustic compatibility landscape, this method analyzes acoustic interference between devices from multiple dimensions, including time, frequency, and space. This method more accurately predicts the acoustic compatibility status between shipboard acoustic devices, supports more rational acoustic device layout within the ship, investigates anti-interference isolation measures, and proposes quantitative anti-interference technical requirements for shipboard acoustic devices. Furthermore, this method supports the demonstration and design of shipboard acoustic devices, more rationally setting operating frequency bands, designing acoustic compatibility technical indicators for transmitters and receivers, and reserving management and control ports. It also more accurately predicts the acoustic compatibility status between shipboard acoustic devices, quantifies the interference between them, and supports the subsequent development of targeted preliminary acoustic compatibility anti-interference measures for saturation and in-band interference. This method effectively supports acoustic compatibility design of both the overall system and the equipment at the initial design stage, optimizing ship-wide acoustic compatibility. Achieve "early detection, early adjustment, and early response" to acoustic compatibility issues, avoid serious problems of "innate deficiencies" in acoustic compatibility between shipborne acoustic equipment, such as co-frequency, overlapping frequency, and small frequency intervals, optimize the acoustic compatibility of the entire ship, and ensure full utilization of the ship's underwater performance.

[0053] Other advantages, objectives and features of the present invention will be described in the following description and will be apparent to those skilled in the art to some extent, or those skilled in the art can be taught from the practice of the present invention. The objectives and other advantages of the present invention can be realized and obtained through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to make the purpose, technical solutions and beneficial effects of the present invention more clear, the present invention provides the following drawings for illustration:

[0055] Figure 1 This is an overall flow chart of a method for predicting full-ship acoustic compatibility based on time-frequency space energy according to the present invention;

[0056] Figure 2 Schematic diagram of the acoustic compatibility technical parameter library of the present invention;

[0057] Figure 3 Schematic diagram of the process of S3 of the present invention;

[0058] Figure 4 Schematic diagram of the process of S5 of the present invention. DETAILED DESCRIPTION

[0059] like Figures 1 to 4 As shown, the present invention discloses a method for predicting the acoustic compatibility of a whole ship based on time-frequency space energy, comprising the following steps:

[0060] S1. Establish an acoustic compatibility technical parameter library based on multi-dimensional characteristics, including the acoustic compatibility technical parameters of the interference source and the acoustic compatibility technical parameters of the interfered acoustic equipment.

[0061] The acoustic compatibility multi-dimensional technical parameters of the interference source include the energy domain, time domain, frequency domain, and spatial domain of the interference source. The energy domain of the interference source represents the acoustic source level of the transmitted signal; the time domain of the interference source represents the pulse width of the transmitted signal; the frequency domain of the interference source represents the frequency, bandwidth, and out-of-band attenuation of the transmitted signal; and the spatial domain of the interference source represents the directivity of the transmitted signal and the main-sidelobe suppression ratio in the direction of the interfered device.

[0062] The acoustic compatibility technical parameters of the interfered acoustic equipment include the acoustic compatibility technical parameters of the receiving sensor, the acoustic compatibility technical parameters of the front-end receiving channel, and the acoustic compatibility technical parameters of the digital signal processing. The acoustic compatibility technical parameters of the receiving sensor include the spatial domain and frequency domain of the receiving sensor; the spatial domain of the receiving sensor represents the spatial directivity of the receiving sensor; and the frequency domain of the receiving sensor represents the sensitivity of the receiving sensor. The acoustic compatibility technical parameters of the front-end receiving channel include the frequency domain and energy domain of the front-end receiving channel; the frequency domain of the front-end receiving channel represents the frequency response characteristics of each stage of amplification and filtering; and the energy domain of the front-end receiving channel represents the limiting voltage of each stage of amplification and filtering. The acoustic compatibility technical parameters of the digital signal processing include the frequency domain, spatial domain, and time domain of the digital signal processing; the frequency domain of the digital signal processing represents the frequency domain interference suppression capability, the spatial domain of the digital signal processing represents the spatial domain interference suppression capability, and the time domain of the digital signal processing represents the time domain interference suppression capability.

[0063] By establishing an acoustic compatibility technical parameter library based on multi-dimensional characteristics, it comprehensively covers the multi-dimensional characteristic parameters of the interference source and interfered acoustic equipment in the energy domain, time domain, frequency domain and spatial domain, providing a more comprehensive and accurate data basis for acoustic compatibility prediction, and effectively making up for the shortcomings of traditional methods in comprehensive consideration of multi-dimensional characteristics.

[0064] S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device based on the energy domain characteristics, spatial domain characteristics, and interference propagation attenuation of the interference source acoustic device.

[0065] Specifically, S2 includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, spatial domain characteristics, and interference propagation attenuation of the interference source acoustic device, using the following formula:

[0066]

[0067] in, is the interference level at the interfered acoustic device, SL is the emission sound source level of the interference source; DIi is the emission directivity of the interference source. When the interfered acoustic device is in the main lobe azimuth of the interference source, DIi=0. When it is in the side lobe azimuth, it is the main lobe to side lobe ratio; TL is the interference propagation attenuation, including distance propagation loss and hull structure shielding. The formula is as follows:

[0068]

[0069] Among them, g represents the acceleration of gravity, R represents the distance between the two devices, represents the seawater absorption coefficient, Indicates the attenuation of interference caused by hull structure obstruction.

[0070] Based on the multi-dimensional characteristics of the interference source and the interference propagation attenuation, the interference level at the interfered acoustic equipment is accurately calculated. Factors such as distance propagation loss, seawater absorption, and hull structure obstruction are comprehensively considered, which improves the accuracy and reliability of the interference level calculation and provides strong support for subsequent interference assessment and equipment design.

[0071] S3. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device based on the interference level at the interfered acoustic device and relevant acoustic compatibility technical parameters, and determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device based on the interference intensity.

[0072] Said S3 specifically includes:

[0073] S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device:

[0074] The intensity of interference picked up by the receiving transducer of the interfered acoustic device The formula is as follows:

[0075]

[0076] Where P represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; N represents the attenuation of the interference source direction compared to the detection direction of the interfered acoustic device;

[0077] S32. Determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device:

[0078] According to the frequency response characteristics of amplification and filtering in the front-end receiving channel of the interfered acoustic device, the output of each level of the front-end receiving channel of the interfered acoustic device is calculated and compared with the corresponding limit voltage. If the limit voltage is exceeded, it indicates saturation; if not, it indicates non-saturation;

[0079] Use the following formula to determine whether the frequency response characteristics of the preamplifier are saturated:

[0080]

[0081] in, is the frequency response characteristic of the preamplifier, Indicates the limiting voltage of the preamplifier; when ≥0, it means saturation occurs. When <0, no saturation occurs;

[0082] Use the following formula to determine whether the pre-filter is saturated:

[0083]

[0084] in, is the frequency response characteristic of the pre-filter, Indicates the limiting voltage of the pre-filter; when ≥0, it means saturation occurs. No saturation occurs.

[0085] The interference intensity entering the front-end receiving channel of the interfered acoustic device is systematically calculated, and saturation judgment is made based on the frequency response characteristics of amplification and filtering. The judgment criteria are clarified, providing clear guidance for the design and optimization of the front-end receiving channel of the device, which helps to improve the stability and reliability of the device.

[0086] S4. Determine whether the interfered acoustic device generates in-band interference based on the multi-dimensional characteristics of acoustic compatibility.

[0087] The S4 specifically includes:

[0088] Based on the multi-dimensional characteristics of the interference source and the interfered acoustic device in the time domain, frequency domain, and spatial domain, the following formula is used to determine whether the interfered acoustic device generates in-band interference:

[0089]

[0090] Where: M represents the background noise level when the interfered acoustic equipment is demonstrated, Y represents the emission bandwidth of the interference source, g represents the acceleration of gravity, Indicates the attenuation of the interference source's transmission band compared to the interference frequency. Indicates the frequency domain interference suppression capability of the interfered acoustic equipment, Indicates the spatial interference suppression capability of the interfered acoustic equipment. Indicates the time domain interference suppression capability of the interfered acoustic device.

[0091] Calculate the frequency domain interference suppression capability based on the frequency domain processing method of the interfered acoustic device , when using line spectrum processing, =0; when broadband processing is used, The calculation method is as follows:

[0092]

[0093] in, represents the processing bandwidth of the interfered acoustic device, Indicates the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic equipment;

[0094] Calculation of spatial interference suppression based on the receiving directivity of the interfered acoustic device , when the interference source is in the detection main lobe direction of the interfered acoustic equipment, =0; when in the side lobe, The main-side lobe ratio of this azimuth is greater than the main lobe;

[0095] Calculates time-domain interference suppression based on the overlap between the interference source's transmission and the victim's acoustic device's processing period, as well as the victim's acoustic device's processing integration time. , The calculation method is shown in the following formula:

[0096]

[0097] in, Indicates the overlap duration, Indicates the integration duration of the disturbed acoustic device.

[0098] S5. Based on the calculation and judgment results of S3 and S4, the comprehensive impact of the interference is quantitatively evaluated to generate the full-ship acoustic compatibility prediction results of time-frequency space energy.

[0099] The correlation between saturation and in-band interference can also be analyzed, including: Saturation-dominated interference: If the receiving channel is saturated (S3 determines yes), but there is no in-band interference (S4 determines no), signal dynamic range compression leads to a degraded signal-to-noise ratio, affecting target detection capabilities. In-band interference-dominated interference: If in-band interference is present (S4 determines yes) but not saturated (S3 determines no), interference energy directly overwhelms the target frequency band, increasing the false alarm rate. Combined interference: If both saturation and in-band interference conditions are met (S3 and S4 both determine yes), signal distortion and spectrum contamination are superimposed, leading to a nonlinear breakdown in system performance.

[0100] The S5 specifically includes: S51, constructing a comprehensive interference impact index: as shown in the following formula:

[0101]

[0102] in, and is the weight coefficient, =0.6, =0.4, is the calculation result of S3, and is the corresponding limiting voltage threshold, K is the calculation result of S4, and NL is the background noise level; and is the nonlinear correction factor, =1.2, =1.5.

[0103] It also includes S52, which divides the interference level based on the CII value. If the CII value is less than 0.3, it is a slight interference, and only the interference event needs to be recorded without active suppression; if 0.3≤CII<0.7, it is moderate interference, and frequency domain notching or spatial domain nulling is enabled; if CII≥0.7, it is severe interference, and the device's working frequency band is forced to switch or time domain asynchronous sampling is started.

[0104] It also includes S53, confidence assessment, which introduces Monte Carlo simulation to verify the reliability of the results, as shown below:

[0105]

[0106] in, represents the mean value of CII, which is obtained by randomly perturbing the input parameters multiple times. is the standard deviation of CII;

[0107] If the confidence level is less than 0.8, the process returns to S51 to recalibrate the parameters.

[0108] By constructing a comprehensive interference impact index, the saturation of the receiving channel and the impact of in-band interference are quantitatively integrated, providing an intuitive evaluation indicator for acoustic compatibility prediction, which helps to fully understand and optimize the performance of acoustic equipment throughout the ship.

[0109] There is feedback optimization for acoustic compatibility: including avoiding saturation between devices: based on the judgment result of S32, if saturation occurs between devices, the acoustic compatibility technical parameters of both the interference source and the interfered device are optimized; for the interference source, the direction of its transmitted signal is optimized to improve the main sidelobe suppression ratio; for the interfered device, including increasing the number of amplifier stages, increasing the out-of-band attenuation of each stage of filter and increasing the limiting or saturation voltage of each stage of amplification and filtering, among which increasing the number of amplifier stages means reducing the total amplification of the front-end receiving channel on the premise that the total amplification of the front-end receiving channel meets the demand, and adjusting the amplification to the back stage.

[0110] Based on the judgment result of S4, if the interfered acoustic device generates in-band interference, an anti-in-band interference method is used to improve the signal quality and signal-to-noise ratio.

[0111] The anti-in-band interference method specifically includes:

[0112] Acquire reference signals and main signals;

[0113] Signal decomposition: Use analysis filter banks to decompose the reference signal and main signal Decomposed into W sub-band signals as follows:

[0114]

[0115]

[0116] in, is the subband signal of the reference signal, is the subband signal of the main signal, is the filter coefficient for the kth subband, calculated using an orthogonal mirror filter or Fourier transform, where I is the order of the analysis filter and k is the subband index;

[0117] Subband adaptive filtering: for each subband signal and Use the NLMS algorithm independently as follows:

[0118]

[0119]

[0120]

[0121] in, is the time domain coefficient of the adaptive filter for the kth subband, is the noise estimate of the kth subband, is the error signal of the kth subband, that is, the difference between the main signal and the noise estimate; is the step size factor, and its value range is ; is the regularization constant, is the energy of the reference signal vector;

[0122] Signal reconstruction: The processed sub-band signals are combined into the final output through the synthesis filter bank as shown below:

[0123]

[0124] in, is the synthesis filter coefficient of the kth subband, is the final output noise reduction signal.

[0125] Finally, it should be noted that the above preferred embodiments are only used to illustrate the technical solutions of the present invention and are not limiting. Although the present invention has been described in detail through the above preferred embodiments, those skilled in the art should understand that various changes can be made in form and details without departing from the scope defined by the claims of the present invention.

Claims

1. A method for predicting the acoustic compatibility of a whole ship based on time-frequency space energy, characterized by: The following steps are involved: S1. Establish an acoustic compatibility technical parameter library based on multi-dimensional characteristics, including the acoustic compatibility technical parameters of the interference source and the acoustic compatibility technical parameters of the interfered acoustic equipment; S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device based on the energy domain characteristics, spatial domain characteristics, and interference propagation attenuation of the interference source acoustic device; S3. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device based on the interference level at the interfered acoustic device and relevant acoustic compatibility technical parameters, and determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device based on the interference intensity; S4. Determine whether the interfered acoustic device generates in-band interference based on the multi-dimensional characteristics of acoustic compatibility; S5. Based on the calculation and judgment results of S3 and S4, quantitatively evaluate the comprehensive impact of interference and generate the full-ship acoustic compatibility prediction results of time-frequency space energy; Said S3 specifically includes: S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device: The formula for the interference intensity H picked up by the receiving transducer of the interfered acoustic device is as follows: H=IL+PN Where P represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; N represents the attenuation of the interference source direction compared to the detection direction of the interfered acoustic device; IL is the interference level at the interfered acoustic device; S32. Determine whether saturation occurs at each level of the front-end receiving channel of the interfered acoustic device: According to the frequency response characteristics of amplification and filtering in the front-end receiving channel of the interfered acoustic device, the output of each level of the front-end receiving channel of the interfered acoustic device is calculated and compared with the corresponding limit voltage. If the limit voltage is exceeded, it indicates saturation; if not, it indicates non-saturation; Use the following formula to determine whether the frequency response characteristics of the preamplifier are saturated: △A=H+X(f)-t1 Where X(f) is the frequency response characteristic of the preamplifier, t1 represents the limiting voltage of the preamplifier; when ΔA ≥ 0, saturation occurs, and when ΔA < 0, no saturation occurs; Use the following formula to determine whether the pre-filter is saturated: △B=H+X(f)+D(f)-t2 Where D(f) is the frequency response characteristic of the prefilter, t2 represents the limiting voltage of the prefilter; when ΔB ≥ 0, saturation occurs, and when ΔB < 0, no saturation occurs; The S4 specifically includes: Based on the multi-dimensional characteristics of the interference source and the interfered acoustic device in the time domain, frequency domain, and spatial domain, the following formula is used to determine whether the interfered acoustic device generates in-band interference: K=IL-101gY-D f -(X f +X s +X t )-M Where: M represents the background noise level when the interfered acoustic equipment is demonstrated, Y represents the interference source emission bandwidth, g represents the acceleration of gravity, D f Indicates the attenuation of the interference source's transmission band compared to the interference frequency, X f It represents the frequency domain interference suppression capability of the interfered acoustic equipment, X s It represents the spatial interference suppression capability of the interfered acoustic equipment, X t Indicates the time domain interference suppression capability of the interfered acoustic equipment; Calculate the frequency domain interference suppression capability X based on the frequency domain processing method of the interfered acoustic device f , when line spectrum processing is used, X f =0; when broadband processing is used, X f The calculation method is as follows: X f =101g(V s / V i ) Among them, V s represents the processing bandwidth of the interfered acoustic device, V i Indicates the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic equipment; Calculate spatial interference suppression X based on the receiving directivity of the interfered acoustic device s , when the interference source is in the detection main lobe direction of the interfered acoustic equipment, X s =0; when in the side lobe, X s The main-side lobe ratio of this azimuth is greater than the main lobe; The time domain interference suppression X is calculated based on the overlapping time of the interference source transmission and the processing period of the interfered acoustic device and the processing integration time of the interfered acoustic device. t , X t The calculation method is shown in the following formula: X t =101g(Z r / Z c ) Among them, Z c Indicates the overlap duration, Z r Indicates the integration time of the disturbed acoustic device; The S5 specifically includes: S51, constructing a comprehensive interference impact index: as shown in the following formula: Wherein, λ and β are weight coefficients, λ = 0.6, β = 0.4, △A and △B are the calculation results of S3, t1 and t2 are the corresponding limiting voltage thresholds, K is the calculation result of S4, NL is the background noise level; γ and k are nonlinear correction factors, γ = 1.2, κ = 1.5, CII is the comprehensive interference impact index; The S5 also includes S52, dividing the interference level based on the CII value. If the CII value is less than 0.3, it is a slight interference, and only the interference event needs to be recorded without active suppression; if 0.3≤CII<0.7, it is moderate interference, and frequency domain notching or spatial domain nulling is enabled; if CII≥0.7, it is severe interference, and the device working frequency band is forced to switch or time domain asynchronous sampling is started.

2. The method for predicting the whole-ship acoustic compatibility based on time-frequency space energy according to claim 1, characterized in that: The acoustic compatibility technical parameters of the interference source described in S1 include the energy domain of the interference source, the time domain of the interference source, the frequency domain of the interference source, and the spatial domain of the interference source; the acoustic compatibility technical parameters of the interfered acoustic device include the acoustic compatibility technical parameters of the receiving sensor, the acoustic compatibility technical parameters of the front-end receiving channel, and the acoustic compatibility technical parameters of the digital signal processing; The acoustic compatibility technical parameters of the receiving sensor include the spatial domain and frequency domain of the receiving sensor; the acoustic compatibility technical parameters of the front-end receiving channel include the frequency domain and energy domain of the front-end receiving channel; and the acoustic compatibility technical parameters of digital signal processing include the frequency domain, spatial domain and time domain of digital signal processing.

3. The method for predicting the whole-ship acoustic compatibility based on time-frequency space energy according to claim 1, characterized in that: Specifically, S2 includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, spatial domain characteristics, and interference propagation attenuation of the interference source acoustic device, using the following formula: IL=SL-DIi-TL Where IL is the interference level at the interfered acoustic device, SL is the emission sound source level of the interference source; DIi is the emission directivity of the interference source. When the interfered acoustic device is in the main lobe azimuth of the interference source, DIi = 0. When it is in the side lobe azimuth, it is the main lobe to side lobe ratio; TL is the interference propagation attenuation, including distance propagation loss and hull structure obstruction. The formula is as follows: TL=201gR+αR+δ Where g represents the acceleration of gravity, R represents the distance between the two devices, α represents the seawater absorption coefficient, and δ represents the attenuation of interference caused by the hull structure.

4. The method for predicting the whole-ship acoustic compatibility based on time-frequency space energy according to claim 1 is characterized by: The S5 also includes S53, confidence assessment, which introduces Monte Carlo simulation to verify the reliability of the results, as shown below: Where μCII represents the mean of CII, which is calculated multiple times by randomly perturbing the input parameters, and σCII is the standard deviation of CII; If the confidence level is less than 0.8, the process returns to S51 to recalibrate the parameters.

Images