Whole-ship acoustic compatibility prediction method based on time-frequency space energy

Abstract

The invention discloses a whole ship acoustic compatibility prediction method based on time-frequency space energy, and belongs to the technical field of acoustic compatibility. Comprising the following steps: S1, establishing an acoustic compatibility technical parameter library based on multi-dimensional characteristics; s2, calculating the interference level at the interfered acoustic equipment; S3, calculating the interference intensity entering the front-end receiving channel of the interfered acoustic equipment, and judging whether each level of the front-end receiving channel of the interfered acoustic equipment is saturated or not based on the interference intensity; s4, judging whether the interfered acoustic equipment generates in-band interference or not based on the acoustic compatible multi-dimensional characteristics; and S5, based on calculation and judgment results of S3 and S4, carrying out quantitative evaluation on the comprehensive influence of interference, and generating a whole-ship acoustic compatibility prediction result of time-frequency space energy. The acoustic compatibility state between the shipborne acoustic equipment can be predicted more accurately, and the ship can be supported to carry out acoustic compatibility design more reasonably on the whole. And the whole ship is effectively supported, the acoustic compatibility design is carried out in the initial design stage of equipment, and the acoustic compatibility of the whole ship is optimized.

Description

Technical Field

[0001] The present invention belongs to the technical field of acoustic compatibility, and particularly relates to a full-ship acoustic compatibility prediction method based on time-frequency-space-energy. Background Art

[0002] With the continuous development and application of full-ship acoustic devices, the problem of acoustic compatibility has become increasingly prominent, becoming a key factor affecting the ship's efficiency and operation safety. Traditional acoustic compatibility analysis methods are mainly based on single-dimensional parameters, such as frequency, time, or space, etc., and it is difficult to comprehensively and accurately evaluate and predict the acoustic interference situation between multiple devices. For example, some methods only focus on the frequency-domain characteristics of the interference source, while ignoring the influence of time-domain and space-domain characteristics on acoustic compatibility, resulting in insufficient prediction accuracy in complex environments. Moreover, in the previous full-ship acoustic compatibility design, it was mainly based on the interference magnitude, and the acoustic compatibility between devices was estimated by using the far-field theory calculation and simulation methods: However, due to the increase in the underwater functions of the whole ship, the number of shipborne acoustic devices has increased, the working frequency bands have expanded, the working modes have become rich, and the problem of working frequency band overlap is prominent. The traditional analysis method based on the interference magnitude can no longer cope with the complex acoustic compatibility form, and it is difficult to comprehensively and accurately predict the acoustic compatibility state between shipborne acoustic devices, and it is difficult to put forward correct and reasonable anti-interference measures and anti-interference technical requirements in the overall ship and equipment design stage, and it cannot support the full-ship acoustic compatibility design.

[0003] Therefore, there is an urgent need 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 full-ship acoustic compatibility prediction method based on time-frequency-space-energy, which can be used to predict the acoustic interference situation between shipborne acoustic devices, support the proposal of quantitative acoustic compatibility technical requirements and preliminary anti-interference measures, support the overall ship and acoustic devices to carry out acoustic compatibility design, is conducive to improving the full-ship acoustic compatibility, solves the problem of no basis for the overall ship, system, and acoustic device to carry out acoustic compatibility design, and improves the underwater efficiency.

[0005] To achieve the above purpose, the present invention provides the following technical solutions: A full-ship acoustic compatibility prediction method based on time-frequency-space-energy, comprising the following steps: 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 device; S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device according to the energy-domain characteristics, space-domain characteristics of the interference source acoustic device, and the interference propagation attenuation; 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 the relevant acoustic compatibility technical parameters, and determine whether saturation occurs at each stage of the front-end receiving channel of the interfered acoustic device based on the interference intensity; S4. Judge 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 the interference, and generate the full-ship acoustic compatibility prediction result of time-frequency-space-energy.

[0006] As a further preference of the present invention, the multi-dimensional acoustic compatibility technical parameters of the interference source 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 space 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 digital signal processing; among them, the acoustic compatibility technical parameters of the receiving sensor include the space 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 digital signal processing include the frequency domain, space domain, and time domain of digital signal processing.

