CN108845290A - A kind of method of the anti-phase ambiguity of ultra-short baseline array - Google Patents
A kind of method of the anti-phase ambiguity of ultra-short baseline array Download PDFInfo
- Publication number
- CN108845290A CN108845290A CN201810261916.0A CN201810261916A CN108845290A CN 108845290 A CN108845290 A CN 108845290A CN 201810261916 A CN201810261916 A CN 201810261916A CN 108845290 A CN108845290 A CN 108845290A
- Authority
- CN
- China
- Prior art keywords
- array
- target direction
- direction angle
- array element
- phase
- 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
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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
Abstract
The present invention is to provide a kind of methods of the anti-phase ambiguity of ultra-short baseline array.According to the formation of ultra-short baseline array, the maximum for obtaining phase between any two array element in same direction obscures periodicity, determines the value range of phase ambiguity periodicity between two array elements；The phase difference and maximum fuzzy periodicity measured according to two array elements calculates all possible target direction angle of this group of array element；The covering of the fan spread factor at target direction angle is calculated, the covering of the fan width at target direction angle is set；According to all possible target direction angle and covering of the fan width, the frequency of occurrence of all target direction angles is counted on the histogram, and the most angle of frequency of occurrence is the target direction angle estimated；The estimated value of the ambiguity period number of each group array element is calculated according to the target direction angle of estimation.The present invention can effectively realize anti-phase ambiguity, breach limitation of the array sizes for basic matrix positioning performance.
Description
Technical field
The present invention relates to a kind of Underwater acoustic signal processing methods.
Background technique
Ultra-short baseline array is because its own size is small, the features such as being easily installed, is widely used in target positioning under water.It passes
Ultra-short baseline array of uniting is less than the triangle battle array of transmitting signal half-wavelength, since array aperture is small, traditional ultrashort base frequently with size
Linear array long distance positioning precision is not often high, can solve the problems, such as this by improving formation or increasing array aperture, but this
The problem of a little means can bring phase ambiguity.So taking into account anti-phase mode while improving the positioning accuracy of ultra-short baseline array
The method of paste is particularly significant for ultra short baseline locating system.
Domestic and foreign scholars used a variety of methods to solve the problems, such as this, and representative research mainly has：Explain quick ([1]
Explain quick, Hui Junying, measurement [J] applied acoustics of Sun Dajun ultra-short baseline array primitive phase shift difference, 2006 (04):229-
233.) method using 8 yuan of battle arrays is proposed, the positioning accuracy of raising system is combined by the maximum array element of spacing in same direction,
Phase compensation is carried out by being less than the array element combination of half-wavelength again, realizes anti-phase ambiguity；Zheng Cuie ([2] Zheng Cuie, Lee
Improved method [J] the Chinese Marine University journal of fine jade, Sun Dajun, a kind of ultra short baseline locating system formation of Zhang Dianlun is (natural
Scientific version), 2009,39 (03):The form for 505-508.) changing transmitting signal, is utilized the transmitting signal of dipulse, uses
Orthogonal 4 yuan of battle array positioning modes solve phase fuzzy problem, improve positioning accuracy while eliminating redundancy array element；Zheng Enming
([3] Zheng Enming, CHEN XINHUA, Sun Changyu, a kind of quaternary ultra-short baseline array of remaining China soldier realize high accuracy positioning [J] applied acoustics,
2013,32(01):The quaternary battle array that spacing does not wait 15-22.) is formed by optimization formation, while reducing element number of array, is simplified
The form of transmitting signal, and ensure that and 8 yuan of comparable positioning accuracies of formation in document [1]；But it is limited to physics
Technique, signal higher for frequency, array element spacing are generally difficult to meet the condition for being less than half-wavelength, and the above method is no longer at this time
It is applicable in.
Summary of the invention
It is able to solve the purpose of the present invention is to provide one kind and is less than half when signal frequency is more high-leveled and difficult to meet array aperture
The method of the anti-phase ambiguity of the ultra-short baseline array of generated phase fuzzy problem when the condition of wavelength.
The object of the present invention is achieved like this：
(1) according to the formation of ultra-short baseline array, the maximum fuzzy week of phase between any two array element in same direction is obtained
Issue determines the value range of phase ambiguity periodicity between two array elements；
(2) phase difference measured according to two array elements and maximum fuzzy periodicity calculate all possible mesh of this group of array element
Mark deflection；
(3) the covering of the fan spread factor for calculating target direction angle, is arranged the covering of the fan width at target direction angle；
(4) according to all possible target direction angle and covering of the fan width, all target direction angles are counted on the histogram
Frequency of occurrence, the most angle of frequency of occurrence be estimate target direction angle；
(5) estimated value of the ambiguity period number of each group array element is calculated according to the target direction angle of estimation.
The invention proposes a kind of methods of the anti-phase ambiguity of ultra-short baseline array, to solve when signal frequency is more high-leveled and difficult with full
Sufficient array aperture is less than generated phase fuzzy problem when the condition of half-wavelength.
The beneficial effects of the invention are as follows：(1) signal frequency received in basic matrix is higher, it is difficult to meet array aperture and be less than
Under conditions of half-wavelength, anti-phase ambiguity can be effectively realized, breach limitation of the array sizes for basic matrix positioning performance；
(2) positioning accuracy can be improved compared with traditional triangle battle array, and takes full advantage of the array element information of redundancy, at the angle of concern
When being θ ∈ [- 60 °, 60 °] in degree range, the correctness approximation 100% of anti-phase ambiguity.
Detailed description of the invention
Fig. 1 is the corresponding basic matrix formation schematic diagram of this method；
Fig. 2 is the statistic histogram of each angle frequency of occurrence when target direction angle is spent for -50；
Fig. 3 is the anti-phase ambiguity accuracy of different target orientation angle.
Specific embodiment
It illustrates below and the present invention is described in more detail.
First in this method using to basic parameter be explained as follows：It illustrates in this method the ultra-short baseline array of use
Formation be 8 yuan of cross battle arrays, formation schematic diagram is as shown in Figure 1.The specific value of array element spacing is d14=0.21m, d12=
0.019m, d23=0.026m.Echo signal form is Single Frequency C W signal, and signal frequency f=75kHz, sound is propagated fast in ocean
Degree is c=1500m/s, and phase difference measurement error is up to 6 °.It is carried out in accordance with the following steps after determining basic parameter：
(1) according to the formation of ultra-short baseline array, the maximum for calculating phase between any two array element in same direction is fuzzy
Periodicity determines the value range of phase ambiguity periodicity between two array elements.
If the spacing of array element i and array element j is dij, then array element i and the maximum ambiguity period N of array element j phase are：
The fuzzy periodicity of maximum is the maximum value between two array elements in all possible ambiguity period number, therefore is calculating direction
All ambiguity period numbers are unable to be more than maximum fuzzy periodicity when angle.
If the spacing of array element i and array element j are dij, all possible ambiguity period number is n between two array elementsij, then nijTake
Value range is：
nij∈[-N,N] (2)
According to known formation and parameter, there are 4 array elements in the same direction so there are 6 kinds of array element combinations, corresponding there are 6
A maximum fuzzy periodicity.Spacing dijIt is as shown in the table with the relationship of maximum cycle number N and value range：
(2) phase difference measured according to two array elements and maximum fuzzy periodicity calculate all possible mesh of this group of array element
Mark deflection.
The calculation formula at target direction angle is：
Wherein,Indicate the phase difference measured between array element i and array element j, nijIndicate own between array element i and array element j
Possible ambiguity period number, dijFor the spacing of array element i and array element j, λ is wavelength, θijThe institute solved between array element i and array element j
Possible target direction angle.
According to phase difference and known array element spacing is measured, all possible values of ambiguity period number are successively substituted into formula
(3) all possible target direction angle can be obtained by being calculated.Such as the phase difference measured between array element 12 isAnd because
It is 1 that maximum between array element, which obscures periodicity, so nijIt can be taken as -1,0,1, it willAnd nijIt brings formula (3) into, can acquire
All possible deflection between array element 1 and array element 2.This calculating process is repeated to all array element combinations, institute can be obtained
Possible target direction angle.
(3) the covering of the fan spread factor for calculating target direction angle, is arranged the covering of the fan width at target direction angle.
The size of target direction angle covering of the fan width is codetermined by covering of the fan spread factor and phase difference measurement tolerance.Covering of the fan is wide
Spending coefficient is target direction angle calculation formula to 2 times of the local derviation numerical value of phase difference, and phase difference measurement tolerance is mainly missed with measurement
The maximum value of difference is related.Target direction angle covering of the fan width is equal to the product of covering of the fan spread factor and phase difference measurement tolerance.Such as
For array element 1 and array element 2, phase difference measurement tolerance is 6 °.If measuring phase difference between array element is -50 °, then calculated local derviation
Number is 0.1712, and covering of the fan spread factor is just 0.3424, and it is 2.0544 ° that target direction angle covering of the fan width, which is calculated,.
(4) according to all possible target direction angle and covering of the fan width, all target direction angles are counted on the histogram
Frequency of occurrence, the most angle of frequency of occurrence be estimate target direction angle.
According to calculated all possible target direction angle and covering of the fan width, unite to the frequency of occurrence of different angle
Meter.Statistical method just adds 11 time for the every appearance of any angle in the frequency of occurrence of record, and record institute is angled on the histogram
There is the statistical result of total degree, the most angle of frequency of occurrence is exactly the estimated value at target direction angle.
It chooses -50 ° to be proved as known target deflection, the estimated result at target direction angle is as shown in Figure 2.Fig. 2
Indicating the number occurred using statistics with histogram all angles, it can be seen that the most angle of frequency of occurrence is -50 °, therefore -50 °
It is the estimated value at target direction angle.
(5) estimated value of the ambiguity period number of each group array element is calculated according to the target direction angle of estimation.
If the target direction angle of estimation isBy the target direction angle of estimationPeriod fuzzy number is sought as known quantity substitution
Formula (4), acquire the ambiguity period number of estimation.Modulus paste periodicity formula be：
Wherein,Indicate the phase difference measured between array element i and array element j,Indicate own between array element i and array element j
The estimated value of possible ambiguity period number, dijFor the spacing of array element i and array element j, λ is wavelength,Between array element i and array element j
The estimated value of the azimuth of target solved.It is found out by formula (4)As required ambiguity period number estimated value.
For the validity for verifying the above method, analysis can be carried out to the anti-phase ambiguity of different target deflection below,
As a result as shown in Figure 3.As it can be seen that the entire minimum accuracy of angle domain is higher than 98.8%, and when deflection θ ∈ [- 60 °, 60 °], inspection
It surveys probability and reaches 100%,.Therefore, this method can accurately complete the task of anti-phase ambiguity.
Finally it should be noted that above embodiments are only to describe technical solution of the present invention rather than to this technology method
It is limited, the present invention can above extend to other modifications, variation, application and embodiment, and therefore, it is considered that institute in application
There are such modification, variation, application, embodiment all within the scope of spirit or teaching of the invention.
Claims (1)
1. a kind of method of the anti-phase ambiguity of ultra-short baseline array, it is characterized in that：
(1) according to the formation of ultra-short baseline array, the maximum for obtaining phase between any two array element in same direction obscures periodicity,
Determine the value range of phase ambiguity periodicity between two array elements；
(2) phase difference measured according to two array elements and maximum fuzzy periodicity calculate all possible target side of this group of array element
To angle；
(3) the covering of the fan spread factor for calculating target direction angle, is arranged the covering of the fan width at target direction angle；
(4) according to all possible target direction angle and covering of the fan width, going out for all target direction angles is counted on the histogram
Occurrence number, the most angle of frequency of occurrence are the target direction angle estimated；
(5) estimated value of the ambiguity period number of each group array element is calculated according to the target direction angle of estimation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810261916.0A CN108845290B (en) | 2018-03-28 | 2018-03-28 | Method for resisting phase ambiguity of ultra-short baseline array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810261916.0A CN108845290B (en) | 2018-03-28 | 2018-03-28 | Method for resisting phase ambiguity of ultra-short baseline array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108845290A true CN108845290A (en) | 2018-11-20 |
| CN108845290B CN108845290B (en) | 2022-08-02 |
Family
ID=64211868
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810261916.0A Active CN108845290B (en) | 2018-03-28 | 2018-03-28 | Method for resisting phase ambiguity of ultra-short baseline array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108845290B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703199A (en) * | 2019-10-22 | 2020-01-17 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN112327280A (en) * | 2020-10-27 | 2021-02-05 | 哈尔滨工程大学 | Method for correcting array element number of underwater sound uniform linear array |
| CN113341390A (en) * | 2021-05-25 | 2021-09-03 | 电子科技大学 | Wide-range linear frequency modulation continuous wave radar angle measurement method |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3935574A (en) * | 1974-04-15 | 1976-01-27 | The United States Of America As Represented By The Secretary Of The Navy | Signal source position-determining process |
| US4639733A (en) * | 1983-05-11 | 1987-01-27 | Racal Communications Equipment Limited | Direction finding |
| CN2854623Y (en) * | 2005-11-17 | 2007-01-03 | 哈尔滨工程大学 | Receiving base distribution structure of ultrashort base-line positioning system |
| JP2011013031A (en) * | 2009-06-30 | 2011-01-20 | Toshiba Corp | Arrival direction estimating apparatus |
| CN102411136A (en) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | Phase interferometer direction finding method for ambiguity resolution by extension baselines |
| CN104330767A (en) * | 2014-11-27 | 2015-02-04 | 中国船舶重工集团公司第七二四研究所 | Interferometer direction-finding method based on phase fuzzy number search and least square fit |
| CN104991573A (en) * | 2015-06-25 | 2015-10-21 | 北京品创汇通科技有限公司 | Locating and tracking method and apparatus based on sound source array |
-
2018
- 2018-03-28 CN CN201810261916.0A patent/CN108845290B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3935574A (en) * | 1974-04-15 | 1976-01-27 | The United States Of America As Represented By The Secretary Of The Navy | Signal source position-determining process |
| US4639733A (en) * | 1983-05-11 | 1987-01-27 | Racal Communications Equipment Limited | Direction finding |
| CN2854623Y (en) * | 2005-11-17 | 2007-01-03 | 哈尔滨工程大学 | Receiving base distribution structure of ultrashort base-line positioning system |
| JP2011013031A (en) * | 2009-06-30 | 2011-01-20 | Toshiba Corp | Arrival direction estimating apparatus |
| CN102411136A (en) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | Phase interferometer direction finding method for ambiguity resolution by extension baselines |
| CN104330767A (en) * | 2014-11-27 | 2015-02-04 | 中国船舶重工集团公司第七二四研究所 | Interferometer direction-finding method based on phase fuzzy number search and least square fit |
| CN104991573A (en) * | 2015-06-25 | 2015-10-21 | 北京品创汇通科技有限公司 | Locating and tracking method and apparatus based on sound source array |
Non-Patent Citations (3)
| Title |
|---|
| ANTHONY GRIFFIN 等: "Real-time multiple speaker DOA estimation in a circular microphone array based on Matching Pursuit", 《2012 PROCEEDINGS OF THE 20TH EUROPEAN SIGNAL PROCESSING CONFERENCE (EUSIPCO)》 * |
| 王燕 等: "超短基线定位系统融合分类解模糊技术研究", 《电子与信息学报》 * |
| 王燕 等: "超短基线水声定位系统相位差解模糊算法研究", 《电子学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703199A (en) * | 2019-10-22 | 2020-01-17 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN110703199B (en) * | 2019-10-22 | 2021-09-24 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN112327280A (en) * | 2020-10-27 | 2021-02-05 | 哈尔滨工程大学 | Method for correcting array element number of underwater sound uniform linear array |
| CN112327280B (en) * | 2020-10-27 | 2022-05-17 | 哈尔滨工程大学 | Method for correcting array element number of underwater sound uniform linear array |
| CN113341390A (en) * | 2021-05-25 | 2021-09-03 | 电子科技大学 | Wide-range linear frequency modulation continuous wave radar angle measurement method |
| CN113341390B (en) * | 2021-05-25 | 2022-06-03 | 电子科技大学 | Wide-range linear frequency modulation continuous wave radar angle measurement method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108845290B (en) | 2022-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108845290A (en) | A kind of method of the anti-phase ambiguity of ultra-short baseline array | |
| CN1996986B (en) | Full phase time shift phase difference spectrum correction method | |
| CN107747931B (en) | For the self-alignment reading head optimization placement method of angular encoder | |
| CN103292701A (en) | Machine-vision-based online dimensional measurement method of precise instrument | |
| CN108802683A (en) | A kind of source localization method based on management loading | |
| CN108548578B (en) | Ultrasonic echo signal characteristic peak identification method based on self-adaptive threshold | |
| CN109782287A (en) | A method of it is surged wave height by high-frequency ground wave radar inverting | |
| CN108152687B (en) | Power transformer partial discharge positioning method utilizing ultrasonic inversion | |
| CN104913830A (en) | Radar level instrument signal processing method based on curve fitting | |
| CN106840247A (en) | A kind of linear sensor bearing calibration and device | |
| CN105891795B (en) | A kind of dual station rcs measurement calibrating method | |
| CN110231620A (en) | A kind of noise correlation system tracking filter method | |
| CN108566677A (en) | A kind of fingerprint positioning method and device | |
| CN103323473A (en) | Method for measuring inverse pole figure | |
| CN107037416A (en) | Time domain door transform method based on frequency sweep rcs measurement | |
| CN104266732B (en) | Dual stage bit stream gauge calibrating installation step-by-step counting precision compensation method | |
| CN103913203A (en) | Ultrasonic water meter coefficient processing method | |
| Rychagov et al. | Multipath flowrate measurements of symmetric and asymmetric flows | |
| CN110057328A (en) | One kind is applied to ward off tin thickness measurement method | |
| CN113343492A (en) | Theoretical spectral data optimization method and system and optical measurement method | |
| CN110360937B (en) | Automatic paper roll width measurement method based on laser range finder | |
| CN108169634B (en) | Method for accurately obtaining partial discharge ultrahigh frequency signal time difference | |
| CN109916361A (en) | A kind of roundness measurement signal processing method without angular position information | |
| CN116359835B (en) | Y-type baseline interferometer direction finding device and direction finding method | |
| Yazdanshenasshad et al. | Reducing the additional error caused by the time‐difference method in transit‐time UFMs |
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 |