CN112782645B - Data fitting angle measurement method for four-arm helical antenna - Google Patents

Data fitting angle measurement method for four-arm helical antenna Download PDF

Info

Publication number
CN112782645B
CN112782645B CN202011631368.XA CN202011631368A CN112782645B CN 112782645 B CN112782645 B CN 112782645B CN 202011631368 A CN202011631368 A CN 202011631368A CN 112782645 B CN112782645 B CN 112782645B
Authority
CN
China
Prior art keywords
angle
estimated value
beams
target
sample data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011631368.XA
Other languages
Chinese (zh)
Other versions
CN112782645A (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.)
Wuxi Guoxin Microelectronics System Co ltd
Original Assignee
Wuxi Guoxin Microelectronics System Co ltd
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 Wuxi Guoxin Microelectronics System Co ltd filed Critical Wuxi Guoxin Microelectronics System Co ltd
Priority to CN202011631368.XA priority Critical patent/CN112782645B/en
Publication of CN112782645A publication Critical patent/CN112782645A/en
Application granted granted Critical
Publication of CN112782645B publication Critical patent/CN112782645B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems 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)

Abstract

The invention discloses a data fitting angle measurement method for a four-arm helical antenna, which relates to the technical field of radar, and comprises the steps of obtaining a first angle estimated value by utilizing a first function relation which is pre-fitted after four original antenna wave beams of a target to be measured are obtained, obtaining a pitch angle estimated value and an azimuth angle estimated value by utilizing a second function relation which is pre-fitted, obtaining a first angle estimated value and a second angle estimated value, and correcting the second angle estimated value to obtain a pitch angle and an azimuth angle of the target to be measured according to the first angle estimated value, the first angle estimated value and the second angle estimated value; the method can conveniently obtain the angle measurement result by adopting a data fitting method, saves the memory overhead and has high angle measurement efficiency.

Description

