CN110703198A - Quaternary cross array envelope spectrum estimation method based on frequency selection - Google Patents
Quaternary cross array envelope spectrum estimation method based on frequency selection Download PDFInfo
- Publication number
- CN110703198A CN110703198A CN201911005242.9A CN201911005242A CN110703198A CN 110703198 A CN110703198 A CN 110703198A CN 201911005242 A CN201911005242 A CN 201911005242A CN 110703198 A CN110703198 A CN 110703198A
- Authority
- CN
- China
- Prior art keywords
- cross
- quaternary
- spectrum
- frequency
- quaternary cross
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 26
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 abstract description 5
- 238000000605 extraction Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 2
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01S5/20—Position of source determined by a plurality of spaced direction-finders
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention provides a quaternary cross array envelope spectrum estimation method based on frequency selection, solves the problem of inaccurate envelope spectrum extraction of ship radiation noise under an interference condition, and belongs to the field of ship target signal parameter estimation. The method is realized based on a quaternary cross array and a compass, the quaternary cross array is distributed at an underwater position, the compass is arranged right below the quaternary cross array, the north direction of the compass points to a first array element of the quaternary cross array from the center, the compass is rigidly connected with the quaternary cross array, the course of the quaternary cross array is monitored in real time, and two cross-power spectrums are respectively corrected according to the horizontal steering angle of the quaternary cross array measured by the compass, so that the high-precision azimuth estimation of long-time observation signals is realized. And carrying out azimuth estimation on the signal by using the quaternary cross array, and estimating the frequency component of the signal by using the histogram statistical result of the azimuth, so that the frequency component of the signal is separated from the frequency components of other interference, and the frequency component of the target signal for envelope spectrum estimation can be accurately obtained.
Description
Technical Field
The invention mainly relates to a quaternary cross-shaped array envelope spectrum estimation method based on frequency selection, and belongs to the field of ship target signal parameter estimation.
Background
The envelope spectrum characteristics are the main basis for ship target identification. The method for extracting the ship envelope spectrum features mainly comprises an absolute value method, a flat method and the like. The performance of the extraction of the envelope spectrum features mainly depends on the signal-to-noise ratio of the radiation noise of the ship target in the analysis frequency band. Therefore, how to optimize the frequency band of the ship radiation noise signal, and designing a corresponding filter to obtain a target signal with high signal-to-noise ratio for envelope spectrum analysis becomes a key for obtaining a clear envelope spectrum. In addition, when the ship radiation noise exists in other directions and is interfered, the envelope spectrum feature extraction performance is sharply reduced or is difficult to correspond to the target.
Disclosure of Invention
The invention provides a quaternary cross array envelope spectrum estimation method based on frequency selection, aiming at solving the problem of inaccurate envelope spectrum extraction of ship radiation noise under the interference condition.
The invention relates to a frequency selection-based quaternary cross array envelope spectrum estimation method, which is realized based on a quaternary cross array and a compass, wherein the quaternary cross array is arranged at an underwater position, the compass is arranged right below the quaternary cross array, the north direction of the compass points to a first array element of the quaternary cross array from the center, the compass is in rigid connection with the quaternary cross array, and four hydrophones of the quaternary cross array are used for synchronously acquiring acoustic signals of a target, and the method comprises the following steps:
s1, the four hydrophones receive acoustic signals, signals output by four channels of the quaternary cross matrix are converted into digital signals, the digital signals of the four channels are divided into M sections, and the length of each section is N;
s2, windowing the digital signals of the four channels, and performing discrete Fourier transform to obtain frequency spectrum values of the four channels;
s3, performing cross-spectrum on the frequency spectrum values of two channels on the same coordinate axis in the quaternary cross matrix in the frequency spectrum values of the four channels to obtain two cross-power spectrums;
s4, respectively correcting the two cross-power spectrums according to the horizontal steering angle of the quaternary cross matrix measured by the compass to obtain corrected cross-power spectrums in the north direction and the east direction;
s5, respectively carrying out integral summation on the north-oriented mutual power spectrum and the east-oriented mutual power spectrum of the M sections to respectively obtain a north-oriented mutual power spectrum and a east-oriented mutual power spectrum after integral summation;
s6, obtaining a north-oriented cross-power spectrum and an east-oriented cross-power spectrum according to the integral summation of S5, calculating the orientation under each frequency value, adding the powers corresponding to the same orientation to obtain the power spectrum of each orientation, selecting the orientation corresponding to the maximum value in the power spectrum of each orientation as a target orientation, and extracting the frequency spectrum of the target orientation;
s7, designing a finite impulse response filter according to the frequency spectrum of the target azimuth;
s8, the signal received by a hydrophone of the four-element cross array is x1(n) filtering the signal x by using the filter designed by S7fir(n) taking the absolute value, the result | x of the absolute value will be takenfir(n) passing through a low-pass filter to obtain a low-frequency component xLPF(n) converting the low frequency component xLPF(n) Fourier transform to obtain a power spectrum Xθ(f)；
S9 according to the power spectrum Xθ(f) Obtaining the envelope spectrum S of the target signal in the theta directionθ(f)。
Preferably, the corrected true north direction cross power spectrum P 'in S4'N(k) And a cross-power spectrum P 'in the east direction'E(k)：
Wherein, P11(k) Representing the self-power spectrum of any channel signal in the quaternary cross matrix, j representing the imaginary part,
is P13(k) Phase of (P)13(k) Cross-power spectra of the quaternary cross-matrix channel 1 and channel 3 signals are shown,
is P24(k) Phase of (P)24(k) Cross power spectra of the quaternary cross-matrix 2 nd and 4 th channel signals are shown, and k represents a frequency value.
Preferably, in S2, the digital signals of the four channels are windowed to obtain:
x1(n)＝s1(n)w(n)，
x2(n)＝s2(n)w(n)，
x3(n)＝s3(n)w(n)，
x4(n)＝s4(n)w(n)，
s1(n)、s2(n)、s3(n) and s4(n) represent the digital signals of four channels, respectively, wherein w (n) is a window function, and the window function is a rectangular window, a hanning window or a hamming window.
Preferably, in S3, two cross-power spectra are obtained, respectively:
in the formula, X1(k) Fourier transformed spectral values, X, of channel 1 signals representing a quaternary cross2(k) Representing the spectral values, X, of the quaternary cross-matrix 2 nd channel signal after Fourier transformation3(k) Representing the spectral values, X, of the quaternary cross-matrix 3 rd channel signal after Fourier transformation4(k) And the frequency spectrum value of the quaternary cross 4 th channel signal after Fourier transformation is represented, k represents a frequency value, and H represents complex conjugate.
As a preference, the first and second liquid crystal compositions are,
wherein,
n denotes the signal length and f denotes the frequency variation.
The method has the advantages that the relative orientation of the target under the matrix coordinate system is estimated by using the quaternary cross array, the compass is rigidly connected with the quaternary cross array, and the course of the quaternary cross array can be monitored in real time, so that the high-precision orientation estimation of long-time observation signals is realized. In order to achieve the purpose of the invention, a quaternary cross matrix is required to carry out azimuth estimation on the signal, and the frequency component of the signal is estimated by utilizing the histogram statistical result of the azimuth, so that the frequency component of the target signal for envelope spectrum estimation can be accurately obtained by separating the frequency component from the frequency components of other interference.
Drawings
Fig. 1 is a diagram of a four-element cross matrix.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The invention is further described with reference to the following drawings and specific examples, which are not intended to be limiting.
The quaternary cross array of the present embodiment is placed at a depth of about 20 meters under water. As shown in fig. 1, the quaternary cross array is composed of 4 hydrophones in the same plane, a left-hand rectangular coordinate system is defined in the plane by using the center of the quaternary cross array as an origin, and the positions of four array elements are respectively: array element 1(α,0), array element 2(0, α), array element 3(- α,0), array element 4(0, - α), 2a is the distance of the diagonal array elements. A compass is arranged right below the quaternary cross array, the north direction of the compass points to the array element 1 from the center, and the compass is rigidly connected with the quaternary cross array. Four hydrophones of the quaternary cross array are used for synchronously acquiring acoustic signals of a target,
the frequency selection-based quaternary cross-shaped envelope spectrum estimation method of the embodiment comprises the following steps:
s1, the four hydrophones receive acoustic signals of a target, and the signals output by the four channels of the quaternary cross array are converted into digital signals; dividing the digital signal of each channel into M sections, wherein the length of each section is N, and M and N are positive integers;
s2, windowing the digital signals of the four channels, and performing discrete Fourier transform, wherein the number of Fourier transform points is determined by the analysis bandwidth, so as to obtain the frequency spectrum values of the four channels;
s3, compensating the cross spectrum on the two axes to the north and east directions by combining compass information:
performing cross-spectrum on the frequency spectrum values of two channels on the same coordinate axis in the quaternary cross matrix in the frequency spectrum values of the four channels to obtain two cross-power spectrums;
s4, respectively correcting the two cross-power spectrums according to the horizontal steering angle of the quaternary cross matrix measured by the compass to obtain corrected cross-power spectrums in the north direction and the east direction;
s5, accumulating the compensated cross spectra of the M sections at different time:
respectively carrying out integral summation on the north-oriented mutual power spectrum and the east-oriented mutual power spectrum of the M sections to respectively obtain a north-oriented mutual power spectrum and a east-oriented mutual power spectrum after integral summation;
s6, obtaining a north-oriented cross power spectrum and a east-oriented cross power spectrum according to the integral summation of S5, calculating the orientation of each frequency value, adding the powers corresponding to the same orientation to target signals on different frequency points, performing histogram statistics to obtain the power spectrum of each orientation, selecting the orientation corresponding to the maximum value in the power spectrum of each orientation as the target orientation, extracting the frequency spectrum of the target orientation, and obtaining a signal passband formed by a plurality of sub-bands;
s7, designing a finite impulse response filter by using a window function method according to a signal passband;
s8, the signal received by a hydrophone of the four-element cross array is x1(n) filtering the signal x by using the filter designed by S7fir(n) taking the absolute value, the result | x of the absolute value will be takenfir(n) passing through a low-pass filter to obtain a low-frequency component xLPF(n) converting the low frequency component xLPF(n) Fourier transform to obtain a power spectrum Xθ(f)；
S9 according to the power spectrum Xθ(f) Obtaining the envelope spectrum S of the target signal in the theta directionθ(f)。
In the embodiment, the power spectrum of the signal after the absolute value is calculated by adopting a Welch method, and the line spectrum of the low-frequency band lower than 50Hz is taken to be judged as the envelope spectrum of the target signal.
The specific calculation process of the embodiment is as follows:
the received signals of the four hydrophones are sampled by an analog-to-digital converter to obtain s1(n)，s2(n)，s3(n)，s4And (n) four paths of signals, wherein n is the serial number of the sample point. The signal is windowed to obtain
x1(n)＝s1(n)w(n) (1)
x2(n)＝s2(n)w(n) (2)
x3(n)＝s3(n)w(n) (3)
x4(n)＝s4(n)w(n) (4)
w (n) is a window function, which may be selected from a rectangular window, a hanning window, a hamming window, and the like. The rectangular window is defined as:
rectangular window:
hanning Window:
hamming window:
n is the length of the window function. The signal length after windowing is also N.
Taking discrete Fourier transform:
X1(k) fourier transformed spectral values, X, of channel 1 signals representing a quaternary cross2(k) Fourier transformed spectral values, X, of a channel 2 signal representing a quaternary cross3(k) Fourier transformed spectral values, X, of a channel 3 signal representing a quaternary cross4(k) And the 4 th channel signal of the quaternary cross matrix is represented as a frequency spectrum value after Fourier transformation, and k represents a frequency value.
Cross-spectra were found in the analysis band:
wherein H represents a complex conjugate. P13(k) Representing cross-power spectra, P, of channel 1 and channel 3 signals24(k) Representing the cross-power spectra of the channel 2 and channel 4 signals.
Combining the horizontal steering angle theta of the matrix obtained by compass measurementcAnd correcting the cross spectrum.
P 'here'N(k) Represents the rotated true north cross-spectrum, where P11(k) Representing the self-power spectrum of the channel 1 signal, j is an imaginary symbol,
is P13(k) The phase of (a) is determined,
is P24(k) Phase of (1), P'E(k) Is the cross-spectrum of the rotated righteast direction.
The signal observation time comprises M segments of signals with length of N, and the compass corrected cross-spectrum obtained from each segment of time signal
And
integral summation is carried out to obtain the north-oriented cross spectrum
Cross spectrum of the east
Showing the corrected north-oriented cross spectrum of the compass obtained from the mth segment of signal,
the corrected righteast cross-spectrum of the compass obtained from the mth signal is shown.
Using cross spectra after integration
Andfinding the bearing at each frequency value k:
here, Arg (-) denotes the argument of the complex number.Representing a complex number.
And
are respectively
And
the argument of (2).
S(θ)＝∑P11(k) (16)
S (theta) represents a power spectrum in theta direction, the power spectrums in all directions form a target signal space spectrum characteristic, and the azimuth theta corresponding to the maximum value of S (theta) is taken as a target azimuth thetaT
Extracting target azimuth theta from histogramTSpectrum of
According to frequency response
And designing the FIR filter with the order of N by adopting a window function method to obtain a filter response of h (N).
As an example, it is not assumed that one hydrophone of a four-element cross array receives a signal of x1(n) is represented by
x1(n)＝C(1+sin 2πf2n)cos 2πf1n(19)
Where C is the signal amplitude, f2Is the frequency of the target signal in the theta direction, f1The frequencies of the signals in the other directions.
Will signal x1(n) filtering through the constructed FIR filter to obtain a filtered result xfir(n)
Filtered result xfirAnd (n) is a target signal in the theta direction. Wherein,representing a convolution, the filtering process can be viewed as a signal x1Convolution of (n) with the filter response h (n), in particular
xfir(n) is the target signal in the theta direction and the noise in that direction, which can be expressed as
xfir(n)＝C(1+sin2πf2n)cos 2πf1n(22)
Where C is the signal amplitude, f2Is the frequency of the target signal in the theta direction, f1Is the frequency of the multiplicative noise in that direction.
Filtering the filtered signal xfir(n) taking the absolute value to obtain | xfir(n)|，
Calculating absolute value of the result | xfir(n) passing through a low-pass filter to obtain a low-frequency component xLPF(n)，
To the filtering result xLPF(n) Fourier transform
Calculating Xθ(f) And power spectrum, obtaining a target signal envelope spectrum in the theta direction:
although the invention herein has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present invention. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims. It should be understood that features described in different dependent claims and herein may be combined in ways different from those described in the original claims. It is also to be understood that features described in connection with individual embodiments may be used in other described embodiments.
Claims (5)
1. The quaternary cross array envelope spectrum estimation method based on frequency selection is characterized by being realized based on a quaternary cross array and a compass, the quaternary cross array is arranged at an underwater position, the compass is arranged right below the quaternary cross array, the north of the compass points to a first array element of the quaternary cross array from the center, the compass is in rigid connection with the quaternary cross array, and four hydrophones of the quaternary cross array are used for synchronously acquiring acoustic signals of a target, and the method comprises the following steps:
s1, the four hydrophones receive acoustic signals, signals output by four channels of the quaternary cross matrix are converted into digital signals, the digital signals of the four channels are divided into M sections, and the length of each section is N;
s2, windowing the digital signals of the four channels, and performing discrete Fourier transform to obtain frequency spectrum values of the four channels;
s3, performing cross-spectrum on the frequency spectrum values of two channels on the same coordinate axis in the quaternary cross matrix in the frequency spectrum values of the four channels to obtain two cross-power spectrums;
s4, measuring the horizontal steering angle theta of the quaternary cross matrix according to the compasscRespectively correcting the two cross power spectrums to obtain a corrected cross power spectrum in the north direction and a corrected cross power spectrum in the east direction;
s5, respectively carrying out integral summation on the north-oriented mutual power spectrum and the east-oriented mutual power spectrum of the M sections to respectively obtain a north-oriented mutual power spectrum and a east-oriented mutual power spectrum after integral summation;
s6, obtaining a north-oriented cross-power spectrum and an east-oriented cross-power spectrum according to the integral summation of S5, calculating the orientation under each frequency value, adding the powers corresponding to the same orientation to obtain the power spectrum of each orientation, selecting the orientation corresponding to the maximum value in the power spectrum of each orientation as a target orientation, and extracting the frequency spectrum of the target orientation;
s7, designing a finite impulse response filter according to the frequency spectrum of the target azimuth;
s8, the signal received by a hydrophone of the four-element cross array is x1(n) filtering the signal x by using the filter designed by S7fir(n) taking the absolute value, the result | x of the absolute value will be takenfir(n)]Obtaining low-frequency component x by low-pass filterLPF(n) converting the low frequency component xLPF(n) Fourier transform to obtain a power spectrum Xθ(f)；
S9 according to the power spectrum Xθ(f) Obtaining the envelope spectrum S of the target signal in the theta directionθ(f)。
2. The frequency-selection-based quaternary cross-matrix envelope spectrum estimation method of claim 1, wherein the corrected true north direction cross-power spectrum P 'in S4'N(k) And a cross-power spectrum P 'in the east direction'E(k)：
Wherein, P11(k) Representing the self-power spectrum of any channel signal in the quaternary cross matrix, j representing the imaginary part,
is P13(k) Phase of (P)13(k) Cross-power spectra of the quaternary cross-matrix channel 1 and channel 3 signals are shown,
is P24(k) Phase of (P)24(k) Cross power spectra of the quaternary cross-matrix 2 nd and 4 th channel signals are shown, and k represents a frequency value.
3. The method for estimating the quaternary cross-matrix envelope spectrum based on frequency selection according to claim 1 or 2, wherein in S2, the digital signals of four channels are windowed to obtain:
x1(n)＝s1(n)w(n)，
x2(n)＝s2(n)w(n)，
x3(n)＝s3(n)w(n)，
x4(n)＝s4(n)w(n)，
s1(n)、s2(n)、s3(n) and s4(n) represent the digital signals of four channels, respectively, wherein w (n) is a window function, and the window function is a rectangular window, a hanning window or a hamming window.
4. The method for estimating the quaternary cross-matrix envelope spectrum based on frequency selection according to claim 3, wherein in step S3, two cross-power spectrums are obtained as follows:
in the formula, X1(k) Fourier transformed spectral values, X, of channel 1 signals representing a quaternary cross2(k) Representing the spectral values, X, of the quaternary cross-matrix 2 nd channel signal after Fourier transformation3(k) Representing the spectral values, X, of the quaternary cross-matrix 3 rd channel signal after Fourier transformation4(k) And the frequency spectrum value of the quaternary cross 4 th channel signal after Fourier transformation is represented, k represents a frequency value, and H represents complex conjugate.
5. The method of estimating quaternary cross-shaped envelope spectrum based on frequency selection according to claim 4,
wherein,
n denotes the signal length and f denotes the frequency variation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911005242.9A CN110703198B (en) | 2019-10-22 | 2019-10-22 | Quaternary cross array envelope spectrum estimation method based on frequency selection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911005242.9A CN110703198B (en) | 2019-10-22 | 2019-10-22 | Quaternary cross array envelope spectrum estimation method based on frequency selection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110703198A true CN110703198A (en) | 2020-01-17 |
| CN110703198B CN110703198B (en) | 2022-03-22 |
Family
ID=69200898
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911005242.9A Active CN110703198B (en) | 2019-10-22 | 2019-10-22 | Quaternary cross array envelope spectrum estimation method based on frequency selection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110703198B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114737448A (en) * | 2022-03-22 | 2022-07-12 | 东南大学 | Ground penetrating radar-based intelligent asphalt pavement compaction system and method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005002070A2 (en) * | 2003-03-20 | 2005-01-06 | Bae Systems Information & Electronic Systems Integration Inc. | Correlation interferometer geolocation |
| US9869752B1 (en) * | 2016-04-25 | 2018-01-16 | Ocean Acoustical Services And Instrumentation Systems, Inc. | System and method for autonomous joint detection-classification and tracking of acoustic signals of interest |
| CN109975813A (en) * | 2019-04-12 | 2019-07-05 | 哈尔滨工程大学 | A kind of ten word array lines spectrum target passive detection method based on underwater gliding machine platform |
| CN110007276A (en) * | 2019-04-18 | 2019-07-12 | 太原理工大学 | A kind of sound localization method and system |
-
2019
- 2019-10-22 CN CN201911005242.9A patent/CN110703198B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005002070A2 (en) * | 2003-03-20 | 2005-01-06 | Bae Systems Information & Electronic Systems Integration Inc. | Correlation interferometer geolocation |
| US9869752B1 (en) * | 2016-04-25 | 2018-01-16 | Ocean Acoustical Services And Instrumentation Systems, Inc. | System and method for autonomous joint detection-classification and tracking of acoustic signals of interest |
| CN109975813A (en) * | 2019-04-12 | 2019-07-05 | 哈尔滨工程大学 | A kind of ten word array lines spectrum target passive detection method based on underwater gliding machine platform |
| CN110007276A (en) * | 2019-04-18 | 2019-07-12 | 太原理工大学 | A kind of sound localization method and system |
Non-Patent Citations (4)
| Title |
|---|
| JACEK DMOCHOWSKI 等: "Direction of Arrival Estimation Using the Parameterized Spatial Correlation Matrix", 《IEEE TRANSACTIONS ON AUDIO, SPEECH, AND LANGUAGE PROCESSING》 * |
| JINYAN DU 等: "Attitude Correction Method in DOA Estimation of Single Vector Hydrophone", 《2018 OCEANS - MTS/IEEE KOBE TECHNO-OCEANS (OTO)》 * |
| SHAO-HUA CHEN 等: "Signal detection and bearing estimation for underwater", 《2012 IEEE 11TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING》 * |
| 牛嗣亮 等: "采用姿态修正的单矢量水听器远场测向技术", 《声学技术》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114737448A (en) * | 2022-03-22 | 2022-07-12 | 东南大学 | Ground penetrating radar-based intelligent asphalt pavement compaction system and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110703198B (en) | 2022-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103454495B (en) | Self-adaptation high-precision fast spectral analytical approach | |
| US9949029B2 (en) | Audio filtering with virtual sample rate increases | |
| CN106019214B (en) | Wide-band coherent signal source DOA estimation method | |
| CN109188017B (en) | Wind direction and wind speed detection device and detection method based on ultrasonic sensor array | |
| CN104768100A (en) | Time domain broadband harmonic region beam former and beam forming method for ring array | |
| CN114428225B (en) | Multi-radiation source arrival angle measuring method and device based on quasi-matched filtering | |
| CN107656237A (en) | The method and its device of a kind of multiple source frequency and DOA joint-detections | |
| CN107561357A (en) | A kind of high-precision instantaneous frequency measurement method and apparatus based on channelizing | |
| CN109597021B (en) | Direction-of-arrival estimation method and device | |
| CN105044769B (en) | The method for improving the resolution ratio of seismic signal | |
| CN110703198B (en) | Quaternary cross array envelope spectrum estimation method based on frequency selection | |
| CN113376576B (en) | Positioning method of sound source positioning sensor based on small-aperture microphone array | |
| US20120162004A1 (en) | Apparatus and method for estimating direction of arrival of signal in communication system | |
| CN114578290A (en) | Microphone array sound source positioning method based on frequency correction | |
| CN108919189B (en) | Array signal processing method for frequency and orientation joint estimation | |
| CN110716172A (en) | Vector hydrophone envelope spectrum estimation method based on frequency selection | |
| CN113627398B (en) | Signal characteristic detection method based on adaptive reconstruction filtering | |
| JP5253026B2 (en) | Signal separation device and signal separation method | |
| US11567162B2 (en) | Device and method for estimating direction of arrival | |
| CN115015942A (en) | Self-adaptive underwater target acoustic speed measurement device and method | |
| CN114036975A (en) | Target signal extraction method based on frequency domain-wavenumber domain deconvolution | |
| CN110703199B (en) | Quaternary cross array high-precision azimuth estimation method based on compass compensation | |
| RU2047278C1 (en) | Method for scanning environment by hydroacoustic system and device for implementation of said method | |
| CN116996137B (en) | Low signal-to-noise ratio broadband linear frequency modulation signal detection method based on weighted superposition | |
| CN109061564A (en) | Simplification near field localization method based on Higher Order Cumulants |
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 |