[0007] As a further preference of the present invention, S2 specifically includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, space domain characteristics of the interfering source acoustic device, and interference propagation attenuation. The formula is as follows: Among them, is the interference level at the interfered acoustic device, SL represents the emission sound source level of the interference source; DIi represents the emission directivity of the interference source. When the interfered acoustic device is in the azimuth of the main lobe of the interference source emission, DIi = 0, and when it is in the azimuth of the side lobe, it is the main side lobe ratio; TL represents the interference propagation attenuation, including distance propagation loss and hull structure shielding. The formula is as follows: Among them, g represents the acceleration of gravity, R represents the distance between the two devices, represents the seawater absorption coefficient, represents the attenuation amount of the hull structure shielding to the interference.

[0008] As a further preference of the present invention, S3 specifically includes, S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device: The interference intensity picked up by the receiving transducer of the interfered acoustic device The formula is as follows: Wherein, P represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; N represents the attenuation amount of the interference source azimuth relative to the detection azimuth of the interfered acoustic device. S32. Determine whether saturation occurs at each stage 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, calculate the output of each stage of the front-end receiving channel of the interfered acoustic device and compare it with the corresponding limiting voltage. If it exceeds the limiting voltage, it indicates saturation; if it does not exceed, it indicates no saturation. Use the following formula to determine whether saturation occurs in the frequency response characteristic of the preamplification: Wherein, is the frequency response characteristic of the preamplification, represents the limiting voltage of the preamplification; when ≥0, it indicates saturation; when <0, no saturation occurs. Use the following formula to determine whether saturation occurs in the pre-filtering: Wherein, is the frequency response characteristic of the pre-filtering, represents the limiting voltage of the pre-filtering; when ≥0, it indicates saturation; when <0, no saturation occurs.

[0009] As a further preference of the present invention, the S4 specifically includes: According to the multi-dimensional characteristics of time domain, frequency domain and space domain of the interference source and the interfered acoustic device, the judgment of whether in-band interference occurs in the interfered acoustic device is shown in the following formula: Wherein: M represents the background noise level during the index demonstration of the interfered acoustic device, Y represents the emission bandwidth of the interference source, g represents the acceleration of gravity, represents the attenuation amount of the interference source emission outside the band relative to the interference frequency, represents the frequency domain interference suppression ability of the interfered acoustic device, represents the space domain interference suppression ability of the interfered acoustic device, represents the time domain interference suppression ability of the interfered acoustic device; Calculate the frequency domain interference suppression ability based on the frequency domain processing method of the interfered acoustic device. When using line spectrum processing, =0; when using broadband processing, The calculation method of Among them, represents the processing bandwidth of the interfered acoustic device, represents the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic device; Calculate the spatial domain interference suppression based on the receiving directivity of the interfered acoustic device , when the interference source is in the azimuth of the detection main lobe of the interfered acoustic device, =0; when in the side lobe, is the main side lobe ratio of this azimuth compared to the main lobe; Calculate the time domain interference suppression based on the overlapping duration of the interference source emission and the processing time period of the interfered acoustic device and the integration time of the interfered acoustic device , The calculation method of is as shown in the following formula: Among them, represents the overlapping duration, represents the integration duration of the interfered acoustic device.

[0010] As a further preference of the present invention, the S5 specifically includes: S51. Construct an interference comprehensive influence index: as shown in the following formula: Among them, and are weight coefficients, =0.6, =0.4, is the calculation result of S3, and are the corresponding limiting voltage thresholds, K is the calculation result of S4, and NL is the background noise level; and are non-linear correction factors, =1.2, =1.5.

[0011] As a further preference of the present invention, the S5 further includes S52. Divide the interference level based on the CII value. If the CII value < 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 a moderate interference, and frequency domain notch filtering or spatial domain nulling is enabled. If CII ≥ 0.7, it is a severe interference, and the device working frequency band is forced to be switched or time domain asynchronous sampling is started.

[0012] As a further preference of the present invention, the S5 further includes S53. Confidence evaluation, introducing Monte Carlo simulation to verify the reliability of the results, as follows: Among them, represents the CII mean value, which is obtained by calculating multiple times by randomly perturbing the input parameters. is the standard deviation of CII; If the confidence level < 0.8, then return to S51 to recalibrate the parameters.

