METHOD FOR DETECTING MOVING TARGET BASED ON SPATIAL SLICES OF TRANSFORMED SPATIO-TEMPORAL FREQUENCY SPACE

Abstract

The present disclosure belongs to the field of underwater target detection, and in particular, to a method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space. The method includes: segmenting a target radiated acoustic signal received by an M-element horizontal line array in an underwater acoustic environment with a low signal-to-noise ratio (SNR); performing N-point discrete Fourier transform (DFT) on the received signal on each array element in each period of time; performing frequency domain beamforming on an array signal after each section of DFT, and performing stacking after compensating a phase difference between arrays brought by an azimuth of each primitive element; performing coordinate transformation on a frequency-azimuth-time three-dimensional (3D) matrix space obtained; taking a slice from the obtained frequency-azimuth-time 3D space subjected to the coordinate transformation; and performing segmented Radon transform on the spatial slice obtained to detect the target.

Claims

1. A method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, comprising the following steps: step 1: segmenting a target radiated acoustic signal s(t) received by an M-element horizontal line array in an underwater acoustic environment with a low signal-to-noise ratio (SNR); step 2: performing N-point discrete Fourier transform (DFT) on the received signal on each array element in each period of time τ.sub.p in step 1, wherein N=T.sub.0.Math.f.sub.s; step 3: performing frequency domain beamforming on an array signal after each section of DFT in step 2, and performing stacking after compensating a phase difference between arrays brought by an azimuth α(τ.sub.p) of each primitive element; step 4: performing coordinate transformation on a frequency-azimuth-time (f-α-t) three-dimensional (3D) matrix space obtained in step 3; step 5: taking a slice from the frequency-azimuth-time (f-cos θ-t) 3D space subjected to the coordinate transformation obtained in step 4; and step 6: performing segmented Radon transform on the spatial slice obtained in step 5 to detect the target.

2. The method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space according to claim 1, wherein step 1 is specifically as follows:
s.sub.m(τ.sub.p,n)=s((τ.sub.p−1)T.sub.b+n) tm (1), where τ.sub.p=1,2, . . . , P, n=1,2, . . . , T.sub.0f.sub.s, and m=1,2, . . . , M, and P represents a number of segments into which data is divided, τ.sub.p represents a p-th segment of signal, T.sub.0 represents a length of each segment of signal, in unit of second, T.sub.b represents a segmentation stride, in unit of second, f.sub.s is a sampling rate of the signal, τ.sub.p represents a slow time, n represents a fast time, and m is an array element number.

3. The method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space according to claim 1, wherein step 2 is specifically as follows: $\begin{array}{cc}s\left({\tau }_p.\omega \right)=\frac{1}{N}{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}\left[\begin{array}{c}1\\ e^{-i\frac{2\pi f_0}{c}d\cos \alpha \left({\tau }_p\right)}\\ .Math.\\ e^{-i\frac{2\pi f_0\left(M-1\right)}{c}d\cos \alpha \left({\tau }_p\right)}\end{array}\right]e^{i\left({\omega }_p-\omega \right)n},& \left(2\right)\end{array}$ wherein within T.sub.0 seconds, a frequency ω.sub.p and azimuth α(τ.sub.p) of the received signal are approximately unchanged, and ω.sub.p=ω.sub.0+{right arrow over (k)}.Math.{right arrow over (v)} is a frequency received by the array at a time τ.sub.p after the Doppler effect caused by target movement, where {right arrow over (v)} is a target velocity vector, and {right arrow over (k)} is a wave number vector; and {right arrow over (r)}.sub.0 is a distance vector between the target and a central receiving array element at an initial time, d is an array element spacing, and both ω.sub.0 and f.sub.0 are frequencies of the target transmitting a single-frequency signal.

