CN113203979A - Phase comparison method angle measurement method after smoothing - Google Patents
Phase comparison method angle measurement method after smoothing Download PDFInfo
- Publication number
- CN113203979A CN113203979A CN202110399034.2A CN202110399034A CN113203979A CN 113203979 A CN113203979 A CN 113203979A CN 202110399034 A CN202110399034 A CN 202110399034A CN 113203979 A CN113203979 A CN 113203979A
- Authority
- CN
- China
- Prior art keywords
- angle
- difference
- antennas
- signal
- azimuth
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000000691 measurement method Methods 0.000 title claims abstract description 14
- 238000009499 grossing Methods 0.000 title claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention relates to the technical field of antenna communication, in particular to a phase comparison method angle measurement method after smoothing, which comprises the following steps: s1, calculating the phase difference corresponding to the time delay difference between the signal reaching each antenna and the origin; s2 determining a received signal si (n) between the multiple antennas; s3 calculating azimuth phase differenceHorizontal direction angle A and pitching phase differenceThe pitch direction angle Ε; s4 passing through the sum of the components of azimuth angle a and elevation angle EThe existing relation, the obtained E/A value is deducedOf the value of (1), whereinThe incident direction and the Z axis form an included angle ofThe included angle between the projection of the incident direction on the XOY plane and the X axis is theta; s5 throughDetermining the position information of the space where the signal is located, and utilizing the self-adaptive beam forming principle to aim at the target userThe method can determine the position information of the space where the signal is located, aims at a target user through a self-adaptive beam forming principle, can adjust the angle of a local transmitting antenna to aim at a target signal through A and E, and solves the problem of inaccurate angle measurement of multiple antennas.
Description
Technical Field
The invention relates to the technical field of antenna communication, in particular to a phase comparison method angle measurement method after smoothing.
Background
A phase comparison antenna system compares the phases of the signals received by the two antennas to determine the coordinates of the object in a coordinate plane. In remote areas, both antennas illuminate the same spatial area, so that the signals reflected from the point targets are of substantially the same amplitude and different phase, in a two-dimensional plane for example
Let the incident angle of the incident signal be theta, the difference of the wave path of the signal reaching the two antennas is
ΔR=d sinθ
Resulting in a phase difference of
Can be calculated by the received signals of two antennas, and the signals arriving at the two antennas are processed by orthogonal down-conversion and filtering to be baseband complex signals which are respectively expressed as s1And s2Then, then
arg () represents finding the phase. To obtainThen, the incident angle theta can be reversely obtained
The method is inaccurate in angle measurement of multiple antennas, so that the invention provides a smoothed phase comparison angle measurement method.
Disclosure of Invention
Aiming at the defects of the prior art, the invention discloses a smoothed phase comparison method angle measurement method, which solves the problem of inaccurate angle measurement of multiple antennas.
The invention is realized by the following technical scheme:
a smoothed phase comparison method angle measurement method, comprising the steps of:
s1, calculating the phase difference corresponding to the time delay difference between the signal reaching each antenna and the origin;
s2 determining a received signal si (n) between the multiple antennas;
s3 calculating azimuth phase differenceHorizontal direction angle A and pitching phase differenceThe pitch direction angle Ε;
s4 passing through the sum of the components of azimuth angle a and elevation angle EThe existing relation, the obtained E/A value is deducedWherein the angle between the signal incidence direction and the Z axis isThe included angle between the projection of the incident direction on the XOY plane and the X axis is theta;
s5 throughAnd determining the position information of the space where the signal is positioned, and aiming at a target user by utilizing the self-adaptive beam forming principle.
Selecting an original point as a reference point, wherein the time delay difference between the signal arriving at each antenna and the original point is
In the formula: c is the speed of light;
a phase difference of
In the formula: f. of0For the frequency of the incident signal, λ0ρ is the radius of the center of the circle for the wavelength of the incident signal, and ρ is d/sqrt (2) when it is a 4-antenna circular array.
Further, in S3, the azimuth angle a component is located on the XOZ plane, signals arriving at the antennas T1 and T4 are simultaneous in the azimuth direction, signals arriving at the antennas T2 and T3 are also simultaneous, the difference in the azimuth direction is the difference in the path length of signals arriving at the antennas T2 and T1, and the difference in the phase due to the difference in the path length is
The estimate of the azimuthal angle A is
Further, in the S3, the pitch direction angle e component is located at the YOZ plane, signals arriving at the antennas T1 and T2 are simultaneous, signals arriving at the antennas T3 and T4 are also simultaneous, the wave path difference in the pitch direction is the wave path difference of signals arriving at the antennas T4 and T1, and the phase difference due to the wave path difference is
The estimated value of the pitch direction angle Ε is
Further, in S4, the azimuth angle a and the pitch angle e are adjusted with the antenna incident angleThe relationship between is derived as follows:
the path difference in the azimuth direction is the path difference between signals arriving at the antennas T2 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be calculated from the azimuthal angular alpha component
LΑ=d sin A
Thus is provided with
Similarly, the path difference in the pitch direction is the path difference between signals arriving at the antennas T4 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be obtained by calculating the E component of the azimuth angle E
LΕ=d sin E
Thus is provided with
further, for the calculatedSmoothing to resist the influence of noise; computation of simultaneous incident signals adaptively adjusts the computation according to SNR(s)1,s2,s3,s4) Length of the sampling point.
The invention has the beneficial effects that:
the invention is provided withThe method can determine the position information of the space where the signal is located, can aim at a target user through a self-adaptive beam forming principle, can adjust the angle of a local transmitting antenna to aim at a target signal through A and E, solves the problem of inaccurate angle measurement of multiple antennas, and has strong creativity.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a diagram of phase relationships in multiple antenna goniometry;
FIG. 2 is a diagram of a four-antenna model of incident angles;
FIG. 4 is a plot of azimuth angle A at different SNR from simulation measurements;
FIG. 5 is a plot of the pitch angle E at different SNR from simulation measurements;
FIG. 6 is a graph of fai angle information for different SNR measured in simulation;
FIG. 7 shows the angle information of theta at different SNR values obtained by simulation measurement.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment discloses a smoothened phase comparison method angle measurement method, as shown in fig. 1, the incident direction is three-dimensional, the antennas are located on an XOY plane, four antennas (T1-T4) are distributed at equal intervals to form a square, the center of the square is located at an origin O, the side length of the square is d, and the coordinates of the antennas are expressed by polar coordinates as (ρ, θ)i) In this case, the square antenna array can be regarded as a circular array,i is 1-4, and the included angle between the signal incidence direction and the Z axis isThe projection of the incident direction on the XOY plane forms an angle theta with the X axis. While also decomposing the spatial angle of incidence into an azimuth angle a (the angle of the component of the plane of the direction of incidence XOZ with the Z-axis) and a pitch angle e (the angle of the component of the plane of the direction of incidence YOZ with the Z-axis).
The multi-antenna, e.g., four-antenna four-channel model, the antenna system is composed of multiple antennas, each having a separate receive channel. The circular arrays of multiple antennas differ only in the polar coordinates (p, theta) of the antennasi) In the same manner as in the above, the difference,selecting an original point as a reference point, wherein the time delay difference between the signal arriving at each antenna and the original point is
In the formula: and c is the speed of light.
A phase difference of
In the formula: f. of0For the frequency of the incident signal, λ0ρ is the radius of the center of the circle for the wavelength of the incident signal, and ρ is d/sqrt (2) when it is a 4-antenna circular array.
Let the 4 antenna receiving channels and the baseband complex signals after the quadrature down-conversion and filtering process be respectively expressed as s1、s2、s3And s4。
The graphical illustration of the wave path difference of the horizontal and the pitching of the four channels of the four antennas is shown in fig. 2:
the azimuth angle a component is located on the XOZ plane, and the upper diagram shows that the signals arriving at the antennas T1 and T4 are simultaneous, the signals arriving at the antennas T2 and T3 are also simultaneous, the azimuth path difference is the path difference between the signals arriving at the antennas T2 and T1, and the phase difference due to the path difference is
The estimate of the azimuthal angle A is
Also, the pitch angle Ε component is located in the YOZ plane, signals arriving at the antennas T1 and T2 are simultaneous, signals arriving at the antennas T3 and T4 are simultaneous, the difference in the pitch direction is the difference in the path lengths of the signals arriving at the antennas T4 and T1, and the difference in the phase due to the difference in the path lengths is
The estimated value of the pitch direction angle Ε is
Azimuth angle a and elevation angle e and antenna incident angleThe relationship between is derived as follows: the path difference in the azimuth direction is the path difference between signals arriving at the antennas T2 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be calculated from the azimuthal angular alpha component
LΑ=d sin A
Thus is provided with
Similarly, the path difference in the pitch direction is the path difference between signals arriving at the antennas T4 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be obtained by calculating the E component of the azimuth angle E
LΕ=d sin Ε
Thus is provided with
So that at this point a and E are calculated,the true angle of incidence can then be deduced as follows:
the angle measurement process of the 4-antenna circular array or square array is as follows:
1) phase difference corresponding to time delay difference of signal arriving at each antenna and origin
3) Calculating the azimuth phase differenceThe estimated value of the horizontal direction angle A is
For more accurate subsequent calculation, the calculated values are comparedSmoothing is performed to resist the influence of noise. Computation of simultaneous incident signals adaptively adjusts the computation according to SNR(s)1,s2,s3,s4) Length of the sampling point.
4) Deriving the sum of the components of azimuth angle A and elevation angle EExisting relationshipAndso that finally the E/A value is deducedThe numerical value of (c). (the included angle between the incident direction of the signal and the Z axis isThe angle between the projection of the incident direction on the XOY plane and the X-axis is θ):
and finally, obtaining the angle position information of the target signal. By passingThe position information of the space where the signal is located can be determined, and then the target user can be aimed through an adaptive beamforming principle, for example, beamforming in MIMO. The angle of the local transmitting antenna can be adjusted by A and E to aim at the target signal, for example, the A/E aiming signal of a parabolic antenna can be adjusted.
This example shows the calculation using 1000 symbols, at different Es/N0The simulation result of the lower angle measurement adopts PSK modulation, 4 antenna incident direction anglesThe 3D graph of the time is shown in FIG. 3, and FIG. 4 is a diagram of an azimuth angle A obtained by simulation measurement under different SNR;
FIG. 5 is a plot of the pitch angle E at different SNR from simulation measurements; FIG. 6 is a graph of fai angle information for different SNR measured in simulation; FIG. 7 shows the angle information of theta at different SNR values obtained by simulation measurement.
The invention is provided withThe method can determine the position information of the space where the signal is located, can aim at a target user through a self-adaptive beam forming principle, can adjust the angle of a local transmitting antenna to aim at a target signal through A and E, solves the problem of inaccurate angle measurement of multiple antennas, and has strong creativity.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (7)
1. A smoothed phase comparison method angle measurement method, characterized in that the angle measurement method comprises the steps of:
s1, calculating the phase difference corresponding to the time delay difference between the signal reaching each antenna and the origin;
s2 determining a received signal si (n) between the multiple antennas;
s3 calculating azimuth phase differenceHorizontal direction angle A and pitching phase differenceThe pitch direction angle Ε;
s4 passing through the sum of the components of azimuth angle a and elevation angle EThe existing relation, the obtained E/A value is deducedWherein the angle between the signal incidence direction and the Z axis isThe included angle between the projection of the incident direction on the XOY plane and the X axis is theta;
2. The smoothed phase ratio method angle measurement method according to claim 1, wherein in S1, the polar coordinates are (ρ, θ)i),
Selecting an original point as a reference point, wherein the time delay difference between the signal arriving at each antenna and the original point is
In the formula: c is the speed of light;
a phase difference of
In the formula: f. of0For the frequency of the incident signal, λ0ρ is the radius of the center of the circle for the wavelength of the incident signal, and ρ is d/sqrt (2) when it is a 4-antenna circular array.
3. The smoothed phase ratio method angle measurement method according to claim 1, wherein in S3, the azimuth angle a component is located on the XOZ plane, signals arriving at the antennas T1 and T4 are simultaneous for channels in the azimuth direction, signals arriving at the antennas T2 and T3 are also simultaneous, the difference in the azimuth direction is the difference in the path of signals arriving at the antennas T2 and T1, and the phase difference due to the difference in the path is the difference in the path of signals arriving at the antennas T2 and T1
The estimate of the azimuthal angle A is
4. The smoothed phase-ratio method angle measurement method according to claim 1, wherein in the S3, the pitch angle e component is located at the YOZ plane, signals arriving at the antennas T1 and T2 are simultaneous, signals arriving at the antennas T3 and T4 are simultaneous, the path difference in the pitch direction is the path difference of signals arriving at the antennas T4 and T1, and the phase difference due to the path difference is the path difference
The estimated value of the pitch direction angle Ε is
5. The smoothed phase ratio method angle measurement method according to claim 1, wherein in S4, the azimuth and elevation direction angles Α and Ε are lower than the antenna incident anglesThe relationship between is derived as follows:
the path difference in the azimuth direction is the path difference between signals arriving at the antennas T2 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be calculated from the azimuthal angular alpha component
LΑ=d sinΑ
Thus is provided with
Similarly, the path difference in the pitch direction is the path difference between signals arriving at the antennas T4 and T1, and is determined by the angle of incidenceCan be calculated to obtain
The wave path difference can be obtained by calculating the E component of the azimuth angle E
LΕ=d sinΕ
Thus is provided with
6. the smoothed phase ratio method of angular measurement according to claim 1, characterized in thatThen, for the calculatedSmoothing to resist the influence of noise; computation of simultaneous incident signals adaptively adjusts the computation according to SNR(s)1,s2,s3,s4) Length of the sampling point.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110399034.2A CN113203979A (en) | 2021-04-14 | 2021-04-14 | Phase comparison method angle measurement method after smoothing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110399034.2A CN113203979A (en) | 2021-04-14 | 2021-04-14 | Phase comparison method angle measurement method after smoothing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113203979A true CN113203979A (en) | 2021-08-03 |
Family
ID=77026804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110399034.2A Pending CN113203979A (en) | 2021-04-14 | 2021-04-14 | Phase comparison method angle measurement method after smoothing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113203979A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120098694A1 (en) * | 2010-10-26 | 2012-04-26 | Hansen Charles T | Beam stabilization for wideband phase comparison monopulse angle estimation with electronically steered antennas |
CN108490428A (en) * | 2018-03-16 | 2018-09-04 | 西安电子科技大学 | A kind of dimensionality reduction submatrix for anti-major lobe suppression is than phase Tracking Angle Measurement method |
CN109669178A (en) * | 2018-12-25 | 2019-04-23 | 中国航天科工集团八五研究所 | A kind of spaceborne three array element pulse bidimensional direction-finding method |
CN110441729A (en) * | 2019-07-31 | 2019-11-12 | 熊军 | A kind of angle-measuring method of large-scale antenna array |
CN111458677A (en) * | 2020-03-05 | 2020-07-28 | 熊军 | Double-channel single-pulse amplitude-comparison angle measurement method and device |
-
2021
- 2021-04-14 CN CN202110399034.2A patent/CN113203979A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120098694A1 (en) * | 2010-10-26 | 2012-04-26 | Hansen Charles T | Beam stabilization for wideband phase comparison monopulse angle estimation with electronically steered antennas |
CN108490428A (en) * | 2018-03-16 | 2018-09-04 | 西安电子科技大学 | A kind of dimensionality reduction submatrix for anti-major lobe suppression is than phase Tracking Angle Measurement method |
CN109669178A (en) * | 2018-12-25 | 2019-04-23 | 中国航天科工集团八五研究所 | A kind of spaceborne three array element pulse bidimensional direction-finding method |
CN110441729A (en) * | 2019-07-31 | 2019-11-12 | 熊军 | A kind of angle-measuring method of large-scale antenna array |
CN111458677A (en) * | 2020-03-05 | 2020-07-28 | 熊军 | Double-channel single-pulse amplitude-comparison angle measurement method and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3254133B1 (en) | Direction finding using signal power | |
CN105445701B (en) | The pulse angle estimating method of DDMA MIMO radar targets | |
CN109946664A (en) | A kind of array radar target seeker Monopulse estimation method under major lobe suppression | |
JP2006270847A (en) | Antenna device | |
CN113252998B (en) | Flatness optimization method for sum and difference beam signal levels of phased array antenna | |
CN103278796B (en) | Conical surface conformal array multi-parameter joint estimation method | |
WO2019200663A1 (en) | Beam control method based on frequency diversity array antenna | |
Jiao et al. | An indoor mmwave joint radar and communication system with active channel perception | |
US11754658B2 (en) | Radio station for client localization in multipath indoor environment | |
CN107167804A (en) | A kind of sane Sidelobe Adaptive beamformer method | |
CN111693979B (en) | Digital array monopulse angle measurement method based on Taylor expansion | |
KR101873273B1 (en) | A 3-D underwater location estimating method using a half power beam width | |
CN113203979A (en) | Phase comparison method angle measurement method after smoothing | |
CN109633600B (en) | DOA estimation method of minimum redundant linear array MIMO-OTHR | |
CN111398907B (en) | Phased array radar space domain two-dimensional digital difference beam forming method | |
CN111812607A (en) | Meter-wave MIMO radar low elevation angle estimation method based on beam space | |
CN114910944A (en) | Phased array antenna tracking precision measuring method | |
CN114545448A (en) | Satellite-borne GNSS-R/S integrated receiving system and ocean intelligent detection method | |
CN108761402A (en) | Steady low-sidelobe beam forming method based on Minimum square error | |
CN114994597A (en) | Interferometer direction finding method based on uniform circular array | |
CN109633563B (en) | Self-adaptive coherent beam forming method based on multipath information | |
CN112965028A (en) | Multi-beam phased array difference sum ratio angle estimation method | |
CN109474549B (en) | Three-dimensional channel estimation method based on three-dimensional beam pattern | |
Prasad et al. | Performance analysis of direction of arrival estimation algorithms in smart antennas | |
CN111967118A (en) | Device and method for constructing circular antenna array |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210803 |