[0013] The beneficial effects of the present invention are as follows: The present invention proposes a full-ship acoustic compatibility prediction method based on time-frequency-space-energy, including the judgment of saturation and in-band interference generated by the front-end receiving channel of the interfered acoustic device, which involves multi-dimensional characteristics such as time, frequency, space, etc. of the interference source, the interfered acoustic device's transmission and reception, and factors such as interference propagation attenuation. To cope with the current complex acoustic compatibility situation, the acoustic interference between devices is analyzed from multiple dimensions such as time, frequency, and space to more accurately predict the acoustic compatibility state between shipborne acoustic devices, support the ship's overall more reasonable arrangement of acoustic devices, study and adopt anti-interference isolation measures, put forward quantitative anti-interference technical requirements for shipborne acoustic devices, etc. At the same time, it can support the demonstration and design of shipborne acoustic devices, more reasonably set the working frequency band, design the acoustic compatibility technical indicators of the transmitting end and the receiving end, reserve management control ports, etc.; more accurately predict the acoustic compatibility state between shipborne acoustic devices, quantify the interference situation between shipborne acoustic devices, and support the subsequent targeted preliminary acoustic compatibility anti-interference measures for saturation and in-band interference. It effectively supports the overall and device acoustic compatibility design in the initial stage of design, realizes the optimization of the full-ship acoustic compatibility. It achieves "early discovery, early adjustment, and early response" to acoustic compatibility problems, avoids serious problems such as "inadequate congenital" acoustic compatibility between shipborne acoustic devices, such as the same frequency, repeated frequency, and small frequency interval, optimizes the full-ship acoustic compatibility, and ensures the full play of the underwater performance of the ship.

[0014] Other advantages, objectives, and features of the present invention will be described in the subsequent specification, and to some extent, they are obvious to those skilled in the art, or those skilled in the art can obtain teachings from the practice of the present invention. The objectives and other advantages of the present invention can be realized and obtained through the following specification. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] To make the objectives, technical solutions, and beneficial effects of the present invention clearer, the present invention provides the following drawings for illustration: Figure 1 is the overall flowchart of a full-ship acoustic compatibility prediction method based on time-frequency-space-energy of the present invention; Figure 2 is a schematic diagram of the acoustic compatibility technical parameter library of the present invention; Figure 3 is the schematic flowchart of S3 of the present invention; Figure 4 is the schematic flowchart of S5 of the present invention. Detailed implementation mode

[0016] As Figures 1 to 4 shown, the present invention discloses a full-ship acoustic compatibility prediction method based on time-frequency-space-energy, including the following steps: S1. Establish an acoustic compatibility technical parameter library based on multi-dimensional characteristics, including the acoustic compatibility technical parameters of interference sources and the acoustic compatibility technical parameters of interfered acoustic devices.

[0017] The multi-dimensional technical parameters of the acoustic compatibility of the interference source include the energy domain of the interference source, the time domain of the interference source, the frequency domain of the interference source, and the space domain of the interference source. The energy domain of the interference source represents the sound 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 transmitted signal frequency, bandwidth, and out-of-band attenuation; the space domain of the interference source represents the directivity of the transmitted signal and the main side lobe suppression ratio in the direction of the interfered device.

[0018] 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 digital signal processing. Among them, the acoustic compatibility technical parameters of the receiving sensor include the space domain and frequency domain of the receiving sensor; the space domain of the receiving sensor represents the spatial directivity of the receiving sensor; 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; 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 digital signal processing include the frequency domain, space domain, and time domain of digital signal processing; the frequency domain of digital signal processing represents the frequency domain interference suppression ability, the space domain of digital signal processing represents the space domain interference suppression ability, and the time domain of digital signal processing represents the time domain interference suppression ability.

[0019] By establishing an acoustic compatibility technical parameter library based on multi-dimensional characteristics, it comprehensively covers multi-dimensional characteristic parameters such as the energy domain, time domain, frequency domain, and space domain of interference sources and interfered acoustic devices, providing a more comprehensive and accurate data basis for acoustic compatibility prediction and effectively making up for the deficiencies of traditional methods in comprehensively considering multi-dimensional characteristics.

[0020] S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device according to the energy domain characteristics, space domain characteristics of the interference source acoustic device, and interference propagation attenuation.

[0021] Specifically, S2 includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, space domain characteristics of the interference source acoustic device, and interference propagation attenuation. The formula is as follows: Wherein, $L_{I}$ is the interference level at the interfered acoustic device, $SL$ represents the transmitting sound source level of the interference source; $DI_{i}$ represents the transmitting directivity of the interference source. When the interfered acoustic device is in the azimuth of the main lobe of the interference source's transmission, $DI_{i}=0$. When it is in the azimuth of the side lobe, it is the main-to-side lobe ratio; $TL$ represents the interference propagation attenuation, including distance propagation loss and hull structure shielding. The formula is as follows: Among them, $g$ represents the acceleration due to gravity, $R$ represents the distance between the two devices, $\alpha$ represents the seawater absorption coefficient, $A_{hull}$ represents the attenuation of the hull structure shielding to the interference.

