CN111323750B - Direct positioning method based on acoustic vector array network - Google Patents
Direct positioning method based on acoustic vector array network Download PDFInfo
- Publication number
- CN111323750B CN111323750B CN202010197246.8A CN202010197246A CN111323750B CN 111323750 B CN111323750 B CN 111323750B CN 202010197246 A CN202010197246 A CN 202010197246A CN 111323750 B CN111323750 B CN 111323750B
- Authority
- CN
- China
- Prior art keywords
- array
- positioning
- sound pressure
- calculating
- spatial spectrum
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
Abstract
The invention provides a direct positioning method based on an acoustic vector array network, and belongs to the field of passive target positioning. Firstly, carrying out segmented frequency domain transformation on sound pressure and vibration speed data received by each array node; then scanning point by point in the area to be measured, and calculating the guiding direction and time delay at the scanning point; weighting and summing the sound pressure and vibration velocity data to obtain weighted data matrixes of all the arrays; calculating equivalent array manifold by using a frequency domain positioning model; then, calculating the spatial spectrum on each frequency by using a spatial spectrum estimation algorithm, and accumulating to obtain the total spatial spectrum of all frequencies; and calculating the spatial spectrum of each point in the region to be measured according to the proper scanning step length, and searching the peak value of the spatial spectrum value to obtain the position of the target to be measured. The invention introduces the acoustic vector array into the multi-array target direct positioning method, obtains higher positioning precision and spatial resolution than a scalar array by utilizing the sound pressure and vibration velocity combined processing technology, and has better adaptability to target positioning under the condition of low signal-to-noise ratio.
Description
Technical Field
The invention relates to a direct positioning method based on an acoustic vector array network, belonging to the field of passive target positioning.
Background
The target passive positioning method utilizes the signal radiated or reflected by the target to position and track, has strong concealment and long action distance, and has wide application in the fields of regional warning, environmental monitoring and the like. The multi-station passive positioning system takes the sensor arrays as nodes to form a distributed array network, utilizes information provided by the arrays to carry out cooperative processing, obtains target position parameters, and has higher positioning precision compared with single-point positioning.
The passive positioning method of an object using an array network is generally classified into a two-step positioning method and a direct positioning method. The two-step positioning method comprises the steps of firstly calculating intermediate positioning parameters such as the target direction of arrival (AOA), the time delay Difference (DTOA), the Doppler frequency shift or the signal intensity, and then calculating the target position by using the intermediate parameters. The method only needs to transmit the intermediate positioning parameters to the resolving center, and has low requirements on communication bandwidth and calculated amount. However, since the estimation of the parameters is completed by each node in isolation, part of information is inevitably lost, the precision is reduced, and meanwhile, the problem of multi-target parameter matching is also brought. Compared with the two-step positioning method, the direct positioning method (DPD) directly utilizes the original data output by the bottom layer array, estimates the position information through single-step calculation, avoids the defects of the two-step positioning method, and has better performance than the two-step positioning method under the condition of low signal-to-noise ratio.
The direct positioning method (DPD) was first proposed by Weiss a J in 2004. According to the method, the processing idea of non-coherent subspace decomposition in broadband spectrum estimation is used for reference, time delay information is separated from a frequency domain through segmented Discrete Fourier Transform (DFT), an equivalent array manifold matrix is defined, and a target position is solved by using a traditional spatial spectrum estimation method. Such as multi-target maximum likelihood direct localization algorithms for known waveforms (Amar A, Weiss A J. direct localization (DPD) of multiple known and unwwn radio-frequency signals [ C ]. 200412 th European Signal Processing Conference, Vienna,2004:1115-1118.), and MUSIC-based direct localization algorithms for unknown waveforms (Weiss AJ, Amar A. direct localization of multiple radio signals [ J ]. EURASIP Journal Advances in Signal Processing,2005, 37-49). In 2015, Tom Tirer and Anthony J.Weiss introduced the minimum variance distortion free response (MVDR) spatial spectrum estimation algorithm into the Direct positioning method (Tirer T, Weiss AJ.high Resolution Direct Position Determination of Radio Frequency Sources [ J ]. IEEE Signal Processing Letters,2015,23(2): 192-. However, most of the array nodes involved in the above direct localization algorithm are based on scalar sensors, and target localization is achieved only by using scalar information, and vector information in a physical field is not considered.
In fact, the physical quantities describing the sound field have a medium particle velocity vector in addition to the sound pressure. Due to the limitation of technical conditions, the traditional sonar only uses scalar sound pressure information and ignores vector vibration velocity information, and the sound vector sensor can synchronously pick up the sound pressure scalar and the vibration velocity vector of a sound field in a space concurrent mode, so that the description of the sound field is more accurate and comprehensive. Compared with an acoustic scalar array, the acoustic vector array has the advantages of high array gain, low detection domain, strong anti-noise interference capability and the like. In 1994, Arye Nehorai first proposed an acoustic vector sensor Signal Processing framework (Nehorai A, paladi E. acoustic vector-sensor array Processing [ J ]. IEEE Transactions on Signal Processing,1994,42(9):2481 and 2491.), and sound pressure and vibration velocity information were processed as independent array elements. The national scholars huijingi et al have set forth the concept and method of sound pressure and vibration velocity combined processing (huijingi. vector acoustic signal processing foundation [ M ]. national defense industry press, 2009.) and theoretically demonstrate the anti-noise advantage of the sound pressure and vibration velocity combined processing technology.
Disclosure of Invention
The invention provides a direct positioning method based on an acoustic vector array network. The invention combines the sound pressure and vibration velocity combined processing technology with the idea of a direct positioning method, and outputs a high-resolution space spectrogram of a target detection area by using a space spectrum estimation method.
The purpose of the invention is realized as follows: the node of the acoustic vector array network is a multi-array element acoustic vector sensor array, each array element can output the same-point sound pressure scalar and the same-point vibration velocity vector information, and the method comprises the following steps:
Step 3, scanning point by point in the area to be measured, and calculating the guiding azimuth psi from the point to the I-th array reference point at the scanning point d l (d) And time delay tau l (d)；
Step 4, using the guide orientation psi l (d) Sound pressure P to the same segment l (k, j) and the vibration velocity V xl (k,j),V yl (k, j) combined treatment to give Y l (k, j), summing Y for all L arrays l (k, j) obtaining a total data matrix Y (k, j);
step 5, defining and calculating equivalent array manifold by using the array positioning frequency domain signal model
Summing all L arrays
Obtaining the total equivalent array manifold matrix
step 7, obtaining a total space spectrum Q (d) of accumulated K frequencies;
and 8, setting a proper scanning step length, repeating the steps 3 to 7, traversing the whole region to be detected, outputting a spatial spectrogram, and searching a spectral value peak value to obtain the position of the target to be detected.
The invention also includes such structural features:
1. step 4 utilizing psi l (d) Sound pressure P to the same segment l (k, j) and the vibration velocity V xl (k,j),V yl (k, j) summing to obtain the vibration velocity of the j-th section sound pressure and the Fourier transform Y of the k-th frequency component l (k,j)：
Y l (k,j)＝P l (k,j)+V xl (k,j)cosψ l (d)+V yl (k,j)sinψ l (d)
Summing Y of all L arrays l (k, j) obtaining a total array data matrix Y (k, j) ═ Y 1 (k,j) T ,…,Y L (k,j) T ] T 。
2. The step 5 specifically comprises the following steps:
step 5-1, outputting Y by the array l in the positioning network l The frequency domain model of (k, j) is:
in the formula, S (k, j), N l (k, j) are frequency domain representations of the source and noise, respectively;
step 5-2, defining and calculating the equivalent array manifold according to the following formula
Summing all L arrays
Is a total equivalent array manifold matrix
Step 5-3, order
Due to b l For unknown complex scalars, it is usual to set 1, b | | | b | |, b |, and b |, are set 1 ,…,b L ] T Then the array positioning frequency domain signal model can be expressed as:
the frequency domain signal model is in accordance with the conventional array signal model.
3. Step 6 relates to a spatial spectrum estimation method including but not limited to CBF, MVDR, MUSIC, etc., taking a MUSIC spatial spectrum estimation algorithm as an example, the method specifically includes the following steps:
step 6-1, extracting fourier transform Y (k, J), J being 1, …, J, constituting Y (k), Y (k) ═ Y (k,1), …, Y (k, J) of J-segment data in k-th frequency component]By the use of R k ＝Y(k)·Y H (k) J calculating frequency f k Upper covariance matrix, R k Is an LM multiplied by LM dimensional covariance matrix;
step 6-2, to covariance matrix R k Performing characteristic decomposition to obtain a signal subspace U consisting of characteristic vectors corresponding to larger D characteristic values S (k) And a noise subspace U consisting of eigenvectors corresponding to other smaller eigenvalues N (k) Wherein D is the number of targets;
step 6-3, at frequency f k The MUSIC spatial spectrum above is:
in the formula (I), the compound is shown in the specification,
I L is L-L dimensional Unit array, J M Is an M multiplied by 1 dimensional full 1 array,
representing the Kronecker product, the above formula is equivalent to:
order to
Since b is 1, G (k, d) is a quadratic form under the constraint, and its eigenvalue is the maximum value of G (k, d) eigenvalue, Q (k, d) is:
Q(k,d)＝max(λ max (G(k,d)))
4. step 7 sums all frequency-space spectra Q (k, d), the total space spectrum Q (d) being:
compared with the prior art, the invention has the beneficial effects that: the invention introduces an acoustic vector sensor into a target direct positioning method of multiple array nodes, and provides a direct positioning method based on an acoustic vector array network. The method combines the superior performance of the acoustic vector array with the advantage of a direct positioning method, and obtains a high-resolution space spectrogram of a target detection area by using a space spectrum estimation method. The method has higher positioning precision and spatial resolution, and has better adaptability to target positioning under the condition of low signal-to-noise ratio.
Drawings
FIG. 1 is a schematic diagram of the distribution of a region to be measured;
FIGS. 2a-b are spatial spectrum pseudocolor maps for single source direct localization at a signal-to-noise ratio of-15 dB, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm;
FIGS. 3a-b are three-dimensional plots of spatial spectra for direct localization by a single source at a signal-to-noise ratio of-15 dB, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm;
FIGS. 4a-b are spatial spectrum pseudocolor maps for single source direct localization at-20 dB, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm;
FIG. 5 is a three-dimensional plot of the spatial spectrum for direct localization by a single source at a signal-to-noise ratio of-20 dB, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm;
FIGS. 6a-b are spatial spectrum pseudocolor maps for direct localization of dual sources at a signal-to-noise ratio of 0dB, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm;
FIGS. 7a-b are spatial spectrum three-dimensional plots for 0dB, dual source direct localization, and (a) scalar MUSIC-DPD algorithm; (b) vector MUSIC-DPD algorithm
FIG. 8 is a graph showing the variation of the positioning error under different SNR conditions.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
The invention combines the sound pressure and vibration velocity combined processing technology with the idea of a direct positioning method, and outputs a high-resolution space spectrogram of a target detection area by using a space spectrum estimation method. The method has higher positioning precision and spatial resolution, and has better adaptability to target positioning under the condition of low signal-to-noise ratio.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps: the node of the acoustic vector array network is a multi-array element acoustic vector sensor array, each array element can output the same-point sound pressure scalar and vibration velocity vector information,
Step 3, gradually adding the materials in the area to be measuredPoint scanning, at scanning point d, the guide orientation psi of the point to the ith array datum position is calculated l (d) And time delay tau l (d)；
Step 4, using the guide orientation psi l (d) Sound pressure P to the same segment l (k, j) and the vibration velocity V xl (k,j),V yl (k, j) performing combined treatment to obtain Y l (k, j), summing Y for all L arrays l (k, j) obtaining a total data matrix Y (k, j);
step 5, defining and calculating equivalent array manifold by using the array positioning frequency domain signal model
Summing all L arrays
Obtaining the total equivalent array manifold matrix
step 7, obtaining a total space spectrum Q (d) of accumulated K frequencies;
and 8, setting a proper scanning step length, repeating the steps 3 to 7, traversing the whole region to be detected, outputting a spatial spectrogram, and searching a spectral value peak value to obtain the position of the target to be detected.
Step 4 of the invention utilizes psi l (d) To the sound pressure P of the same segment l (k, j) and the vibration velocity V xl (k,j),V yl (k, j) summing to obtain the vibration velocity of the j-th section sound pressure and the Fourier transform Y of the k-th frequency component l (k,j)：
Y l (k,j)＝P l (k,j)+V xl (k,j)cosψ l (d)+V yl (k,j)sinψ l (d)
Summing Y of all L arrays l (k, j) obtaining a total array data matrix Y (k, j) [ [ solution ] ]Y 1 (k,j) T ,…,Y L (k,j) T ] T 。
The step 5 of the invention specifically comprises the following steps:
step 5-1, outputting Y by the array l in the positioning network l The frequency domain model of (k, j) is:
in the formula, S (k, j), N l (k, j) are frequency domain representations of the source and noise, respectively;
step 5-2, defining and calculating the equivalent array manifold according to the following formula
Summing all L arrays
Is a total equivalent array manifold matrix
Step 5-3, order
Due to b l For unknown complex scalars, it is usual to set 1, b | | | b | |, b |, and b |, are set 1 ,…,b L ] T Then the array positioning frequency domain signal model can be expressed as:
the frequency domain signal model is in accordance with the conventional array signal model.
The step 6 of the invention relates to a spatial spectrum estimation method which comprises, but is not limited to CBF, MVDR, MUSIC and the like, taking a MUSIC spatial spectrum estimation algorithm as an example, the method specifically comprises the following steps:
step 6-1, extracting fourier transform Y (k, J), J being 1, …, J, constituting Y (k), Y (k) ═ Y (k,1), …, Y (k, J) of J-segment data in k-th frequency component]By the use of R k ＝Y(k)·Y H (k) J calculating frequency f k Upper covariance matrix, R k Is an LM multiplied by LM dimensional covariance matrix;
step 6-2, to covariance matrix R k Performing characteristic decomposition to obtain a signal subspace U consisting of characteristic vectors corresponding to larger D characteristic values S (k) And a noise subspace U consisting of eigenvectors corresponding to other smaller eigenvalues N (k) Wherein D is the target number;
step 6-3, calculating the frequency f by using the following formula k MUSIC spatial spectrum above:
in the formula (I), the compound is shown in the specification,
I L is a L × L dimensional unit array, J M Is an M multiplied by 1 dimensional full 1 array,
representing the Kronecker product, the above formula is equivalent to:
order to
Since | | | b | | | 1, G (k, d) is a quadratic form under the constraint, whose characteristic value is the maximum value of the characteristic value of G (k, d), Q (k, d) can be expressed as:
Q(k,d)＝max(λ max (G(k,d)))
step 7 of the present invention sums all frequency space spectra Q (k, d), and the total space spectrum Q (d) can be expressed as:
the invention is described in detail with reference to the accompanying drawings:
Sound pressure data P outputted from the ith array node l (t) and vector vibration velocity data V xl (t),V yl (t) can be represented by
In the formula, P l (t),V xl (t),V yl (t) respectively representing sound pressure received by the ith array and vector vibration velocity signals in x and y directions, S (t) is an information source time domain snapshot matrix, N pl (t),N xl (t),N yl (t) receiving noise of sound pressure and x, y direction vibration velocity channels respectively; b l To the attenuation coefficient, α l (d) For an M x D dimensional array manifold, tau l (d) The time delay from the source at the d position to the l array is obtained;
the d-th position source direction.
In the formula, S (k, j), N pl (k,j),N xl (k,j),N yl (k, j) are S (t), N pl (t),N xl (t),N yl (t) discrete fourier transform of the jth segment of data at the kth frequency;
step 3, scanning point by point in the area to be measured, and scanning the area in the coordinate (x) d ,y d ) At the scanning point d, the guiding orientation ψ from the scanning point d to the ith array reference point position is calculated by the following formula l (d) And time delay tau l (d)：
In the formula (x) d ,y d ) As the position coordinates of the scanning point, (x) l ,y l ) Is the coordinate of the position of the datum point of the array l, and C is the sound velocity.
Step 4, using the guide orientation psi l (d) For P of the same segment l (k, j) and V xl (k,j),V yl (k, j) summing to obtain the vibration velocity of the j-th section sound pressure and the Fourier transform Y of the k-th frequency component l (k,j)：
Y l (k,j)＝P l (k,j)+V xl (k,j)cosψ l (d)+V yl (k,j)sinψ l (d) (4)
And Y is l (k, j) can also be expressed as:
in the formula, N l (k, j) is N pl (k,j),N xl (k,j),N yl The sum of (k, j);
summing Y of all L arrays l (k, j) obtaining a total array data matrix Y (k, j) ═ Y 1 (k,j) T ,…,Y L (k,j) T ] T 。
Step 5, defining and calculating equivalent array manifold
Summing all L arrays
Obtaining the total equivalent array manifold matrix
Order to
b l For unknown complex scalars, it is usual to set 1, b | | | b | |, b |, and b |, are set 1 ,…,b L ] T Then the array positioning frequency domain signal model can be expressed as:
the frequency domain signal model represented by equation (6) is consistent with the conventional array signal model form.
The frequency f is calculated using the following formula k MUSIC spatial spectrum above:
in the formula (I), the compound is shown in the specification,
I L is a L × L dimensional unit array, J M Is an M multiplied by 1 dimensional full 1 array,
representing the Kronecker product, the above formula is equivalent to:
order to
Since | | | b | | | 1, G (k, d) is a quadratic form under the constraint, whose characteristic value is the maximum value of the characteristic value of G (k, d), Q (k, d) can be expressed as:
Q(k,d)＝max(λ max (G(k,d))) (9)
step 7, the total spatial spectrum q (d) of all frequencies can be expressed as:
and 8, setting a proper scanning step length, repeating the steps 3 to 7, traversing the whole region to be detected, outputting a spatial spectrogram, and searching a spectral value peak value to obtain the position of the target to be detected.
The above description is directed to the embodiments of the present invention, and the following description is directed to the simulation examples.
Example one: single source processing effect analysis
Example parameter settings are as follows: as shown in fig. 1, in a 500m × 500m area to be measured, a narrowband signal source to be measured is located at (200,300) m, the number of enabled array nodes is 2, and the enabled array nodes are respectively located at (0,0) m and (500,0) m positions, a network node is an 8-element uniform linear array, the array element interval of the array is a half wavelength, the horizontal direction of the linear array is consistent with the set x direction, and environmental noise is stable narrowband gaussian white noise. The sound velocity is set to be 1500m/s, the sampling rate of the system is 2kHz, the fast beat number is 320, and the received data is evenly divided into 8 sections, and each section has 40 points. And traversing the whole area to be measured, wherein the scanning step length is 1 m.
And analyzing the positioning performance of the scalar MUSIC direct positioning algorithm and the vector MUSIC direct positioning algorithm on the single information source under the conditions of different signal-to-noise ratios. FIG. 2 and FIG. 3 show a spatial spectrum pseudo-color map and a three-dimensional map of single source positioning by two direct positioning algorithms when the signal-to-noise ratio is-15 dB, respectively; FIGS. 4 and 5 show the spatial spectrum pseudo-color map and three-dimensional map of single source positioning by two direct positioning algorithms when the signal-to-noise ratio is-20 dB, respectively;
comparing fig. 2-5, it can be seen that under the signal-to-noise ratio-15 dB condition, the scalar MUSIC direct positioning algorithm is affected by noise, and the background fluctuation is obvious, compared with the vector MUSIC direct positioning method, the space spectrum estimation background is smooth, the spectrum peak is more obvious and prominent, and the positioning accuracy is better; when the signal-to-noise ratio is-20 dB, the scalar MUSIC direct positioning algorithm basically fails, and the vector MUSIC direct positioning method using the vector array has stronger noise suppression capability, and can still effectively position the target with high precision under the condition of lower signal-to-noise ratio.
Example two: analysis of dual target resolution
The parameter setting is the same as in the first example, two narrowband source positions are set to be (250,245) m and (250,255) m, and the positioning performance and the resolution effect of the three direct positioning algorithms on two sources which are very close to each other are analyzed.
FIGS. 6 and 7 show the spatial spectrum pseudo-color and three-dimensional plots of the dual-signal-source resolution effect of the three direct positioning algorithms when the signal-to-noise ratio is 0dB, respectively;
through comparative analysis, when the two information sources are 10m apart, under the same condition, the space spectrum peaks corresponding to the two information sources are partially overlapped by the scalar MUSIC direct positioning algorithm, and the existence of the two targets is difficult to identify; the vector MUSIC direct positioning algorithm provided by the invention has higher spatial resolution, two adjacent targets are completely separated, and the positions of the two targets can be obviously distinguished in the pseudo-color image and the three-dimensional image.
Example three: performance analysis of different signal-to-noise ratio methods
The parameter setting is the same as that in the first example, the signal-to-noise ratio is increased from-30 dB to 20dB, 50 Monte Carlo experiments are carried out at each signal-to-noise ratio, and the positioning accuracy of the two direct positioning algorithms under different signal-to-noise ratios is analyzed and compared. Fig. 8 shows the positioning accuracy of two direct positioning algorithms as a function of the signal to noise ratio.
It can be clearly seen from the figure that when the signal-to-noise ratio is increased from-30 dB to 20dB, the positioning accuracy of the vector MUSIC direct positioning is better, the positioning error is lower, especially under the low signal-to-noise ratio condition below-10 dB, the positioning performance of the vector processing on the target is obviously improved compared with a scalar, that is, the vector MUSIC direct positioning algorithm has better adaptability to the target positioning under the low signal-to-noise ratio condition.
In summary, the invention discloses a direct positioning method based on an acoustic vector array network, and belongs to the field of target passive positioning. Firstly, carrying out segmented frequency domain transformation on sound pressure and vibration speed data received by each array node; then scanning point by point in the area to be measured, and calculating the guiding direction and time delay at the scanning point; weighting and summing the sound pressure and vibration velocity data to obtain weighted data matrixes of all the arrays; calculating equivalent array manifold by using a frequency domain positioning model; then, calculating the spatial spectrum on each frequency by using a spatial spectrum estimation algorithm, and accumulating to obtain the total spatial spectrum of all frequencies; and calculating the spatial spectrum of each point in the region to be measured according to the proper scanning step length, and searching the peak value of the spatial spectrum value to obtain the position of the target to be measured. The invention introduces the acoustic vector array into the multi-array target direct positioning method, and utilizes the sound pressure and vibration velocity combined processing technology to obtain higher positioning precision and spatial resolution than a scalar array, and has better adaptability to target positioning under the condition of low signal-to-noise ratio.
Claims (4)
1. A direct positioning method based on an acoustic vector array network is characterized in that: the node of the acoustic vector array network is a multi-array element acoustic vector sensor array, each array element can output the same-point sound pressure scalar and the same-point vibration velocity vector information, and the method comprises the following steps:
step 1, establishing a narrow-band signal array positioning model of the array network, and obtaining sound pressure output P of the ith array node l (t) sum velocity vector output V xl (t)，V yl (t)；
Step 2, dividing the sound pressure and vibration velocity data into J sections, and respectively carrying out frequency domain transformation on each section of data to obtain Fourier transformation P of the jth section of observation signal of the ith array in the kth frequency component l (k，j)，V xl (k，j)，V yl (k，j)；
Step 3, scanning point by point in the area to be measured, and calculating the guiding azimuth psi from the point to the I-th array reference point at the scanning point d l (d) And time delay tau l (d)；
Step 4, using the guide orientation psi l (d) Sound pressure P to the same segment l (k, j) and the vibration velocity V xl (k，j)，V yl (k, j) combined treatment to give Y l (k, j), summing Y for all L arrays l (k, j) obtaining a total data matrix Y (k, j);
by making use of psi l (d) Sound pressure P to the same segment l (k, j) and the vibration velocity V xl (k，j)，V yl (k, j) summing to obtain the vibration velocity of the j-th section sound pressure and the Fourier transform Y of the k-th frequency component l (k，j)：
Y l (k，j)＝P l (k，j)+V xl (k，j)cosψ l (d)+V yl (k，j)sinψ l (d)
Summing Y of all L arrays l (k, j) obtaining a total array data matrix Y (k, j) ═ Y 1 (k，j) T ，…，Y L (k，j) T ] T ；
Step 5, positioning the frequency domain signal model by using the array, defining and calculating and the likeArray manifold
Summing all L arrays
Obtaining the total equivalent array manifold matrix
Step 5-1, the array l in the array network outputs Y l The frequency domain model of (k, j) is:
in the formula, S (k, j), N l (k, j) are frequency domain representations of the source and noise, respectively;
step 5-2, defining and calculating the equivalent array manifold according to the following formula
Summing all L arrays
Is a total equivalent array manifold matrix
Step 5-3, order
Due to b l For unknown complex scalar quantities, | | b | | | 1, b ═ b | > [ b | 1 ，…，b L ] T Then the array positioning frequency domain signal model can be expressed as:
the frequency domain signal model is consistent with the array signal model in form;
step 6, extracting the fourier transform Y (k, J) of the J-segment data under the k-th frequency component, wherein J is 1, … and J to form Y (k), and calculating the frequency f by using a spatial spectrum estimation algorithm k Upper spatial spectral values Q (k, d);
step 7, obtaining a total space spectrum Q (d) of accumulated K frequencies;
and 8, setting a proper scanning step length, repeating the steps 3 to 7, traversing the whole region to be detected, outputting a spatial spectrogram, and searching a spectral value peak value to obtain the position of the target to be detected.
2. The direct positioning method based on the acoustic vector array network as claimed in claim 1, wherein: step 6 relates to spatial spectrum estimation methods including CBF, MVDR or MUSIC.
3. The direct positioning method based on the acoustic vector array network as claimed in claim 1, wherein: step 6 relates to a spatial spectrum estimation method, and the MUSIC spatial spectrum estimation algorithm used in the spatial spectrum estimation method specifically comprises the following steps:
step 6-1, extracting fourier transform Y (k, J), J being 1, …, J, constituting Y (k), Y (k) ═ Y (k,1), …, Y (k, J) of J-segment data in k-th frequency component]By the use of R k ＝Y(k)·Y H (k) Calculating frequency f k Upper covariance matrix, R k Is an LM multiplied by LM dimensional covariance matrix;
step 6-2, to covariance matrix R k Performing characteristic decomposition to obtain a signal subspace U consisting of characteristic vectors corresponding to larger D characteristic values S (k) And a noise subspace U consisting of eigenvectors corresponding to other smaller eigenvalues N (k) Wherein, in the step (A),d is the target number;
step 6-3, at frequency f k The MUSIC spatial spectrum above is:
in the formula (I), the compound is shown in the specification,
I L is a L × L dimensional unit array, J M Is an M multiplied by 1 dimensional full 1 array,
representing the Kronecker product, the above formula is equivalent to:
order to
Since | | | b | | | 1, G (k, d) is a quadratic form under the constraint, whose eigenvalue is the maximum of the eigenvalues of G (k, d), Q (k, d) is:
Q(k，d)＝max(λ max (G(k，d)))。
4. a direct positioning method based on acoustic vector array network according to claim 2 or 3, characterized in that: step 7 sums all frequency-space spectra Q (k, d), the total space spectrum Q (d) being:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010197246.8A CN111323750B (en) | 2020-03-19 | 2020-03-19 | Direct positioning method based on acoustic vector array network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010197246.8A CN111323750B (en) | 2020-03-19 | 2020-03-19 | Direct positioning method based on acoustic vector array network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111323750A CN111323750A (en) | 2020-06-23 |
| CN111323750B true CN111323750B (en) | 2022-08-02 |
Family
ID=71171826
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010197246.8A Active CN111323750B (en) | 2020-03-19 | 2020-03-19 | Direct positioning method based on acoustic vector array network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111323750B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112578335B (en) * | 2020-11-29 | 2023-05-02 | 西北工业大学 | Direct positioning and source number estimation method for strong and weak signals |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166120A (en) * | 2014-07-04 | 2014-11-26 | 哈尔滨工程大学 | Acoustic vector circular matrix steady broadband MVDR orientation estimation method |
| JP2015079080A (en) * | 2013-10-16 | 2015-04-23 | 日本電信電話株式会社 | Sound source position estimation device, method, and program |
| CN106066468A (en) * | 2016-05-25 | 2016-11-02 | 哈尔滨工程大学 | A kind of based on acoustic pressure, the vector array port/starboard discrimination method of vibration velocity Mutual spectrum |
| CN106199505A (en) * | 2016-06-28 | 2016-12-07 | 哈尔滨工程大学 | The sane direction estimation method in a kind of acoustic vector circle battle array mode territory |
| CN107290717A (en) * | 2017-05-19 | 2017-10-24 | 中国人民解放军信息工程大学 | For the direct localization method of multiple target of not rounded signal |
| US9829565B1 (en) * | 2016-02-19 | 2017-11-28 | The United States Of America As Represneted By The Secretary Of The Navy | Underwater acoustic beacon location system |
| US10042038B1 (en) * | 2015-09-01 | 2018-08-07 | Digimarc Corporation | Mobile devices and methods employing acoustic vector sensors |
| JP2018146948A (en) * | 2017-03-03 | 2018-09-20 | パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカＰａｎａｓｏｎｉｃ Ｉｎｔｅｌｌｅｃｔｕａｌ Ｐｒｏｐｅｒｔｙ Ｃｏｒｐｏｒａｔｉｏｎ ｏｆ Ａｍｅｒｉｃａ | Sound source probing apparatus, sound source probing method, and program therefor |
| CN109283492A (en) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | Multi-target DOA estimation method and underwater sound vertical vector array system |
| CN109581291A (en) * | 2018-12-11 | 2019-04-05 | 哈尔滨工程大学 | A kind of direct localization method based on artificial bee colony |
| CN110501669A (en) * | 2019-09-25 | 2019-11-26 | 哈尔滨工程大学 | A kind of quick spatial spectrum of central symmetry acoustic vector circle battle array compresses super-resolution direction estimation method |
| CN110764055A (en) * | 2019-10-25 | 2020-02-07 | 哈尔滨工程大学 | Virtual plane array underwater moving target radiation noise vector measurement system and measurement method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9264809B2 (en) * | 2014-05-22 | 2016-02-16 | The United States Of America As Represented By The Secretary Of The Navy | Multitask learning method for broadband source-location mapping of acoustic sources |
| CN107132503B (en) * | 2017-03-23 | 2019-09-27 | 哈尔滨工程大学 | Acoustic vector circle battle array broadband coherent source direction estimation method based on vector singular value decomposition |
| CN110082712B (en) * | 2019-03-14 | 2022-12-13 | 哈尔滨工程大学 | Acoustic vector circular array coherent target azimuth estimation method |
-
2020
- 2020-03-19 CN CN202010197246.8A patent/CN111323750B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015079080A (en) * | 2013-10-16 | 2015-04-23 | 日本電信電話株式会社 | Sound source position estimation device, method, and program |
| CN104166120A (en) * | 2014-07-04 | 2014-11-26 | 哈尔滨工程大学 | Acoustic vector circular matrix steady broadband MVDR orientation estimation method |
| US10042038B1 (en) * | 2015-09-01 | 2018-08-07 | Digimarc Corporation | Mobile devices and methods employing acoustic vector sensors |
| US9829565B1 (en) * | 2016-02-19 | 2017-11-28 | The United States Of America As Represneted By The Secretary Of The Navy | Underwater acoustic beacon location system |
| CN106066468A (en) * | 2016-05-25 | 2016-11-02 | 哈尔滨工程大学 | A kind of based on acoustic pressure, the vector array port/starboard discrimination method of vibration velocity Mutual spectrum |
| CN106199505A (en) * | 2016-06-28 | 2016-12-07 | 哈尔滨工程大学 | The sane direction estimation method in a kind of acoustic vector circle battle array mode territory |
| JP2018146948A (en) * | 2017-03-03 | 2018-09-20 | パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカＰａｎａｓｏｎｉｃ Ｉｎｔｅｌｌｅｃｔｕａｌ Ｐｒｏｐｅｒｔｙ Ｃｏｒｐｏｒａｔｉｏｎ ｏｆ Ａｍｅｒｉｃａ | Sound source probing apparatus, sound source probing method, and program therefor |
| CN107290717A (en) * | 2017-05-19 | 2017-10-24 | 中国人民解放军信息工程大学 | For the direct localization method of multiple target of not rounded signal |
| CN109283492A (en) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | Multi-target DOA estimation method and underwater sound vertical vector array system |
| CN109581291A (en) * | 2018-12-11 | 2019-04-05 | 哈尔滨工程大学 | A kind of direct localization method based on artificial bee colony |
| CN110501669A (en) * | 2019-09-25 | 2019-11-26 | 哈尔滨工程大学 | A kind of quick spatial spectrum of central symmetry acoustic vector circle battle array compresses super-resolution direction estimation method |
| CN110764055A (en) * | 2019-10-25 | 2020-02-07 | 哈尔滨工程大学 | Virtual plane array underwater moving target radiation noise vector measurement system and measurement method |
Non-Patent Citations (6)
| Title |
|---|
| DEMON feature extraction of acoustic vector signal based on 3/2-D spectrum;Sichun, Li等;《ICIEA 2007: 2ND IEEE CONFERENCE ON INDUSTRIAL ELECTRONICS AND APPLICATIONS, VOLS 1-4, PROCEEDINGS》;20071231;2239-2243 * |
| Direct position determination(DPD)of multiple known and unknown radio-frequency signals;Alon Amar等;《2004 12th European Signal Processing Conference》;20150406;1115-1118 * |
| Sequential Maximum-Likelihood Source Localization of a Near Field Emitter of Unknown Spectrum, Using an Acoustic Vector Sensor;Song, Yang等;《2016 IEEE INTERNATIONAL CONFERENCE ON CONSUMER ELECTRONICS-CHINA (ICCE-CHINA)》;20161231;全文 * |
| 中心对称声矢量圆阵的相干双声源方位估计方法;时胜国等;《哈尔滨工程大学学报》;20190731;第40卷(第7期);1187-1193 * |
| 基于广义逆波束形成的噪声源定位识别方法研究;田德艳;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20190115;C036-31 * |
| 基于矢量水听器的声源定位算法研究;韩金金;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20190115;C028-242 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111323750A (en) | 2020-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108375763B (en) | Frequency division positioning method applied to multi-sound-source environment | |
| CN109188344B (en) | Estimation method for source number and incoming wave direction angle based on mutual cyclic correlation MUSIC algorithm in impulse noise environment | |
| CN111123192B (en) | Two-dimensional DOA positioning method based on circular array and virtual extension | |
| CN111798869B (en) | Sound source positioning method based on double microphone arrays | |
| Bosse et al. | A spatio-temporal array processing for passive localization of radio transmitters | |
| CN111413666A (en) | Design method of array frequency and direction measurement combined receiver | |
| CN108398659B (en) | Direction-of-arrival estimation method combining matrix beam and root finding MUSIC | |
| Cao et al. | Direction of arrival estimation for monostatic multiple-input multiple-output radar with arbitrary array structures | |
| CN109696657A (en) | A kind of coherent sound sources localization method based on vector hydrophone | |
| CN108828586B (en) | Bistatic MIMO radar angle measurement optimization method based on beam domain | |
| CN111323750B (en) | Direct positioning method based on acoustic vector array network | |
| CN111273269A (en) | IPSO-BP-based radar target positioning method of frequency diversity array | |
| CN108594165B (en) | Narrow-band signal direction-of-arrival estimation method based on expectation maximization algorithm | |
| CN102062851A (en) | Direction finding method based on improved L array star-carrying broadband multipurpose | |
| CN115826004B (en) | Three-star cooperative direct positioning method based on two-dimensional angle and time difference combination | |
| CN114563760B (en) | Second-order super-beam forming method, equipment and medium based on SCA array | |
| CN114755628A (en) | Method for estimating direction of arrival of acoustic vector sensor array under non-uniform noise | |
| CN114167346B (en) | DOA estimation method and system based on covariance matrix fitting array element expansion | |
| CN114740434A (en) | Equidistance distribution subarray system and method for resisting accompanying interference | |
| Wan et al. | A range-Doppler-angle estimation method for passive bistatic radar | |
| CN114371441A (en) | Virtual array direction of arrival estimation method, device, product and storage medium | |
| CN114325560A (en) | Super-resolution target direction finding method for beam scanning radar | |
| CN112698263A (en) | Orthogonal propagation operator-based single-basis co-prime MIMO array DOA estimation algorithm | |
| CN113589223A (en) | Direction finding method based on nested array under mutual coupling condition | |
| CN112649787B (en) | Target azimuth estimation method based on low-frequency circular array |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |