CN112162233A - Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment - Google Patents
Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment Download PDFInfo
- Publication number
- CN112162233A CN112162233A CN202011040289.1A CN202011040289A CN112162233A CN 112162233 A CN112162233 A CN 112162233A CN 202011040289 A CN202011040289 A CN 202011040289A CN 112162233 A CN112162233 A CN 112162233A
- Authority
- CN
- China
- Prior art keywords
- angle
- radio frequency
- radiation source
- port
- baseline
- 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
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/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
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/48—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention discloses a two-dimensional wide-angle high-precision angle measuring method based on eight-port four-baseline radio frequency equipment, which is characterized by comprising the following steps of: step 1, building eight-port four-baseline radio frequency equipment; step 2, determining a baseline interval and a frequency factor according to the angle range and the wavelength of the radiation source; step 3, determining the number of steps where the pitch angle of the radiation source is located; step 4, restoring the phase value of the actual pitch angle; step 5, calculating a radiation source pitch angle; step 6, determining the number of steps where the azimuth angle of the radiation source is located; step 7, restoring the actual azimuth angle phase value; and 8, calculating the azimuth angle of the radiation source. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency equipment is simple and convenient, has a wider angle measurement range and higher precision under the condition of the same cost or volume or number of antenna units, and can be used for engineering practice.
Description
Technical Field
The invention belongs to the technical field of radar radiation source direction finding methods, and particularly relates to a two-dimensional wide-angle high-precision angle finding method based on eight-port four-baseline radio frequency equipment.
Background
The radar radiation source direction finding technology can demodulate the direction of electromagnetic waves by utilizing the principle that the electromagnetic waves in different directions reach different amplitudes or phase responses generated by a direction finding antenna system. Depending on different factors, it can be classified into an amplitude method and a phase method. The phase method angle measurement utilizes the wave path difference between the electromagnetic wave signals radiated by a target and the antenna base line to measure the angle, namely, the direction of the signals is determined according to the relative phase difference of the same signals detected and received by a direction-finding antenna system, and then angle error signals are demodulated through the phase difference, and the antenna is driven to perform passive tracking on a radiation source. The relative phase difference is derived from the ratio of the relative wave path difference to the wavelength, and the principle is simpler. However, the phase method has a great disadvantage that measurement is blurred when the baseline width between the two antennas is too large, and a large measurement error is caused when the baseline width is too small. The novel wide-angle and high-precision angle measurement method researched by the technical scheme is based on phase method direction measurement, a set of simple and feasible novel angle measurement method is researched, and the problem of measurement ambiguity existing in the traditional phase method angle measurement is solved.
One, two base line angle measuring principle
At present, the angle measuring device based on the phase method direction measurement is mainly an interferometer, which uses the angle measuring principle of the double-baseline phase method, that is, uses the phase difference between echo signals received by a plurality of antennas to measure the angle. As shown in fig. 1, when the target radiates an electromagnetic wave signal in the θ direction, the electric wave reflected by the target reaching the receiving point is approximated to a plane wave. Because the base line interval between the two antennas is d, the signals received by the two antennas reach the difference of the wave path lengths of the two base lines, namely Delta R, so as to generate phase differenceThe relationship between the phase difference and the baseline interval is:
where λ is the wavelength of the electromagnetic wave signal radiated by the target. The phase difference resulting from the wave path difference can be measured by a phase meter. Therefore, the angle θ of the electromagnetic wave signal radiated from the target can be derived from equation (1)
Knowing the wavelength of the electromagnetic wave signal radiated from the target, the azimuth of the target signal can be calculated from the equation (2) by using the phase difference measured by the phase meter.
Second, the problems of angle measurement blur and precision
The simple two-baseline phase method direction finding actually has a great problem, namely the measurement ambiguity problem.
From the equation (2), if the phase difference is smallInaccurate value measurement will result in angle measurement errors. In order to research the relevant factors influencing the angle measurement precision, the two sides of the formula (1) are differentiated, namely
As can be seen from the formula (3), the reading accuracy is highThe angle measurement precision can be improved by the phase meter or by reducing the lambda/d value. In addition, when θ is 0, that is, when the target is in the antenna normal direction, the angle measurement error d θ is minimum, and when θ increases, d θ also increases, so that the range of θ is also limited to a certain extent to ensure a certain angle measurement accuracy. Although the angle measurement accuracy can be improved by reducing the value of λ/d, when the value of λ/d is reduced to a certain extent in a certain angle measurement range θ,the value may exceed 2 pi, in which caseWhere N is an integer, psi < 2 pi, and the actual reading of the phase meter is psi. Since the value of N is unknown, it is trueThe value cannot be determined and a blurring problem (multivalue) occurs.
Third, the current method for solving the fuzzy problem
The measurement range is reduced when the measurement range is reduced, which means that the measurement accuracy and the measurement range are contradictory. The key to solve the contradiction is to solve the ambiguity problem, so how to deblur becomes a hot spot to be considered by applying the phase method to direction finding, and a plurality of methods for solving the ambiguity problem are researched and developed.
1. Stagger baseline deblurring
In order to solve the problem of direction finding ambiguity in the phase method direction finding, a direction finding method of a staggered baseline interferometer for solving the ambiguity problem by using the Chinese remainder theorem is provided by imitating the multi-frequency continuous wave distance measuring technology, the remainder theorem is applied to the direction finding of the interferometer, and a basic angle measuring schematic diagram of the ambiguity solving method is shown in FIG. 2.
An M-dimensional baseline interferometer with length of li(i-1, 2.., M-1), taking a base baseline l0≤λ min2, all base lines are l0Integer multiple of (a) of
l1:l2:…:lM-1=m1:m2:…:mM-1 (4)
Wherein m isi(i ═ 1, 3.., M-1) is an integer
The interferometer measures the spacing between the base lines as liWhen corresponding to a phase difference ofAnd the actual phase difference is 2 pi lisin theta/lambda in a relationship of
Wherein k isiIndicates that the spacing between the base lines is liThe number of direction finding ambiguities in time.
Equation (5) is a same remainder with the same remainder in the real number domain with the divisor being an integerThe equation set, if two are selected to be co-prime, can be known according to the Chinese remainder theoremHaving a unique set of solutions k within the determined maximum unambiguous direction finding rangei. However, in this method, phase errors caused by antenna elements, microwave channels, receivers, and the like easily cause ambiguity resolution failure, and the calculation amount is large.
2. Virtual baseline disambiguation
The term "virtual baseline" refers to the difference in length between two different baselines. When the length difference is smaller than the half wavelength of the highest frequency of the broadband signal, the phase difference of the virtual baselines is the unambiguous phase. The schematic diagram is shown in FIG. 3, where the base line intervals of base lines 1 and 2 and base lines 2 and 3 are respectively l1,l2(l2>l1) Subtracting two different baseline intervals can yield an interval of l2-l1The baseline interval and the corresponding phase difference of the virtual short baselineIn a relationship of
However, in the wide-angle direction finding, the virtual baseline method causes ambiguity resolution errors due to the influence of system errors and random errors, and even cannot resolve ambiguity.
3. Long and short baseline ambiguity resolution
The long and short baseline method is also called a three-baseline angle measurement method, and is implemented by using three baselines with proper two different baseline intervals, wherein one baseline is long and the other baseline is short. The schematic diagram is shown in fig. 4, 1 and 3 antennas with large intervals are used for obtaining high-precision measurement, and 1 and 2 antennas with small intervals are used for solving the measurement multivalue. Assuming that the electromagnetic wave signal with the target radiating direction theta is outwards, the distance between the antennas 1 and 2 is d12The distance between the antennas 1,3 is d13. By selecting d appropriately12So that the phase difference between the signals received by the antennas 1 and 2 can satisfy the requirement in the angle measurement range
According to requirements, larger d is selected13The phase difference of the signals received by the antennas 1,3 is
Where the phase meter 2 reads psi less than 2 pi, the following relationship may be used to determine the value of N
When the error of the phase meter 1 is within the acceptable range, the reading of the phase meter 1 is takenAnd equation (10) can calculateThen, according to the formula (9), the value of N and the value of theta can be determined. d13The/lambda value is larger, and the required precision is ensured.
Although the method for measuring the angle of the long and short baselines can solve the problem of ambiguity of direction measurement by using the short baselines and the problem of direction measurement range by using the long baselines, in the broadband direction measurement, when a target signal is a high-frequency signal, high requirements are required on the short baselines, and the method cannot be widely applied due to engineering limitation of physical realization of the short baselines.
Disclosure of Invention
The invention aims to provide a two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment, and solves the problem that the existing angle measurement method cannot realize high measurement precision and wide measurement range at the same time.
The technical scheme adopted by the invention is as follows: the two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency equipment comprises the following steps of:
and 8, calculating the azimuth angle of the radiation source.
The invention has the beneficial effects that: the two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency equipment is simple and convenient, has a wider angle measurement range (can measure all angles) and higher precision under the condition of the same cost or volume or the number of antenna units, and can be used for engineering practice.
Drawings
FIG. 1 is a schematic diagram of a dual baseline goniometry;
FIG. 2 is a ragged baseline deblurring diagram;
FIG. 3 is a diagram of a virtual baseline deblurring principle;
FIG. 4 is a long and short baseline deblurring schematic;
FIG. 5 is a schematic diagram of the principle of eight-port angle measurement;
FIG. 6 is a test chart of two baselines of an antenna in an eight-port goniometric system;
FIG. 7 is a diagram showing the positional relationship between the electromagnetic wave to be measured and the antenna in the eight-port angle measurement system;
FIG. 8 is a schematic diagram of an eight port junction;
FIG. 9 shows the fuzzy phase value θ measured by an eight-port device when q is 0.8EAnd thetaH;
FIG. 10 shows an angle value θ obtained by resolving the ambiguity phase values measured by an eight-port device at two frequencies when q is 0.8E;
Fig. 11 shows angle values θ calculated at two frequencies when q is 0.8ESubtracting to obtain a step map;
fig. 12 shows the difference θ between two frequencies when q is 0.8EThe value is calculated by using the fuzzy phase value measured by the eight-port device to obtain an angle value thetaH;
Fig. 13 shows the difference θ between two frequencies when q is 0.8EValue of fuzzy phase value theta measured by using eight-port deviceHSubtracting the calculated angle value to obtain a step map;
fig. 14 shows a blur-free angle value θ obtained by reducing the measured angle value by the step method when q is 0.8E;
In fig. 15, q is 0.8 and θ is differentEReducing the measured angle value by a step method under the value to obtain a non-fuzzy angle value thetaH;
Fig. 16 shows an unambiguous angle value θ obtained by reducing the measured angle value by the step method in the range of plus and minus 90 degrees when q is 0.8E;
In fig. 17, q is 0.8 and θ is differentEWithin the range of plus or minus 90 degrees, the step method is utilized to reduce the measured angle value to obtain the non-fuzzy angle value thetaH;
Fig. 18 shows angle values θ obtained by solving for two frequencies at plus and minus 90 degrees when q is 0.6E;
Fig. 19 shows angle values θ obtained by solving for two frequencies at plus and minus 90 degrees when q is 0.6ESubtracting to obtain a step map;
fig. 20 shows an unambiguous angle value θ obtained by restoring the measured angle value by a step method at plus or minus 90 degrees when q is 0.6E;
Fig. 21 shows θ at two frequencies of plus and minus 90 degrees when q is 0.6EThe value is calculated by using the fuzzy phase value measured by the eight-port device to obtain an angle value thetaH;
Fig. 22 shows θ at two frequencies of plus or minus 90 degrees when q is 0.6EThe value is calculated by using the fuzzy phase value measured by the eight-port device to obtain an angle value thetaHSubtracting to obtain a step map;
fig. 23 shows θ at two frequencies of plus and minus 90 degrees when q is 0.6EThe value utilizes a step method to reduce the measured angle value to obtain a non-fuzzy angle value thetaH;
Fig. 24 shows an angle value θ obtained by solving for two frequencies at plus and minus 90 degrees when q is 0.3ESubtracting to obtain a step map;
fig. 25 shows an unambiguous angle value θ obtained by reducing the measured angle value by a step method at plus or minus 90 degrees when q is 0.3E;
Fig. 26 shows θ at two frequencies of plus and minus 90 degrees when q is 0.3EThe value is calculated by using the fuzzy phase value measured by the eight-port device to obtain an angle value thetaHSubtracting to obtain a step map;
fig. 27 shows θ at two frequencies of plus and minus 90 degrees when q is 0.3EThe value utilizes a step method to reduce the measured angle value to obtain a non-fuzzy angle value thetaH;
Fig. 28 shows angle values θ obtained by solving for two frequencies at plus and minus 90 degrees when q is 0.1ESubtracting to obtain a step map;
fig. 29 is a non-blur angle value θ obtained by restoring the measured angle value by the step method at plus or minus 90 degrees when q is 0.1E;
Fig. 30 shows θ at two frequencies of plus and minus 90 degrees when q is 0.1EThe value is calculated by using the fuzzy phase value measured by the eight-port device to obtain an angle value thetaHSubtracting to obtain a step map;
fig. 31 shows θ at two frequencies of plus and minus 90 degrees when q is 0.1EThe value utilizes a step method to reduce the measured angle value to obtain a non-fuzzy angle value thetaH。
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The invention provides a two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment, which comprises the following steps of firstly building the eight-port four-baseline radio frequency equipment, namely providing the radio frequency equipment, wherein two ports are arranged in each direction around the radio frequency equipment, baselines are arranged on four ports of the radio frequency equipment opposite to the two directions, and the other four ports of the radio frequency equipment are respectively connected with a Schottky diode detector:
firstly, determining the length and frequency factor of the base line according to the angle range and wavelength of the radiation source
Firstly, a schematic diagram of the method is given, as shown in fig. 5, four large dots are a top view of the base line position of the antenna, the interval between the base lines 5 and 7 and the interval between the base lines 6 and 8 are d, the same four schottky diode detectors are respectively connected to the ports 1,2, 3 and 4, and the radiation source angle can be calculated through the reading of the detectors. Fig. 6 is a side view of two base wires of the antenna, and the other two base wires are the same as the side view of the antenna in the figure, but are perpendicular to the base wires.
Assuming an angle theta of an electromagnetic wave signal radiated from a targetEAnd thetaHWithin (-80 degrees and 80 degrees), d is 18mm for convenient engineering realization. Unlike the existing methods, the method of the invention only requires the ratio q to f of the two frequencies at which the system operates1/f2The frequency factor is fixed, so on the other hand, except some special values, the frequency factor can be arbitrarily taken, where q is 0.8, 24GHz is taken as the center frequency, two working frequencies are 21.3GHz and 26.7GHz, and the corresponding wavelength is λ114.1mm and λ2=11.2mm。
Second, determining the number of steps of the radiation source
The angle value with fuzzy value of two base line intervals at two frequencies can be measured by using an eight-port deviceAndn is 1,2 and indicates the nth frequency. The principle of the steps proposed by the method of the invention is explained when incidentAngle of rotationAndin the range of (-80 deg., 80 deg.), at different timesAndwithin the value range of (2), the measured angle difference values under different frequencies only contain limited fixed values called step values, and the steps are numbered in the sequence from left to right. So that the second step is to determine from the actually measured valuesAndthe number of the step.
Thirdly, restoring the actual phase value
Because each step corresponds to an angle interval, the step value can be reduced to an actual phase value by adopting a certain reduction criterion. For this purpose, a reduction criterion is established:
let LmWhere m is E, H is the number of steps, XmThe step position is numbered, the interval factor t is selected according to the specific condition of angle measurement, the manual adjustment of the 360-degree interval is carried out, and the phase of any one frequency measurement of two frequencies is selectedHaving a real phase valueIs composed of
Simulation verification is performed according to the above-mentioned reduction process, and in this example, when t is equal to 1, the blurred phase value can be reduced to the actual phase value.
Fourthly, calculating the angle of the radiation source
After the actual phase value is obtained, the actual incident angle theta can be easily calculated according to the formula (2-1)m。
The invention relates to a two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment, which is based on the principle that:
one, eight port angle measuring principle
The antenna layout represented by the red dots in fig. 5 is drawn as a perspective view, as in fig. 7, with red linesFor plane waves arriving at the antenna, thetaEIs the pitch angle of the incoming wave, thetaHIs the azimuth angle of the incoming wave, betax、βy、βzThe projections of the incoming wave on the x, y and z axes respectively.
The plane wave that is set to reach the antenna can be expressed as:
wherein t represents the time of day in which,a position vector is represented by a vector of positions,represents a wave propagation vector, so the wave propagation vector can be decomposed into:
βx=β0cosθEsinθH (13)
βy=β0cosθEcosθH (14)
βz=β0sinθE (15)
it can be seen that, assuming the receiving unit 5 as the center of the coordinate system, it can be derived:
so in the coordinate system with the receiving unit 5 as the origin, the phase of the incident plane wave at each receiving unit can be expressed as:
now assume phiE=β0dsinθEIs the phase difference between the antennas 5 and 7H=β0dcosθEsinθHIs the phase difference between the antennas 6 and 8, i.e. + -.)ECorresponding to the elevation angle phiHFor azimuth. The above equation can be written as:
the eight-port network proposed by Zhanxuchun is composed of 4 180-degree directional couplers (ring bridge) and 1 90-degree item shifter, and the concept of incident wave and reflected wave in S parameter of microwave network is utilized, and the schematic diagram of eight ports is separately drawn as shown in FIG. 8, and then the eight ports have
According to the S parameter characteristics of the eight-port network
The combination of (28) and (29) has
The square of the modulus value is taken on both sides of each equation in equation (31), the ratio of the left-hand reflected voltage to the entrance voltage of the equation becomes the ratio of the reflected power to the incident power of each port, and the converted result is simplified as follows:
and due to Pi/Pk=|Sik|2Can be combined with (32)
So that x ═ phi can be solvedE+φH) (phi) 2 and y ═ yE-φH) /2, according to phiE=β0dsinθEAnd phiH=β0dcosθEsinθHThe two angles can be solved as
θE、θHThe phase difference measured for the two baseline intervals respectively shows that the eight-port can measure the phase difference between the two baseline intervals, but the two phase differences have a fuzzy problem as the existing method at present, and therefore, a step concept is provided.
Second, step principle
Firstly, when the radiation source reaches the receiving device, the corresponding wavelengths lambda of the two frequencies1And λ2Next, the actual phase difference in both directions and the fuzzy phase value φ measured using an eight port deviceEAnd phiHAs shown in fig. 9. At this time can seeFor example, the two baseline intervals have phase ambiguity problems, and the ambiguity-free phase regions of different frequencies are different.
As can be seen from the observation formula (34), the incident angle θ is calculatedEAnd thetaHMust first find the unambiguous thetaEThe value is obtained. For this purpose, first, the angle theta is alignedEAnd performing deblurring.
First, the measured ambiguity phase value is resolved into an angle value according to equation (34), and the relationship with the actual incident angle is shown in fig. 10. It can be seen that theta is measured at two frequenciesEThe value is only equal to the true value near 0 degrees, and the fuzzy characteristic exists after a certain range is exceeded, so that the angle measurement is inaccurate.
However, when the angle values measured at the two frequencies are subtracted from each other, a step diagram as shown in fig. 11 is obtained, and the value of each step is different, that is, the correspondence between the angle difference value measured at the two different frequencies and the real angle value of the callback signal is the step shape. However, each step is not very flat and has a small deviation, that is, the actual subtracted angular difference is not a constant value in the corresponding partial interval, but the number of the steps is approximately equal, and in order to make the step effect more obvious and to determine the steps more easily, the angular difference is allowed to fluctuate within the range of ± 1, so as to obtain 7 determined steps, the heights of which are-27, 9, -36, 0, 36, -9 and 27 respectively, and the steps are numbered from left to right from 1 to 7.
The second step of the process of the present invention is therefore based on θEThe step number is determined. By means of the different values of the steps, a wide angle theta of (-80 DEG, 80 DEG) can be achievedEThe measurement is even the omnidirectional measurement of (-90 degrees and 90 degrees), and the high-precision advantage is ensured because the high-frequency test result can be adopted.
Now, the unambiguous elevation angle theta is obtainedEThe measured azimuth angle theta can be calculated according to the formula (34)H. Note θHMeasured value of (a) is specified byEInfluence of the value, so different θECorresponding to different thetaHIs measured. Due to the function cos θEAbout thetaESymmetric about 0, soHIs also related to thetaESymmetric about 0, so the equal spacing is chosen at this point toE=[-80,0]Of 17, solving for the measured theta using the second equation of equation (34)HThe values are shown in fig. 12. The corresponding difference is shown in fig. 13, where θ can be seenEThe difference maps are relatively straight step maps (or a straight line) and contain limited different step values.
At this time, though different thetaEThe value will result in a solution of θ using equation (34)HDifferent measured values are obtained, but FIG. 12 tells that a step phenomenon still exists at the moment, so that the ambiguity resolution process is the same, and the same calculation method can be adopted.
The third step of the process according to the invention is therefore based on θHThe step number is determined. By different values of the step and the ambiguity-free theta calculated in the second stepEValue, a wide angle theta of (-80 deg., 80 deg.) can be achievedEThe measurement is even the omnidirectional measurement of (-90 degrees and 90 degrees), and the high-precision advantage is ensured because the high-frequency test result can be adopted.
The two-dimensional wide-angle high-precision angle measuring method based on the eight-port four-baseline radio frequency equipment has the beneficial effects that:
high-precision merit evidence
The fuzzy angle value theta after the second step of reductionEAngle with actual value thetaEThe reduction degree was observed by comparison, and the simulation result is shown in fig. 14. It can be seen that the reduced angle curve completely coincides with the actual angle curve, indicating the angle θ reduced according to this methodEAngle theta with respect to the actualEThe method is completely equal, and the high-precision advantage is ensured due to the adoption of a high-frequency test result.
Similarly, the fuzzy angle value theta restored in the third step isHAngle with actual value thetaHThe reduction degree was observed by comparison, and the simulation result is shown in fig. 15. It can be seen that the reduced angle curve completely coincides with the actual angle curve, indicating the angle θ reduced according to this methodHAngle theta with respect to the actualHThe method is completely equal, and the high-precision advantage is ensured due to the adoption of a high-frequency test result.
Two, wide angle merit evidence
The measurement of the incident angle in the range of (-80 deg., 80 deg.) was discussed above, but it has not been proven that the method can measure the incident angle in the range of (-90 deg., 80 deg.) and (80 deg., 90 deg.). The maximum goniometric range of the method of the invention is investigated next.
Therefore, only the incident angle needs to be expanded to (-90 degrees and 90 degrees), and then the angle measurement principle is simulated, so that the incident angle in the interval of (-90 degrees and 90 degrees can be measured under a feasible frequency factor q, that is, the method can measure targets with all angles in the space, and the simulation result is shown in fig. 16 and 17.
In conclusion, the method can be realized in physical practice, and the incident angle in a wide angle range can be measured with high precision.
Examples
This section illustrates the selection of the frequency factor q and the interval factor t.
The simulation condition of the condition that the frequency factor q is 0.8 is given above, the frequency factors are assigned one by one, and the angle measurement principle of different frequency factors is simulated and verified one by one.
When q is 0.6, let λ116.6mm and λ2=10mm
The fuzzy phase theta measured at two frequencies is obtained from equation (1)EAnd the fuzzy phase theta is corrected by using the first expression of the expression (34)EResolving to obtain the angle value theta with ambiguity tested at the momentEAs shown in fig. 18.
The two measured angle values are subtracted to find the step, which is shown in fig. 19. The number of steps is 7, namely-11, 21, -32, 0, 32, -21 and 11. To restore the actual phase thetaEWhen t is 1, the reduction result is shown in fig. 20. As a result of observation, the phase θ after reduction was foundEThe two curves are completely consistent with the actual phase, and the reduction phase theta is illustratedEWith the actual phase thetaEExactly equal, ambiguous phase θEThe reduction is successful.
Reuse of reduced thetaEAnd measured thetaHValue, according to the second expression of expression (34), on the blur phase θHResolving to obtain the angle value theta with ambiguity tested at the momentHAs shown in fig. 21. The two measured angle values are subtracted to find the step, which is shown in fig. 22. At this time, different thetaEAlthough the values correspond to different orders, it is clear from the foregoing that it is entirely possible to deblur these cases with one program. To restore the actual phase thetaHWhen t is 1, the reduction result is shown in fig. 23. As a result of observation, the phase θ after reduction was foundHThe two curves are completely consistent with the actual phase, and the reduction phase theta is illustratedHWith the actual phase thetaHExactly equal, ambiguous phase θHThe reduction is successful.
Therefore, the method of the present invention is effective when the frequency factor q is 0.6.
When q is 0.3, let λ127.1mm and λ2=8.1mm
Still according to the above method, a step map is obtained as shown in fig. 24 and 26, and when t is 2, the angle value calculated by the reduction phase is also completely equal to the actual angle value of incidence, and the result is shown in fig. 25 and 27. When q is 0.3, the method of the present invention is effective.
When q is 0.1, let λ168.8mm and λ2=6.9mm
Still according to the above method, a step map is obtained as shown in fig. 28 and fig. 30, and when t is taken to be 3, the angle value calculated by the reduction phase is also completely equal to the actual incident angle value, and the result is shown in fig. 29 and fig. 31. When q is 0.1, the method of the present invention is effective.
Fourthly, summarize
And then, taking other values for q, and simulating one by one according to the process, wherein simulation results are not listed one by one.
And determining the value range of k as (0,1), scanning k at intervals of 0.1, and performing multiple times of simulation to obtain the following conclusion.
When q is equal to (0,1), the frequency factor can be measured by the method of the invention.
Further, by the generalization, the interval factor t is 3 when q ∈ (0,0.2), 2 when q ∈ (0.2,0.3), and 1 when q ∈ (0.3, 1).
Claims (9)
1. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency equipment is characterized by comprising the following steps of:
step 1, building eight-port four-baseline radio frequency equipment;
step 2, determining a baseline interval and a frequency factor according to the angle range and the wavelength of the radiation source;
step 3, determining the number of steps where the pitch angle of the radiation source is located;
step 4, restoring the phase value of the actual pitch angle;
step 5, calculating a radiation source pitch angle;
step 6, determining the number of steps where the azimuth angle of the radiation source is located;
step 7, restoring the actual azimuth angle phase value;
and 8, calculating the azimuth angle of the radiation source.
2. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 1, wherein the step 1 specifically comprises: providing radio frequency equipment, wherein two ports are arranged in each direction around the radio frequency equipment, base lines are arranged on four ports of the radio frequency equipment relative to the two directions, and the other four ports of the radio frequency equipment are respectively connected with a Schottky diode detector.
3. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 2, wherein the step 2 specifically comprises: according to the angle range and the wavelength of the dual-frequency radiation source, the ratio of two frequencies of the radiation source is determined as a frequency factor q, and any distance larger than the size of an antenna receiving the radiation source is selected as a base line interval d between every two opposite base lines.
4. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 3, wherein the step 3 specifically comprises:
step 3.1, resolving all fuzzy pitch angle values theta which can be measured under two frequencies when the radiation source is in all angles in the angle range by utilizing a first expression of the expression (34)E1And thetaE2The subscripted numbers 1 and 2 denote a first frequency and a second frequency;
in the formula (34), x is (phi)E+φH) (phi) 2 and y ═ yE-φH) [ 2 ] can be calculated from the following equation (33):
in the formula (33), Si5=Pi/P5,i=1,2,3,4;PiThe signal power output by the Schottky diode detector connected with the port i of the radio frequency equipment is P5, and the signal power received by the base line connected with the port 5 of the radio frequency equipment is P5;
step 3.2, subtracting the fuzzy pitch angle values under the two frequencies, namely thetaE1-θE2Obtaining a step map of a pitch angle, wherein the total number of steps is LEThe steps are numbered according to the relationship from small to large of the corresponding radiation source angle and are 1,2, … and LE;
Step 3.3, measuring and calculating fuzzy pitch angle difference theta 'corresponding to the actual radiation source under two baseline intervals by adopting eight-port four-baseline dual-frequency radio frequency equipment and the first equation of equation (34)'E1-θ′E2The value is the actual step value;
step 3.4, mapping the actual step value obtained in the step 3.3 and the step map obtained in the step 3.2Obtaining the actual step number XE。
5. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 4, wherein the step 4 specifically comprises:
step 4.1, selecting the phase measured by any one of the two frequenciesn is 1,2 denotes the nth frequency;
In the formula (11-1), t is an interval factor, and when q ∈ (0,0.2), the interval factor t becomes 3, when q ∈ (0.2,0.3), the interval factor t becomes 2, and when q ∈ (0.3,1), the interval factor t becomes 1.
6. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 5, wherein the step 5 specifically comprises: according to the actual phase value obtained in step 4The radiation source pitch angle theta is obtained through the formula (2-1-1)E:
7. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 6, wherein the step 6 specifically comprises:
step 6.1, utilizing the pitch angle theta of the non-fuzzy radiation source obtained in the step 5EAnd a second equation (34) for solving all ambiguous azimuth values θ that are possible to measure at two frequencies at all angles of the radiation source angular rangeH1And thetaH2The subscripted numbers 1 and 2 denote a first frequency and a second frequency;
step 6.2, subtracting the fuzzy azimuth angle values under the two frequencies, namely thetaH1-θH2To obtain the actual pitch angle thetaEStep diagram of lower azimuth angle, total number of steps is LHThe steps are numbered according to the relationship from small to large of the corresponding radiation source angle and are 1,2, … and LH;
And 6.3, measuring and calculating fuzzy azimuth angle difference theta 'corresponding to the actual radiation source under two baseline intervals by adopting an eight-port four-baseline dual-frequency radio frequency device and a second equation of the equation (34)'H1-θ′H2The value is the actual step value;
step 6.4, mapping the actual step value obtained in the step 6.3 and the step map obtained in the step 6.2 to obtain the actual step number XH。
8. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 7, wherein the step 7 specifically comprises:
step 7.1, selecting the phase of any one of the two frequencies for measurementn is 1,2 denotes the nth frequency;
In the formula (11-2), t is a section factor, the magnitude of which is determined by the frequency factor q, and the section factor t is 3 when q is equal to (0,0.2), 2 when q is equal to (0.2,0.3), and 1 when q is equal to (0.3, 1).
9. The two-dimensional wide-angle high-precision angle measurement method based on the eight-port four-baseline radio frequency device according to claim 8, wherein the step 8 specifically comprises: according to the actual phase value obtained in step 7The radiation source angle theta is obtained through the formula (2-2)H:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011040289.1A CN112162233B (en) | 2020-09-28 | 2020-09-28 | Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011040289.1A CN112162233B (en) | 2020-09-28 | 2020-09-28 | Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112162233A true CN112162233A (en) | 2021-01-01 |
CN112162233B CN112162233B (en) | 2023-07-25 |
Family
ID=73862297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011040289.1A Active CN112162233B (en) | 2020-09-28 | 2020-09-28 | Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112162233B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4845502A (en) * | 1988-04-07 | 1989-07-04 | Carr James L | Direction finding method and apparatus |
JPH02231589A (en) * | 1989-03-03 | 1990-09-13 | Koden Electron Co Ltd | Wide frequency band type direction finding device |
US5657027A (en) * | 1996-06-02 | 1997-08-12 | Hughes Electronics | Two dimensional interferometer array |
US5724047A (en) * | 1996-11-27 | 1998-03-03 | Hughes Electronics | Phase and time-difference precision direction finding system |
JP2009300161A (en) * | 2008-06-11 | 2009-12-24 | Mitsubishi Electric Corp | Radio detector |
CN105510873A (en) * | 2015-12-01 | 2016-04-20 | 中电科航空电子有限公司 | S-mode response system in amplitude comparison direction-finding system, and signal processing method |
CN109828232A (en) * | 2019-03-26 | 2019-05-31 | 北京遥感设备研究所 | A kind of array antenna angle measurement fuzzy believable degree computing device and method |
US20190302217A1 (en) * | 2018-03-28 | 2019-10-03 | Kabushiki Kaisha Toshiba | 3d direction finding method and device therefore |
CN111044971A (en) * | 2019-12-23 | 2020-04-21 | 南京长峰航天电子科技有限公司 | Two-dimensional interferometer direction finding system |
-
2020
- 2020-09-28 CN CN202011040289.1A patent/CN112162233B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4845502A (en) * | 1988-04-07 | 1989-07-04 | Carr James L | Direction finding method and apparatus |
JPH02231589A (en) * | 1989-03-03 | 1990-09-13 | Koden Electron Co Ltd | Wide frequency band type direction finding device |
US5657027A (en) * | 1996-06-02 | 1997-08-12 | Hughes Electronics | Two dimensional interferometer array |
US5724047A (en) * | 1996-11-27 | 1998-03-03 | Hughes Electronics | Phase and time-difference precision direction finding system |
JP2009300161A (en) * | 2008-06-11 | 2009-12-24 | Mitsubishi Electric Corp | Radio detector |
CN105510873A (en) * | 2015-12-01 | 2016-04-20 | 中电科航空电子有限公司 | S-mode response system in amplitude comparison direction-finding system, and signal processing method |
US20190302217A1 (en) * | 2018-03-28 | 2019-10-03 | Kabushiki Kaisha Toshiba | 3d direction finding method and device therefore |
CN109828232A (en) * | 2019-03-26 | 2019-05-31 | 北京遥感设备研究所 | A kind of array antenna angle measurement fuzzy believable degree computing device and method |
CN111044971A (en) * | 2019-12-23 | 2020-04-21 | 南京长峰航天电子科技有限公司 | Two-dimensional interferometer direction finding system |
Non-Patent Citations (3)
Title |
---|
TAKASHI MARUYAMA等: ""Four-Element Array Antenna based on Pattern Reconfigurable Yagi-Uda Antenna with Complementary Parasitic Elements "", 《2015 IEEE 4TH ASIA-PACIFIC CONFERENCE ON ANTENNAS AND PROPAGATION (APCAP)》, pages 183 - 184 * |
张旭春等: ""基于对称性的六端口测相***校准方法研究"", 《电子测量与仪器学报》, vol. 32, no. 6, pages 17 * |
田德民等: ""宽带二维多基线模拟干涉仪测向技术研究"", 《舰船电子对抗》, vol. 34, no. 5, pages 9 - 12 * |
Also Published As
Publication number | Publication date |
---|---|
CN112162233B (en) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | DOA estimation with sub-array divided technique and interporlated esprit algorithm on a cylindrical conformal array antenna | |
US6061022A (en) | Azimuth and elevation direction finding system based on hybrid amplitude/phase comparison | |
US5657027A (en) | Two dimensional interferometer array | |
CN103235292B (en) | Full-dimension and difference angle measurement method for zero setting conformal calibration of a planar phased array | |
US20220146623A1 (en) | Radar System and Vehicle | |
KR101357690B1 (en) | The calculation method of interferometer array antenna spacing ratios for direction finder | |
CN108875099B (en) | Baseline selection method based on long and short baseline interferometer direction-finding system | |
Van Doan et al. | Optimized algorithm for solving phase interferometer ambiguity | |
CN110673086A (en) | Two-dimensional angle super-resolution method based on digital array radar | |
CN108037374A (en) | A kind of array antenna near field scaling method | |
US6407702B1 (en) | Method and system for obtaining direction of an electromagnetic wave | |
CN109507635A (en) | Utilize the array amplitude phase error evaluation method of two unknown orientation auxiliary sources | |
CN106970348A (en) | Electromagnetic Vector Sensor Array decorrelation LMS two dimension MUSIC method for parameter estimation | |
Chen et al. | DOA estimation with non-uniform linear arrays: A phase-difference projection approach | |
IL257479A (en) | Radar system and method for determining direction to an object | |
CN215579031U (en) | Rotary multi-baseline phase interferometer direction-finding antenna array | |
RU2598648C1 (en) | Method for radio direction-finding and radio direction finder therefor | |
CN112162234B (en) | Wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment | |
CN112162233A (en) | Two-dimensional wide-angle high-precision angle measurement method based on eight-port four-baseline radio frequency equipment | |
Boiko et al. | Design Concepts for Mobile Computing Direction Finding Systems | |
CN114265005A (en) | Polarization phase interferometer direction finding method and device | |
Shcherbyna et al. | Accuracy characteristics of radio monitoring antennas | |
Mollai et al. | Wideband two dimensional interferometric direction finding algorithm using base-triangles and a proposed minimum planar array | |
RU2686113C1 (en) | Method of amplitude two-dimensional direction-finding | |
Ryan et al. | Plane wave spectrum scattering analysis of near-field obstacle effects on directive antenna patterns |
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 |