4. The method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space according to claim 3, wherein step 3 is specifically as follows: $\begin{array}{cc}\begin{array}{c}S\left({\tau }_p,\omega ,\alpha \right)=\underset{m=1}{\stackrel{M}{.Math.}}{s}_m\left({\tau }_p,\omega \right)e^{i\frac{2\pi f_0\left(m-1\right)}{c}d\cos \alpha }\\ =\frac{1}{N}{.Math.}_{m=1}^M{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}{e}^{i\left({\omega }_p-\omega \right)n}{e}^{-i\frac{2\pi f_0\left(m-1\right)}{c}d\left(\cos \alpha -\cos \alpha \left({\tau }_p\right)\right)}\end{array},& \left(3\right)\end{array}$ where α=α(τ.sub.p) represents the azimuth, M represents a number of array elements, m represents an m-th array element, τ.sub.p represents a p-th segment of signal, i represents an imaginary unit, cos α represents cosine of the azimuth, c represents an acoustic velocity, N represents a number of points in each segment of the segmented signal, N−1 represents a point position, and {right arrow over (k)} is the wave number vector; and {right arrow over (r)}.sub.0 is the distance vector between the target and the central receiving array element at the initial time.

5. The method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space according to claim 1, wherein step 3 is specifically as follows: according to a mathematical expression of space transformation: $\begin{array}{cc}\{\begin{array}{c}{f}^{\prime }=f\\ t^{\prime }=t\\ \cos \theta =-\sin \left(\alpha +{\gamma }_0\right)\end{array},& \left(4\right)\end{array}$ wherein γ.sub.0∈[0,360), and γ.sub.0 is an angle measured clockwise from a true north line to a target course line, searching the angle γ.sub.0, and performing coordinate transformation as shown in Formula (4) on the (f-α-t) 3D space according to a parameter to obtain a new 3D space (f-cos θ-t): $S\left(f,\cos \theta ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\cos \theta -\cos \theta \left(t\right)\right)\delta \left(f-\left[f_0+f_0\frac{v}{c}\cos \theta \left(t\right)\right]\right),$ wherein the signal is represented as a curve on a plane $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ in the new 3D space.

6. The method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space according to claim 1, wherein step 5 is specifically as follows: searching a transmitting signal frequency f.sub.0 and a velocity v, and determining a slice in the 3D space A: $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ according to the parameters.

7. A computer readable storage medium, storing a computer program, wherein when executed by a processor, the computer program implements steps of the method according to claim 1.

8. The computer readable storage medium according to claim 7, wherein step 1 is specifically as follows:
s.sub.m(τ.sub.p,n)=s((τ.sub.p−1)T.sub.b+n) tm (1), where τ.sub.p=1,2, . . . , P, n=1,2, . . . , T.sub.0f.sub.s, and m=1,2, . . . , M, and P represents a number of segments into which data is divided, τ.sub.p represents a p-th segment of signal, T.sub.0 represents a length of each segment of signal, in unit of second, T.sub.b represents a segmentation stride, in unit of second, f.sub.s is a sampling rate of the signal, τ.sub.p represents a slow time, n represents a fast time, and m is an array element number.

9. The computer readable storage medium according to claim 7, wherein step 2 is specifically as follows: $\begin{array}{cc}s\left({\tau }_p,\omega \right)=\frac{1}{N}{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}\left[\begin{array}{c}1\\ e^{-i\frac{2\pi f_0}{c}d\cos \alpha \left({\tau }_p\right)}\\ .Math.\\ e^{-i\frac{2\pi f_0\left(M-1\right)}{c}d\cos \alpha \left({\tau }_p\right)}\end{array}\right]e^{i\left({\omega }_p-\omega \right)n},& \left(2\right)\end{array}$ wherein within T.sub.0 seconds, a frequency ω.sub.p and azimuth α(τ.sub.p) of the received signal are approximately unchanged, and ω.sub.p=ω.sub.0+{right arrow over (k)}.Math.{right arrow over (v)} is a frequency received by the array at a time τ.sub.p after the Doppler effect caused by target movement, where {right arrow over (v)} is a target velocity vector, and {right arrow over (k)} is a wave number vector; and {right arrow over (r)}.sub.0 is a distance vector between the target and a central receiving array element at an initial time, d is an array element spacing, and both ω.sub.0 and f.sub.0 are frequencies of the target transmitting a single-frequency signal.

10. The computer readable storage medium according to claim 9, wherein step 3 is specifically as follows: $\begin{array}{cc}\begin{array}{c}S\left({\tau }_p,\omega ,\alpha \right)=\stackrel{M}{\underset{m=1}{.Math.}}{s}_m\left({\tau }_p,\omega \right)e^{i\frac{2\pi f_0\left(m-1\right)}{c}d\cos \alpha }\\ =\frac{1}{N}{.Math.}_{m=1}^M{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}{e}^{i\left({\omega }_p-\omega \right)n}{e}^{-i\frac{2\pi f_0\left(m-1\right)}{c}d\left(\cos \alpha -\mathrm{co}s\alpha \left({\tau }_p\right)\right)}\end{array},& \left(3\right)\end{array}$ where α=α(τ.sub.p) represents the azimuth, M represents a number of array elements, m represents an m-th array element, τ.sub.p represents a p-th segment of signal, i represents an imaginary unit, cos α represents cosine of the azimuth, c represents an acoustic velocity, N represents a number of points in each segment of the segmented signal, N−1 represents a point position, and {right arrow over (k)} is the wave number vector; and {right arrow over (r)}.sub.0 is the distance vector between the target and the central receiving array element at the initial time.

11. The computer readable storage medium according to claim 7, wherein step 3 is specifically as follows: according to a mathematical expression of space transformation: $\begin{array}{cc}\{\begin{array}{c}{f}^{\prime }=f\\ t^{\prime }=t\\ \cos \theta =-\sin \left(\alpha -{\gamma }_0\right)\end{array},& \left(4\right)\end{array}$ wherein γ.sub.0∈[0,360), and γ.sub.0 is an angle measured clockwise from a true north line to a target course line, searching the angle γ.sub.0, and performing coordinate transformation as shown in Formula (4) on the (f-α-t) 3D space according to the parameter to obtain a new 3D space (f-cos θ-t): $S\left(f,\cos \theta ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\cos \theta -\cos \theta \left(t\right)\right)\delta \left(f-\left[f_0+f_0\frac{v}{c}\cos \theta \left(t\right)\right]\right),$ wherein the signal is represented as a curve on a plane $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ in the new 3D space.

12. The computer readable storage medium according to claim 7, wherein step 5 is specifically as follows: searching a transmitting signal frequency f.sub.0 and a velocity v, and determining a slice in the 3D space A: $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ according to the parameters.

13. An electronic device, comprising a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; the memory is configured to store a computer program; and the processor is configured to implement steps of the method according to claim 1 when executing the program stored in the memory.

14. The electronic device according to claim 13, wherein step 1 is specifically as follows:
s.sub.m(τ.sub.p,n)=s((τ.sub.p−1)T.sub.b+n) tm (1), where τ.sub.p=1,2, . . . , P, n=1,2, . . . , T.sub.0f.sub.s, and m=1,2, . . . , M, and P represents a number of segments into which data is divided, τ.sub.p represents a p-th segment of signal, T.sub.0 represents a length of each segment of signal, in unit of second, T.sub.b represents a segmentation stride, in unit of second, f.sub.s is a sampling rate of the signal, τ.sub.p represents a slow time, n represents a fast time, and m is an array element number.

15. The electronic device according to claim 13, wherein step 2 is specifically as follows: $\begin{array}{cc}s\left({\tau }_p,\omega \right)=\frac{1}{N}{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}\left[\begin{array}{c}1\\ e^{-i\frac{2\pi f_0}{c}d\cos \alpha \left({\tau }_p\right)}\\ .Math.\\ e^{-i\frac{2\pi f_0\left(M-1\right)}{c}d\cos \alpha \left({\tau }_p\right)}\end{array}\right]e^{i\left({\omega }_p-\omega \right)n},& \left(2\right)\end{array}$ within T.sub.0 seconds, a frequency ω.sub.p and azimuth α(τ.sub.p) of the received signal are approximately unchanged, and ω.sub.p=ω.sub.0+{right arrow over (k)}.Math.{right arrow over (v)} is a frequency received by the array at a time τ.sub.p after the Doppler effect caused by target movement, where {right arrow over (v)} is a target velocity vector, and {right arrow over (k)} is a wave number vector; and {right arrow over (r)}.sub.0 is a distance vector between the target and a central receiving array element at an initial time, d is an array element spacing, and both ω.sub.0 and f.sub.0 are frequencies of the target transmitting a single-frequency signal.

16. The electronic device according to claim 15, wherein step 3 is specifically as follows: $\begin{array}{cc}\begin{array}{c}S\left({\tau }_p,\omega ,\alpha \right)=\stackrel{M}{\underset{m=1}{.Math.}}{s}_m\left({\tau }_p,\omega \right)e^{i\frac{2\pi f_0\left(m-1\right)}{c}d\cos \alpha }\\ =\frac{1}{N}{.Math.}_{m=1}^M{.Math.}_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}{e}^{i\left({\omega }_p-\omega \right)n}{e}^{-i\frac{2\pi f_0\left(m-1\right)}{c}d\left(\cos \alpha -\mathrm{co}s\alpha \left({\tau }_p\right)\right)}\end{array},& \left(3\right)\end{array}$ where α=α(τ.sub.p) represents the azimuth, M represents a number of array elements, m represents an m-th array element, τ.sub.p represents a p-th segment of signal, i represents an imaginary unit, cos α represents cosine of the azimuth, c represents an acoustic velocity, N represents a number of points in each segment of the segmented signal, N−1 represents a point position, and {right arrow over (k)} is the wave number vector; and {right arrow over (r)}.sub.0 is the distance vector between the target and the central receiving array element at the initial time.

17. The electronic device according to claim 13, wherein step 3 is specifically as follows: according to a mathematical expression of space transformation: $\begin{array}{cc}\{\begin{array}{c}{f}^{\prime }=f\\ t^{\prime }=t\\ \cos \theta =-\sin \left(\alpha -{\gamma }_0\right)\end{array},& \left(4\right)\end{array}$ wherein γ.sub.0∈[0,360), and γ.sub.0 is an angle measured clockwise from a true north line to a target course line, searching the angle γ.sub.0, and performing coordinate transformation as shown in Formula (4) on the (f-α-t) 3D space according to the parameter to obtain a new 3D space (f-cos θ-t): $S\left(f,\cos \theta ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\cos \theta -\cos \theta \left(t\right)\right)\delta \left(f-\left[f_0+f_0\frac{v}{c}\cos \theta \left(t\right)\right]\right),$ wherein the signal is represented as a curve on a plane $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ in the new 3D space.

18. The electronic device according to claim 13, wherein step 5 is specifically as follows: searching a transmitting signal frequency f.sub.0 and a velocity v, and determining a slice in the 3D space A: $f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$ according to the parameters.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0040] FIG. 1 is a schematic diagram of a target movement model of the present disclosure;

[0041] FIG. 2 is a flowchart of a method of the present disclosure;

[0042] FIGS. 3A-D are schematic diagrams of a position of a signal in a 3D space of the present disclosure, where FIG. 3A shows a curve diagram of the signal in a 3D space (f-α-t), FIG. 3B shows a top view of the 3D space (f-α-t), FIG. 3C shows a side view of the 3D space (f-α-t), and FIG. 3D shows a front view of the 3D space (f-α-t);

[0043] FIGS. 4A-D are schematic diagrams of a position of a signal in a 3D space of the present disclosure, where FIG. 4A shows a curve diagram of the signal in a 3D space (f-cos θ-t), FIG. 4B shows a top view of the 3D space (f-cos θ-t), FIG. 4C shows a side view of the 3D space (f-cos θ-t), and FIG. 4D shows a front view of the 3D space (f-cos θ-t);

[0044] FIG. 5 is a schematic diagram of taking a slice in the 3D space of the present disclosure;

[0045] FIGS. 6A-D are LOFAR diagrams of the present disclosure, where FIG. 6A shows a 64 Hz LOFAR diagram, FIG. 6B shows a 67 Hz LOFAR diagram, FIG. 6C shows a 70 Hz LOFAR diagram, and FIG. 6D shows a 73 Hz LOFAR diagram;

[0046] FIGS. 7A-D are azimuth history diagrams of the present disclosure, where FIG. 7A shows a 64 Hz azimuth history diagram, FIG. 7B shows a 67 Hz azimuth history diagram, FIG. 7C shows a 70 Hz azimuth history diagram, and FIG. 7D shows a 73 Hz azimuth history diagram; and

[0047] FIGS. 8A-D are effect diagrams of the present disclosure, where FIG. 8A shows a 64 Hz effect diagram, FIG. 8B shows a 67 Hz effect diagram, FIG. 8C shows a 70 Hz effect diagram, and FIG. 8D shows a 73 Hz effect diagram.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0048] The technical solutions of the embodiments of the present disclosure are clearly and completely described below with reference to the accompanying drawings. Apparently, the described embodiments are only illustrative ones, and are not all possible ones of the present disclosure. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

