CN112114287B - Outlier real-time eliminating method for azimuth observation data - Google Patents
Outlier real-time eliminating method for azimuth observation data Download PDFInfo
- Publication number
- CN112114287B CN112114287B CN202010994478.6A CN202010994478A CN112114287B CN 112114287 B CN112114287 B CN 112114287B CN 202010994478 A CN202010994478 A CN 202010994478A CN 112114287 B CN112114287 B CN 112114287B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- array
- outliers
- data
- observation data
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000001228 spectrum Methods 0.000 claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 230000008030 elimination Effects 0.000 description 6
- 238000003379 elimination reaction Methods 0.000 description 6
- 238000013213 extrapolation Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003672 processing method Methods 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a method for eliminating outliers of azimuth observation data in real time, which comprises the steps of firstly obtaining broadband array receiving signals of an underwater moving target, estimating the azimuth of the target by calculating the space spectrum of the receiving signals at different moments, obtaining the azimuth observation data, and distinguishing and correcting the outliers of the azimuth observation data by utilizing the uniform motion characteristic of the target. The invention is suitable for eliminating isolated outliers, spot outliers and mixed outliers, has simple method, improves the accuracy of dynamic outlier point detection and azimuth tracking performance of the moving target azimuth measurement data, can be popularized and applied to the discrimination and correction of isolated outliers, spot outliers and mixed outliers of measurement data of other moving target measurement systems, and has good theoretical and engineering application values.
Description
Technical Field
The invention belongs to the field of underwater acoustic signal processing, and particularly relates to a method for removing outliers of underwater moving target azimuth observation data in real time.
Background
Sonar is equipment for underwater detection by using sound waves, and orientation estimation is a basic function of sonar. In the azimuth estimation process, due to the influence of various factors such as noise, interference and ocean multipath propagation channels, abnormal values, i.e. outliers, which are seriously deviated from the normal azimuth change rule, may exist in the azimuth measurement data. The existence of the outlier can seriously affect the performance of subsequent sonar signal processing. The outliers are generally classified into isolated, spot and mixed outliers. For isolated outliers, the widely adopted methods mainly comprise an alpha-beta-gamma filtering method, a least square estimation method, a matching element measurement method and the like, and the elimination effect is good. For speckle outliers, the elimination method mainly comprises a polynomial extrapolation fitting method, a Kalman filtering method, a robust filtering M estimation method, a least square B spline approximation method and the like, and is suitable for outlier elimination under different conditions, wherein only the robust filtering M estimation method has a good elimination effect, but the processing process is complex. For the mixed type outlier with both the isolated outlier and the speckle outlier, the common method is to distinguish the isolated outlier and the speckle outlier first, then to process them separately, the process is complicated, and further research is needed for the accurate determination of the starting point and the ending point of the speckle outlier. In addition, for a moving object, the azimuth changes with time, and it is difficult to effectively correct the speckle pattern outlier by methods such as polynomial extrapolation fitting, and further research and study are also required.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects of the prior art, the invention provides a field value real-time eliminating method of azimuth observation data by distinguishing and correcting the field value of azimuth measurement data through course data by utilizing the characteristic of uniform motion of a target.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the technical scheme that:
a method for eliminating outliers of azimuth observation data in real time comprises the following steps:
step 1, acquiring a broadband array receiving signal of an underwater moving target by using a base array.
And 2, calculating the spatial spectrum of the broadband array received signals at different moments, wherein the spatial spectrum is the energy distribution of all possible incoming wave directions in the space, estimating the target position through the spatial spectrum, and acquiring target position observation data.
And 3, for the uniform motion target or the segmented uniform motion target, taking three adjacent point azimuth measurement data with the time interval T in the uniform motion time period, and calculating the course:
where phi is the heading and alpha is1,α2,α3For adjacent three-point azimuth measurement data of time interval T, theta1=α2-α1,θ2=α3-α2。
Step 4, repeating the step 3 to obtain a group of course dataφiIndicating the heading at the ith time, and n indicating the n time.
Step 6, carrying out outlier discrimination on the course data, and calculating the residual error of the course data at each momentComparing the residual with a 3 sigma threshold, if viIf | ≧ 3 σ, heading data φ corresponding to the timeiAnd the azimuth data is a outlier, the outlier is removed, course data is utilized to correct the outlier, and sigma is the variance of the course data. Other methods may be used to determine the outlier of the heading data.
Further: for the plane wave incident signal in step 1, the frequency domain of the broadband array received signal is represented as:
X(ω)=AS(ω)+N(ω)
where X (ω) is a frequency spectrum of an M × 1 dimensional array received signal, ω is an angular frequency, S (ω) is a frequency spectrum of an N × 1 dimensional spatial source signal, N (ω) is a frequency spectrum of an M × 1 dimensional noise signal, and a is an M × N dimensional array flow pattern matrix:
A=[a0(ω),a1(ω),…,aN-1(ω)]
wherein, aiAnd (omega) is a guide vector of the ith signal, N is the number of information sources, and M is the number of array elements.
For a uniform linear array:
wherein i is 0,1iIs the incident angle of the ith signal, d is the array element spacing, j is the complex factor, c is the acoustic velocity [ ·]TIs a matrix transposition operation.
Further: in step 2, the spatial spectrum can be calculated by methods such as adaptive beam forming (MVDR), high-resolution spatial spectrum estimation (MUSIC) and the like, and the target azimuth estimation can be obtained by searching spatial spectrum peaks.
Further: calculating the spatial spectrum of the array received signals at different moments by adopting a self-adaptive beam forming MVDR method:
wherein, PMVDR(theta) is the spatial spectrum, theta is the scan orientation, a (theta) is the steering vector, [ ·]HFor conjugate transposition, R-1Is the inverse of the received signal covariance matrix. Searching a spatial spectrum peak to obtain a target azimuth, and repeating the processing process on array received signals at different moments to obtain a group of azimuth observation data.
Further: in the step 1, the array is a linear array or a cylindrical array or any other array type.
Compared with the prior art, the invention has the following beneficial effects:
the invention utilizes the uniform motion characteristic of the target, judges and corrects the outlier of the azimuth measurement data through the course data, does not need to distinguish the isolated outlier and the spot outlier, can directly process the mixed outlier, has simple method, improves the detection accuracy and the azimuth tracking performance of the dynamic outlier point of the azimuth measurement data of the moving target, can be popularized and applied to the judgment and correction of the isolated outlier, the spot outlier and the mixed outlier of the measurement data of other moving target measurement systems, and has good theoretical and engineering application values.
Drawings
FIG. 1 shows the elimination of mixed outliers in azimuth measurement data by a polynomial extrapolation fitting method.
FIG. 2 shows the elimination result of the mixed outlier in the azimuth measurement data according to the method of the present invention.
Detailed Description
The present invention is further illustrated by the following description in conjunction with the accompanying drawings and the specific embodiments, it is to be understood that these examples are given solely for the purpose of illustration and are not intended as a definition of the limits of the invention, since various equivalent modifications will occur to those skilled in the art upon reading the present invention and fall within the limits of the appended claims.
A method for removing outliers of azimuth observation data in real time takes outliers of a group of azimuth measurement data to be removed in real time as an example, and the specific implementation process is as follows:
step (1) obtaining a broadband array receiving signal: for plane wave incident signals, the frequency domain of the broadband array received signal is represented as
X(ω)=AS(ω)+N(ω)
Where X (ω) is a frequency spectrum of an M × 1 dimensional array received signal, S (ω) is a frequency spectrum of an N × 1 dimensional spatial source signal, N (ω) is a frequency spectrum of an M × 1 dimensional noise signal, a is an M × N dimensional array flow pattern matrix:
A=[a0(ω),a1(ω),…,aN-1(ω)]
wherein a isiAnd (omega) is a guide vector of the ith signal, N is the number of information sources, and M is the number of array elements. For a uniform linear array of the wires,
wherein i is 0,1iIs the incident angle of the ith signal, d is the array element spacing [. degree]TIs a matrix transposition operation.
The given target signal is band-limited noise, the frequency band is 200 Hz-4000 Hz, the sampling frequency is 48kHz, the target moves at a constant speed, the speed is 20m/s, and the heading is 150 degrees. The receiving array is a uniform linear array, the number of the array elements is 10, and the spacing between the array elements is 0.5 m. The noise is white Gaussian noise with zero mean value, the signal-to-noise ratio is-10 dB, the sampling fast beat number is 2048, and the total observation time is 5.12 s.
Acquiring azimuth observation data: computing spatial spectrum of array received signals at different moments by adopting MVDR algorithm
Where a (θ) is the steering vector. And searching a spatial spectrum peak to obtain target azimuth estimation, and repeating the processing process on array received signals at different moments to obtain a group of azimuth observation data.
Step (3), calculating course: taking the azimuth measurement data of adjacent three points with time interval T, and calculating course by the following formula
Wherein alpha is1,α2,α3For adjacent three-point azimuth measurement data of time interval T, theta1=α2-α1,θ2=α3-α2. Other methods of calculating heading may also be employed.
And (6) carrying out wild value discrimination and correction on the azimuth measurement data: calculating residual error of course data at each momentComparing the residual with a 3 sigma threshold, if viIf | ≧ 3 σ, heading data φ corresponding to the timeiAnd the orientation data is the outlier, the outlier is removed, and the outlier is corrected by utilizing the course data.
FIG. 1 is a diagram illustrating a mixed outlier in azimuth measurement data processed by a polynomial extrapolation fitting method, and FIG. 2 is a diagram illustrating a mixed outlier in azimuth measurement data processed by the method of the present invention. Comparing fig. 1 and fig. 2, it can be seen that, for the mixed outlier, the correction of the outlier point by the polynomial extrapolation fitting method, especially the correction error of the speckle outlier point is larger, and the isolated outlier and the speckle outlier in the outlier can be better distinguished and corrected by the present invention.
The method for eliminating the outlier of the orientation measurement data in real time provides an effective solution to the problems that the existing mixed outlier processing method is complex, the positions of the starting point and the ending point of the spot outlier are difficult to accurately determine, the dynamic spot outlier is difficult to correct and the like, improves the detection accuracy and the orientation tracking performance of the dynamic outlier of the orientation measurement data of the moving target, can be popularized and applied to the discrimination and correction of the isolated outlier, the spot outlier and the mixed outlier of the measurement data of other moving target measurement systems, and has good theoretical and engineering application values.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (5)
1. A method for eliminating outliers of azimuth observation data in real time is characterized by comprising the following steps:
step 1, acquiring a broadband array receiving signal of an underwater moving target by using a base array;
step 2, calculating the spatial spectrum of the broadband array received signal at different moments, wherein the spatial spectrum is the energy distribution of all possible incoming wave directions in the space, estimating the target position through the spatial spectrum, and acquiring target position observation data;
and 3, for the uniform motion target or the segmented uniform motion target, taking three adjacent point azimuth measurement data with the time interval T in the uniform motion time period, and calculating the course:
where phi is the heading and alpha is1,α2,α3For adjacent three-point azimuth measurement data of time interval T, theta1=α2-α1,θ2=α3-α2;
Step 4, repeating the step 3 to obtain a group of course dataφiIndicating the course of the time i, and n indicates the time;
step 6, carrying out outlier discrimination on the course data, and calculating the residual error of the course data at each momentComparing the residual with a 3 sigma threshold, if viIf | ≧ 3 σ, heading data φ corresponding to the timeiAnd the azimuth observation data are outliers, the outliers are removed, course data are used for carrying out outlier correction, and sigma is the variance of the course data.
2. The method for eliminating outliers of the azimuth observation data in real time according to claim 1, which is characterized in that: for the plane wave incident signal in step 1, the frequency domain of the broadband array received signal is represented as:
X(ω)=AS(ω)+N(ω)
where X (ω) is a frequency spectrum of an M × 1 dimensional array received signal, ω is an angular frequency, S (ω) is a frequency spectrum of an N × 1 dimensional spatial source signal, N (ω) is a frequency spectrum of an M × 1 dimensional noise signal, and a is an M × N dimensional array flow pattern matrix:
A=[a0(ω),a1(ω),…,ak(ω),…,aN-1(ω)]
wherein, ak(omega) is a steering vector of a kth signal, N is an information source number, and M is an array element number;
for a uniform linear array:
3. The method for eliminating outliers of the azimuth observation data in real time according to claim 1, which is characterized in that: and 2, calculating the spatial spectrum of the array received signals at different moments by adopting a self-adaptive beam forming MVDR method or a high-resolution spatial spectrum estimation MUSIC method, wherein the position corresponding to the spatial spectrum peak is the target position.
4. The method for eliminating outliers of the azimuth observation data in real time according to claim 1, which is characterized in that: calculating the spatial spectrum of the array received signals at different moments by adopting a self-adaptive beam forming MVDR method:
wherein, PMVDR(theta) is the spatial spectrum, theta is the scan orientation, a (theta) is the steering vector, [ ·]HFor conjugate transposition, R-1Is the inverse of the received signal covariance matrix; searching a spatial spectrum peak to obtain a target azimuth, and repeating the processing process on array received signals at different moments to obtain a group of azimuth observation data.
5. The method for eliminating outliers of the azimuth observation data in real time according to claim 1, which is characterized in that: in the step 1, the array is a linear array or a cylindrical array or any other array type.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010994478.6A CN112114287B (en) | 2020-09-21 | 2020-09-21 | Outlier real-time eliminating method for azimuth observation data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010994478.6A CN112114287B (en) | 2020-09-21 | 2020-09-21 | Outlier real-time eliminating method for azimuth observation data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112114287A CN112114287A (en) | 2020-12-22 |
CN112114287B true CN112114287B (en) | 2021-04-20 |
Family
ID=73801422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010994478.6A Expired - Fee Related CN112114287B (en) | 2020-09-21 | 2020-09-21 | Outlier real-time eliminating method for azimuth observation data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112114287B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108335A (en) * | 2017-12-26 | 2018-06-01 | 北京邮电大学 | A kind of method of abnormal value removing and correction and device |
CN109447133A (en) * | 2018-10-11 | 2019-03-08 | 东南大学 | A kind of elimination method of the location information outlier based on SVR algorithm |
CN111307160A (en) * | 2020-03-20 | 2020-06-19 | 中国人民解放军军事科学院国防科技创新研究院 | Correction method and device for long-endurance flight path of autonomous underwater vehicle |
CN111474536A (en) * | 2020-05-06 | 2020-07-31 | 智慧航海(青岛)科技有限公司 | Intelligent ship autonomous positioning system and method based on shore-based radar system |
CN111679246A (en) * | 2020-06-04 | 2020-09-18 | 哈尔滨工程大学 | Small-scale array high-resolution direction finding method carried by three-dimensional motion platform |
-
2020
- 2020-09-21 CN CN202010994478.6A patent/CN112114287B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108335A (en) * | 2017-12-26 | 2018-06-01 | 北京邮电大学 | A kind of method of abnormal value removing and correction and device |
CN109447133A (en) * | 2018-10-11 | 2019-03-08 | 东南大学 | A kind of elimination method of the location information outlier based on SVR algorithm |
CN111307160A (en) * | 2020-03-20 | 2020-06-19 | 中国人民解放军军事科学院国防科技创新研究院 | Correction method and device for long-endurance flight path of autonomous underwater vehicle |
CN111474536A (en) * | 2020-05-06 | 2020-07-31 | 智慧航海(青岛)科技有限公司 | Intelligent ship autonomous positioning system and method based on shore-based radar system |
CN111679246A (en) * | 2020-06-04 | 2020-09-18 | 哈尔滨工程大学 | Small-scale array high-resolution direction finding method carried by three-dimensional motion platform |
Non-Patent Citations (1)
Title |
---|
一种长基线中野值点判定、剔除及修正方法;李钉云 等;《声学技术》;20200229;第39卷(第1期);第117-120页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112114287A (en) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110764053B (en) | Multi-target passive positioning method based on underwater sensor network | |
CN109799495B (en) | Broadband time delay estimation method for high fidelity array processing | |
CN111025273B (en) | Distortion drag array line spectrum feature enhancement method and system | |
CN111537982B (en) | Distortion drag array line spectrum feature enhancement method and system | |
CN109540172A (en) | A kind of Target moving parameter estimation method for submarine mine platform | |
CN111896913B (en) | High-frequency radar monopole/crossed loop antenna channel gain calibration method and device | |
CN111580079B (en) | Underwater moving target nearest distance estimation method based on single-vector hydrophone | |
CN114089279A (en) | Sound target positioning method based on uniform concentric circle microphone array | |
CN111624573A (en) | Time domain self-adaptive target detection method under sea clutter background | |
CN115494497A (en) | Phase coding pulse radar signal direction finding method | |
CN111175727B (en) | Method for estimating orientation of broadband signal based on conditional wave number spectral density | |
CN108957389A (en) | A kind of real number field multi channel signals method for estimating target azimuth | |
CN109814065B (en) | Beam forming method based on phase factor weighting | |
CN112114287B (en) | Outlier real-time eliminating method for azimuth observation data | |
CN105652256B (en) | A kind of high-frequency ground wave radar TBD methods based on polarization information | |
CN115656994B (en) | Real-time calibration method for bistatic active detection towing array shape | |
CN109541572B (en) | Subspace orientation estimation method based on linear environment noise model | |
CN115201821B (en) | Small target detection method based on strong target imaging cancellation | |
CN111722178B (en) | Far-field narrow-band signal incoming wave direction estimation method based on numerical solution of directivity model | |
CN113687364B (en) | Method for inhibiting false bottom of three-dimensional imaging sonar image | |
CN113075645B (en) | Distorted matrix line spectrum enhancement method based on principal component analysis-density clustering | |
CN112415469B (en) | Rapid interference direction finding method for two-dimensional digital array radar | |
CN110632579B (en) | Iterative beam forming method using subarray beam domain characteristics | |
CN117991245B (en) | Ship azimuth estimation method based on combination of multi-channel secondary DFT and interpolation | |
CN113532617B (en) | Line spectrum detection method for long-term beam phase statistical characteristics |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210420 |