Data fitting angle measurement method for four-arm helical antenna
Technical Field
The invention relates to the technical field of radars, in particular to a data fitting angle measurement method for a four-arm helical antenna.
Background
Angle measurement is the most critical application of radar signal processing, angle measurement provides basic information for further processing and utilization of rear-end radar signals, and the two most commonly used angle measurement systems at present are single-pulse angle measurement and phase interferometer angle measurement. The phase interferometer angle measurement system obtains error signals related to the magnitude and direction of a target deviation signal axis by using the signal phase difference of continuous wave carriers of paired antenna receiving templates, and although the phase interferometer angle measurement system is high in tracking accuracy, the debugging difficulty is high, and the ambiguity resolution problem exists. Therefore, the single-pulse angle measuring system which has a simple structure and is widely applied is more commonly used at present, although the angle measuring precision is slightly worse than that of a phase interferometer.
In the single-pulse angle measurement theory, target angle error information can be determined only by one echo pulse, and in engineering, a plurality of echo pulses are usually accumulated to improve the detection probability and the angle error measurement precision. The single-pulse angle measuring system utilizes the paired wave beams to receive the amplitude or phase of a target signal and simultaneously compares the amplitude or phase to obtain an error signal related to the magnitude and direction of a target deviation equal signal axis. Therefore, the single-pulse angle measurement is further divided into a phase-comparison single-pulse angle measurement and a amplitude-comparison single-pulse angle measurement: compared with monopulse angle measurement, the method is more applied to an array antenna system with scattered feed sources, and phase differences of signals are compared to obtain an angle error value; the amplitude-comparison single-pulse angle measurement is mostly used for a surface antenna system with a concentrated feed source, and the amplitude difference of signals is compared to obtain an angle error value.
Amplitude-contrast monopulse goniometry uses multiple independent antennas to generate multiple independent contiguous beams covering 360 degrees of azimuth, these antennas using the same directional pattern function F (theta) and being evenly distributed, the field angles theta of adjacent antennass360 DEG/M, each antenna has an azimuth Fi(θ)=F(θ-iθs) i is 0,1, …, M-1, K is the number of antennas. Four antennas, six antennas and eight antennas are commonly used, the working principles of the antennas are basically similar, and the direction finding precision is improved along with the increase of the number of the antennas. Taking a four-antenna omnidirectional amplitude monopulse angle measurement system including four antennas as an example, referring to fig. 1, a four-antenna directional diagram and a corresponding digital receiver system block diagram are shown, where the field angle of each antenna is 90 °, and the respective amplitude response of the signal received by each antenna is KiThe logarithmic envelope signal of the output pulse is si=log[KiF(θ-iθs)A(t)]i is 0,1, …, M-1, where a (t) is the amplitude modulation of the radar signal, and the logarithmic envelope signal is fed to a signal processor, which generates a corresponding angle estimate. Specifically, the method comprises the following steps: in the system, a pair of adjacent beams respectively output strongest and second strongest signals, and the radar azimuth can be determined by comparing the relative magnitude of the envelope amplitudes of the signals output by the pair of adjacent beams. As shown in fig. 2, assuming that the antenna pattern satisfies the symmetry of the amplitude direction shown in fig. 2, i.e., F (θ) ═ F (- θ), the radar direction is located between the two antennas and is offset from the signal direction of the antennas, etc., by the angle
Figure BDA0002876627880000021
When the corresponding channel output signals are respectively
Figure BDA0002876627880000022
I.e. after subtraction, the logarithmic voltage ratio in decibels R is
Figure BDA0002876627880000023
If the F (theta) function is in the interval [ -theta ]ss]Has monotonicity in the inside, namely satisfies F (theta)1)<F(θ2)
Figure BDA0002876627880000024
θ12∈[-θss]Then the logarithmic voltage ratio R to
Figure BDA0002876627880000025
The angle information of the target can be obtained according to the size of R.
The planar helical antenna can determine angle information through a beam relationship between the symmetrical arms, the four-arm helix can determine the azimuth and the elevation through the relationship between the two symmetrical arms so as to achieve the purpose of angle measurement, and the planar helical antenna is widely applied to the fields of GPS, PCS and the like due to the characteristics of wide beams, wide frequency bands, circular polarization and the like. In a passive direction-finding system based on a quadrifilar helical antenna, angle measurement can be performed by adopting the method, but for an actual angle-finding system, the specific calculation R and R are caused by the fact that the directional diagram function of the antenna is complicated
Figure BDA0002876627880000026
The function between the two is not an analytic solution, the relation between the angle of each direction and the beam size of the radar antenna is obtained through experiments in engineering application, a specific lookup table is manufactured, and the corresponding angle is obtained through the lookup table in the actual use process. However, the table lookup method has two disadvantages: 1, table data occupies a large number of memory units; 2, the comparison of the lookup method wastes the lookup time, the lookup algorithm is not well executed, and the lookup time is differentStrain is slow which also presents certain challenges to real-time processing systems.
Disclosure of Invention
The invention provides a data fitting angle measurement method for a four-arm helical antenna aiming at the problems and the technical requirements, and the technical scheme of the invention is as follows:
a data fitting angle measurement method for a four-arm helical antenna comprises the following steps:
carrying out direction finding on sample targets at different angles by using a passive direction finding system based on a four-arm helical antenna, and determining the target direction of the sample targets in an antenna coordinate system and four corresponding original antenna beams acquired through the four-arm helical antenna in each direction finding to obtain a group of sample data;
fitting based on sample data to obtain a first functional relation, wherein the first functional relation is a functional relation between a first angle and four original antenna beams, and the first angle is an included angle between a target incoming direction and a horizontal plane of an antenna coordinate system;
fitting based on sample data to obtain a second functional relation, wherein the second functional relation is a functional relation between the pitch angle and the azimuth angle which are determined based on the target direction and the four original antenna beams;
acquiring four original antenna beams of a target to be measured by using a passive direction finding system based on a four-arm helical antenna;
obtaining a first angle first estimated value based on a first functional relation by utilizing four original antenna beams of a target to be measured;
obtaining a pitch angle estimated value and an azimuth angle estimated value based on a second function relation by utilizing four original antenna beams of the target to be detected, and obtaining a first angle second estimated value and a second angle estimated value based on the pitch angle estimated value and the azimuth angle estimated value, wherein the second angle estimated value is an included angle between the projection of the target to be detected on the horizontal plane and a preset direction;
and correcting the pitch angle pre-estimated value and the azimuth angle pre-estimated value according to the first angle first pre-estimated value, the first angle second pre-estimated value and the second angle pre-estimated value to obtain the pitch angle and the azimuth angle of the target to be detected.
The further technical scheme is that the pitch angle and azimuth angle of the target to be measured are obtained according to the first angle estimated value, the second angle estimated value and the second angle estimated value, and the pitch angle estimated value and the azimuth angle estimated value are corrected, and the method comprises the following steps:
when the difference value between the first angle first estimated value and the first angle second estimated value reaches a preset threshold value, correcting the pitch angle estimated value and the azimuth angle estimated value by using the first angle first estimated value and the second angle estimated value to obtain a pitch angle and an azimuth angle of the target to be measured;
and when the difference value between the first angle first estimated value and the first angle second estimated value does not reach a preset threshold value, taking the pitch angle estimated value as the pitch angle of the target to be detected, and taking the azimuth angle estimated value as the azimuth angle of the target to be detected.
The further technical scheme is that the pitch angle estimation value and the azimuth angle estimation value are corrected by utilizing the first angle estimation value and the second angle estimation value to obtain the pitch angle and the azimuth angle of the target to be measured, and the method comprises the following steps:
corrected pitch angle estimate as β tan-1(cosφ0tanθ1) Obtaining the pitch angle of the target to be measured, and correcting the estimated value of the azimuth angle to alpha sin-1(sinφ0sinθ1) Obtaining the azimuth angle phi of the target to be measured0For a second angle estimate, theta1A first estimate is made for the first angle.
The further technical scheme is that a first angle second estimated value and a second angle estimated value are obtained based on a pitch angle estimated value and an azimuth angle estimated value, and the method comprises the following steps:
determining the first angle and the second estimated value as theta2=cos-1(cos(β0)cos(α0) Determine a second angle estimate as phi0=tan-1(tan(β0)/sin(α0) Wherein α) is0For azimuthal prediction, beta0Is an estimate of pitch angle.
According to a further technical scheme, a first angle of a sample target is indicated to a target of the sample target in each group of sample data, and a first functional relation is obtained based on sample data fitting, wherein the method comprises the following steps:
for each sample data, calculating sum beams and difference beams of four original antenna beams of the sample data, and calculating the ratio of the difference beams to the sum beams to obtain difference mode sum mode ratios;
and fitting to obtain a first functional relation based on the corresponding relation between the first angle and the difference mode and the mode ratio in each group of sample data, wherein the first function reflects the functional relation between the first angle and the difference mode and the mode ratio.
According to a further technical scheme, the pitch angle and the azimuth angle of the sample target are indicated to the target of the sample target in each group of sample data, and a second functional relation is obtained based on sample data fitting, and the method comprises the following steps:
for each sample data, calculating sum beams and difference beams of four original antenna beams of the sample data, and performing different beam operations on the sum beams and the difference beams respectively to obtain four off-axis orthogonal beams;
determining a first parameter according to one set of two diagonal off-axis orthogonal beams and a sum beam, and determining a second parameter according to the other set of two diagonal off-axis orthogonal beams and the sum beam;
and fitting based on the sample data to obtain a functional relation between the azimuth angle and the first parameter, and fitting to obtain a functional relation between the pitch angle and the second parameter to obtain a second functional relation.
The further technical scheme is that a functional relation between the azimuth angle and the first parameter is obtained through fitting based on sample data, and a functional relation between the pitch angle and the second parameter is obtained through fitting, and the method comprises the following steps: according to a ═ b0+b1·k1Is fitted to obtain a functional relationship between the azimuth angle and the first parameter, in terms of β ═ c0+c1·k2Is fitted to obtain a functional relationship between the pitch angle and the second parameter, where α is the azimuth angle, b0、b1Is a coefficient, k1Is a first parameter, beta is the pitch angle, c0、c1Is a coefficient, k2Is the second parameter.
A further technical solution is that, determining a first parameter according to one set of two diagonal off-axis orthogonal beams and a sum beam, and determining a second parameter according to the other set of diagonal two off-axis orthogonal beams and the sum beam, including:
determining a first parameter as
Figure BDA0002876627880000041
Determining the second parameter as
Figure BDA0002876627880000042
Wherein, B1And B3Is a set of two diagonal off-axis orthogonal beams, B2And B4Is another set of two diagonally opposite, off-axis orthogonal beams, M1Is a sum beam, the symbol | | | represents the amplitude of the beam.
The method further adopts the technical scheme that the method comprises the steps of calculating sum beams and difference beams of four original antenna beams of sample data and respectively carrying out different beam operations on the sum beams and the difference beams to obtain four off-axis orthogonal beams, and comprises the following steps:
determining a sum beam as M1P1-j P2-P3+ j P4, difference beam M2P1-P2+ P3-P4, P1, P2, P3, P4 are four original antenna beams, respectively;
determining four off-axis orthogonal beams B1、B2、B3、B4Are respectively as
Figure BDA0002876627880000043
Wherein, B1And B3As a set of two diagonal off-axis orthogonal beams, B2And B4Is another set of two off-axis orthogonal beams that are diagonal.
The beneficial technical effects of the invention are as follows:
the application discloses a data fitting angle measurement method for a four-arm helical antenna, which adopts a data fitting algorithm to reduce the problems of overlarge memory occupation and low table look-up efficiency caused by a table look-up method; the fitting algorithm can conveniently obtain angle measurement results, meanwhile, the stored parameters are few, the memory cost is saved, the two major defects caused by the current table look-up method structure are overcome, the method can be widely applied to passive seeker antennas, the four-arm helical antenna is adopted for angle measurement, and the frequency band range can be 580M-2G.
Drawings
Fig. 1 is a four-antenna omnidirectional amplitude monopulse measurement schematic block diagram.
Fig. 2 is an amplitude pattern of adjacent antennas.
Fig. 3 is a flowchart of a method of the present application for a four-arm helical antenna data fitting goniometry method.
Detailed Description
The following further describes the embodiments of the present invention with reference to the drawings.
The application discloses a data fitting angle measurement method for a quadrifilar helix antenna, please refer to fig. 3, the method comprises the following steps:
step S1, performing direction finding on the sample target at different angles by using a passive direction finding system based on the quadrifilar helix antenna, and determining a target arrival direction of the sample target in the antenna coordinate system and corresponding four original antenna beams P1, P2, P3, and P4 acquired through the quadrifilar helix antenna in each direction finding to obtain a set of sample data. In the present application, the target direction of the sample target indicates the pitch angle β and the azimuth angle α of the sample target, and also indicates a first angle θ and a second angle Φ of the sample target, the first angle is an angle between the target direction and a horizontal plane of the antenna coordinate system, that is, an XY plane, and the second angle is an angle between a projection of the target direction onto the horizontal plane, that is, the XY plane, and a predetermined direction, which is predefined in advance, and usually a direction of an X axis or a Y axis of the antenna coordinate system is taken as the predetermined direction.
The (α, β) and (θ, Φ) can be converted into each other by coordinate conversion, so that the target orientation of the sample target can actually include only one set of angles of (α, β) and (θ, Φ), and the other set of angles can be converted accordingly. Or the target direction of arrival may actually include both (α, β) and (θ, Φ), for example, a set of sample data obtained at each direction finding may be represented as { (α, β), (θ, Φ), (P1, P2, P3, P4) }.
Step S2, fitting to obtain a first functional relation based on the sample data, wherein the first functional relation is a first angle thetaAs a function of the four original antenna beams. Alternatively, in the present application, the first functional relationship does not directly reflect the functional relationship between the first angle θ and the four original antenna beams, but reflects the first angle θ and the differential mode-sum mode ratio
Figure BDA0002876627880000061
Functional relationship between them, i.e.
Figure BDA0002876627880000062
f1() Is a first functional relationship, wherein M1Is a sum beam, M2For a difference beam, the symbol | | | represents the amplitude of the beam. Thus, in fitting to obtain the first functional relationship, first for each sample data, the sum beam M of the four original antenna beams P1, P2, P3, P4 of the sample data is calculated1Sum and difference beam M2The calculation method is M1=P1-j*P2-P3+j*P4,M2P1-P2+ P3-P4, and then calculate the difference beam M2And beam M1The ratio of (A) to (B) yields a difference mode and a mode ratio
Figure BDA0002876627880000063
And then based on the first angle theta and the difference mode and the mode ratio in each group of sample data
Figure BDA0002876627880000064
Fitting the corresponding relation to obtain a first functional relation f1()。
And step S3, fitting based on the sample data to obtain a second functional relation, wherein the second functional relation is a functional relation between the pitch angle beta and the azimuth angle alpha which are determined based on the target direction and the four original antenna beams. Similarly, the present application does not directly determine the functional relationship between β, α and the four original antenna beams, but first computes, for each sample data, the sum beam M of the four original antenna beams P1, P2, P3, P4 of the sample data1Sum and difference beam M2The calculation method is as described in step S2 above. Then respectively to sum beams M1Sum and difference beam M2Carrying out different beam operations to obtain four off-axis orthogonal beams B1、B2、B3、B4The calculation method comprises the following steps:
Figure BDA0002876627880000065
wherein, B1And B3As a set of two diagonal off-axis orthogonal beams, B2And B4Is another set of two off-axis orthogonal beams that are diagonal. Two off-axis orthogonal beams B from one set of diagonals1And B3And beam M1Determining a first parameter k1From another set of diagonal two off-axis orthogonal beams B2And B4And beam M1Determining a second parameter k2Specifically, the method comprises the following steps: the first parameter is
Figure BDA0002876627880000066
Determining the second parameter as
Figure BDA0002876627880000067
Therefore, each group of sample data can be processed to obtain a target to determine the pitch angle beta and the azimuth angle alpha and the corresponding first parameter k1And a second parameter k2
Then fitting based on sample data to obtain an azimuth angle alpha and a first parameter k1The pitch angle beta and the second parameter k are obtained by functional relation between the two and fitting2And obtaining a second functional relation through the functional relation between the first and second groups. The application uses polynomial fitting, i.e. according to α ═ b0+b1·k1The azimuth angle alpha and the first parameter k are obtained by form fitting1According to β ═ c0+c1·k2Form fitting to obtain a pitch angle beta and a second parameter k2Functional relationship between b0、b1Is a coefficient of, c0、c1Is a coefficient.
And step S4, acquiring four original antenna beams of the target to be measured by using the passive direction finding system based on the four-arm helical antenna. Using four original antenna beams of the target to be measured based on a first functionThe relationship is used to obtain a first estimated value of the first angle, and similar to the fitting process, the difference mode and the mode ratio are obtained by processing four original antenna beams
Figure BDA0002876627880000071
Then substituting the first function relation to obtain a corresponding first angle first estimated value theta1
Step S5, obtaining a pitch angle estimated value beta based on a second functional relation by using four original antenna wave beams of the target to be measured0Sum azimuth angle estimate α0Similar to the fitting process described above, the first parameter k is obtained by processing four original antenna beams1And a second parameter k2Then substituting the second function relationship to obtain a pitch angle estimated value beta0Sum azimuth angle estimate α0
Estimation value beta based on pitch angle0Sum azimuth angle estimate α0Obtaining a second estimated value theta of the first angle2And a second angle estimate phi0Obtained according to the following formula: determining the first angle and the second estimated value as theta2=cos-1(cos(β0)cos(α0) Determine a second angle estimate as phi0=tan-1(tan(β0)/sin(α0))
Step S6, a first estimated value theta is calculated according to the first angle1A second estimated value theta of the first angle2And a second angle estimate phi0Corrected pitch angle estimate beta0Sum azimuth angle estimate α0And obtaining the pitch angle and the azimuth angle of the target to be detected. Specifically, the method comprises the following steps:
(1) when the first angle is the first estimated value theta1And a second estimated value theta of the first angle2When the difference between the first and second angles reaches a predetermined threshold T, a first predicted value theta is used1And a second angle estimate phi0Corrected pitch angle estimate beta0Sum azimuth angle estimate α0And obtaining the pitch angle and the azimuth angle of the target to be detected. The predetermined threshold is a user-defined value, and the value of the predetermined threshold is typically 0.02.
Using theta1And phi0The method for obtaining the pitch angle and the azimuth angle through correction comprises the following steps: corrected pitch angle estimate as β tan-1(cosφ0tanθ1) Obtaining the pitch angle of the target to be measured, and correcting the estimated value of the azimuth angle to alpha sin-1(sinφ0sinθ1) And obtaining the azimuth angle of the target to be measured.
(2) When the first angle is the first estimated value theta1And a second estimated value theta of the first angle2When the difference between the pitch angle and the pitch angle does not reach the predetermined threshold value T, the pitch angle estimated value beta obtained in step S5 is directly used0As the pitch angle of the target to be measured, the azimuth angle estimated value alpha is estimated0As the azimuth of the object to be measured.
What has been described above is only a preferred embodiment of the present application, and the present invention is not limited to the above embodiment. It is to be understood that other modifications and variations directly derivable or suggested by those skilled in the art without departing from the spirit and concept of the present invention are to be considered as included within the scope of the present invention.

Claims (6)

1. A method for measuring an angle by data fitting of a four-arm helical antenna is characterized by comprising the following steps:
carrying out direction finding on sample targets at different angles by using a passive direction finding system based on a four-arm helical antenna, and determining the target direction of the sample targets in an antenna coordinate system and four corresponding original antenna beams acquired through the four-arm helical antenna in each direction finding to obtain a group of sample data;
fitting based on sample data to obtain a first functional relation, wherein the first functional relation is a functional relation between a first angle and four original antenna beams, and the first angle is an included angle between a target incoming direction and a horizontal plane of the antenna coordinate system;
fitting based on sample data to obtain a second functional relation, wherein the second functional relation is a functional relation between the pitch angle and the azimuth angle which are determined based on the target direction and the four original antenna beams;
acquiring four original antenna beams of a target to be measured by using a passive direction finding system based on a four-arm helical antenna;
obtaining a first angle first estimated value based on the first functional relation by utilizing four original antenna beams of the target to be measured;
obtaining a pitch angle estimated value and an azimuth angle estimated value based on the second functional relation by utilizing four original antenna beams of the target to be detected, and obtaining a first angle second estimated value and a second angle estimated value based on the pitch angle estimated value and the azimuth angle estimated value, wherein the second angle estimated value is an included angle between the projection of the target to be detected on the horizontal plane and a preset direction; wherein the second estimated value of the first angle is theta2=cos-1(cos(β0)cos(α0) The second angle estimate is phi)0=tan-1(tan(β0)/sin(α0) Wherein α) is0For the estimation of said azimuth angle, β0Estimating a pitch angle for the vehicle;
correcting the pitch angle estimated value and the azimuth angle estimated value according to the first angle estimated value, the second angle estimated value and the second angle estimated value to obtain a pitch angle and an azimuth angle of the target to be detected, and the method comprises the following steps: correcting the pitch angle estimate to β tan when a difference between the first angle first estimate and the first angle second estimate reaches a predetermined threshold-1(cosφ0tanθ1) Obtaining the pitch angle of the target to be measured, and correcting the estimated azimuth angle to alpha sin-1(sinφ0sinθ1) Obtaining the azimuth angle theta of the target to be measured1A first estimated value for the first angle; and when the difference value between the first angle first estimated value and the first angle second estimated value does not reach the preset threshold value, taking the pitch angle estimated value as the pitch angle of the target to be detected, and taking the azimuth angle estimated value as the azimuth angle of the target to be detected.
2. The method of claim 1, wherein the target of the sample target in each set of sample data is at a first angle indicative of the sample target, then said fitting based on the sample data results in a first functional relationship comprising:
for each sample data, calculating sum beams and difference beams of four original antenna beams of the sample data, and calculating the ratio of the difference beams to the sum beams to obtain difference mode sum mode ratios;
and fitting to obtain the first functional relation based on the corresponding relation between the first angle and the difference mode and the mode ratio in each group of sample data, wherein the first function reflects the functional relation between the first angle and the difference mode and the mode ratio.
3. The method of claim 1, wherein the target of the sample target in each set of sample data indicates a pitch angle and an azimuth angle of the sample target, and said fitting based on the sample data yields a second functional relationship comprising:
for each sample data, calculating sum beams and difference beams of four original antenna beams of the sample data, and performing different beam operations on the sum beams and the difference beams respectively to obtain four off-axis orthogonal beams;
determining a first parameter from one set of two diagonal off-axis orthogonal beams and the sum beam, and determining a second parameter from the other set of two diagonal off-axis orthogonal beams and the sum beam;
and fitting based on sample data to obtain a functional relation between the azimuth angle and the first parameter and fitting to obtain a functional relation between the pitch angle and the second parameter to obtain the second functional relation.
4. The method of claim 3, wherein fitting a functional relationship between an azimuth angle and a first parameter and fitting a functional relationship between a pitch angle and a second parameter based on sample data comprises: according to a ═ b0+b1·k1Is fitted to obtain a functional relationship between the azimuth angle and the first parameter, in terms of β ═ c0+c1·k2Is fitted to obtain a functional relationship between the pitch angle and the second parameter, where α is the azimuth angle, b0、b1Is a coefficient, k1Is the firstParameter, beta is pitch angle, c0、c1Is a coefficient, k2Is the second parameter.
5. The method of claim 3, wherein determining a first parameter from one of the two diagonal off-axis orthogonal beams and the sum beam and determining a second parameter from the other of the two diagonal off-axis orthogonal beams and the sum beam comprises:
determining the first parameter as
Figure FDA0003283469450000021
Determining the second parameter as
Figure FDA0003283469450000022
Wherein, B1And B3Is a set of two diagonal off-axis orthogonal beams, B2And B4Is another set of two diagonally opposite, off-axis orthogonal beams, M1Is the sum beam, the symbol | | | represents the amplitude of the beam.
6. The method of claim 3, wherein said calculating a sum beam and a difference beam of four original antenna beams of the sample data and performing different beam operations on the sum beam and the difference beam to obtain four off-axis orthogonal beams comprises:
determining a sum beam as M1P1-j P2-P3+ j P4, difference beam M2P1-P2+ P3-P4, P1, P2, P3, P4 are four original antenna beams, respectively;
determining four off-axis orthogonal beams B1、B2、B3、B4Are respectively as
Figure FDA0003283469450000031
Wherein, B1And B3As a set of two diagonal off-axis orthogonal beams, B2And B4Is another set of two off-axis orthogonal beams that are diagonal.
CN202011631368.XA 2020-12-30 2020-12-30 Data fitting angle measurement method for four-arm helical antenna Active CN112782645B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011631368.XA CN112782645B (en) 2020-12-30 2020-12-30 Data fitting angle measurement method for four-arm helical antenna

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011631368.XA CN112782645B (en) 2020-12-30 2020-12-30 Data fitting angle measurement method for four-arm helical antenna

Publications (2)

Publication Number Publication Date
CN112782645A CN112782645A (en) 2021-05-11
CN112782645B true CN112782645B (en) 2021-11-19

Family

ID=75753325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011631368.XA Active CN112782645B (en) 2020-12-30 2020-12-30 Data fitting angle measurement method for four-arm helical antenna

Country Status (1)

Country Link
CN (1) CN112782645B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640755A (en) * 2021-05-24 2021-11-12 中国南方电网有限责任公司超高压输电公司广州局 Target pitch angle acquisition method and device based on radar photoelectric linkage system
CN113884977B (en) * 2021-08-17 2023-09-08 中国电子科技集团公司第二十九研究所 Correction method for one-dimensional interferometer direction-finding conical effect
CN114234904B (en) * 2021-12-16 2023-01-24 四川九洲电器集团有限责任公司 Angle measuring method, system and medium based on lens multi-beam antenna
CN116540219B (en) * 2023-07-04 2023-09-22 北醒(北京)光子科技有限公司 Method and device for correcting radar emergent light angle, storage medium and electronic equipment

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104198985A (en) * 2014-09-24 2014-12-10 中国电子科技集团公司第五十四研究所 Single antenna direction-finding method
CN105445701A (en) * 2015-11-11 2016-03-30 西安电子科技大学 Mono-pulse angle estimation method for DDMA-MIMO radar target
CN107656238A (en) * 2017-09-12 2018-02-02 中国海洋大学 A kind of high-frequency ground wave radar direction finding new method based on double LUT Methods
CN109581279A (en) * 2018-11-06 2019-04-05 中国电子科技集团公司第三十八研究所 A kind of bearing calibration of ultra wide band multiple-beam system direction finding precision and device
CN109856605A (en) * 2019-01-11 2019-06-07 中国船舶重工集团公司第七二四研究所 A kind of while formation of the digital multiple beam quadratic fit curve is directed toward modification method
CN111398914A (en) * 2020-04-23 2020-07-10 中国北方工业有限公司 Radar tracking beam angle measurement method
CN111693979A (en) * 2020-05-25 2020-09-22 北京理工大学 Digital array monopulse angle measurement method based on Taylor expansion

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140002306A1 (en) * 2012-03-07 2014-01-02 Son Thach Nguyen Directional radio signal detection apparatus and methods of use

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104198985A (en) * 2014-09-24 2014-12-10 中国电子科技集团公司第五十四研究所 Single antenna direction-finding method
CN105445701A (en) * 2015-11-11 2016-03-30 西安电子科技大学 Mono-pulse angle estimation method for DDMA-MIMO radar target
CN107656238A (en) * 2017-09-12 2018-02-02 中国海洋大学 A kind of high-frequency ground wave radar direction finding new method based on double LUT Methods
CN109581279A (en) * 2018-11-06 2019-04-05 中国电子科技集团公司第三十八研究所 A kind of bearing calibration of ultra wide band multiple-beam system direction finding precision and device
CN109856605A (en) * 2019-01-11 2019-06-07 中国船舶重工集团公司第七二四研究所 A kind of while formation of the digital multiple beam quadratic fit curve is directed toward modification method
CN111398914A (en) * 2020-04-23 2020-07-10 中国北方工业有限公司 Radar tracking beam angle measurement method
CN111693979A (en) * 2020-05-25 2020-09-22 北京理工大学 Digital array monopulse angle measurement method based on Taylor expansion

Also Published As

Publication number Publication date
CN112782645A (en) 2021-05-11

Similar Documents

Publication Publication Date Title
CN112782645B (en) Data fitting angle measurement method for four-arm helical antenna
CN109669178B (en) Satellite-borne three-array-element single-pulse two-dimensional direction finding method
US20220113363A1 (en) Direction of arrival estimation
CN109856605B (en) Quadratic fitting curve pointing correction method for simultaneous digital multi-beam forming
US20050270229A1 (en) Positioning system with a sparse antenna array
CN112083387B (en) Radar calibration method and device
US6169519B1 (en) TCAS bearing measurement receiver apparatus with phase error compensation method
CN104007413A (en) Array position error correction method taking information source azimuth error into account
CN113189625B (en) Error correction method and device based on single-satellite interferometer direction finding system and satellite
CN113884977A (en) Correction method for one-dimensional interferometer direction finding cone effect
CN116087871A (en) Direction finding method, device and system for circular array interferometer and storage medium
CN110927751A (en) Array antenna self-adaptive correction implementation method based on carrier phase measurement
CN114236465A (en) Search type two-dimensional high-precision amplitude comparison direction finding method and device
KR101796472B1 (en) Radar apparatus and DOA estimation method using the same
CN110515066B (en) Vehicle-mounted millimeter wave radar and target height measuring method thereof
RU2711341C1 (en) Two-dimensional direction finding method
JP5722026B2 (en) Direction of arrival estimation method
CN116091782A (en) Three-dimensional feature extraction method for space target component
KR20160019673A (en) A radio direction finder
US11949488B2 (en) Method and apparatus for realizing beam alignment
CN112986701B (en) Holographic measurement method and system based on radio frequency power supply broadband signal
JP7258264B2 (en) Angular Width Estimating Device, Angular Width Estimating Method, and Target Angle Measuring Device
RU2686113C1 (en) Method of amplitude two-dimensional direction-finding
CN112946615A (en) Amplitude-phase error correction method for phased array system
CN111123250A (en) Pulse Doppler radar based on pattern search algorithm and beam forming method

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