Shallow sea surface underwater target distinguishing method based on correlation matrix characteristic value attenuation coefficient
Technical Field
The present invention relates to computing; the technical field of calculation or counting, in particular to a shallow sea water surface underwater target distinguishing method based on a correlation matrix eigenvalue attenuation coefficient in the field of underwater target detection.
Background
With the development of marine resources and the need of military combat, small target detection in shallow sea becomes a hot spot of underwater acoustic signal processing research. In the waveguide environment of the offshore area, the detection of the small-delay target of the propagation of the acoustic signal is a difficult problem due to the multiple reflections at the sea bottom and the sea surface. Ocean background noise level increases; in addition, in a specific sea area (such as a muddy water sea area with a high sediment content), active detection faces strong reverberation. Therefore, small target detection in shallow sea is a difficult problem.
Disclosure of Invention
The invention solves the problems in the prior art, provides an optimized shallow sea surface and underwater target resolution method based on the attenuation coefficient of the characteristic value of the related matrix, aims at the problem of autonomous resolution of a surface and underwater sound source target, and can effectively realize shallow sea surface and underwater target resolution by relying on the theoretical support of normal waves and the like and utilizing the correlation characteristics of different target vertical sound fields; the vertical submerged buoy array is used as a support, and technical support can be provided for realizing the full-coverage underwater warning capability of forming a specific sea area in a short time.
The method adopts the technical scheme that the method for distinguishing the underwater target on the shallow sea surface based on the attenuation coefficient of the characteristic value of the correlation matrix characterizes the vertical correlation of a sound field according to the size of the attenuation coefficient eta (rho) of the characteristic value of the correlation matrix of a matrix signal; the attenuation coefficient of the matrix signal correlation matrix eigenvalue on the threshold value corresponds to a water surface sound source, otherwise, the attenuation coefficient is an underwater sound source.
The invention solves the problem of positioning and/or distinguishing different targets (water surface and underwater targets) in the vertical direction by adopting the fading coefficient of the vertical array signal.
In the invention, the larger the fading coefficient of the eigenvalue of the correlation matrix of the array signal is, the more concentrated the distribution of the eigenvalue is, the weaker the corresponding oscillation of the correlation matrix is, the stronger the vertical correlation of the sound field is, and the corresponding sound source on the water surface is; the smaller the fading coefficient of the eigenvalue of the matrix signal correlation matrix, the more dispersed the eigenvalue distribution, the stronger the corresponding correlation matrix oscillation, and the weaker the vertical correlation of the sound field, corresponding to the underwater sound source.
Preferably, the method comprises the steps of:
step 1: calculating a correlation coefficient matrix rho of a shallow sea vertical array receiving target signal;
step 2: obtaining the distribution of the eigenvalue lambda (rho) of the sound field correlation coefficient matrix rho;
and step 3: processing the characteristic value lambda (rho);
and 4, step 4: obtaining a fading coefficient eta (rho) of a correlation matrix eigenvalue based on the eigenvalue, wherein the fading coefficient eta (rho) is used for representing the correlation of the sound field;
and 5: performing sound field modeling by using hydrological data of measuring environment, and calculating a fading coefficient eta of a correlation matrix eigenvalue of a vertical array copy field by using simulated sound field dataref;
Step 6: utilizing actual correlation matrix eigenvalue fading coefficient eta (rho) and simulated correlation matrix eigenvalue fading coefficient eta (rho)refAnd comparing, realizing the estimation of the sound source depth, and distinguishing the underwater target on the shallow sea surface.
Preferably, in the step 2, λ (ρ) ═ λ (λ [)1(ρ),λ2(ρ),…,λN(ρ)]And N is the number of the vertical array elements, namely the dimension of a correlation coefficient matrix.
Preferably, after the treatment of step 3, lambda1(ρ)≥λ2(ρ)≥…≥λN(ρ)。
Preferably, in the
step 4,
preferably, in said
step 6, defining
With H
0To representSound source from the surface H
1Representing a sound source from underwater; that is, when eta (rho) is less than or equal to eta
refWhen the target is a sound source from the water surface, and when η (ρ) > η
refWhen the target is a sound source from underwater.
Preferably, the critical depth is 20-30 m underwater; specifically, this is the critical depth under shallow sea conditions in the simulation of the present invention.
The invention relates to an optimized shallow sea surface underwater target resolution method based on correlation matrix eigenvalue decay coefficient, which characterizes the vertical correlation of a sound field by the size of matrix signal correlation matrix eigenvalue decay coefficient eta (rho); the fading coefficient of the matrix signal correlation matrix eigenvalue on the threshold value corresponds to a water surface sound source, otherwise, the fading coefficient is an underwater sound source; the larger the fading coefficient of the eigenvalue of the correlation matrix of the matrix signal is, the more concentrated the distribution of the eigenvalue is, the weaker the corresponding oscillation of the correlation matrix is, the stronger the vertical correlation of the sound field is, and the corresponding sound source on the water surface is; the smaller the fading coefficient of the eigenvalue of the matrix signal correlation matrix, the more dispersed the eigenvalue distribution, the stronger the corresponding correlation matrix oscillation, the weaker the vertical correlation of the sound field, and the corresponding underwater sound source.
According to the method, a dimensionless physical quantity, namely a characteristic value fading coefficient of a matrix signal correlation matrix is provided based on the influence of sound source depth on the vertical correlation of a sound field according to a normal wave theory, wherein the fading coefficient can effectively represent the vertical correlation of the sound field, and the size of the fading coefficient is closely related to the sound source depth; simulation data show that the underwater target on the surface of the shallow sea can be effectively judged by distinguishing the attenuation coefficient.
The method can be applied to deep sea vertical array anti-submarine warning detection, and has certain innovative significance.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 shows three typical hydrological conditions of a subsequent simulation in shallow sea;
FIG. 3 is a diagram of two-dimensional distribution of fading coefficients of characteristic values of correlation matrices for a matrix with sound source distances of 10km and 20km respectively in good hydrology at different depths of first matrix elements and different sound sources;
FIG. 4 is the distribution of fading coefficients of eigenvalues of matrix signal correlation matrices with sound source distances of 10km and 20km respectively under good hydrology and with a matrix first array element depth of 30m, and a black line is a water surface and underwater target resolution threshold;
FIG. 5 shows two-dimensional distribution of fading coefficients of characteristic values of correlation matrices for matrixes with sound source distances of 10km and 20km respectively in medium hydrology at different depths of initial elements and different sound sources;
FIG. 6 shows the distribution of fading coefficients of eigenvalues of matrix signal correlation matrices with sound source distances of 10km and 20km respectively under medium hydrology and with a matrix first array element depth of 30m, and a black line is a water surface underwater target resolution threshold;
FIG. 7 shows two-dimensional distribution of fading coefficients of characteristic values of correlation matrices for matrixes with sound source distances of 10km and 20km respectively in severe hydrology at different depths of initial elements and different sound sources;
fig. 8 shows the distribution of fading coefficients of eigenvalues of matrix signal correlation matrices of which the sound source distances are 10km and 20km respectively under severe hydrology and the depth of the first array element of the matrix is 30m, and a black line is a discrimination threshold of underwater targets on the water surface.
Detailed Description
The present invention is described in further detail with reference to the following examples, but the scope of the present invention is not limited thereto.
As shown in FIG. 1, the invention relates to a shallow sea water surface underwater target resolution method based on correlation matrix eigenvalue decay coefficient, which characterizes the vertical correlation of sound field by the magnitude of matrix signal correlation matrix eigenvalue decay coefficient eta (rho); the fading coefficient of the matrix signal correlation matrix eigenvalue on the threshold corresponds to the water surface sound source, otherwise, the matrix signal correlation matrix eigenvalue is an underwater sound source.
In the invention, the larger the fading coefficient of the eigenvalue of the correlation matrix of the array signal is, the more concentrated the distribution of the eigenvalue is, the weaker the corresponding oscillation of the correlation matrix is, the stronger the vertical correlation of the sound field is, and the corresponding sound source on the water surface is; the smaller the fading coefficient of the eigenvalue of the matrix signal correlation matrix, the more dispersed the eigenvalue distribution, the stronger the corresponding correlation matrix oscillation, and the weaker the vertical correlation of the sound field, corresponding to the underwater sound source.
In the actual shallow sea sound transmission process, under the condition of neglecting the side wave, the sound signal can be represented by the sum of a plurality of transmittable normal waves; the order and amplitude of the normal wave mode have close relation with the depth of a target sound source and a receiver and the like.
The method comprises the following steps:
step 1: calculating a correlation coefficient matrix rho of a shallow sea vertical array receiving target signal;
step 2: obtaining the distribution of the eigenvalue lambda (rho) of the sound field correlation coefficient matrix rho;
in the step 2, λ (ρ) ═ λ1(ρ),λ2(ρ),…,λN(ρ)]And N is the number of the vertical array elements.
And step 3: processing the characteristic value lambda (rho);
after treatment in step 3,. lambda.1(ρ)≥λ2(ρ)≥…≥λN(ρ)。
And 4, step 4: obtaining a fading coefficient eta (rho) of a correlation matrix eigenvalue based on the eigenvalue;
in the
step 4, the process of the step,
and 5: performing sound field modeling by using hydrological data of measuring environment, and calculating a fading coefficient eta of a correlation matrix eigenvalue of a vertical array copy field by using simulated sound field dataref;
In the simulation of the invention, the critical depth of shallow sea is 20-30 m underwater.
In the invention, the critical depth is usually used for simulation, namely, the vertical array is usually set to be the critical depth of 20-30 m under water during simulation; on the basis of the above, with etarefAnd (3) representing the simulation attenuation coefficient of the sound field of the vertical array, wherein the vertical array is usually set to be the critical depth of 20-30 m underwater in the simulation process, and the sound source is also the critical depth of 20-30 m.
Step 6: utilizing actual correlation matrix eigenvalue fading coefficient eta (rho) and simulated correlation matrix eigenvalue fading coefficient eta (rho)refContrast to realize the depth of sound sourceAnd (4) identifying the underwater target on the shallow sea surface.
In said
step 6, defining
With H
0Representing sound sources from the surface, H
1Representing a sound source coming from underwater.
Fig. 2 shows three typical hydrological conditions of the shallow sea for subsequent simulation, including good hydrology, medium hydrology and severe hydrology.
Correspondingly, as shown in fig. 3, for a good underwater sound source distance of 10km and 20km, the fading coefficients of the eigenvalues of the correlation matrix are distributed in two dimensions at different depths of the first array element and at different depths of the sound source; as shown in fig. 4, in good hydrology, when the depth of the first array element of the matrix is 30m, the distances of the sound sources are respectively 10km and 20km, the characteristic values of the matrix signal correlation matrix decay coefficient are distributed, and the black line is the resolution threshold of the underwater target on the water surface.
Correspondingly, as shown in fig. 5, the two-dimensional distribution of the fading coefficients of the eigenvalues of the correlation matrix is obtained for the matrixes with the sound source distances of 10km and 20km respectively in the medium hydrology at different first array element depths and different sound source depths; as shown in fig. 6, under the condition that the depth of the first array element of the matrix is 30m in medium hydrology, the distances of the sound sources are respectively 10km and 20km, the characteristic values of the correlation matrix of the matrix signals are distributed with fading coefficients, and the black line is the underwater target resolution threshold of the water surface.
Correspondingly, as shown in fig. 7, the two-dimensional distribution of the fading coefficients of the eigenvalues of the correlation matrix at different depths of the first array element and at different depths of the sound source is a matrix with sound source distances of 10km and 20km respectively in the severe hydrology; as shown in fig. 8, when the depth of the first array element of the matrix is 30m in severe water, the distances of the sound sources are 10km and 20km respectively, the characteristic values of the matrix signal correlation matrix decay coefficient are distributed, and the black line is the resolution threshold of the underwater target on the water surface.