[0049] The present disclosure provides a method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, including the following steps.

[0050] Step 1: A target radiated acoustic signal s(t) received by an M-element horizontal line array in an underwater acoustic environment with a low SNR is segmented.

[0051] Step 2: N-point DFT is performed on the received signal on each array element in each period of time τ.sub.p in step 1, where N=T.sub.0.Math.f.sub.s.

[0052] For the azimuth α(τ.sub.p), each period of time τ.sub.p has an azimuth α(τ.sub.p), and step S3 compensates this azimuth for each segment of signal. However, during compensation, the azimuth is unknown, so it is compensated once in sequence from 0 to 180 degrees. In the formula, α is a variable, and α(τ.sub.p) is a constant.

[0053] Step 3: Frequency domain beamforming is performed on an array signal after each section of DFT in step 2, and stacking is performed after a phase difference between arrays brought by an azimuth α(τ.sub.p) of each primitive element is compensated.

[0054] Step 4: Coordinate transformation is performed on a frequency-azimuth-time (f-α-t) 3D matrix space obtained in step 3.

[0055] Step 5: A slice is taken from the frequency-azimuth-time (f-cos θ-t) 3D space subjected to the coordinate transformation obtained in step 4.

[0056] Step 6: Segmented Radon transform is performed on the spatial slice obtained in step 5 to detect the target.