[0022] According to the multi-dimensional characteristics of the interference source and the interference propagation attenuation, the interference level at the interfered acoustic device is accurately calculated, comprehensively considering factors such as distance propagation loss, seawater absorption, and hull structure shielding, improving the accuracy and reliability of the interference level calculation, and providing strong support for subsequent interference assessment and equipment design.

[0023] 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 judge whether each stage of the front-end receiving channel of the interfered acoustic device is saturated based on the interference intensity.

[0024] The specific content of S3 includes, S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device: The interference intensity picked up by the receiving transducer of the interfered acoustic device The formula is as follows: Among them, $P$ represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; $N$ represents the attenuation amount of the interference source azimuth relative to the detection azimuth of the interfered acoustic device; S32. Judge whether each stage of the front-end receiving channel of the interfered acoustic device is saturated: According to the frequency response characteristics of amplification and filtering in the front-end receiving channel of the interfered acoustic device, calculate the output of each stage of the front-end receiving channel of the interfered acoustic device and compare it with the corresponding limiting voltage. If it exceeds the limiting voltage, it means saturation occurs. If it does not exceed, it means unsaturation; The following formula is used to judge whether the frequency response characteristics of the pre-amplification are saturated: Among them, $G_{pre}$ is the frequency response characteristic of the pre-amplification, $V_{lim}$ represents the limiting voltage of the pre-amplification; when $G_{pre}\geq0$, it means saturation occurs. When $G_{pre}<0$, no saturation occurs; The following formula is used to determine whether pre-filtering causes saturation: where, is the frequency response characteristic of the pre-filter, represents the limiting voltage of the pre-filter; when ≥0, it indicates saturation, and when , there is no saturation.

[0025] The interference intensity entering the front-end receiving channel of the interfered acoustic device is systematically calculated, and saturation judgment is carried out based on the frequency response characteristics of amplification and filtering, clarifying the judgment criteria, providing clear guidance for the design and optimization of the front-end receiving channel of the device, and helping to improve the stability and reliability of the device.

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

[0027] The specific steps of S4 include: According to the multi-dimensional characteristics of the time domain, frequency domain, and space domain of the interference source and the interfered acoustic device, the judgment of whether the interfered acoustic device generates in-band interference is shown in the following formula: where: M represents the background noise level during the index demonstration of the interfered acoustic device, Y represents the emission bandwidth of the interference source, g represents the acceleration of gravity, represents the attenuation of the interference source outside the emission band compared to the interference frequency, represents the frequency domain interference suppression ability of the interfered acoustic device, represents the space domain interference suppression ability of the interfered acoustic device, represents the time domain interference suppression ability of the interfered acoustic device.

[0028] Calculate the frequency domain interference suppression ability based on the frequency domain processing method of the interfered acoustic device , when using line spectrum processing, =0; when using broadband processing, The calculation method of is shown in the following formula: where, represents the processing broadband of the interfered acoustic device, represents the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic device; Calculate the space domain interference suppression based on the receiving directivity of the interfered acoustic device , when the interference source is in the azimuth of the detection main lobe of the interfered acoustic device, =0; when in the side lobe, For the main side lobe ratio of this azimuth to the main lobe; Calculate the time-domain interference suppression based on the overlapping duration of the interference source emission and the processing period of the interfered acoustic device and the processing integration time of the interfered acoustic device , The calculation method is as shown in the following formula: Wherein, represents the overlapping duration, represents the integration duration of the interfered acoustic device.

[0029] 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.

[0030] It is also possible to analyze the correlation between saturation and in-band interference, including: saturation-dominated interference: If the receiving channel is saturated (S3 determines yes), that is, there is no in-band interference (S4 determines no), the signal dynamic range compression will cause the signal-to-noise ratio to deteriorate, affecting the target detection ability. In-band interference-dominated: If there is in-band interference (S4 determines yes) but not saturated (S3 determines no), the interference energy directly covers the target frequency band, resulting in an increase in the false alarm rate. Composite interference mode: If both the saturation and in-band interference conditions are met (both S3 and S4 are yes), the signal distortion and spectrum pollution are superimposed, and the system performance shows a non-linear collapse trend.

