CN114070403A - Feedforward tracking control method and system for inter-satellite laser communication system - Google Patents
Feedforward tracking control method and system for inter-satellite laser communication system Download PDFInfo
- Publication number
- CN114070403A CN114070403A CN202111346355.2A CN202111346355A CN114070403A CN 114070403 A CN114070403 A CN 114070403A CN 202111346355 A CN202111346355 A CN 202111346355A CN 114070403 A CN114070403 A CN 114070403A
- Authority
- CN
- China
- Prior art keywords
- satellite
- tracking
- coordinate system
- feedforward
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000007246 mechanism Effects 0.000 claims abstract description 65
- 230000004044 response Effects 0.000 claims abstract description 14
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 88
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000013459 approach Methods 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 6
- 230000000630 rising effect Effects 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/11—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
- H04B10/118—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum specially adapted for satellite communication
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18519—Operations control, administration or maintenance
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Astronomy & Astrophysics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Aviation & Aerospace Engineering (AREA)
- Optical Communication System (AREA)
Abstract
The invention provides a feedforward tracking control method and a feedforward tracking control system for an inter-satellite laser communication system. The method comprises the following steps: acquiring double-star orbit information, and acquiring light spot centroid position information according to real-time feedback information of a detection camera; if the light spot center of mass is in the fine tracking window, calculating the tracking quantity of the fine aiming mechanism according to the center of mass miss distance, and if the light spot center of mass is outside the fine tracking window, calculating the tracking quantity of the coarse aiming mechanism according to the center of mass miss distance; calculating a tracking error caused by the relative motion of the satellite under the condition of considering the system time delay according to the response control bandwidth of the tracking system and the relative motion speed between the satellites to obtain a feedforward tracking compensation quantity; and if the feedforward tracking compensation quantity is smaller than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the tracking quantity of the fine aiming mechanism to obtain the tracking angle control quantity. The method effectively improves the tracking control precision and improves the steady-state tracking capability of the link between the laser satellites.
Description
Technical Field
The invention relates to the field of optical communication, in particular to a feedforward tracking control method and a feedforward tracking control system for an inter-satellite laser communication system.
Background
The inter-satellite laser communication system has the advantages of high communication speed, large communication capacity, strong anti-interference capability, good confidentiality, light and small terminal and the like, and is one of important development trends of space networks in the future. At present, the laser communication test among the satellites is carried out successively in China, the United states, Europe and the like. In China, satellite laser communication tests are carried out on ocean No. 2 satellites at the earliest time in 2011. Thereafter, the satellite laser communication experiment is successfully carried out in China on the practice of the No. thirteen satellite and the Beidou system. In 2020, the aerospace science and technology group eight hospitals and the aerospace science and technology group cloud engineering successively performed low-orbit satellite laser communication tests. Because the inter-satellite laser communication distance is long, the relative movement speed between the satellites is high, and the divergence angle of signal beams is small, the accuracy and the frequency response characteristic of tracking control directly influence the stability of a link between the laser satellites, which puts high requirements on the tracking technology of an inter-satellite laser communication system.
The traditional tracking method only calculates the tracking control quantity according to the spot centroid position fed back by the camera, and has the problems of high steady-state tracking difficulty of an inter-satellite laser communication link with high-speed relative motion due to tracking response lag, difficulty in long-time stable communication of the inter-satellite laser link and the like, and has adverse effect on high-speed networking of a space laser network system.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention aims to provide a feedforward tracking control method and a feedforward tracking control system for an inter-satellite laser communication system.
In order to achieve the above object, the present invention provides a feedforward tracking control method for an inter-satellite laser communication system, comprising the following steps:
acquiring double-star orbit information, and acquiring light spot centroid position information according to real-time feedback information of a detection camera;
judging whether the light spot mass center of each satellite in the double satellites is in the respective precise tracking window, and if the light spot mass center of any satellite is in the precise tracking window of the satellite, calculating the precise sighting mechanism tracking quantity of the satellite according to the mass center miss quantity; if the light spot center of mass of any satellite is outside the fine tracking window of the satellite, calculating the tracking quantity of the coarse aiming mechanism of the satellite according to the center of mass miss quantity;
calculating a tracking error caused by relative motion of the satellites under the condition of considering system time delay according to the response control bandwidth of the tracking system and the relative motion speed between the satellites to obtain respective feedforward tracking compensation quantities of the double satellites;
judging the magnitude relation between the respective feedforward tracking compensation quantity of the double satellites and the control precision of the coarse aiming mechanism, if the feedforward tracking compensation quantity is greater than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity of the satellite and the corresponding tracking quantity of the coarse aiming mechanism to obtain the tracking angle control quantity of the satellite, and controlling the coarse aiming mechanism to perform closed-loop tracking control on the satellite according to the tracking angle control quantity; and if the feedforward tracking compensation quantity is smaller than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the tracking quantity of the fine aiming mechanism corresponding to the feedforward compensation quantity to obtain a tracking angle control quantity, and controlling the fine aiming mechanism to perform closed-loop tracking control on the satellite according to the tracking angle control quantity.
The method effectively improves the tracking control precision and improves the steady-state tracking capability of the link between the laser satellites.
The optimal scheme of the feedforward tracking control method of the inter-satellite laser communication system is that the calculation formula of the tracking quantity of the coarse aiming mechanism is
The tracking quantity calculation formula of the fine aiming mechanism is as follows
Wherein, thetaTracking_Az，θTracking_ElRespectively tracking the azimuth and the pitch angle; x and y are respectively the miss distance of the spot centroid in the azimuth axis and the pitch axis of the camera; alpha is alphacoarse_x，αcoarse_yThe multiplication coefficient of the miss distance and the tracking distance of the coarse aiming mechanism is obtained; alpha is alpharefined_x，αrefined_yThe multiplication coefficient of the miss distance and the tracking distance of the fine aiming mechanism is obtained; beta is acoarse_x，βcoarse_yAdding coefficients of the miss distance and the tracking distance of the coarse aiming mechanism; beta is arefined_x，βrefined_yFor miss distance and fine aimingAnd adding the coefficient by the mechanism tracking amount.
According to the optimal scheme of the feedforward tracking control method of the inter-satellite laser communication system, the calculation steps of the feedforward tracking compensation amount are as follows:
the method comprises the steps of extrapolating data in real time based on double-star orbit information data, calculating double-star position vectors of a near focus coordinate system at the time delay caused by tracking delay response, converting the double-star position vectors from the near focus coordinate system to a geocentric equatorial coordinate system, subtracting the position vectors under the geocentric equatorial coordinate system to obtain feedforward tracking vectors, converting the feedforward tracking vectors from the geocentric equatorial coordinate system to a satellite orbit coordinate system, and calculating to obtain respective feedforward tracking azimuth/pitching angle compensation quantities of double stars.
According to the optimal scheme of the feedforward tracking control method of the inter-satellite laser communication system, the double-satellite position vector solving step of the near focus coordinate system at the time delay caused by tracking lag response comprises the following steps:
mean angle of approach of two stars at time t:
MA(t)＝EA(t)-eAsin(EA(t))，EA(t) is the angle of approach of satellite A at time t, eAIs the orbital eccentricity of satellite a; mB(t)＝EB(t)-eBsin(EB(t))，EB(t) is the angle of approach of satellite B at time t, eBIs the orbital eccentricity of satellite B;
and (3) calculating the position vectors of the double stars at the current moment t under the near focus coordinate system according to the formula:
position vector r of satellite A at current moment t in near-focus coordinate systemP A(t)＝(rA(t)·cos(fA(t)))Xp A+(rA(t)·sin(fA(t)))Yp A；
Position vector r of satellite B under near focus coordinate system at current moment tP B(t)＝(rB(t)·cos(fB(t)))Xp B+(rB(t)·sin(fB(t)))Yp B；
Wherein r isA(t) is the modulus of the satellite A position vector at time t, rB(t) is the modulus of the satellite B position vector at time t, and the expression is as follows:
rA(t)＝aA(1-eAcosEA(t))，rB(t)＝aB(1-eBcosEB(t))；
wherein f isA(t) true anomaly of satellite A at time t, fB(t) is the true anomaly of satellite B at time t, and is expressed as follows:
the position vector of the satellite in the near-focus coordinate system at time t is simplified as follows:
XA p、YA pare unit vectors, X, of the satellite A in a near focus coordinate systemB p、YB pAre unit vectors, X, of the satellite B in a near focus coordinate systemA p＝XB p＝[1 0 0]T，YA p＝YB p＝[0 1 0]T。
Position vectors of the two stars under the geocentric equatorial coordinate system at the moment t:
position vector of satellite A under geocentric equatorial coordinate system at time t
Position vector of satellite B in geocentric equatorial coordinate system at time t
Wherein, the transformation matrix of the near focus coordinate system and the geocentric equatorial coordinate system of the satellite A
Transformation matrix of near-focus coordinate system and geocentric equatorial coordinate system of satellite B
ΩAIs the rising intersection right ascension and omega of the satellite AAIs the orbital perigee angle, i, of satellite AAIs the orbital inclination of the satellite A, omegaBIs the rising intersection right ascension and omega of the satellite BBOrbital perigee angle, i, for satellite BBOrbital inclination of satellite B, then Xp A、Yp ATransformation to the geocentric equatorial coordinate system is denoted Xp lA、Yp lAThe expression is as follows:
Xp B、Yp Btransformation to the geocentric equatorial coordinate system is denoted Xp lB、Yp lBThe expression is as follows:
feedforward tracking vector v of satellite A to satellite B under geocentric equatorial coordinate system at time tI AThe calculation expression of (a) is:
respectively are position vectors of the satellite A and the satellite B under the geocentric equatorial coordinate system at the moment t;
feedforward tracking vector v of satellite B to satellite A under geocentric equatorial coordinate system at time tI BThe calculation expression of (a) is:
the position vectors of the satellite A and the satellite B in the equatorial coordinate system of the geocentric at the moment t are respectively.
Transforming the feedforward tracking vector of the satellite A from the earth center equatorial coordinate system to the satellite orbit coordinate system to obtain the feedforward tracking vector of the satellite A in the orbit coordinate system
Wherein,
，fA(t) is the true anomaly of satellite A at time t;
according to the feedforward tracking vector v of the satellite A under the orbit coordinate systemO AAnd calculating the feedforward tracking pitch angle compensation quantity of the satellite A to the satellite B under the satellite orbit coordinate system at the time t
Feedforward tracking azimuth angle compensation quantity
transforming the feedforward tracking vector of the satellite B from the earth center equatorial coordinate system to the satellite orbit coordinate system to obtain the feedforward tracking vector of the satellite B in the orbit coordinate system
Wherein,
，fB(t) is the true anomaly of satellite B at time t;
according to the feedforward tracking vector v of the satellite B in the orbit coordinate systemO BAnd calculating the feedforward tracking pitch angle compensation quantity of the satellite B to the satellite A under the satellite orbit coordinate system at the time t
Feedforward tracking azimuth angle compensation quantity
The invention further provides a feedforward tracking control system of the inter-satellite laser communication system, which comprises a control unit and a storage unit, wherein the control unit is in communication connection with the storage unit, the storage unit is used for storing at least one executable instruction, and the executable instruction enables the processing unit to execute the operation corresponding to the high-precision aiming pointing method of the inter-satellite laser communication system.
The invention has the beneficial effects that: compared with the existing inter-satellite tracking technology, the method can be applied to low-orbit, medium-orbit and high-orbit satellite laser communication links, integrates factors such as satellite relative motion speed, tracking control time delay, facula mass center miss distance and the like, sets feed-forward tracking compensation quantity, effectively improves tracking control precision, and further improves steady-state tracking capability of the laser inter-satellite links.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow chart of a feedforward tracking control method of an inter-satellite laser communication system;
FIG. 2 is a schematic view of a camera tracking window;
FIG. 3 is a schematic diagram of coordinate system conversion;
fig. 4 is a diagram showing simulation results of the inter-satellite feedforward tracking control amount.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
In the description of the present invention, unless otherwise specified and limited, it is to be noted that the terms "mounted," "connected," and "connected" are to be interpreted broadly, and may be, for example, a mechanical connection or an electrical connection, a communication between two elements, a direct connection, or an indirect connection via an intermediate medium, and specific meanings of the terms may be understood by those skilled in the art according to specific situations.
As shown in fig. 1, the present invention provides a feedforward tracking control method for an inter-satellite laser communication system, and as shown in fig. 1, the feedforward tracking control method for the inter-satellite laser communication system specifically includes the following steps:
firstly, acquiring double-satellite orbit information according to a satellite-borne GPS or GNSS, and acquiring light spot centroid position information according to real-time feedback information of a detection camera.
In this embodiment, the two satellites are respectively referred to as a satellite a and a satellite B, specifically, the two-satellite orbit information acquired by the satellite-borne GPS and the GNSS includes six double-satellite orbits, including 12 parameters in total, of the two-satellite orbit semimajor axis a, the eccentricity e, the near location angle ω, the ascent point right ascension Ω, the inclination angle i, and the time T of passing through the near location, and the two-satellite orbit information is subjected to high-frequency real-time extrapolation. The detection camera is preferably an InGaAs camera, the coarse aiming mechanism is preferably an aerospace high-precision torque motor, and the fine aiming mechanism is preferably a piezoelectric deflection mirror (PZT).
And then judging whether the light spot mass center of each satellite in the double satellites is in a fine tracking window, wherein the fine tracking window is preset in a system, as shown in fig. 2, the upper surface of the fig. 2 is provided with a window, the fine tracking window and a coarse tracking window are divided, light can be detected by the terminal as long as the light enters the window, and the window in which the light is located can be judged, if the light spot mass center of any satellite is in the fine tracking window of the satellite, the tracking quantity of a fine aiming mechanism of the satellite is calculated according to the mass center miss distance, and if the light spot mass center of any satellite is outside the fine tracking window of the satellite, the tracking quantity of the coarse aiming mechanism of the satellite is calculated according to the mass center miss distance. When the center of mass of the light spot is detected, the center of mass miss distance can be obtained.
The calculation formula of the tracking quantity of the coarse aiming mechanism is as follows
The tracking quantity calculation formula of the fine aiming mechanism is as follows
Wherein, thetaTracking_Az，θTracking_ElRespectively tracking an azimuth angle and a pitch angle; x and y are respectively the miss distance of the spot centroid in the azimuth axis and the miss distance of the pitch axis in the camera；αcoarse_x，αcoarse_yThe multiplication coefficient of the miss distance and the tracking distance of the coarse aiming mechanism is obtained; alpha is alpharefined_x，αrefined_yThe multiplication coefficient of the miss distance and the tracking distance of the fine aiming mechanism is obtained; beta is acoarse_x，βcoarse_yAdding coefficients of the miss distance and the tracking distance of the coarse aiming mechanism; beta is arefined_x，βrefined_yAnd adding coefficients of the miss distance and the tracking distance of the fine aiming mechanism.
If the tracking quantity of the fine and coarse aiming mechanisms of the satellite A in the double satellites is calculated, substituting the tracking quantity of the spot centroid of the satellite A in the azimuth axis and the tracking quantity of the pitch axis into the camera to calculate; if the tracking quantity of the fine and coarse aiming mechanisms of the satellite B in the double-star is calculated, the tracking quantity is substituted into the miss distance of the spot centroid of the satellite B in the azimuth axis and the miss distance of the pitch axis in the camera for calculation.
And calculating the tracking error caused by the relative motion of the satellite under the condition of considering the system time delay according to the response control bandwidth of the tracking system and the relative motion speed between the satellites to obtain respective feedforward tracking compensation quantities of the double satellites. The response control bandwidth and the relative motion speed between the satellites can be directly obtained from the tracking system.
The step of calculating the feedforward tracking compensation amount mainly comprises the following steps: the method comprises the steps of extrapolating data in real time based on double-star orbit information data, calculating double-star position vectors of a near focus coordinate system at the time delay caused by tracking delay response, converting the double-star position vectors from the near focus coordinate system to a geocentric equatorial coordinate system, subtracting the position vectors under the geocentric equatorial coordinate system to obtain feedforward tracking vectors, converting the feedforward tracking vectors from the geocentric equatorial coordinate system to a satellite orbit coordinate system, and calculating to obtain respective feedforward tracking azimuth/pitching angle compensation quantities of double stars. The coordinate relationships are shown in fig. 3.
Specifically, the time delay time t is the position vector r of two satellites under the near-focus coordinate systemp A，rp BThe solving steps are as follows:
the average near point angle M of two satellite time t can be iteratively solved according to the Keplerian third lawA，MBThe expression is as follows:
where μ is 398601.2 as an attractive force constant, aAIs the orbital semi-major axis, T, of satellite AAIs the time of the past location of satellite A, aBIs the orbital semi-major axis, T, of the satellite BBIs the time of the past location of satellite B.
According to the Kepler equation, the approximate point angle E of the two satellites is calculated by a Newton iteration methodA，EBThe expression is as follows:
MA(t)＝EA(t)-eAsin(EA(t))，MB(t)＝EB(t)-eBsin(EB(t))，eAis the orbital eccentricity of the satellite A, eBIs the orbital eccentricity of satellite B.
Calculating the position vectors r of the two satellites at the current time t under the near-focus coordinate system according to the formulap A，rp B。
rP A(t)＝(rA(t)·cos(fA(t)))Xp A+(rA(t)·sin(fA(t)))Yp A，
rP B(t)＝(rB(t)·cos(fB(t)))Xp B+(rB(t)·sin(fB(t)))Yp B；
Wherein r isA，rBThe expression, modulo the two satellite position vectors, is as follows:
rA(t)＝aA(1-eAcosEA(t))，rB(t)＝aB(1-eBcosEB(t))；
wherein f isA，fBThe true near point angle for two satellites can be obtained from the approximate near point angle, and the expression is as follows:
and finally, the position vectors r of the two satellites at the time t in the near-focus coordinate systemp A，rp BThe simplification can be as follows:
position vector of satellite A in double stars under near focus coordinate system at moment t
Position vector of satellite B
Wherein, Xp A，Yp AIs the unit vector, X, of satellite A in its near focus coordinate systemp B，Yp BIs the unit vector, X, of satellite B in its near-focus coordinate systemp A＝[1 0 0]T，Yp A＝[0 1 0]T，Xp B＝[1 0 0]T，Yp B＝[0 1 0]T。
The expression of the position vector of the satellite a in the two-star at time t under the geocentric equatorial coordinate system is as follows:
the expression of the position vector of the satellite B in the isocenter equatorial coordinate system at time t in the two stars is as follows:
the expression of the transformation matrix of the near-focus coordinate system and the geocentric equatorial coordinate system of the satellite A is as follows:
The expression of the transformation matrix of the near-focus coordinate system and the geocentric equatorial coordinate system of the satellite B is as follows:
Xp B、Yp Btransformation to the geocentric equatorial coordinate system is denoted Xp lB、Yp lBThe expression is as follows:
feedforward tracking vector of satellite A to satellite B under geocentric equatorial coordinate system at time t
Is calculated as
The position vectors of the satellite A and the satellite B in the equatorial coordinate system of the earth center at the moment t respectively, and similarly, the feedforward tracking vector of the satellite B to the satellite A
To avoid repeating the description, the feedforward tracking vector of satellite A to satellite B is used
The process of the present invention is described in detail for the purpose of example.
The feedforward tracking vector v of the satellite A under the orbit coordinate system can be obtained by transforming the feedforward tracking vector of the satellite A into the orbit coordinate system of the satellite from the earth center equatorial coordinate systemO AThe calculation expression is
Wherein,
according to the feedforward tracking vector v of the satellite A under the orbit coordinate systemO AAnd calculating the feedforward tracking pitch angle compensation quantity of the satellite A to the satellite B under the satellite orbit coordinate system at the time t
And feedforward tracking azimuth angle compensation quantity
Respectively mean VoA(t) the first, second and third terms of the vector, VoAThe (t) vector is a matrix of 3 x 1 size.
Judging the magnitude relation between feedforward tracking compensation quantity (feedforward tracking pitch angle compensation quantity and feedforward tracking azimuth angle compensation quantity) of the satellite A to the satellite B at the time delay (t time) caused by considering tracking lag response and the control precision of the coarse aiming mechanism, if the feedforward tracking compensation quantity is greater than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the calculated tracking quantity of the coarse aiming mechanism to obtain the tracking angle control quantity (feedforward tracking control quantity of a terminal azimuth axis and a pitch axis) of the satellite A at the time delay t, and controlling the coarse aiming mechanism to perform closed-loop tracking control on the satellite A according to the tracking angle control quantity; and if the feedforward tracking compensation quantity is smaller than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the tracking quantity of the fine aiming mechanism to obtain the tracking angle control quantity of the satellite A at the time delay t, and controlling the fine aiming mechanism to perform closed-loop tracking control on the satellite A according to the tracking angle control quantity.
When the relation between the feedforward tracking compensation amount (feedforward tracking pitch angle compensation amount, feedforward tracking azimuth angle compensation amount) and the control accuracy of the coarse aiming mechanism is judged, if one of the feedforward tracking pitch angle compensation amount and the feedforward tracking azimuth angle compensation amount is larger than the control accuracy of the coarse aiming mechanism, that is, if the feedforward tracking compensation amount is larger than the control accuracy of the coarse aiming mechanism, and both the feedforward tracking pitch angle compensation amount and the feedforward tracking azimuth angle compensation amount are smaller than the control accuracy of the coarse aiming mechanism, the feedforward tracking compensation amount is smaller than the control accuracy of the coarse aiming mechanism.
Similarly, a feedforward tracking vector v of the satellite B at the time t in the orbit coordinate system can be obtainedO BThe feedforward tracking pitch angle compensation amount theta of the satellite B to the satellite AEl BFeedforward tracking azimuth angle compensation amount thetaAz B。
Judging the magnitude relation between feedforward tracking compensation quantity (feedforward tracking pitch angle compensation quantity and feedforward tracking azimuth angle compensation quantity) of the satellite B to the satellite A at a time delay (time t) caused by considering tracking lag response and the control precision of the coarse aiming mechanism, if the feedforward tracking compensation quantity is greater than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the calculated tracking quantity of the coarse aiming mechanism to obtain the tracking angle control quantity (feedforward tracking control quantity of a terminal azimuth axis and a pitch axis) of the satellite B at the time delay t, and controlling the coarse aiming mechanism to perform closed-loop tracking control on the satellite B according to the tracking angle control quantity; and if the feedforward tracking compensation quantity is smaller than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the tracking quantity of the fine aiming mechanism to obtain the tracking angle control quantity of the satellite B at the time delay t, and controlling the fine aiming mechanism to perform closed-loop tracking control on the satellite B according to the tracking angle control quantity.
The feedforward tracking control quantity of the terminal azimuth axis and the terminal pitch axis of the time delay time t satellite A is
The feedforward tracking control quantity of the terminal azimuth axis and the terminal pitch axis of the time delay time t satellite B is
According to the method, feedforward tracking compensation quantity simulation is carried out on the different-orbit satellite, wherein the height of the orbit of the satellite is about 1100km, the inclination angle of the orbit is 86.5 degrees, the phase difference of the satellite is 7.5 degrees, and the right ascension of the orbit of the satellite is 0 degree and 15 degrees respectively. As shown in fig. 4, the amount of feed forward tracking compensation for the terminal azimuth axis and the pitch axis.
The invention also provides an embodiment of the feedforward tracking control system of the inter-satellite laser communication system, which comprises a control unit and a storage unit, wherein the control unit is in communication connection with the storage unit, and the storage unit is used for storing at least one executable instruction, and the executable instruction enables the processing unit to execute the operation corresponding to the high-precision aiming pointing method of the inter-satellite laser communication system.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the invention have been shown and described, it will be understood by those of ordinary skill in the art that: various changes, modifications, substitutions and alterations can be made to the embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.
Claims (8)
1. A feedforward tracking control method for an inter-satellite laser communication system is characterized by comprising the following steps:
acquiring double-star orbit information, and acquiring light spot centroid position information according to real-time feedback information of a detection camera;
judging whether the light spot mass center of each satellite in the double satellites is in the respective precise tracking window, and if the light spot mass center of any satellite is in the precise tracking window of the satellite, calculating the precise sighting mechanism tracking quantity of the satellite according to the mass center miss quantity; if the light spot center of mass of any satellite is outside the fine tracking window of the satellite, calculating the tracking quantity of the coarse aiming mechanism of the satellite according to the center of mass miss quantity;
calculating a tracking error caused by relative motion of the satellites under the condition of considering system time delay according to the response control bandwidth of the tracking system and the relative motion speed between the satellites to obtain respective feedforward tracking compensation quantities of the double satellites;
judging the magnitude relation between the respective feedforward tracking compensation quantity of the double satellites and the control precision of the coarse aiming mechanism, if the feedforward tracking compensation quantity is greater than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity of the satellite and the corresponding tracking quantity of the coarse aiming mechanism to obtain the tracking angle control quantity of the satellite, and controlling the coarse aiming mechanism to perform closed-loop tracking control on the satellite according to the tracking angle control quantity; and if the feedforward tracking compensation quantity is smaller than the control precision of the coarse aiming mechanism, adding the feedforward compensation quantity and the tracking quantity of the fine aiming mechanism corresponding to the feedforward compensation quantity to obtain a tracking angle control quantity, and controlling the fine aiming mechanism to perform closed-loop tracking control on the satellite according to the tracking angle control quantity.
2. The feedforward tracking control method of the inter-satellite laser communication system according to claim 1, wherein the calculation formula of the tracking amount of the coarse aiming mechanism is
The tracking quantity calculation formula of the fine aiming mechanism is as follows
Wherein, thetaTracking_Az，θTracking_ElRespectively tracking the azimuth and the pitch angle; x and y are respectively the miss distance of the spot centroid in the azimuth axis and the pitch axis of the camera; alpha is alphacoarse_x，αcoarse_yThe multiplication coefficient of the miss distance and the tracking distance of the coarse aiming mechanism is obtained; alpha is alpharefined_x，αrefined_yThe multiplication coefficient of the miss distance and the tracking distance of the fine aiming mechanism is obtained; beta is acoarse_x，βcoarse_yAdding coefficients of the miss distance and the tracking distance of the coarse aiming mechanism; beta is arefined_x，βrefined_yAnd adding coefficients of the miss distance and the tracking distance of the fine aiming mechanism.
3. The feedforward tracking control method of the inter-satellite laser communication system according to claim 1, wherein the feedforward tracking compensation amount calculating step is:
the method comprises the steps of extrapolating data in real time based on double-star orbit information data, calculating double-star position vectors of a near focus coordinate system at the time delay caused by tracking delay response, converting the double-star position vectors from the near focus coordinate system to a geocentric equatorial coordinate system, subtracting the position vectors under the geocentric equatorial coordinate system to obtain feedforward tracking vectors, converting the feedforward tracking vectors from the geocentric equatorial coordinate system to a satellite orbit coordinate system, and calculating to obtain respective feedforward tracking azimuth/pitching angle compensation quantities of double stars.
4. A feedforward tracking control method of an inter-satellite laser communication system according to claim 3, wherein the two-star position vector solving step of the near-focus coordinate system at the time delay caused by the tracking lag response is:
mean angle of approach of two stars at time t:
MA(t)＝EA(t)-eAsin(EA(t))，EA(t) is the angle of approach of satellite A at time t, eAIs the orbital eccentricity of satellite a; mB(t)＝EB(t)-eBsin(EB(t))，EB(t) is the angle of approach of satellite B at time t, eBIs the orbital eccentricity of satellite B;
and (3) calculating the position vectors of the double stars at the current moment t under the near focus coordinate system according to the formula:
position vector r of satellite A at current moment t in near-focus coordinate systemP A(t)＝(rA(t)·cos(fA(t)))Xp A+(rA(t)·sin(fA(t)))Yp A；
Position vector r of satellite B under near focus coordinate system at current moment tP B(t)＝(rB(t)·cos(fB(t)))Xp B+(rB(t)·sin(fB(t)))Yp B；
Wherein r isA(t) is the modulus of the satellite A position vector at time t, rB(t) is the modulus of the satellite B position vector at time t, and the expression is as follows:
rA(t)＝aA(1-eAcosEA(t))，rB(t)＝aB(1-eBcosEB(t))；
wherein f isA(t) true anomaly of satellite A at time t, fB(t) is the true anomaly of satellite B at time t, and is expressed as follows:
the position vector of the satellite in the near-focus coordinate system at time t is simplified as follows:
XA p、YA pare unit vectors, X, of the satellite A in a near focus coordinate systemB p、YB pAre unit vectors, X, of the satellite B in a near focus coordinate systemA p＝XB p＝[1 0 0]T，YA p＝YB p＝[0 1 0]T。
5. A feed-forward tracking control method for an inter-satellite laser communication system as claimed in claim 3, wherein the position vectors of the double stars in the geocentric equatorial coordinate system at time t are:
position vector of satellite A under geocentric equatorial coordinate system at time t
Position vector of satellite B in geocentric equatorial coordinate system at time t
Wherein, the transformation matrix of the near focus coordinate system and the geocentric equatorial coordinate system of the satellite A
Transformation matrix of near-focus coordinate system and geocentric equatorial coordinate system of satellite B
ΩAIs the rising intersection right ascension and omega of the satellite AAIs the orbital perigee angle, i, of satellite AAIs the orbital inclination of the satellite A, omegaBIs the rising intersection right ascension and omega of the satellite BBOrbital perigee angle, i, for satellite BBOrbital inclination of satellite B, then Xp A、Yp ATransformation to the geocentric equatorial coordinate system is denoted Xp lA、Yp lAThe expression is as follows:
Xp B、Yp Btransformation to the geocentric equatorial coordinate system is denoted Xp lB、Yp lBThe expression is as follows:
6. the feedforward tracking control method for the inter-satellite laser communication system according to claim 3, wherein the feedforward tracking vector v of the satellite A to the satellite B under the geocentric equatorial coordinate system at the time t is the vectorI AThe calculation expression of (a) is:
respectively are position vectors of the satellite A and the satellite B under the geocentric equatorial coordinate system at the moment t;
red earth in the heart at time tFeedforward tracking vector v of satellite B to satellite A under orbit coordinate systemI BThe calculation expression of (a) is:
the position vectors of the satellite A and the satellite B in the equatorial coordinate system of the geocentric at the moment t are respectively.
7. A feedforward tracking control method of an inter-satellite laser communication system according to claim 3, wherein the feedforward tracking vector of the satellite A is transformed from the earth center equatorial coordinate system to the satellite orbit coordinate system to obtain the feedforward tracking vector v of the satellite A in the orbit coordinate systemO A，
Wherein,
，fA(t) is the true anomaly of satellite A at time t;
according to the feedforward tracking vector v of the satellite A under the orbit coordinate systemO AAnd calculating the feedforward tracking pitch angle compensation quantity of the satellite A to the satellite B under the satellite orbit coordinate system at the time t
Feedforward tracking azimuth angle compensation quantity
transforming the feedforward tracking vector of the satellite B from the earth center equatorial coordinate system to the satellite orbit coordinate system to obtain the feedforward tracking vector of the satellite B in the orbit coordinate system
Wherein,
，fB(t) is the true anomaly of satellite B at time t;
according to the feedforward tracking vector v of the satellite B in the orbit coordinate systemO BAnd calculating the feedforward tracking pitch angle compensation quantity of the satellite B to the satellite A under the satellite orbit coordinate system at the time t
Feedforward tracking azimuth angle compensation quantity
8. A feed-forward tracking control system of an inter-satellite laser communication system, which is characterized by comprising a control unit and a storage unit, wherein the control unit is connected with the storage unit in a communication way, and the storage unit is used for storing at least one executable instruction which causes the processing unit to execute the operation corresponding to the high-precision aiming pointing method of the inter-satellite laser communication system according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111346355.2A CN114070403B (en) | 2021-11-15 | 2021-11-15 | Feedforward tracking control method and system for inter-satellite laser communication system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111346355.2A CN114070403B (en) | 2021-11-15 | 2021-11-15 | Feedforward tracking control method and system for inter-satellite laser communication system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114070403A true CN114070403A (en) | 2022-02-18 |
| CN114070403B CN114070403B (en) | 2022-12-02 |
Family
ID=80271759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111346355.2A Active CN114070403B (en) | 2021-11-15 | 2021-11-15 | Feedforward tracking control method and system for inter-satellite laser communication system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114070403B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115267854A (en) * | 2022-09-22 | 2022-11-01 | 鹏城实验室 | Advanced alignment method and device based on satellite trajectory prediction |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002314487A (en) * | 2001-04-16 | 2002-10-25 | Telecommunication Advancement Organization Of Japan | Free-space optical communications unit |
| CN107707297A (en) * | 2017-11-03 | 2018-02-16 | 潘运滨 | A kind of airborne laser communication system and its communication means |
| CN113568133A (en) * | 2021-07-20 | 2021-10-29 | 西安工业大学 | Advanced aiming and precise tracking two-in-one system and light beam pointing method |
| CN113641182A (en) * | 2021-07-27 | 2021-11-12 | 东方红卫星移动通信有限公司 | High-precision aiming and pointing method and system for inter-satellite laser communication system |
-
2021
- 2021-11-15 CN CN202111346355.2A patent/CN114070403B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002314487A (en) * | 2001-04-16 | 2002-10-25 | Telecommunication Advancement Organization Of Japan | Free-space optical communications unit |
| CN107707297A (en) * | 2017-11-03 | 2018-02-16 | 潘运滨 | A kind of airborne laser communication system and its communication means |
| CN113568133A (en) * | 2021-07-20 | 2021-10-29 | 西安工业大学 | Advanced aiming and precise tracking two-in-one system and light beam pointing method |
| CN113641182A (en) * | 2021-07-27 | 2021-11-12 | 东方红卫星移动通信有限公司 | High-precision aiming and pointing method and system for inter-satellite laser communication system |
Non-Patent Citations (1)
| Title |
|---|
| 李少辉 等: "高精度卫星激光通信地面验证***", 《光学精密工程》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115267854A (en) * | 2022-09-22 | 2022-11-01 | 鹏城实验室 | Advanced alignment method and device based on satellite trajectory prediction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114070403B (en) | 2022-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kobayashi et al. | The high-resolution coronal imager (Hi-C) | |
| CN111897357B (en) | Attitude tracking control method for satellite earth scanning | |
| US6184825B1 (en) | Method and apparatus for radio frequency beam pointing | |
| CN113641182B (en) | High-precision aiming and pointing method and system for inter-satellite laser communication system | |
| Baister et al. | Pointing, acquisition and tracking for optical space communications | |
| CN113296128A (en) | System and method for establishing laser communication link between high-capture-rate low-orbit satellites | |
| Jono et al. | Acquisition, tracking, and pointing systems of OICETS for free space laser communications | |
| CN114070403B (en) | Feedforward tracking control method and system for inter-satellite laser communication system | |
| CN112769466B (en) | Low-orbit satellite constellation configuration keeping method | |
| CN114889849A (en) | Estimation method for multi-constraint flying-over optimal speed of spacecraft | |
| CN115562378B (en) | Photoelectric stable platform, angular velocity compensation method and storage medium | |
| Marshall et al. | Investigation of equatorial medium earth orbits for space solar power | |
| CN113489540B (en) | Light beam tracking method based on satellite platform fixed frequency decoupling | |
| Helvajian et al. | Mission Architecture to Reach and Operate at the Focal Region of the Solar Gravitational Lens | |
| CN112093081B (en) | Solar wing microwave remote sensing integrated load and control method thereof | |
| JP2000001200A (en) | Improved method and system for determining posture of space craft | |
| CN116184536A (en) | Method for detecting and inhibiting cold air view field interference of moon of satellite-borne microwave radiometer | |
| CN113237439B (en) | Decoupling tracking method of periscopic laser communication terminal | |
| CN109213190A (en) | Electronic sail Spacecraft Formation Flying cooperative control method on a kind of day heart suspension railway | |
| CN111891402B (en) | Mars detection ground antenna pointing recovery method based on autonomous maneuvering | |
| Yen | Nuclear electric propulsion for future NASA space science missions | |
| CN117811661A (en) | High-precision on-orbit pointing method of inter-satellite laser communication system | |
| Liu et al. | Full-field-of-view image quality analysis based on residual distribution of drift angle in agile satellite | |
| Morine | Zenith Star: A Structural Control Challenge | |
| Zhang et al. | Adaptive staring imaging control of micro-satellite based on target information |
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 | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 401120 39-1, Building 1, No. 64, Middle Huangshan Avenue, Yubei District, Chongqing Patentee after: China Star Network Application Co.,Ltd. Address before: 618 Liangjiang Avenue, Longxing Town, Yubei District, Chongqing Patentee before: Dongfanghong Satellite Mobile Communication Co.,Ltd. |