[0057] In the method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, step 1 is specifically as follows:

(τ.sub.p,n)=s((τ.sub.p−1)T.sub.b+n) tm (1),

where [0058] τ.sub.p=1,2, . . . , P, n=1,2, . . . , T.sub.0f.sub.s, and m=1,2, . . . , M, and

[0059] P represents a number of segments into which data is divided, τ.sub.p represents a p-th segment of signal, T.sub.0 represents a length of each segment of signal, in unit of second, T.sub.b represents a segmentation stride, in unit of second, f.sub.s is a sampling rate of the signal, τ.sub.p represents a slow time, n represents a fast time, and m is an array element number.

[0060] In the method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, step 2 is specifically as follows:

[00007]$\begin{array}{cc}s\left({\tau }_p,\omega \right)=\frac{1}{N}{.Math.}_{n=0}^{N-1}{\mathrm{Ae}}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}\left[\begin{array}{c}\begin{array}{c}\begin{array}{c}1\\ e^{-i\frac{2\pi f_0}{c}\mathrm{dcos}\alpha \left({\tau }_p\right)}\end{array}\\ .Math.\end{array}\\ e^{-i\frac{2\pi f_0\left(M-1\right)}{c}\mathrm{dcos}\alpha \left({\tau }_p\right)}\end{array}\right]e^{i\left({\omega }_p-\omega \right)n},& \left(2\right)\end{array}$