[0031] The specific steps of S5 include: S51. Construct an interference comprehensive impact index: as shown in the following formula: Wherein, and are weight coefficients, = 0.6, = 0.4, is the calculation result of S3, and are the corresponding limiting voltage thresholds, K is the calculation result of S4, and NL is the background noise level; and are non-linear correction factors, = 1.2, = 1.5.

[0032] It also includes S52. Divide the interference level based on the CII value. If the CII value < 0.3, it is a minor interference, and only the interference event needs to be recorded without active suppression; if 0.3 ≤ CII < 0.7, it is a moderate interference, and frequency-domain notch or spatial nulling is enabled. If CII ≥ 0.7, it is a severe interference, and the device working frequency band is forced to be switched or time-domain asynchronous sampling is started.

[0033] It also includes S53, confidence evaluation, and introduces Monte Carlo simulation to verify the reliability of the results, as follows: Among them, represents the mean value of CII, which is obtained by calculating multiple times with randomly perturbed input parameters, is the standard deviation of CII; If the confidence level < 0.8, then return to S51 to recalibrate the parameters.

[0034] By constructing the interference comprehensive influence index, the saturation situation of the receiving channel and the influence of in-band interference are quantitatively integrated, providing an intuitive evaluation index for acoustic compatibility prediction, which helps to comprehensively understand and optimize the performance of the ship's acoustic equipment.

[0035] There is feedback optimization for acoustic compatibility: including avoiding saturation between devices: based on the judgment result of S32, if saturation occurs between devices, then optimize the acoustic compatibility technical parameters for both the interference source and the interfered device; for the interference source, optimize the direction of its transmitted signal to improve the main side lobe suppression ratio; for the interfered device, it includes increasing the number of amplifier stages, improving the out-of-band attenuation of each stage of the filter, and increasing the limiting or saturation voltage of each stage of amplification and filtering, where 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 requirements, and adjusting the amplification to the subsequent stage.

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

[0037] The anti-in-band interference method specifically includes: Collect the reference signal and the main signal; Signal decomposition: Use the analysis filter bank to decompose the reference signal and the main signal into W sub-band signals, as follows: Among them, is the sub-band signal of the reference signal, is the sub-band signal of the main signal, is the filter coefficient of the k-th sub-band, which is calculated using the quadrature mirror filter or the Fourier transform, I is the order of the analysis filter, and k is the sub-band index; Sub-band adaptive filtering: Independently use the NLMS algorithm for each sub-band signal and , as follows: wherein, are the time-domain coefficients of the adaptive filter for the k-th sub-band, is the noise estimation value for the k-th sub-band, is the error signal for the k-th sub-band, i.e., the difference between the main signal and the noise estimation; is the step size factor, and its value range is ; is the regularization constant, is the energy of the reference signal vector; Signal reconstruction: The processed sub-band signals are combined into the final output through a synthesis filter bank as follows: wherein, are the synthesis filter coefficients for the k-th sub-band, is the noise-reduced signal of the final output.

[0038] Finally, it should be noted that the above preferred embodiments are only used to illustrate the technical solutions of the present invention rather than to limit them. 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 to it in form and details without departing from the scope defined by the claims of the present invention.

Claims

1. A full-ship acoustic compatibility prediction method based on time-frequency-space-energy, characterized in that: It includes the following steps: S1. Establish an acoustic compatibility technical parameter library based on multi-dimensional characteristics, including the acoustic compatibility technical parameters of interference sources and the acoustic compatibility technical parameters of the interfered acoustic devices; S2. Calculate the interference level at the interfered acoustic device: Calculate the interference level at the interfered acoustic device according to the energy domain characteristics, spatial domain characteristics of the interference source acoustic device and the interference propagation attenuation; 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 the relevant acoustic compatibility technical parameters, and judge whether each stage of the front-end receiving channel of the interfered acoustic device is saturated based on the interference intensity; S4. Judge 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 influence of the interference, and generate the full-ship acoustic compatibility prediction result of time-frequency-space-energy.

2. The full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 1, characterized in that: The multi-dimensional technical parameters of acoustic compatibility of the interference source 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 digital signal processing; Among them, 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 digital signal processing include the frequency domain, spatial domain and time domain of digital signal processing.

3. A full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 1, characterized in that: The specific content of S2 includes calculating the interference level at the interfered acoustic device according to the energy domain characteristics, spatial domain characteristics of the interference source acoustic device and the interference propagation attenuation. The formula is as follows: Among them, is the interference level at the interfered acoustic device, SL represents the emission sound source level of the interference source; DIi represents the emission directivity of the interference source. When the interfered acoustic device is in the azimuth of the main lobe of the interference source's emission, DIi = 0, and when it is in the azimuth of the side lobe, it is the main-to-side lobe ratio; TL represents the interference propagation attenuation, including distance propagation loss and hull structure shielding. The formula is as follows: where g represents the acceleration due to gravity and R represents the distance between the two devices. represents the seawater absorption coefficient. represents the attenuation of interference caused by the shielding of the hull structure.

