SPARSE OPTIMIZATION METHOD BASED ON CROSS-SHAPED THREE-DIMENSIONAL IMAGING SONAR ARRAY

Abstract

The present invention a sparse optimization method based on cross-shaped three-dimensional imaging sonar array, comprising the following steps: first, constructing a beam pattern simultaneously applicable to a near field and a far field based on a cross-shaped array; then, constructing an energy function required by sparse optimization according to the beam pattern; then, introducing an array element position disturbance into a simulated annealing algorithm to increase the degree of freedom of the sparse process and increase the sparse rate of the sparse array, and using the simulated annealing algorithm to sparse optimization of the energy function; finally, after optimization, a sparse optimization cross-shaped array is obtained. The present invention ensures that the three-dimensional imaging sonar system has the desired performance at any distance, and greatly reduces the hardware complexity of the system. It provides an effective method to achieve high performance and ultra-low complexity 3D imaging sonar system.

Claims

1. A sparse optimization method based on cross-shaped three-dimensional imaging sonar array, comprising the following steps: (1) constructing a beam pattern BP(W,u,v,δ,f.sub.j) simultaneously applicable to a near field and a far field based on a cross-shaped array, the beam pattern BP(W,u,v,δ,f.sub.j) being: $\mathrm{BP}\left(W,u,v,\delta ,f_j\right)=.Math.\underset{n=0}{\overset{N-1}{.Math.}}.Math.{\omega }_n.Math.\exp \left[-j.Math.\frac{2.Math.\pi .Math..Math.f_j}{c}.Math.y_n.Math.v+\delta .Math.\frac{y_n^2}{2}\right].Math.\times .Math.\underset{m=1}{\overset{M}{.Math.}}.Math.{\omega }_m.Math.\exp \left[-j.Math.\frac{2.Math.\pi .Math..Math.f_j}{c}.Math.x_m.Math.u+\delta .Math.\frac{x_m^2}{2}\right].Math.$ wherein, W is a weight coefficient of the array, including a weight coefficient ω.sub.n of the vertical transmitting array and a weight coefficient ω.sub.m of the horizontal receiving array; f.sub.j is a transmitting frequency in the vertical beam j direction; x.sub.m is a position of the m-th element of the horizontal receiving array; y.sub.n is a position of the nth element of the vertical transmitting array; c is a speed at which sound waves propagate in water;
δ=1/r−1/r.sub.0; r is a target distance; r.sub.0 is a beam focusing distance;
u=sin β.sub.a−sin θ.sub.a;
v=sin β.sub.e−sin θ.sub.e; β.sub.a is a horizontal beam arrival direction; θ.sub.a is a horizontal beam focusing direction; β.sub.e is a vertical beam arrival direction; θ.sub.e is a vertical beam focusing direction; when δ=0, BP is a far field beam pattern; when δ≠0, BP is a near field beam pattern; (2) constructing an energy function E(W,A) required by sparse optimization according to the beam pattern BP(W,u,v,δ,f.sub.j), the energy function E(W,A) being: $E\left(W,A\right)={k_1\left({\int }_{\delta .Math.\min }^{\delta .Math..Math.\max }.Math.\left(\underset{\left(u,v\right)\in \Omega }{.Math.}.Math.\left(\frac{B.Math.P\left(W,\delta ,u,v,f_j\right)}{B.Math.P_{M.Math.A.Math.X}}-b_d\right)\right).Math.d.Math.\delta \right)}^2+k_2.Math.A^2+{k_3\left(R_o-R_d\right)}^2$ wherein, k.sub.1, k.sub.2 and k.sub.3 are the weight coefficients of the corresponding items; b.sub.d is the desired beam pattern sidelobe peak; R.sub.d is the ratio of the maximum weight coefficient to the minimum weight coefficient in the weight coefficient matrix W; Ro is the ratio of the desired maximum weight coefficient to the minimum weight coefficient; the value range Ω of u and v corresponds to the part of the sidelobe beam whose intensity is greater than b.sub.d; δ.sub.min and δ.sub.max represent the minimum and maximum values of S, respectively; (3) introducing an array element position disturbance into a simulated annealing algorithm, and using the simulated annealing algorithm to sparse optimization of the energy function E(W,A); and (4) after step (3) optimization, obtaining a sparse optimization cross-shaped array.

