CN110763327A - Single-vector hydrophone signal detection method under condition of plane baffle - Google Patents
Single-vector hydrophone signal detection method under condition of plane baffle Download PDFInfo
- Publication number
- CN110763327A CN110763327A CN201910938268.2A CN201910938268A CN110763327A CN 110763327 A CN110763327 A CN 110763327A CN 201910938268 A CN201910938268 A CN 201910938268A CN 110763327 A CN110763327 A CN 110763327A
- Authority
- CN
- China
- Prior art keywords
- vector hydrophone
- signal
- received
- vibration velocity
- under
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to the technical field of underwater sound, and discloses a single-vector hydrophone signal detection method under the condition of a plane baffle, which comprises the following steps: step (1): establishing a transmission matrix of the planar cavity baffle, and obtaining a reflection coefficient R of the planar cavity baffle according to the transmission matrix and boundary conditions; step (2): establishing a single-frequency signal model received by the vector hydrophone and a broadband signal model received by the vector hydrophone under the condition of a plane cavity baffle according to the reflection coefficient R; and (3): solving a single-frequency signal received by the vector hydrophone under the isotropic noise interference background and a broadband signal received by the vector hydrophone under the isotropic noise interference background; and (4): converting the signals obtained in the step (3) to a frequency domain to obtain corresponding spectrums, averaging after sound pressure vibration velocity cross spectrums to obtain complex sound intensity output of an average periodogram, and then detecting the signals by using a complex sound intensity device; and (5): and evaluating the performance of the complex sound intensity device real part detection and the modulus detection by using an ROC curve.
Description
Technical Field
The invention belongs to the technical field of underwater sound, and particularly relates to a single-vector hydrophone signal detection method under the condition of a plane baffle.
Background
Water acoustics is used as a branch of acoustics, and the theories of generation, radiation, propagation and reception of sound waves underwater are mainly studied to solve the acoustics problems related to underwater target detection, identification and information transmission processes. In sea warfare, sonar is the five sense organs of individual sea warfare, and all underwater battlefield reconnaissance needs sonar as a medium, but the sonar is not necessary. The underwater acoustic transducer is one of important parts of a sonar system and is an important research direction of underwater acoustics, and the research of a novel underwater acoustic transducer is a key content of the development of naval sonar technology. The underwater acoustic transducer is a general term for various underwater transmitting and receiving measuring sensors, converts underwater acoustic signals into electric signals, or converts electric signals into underwater acoustic signals, and is an important component of sonar. The performance of a sonar is directly related to the performance of an underwater acoustic transducer. In underwater acoustic engineering, transducer technology is at a fundamental position. The progress of transducer technology can drive the improvement of the technical level of sonar systems, so that the research work of novel underwater acoustic transducers has important significance.
Underwater target detection is an important content of water acoustics, is a premise of target parameter estimation and positioning, and is one of main functions of a sonar system. In more than half a century, from the initial analog amplitude limiting related detection method to the combination of digital signal processing methods such as DEMON and LOFAR spectrum analysis to detect the target signal, passive sonar signal detection has achieved many excellent research results, and new vitality is injected into the passive sonar detection technology.
In passive sonar, target detection using the line spectrum or broadband continuum spectrum of the target radiation is an important means of finding the target. As early as the twenty-first century, the students of Junying and Liuhong have proposed a simple and effective signal detector, namely an average intensity detector. The monograph "huijing" vector sound signal processing foundation "national defense industry press, 2009" indicates that the average sound intensity device has good anti-interference capability because the sound pressure and the vibration velocity of the target signal in the marine waveguide are correlated, and the sound pressure and the vibration velocity of the isotropic environmental interference are uncorrelated or correlated very weakly. Averaging the intensity of the sound is very effective when detecting a single target of radiation continuum noise. However, when multiple targets exist, the average intensity speaker can only measure the direction of the synthesized sound intensity flow of the multiple targets, and cannot distinguish the directions of the multiple targets.
In the nineties of the last century, by means of technical introduction and international cooperative research, the underwater acoustic vector signal sensor, also called a vector hydrophone, is successfully developed in China, and can output sound pressure and vibration velocity signals at the same point, so that signal detection under the condition of a baffle can be realized only by means of a single vector hydrophone. However, in practical engineering applications, the sonar basic array is provided with an acoustic baffle, the sound field is changed due to the existence of the baffle, the above traditional underwater acoustic signal detection methods are all based on the free field assumption, and relevant documents of the signal complex sound intensity detection method under the non-free field condition are not reported yet.
Disclosure of Invention
The invention aims to disclose a single-vector hydrophone signal detection method under the condition of a plane baffle, which has strong anti-interference capability and strong adaptability.
The purpose of the invention is realized as follows:
a single-vector hydrophone signal detection method under the condition of a plane baffle comprises the following steps:
step (1): establishing a transmission matrix of a planar cavity baffle with a thin steel plate-air-thin steel plate three-layer structure, and obtaining a reflection coefficient R of the planar cavity baffle according to the transmission matrix and boundary conditions between layers of the planar cavity baffle;
step (2): establishing a single-frequency signal model received by the vector hydrophone and a broadband signal model received by the vector hydrophone under the condition of a plane cavity baffle according to the reflection coefficient R, and obtaining a single-frequency signal v received by a vibration velocity x channel of the vector hydrophone according to the single-frequency signal model received by the vector hydrophonexSingle-frequency signal v received by y channel of sum vector hydrophone vibration velocityyWideband signal model received by vector hydrophoneObtaining a vibration velocity signal S received by a vibration velocity x channel of the vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of vector hydrophonevy；
And (3): carrying out single-frequency signal receiving by the vector hydrophone, broadband signal receiving by the vector hydrophone and single-frequency signal v receiving by the vibration speed x channel of the vector hydrophone in the step (2)xSingle-frequency signal v received by y channel of vibration velocity of vector hydrophoneyVibration velocity signal S received by vibration velocity x channel of vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of sum vector hydrophonevyAdding isotropic noise to obtain a single-frequency signal received by the vector hydrophone under the isotropic noise interference background and a broadband signal received by the vector hydrophone under the isotropic noise interference background;
and (4): converting the single-frequency signal received by the vector hydrophone under the isotropic noise interference background and the broadband signal received by the vector hydrophone under the isotropic noise interference background obtained in the step (3) to a frequency domain to obtain corresponding spectrums, averaging after mutual spectrum of sound pressure vibration velocity to obtain complex sound intensity output of an average periodogram, and then respectively adding real sound intensity parts in x and y directions and adding a module value by using a complex sound intensity device to detect signals;
and (5): and evaluating the performance of the complex sound intensity device real part detection and the modulus detection by using an ROC curve.
Further, the step (2) comprises:
step (2.1): according to the reflection coefficient R, establishing a single-frequency signal model received by the vector hydrophone under the condition of a plane cavity baffle:
in the above formula, kx、kyIs the wavenumber, ω is the angular frequency, t is time;
step (2.2): obtaining single-frequency signals received by x and y channels of the vibration speed of the vector hydrophone according to the single-frequency signal model received by the vector hydrophone:
in the above formula, ρ is the density of the propagation medium; v. ofxIs a single-frequency signal received by a vector hydrophone vibration velocity x channelyThe signal is a single-frequency signal received by a vibration velocity y channel of the vector hydrophone;
step (2.3): establishing a broadband signal model received by the vector hydrophone:
s(t)＝si(t)+sr(t)；
in the above formula, si(t) is a band-limited white plane wave noise broadband incident sound pressure signal; sr(t) is si(t) signals reflected by the planar cavity baffles, sr(t)＝R·si(t- τ); tau is the time delay difference existing between the incident sound pressure signal and the reflected sound pressure;
step (2.4): obtaining a vibration velocity signal S received by a vibration velocity x channel of the vector hydrophone according to a broadband signal model received by the vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of vector hydrophonevy：
Further, the single-frequency signal received by the vector hydrophone under the isotropic noise interference background in step (3) is: (ii) a
In the above formula, np(t) is the acoustic pressure signal noise received by the vector hydrophone, nvx(t) is the vector hydrophone vibration velocity x channel noise, nvy(t) is the vector hydrophone vibration velocity y-channel noise;
the broadband signal received by the vector hydrophone under the isotropic noise interference background is as follows:
the invention has the beneficial effects that:
the detection method has good anti-interference capability, can detect signals only by a single hydrophone, has a simple structure, reduces the equipment volume and improves the ship suitability.
Drawings
FIG. 1 is a method for detecting a single-vector hydrophone signal under a planar mask condition;
FIG. 2 is a model of a vector hydrophone received signal under shadow mask conditions;
FIG. 3 is a block diagram of a complex sound intensifier;
FIG. 4 is [ H ]0]、[H1]Assuming a probability density curve of a lower random process x;
FIG. 5 is a ROC curve for complex sound intensity detection of single-frequency signals under the condition of a planar baffle with a signal-to-noise ratio of-26 dB;
FIG. 6 is a ROC curve of the average sound intensity detection of broadband signals under the condition of a planar baffle with a signal-to-noise ratio of-18 dB;
FIG. 7 is a complex sound intensity detection ROC curve of broadband signals under the condition of a planar baffle at a signal-to-noise ratio of-18 dB.
Detailed Description
The invention is further described below with reference to the accompanying drawings:
a single-vector hydrophone signal detection method under the condition of a plane baffle comprises the following steps:
step (1): establishing a transmission matrix of a planar cavity baffle plate with a thin steel plate-air-thin steel plate three-layer structure, and obtaining a reflection coefficient R of the baffle plate according to the transmission matrix and boundary conditions between layers of the cavity baffle plate;
the baffle plate is a thin steel plate-air-thin steel plate 3-layer cavity structure, and as boundary conditions between layers are that normal particle vibration speed is continuous and sound pressure is continuous, sound pressure and vibration speed of the upper and lower boundaries of the planar cavity baffle plate can be related through a transfer matrix to be expressed as a formula. Wherein the transfer matrix is determined by the structure and parameters of the planar cavity baffle itself. And solving the formula to obtain the complex reflection coefficient R of the baffle.
Step (2): obtaining reflected sound pressure of a single-frequency incident sound pressure signal and a broadband incident sound pressure signal according to the reflection coefficient R, and superposing the incident sound pressure and the reflected sound pressure to obtain a sound pressure channel receiving signal model under the condition of a plane baffle, so as to obtain a vibration velocity signal model received by a vector hydrophone vibration velocity x and y channel;
specifically, according to the baffle reflection coefficient R, obtaining a single-frequency sound pressure receiving signal model of the single-vector hydrophone under the baffle condition:
wherein k isx、kyIs the wave number, if the signal incidence direction is theta, kx＝kcosθ，ky＝ksinθ。
Therefore, single-frequency signals received by x and y channels of the vector hydrophone vibration speed can be obtained:
where ρ is the density of the propagation medium (water).
It can be seen from the above that there is a fixed phase difference e between the reflected sound pressure and the incident sound pressure-j2kydFrom this, it can be expressed that the delay difference τ between the incident sound pressure signal and the reflected sound pressure signal of the baffle is:
if the band-limited white plane wave noise broadband incident sound pressure signal is si(t) then si(t) signal s reflected by the baffler(t)＝R·si(t- τ), the broadband sound pressure signal received by the vector hydrophone is:
s(t)＝si(t)+sr(t)；
similarly, the broadband vibration velocity signal received by the single-vector hydrophone under the condition of the plane baffle is as follows:
and (3): adding isotropic noise to the sound pressure signal and the vibration velocity signals in the x and y directions in the step (2);
according to the model that a single-vector hydrophone receives single-frequency and broadband signals under the condition of the baffle, the single-frequency signals received by the vector hydrophone under the isotropic noise interference background can be represented as follows:
the broadband signal received by the vector hydrophone under the isotropic noise interference background is as follows:
and (4): converting the single-vector hydrophone non-free field received sound pressure signals and x and y direction vibration velocity signals obtained in the step (3) under the isotropic noise background into a frequency domain to obtain corresponding spectrums, averaging after sound pressure vibration velocity cross-spectrums to obtain complex sound intensity output of an average periodogram, and then adding real parts of the x and y direction sound intensity and adding a mode value respectively by using a complex sound intensity device to carry out signal detection;
firstly, FFT conversion is carried out on the sound pressure signal and the vibration velocity signal in the x and y directions corresponding to the sound pressure signal, and a corresponding spectrum is obtained. The same way yields a spectrum of isotropic noise. And then averaging the spectrum by a sliding time window to obtain an average periodogram output. When the complex sound intensity real part is used for detection, the complex sound intensity real parts in the x and y directions are added and output; when the complex sound intensity device is used for conducting the module taking detection, the complex sound intensities in the x direction and the y direction are subjected to the module taking addition output.
And (5): and evaluating the performance of the complex sound intensity device real part detection and the modulus detection by using an ROC curve.
The output of the signal containing noise through the sound intensity device is larger than the threshold yTIs the detection probability P of the target signald(ii) a The output of the ambient noise passing through the sound intensity device is larger than the threshold yTThe probability of occurrence of an event in the set is a false alarm probability Pf. Can utilize the detection probability PdAnd false alarm probability PfConfidence limits are described to evaluate the performance of the detection methods proposed in the present invention.
With reference to fig. 1, a specific embodiment of the present invention is as follows:
first, the reflection coefficient R of the plane shadow mask is calculated as follows:
the flat shadow mask shown in FIG. 2 is simplified to a sheet-air-sheet structure with a thickness of 0.3cm for the upper and lower sheets and a thickness of 3cm for the air layer. And establishing a mathematical model according to the baffle structure and the baffle parameters, solving a transmission matrix of the sound wave incident baffle, solving the transmission matrix according to boundary conditions between layers, and solving a reflection coefficient R of the sound wave incident baffle. The reflection coefficient of the planar cavity mask is similar to that of a soft mask, and in order to highlight the key problem, only R-1 is taken for analysis.
It must be noted here that the process of solving is not illustrated in the present invention. The unexplained parts are well known to those skilled in the art, and the literature on the relevant records can be found by those skilled in the art according to the names or functions of the present invention, and thus are not further described.
Secondly, obtaining a single-frequency and broadband signal model received by the single-vector hydrophone under the condition of the baffle:
for a single-frequency signal, considering the two-dimensional case, taking the vector hydrophone coordinates as (x, y), without loss of generality, the incident sound pressure signal can be expressed as:
in the formulaθ is the plane wave incident angle, k ═ ω/c is the wave number, c is the speed of sound wave propagating in water, ω is the angular frequency, t is the time of sound wave propagation, and a is the sound pressure amplitude.
The reflected sound pressure received by the vector hydrophone is:
wherein R is the complex reflection coefficient of the mask, and it can be seen from the first step that R in the present invention is equal to-1.
Under the premise of not considering propagation loss and other noise interference, the sound pressure signal received by a single vector hydrophone is as follows:
p(t)＝pi(t)+pr(t)(3)；
the signal received by the x and y channels of the vibration velocity of the vector hydrophone can be obtained:
where ρ is the density of the propagation medium (water).
For a broadband signal, obtaining a band-limited white noise signal s satisfying Gaussian distribution through an I-type Chebyshev filteri(t) of (d). From the formulas (1) and (2), it is known that:
and d is the spacing between the vector hydrophone and the baffle. The equation (5) shows that there is a fixed phase difference between the reflected sound pressure and the incident sound pressure
Accordingly, the time delay difference tau existing between the incident sound pressure signal and the reflected sound pressure signal of the baffle can be represented:
therefore si(t) signal s reflected by the baffler(t)＝R·si(t- τ), the broadband sound pressure signal received by the vector hydrophone is:
s(t)＝si(t)+sr(t) (7)
similarly, the broadband vibration velocity signal received by the single-vector hydrophone under the condition of the plane baffle is as follows:
thirdly, adding isotropic noise to the sound pressure signal and the vibration velocity signal in the x and y directions under the condition of the non-free field obtained in the second step:
for isotropic noise fields, the spectral density distribution is independent of direction and uniformly distributed and non-directional, i.e. plane waves from all directions in space (which can be regarded as a sphere with a large radius) have the same power, so that the isotropic noise is generally gaussian noise. Setting a SNR value to obtain the noise signal n in the present inventionp(t)、nvx(t)、nvy(t)。
Considering noise interference, the single-frequency signal received by the vector hydrophone is:
the vector hydrophone receives a broadband signal of:
fourthly, respectively taking the sound intensity real part and the modulus value by using a complex sound intensity device to carry out signal detection:
a block diagram of a complex sound intensifier is shown in fig. 3.
For single frequency signals, first, P (t) and V are comparedi(t) (i ═ x, y) is FFT transformed to obtain corresponding spectra P (ω) and Vi(ω), the sound pressure vibration velocity cross spectrum is:
in the formula: symbol denotes a conjugate operation.
The FFT is used for replacing the Fourier transform in operation.
Using sliding time window averaging, for Ipvi(ω) average, resulting in an average periodogram output of:
in the formula: < > represents a sliding window average periodogram.
Summing the sound impositions in x, y directions obtained in formula (12) to obtain the total sound intensity Ipv(ω). It should be noted that I is the actual value of the actual valuepvx(omega) and Ipvy(ω) addition of real parts; when taking the model to examine Ipvx(omega) and IpvyAnd (omega) adding the modulus values. Corresponding to two detection methods of real part and modulus value of the complex sound intensity device in fig. 3.
In the same way, the complex sound intensifier detection output of the broadband signal and the isotropic noise can be obtained.
And fifthly, evaluating the performance of real part detection and modulus detection of the complex sound intensity device by using an ROC curve:
the ROC curve, i.e., the receiver operating characteristic curve, can describe the detection performance of the detector under a given signal-to-noise ratio condition. Setting the threshold of the detector to yTFig. 4 shows a probability density curve of noise and a probability density function curve of a signal containing noise.
The thin shaded portion of FIG. 4 represents the output of the noise-containing signal through the beamformer being greater than the threshold yTIs the detection probability P of the target signald(ii) a FIG. 4 shows the thick slope portion when the output of the ambient noise passing through the sound intensity device is greater than the threshold yTThe probability of occurrence of an event in the set is a false alarm probability Pf. Can utilize the detection probability PdAnd false alarm probability PfConfidence limits are described.
The signal detection process can be regarded as a random process, denoted by x, with a detection probability P for the detected sound intensity x, i.e. the output of the sound intensity generatordAnd false alarm probability PfCan be represented by the following formula:
in the formula: [ H ]0]Indicates only interference, [ H ]1]Representing the interference plus target signal; f0(x) Is [ H ]0]A probability density function of x under the assumed conditions; f1(x) Is [ H ]1]A conditional probability density function of x is assumed.
The above is a description of the specific embodiments of the present invention, and the following is a further description of the specific embodiments of the present invention by using simulation examples.
Simulation example 1: detecting single-frequency signals under the condition of a plane baffle:
the simulation parameters are set as follows: the plane wave line spectrum frequency f is 1000Hz, the isotropic noise background, the sampling interval T is 0.0002s, the integration time T is 2s, and the signal-to-noise ratio SNR of the detection signal is-26 dB. Wherein the parameters of the baffle are as follows: density of steel sheet rho17800kg/m3Velocity of longitudinal wave cl5900 m/s; air density ρ2Is 1.29kg/m3The sound velocity c is 1500m/s, and the distance d between the vector hydrophone and the baffle is 0.1 m. 1000 Monte Carlo simulation experiments were performed. Fig. 5 shows the detection results of two different detection modes of complex sound intensity, and it can be seen that the effect of cross-spectrum output mode-taking detection of single-vector hydrophone detection is better than that of cross-spectrum output real-taking part detection.
Simulation example 2: detecting a broadband signal average sound intensity device under the condition of a plane baffle:
the simulation parameters are set as follows: the incident sound wave is band-limited white plane wave noise, the frequency range f is 500 Hz-1300 Hz, the sampling interval T is 0.0002s, the integration time T is 2s, and the signal-to-noise ratio SNR of the detection signal is-18 dB. Wherein the parameters of the baffle are as follows: density of steel sheet rho17800kg/m3Velocity of longitudinal wave cl5900 m/s; air density ρ2Is 1.29kg/m3The sound velocity c is 1500m/s, and the distance d between the vector hydrophone and the baffle is 0.1 m. 1000 Monte Carlo simulation experiments were performed. Fig. 6 shows the detection result of the average intensity detector.
Simulation example 3: detecting a broadband signal complex intensity detector under the condition of a plane baffle:
the simulation parameter settings were the same as those of simulation example 2. Fig. 7 shows the detection results of two different detection modes of complex sound intensity, which can show that the detection of broadband signals under the condition of a planar mask, the detection performance of the modulus taking part is obviously better than that of the real part, and the detection has the same characteristics as the detection of narrowband signals under the condition of the mask. Comparing with the detection result of the average sound intensity device in fig. 6, it can be seen that the complex sound intensity device has better detection performance under the condition of the same signal to noise ratio.
In the ocean channel, acoustic ohm's law is approximately satisfied, and the sound pressure and the vibration velocity are in phase. According to the basic characteristic of Fourier transform, the energy of two same-phase inputs is concentrated on the real part of a cross spectrum, so that the energy of a target signal under the condition of a free field is concentrated on the real part of a cross spectrum output of a complex intensity amplifier, and the imaginary part is mainly interference energy. The signal obtained by reflection from the baffle has one incident signal
If the phase difference of (2) is taken as a real part only, a part of the target signal information is lost. Therefore, for complex sound intensity detection under the condition of a plane baffle, a real part is preferably not simply taken, but a mode of cross-spectrum output is taken, so that the energy of a target signal can be reserved, and the target signal can be better detected.
The analysis result of the simulation example shows that: (1) the invention relates to a signal detection method based on a free field, which realizes the detection of signals under the condition of a non-free field by constructing single-frequency and broadband signal models under the condition of a plane baffle, performing cross-spectrum processing on corresponding sound pressure and vibration velocity and using two complex sound intensity detection methods of real part detection and modulus detection. (2) Although the operation amount of the complex sound intensity detector is slightly larger than that of the average sound intensity device, the detection performance of the complex sound intensity detector is obviously superior to that of the average sound intensity device. (3) The invention popularizes the signal detection method under the condition of the free field to the non-free field, and theoretically analyzes the detection probability of the complex sound intensifier modulus detection higher than the detection probability of the real part. The method can realize signal detection with high detection probability under the condition of non-free field only by using a small-sized vector hydrophone, is simple and effective, has strong anti-interference capability and huge application potential in engineering, and can be widely applied to the technical field of underwater sound signal detection.
In conclusion, the invention provides a single-vector hydrophone signal detection method under the condition of a plane baffle. Firstly, establishing a transfer matrix according to the structure and parameters of a baffle, and solving the reflection coefficient R of a plane baffle according to the transfer matrix and the boundary conditions between layers of the plane cavity baffle; determining a single-frequency and broadband signal model received by the vector hydrophone under the condition of the baffle by using the reflection coefficient R; adding isotropic noise to the two obtained signal models; the corresponding sound pressure and vibration velocity signals are processed with cross spectrum to respectively obtain the complex sound intensity of the signals in the x and y directions, the detection method used in the invention is the module taking detection of the complex sound intensity device, and the output of the complex sound intensity device is the addition of the module values of the complex sound intensity in the x and y directions during the module taking detection; and finally, evaluating the detection performance of the complex sound intensity device modulus detection method by using an ROC curve. Because a fixed phase difference exists between the incident sound signal and the reflected sound signal, the complex sound intensity modular detection can fully utilize signal information under the condition of the baffle, and the cross-spectrum output only ignores the imaginary part and loses part of the signal information, so that the complex sound intensity modular detection has higher detection gain. The conventional signal detection method is based on the free field hypothesis, the signal detection method provided by the invention can be used for detecting signals under the non-free field condition, is simple and effective, has strong anti-interference capability, and can be widely applied to the technical field of underwater sound signal detection.
Claims (3)
1. A single-vector hydrophone signal detection method under the condition of a plane baffle is characterized by comprising the following steps: comprises the following steps:
step (1): establishing a transmission matrix of a planar cavity baffle with a thin steel plate-air-thin steel plate three-layer structure, and obtaining a reflection coefficient R of the planar cavity baffle according to the transmission matrix and boundary conditions between layers of the planar cavity baffle;
step (2): establishing a single-frequency signal model received by the vector hydrophone and a wide-band signal received by the vector hydrophone under the condition of a plane cavity baffle according to the reflection coefficient RA signal model, which is used for obtaining a single-frequency signal v received by a vibration velocity x channel of the vector hydrophone from the single-frequency signal model received by the vector hydrophonexSingle-frequency signal v received by y channel of sum vector hydrophone vibration velocityyObtaining a vibration velocity signal S received by a vibration velocity x channel of the vector hydrophone through a broadband signal model received by the vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of vector hydrophonevy；
And (3): carrying out single-frequency signal receiving by the vector hydrophone, broadband signal receiving by the vector hydrophone and single-frequency signal v receiving by the vibration speed x channel of the vector hydrophone in the step (2)xSingle-frequency signal v received by y channel of vibration velocity of vector hydrophoneyVibration velocity signal S received by vibration velocity x channel of vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of sum vector hydrophonevyAdding isotropic noise to obtain a single-frequency signal received by the vector hydrophone under the isotropic noise interference background and a broadband signal received by the vector hydrophone under the isotropic noise interference background;
and (4): converting the single-frequency signal received by the vector hydrophone under the isotropic noise interference background and the broadband signal received by the vector hydrophone under the isotropic noise interference background obtained in the step (3) to a frequency domain to obtain corresponding spectrums, averaging after mutual spectrum of sound pressure vibration velocity to obtain complex sound intensity output of an average periodogram, and then respectively adding real sound intensity parts in x and y directions and adding a module value by using a complex sound intensity device to detect signals;
and (5): and evaluating the performance of the complex sound intensity device real part detection and the modulus detection by using an ROC curve.
2. The method for detecting the signal of the single-vector hydrophone under the condition of the planar mask as claimed in claim 1, wherein the method comprises the following steps: the step (2) comprises the following steps:
step (2.1): according to the reflection coefficient R, establishing a single-frequency signal model received by the vector hydrophone under the condition of a plane cavity baffle:
in the above formula, kx、kyIs the wavenumber, ω is the angular frequency, t is time;
step (2.2): obtaining single-frequency signals received by x and y channels of the vibration speed of the vector hydrophone according to the single-frequency signal model received by the vector hydrophone:
in the above formula, ρ is the density of the propagation medium; v. ofxIs a single-frequency signal received by a vector hydrophone vibration velocity x channelyThe signal is a single-frequency signal received by a vibration velocity y channel of the vector hydrophone;
step (2.3): establishing a broadband signal model received by the vector hydrophone:
s(t)＝si(t)+sr(t)；
in the above formula, si(t) is a band-limited white plane wave noise broadband incident sound pressure signal; sr(t) is si(t) signals reflected by the planar cavity baffles, sr(t)＝R·si(t- τ); tau is the time delay difference existing between the incident sound pressure signal and the reflected sound pressure;
step (2.4): obtaining a vibration velocity signal S received by a vibration velocity x channel of the vector hydrophone according to a broadband signal model received by the vector hydrophonevxVibration velocity signal S received by vibration velocity y channel of vector hydrophonevy：
3. The method for detecting the signal of the single-vector hydrophone under the condition of the planar mask as claimed in claim 2, wherein the method comprises the following steps: the single-frequency signal received by the vector hydrophone under the isotropic noise interference background in the step (3) is:
in the above formula, np(t) is the acoustic pressure signal noise received by the vector hydrophone, nvx(t) is the vector hydrophone vibration velocity x channel noise, nvy(t) is the vector hydrophone vibration velocity y-channel noise;
the broadband signal received by the vector hydrophone under the isotropic noise interference background is as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910938268.2A CN110763327A (en) | 2019-09-30 | 2019-09-30 | Single-vector hydrophone signal detection method under condition of plane baffle |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910938268.2A CN110763327A (en) | 2019-09-30 | 2019-09-30 | Single-vector hydrophone signal detection method under condition of plane baffle |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110763327A true CN110763327A (en) | 2020-02-07 |
Family
ID=69330274
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910938268.2A Pending CN110763327A (en) | 2019-09-30 | 2019-09-30 | Single-vector hydrophone signal detection method under condition of plane baffle |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110763327A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111273265A (en) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | Single-vector hydrophone autonomous ranging device and method based on DEMON spectrum |
| CN111580078A (en) * | 2020-04-14 | 2020-08-25 | 哈尔滨工程大学 | Single hydrophone target identification method based on fusion mode flicker index |
| CN112099019A (en) * | 2020-07-16 | 2020-12-18 | 中国海洋石油集团有限公司 | Underwater sound vector detection device |
| CN112964349A (en) * | 2020-12-03 | 2021-06-15 | 中国船舶重工集团公司第七一五研究所 | Noise measurement hydrophone sensitivity absolute calibration device and method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2312573A2 (en) * | 2009-10-14 | 2011-04-20 | Thales Holdings UK Plc | Electronic baffling of sensor arrays |
| CN102353937A (en) * | 2011-07-06 | 2012-02-15 | 哈尔滨工程大学 | Single-vector active acoustic intensity averager |
| CN109001297A (en) * | 2018-06-05 | 2018-12-14 | 哈尔滨工程大学 | Large sample underwater acoustic materials acoustical reflection factor measurement method based on single vector hydrophone |
-
2019
- 2019-09-30 CN CN201910938268.2A patent/CN110763327A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2312573A2 (en) * | 2009-10-14 | 2011-04-20 | Thales Holdings UK Plc | Electronic baffling of sensor arrays |
| CN102353937A (en) * | 2011-07-06 | 2012-02-15 | 哈尔滨工程大学 | Single-vector active acoustic intensity averager |
| CN109001297A (en) * | 2018-06-05 | 2018-12-14 | 哈尔滨工程大学 | Large sample underwater acoustic materials acoustical reflection factor measurement method based on single vector hydrophone |
Non-Patent Citations (1)
| Title |
|---|
| 秦松磊: "矢量舷侧阵信号检测和目标方位估计研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111273265A (en) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | Single-vector hydrophone autonomous ranging device and method based on DEMON spectrum |
| CN111273265B (en) * | 2020-03-02 | 2020-09-29 | 中国人民解放军国防科技大学 | Single-vector hydrophone autonomous ranging device and method based on DEMON spectrum |
| CN111580078A (en) * | 2020-04-14 | 2020-08-25 | 哈尔滨工程大学 | Single hydrophone target identification method based on fusion mode flicker index |
| CN111580078B (en) * | 2020-04-14 | 2022-09-09 | 哈尔滨工程大学 | Single hydrophone target identification method based on fusion modal flicker index |
| CN112099019A (en) * | 2020-07-16 | 2020-12-18 | 中国海洋石油集团有限公司 | Underwater sound vector detection device |
| CN112964349A (en) * | 2020-12-03 | 2021-06-15 | 中国船舶重工集团公司第七一五研究所 | Noise measurement hydrophone sensitivity absolute calibration device and method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110763327A (en) | Single-vector hydrophone signal detection method under condition of plane baffle | |
| CN105301114B (en) | A kind of acoustic stimulation insertion loss measuring method based on multichannel space-time liftering technology | |
| CN109541548B (en) | Air sonar positioning method based on matching field | |
| CN102645265A (en) | Ship radiated noise level measuring method based on virtual time reversal mirror | |
| CN113884986B (en) | Beam focusing enhanced strong impact signal space-time domain joint detection method and system | |
| CN111024209A (en) | Line spectrum detection method suitable for vector hydrophone | |
| CN105301580A (en) | Passive detection method based on split array cross-spectrum phase difference variance weighing | |
| Ma et al. | Underwater broadband source detection using an acoustic vector sensor with an adaptive passive matched filter | |
| CN109001297B (en) | Method for measuring acoustic reflection coefficient of large-sample underwater acoustic material based on single-vector hydrophone | |
| CN112987003A (en) | HFM signal separation method and system in active sonar | |
| CN113281727A (en) | Output enhanced beam forming method and system based on horizontal line array | |
| CN104180891A (en) | Method for measuring ocean sound transmission loss in real time based on sound matrix | |
| Liu et al. | High-resolution and low-sidelobe forward-look sonar imaging using deconvolution | |
| Zou et al. | Vector hydrophone array development and its associated DOA estimation algorithms | |
| CN110361690B (en) | Single-vector hydrophone direction finding method under condition of plane baffle | |
| Chen et al. | Insight into split beam cross-correlator detector with the prewhitening technique | |
| CN116106879A (en) | Linear array line spectrum coherent accumulation detection method in multi-path environment | |
| CN113126029B (en) | Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment | |
| CN112684437B (en) | Passive ranging method based on time domain warping transformation | |
| CN113704685B (en) | Deep sea blind deconvolution method based on vertical line array | |
| CN112834623B (en) | Acoustic covering layer echo reduction measurement method based on compressed sensing technology | |
| CN117849882B (en) | Submerged oil detection method based on sound vortex technology | |
| CN115586493A (en) | Noise positioning method and system based on GCC-PATH method | |
| Waters et al. | Isolating scattering resonances of an air-filled spherical shell using iterative, single-channel time reversal | |
| Song et al. | Tow-ship interference suppression based on blind source separation for passive sonar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200207 |
|
| RJ01 | Rejection of invention patent application after publication |