4. The full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 3, characterized in that: The specific content of S3 includes: S31. Calculate the interference intensity entering the front-end receiving channel of the interfered acoustic device: The interference intensity picked up by the receiving transducer of the interfered acoustic device is given by the following formula: Wherein, P represents the receiving sensitivity of the receiving transducer of the interfered acoustic device; N represents the attenuation amount of the interference source azimuth relative to the detection azimuth of the interfered acoustic device; S32. Judge whether each stage of the front-end receiving channel of the interfered acoustic device is saturated: According to the frequency response characteristics of amplification and filtering in the front-end receiving channel of the interfered acoustic device, calculate the output of each stage of the front-end receiving channel of the interfered acoustic device, and compare it with the corresponding limiting voltage. If it exceeds the limiting voltage, it means saturation occurs; if it does not exceed, it means unsaturation; Use the following formula to judge whether the frequency response characteristics of pre-amplification are saturated: Among them, is the frequency response characteristic of pre-amplification, represents the limiting voltage of pre-amplification; when ≥0, it indicates saturation occurs, and when <0, saturation does not occur; Use the following formula to judge whether pre-filtering is saturated: Among them, is the frequency response characteristic of the pre-filter, represents the limiting voltage of the pre-filter; when ≥0, it indicates saturation occurs, and when , saturation does not occur.

5. A full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 4, characterized in that: The specific content of S4 includes: According to the time domain, frequency domain and spatial domain multi-dimensional characteristics of the interference source and the interfered acoustic device, the judgment of whether the interfered acoustic device generates in-band interference is shown in the following formula: Where: M represents the background noise level during the demonstration of the indicators of the interfered acoustic device, Y represents the emission bandwidth of the interference source, g represents the acceleration due to gravity, represents the attenuation of the interference source outside the emission band compared to the interference frequency, represents the frequency-domain interference suppression ability of the interfered acoustic device, represents the spatial-domain interference suppression ability of the interfered acoustic device, represents the time-domain interference suppression ability of the interfered acoustic device; Calculating the frequency-domain interference suppression ability based on the frequency-domain processing method of the interfered acoustic device , when line spectrum processing is adopted, = 0; when broadband processing is adopted, The calculation method is shown in the following formula: Among them, represents the processing broadband of the interfered acoustic device, represents the overlapping bandwidth between the interference frequency of the interference source and the processing frequency band of the interfered acoustic device; Calculation of Spatial Interference Suppression Based on the Reception Directivity of Interfered Acoustic Devices When the interference source is in the azimuth of the detection main lobe of the interfered acoustic device, = 0; when it is in the side lobe, This is the main-to-side lobe ratio of this azimuth compared to the main lobe; Calculating time-domain interference suppression based on the overlapping duration of the interference source emission and the processing period of the interfered acoustic device and the processing integration time of the interfered acoustic device , The calculation method is as shown in the following formula: Among them, represents the overlapping duration, represents the integration duration of the interfered acoustic device.

6. The full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 5, wherein: The specific content of S5 includes: S51. Construct an interference comprehensive influence index: As shown in the following formula: Among them, and are weighting coefficients, = 0.6, = 0.4, is the calculation result of S3, and are the corresponding clipping voltage thresholds, K is the calculation result of S4, and NL is the background noise level; and are non - linear correction factors, = 1.2, = 1.

5.

7. A full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 6, characterized in that: The S5 further includes S52, which divides the interference level based on the CII value. If the CII value < 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 a moderate interference, and frequency domain notch or spatial nulling is enabled. If CII ≥ 0.7, it is a severe interference, and the device operating frequency band is forcibly switched or time domain asynchronous sampling is started.

8. A full-ship acoustic compatibility prediction method based on time-frequency-space-energy according to claim 7, characterized in that: The S5 further includes S53, which is confidence evaluation. Monte Carlo simulation is introduced to verify the reliability of the results, as follows: Among them, represents the CII mean value, which is obtained by calculating multiple times with randomly perturbed input parameters. is the standard deviation of CII; If the confidence level < 0.8, return to S51 to recalibrate the parameters.

Images