[0061] Within T.sub.0 seconds, a frequency ω.sub.p and azimuth α(τ.sub.p) of the received signal are approximately unchanged, and ω.sub.p=ω.sub.0+{right arrow over (k)}.Math.{right arrow over (v)} is a frequency received by the array at a time τ.sub.p after the Doppler effect caused by target movement, where {right arrow over (v)} is a target velocity vector, and {right arrow over (k)} is a wave number vector; and {right arrow over (r)}.sub.0 is a distance vector between the target and a central receiving array element at an initial time, d is an array element spacing, and both ω.sub.0 and f.sub.0 are frequencies of the target transmitting a single-frequency signal.

[0062] In the method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, step 3 is specifically as follows:

[00008]$\begin{array}{cc}S\left({\tau }_p,\omega ,\alpha \right)=\underset{m=1}{\overset{M}{.Math.}}{s}_m\left({\tau }_p,\omega \right)e^{i\frac{2\pi f_0\left(m-1\right)}{c}\mathrm{dcos}\alpha }=\frac{1}{N}{\Sigma }_{m=1}^M{\Sigma }_{n=0}^{N-1}A{e}^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}{e}^{i{\omega }_p{\tau }_p}{e}^{i\left({\omega }_p-\omega \right)n}{e}^{-i\frac{2\pi f_0\left(m-1\right)}{c}d\left(\cos \alpha -\cos \alpha \left({\tau }_p\right)\right)},& \left(3\right)\end{array}$ [0063] where α=α(τ.sub.p) represents the azimuth, M represents a number of array elements, m represents an m-th array element, τ.sub.p represents a p-th segment of signal, i represents an imaginary unit, cos α represents cosine of the azimuth, c represents an acoustic velocity, N represents a number of points in each segment of the segmented signal, N−1 represents a point position, and {right arrow over (k)} is the wave number vector; and {right arrow over (r)}.sub.0 is the distance vector between the target and the central receiving array element at the initial time.