2. The sparse optimization method based on cross-shaped three-dimensional imaging sonar array according to claim 1, wherein, the introduction of an array element position disturbance into the simulated annealing algorithm comprises: a. if the selected array element weight coefficient is not 0, that is, the selected array is in the on state, caching (ω.sub.temp) the array element weight coefficient and the current array element position p.sub.xy=(x.sub.m, y.sub.n); b. closing the array element, updating the array element weight coefficient matrix W and the number of array elements A, and calculating the energy function; c. when the energy function decreases, accepting the state and selecting the next random array element; and d. when the energy function increases, turning on the array element again, and the weight coefficient adds a random disturbance within a certain range, and a disturbance is added to the original position p.sub.xy of the array element at the same time.

3. The sparse optimization method based on cross-shaped three-dimensional imaging sonar array according to claim 2, wherein, in step (d), the formula for adding disturbance is: $\left(x_m,y_n\right)=\{\begin{array}{cc}{p}_{\mathrm{xy}}+\left(\mathrm{unifrnd}\left(-0.1\times {\lambda }_{\min },0.1\times {\lambda }_{\min }\right),0\right),& \mathrm{if}.Math..Math.\mathrm{receiving}.Math..Math.\mathrm{array}\\ p_{\mathrm{xy}}+\left(\mathrm{unifrnd}\left(-0.1\times {\lambda }_{\min },0.1\times {\lambda }_{\min }\right)\right),& \mathrm{if}.Math..Math.\mathrm{transmitting}.Math..Math.\mathrm{array}\end{array}$ wherein, unifrnd(−0.1×λ.sub.min, 0.1×λ.sub.min) is a function of the software MATLAB, which means that a random number is generated in the interval (−0.1×λ.sub.min, 0.1×λ.sub.min).

4. The sparse optimization method based on cross-shaped three-dimensional imaging sonar array according to claim 1, wherein, the calculation formula of δ.sub.min and δ.sub.max is: ${\delta }_{\min }=-\frac{2.Math.{\lambda }_{\min }}{D^2}.Math..Math.{\delta }_{\max }=\frac{2.Math.{\lambda }_{\min }}{D^2}$ wherein, λ.sub.min is the wavelength corresponding to the highest frequency among all emitted sound waves; D is the array aperture.

5. The sparse optimization method based on cross-shaped three-dimensional imaging sonar array according to claim 1, wherein, the transmitting frequency in the vertical beam j direction is 205 kHz˜300 kHz, and the step is 5 kHz.

6. The sparse optimization method based on cross-shaped three-dimensional imaging sonar array according to claim 1, wherein, in step (2), setting k.sub.1=10000, k.sub.2=20000, k.sub.3=1, b.sub.d=−22 dB.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings may be obtained from these drawings without creative work.

[0034] FIG. 1 is a flowchart of the array sparse optimization process of the present invention;

[0035] FIG. 2 is a schematic diagram of dividing the far field and near field focusing intervals of the present invention;

[0036] FIG. 3 is a schematic diagram of array disturbance in the sparse optimization process of the present invention;

[0037] FIG. 4 is a schematic diagram of the sparse cross-shaped array of the present invention;

[0038] FIG. 5 is a near field beam pattern of the obtained cross-shaped array of the present invention;

[0039] FIG. 6 is a far field beam pattern of the obtained cross-shaped array of the present invention.

DETAILED DESCRIPTION OF THE INVENTION

[0040] In order to make the objectives, technical solutions, and advantages of the present invention clearer, the following further describes the present invention in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, and do not limit the protection scope of the present invention.

[0041] In this embodiment, the initial array is a 100-element vertical transmitting array and a 100-element horizontal receiving array. The transducers are evenly distributed in a rectangular plane at half-wavelength spacing. The horizontal spacing and vertical spacing of the transducers are equal. The transmitting frequency f.sub.j is 205 kHz-300 kHz, and the step is 5 kHz, and the sound velocity is c=1500 m/s.

[0042] As shown in FIG. 1, the sparse optimization method based on cross-shaped three-dimensional imaging sonar array provided by this embodiment comprises the following steps:

[0043] S101, based on the multi-frequency transmission algorithm and the cross-shaped array, a beam pattern BP(W,u,v,δ,f.sub.j) simultaneously applicable to a near field and afar field is provided, the beam pattern BP(W,u,v,δ,f.sub.j) is:

[00003]$\mathrm{BP}\left(W,u,v,\delta ,f_j\right)=.Math.\underset{n=0}{\overset{N-1}{.Math.}}.Math.{\omega }_n.Math.\exp \left[-j.Math.\frac{2.Math.\pi .Math..Math.f_j}{c}.Math.y_n.Math.v+\delta .Math.\frac{y_n^2}{2}\right].Math.\times .Math.\underset{m=1}{\overset{M}{.Math.}}.Math.{\omega }_m.Math.\exp \left[-j.Math.\frac{2.Math.\pi .Math..Math.f_j}{c}.Math.x_m.Math.u+\delta .Math.\frac{x_m^2}{2}\right].Math.$

[0044] wherein, W is the weight coefficient of the array, including a weight coefficient ω.sub.n of the vertical transmitting array and a weight coefficient ω.sub.m of the horizontal receiving array;

[0045] f.sub.j is the transmitting frequency 205 kHz-300 kHz in the vertical beam j direction, and the step is 5 kHz, j=1, 2, . . . , 20;

[0046] x.sub.m is the position of the m-th element of the horizontal receiving array, x.sub.m=mλ.sub.min/2, m=1, 2, . . . , 100; λ.sub.min is the wavelength corresponding to the highest frequency among all emitted sound waves.

[0047] y.sub.n is the position of the n-th element of the vertical transmitting array, y.sub.n=nλ.sub.min/2, n=1, 2, . . . , 100; [0048] c is the speed at which sound waves propagate in water 1500 m/s;

δ=1/r−1/r.sub.0; [0049] r is the target distance, detecting a target within 50 meters; [0050] r.sub.0 is the beam focusing distance, the focusing range is 50 meters;

u=sin β.sub.a−sin θ.sub.a, the value range u∈[0,1];

v=sin β.sub.e−sin θ.sub.e, the value range v∈[0,1];

[0051] β.sub.a is the horizontal beam arrival direction, the viewing angle range is 60°;

[0052] θ.sub.a is the horizontal beam focusing direction, the viewing angle range is 60°;

[0053] β.sub.e is the vertical beam arrival direction, the viewing angle range is 60°;

[0054] θ.sub.e is the vertical beam focusing direction, the viewing angle range is 60°;

[0055] S102, an energy function E(W,A) required by sparse optimization according to the beam pattern BP(W,u,v,δ,f.sub.j) is constructed.

[0056] As shown in FIG. 2, the entire observation scene is divided into two parts, the near field and the far field; when δ=0, BP(W,u,v,δ,f.sub.j) is a far field beam pattern; when δ≠0, BP(W,u,v,δ,f.sub.j) is a near field beam pattern.

[0057] The near field is composed of multiple focus intervals, and each focus interval selects a beam focus distance r.sub.0, then the boundaries of the focus interval are r.sub.0− and r.sub.0+, that is, the Depth of Field (DOF) is [r.sub.0−,r.sub.0+]. r.sub.0− corresponds to δ.sub.max, r.sub.0+ corresponds to δ.sub.min.

[0058] In the depth of field, the main lobe attenuation is less than 3 dB. An energy function E(W,A) required by sparse optimization according to the beam pattern is constructed, the energy function E(W,A) is:

[00004]$E\left(W,A\right)={k_1\left({\int }_{\delta .Math.\min }^{\delta .Math..Math.\max }.Math.\left(\underset{\left(u,v\right)\in \Omega }{.Math.}.Math.\left(\frac{B.Math.P\left(W,\delta ,u,v,f_j\right)}{B.Math.P_{M.Math.A.Math.X}}-b_d\right)\right).Math.d.Math.\delta \right)}^2+k_2.Math.A^2+{k_3\left(R_o-R_d\right)}^2$

[0059] wherein, k.sub.1, k.sub.2 and k.sub.3 are the weight coefficients of the corresponding items, setting k.sub.1=10000, k.sub.2=20000, k.sub.3=1; b.sub.d is the desired beam pattern sidelobe peak, b.sub.d is −22 dB; Ro is the ratio of the maximum weight coefficient to the minimum weight coefficient in the weight coefficient matrix W; R.sub.d is the ratio of the desired maximum weight coefficient to the minimum weight coefficient, R.sub.d is 3; the value range Ω of u and v corresponds to the part of the sidelobe beam whose intensity is greater than b.sub.d; δ.sub.min and δ.sub.max represent the minimum and maximum values of δ respectively; the calculation formula of δ.sub.min and δ.sub.max is:

[00005]${\delta }_{\min }=-\frac{2.Math.{\lambda }_{\min }}{D^2}.Math..Math.{\delta }_{\max }=\frac{2.Math.{\lambda }_{\min }}{D^2}$

[0060] wherein, λ.sub.min is the wavelength corresponding to the highest frequency among all emitted sound waves, λ.sub.min is 2.5 mm; D is the array aperture, D is 25 cm.

[0061] S103, introducing an array element position disturbance into a simulated annealing algorithm, and using the simulated annealing algorithm to sparse optimization of the energy function E(W,A).

[0062] An array element position disturbance is introduced into the simulated annealing algorithm to increase the degree of freedom of the sparse process and increase the sparse rate of the sparse array. As shown in FIG. 3, the distance between the array elements is d.sub.i(i=1, 2, . . . 100), and the array elements deviate from their original positions due to position disturbances. The horizontal and vertical deviation distances are Δx and Δy respectively. Specifically, the formula for adding disturbance is:

[0063] a. if the selected array element weight coefficient is not 0, that is, the selected array is in the on state, the array element weight coefficient and the current array element position p.sub.xy=(x.sub.m, y.sub.n) are cached (ω.sub.temp).

[0064] b. closing the array element, updating the array element weight coefficient matrix W and the number of array elements A, and calculating the energy function.

[0065] c. when the energy function decreases, accepting the state and selecting the next random array element.

[0066] d. when the energy function increases, turning on the array element again, and the weight coefficient adds a random disturbance within a certain range, and a disturbance is added to the original position p.sub.xy of the array element at the same time. the formula for adding disturbance is:

[00006]$\left(x_m,y_n\right)=\{\begin{array}{cc}{p}_{\mathrm{xy}}+\left(\mathrm{unifrnd}\left(-0.1\times {\lambda }_{\min },0.1\times {\lambda }_{\min }\right),0\right),& \mathrm{if}.Math..Math.\mathrm{receiving}.Math..Math.\mathrm{array}\\ p_{\mathrm{xy}}+\left(\mathrm{unifrnd}\left(-0.1\times {\lambda }_{\min },0.1\times {\lambda }_{\min }\right)\right),& \mathrm{if}.Math..Math.\mathrm{transmitting}.Math..Math.\mathrm{array}\end{array}$

[0067] wherein, unifrnd(−0.1×λ.sub.min, 0.1×λ.sub.min) is a function of the software MATLAB, which means that a random number is generated in the interval (−0.1×λ.sub.min, 0.1×λ.sub.min).

[0068] S104, after S103 optimization, a sparse optimization cross-shaped array can be obtained, the number of array elements is 120, and the weight coefficient is between 0.5 and 1.5, as shown in FIG. 4. In addition, the beam pattern of the obtained array in the near field is shown in FIG. 5, and the sidelobe is suppressed at −22 dB; the beam pattern in the far field is shown in FIG. 6, and the sidelobe is also suppressed at −22 dB.

[0069] A cross-shaped array that satisfies the three-dimensional imaging sonar system with desired performance at any distance can be constructed by using the sparse optimization method of the cross-shaped three-dimensional imaging sonar array provided by this embodiment.

[0070] The specific implementations described above describe the technical solutions and beneficial effects of the present invention in detail. It should be understood that the above descriptions are only the most preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, additions, and equivalent replacements within the scope shall be included in the protection scope of the present invention.