CN105301564A - Phased array antenna beam pointing error compensation method based on bilinear fitting - Google Patents
Phased array antenna beam pointing error compensation method based on bilinear fitting Download PDFInfo
- Publication number
- CN105301564A CN105301564A CN201510581867.5A CN201510581867A CN105301564A CN 105301564 A CN105301564 A CN 105301564A CN 201510581867 A CN201510581867 A CN 201510581867A CN 105301564 A CN105301564 A CN 105301564A
- Authority
- CN
- China
- Prior art keywords
- beam pointing
- wave
- pointing error
- error
- phased array
- 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 35
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000000007 visual effect Effects 0.000 claims description 6
- 208000035126 Facies Diseases 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000013139 quantization Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a phased array antenna beam pointing error compensation method based on bilinear fitting. A bilinear fitting method which is easily realized by a wave control machine is used to divide a curved surface formed by four adjacent beam position errors in a space into two planes. In each plane, two space points are used as a straight line to calculate error slope. Beam pointing error slopes in azimuth and elevation are respectively calculated. According to the slopes, the beam pointing error of a beam position between given beam position points can be calculated. The method is used to calculate the beam pointing error of the whole two-dimensional visible area. According to the calculated beam pointing error, a phased array antenna beam compensation table is generated and is written into the wave control machine to realize phased array antenna online beam pointing error compensation. The phased array antenna beam pointing accuracy is improved. According to the bilinear interpolation fitting method, the number of the beam positions, which need to be measured, of a phased array antenna can be reduced; the measurement test time is saved; the efficiency is significantly improved; and the cost is reduced.
Description
Technical Field
The invention relates to the technical field of phased array antennas, in particular to a beam pointing error compensation method of a phased array antenna based on bilinear fitting.
Background
The phased array radar seeker adopts a phased array antenna system, and can transfer a wave beam from one direction to the other direction within microsecond-level time by means of a digital phase shifter, so that fast scanning of an airspace and effective tracking of a target are realized, and a missile is guided to fly to the target. Because the digital phase shifter cannot continuously adjust the phase, and errors are introduced due to design processes and the like, the antenna beam cannot accurately point to a desired direction, and pointing deviation is caused. The size of the pointing error of the phased array antenna has very important influence on a phased array radar seeker system, and the accurate hitting capability of the phased array radar seeker is directly limited.
The phased array antenna has the inherent characteristic of beam pointing error, and the beam pointing error is different at different pointing angles, so that the precision of a subsequent guidance control instruction is seriously influenced. Because the factors influencing the phased array antenna are many, including system errors and random errors, the improvement of the beam pointing accuracy is a difficult point of the phased array antenna, and in order to improve the beam pointing accuracy of the phased array antenna and improve the technical and combat performance of a phased array radar seeker system, the phased array antenna pointing error needs to be effectively compensated.
Research on improving the beam pointing accuracy of phased array antennas is focused on by researchers in the technical field of phased arrays, and researchers in this technical field at home and abroad propose various methods for improving the beam pointing accuracy, and table 1 is related domestic and foreign patents and papers.
Table 1 phased array antenna beam pointing accuracy related patent and thesis situation table
FRAZITARICHARDF et al propose a method for reducing antenna beam pointing errors due to phase quantization of phase shifters using a coupling network structure, which requires an antenna array structure having a corresponding structure with a phase shift that is an odd integer multiple of half of the minimum phase step of the phase shifters of the antenna elements located on both sides of the center of the array, and have filed patent application with US4188633 (a). The method has strict requirements on the antenna array structure and is not suitable for common phased array antennas. CHANGKAICHIANG discloses a method for quantizing phase associated with phased array radar antenna sections, which reduces beam pointing error by performing digital randomized assignment on each element, and simultaneously reduces beam sidelobe peak value, and the method requires the use of distributed beam controllers, designs a random phase adjustment algorithm in each beam controller, and applies for US5103232 (a). The phase error method is used for improving the beam pointing accuracy, and the mean value of the phase errors which can not be accurately phase-shifted in the digital phase shifter is equal to zero, so that a group of carry probabilities and truncation probabilities of random phase feeds are obtained. The set of probabilities is used to perform random phase feeding for the phased array antenna. The saphen saki et al developed phased array antenna pointing accuracy test studies using statistical methods. Liuxuan et al utilize the method of the cyclic extrapolation array element phase feed, determine the position of the central array element according to the difference between the ideal phase feed value and the quantization step value, and feed the phase quantization step value to the phase shifter with smaller absolute value, and extrapolate in turn, and complete the phase feed of all the array elements under the condition of ensuring the minimum accumulated phase error value of all the array elements.
The methods either put requirements on the antenna array structure or carry out random phase feeding by using a probabilistic method, and the methods have certain probability, and require random phase feeding calculation on each antenna unit, so that algorithm divergence sometimes occurs, beam pointing error is increased, and the calculation amount is large, and actual engineering application is difficult.
Disclosure of Invention
The invention provides a bilinear fitting-based beam pointing error compensation method for a phased array antenna, which can calculate pointing error values of all wave positions according to beam pointing errors of limited wave positions of the phased array antenna, and can compensate the pointing error of the phased array antenna according to the pointing error values so as to improve the pointing accuracy of the phased array antenna beam.
In order to achieve the above object, the present invention provides a method for compensating a beam pointing error of a phased array antenna based on bilinear fitting, comprising the following steps:
step S1, discretizing a visual angle area;
given a scan angle region of a phased array antennaAt azimuthal angular intervalsDiscretizing the visual angle region by the discrete angle formed by the pitch angle interval delta theta to obtain a spatial discrete anglei＝1,2,…,I,k＝1,2,…,K；
Taking the discrete angle as an ideal wave site, and forming a wave site grid by any 4 adjacent wave sites;
step S2, obtaining an ideal beam pointing angle of a wave position point on the wave position grid;
obtaining ideal wave beam pointing angles of all wave sites according to theoretical calculation, randomly selecting 4 adjacent space wave sites, and pointing the ideal wave beam of the wave site 1 to beThe ideal beam pointing at wave position 2 isThe ideal beam pointing at wave position point 3 isThe ideal beam pointing at wave position 4 is
Step S3, obtaining the actual beam pointing angle of the wave position point on the wave position network;
combining with an actual test system, scanning antenna beams through a wave control machine to obtain actual beam pointing angles of all wave positions, wherein the actual beam pointing direction of the wave position 1 is (α)1,β1) The actual beam pointing at wave position 2 is (α)2,β1) The actual beam pointing at wave position 3 is (α)1,β2) The actual beam pointing direction 4 of the wave site is (α)2,β2)；
Step S4, calculating the beam pointing error of the wave position on the wave position grid;
ideal beam according to all wave positionsThe pointing angle and the corresponding actual beam pointing angle thereof, and beam pointing errors of all the wave sites are calculated, wherein the beam pointing errors comprise azimuth beam pointing errorsAnd elevation beam pointing errorθ1,θ2,...,θKPointing azimuth beam with errorAnd elevation beam pointing errorθ1,θ2,...,θKMaking table data and storing the table data in a wave control machine hardware system;
wherein the wave position 1 pointing error isWave position 2 pointing error isWave position 3 pointing error isWave position 4 pointing error is
Step S5, dividing a wave beam pointing error plane of a wave position point on the wave position network;
dividing 4 wave sites in the wave position network into a first wave beam pointing error plane A and a second wave beam pointing error plane B by taking any wave site in the wave position network as a reference point, and decomposing the slope of each wave beam pointing error plane into an azimuth direction wave beam pointing error slope and a pitch direction wave beam pointing error slope;
step S6, calculating the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope of the beam pointing error plane according to the beam pointing error of the wave point on the wave position grid, and making the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope into table data and storing the table data in the wave controller hardware system;
wherein,the wave position interval in the azimuth direction is 4 wave position pointsDelta theta is the wave position interval in the pitching direction theoretically of 4 wave positions, namely Delta theta2-θ1；Error slope, k, for azimuth beam pointingθ1,kθ2,...,kθKA beam pointing error slope in a pitch direction;
step S7, selecting a corresponding beam pointing error plane according to the angle difference between the wave position point to be calculated and the reference point to calculate the beam pointing error on line, generating a phased array antenna beam compensation table by the beam pointing errors of all the wave position points in the visible angle area, and storing the phased array antenna beam compensation table in a wave control machine hardware system to realize the compensation of the phased array antenna beam pointing error;
wherein,andthe selection of the plane is shown, the wave position point 1 is taken as an example as a reference point,indicating that a second beam pointing error plane B is selected,indicating selection of a first beam pointing error plane a;
and the formula (3) is written into a wave control machine hardware system, so that the beam pointing error value can be calculated on line conveniently.
The step S7 includes the following steps:
step S701, calculating an angle difference value between a wave position point to be calculated and a reference point;
taking wave position point 1 as an example of a reference point, the azimuth direction angle difference is calculatedAnd a difference in pitch direction angle (theta-theta)1)；
Step S702, selecting a corresponding beam pointing error plane to calculate a beam pointing error according to the size of the azimuth direction angle difference and the pitch direction angle difference;
comparing the difference of azimuth anglesAngle difference of the sum and azimuthThe size of (1) whenIn the second beam pointing error plane BCalculating the beam pointing error of a wave position point to be calculated; when in useIn the first beam pointing error plane AAnd calculating the beam pointing error of the wave position to be calculated.
High-order interception is carried out on the dispersion angle of the input wave control machine hardware system, low 8-order data are abandoned, and a wave beam control command is obtained according to the intercepted high orderi＝1,2,…,I、θkK-1, 2, …, M value, used to generate the lookup angle slopekθ1,kθ2,...,kθKSum beam pointing error value θ1,θ2,...,θKThe address of the data table, and the lower 8 of the discard is the data representing the azimuth beam pointing error valueAnd elevation beam pointing error theta-theta1,θ-θ2,...θ-θMA value for calculating a beam pointing error.
Said azimuth angle intervalAnd the pitch angle interval Δ θ are both 2.56 °.
The invention realizes the calculation of the phased array antenna beam pointing error compensation value by using the bilinear fitting method, provides technical support for improving the phased array antenna beam pointing accuracy, calculates the beam pointing error value between given wave positions by using the bilinear fitting method according to the beam pointing error of the phased array radar seeker antenna at the limited wave positions, writes the error compensation value into a beam control system, realizes the effective compensation of the phased array antenna beam pointing error, thereby improving the phased array antenna beam pointing accuracy, reduces the wave positions of the phased array antenna to be measured by using the bilinear interpolation fitting method, saves the measurement test time, greatly improves the efficiency and reduces the cost.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a schematic diagram of beam pointing angles for a given ideal wave position in an embodiment of the present invention.
Fig. 3 is a schematic diagram of an actual beam pointing angle of a phased array antenna obtained by scanning with a wave control machine according to an embodiment of the present invention.
FIG. 4 is a diagram illustrating beam pointing errors between an ideal wave position and an actual beam pointing according to an embodiment of the present invention.
FIG. 5 is a schematic diagram of the curved surface division and bilinear beam pointing error slope calculation according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of a Ku-band planar phased array antenna according to an embodiment of the present invention.
Figure 7a is a plot of beam pointing error in a azimuth direction before beam pointing error compensation in an embodiment of the present invention.
Figure 7b is a diagram of beam pointing errors before beam pointing error compensation in an embodiment of the present invention.
FIG. 7c is a diagram of the azimuth beam pointing error after beam pointing error compensation in an embodiment of the present invention.
Figure 7d is a plot of the elevation beam pointing error after beam pointing error compensation in an embodiment of the present invention.
Detailed Description
The preferred embodiment of the present invention will be described in detail with reference to fig. 1 to 7 d.
The most effective method for compensating the antenna beam pointing error is to measure the error value of each antenna wave position and then compensate the error value. Because of the large number of wave bits of the phased array antenna, the measurement of the antenna pointing error by performing experiments on each wave bit causes heavy test workload and time inadmissible. Therefore, effective fitting calculation needs to be performed according to actually measured data of antenna pointing errors of the phased array antenna on the finite wave positions by combining an effective fitting calculation method, so that beam pointing errors on non-test wave positions are obtained and compensated. Because the wave positions of the wave beams of the phased array antenna form a two-dimensional space grid, other wave positions needing to be calculated are positioned in the middle of four adjacent wave positions, and a two-dimensional fitting calculation method is needed. Because the fitting calculation methods are more, but most of the fitting calculation methods are limited in the category of theoretical research and data analysis, and a complex fitting calculation formula is not easy to implement in engineering, therefore, an easy-to-implement bilinear fitting calculation method is provided by combining with an FPGA (field programmable gate array) fast wave control hardware system, the data of any curved surface can be obtained by the provided bilinear fitting method, the calculation method is simple and convenient to implement, and an effective way is provided for the pointing error compensation of each angle in the visible area of the phased array antenna.
Given a two-dimensional spatial viewing angle ofAt certain discrete angular intervalsDiscretizing the viewing angle into wave positionsi＝1,2,L,I,k＝1,2,L,K is whereinAnd V theta are respectively an azimuth angle discrete interval and a pitch angle discrete interval, the angle step of the current wave control machine is 0.01 DEG, and the angle intervalAnd V theta should be 0.01 degrees, scanning is carried out within the range of (-60 degrees and 60 degrees) according to the orientation and the pitching requirement of the phased array antenna, the number of azimuth wave bits is 12001, the number of pitching wave bits is 12001, and the number of spatial two-dimensional angle wave bits is 144024001. If the beam pointing error measurement is performed on all the wave positions and then compensation is performed, 144024001 wave position beam pointing error measurements are required, 1s is required for 1 wave position pointing error measurement, 144024001s is required for completing beam pointing error compensation of one phased array antenna, 1667 days are required for continuous measurement of 24 hours in 1 whole day, and the time is unacceptable. Therefore, the discrete angle interval needs to be enlarged, a limited number of wave positions are selected for measuring the beam pointing error, in order to facilitate the simplification of wave control software,and Δ θ can be set to 2 of 0.01 °NIn this embodiment, the measurement time is greatly reduced to the original (1/65536) s by measuring the beam pointing error at an angle interval of 2.56 °. And calculating the beam pointing error value of the intermediate wave position by using a bilinear interpolation fitting method according to the beam pointing error at the interval of 2.56 degrees, and further performing effective compensation.
As shown in fig. 1, the present invention provides a method for compensating a beam pointing error of a phased array antenna based on bilinear fitting, which comprises the following steps:
step S1, discretizing a visual angle area;
given a scan angle region of a phased array antennaAt azimuthal angular intervalsDiscretizing the visual angle area by the pitch angle interval delta theta to obtain a space discrete anglei＝1,2,…,I,k＝1,2,…,K；
Randomly selecting 4 adjacent reference wave sites in space, wherein the wave site 1 isWave site 2 isWave site 3 isWave site 4 is
4 adjacent reference wave sites are randomly selected from all the wave sites, 4 adjacent wave sites form a wave position grid, and the wave beam pointing error compensation of the wave positions in the wave position grid formed by the 4 wave sites can be realized by the method;
step S2, obtaining an ideal beam pointing angle of a wave position point on the wave position grid;
obtaining the ideal beam pointing angles of all wave positions according to theoretical calculation, and randomly selecting 4 adjacent space wave positions, as shown in fig. 2, the ideal beam pointing direction of the wave position 1 isThe ideal beam pointing at wave position 2 isThe ideal beam pointing at wave position point 3 isThe ideal beam pointing at wave position 4 is
Step S3, obtaining the actual beam pointing angle of the wave position point on the wave position network;
combining with the actual test system, the wave beam scanning is performed by the wave control machine to obtain the actual beam pointing angles of all wave positions, as shown in fig. 3, the wave position 1 is (α)1,β1) Wave position 2 is (α)2,β1) Wave position 3 is (α)1,β2) Wave position 4 is (α)2,β2)；
Step S4, calculating the beam pointing error of the wave position on the wave position grid;
calculating beam pointing errors of all wave points according to the ideal beam pointing angles of all the wave points and the corresponding actual beam pointing angles, wherein the beam pointing errors comprise azimuth beam pointing errorsAnd elevation beam pointing errorθ1,θ2,...,θKPointing azimuth beam with errorAnd elevation beam pointing errorθ1,θ2,...,θKMaking table data and storing the table data in a wave control machine hardware system;
in equation (1), the beam pointing errors due to different wave sites are not necessarily the same, i.e. Respectively are error values of the azimuth directions of four wave position points on any complex curved surface,θ1、θ2、θ3、θ4respectively, the error values of the pitching directions of the four wave positions on the arbitrary complex curved surface are shown in fig. 4, and the pointing error of the wave position 1 isWave position 2 pointing error isWave position 3 pointing error isWave position 4 pointing error is
Step S5, dividing a wave beam pointing error plane of a wave position point on the wave position network;
as shown in fig. 5, a plane is determined according to 3 points in space, and 4 wave sites in the wave site network are divided into 2 beam pointing error planes (plane a and plane B) by using any one wave site in the wave site network as a reference point, and the slope of each plane is decomposed into an azimuth beam pointing error slope and a pitch beam pointing error slope;
the wave point 1 is used as a reference point to realize beam pointing error plane division, the wave point 3 can also be used as a reference point to carry out beam pointing error plane division, and the calculation of the beam pointing error slope in the azimuth direction and the beam pointing error slope in the elevation direction is completed in each plane;
step S6, calculating the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope of the beam pointing error plane according to the beam pointing error of the wave position on the wave position grid;
as shown in fig. 5, according to the linear relationship, the pointing error slope calculation formula corresponding to the 4 sides of the two beam pointing error planes is as follows:
wherein,the wave position interval in the azimuth direction is 4 wave position pointsDelta theta is the wave position interval in the pitching direction theoretically of 4 wave positions, namely Delta theta2-θ1；Error slope, k, for azimuth beam pointingθ1,kθ2,...,kθKA beam pointing error slope in a pitch direction;
the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope are made into table data and stored in a wave control machine hardware system;
step S7, selecting a corresponding beam pointing error plane according to the angle difference between the wave position point to be calculated and the reference point to calculate the beam pointing error on line, generating a phased array antenna beam compensation table by the beam pointing errors of all the wave position points in the visible angle area, and storing the phased array antenna beam compensation table in a wave control machine hardware system to realize the compensation of the phased array antenna beam pointing error;
wherein,andthe selection of the representation plane is shown,it is shown that the selection plane B is selected,represents the selection plane a;
the formula (3) is written into a wave control machine hardware system, so that the wave beam pointing error value can be calculated on line conveniently;
the step S7 includes the following steps:
step S701, calculating an angle difference value between a wave position point to be calculated and a reference point;
calculating the angle difference of the azimuth direction by taking the wave position point 1 as a reference pointAnd a difference in pitch direction angle (theta-theta)1)；
Step S702, selecting a corresponding beam pointing error plane to calculate a beam pointing error according to the size of the azimuth direction angle difference and the pitch direction angle difference;
comparing the difference of azimuth anglesAnd azimuth angle difference (theta-theta)1) The size of (1) whenIn plane B according toCalculating the beam pointing error of a wave position point to be calculated; when in useIn plane A according toCalculating the beam pointing error of a wave position point to be calculated;
due to the azimuthal interval in this embodimentAnd the pitch angle interval delta theta are both adopted to carry out the measurement of the beam pointing error of the phased array antenna by 2.56 degrees, while the angle stepping of the wave control machine is 0.01 degree, and the angle interval is 256 times of the angle stepping, namely 28Therefore, the high order is directly intercepted from the input angle value in the wave control machine, the low order 8 orders are directly abandoned, and the wave beam control instruction can be obtainedi＝1,2,…,I、θkK is 1,2, …, M value, used to generate the address for looking up the angle slope and error angle table, and the discarded lower 8 bits of data areθ-θ1,θ-θ2,...θ-θKIt can be used for wave control formula calculation and making table to store data θ1,θ2,...,θK、kθ1,kθ2,...,kθKThe beam pointing command of each wave bit can be obtained by using the formula (3).
According to the invention, for different phased array antennas, the azimuth and elevation beam pointing errors of the wave point points are measured through tests, and the error slope between all adjacent wave points can be obtained according to the azimuth and elevation angle intervals, so that the beam pointing error compensation of the phased array antenna can be realized on line, and the method is effective, simple and easy to realize.
In a preferred embodiment of the invention, a Ku-band planar phased array radar seeker is provided, wherein a phased array antenna array of the Ku-band planar phased array radar seeker is composed of 512 antenna units, the line spacing and the column spacing of the antenna units are both half wavelengths, and the phased array antenna units are distributed as shown in fig. 6. The beam pointing error measurement test is carried out through a phased array antenna beam pointing accuracy test system, the pointing accuracy of the phased array antenna in the range of the azimuth direction (-60 degrees, 60 degrees) and the pitching direction (-60 degrees, 60 degrees) is tested, and the beam pointing accuracy after the beam pointing error is corrected is verified at the same time. Fig. 7a and 7b are beam pointing error diagrams of the Ku-band planar phased array antenna in the azimuth direction and the elevation direction, respectively, as can be seen from fig. 7a and 7b, wherein the azimuth pointing error is positive or negative, and the elevation pointing error is positive; the pointing error is gradually increased along with the increase of the pointing angle, but the pointing jump basically does not exist; within ± 30 °, the pointing error is less than 0.1 °, and at 60 ° the pointing error is maximal, about 0.6 °. According to the data of the figures 7a and 7b, the antenna beam pointing error compensation is carried out, the figures 7c and 7d are beam pointing error graphs in the azimuth direction and the elevation direction respectively after the antenna beam pointing error compensation, as seen from the graphs, the pointing errors in the azimuth direction and the elevation direction are greatly reduced, the pointing errors are gradually increased along with the increase of the pointing angle, no pointing jump exists, the pointing error is not more than 0.02 degrees within +/-30 degrees, the pointing error is maximum at 60 degrees and is about 0.1 degrees, the pointing accuracy is greatly improved, and the verification method is very effective.
Phased array radar seeker needs to be developed in batches, phased array antennas need to carry out antenna beam pointing error compensation, if each phased array antenna is required to carry out pointing error measurement at a discrete angle of a beam stepping angle level, then error compensation is carried out, and the phased array radar seeker is a huge project and cannot complete. The method for compensating the beam pointing error of the phased array antenna by utilizing bilinear fitting can greatly improve the beam pointing accuracy of the phased array antenna, is simple and easy to realize in engineering, and provides a feasible, efficient and high-precision technical approach for compensating the beam pointing error of the phased array antenna. Through a large number of phased array antenna pointing accuracy tests, test results prove that the phased array antenna beam pointing error compensation method based on bilinear fitting provided by the invention is very effective.
While the present invention has been described in detail with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.
Claims (4)
1. A phased array antenna beam pointing error compensation method based on bilinear fitting is characterized by comprising the following steps:
step S1, discretizing a visual angle area;
given a scan angle region of a phased array antennaAt azimuthal angular intervalsDiscretizing the visual angle region by the discrete angle formed by the pitch angle interval delta theta to obtain a spatial discrete anglei＝1,2,…,I,k＝1,2,…,K；
Taking the discrete angle as an ideal wave site, and forming a wave site grid by any 4 adjacent wave sites;
step S2, obtaining an ideal beam pointing angle of a wave position point on the wave position grid;
obtaining ideal wave beam pointing angles of all wave sites according to theoretical calculation, randomly selecting 4 adjacent space wave sites, and pointing the ideal wave beam of the wave site 1 to beThe ideal beam pointing at wave position 2 isThe ideal beam pointing at wave position point 3 isThe ideal beam pointing at wave position 4 is
Step S3, obtaining the actual beam pointing angle of the wave position point on the wave position network;
combining with an actual test system, scanning antenna beams through a wave control machine to obtain actual beam pointing angles of all wave positions, wherein the actual beam pointing direction of the wave position 1 is (α)1,β1) The actual beam pointing at wave position 2 is (α)2,β1) The actual beam pointing at wave position 3 is (α)1,β2) The actual beam pointing direction 4 of the wave site is (α)2,β2)；
Step S4, calculating the beam pointing error of the wave position on the wave position grid;
calculating beam pointing errors of all wave points according to the ideal beam pointing angles of all the wave points and the corresponding actual beam pointing angles, wherein the beam pointing errors comprise azimuth beam pointing errorsAnd elevation beam pointing errorθ1,θ2,...,θKPointing azimuth beam with errorAnd elevation beam pointing errorθ1,θ2,...,θKMaking table data and storing the table data in a wave control machine hardware system;
wherein the wave position 1 pointing error isWave position 2 pointing error isWave position 3 pointing error isWave position 4 pointing error is
Step S5, dividing a wave beam pointing error plane of a wave position point on the wave position network;
dividing 4 wave sites in the wave position network into a first wave beam pointing error plane A and a second wave beam pointing error plane B by taking any wave site in the wave position network as a reference point, and decomposing the slope of each wave beam pointing error plane into an azimuth direction wave beam pointing error slope and a pitch direction wave beam pointing error slope;
step S6, calculating the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope of the beam pointing error plane according to the beam pointing error of the wave point on the wave position grid, and making the azimuth direction beam pointing error slope and the elevation direction beam pointing error slope into table data and storing the table data in the wave controller hardware system;
wherein,the wave position interval in the azimuth direction is 4 wave position pointsDelta theta is the wave position interval in the pitching direction theoretically of 4 wave positions, namely Delta theta2-θ1；Error slope, k, for azimuth beam pointingθ1,kθ2,...,kθKA beam pointing error slope in a pitch direction;
step S7, selecting a corresponding beam pointing error plane according to the angle difference between the wave position point to be calculated and the reference point to calculate the beam pointing error on line, generating a phased array antenna beam compensation table by the beam pointing errors of all the wave position points in the visible angle area, and storing the phased array antenna beam compensation table in a wave control machine hardware system to realize the compensation of the phased array antenna beam pointing error;
wherein,andthe selection of the plane is shown, the wave position point 1 is taken as an example as a reference point,indicating that a second beam pointing error plane B is selected,indicating selection of a first beam pointing error plane a;
and the formula (3) is written into a wave control machine hardware system, so that the beam pointing error value can be calculated on line conveniently.
2. The method for compensating for beam pointing error of a phased array antenna based on bilinear fitting as claimed in claim 1, wherein said step S7 comprises the steps of:
step S701, calculating an angle difference value between a wave position point to be calculated and a reference point;
taking wave position point 1 as an example of a reference point, the azimuth direction angle difference is calculatedAnd a difference in pitch direction angle (theta-theta)1)；
Step S702, selecting a corresponding beam pointing error plane to calculate a beam pointing error according to the size of the azimuth direction angle difference and the pitch direction angle difference;
comparing the difference of azimuth anglesAnd azimuth angle difference (theta-theta)1) The size of (1) whenIn the second beam pointing error plane BCalculating the beam pointing error of a wave position point to be calculated;
when in useIn the first beam pointing error plane AAnd calculating the beam pointing error of the wave position to be calculated.
3. The phased array antenna beam pointing error compensation method based on bilinear fitting of claim 2, characterized in that high-order interception is performed on the discrete angle input to the wave controller hardware system, low-order 8-order data is discarded, and a beam control command is obtained according to the intercepted high orderi＝1,2,…,I、θkK-1, 2, …, M value, used to generate the lookup angle slopekθ1,kθ2,...,kθKSum beam pointing error value θ1,θ2,...,θKThe address of the data table, and the lower 8 of the discard is the data representing the azimuth beam pointing error valueAnd elevation beam pointing error theta-theta1,θ-θ2,...θ-θMA value for calculating a beam pointing error.
4. Bilinear fitting based facies as claimed in claim 3The method for compensating the beam pointing error of the steering array antenna is characterized in that the azimuth angle intervalAnd the pitch angle interval Δ θ are both 2.56 °.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510581867.5A CN105301564A (en) | 2015-09-14 | 2015-09-14 | Phased array antenna beam pointing error compensation method based on bilinear fitting |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510581867.5A CN105301564A (en) | 2015-09-14 | 2015-09-14 | Phased array antenna beam pointing error compensation method based on bilinear fitting |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105301564A true CN105301564A (en) | 2016-02-03 |
Family
ID=55199046
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510581867.5A Pending CN105301564A (en) | 2015-09-14 | 2015-09-14 | Phased array antenna beam pointing error compensation method based on bilinear fitting |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105301564A (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107949066A (en) * | 2017-11-21 | 2018-04-20 | 西安空间无线电技术研究所 | A kind of ripple position resource flexible scheduling system and dispatching method towards beam-hopping |
| CN107991668A (en) * | 2017-11-20 | 2018-05-04 | 中国电子科技集团公司第十四研究所 | A kind of vacant lot deviational survey suitable for phased-array radar is away from method |
| CN108205137A (en) * | 2016-12-20 | 2018-06-26 | 北京行易道科技有限公司 | Lens radar and the vehicles |
| CN109239736A (en) * | 2018-07-28 | 2019-01-18 | 上海航天电子有限公司 | A kind of phase array antenna beam error in pointing modification method |
| CN109541558A (en) * | 2018-09-30 | 2019-03-29 | 航天恒星科技有限公司 | A kind of calibration method of whole process total system Active Phase-Array Radar target seeker |
| CN109639329A (en) * | 2018-11-16 | 2019-04-16 | 上海无线电设备研究所 | The only quick shaping method of phase weighting wave beam |
| CN110426690A (en) * | 2019-07-02 | 2019-11-08 | 中国航空工业集团公司雷华电子技术研究所 | A kind of airborne weather radar beam position automatic calibrating method |
| CN110531328A (en) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | A kind of spaceborne Circular scanning radar beam center error in pointing estimation method |
| CN110556630A (en) * | 2019-09-18 | 2019-12-10 | 上海无线电设备研究所 | Special-shaped radome aiming line error compensation method |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN113419208A (en) * | 2021-05-11 | 2021-09-21 | 四川九洲防控科技有限责任公司 | Interference source real-time direction finding method and device, storage medium and electronic equipment |
| CN113899956A (en) * | 2021-09-18 | 2022-01-07 | 中国电子科技集团公司第十四研究所 | Phased array antenna array element calibration and directional diagram recording method |
| CN114553334A (en) * | 2022-04-28 | 2022-05-27 | 浩泰智能(成都)科技有限公司 | Phased array antenna pointing error measurement method, system, terminal and device |
| CN116826381A (en) * | 2023-08-23 | 2023-09-29 | 成都天成电科科技有限公司 | Phased array antenna pointing accuracy correction method, device and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS597277A (en) * | 1982-07-06 | 1984-01-14 | Mitsubishi Electric Corp | Phased array radar |
| GB2435129B (en) * | 2006-02-10 | 2009-11-11 | Thales Holdings Uk Plc | Antenna signal processing apparatus |
| CN102565756A (en) * | 2010-12-17 | 2012-07-11 | 上海无线电设备研究所 | Single-channel radiolocation direction-finding system direction-finding method thereof |
| US20120206291A1 (en) * | 2011-02-11 | 2012-08-16 | Src, Inc. | Bench-top measurement method, apparatus and system for phased array radar apparatus calibration |
-
2015
- 2015-09-14 CN CN201510581867.5A patent/CN105301564A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS597277A (en) * | 1982-07-06 | 1984-01-14 | Mitsubishi Electric Corp | Phased array radar |
| GB2435129B (en) * | 2006-02-10 | 2009-11-11 | Thales Holdings Uk Plc | Antenna signal processing apparatus |
| CN102565756A (en) * | 2010-12-17 | 2012-07-11 | 上海无线电设备研究所 | Single-channel radiolocation direction-finding system direction-finding method thereof |
| US20120206291A1 (en) * | 2011-02-11 | 2012-08-16 | Src, Inc. | Bench-top measurement method, apparatus and system for phased array radar apparatus calibration |
Non-Patent Citations (1)
| Title |
|---|
| 李毓琦等: "相控阵天线波束指向误差的线性插值补偿法", 《制导与引信》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108205137A (en) * | 2016-12-20 | 2018-06-26 | 北京行易道科技有限公司 | Lens radar and the vehicles |
| CN108205137B (en) * | 2016-12-20 | 2024-03-01 | 北京行易道科技有限公司 | Lens radar and vehicle |
| CN107991668A (en) * | 2017-11-20 | 2018-05-04 | 中国电子科技集团公司第十四研究所 | A kind of vacant lot deviational survey suitable for phased-array radar is away from method |
| CN107949066A (en) * | 2017-11-21 | 2018-04-20 | 西安空间无线电技术研究所 | A kind of ripple position resource flexible scheduling system and dispatching method towards beam-hopping |
| CN109239736A (en) * | 2018-07-28 | 2019-01-18 | 上海航天电子有限公司 | A kind of phase array antenna beam error in pointing modification method |
| CN109541558A (en) * | 2018-09-30 | 2019-03-29 | 航天恒星科技有限公司 | A kind of calibration method of whole process total system Active Phase-Array Radar target seeker |
| CN109639329A (en) * | 2018-11-16 | 2019-04-16 | 上海无线电设备研究所 | The only quick shaping method of phase weighting wave beam |
| CN109639329B (en) * | 2018-11-16 | 2022-03-29 | 上海无线电设备研究所 | Phase-only weighted beam fast shaping method |
| CN110531328B (en) * | 2019-06-25 | 2021-09-07 | 西安空间无线电技术研究所 | Satellite-borne ring-scan radar beam center pointing error estimation method |
| CN110531328A (en) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | A kind of spaceborne Circular scanning radar beam center error in pointing estimation method |
| CN110426690A (en) * | 2019-07-02 | 2019-11-08 | 中国航空工业集团公司雷华电子技术研究所 | A kind of airborne weather radar beam position automatic calibrating method |
| CN110556630B (en) * | 2019-09-18 | 2021-06-04 | 上海无线电设备研究所 | Special-shaped radome aiming line error compensation method |
| CN110556630A (en) * | 2019-09-18 | 2019-12-10 | 上海无线电设备研究所 | Special-shaped radome aiming line error compensation method |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN113419208A (en) * | 2021-05-11 | 2021-09-21 | 四川九洲防控科技有限责任公司 | Interference source real-time direction finding method and device, storage medium and electronic equipment |
| CN113899956A (en) * | 2021-09-18 | 2022-01-07 | 中国电子科技集团公司第十四研究所 | Phased array antenna array element calibration and directional diagram recording method |
| CN113899956B (en) * | 2021-09-18 | 2023-10-24 | 中国电子科技集团公司第十四研究所 | Phased array antenna array element calibration and directional diagram recording method |
| CN114553334A (en) * | 2022-04-28 | 2022-05-27 | 浩泰智能(成都)科技有限公司 | Phased array antenna pointing error measurement method, system, terminal and device |
| CN114553334B (en) * | 2022-04-28 | 2022-07-22 | 浩泰智能(成都)科技有限公司 | Phased array antenna pointing error measurement method, system, terminal and device |
| CN116826381A (en) * | 2023-08-23 | 2023-09-29 | 成都天成电科科技有限公司 | Phased array antenna pointing accuracy correction method, device and system |
| CN116826381B (en) * | 2023-08-23 | 2023-11-14 | 成都天成电科科技有限公司 | Phased array antenna pointing accuracy correction method, device and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105301564A (en) | Phased array antenna beam pointing error compensation method based on bilinear fitting | |
| CN106407723B (en) | The determination method of sparse arrangement array antenna exciting current amplitude towards Sidelobe | |
| CN109669178A (en) | A kind of spaceborne three array element pulse bidimensional direction-finding method | |
| CN109633526B (en) | Direction finding ambiguity resolving method of non-uniform circular array phase interferometer based on direction function | |
| CN111624565B (en) | Multi-area joint beam tracking method for large conformal phased array system | |
| CN104063544B (en) | Near-field effect error analysis method | |
| CN104007413B (en) | Consider the array position error calibration method of information source azimuthal error | |
| CN102683899A (en) | Other vector rotation calibration method for phased array antenna | |
| CN105606906B (en) | A kind of millimeter wave phased array test calibration method | |
| CN104375129B (en) | A kind of distributive array coherent synthesis radar emission coherent parameter calibrating method | |
| CN104200110A (en) | Decoupling-based two-dimensional beam scanning angle measurement method | |
| CN104409853A (en) | Method for controlling electric scanning beam shape of planar array antenna | |
| CN109343004A (en) | The iteration for improving Planar Phased Array Antenna beam-pointing accuracy presents phase calculation method | |
| CN109818689B (en) | Array antenna calibration method, device, system and computer readable storage medium | |
| CN104852775A (en) | Rapid evaluation method for large phased-array antenna wave beam performance | |
| CN106990385B (en) | Non-homogeneous face battle array design method based on minimum-redundancy linear arrays | |
| CN111521988B (en) | Radar angle measurement method and device based on beam forming, radar and vehicle | |
| CN109818688B (en) | Array antenna calibration method, device, system and computer readable storage medium | |
| CN103955602A (en) | Array factor molding method of synthesis aperture microwave radiometer | |
| CN107991668A (en) | A kind of vacant lot deviational survey suitable for phased-array radar is away from method | |
| CN104199005B (en) | The distributed meter wave array radar optimum length of base determines method | |
| CN116908805B (en) | Combined sum and difference angle measurement method, system and storage medium of two-dimensional radar | |
| CN106092047B (en) | Beam position deviation computational methods when large-scale antenna free-position based on track unevenness | |
| CN106017417B (en) | Large-scale wheel-track type antenna based on track unevenness is directed toward the computational methods of root-mean-square error | |
| CN113504518B (en) | External field calibration method for phased array electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into 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: 20160203 |