[0064] In the method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, step 3 is specifically as follows.

[0065] According to a mathematical expression of space transformation:

[00009]$\begin{array}{cc}\{\begin{array}{c}{f}^{\prime }=f\\ t^{\prime }=t\\ \cos \theta =-\sin \left(\alpha -{\gamma }_0\right)\end{array},& \left(4\right)\end{array}$ [0066] where γ.sub.0∈[0,360), and γ.sub.0 is an angle measured clockwise from a true north line to a target course line, and is a constant, [0067] the angle γ.sub.0 is searched, and coordinate transformation as shown in Formula (4) is performed on the (f-α-t) 3D space according to the parameter to obtain a new 3D space (f-cos θ-t):

[00010]$S\left(f,\cos \theta ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\cos \theta -\cos \theta \left(t\right)\right)\delta \left(f-\left[f_0+f_0\frac{v}{c}\cos \theta \left(t\right)\right]\right).$

[0068] The signal is represented as a curve on a plane

[00011]$f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$

in the new 3D space.

[0069] In the method for detecting a moving target based on spatial slices of transformed spatio-temporal frequency space, step 5 is specifically as follows: searching a transmitting signal frequency f.sub.0 and a velocity v, and determining a slice in the 3D space A:

[00012]$f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$

according to the parameters.

[0070] The present disclosure provides a computer readable storage medium, storing a computer program. When executed by a processor, the computer program implements steps of the method described above.

[0071] The present disclosure provides an electronic device, including a processor, a communication interface, a memory, and a communication bus. The processor, the communication interface, and the memory communicate with each other through the communication bus.

[0072] The memory is configured to store a computer program.

[0073] The processor is configured to implement steps of the method described above when executing the program stored in the memory.

[0074] FIG. 1 shows a target movement model established herein. First of all, it is assumed that the receiving device placed underwater is a uniform horizontal line array with M array elements. The reference array and the equivalent acoustic center of the line array are located at the origin of coordinates, and the moving target moves uniformly in a straight line from the initial point A(x.sub.0, y.sub.0) at {right arrow over (v)}, radiating a single frequency signal with an angular frequency of ω.sub.0. The course γ.sub.0, γ.sub.0∈[0,360) is an angle measured clockwise from a true north line to a target course line, and is a constant. At a time t, the target moves to the point B, and the relative bearing θ, θ∈[0,180] is an included angle between the moving direction of the target and the distance vector {right arrow over (BO)} of the target position. α is the target azimuth, that is, the angle between the vector {right arrow over (OB)} and the x-axis of the line between the target position and the equivalent acoustic center of the array. θ and α change continuously with the target movement or time, which can be expressed as the function θ(t) and α(t) of the time t. Only considering the movement of the target in the first and second quadrants, the position relation between the three angles cos θ, α, and γ.sub.0 is:

cos θ(t)=−sin(α(t)+γ.sub.0).

[0075] Assuming that the acoustic field is an infinite free plane and the receiving array is located in the far field of the signal source, the received signal can be approximately considered as a plane wave. Within a period of time, due to the Doppler effect brought by the target movement, the signal received by the array element at the origin of coordinates is:

[00013]$s_1\left(t\right)=e^{i\left(\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0-\left({\omega }_0+\stackrel{.fwdarw.}{k}.Math.\stackrel{.fwdarw.}{v}\right)t\right)},$ [0076] where {right arrow over (k)} is a wave number vector. The coordinate of any array element located on the x-axis can be expressed as (x, 0). According to the plane wave incidence principle, the phase difference between the received signal of this array element and the central array element is

[00014]$e^{-i\frac{{\omega }_0}{c}x\cos \alpha },$

then the signal received by any array element can be written as:

[00015]$s\left(t,x\right)=e^{i\left(\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0-\left({\omega }_0+\frac{{\omega }_0}{c}v\cos \theta \right)t-\frac{{\omega }_0}{c}x\cos \alpha \left(t\right)\right)}.$

[0077] According to the moving target model in FIG. 1, it can be seen that the azimuth of the target changes with time, and changes with the radial velocity of the target movement. Due to the Doppler effect, the received frequency also changes with time, and a non-stationary signal is received. Therefore, the above array signals are subjected to the spatio-temporal two-dimensional (2D) short-time Fourier transform, that is, the time window function γ* (t′−t) is multiplied by the source signal s(t′, x) to achieve the interception and translation of the source signal, and the signal in the spatio-temporal domain is subjected to the 2D Fourier transform to obtain the frequency-azimuth domain of the signal. The result after transformation can be expressed by the following formula:

S(f,α,t)=∫.sub.−∞.sup.+∞∫.sub.−∞.sup.+∞s(t′,x)γ*(t′−t)e.sup.−i2π(ft′+x cos α)dtdx.

[0078] By substituting Formula (3), the formula is simplified into:

[00016]$S\left(f,\alpha ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\alpha -\alpha \left(t\right)\right)\delta \left(f-f\left({\tau }_p\right)\right).$

[0079] By sliding the window function and performing the 2D Fourier transform, the function expression of the received signal in the 3D space represented by the three axes of frequency-azimuth-time (f-α-t) is obtained. As shown in FIGS. 3A-D, it is a curve on the surface in the 3D space.

[0080] According to the above derivation, the frequency of the received signal at any time is:

[00017]$f\left(t\right)=f_0+f_0\frac{v}{c}\cos \theta \left(t\right).$

[0081] Since f.sub.0, v, and c are all constants, the relationship between f(t) and cos θ(t) is linear. Therefore, according to the mathematical expression of space transformation:

[00018]$\{\begin{array}{c}{f}^{\prime }=f\\ t^{\prime }=t\\ \cos \theta =-\sin \left(\alpha +{\gamma }_0\right)\end{array},$

[0082] the angle γ.sub.0 is searched, and coordinate transformation as shown in Formula (6) is performed on the (f-α-t) 3D space according to the parameter to obtain a new 3D space (f-cos θ-t). The position of the signal in (f-cos θ-t) is shown in FIGS. 4A-D. Since the relationship between f(t) and cos θ(t) is linear, it can be seen that the points representing the signal in the top view of FIG. 4 (B) form a straight line, that is, the signal can be represented as a curve on the plane

[00019]$f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$

in the new 3D space, whose equation is:

[00020]$S\left(f,\cos \theta ,t\right)=e^{i\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{r}}_0}\delta \left(\cos \theta -\cos \theta \left(t\right)\right)\delta \left(f-\left[f_0+f_0\frac{v}{c}\cos \theta \left(t\right)\right]\right).$

[0083] A transmitting signal frequency f.sub.0 and a velocity v are searched, and a slice in the 3D space A:

[00021]$f-\left(f_0+f_0\frac{v}{c}\cos \theta \right)=0$

can be determined according to the parameters, that is, the plane represented by the blue block diagram in FIG. 5.

[0084] As shown in FIGS. 8A-D, 0-32 min data of the Swellex96 S5 experiment is processed by the present disclosure, and compared with the LOFAR diagram in FIGS. 6A-D and azimuth history diagram in FIGS. 7A-D. The left figure is the slice taken from the space, and the right figure is the result after Radon transform of the slice.