CN105954707B - A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid - Google Patents

A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid Download PDF

Info

Publication number
CN105954707B
CN105954707B CN201610248569.9A CN201610248569A CN105954707B CN 105954707 B CN105954707 B CN 105954707B CN 201610248569 A CN201610248569 A CN 201610248569A CN 105954707 B CN105954707 B CN 105954707B
Authority
CN
China
Prior art keywords
isosceles trapezoid
battle array
array
receives
array element
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
Application number
CN201610248569.9A
Other languages
Chinese (zh)
Other versions
CN105954707A (en
Inventor
韩宁
方世良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201610248569.9A priority Critical patent/CN105954707B/en
Publication of CN105954707A publication Critical patent/CN105954707A/en
Application granted granted Critical
Publication of CN105954707B publication Critical patent/CN105954707B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/02—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 radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Preparing Plates And Mask In Photomechanical Process (AREA)

Abstract

The invention discloses a kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid, including step：Step (1)：Array element construction isosceles trapezoid is received with four and receives battle array, it is made to meet L1＞ L2Cos Φ and L1≠2L2cosΦ；Wherein, L1Upper bottom edge length, the L of battle array are received for isosceles trapezoid2The waist length of battle array is received for isosceles trapezoid, Φ is that isosceles trapezoid receives the waist of battle array and the angle of bottom；Step (2)：Measuring signal reaches the phase difference Ω that isosceles trapezoid receives two reception array element above and below waist on the left of battle array respectively1, signal reach isosceles trapezoid receive battle array on the right side of waist up and down two reception array element phase difference Ω2, signal reaches isosceles trapezoid and receives receives array element on the right side of bottom and left side receives the phase difference Ω of array element in battle array3；Step (3)：The combination calculating signal arrival bearing of phase difference between array element is received using isosceles trapezoid.The present invention reduces multivalue problem present in phase direction-finding method using the combination of phase difference between array element.

Description

A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid
Technical field
It is particularly a kind of that subtracting for battle array is received based on isosceles trapezoid the present invention relates to a kind of method for reducing direction finding multivalue problem The method of small direction finding multivalue problem.
Background technology
Direction finding technology has ten in numerous application fields such as radar, sonar, communication, navigation, geo-location, military weapon Divide important role.Traditional direction finding technology includes amplitude of vibration method and phase method, [Reuben Shar, Simulation of Polarizer Impact on Circular Interferometer performance.The 11th IEEE International Conference on Computational Science and Engineering,IEEE, 2008: 334-339], wherein phase method direction finding is to judge radiation source by measuring the phase difference of plane wave signal arrival reception array element Incident direction, the stability of algorithm realization is high, system structure is simple, convenient for Project Realization.In practical applications, phase method by The limitation of radiation wavelength and array element spacing：When array element spacing is more than radiation signal half-wavelength, the multivalue in direction to be measured will be caused Problem [Ernest Jacobs, Elizabeth W.Ralston. Ambiguity resolution in interferometry.IEEE Transactions on Aerospace and Electronic Systems,1981,17 (6):766-779].Therefore the research for carrying out reduction direction finding multivalue problem is of great significance.
The method for solving direction finding multivalue problem at present mainly has：Long-short baselines method, virtual baseline method, is stood remainder theorem method Body base-line method etc..There is one group of baseline for being less than radiation signal half-wavelength in long-short baselines method requirement receiving array, often by physics The limitation of condition can not be realized；Remainder theorem method and virtual baseline method each group baseline length need not be less than signal half-wavelength, but real It is often limited in the application of border by calculation method, as remainder theorem method needs each group baseline all to meet stringent relatively prime relationship；It is vertical The spatial relation of each array element in receiving array is rationally utilized in body baseline direction finding method, and algorithm is to each group array element spacing without spy Different requirement, but required array element number is more.With the development and application of array technique, how rationally to lay in a limited space The problem of receiving array element, and effectively reduce multivalue problem present in phase direction-finding method, being still current constantly research and improve. On the other hand, it on the basis of keeping existing reception formation constant, constructs new formation using array elements several in formation, be aided with newly Direction finding calculation method, can also effectively reduce existing multivalue problem during phase direction finding.
Invention content
The object of the present invention is to provide a kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid, the party Method can effectively reduce multivalue problem present in phase direction-finding method.The technical solution of use is as follows：
A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid, including step：
Step (1)：Array element construction isosceles trapezoid is received with four and receives battle array, it is made to meet L1＞ L2Cos Φ and L1≠ 2L2cosΦ；Wherein, L1Upper bottom edge length, the L of battle array are received for isosceles trapezoid2The waist length of battle array is received for isosceles trapezoid, Φ is isosceles ladder Shape receives the waist of battle array and the angle of bottom；
Step (2)：Measuring signal reaches the phase difference that isosceles trapezoid receives two reception array element above and below waist on the left of battle array respectively Ω1, signal reach isosceles trapezoid receive battle array on the right side of waist up and down two reception array element phase difference Ω2, signal reach isosceles trapezoid connect Receive the phase difference Ω for receiving array element and left side reception array element in battle array on the right side of bottom3；
Step (3)：Calculate signal arrival bearing
Wherein f0For signal frequency, c is the spread speed of signal, and n takes all satisfactions Integer value.
The advantages of the present invention are：The present invention receives the combination progress of phase difference between array element using isosceles trapezoid Orientation measurement can effectively reduce direction finding multivalue problem present in original technology.
Description of the drawings
Fig. 1 is the schematic diagram for the reduction direction finding multivalue methods that battle array is received based on isosceles trapezoid；
Fig. 2 is the specific implementation process of the embodiment of the present invention；
Fig. 3 is orientation measurement result of the inventive embodiments 1 to different directions signal.
Fig. 4 is orientation measurement result of the inventive embodiments 2 to different directions signal.
Specific embodiment
Technical scheme of the present invention is described in detail below in conjunction with the accompanying drawings：
As shown in Figure 1, four receive array element and are located at four vertex that the isosceles trapezoid receives battle array, respectively a, b, c, d, Length between wherein ab is L1, the length of ad and bc are L2, the length between cd is L3=L1The height of+2*l, the then isosceles trapezoid ForAngle, φ=arctan (H/l) that ad lines and cd lines are formed.If signal arrival bearing and level side To angle be θ, then
D1=L2* cos (Φ-θ)=τ1C, D2=L2* cos (π-Φ-θ)=τ2C, D3=L1Cos θ=τ3c (1)
In above formula, spread speeds of the c for signal, frequency f0Signal reach 2 points of a, d delay inequality be τ1, phase difference For Ω1, then Ω1=2 π f0τ1+2πn1,Wherein n1It is an integer, makes-π ＜ Ω1≤π；Signal reach 2 points of b, c delay inequality be τ2, phase difference Ω2, then Ω2=2 π f0τ2+2πn2, wherein n2It is an integer, makes-π ＜ Ω2≤π；Signal reaches 2 points of a, b's Delay inequality is τ3, phase difference Ω3, then Ω3=2 π f0τ3+2πn3,Wherein n3It is an integer, makes-π ＜ Ω3≤π。
If carrying out direction finding using ab array elements, direction to be measured is calculated
During practical direction finding, n3It is unknown, so the n in formula (2)3All satisfactions need to be chosen's Integer value carries out againSolution, therefore multiple directions to be measured can be calculated, that is, the multivalue problem in phase direction-finding method occur.
If receiving battle array using the isosceles trapezoid carries out direction finding, because
So
D3-(D1-D2)=L1cosθ-2L2* cos θ * cos Φ=τ3c-(τ1-τ2)c (4)
Thus direction to be measured can be calculated
Wherein n=n3-(n1-n2) it is integer.During practical direction finding, n is unknown, so the n in formula (5) need to choose institute There is satisfactionInteger value carry out againSolution, can also calculate multiple directions.It is if such The trapezoidal battle array of waist meets L1＞ L2Cos Φ and L1≠2L2Cos Φ, contrast equation (2) and formula (5), then | 2 π f0(L1-2L2*cos Φ)/c | ＜ | 2 π f0L1/ c |, then Ω3-(Ω1-Ω2) value range of -2 π n compares Ω3-2πn3It is small, so n in formula (5) Value number can be less than formula (2) in n3Value number, so as to reduce direction finding multivalue problem present in original technology.It is special It is other, if the ＜ L of-λ/21-2L2* cos Φ ＜ λ/2, at this time n take unique value, unique signal arrival bearing can be calculated.
Fig. 2 is the specific implementation process for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid.
As shown in Fig. 2, the reduction direction finding multivalue problem that battle array is received the present invention is based on isosceles trapezoid is as follows：
(1) it receives array element construction isosceles trapezoid with four and receives battle array, it is made to meet L1＞ L2Cos Φ and L1≠2L2cos Φ, the wherein isosceles trapezoid receive a length of L in upper bottom edge of battle array1, a length of L of waist2, a length of L in bottom3=L1+ 2l, then the isosceles are terraced The height of shapeWaist and the angle Φ=arctan (H/l) on base.Such isosceles trapezoid reception battle array, which has, to be subtracted The advantage of small direction finding multivalue.
(2) measuring signal reaches the phase difference Ω of two array elements above and below left waist1, reach right waist up and down two array elements phase difference Ω2, the phase difference Ω of right left two array element in bottom in arrival3。
(3) signal arrival bearing is calculated
Wherein λ be signal wavelength, f0For signal frequency, c is the spread speed of signal, and n takes all satisfactionsInteger value.
Fig. 3 is a specific embodiment 1 according to Fig. 2.In the isosceles trapezoid receives battle array, the length between ab is L1= The length of 2m, ad and bc are L2Length between=2.06m, cd is L3=3m, trapezoidal height are H=2m, ad lines and cd lines Angle, φ=76 ° of composition.The isosceles trapezoid receives battle array and meets L1＞ L2Cos Φ and L1≠2L2cosΦ.The radiation of signal source Frequency f0=500Hz, signal velocity c=1500m/s, then signal wavelength lambda=3m.The incident direction in setting signal source is surveyed Measure the phase difference Ω that signal reaches two array element of a, d1, the phase difference Ω of arrival two array element of b, c2, the phase difference of arrival two array element of b, a Ω3.Because of-λ/2 ＜ L1-2L2* cos Φ ＜ λ/2, then unique signal incoming wave side can be calculated by receiving battle array using the isosceles trapezoid To.
The actual direction of signal source is to be spaced apart with 10 degree from 0 degree to 180 degree in Fig. 3；The result measured with array element ab There is multivalue in 0~40 degree and 140~180 degree, angle and practical incident direction are basically identical wherein shown in black matrix, are considered as just True measurement result, angle shown in italic are erroneous measurements caused by multivalue phenomenon；The direction measured with array element cd 0~ There is multivalue phenomenon in 180 degree, and three measurement results occurs at 0,90,180 degree, angle and reality wherein shown in black matrix Border incident direction is basically identical, is considered as correct measurement result, and mistake caused by angle shown in italic is multivalue phenomenon measures knot Fruit；When carrying out orientation measurement using the isosceles trapezoid battle array, do not occur multivalue phenomenon, and measurement result and practical incident direction are basic Unanimously.Therefore the invention can effectively reduce the multivalue problem occurred during direction finding, and then obtain correct azimuthal measurements. Fig. 4 is a specific embodiment 2 according to Fig. 2.Except the radiation frequency f of signal source0Outside=800Hz, remaining parameter and Fig. 3 phases Together.The actual direction of signal source is to be spaced apart with 10 degree from 0 degree to 180 degree in Fig. 4；It is surveyed with array element ab, cd, isosceles trapezoid battle array There is phase multivalue phenomenon in 0~180 degree in the result of amount, angle and practical incident direction basic one wherein shown in black matrix It causes, is considered as correct measurement result, angle shown in italic is erroneous measurements caused by multivalue phenomenon.As seen from the figure, although Direction finding, which is carried out, using the isosceles trapezoid battle array still has phase multivalue problem, but its obtained number of measurements is minimum.Therefore should Invention can effectively reduce the multivalue problem occurred during phase direction finding.

Claims (1)

  1. A kind of 1. method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid, which is characterized in that including step：
    Step (1)：Array element construction isosceles trapezoid is received with four and receives battle array, it is made to meet L1＞ L2Cos Φ and L1≠2L2cosΦ； Wherein, L1Upper bottom edge length, the L of battle array are received for isosceles trapezoid2The waist length of battle array is received for isosceles trapezoid, Φ is that isosceles trapezoid receives battle array Waist and bottom angle；
    Step (2)：Measuring signal reaches the phase difference Ω that isosceles trapezoid receives two reception array element above and below waist on the left of battle array respectively1, letter Number reach isosceles trapezoid receive battle array on the right side of waist up and down two reception array element phase difference Ω2, signal arrival isosceles trapezoid reception battle array Array element is received on the right side of upper bottom and left side receives the phase difference Ω of array element3；
    Step (3)：Calculate signal arrival bearing
    Wherein f0For signal frequency, c is the spread speed of signal, and n takes all satisfactions
    Integer value.
CN201610248569.9A 2016-04-20 2016-04-20 A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid Active CN105954707B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610248569.9A CN105954707B (en) 2016-04-20 2016-04-20 A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610248569.9A CN105954707B (en) 2016-04-20 2016-04-20 A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid

Publications (2)

Publication Number Publication Date
CN105954707A CN105954707A (en) 2016-09-21
CN105954707B true CN105954707B (en) 2018-06-12

Family

ID=56917823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610248569.9A Active CN105954707B (en) 2016-04-20 2016-04-20 A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid

Country Status (1)

Country Link
CN (1) CN105954707B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108480921A (en) * 2018-04-11 2018-09-04 西安电子工程研究所 A kind of method of rapid processing corner cut waveguide

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6839025B1 (en) * 2002-06-03 2005-01-04 Ken Reigle Precision direction finding sensing systems and methods
CN101105525A (en) * 2007-07-06 2008-01-16 哈尔滨工程大学 Pure phase type broad frequency band microwave radiation source direction finding system and method
CN101109799A (en) * 2007-07-06 2008-01-23 哈尔滨工程大学 Sense finding device for solid base line broad band microwave radiation source and sense finding method thereof
CN101403791A (en) * 2008-11-06 2009-04-08 哈尔滨工程大学 Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof
RU2419805C1 (en) * 2009-09-23 2011-05-27 Государственное образовательное учреждение высшего профессионального образования "Военная академия связи имени С.М. Буденного" Министерства обороны Российской Федерации Method of radio signal direction finding and direction finder to this end
CN102411136A (en) * 2011-08-09 2012-04-11 电子科技大学 Phase interferometer direction finding method for ambiguity resolution by extension baselines

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6839025B1 (en) * 2002-06-03 2005-01-04 Ken Reigle Precision direction finding sensing systems and methods
CN101105525A (en) * 2007-07-06 2008-01-16 哈尔滨工程大学 Pure phase type broad frequency band microwave radiation source direction finding system and method
CN101109799A (en) * 2007-07-06 2008-01-23 哈尔滨工程大学 Sense finding device for solid base line broad band microwave radiation source and sense finding method thereof
CN101403791A (en) * 2008-11-06 2009-04-08 哈尔滨工程大学 Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof
RU2419805C1 (en) * 2009-09-23 2011-05-27 Государственное образовательное учреждение высшего профессионального образования "Военная академия связи имени С.М. Буденного" Министерства обороны Российской Федерации Method of radio signal direction finding and direction finder to this end
CN102411136A (en) * 2011-08-09 2012-04-11 电子科技大学 Phase interferometer direction finding method for ambiguity resolution by extension baselines

Also Published As

Publication number Publication date
CN105954707A (en) 2016-09-21

Similar Documents

Publication Publication Date Title
CN104122527B (en) A kind of round battle array phase-interferometer broadband based on look-up table instantaneous direction finding method
CN108875099B (en) Baseline selection method based on long and short baseline interferometer direction-finding system
CN102175988B (en) Direction finding method of correlative interferometer based on dimension division
CN105403904B (en) A kind of satellite navigation single-frequency based on aerial array surveys attitude positioning method
Van Doan et al. Optimized algorithm for solving phase interferometer ambiguity
CN103091660A (en) Virtual base line direction-finding method based on phase jump correction
CN106249057B (en) It is a kind of large-scale around focus rotation beam scanning antenna radiation characteristics equivalent detecting method
CN113504503B (en) Ambiguity resolution method for rotary multi-baseline phase interferometer
CN106199503A (en) Correlation interferometer direction-finding method and device
CN104699965A (en) Estimation method of parameters of near-field source based on rotating interferometer
KR100902559B1 (en) Radar interferometer and method of target position estimation using the same
CN104614736A (en) Calibration method of GPS receiver
CN105954707B (en) A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid
CN109215072A (en) Chaff cloud RCS acquisition methods based on tanh-sinh integration method
CN215579031U (en) Rotary multi-baseline phase interferometer direction-finding antenna array
CN109029744B (en) A kind of detection method of vortex beams topological charge number
Kang et al. An Efficient Hybrid Multi-Station TDOA and Single-Station AOA Localization Method
CN102735207B (en) Two-dimension angle calculation method based on distance and carrier phase difference
CN105334490B (en) A kind of method of broadband signal interference direction-finding
CN109407089B (en) A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view
CN108205128A (en) A kind of passive measurement distance method based on long baseline interferometer
KR102022589B1 (en) The method for calculating ambiguity probability of the direction of arrival
Bondesson et al. Standard errors of area estimates obtained by traversing and GPS
CN102736057B (en) Two-dimensional angle calculating method based on distance and tracking regeneration carrier phase
KR101241926B1 (en) Conformity evaluation method of radome for 3 dimensional interferometric array by phase difference error estimation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant