CA2000214C - Autonomous orbit control method and system for a geostationary satellite - Google Patents
Autonomous orbit control method and system for a geostationary satelliteInfo
- Publication number
- CA2000214C CA2000214C CA002000214A CA2000214A CA2000214C CA 2000214 C CA2000214 C CA 2000214C CA 002000214 A CA002000214 A CA 002000214A CA 2000214 A CA2000214 A CA 2000214A CA 2000214 C CA2000214 C CA 2000214C
- Authority
- CA
- Canada
- Prior art keywords
- satellite
- vector
- delay
- bias
- alpha
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 12
- 239000013598 vector Substances 0.000 claims abstract description 59
- 238000005259 measurement Methods 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000035945 sensitivity Effects 0.000 claims abstract description 8
- 235000019892 Stellar Nutrition 0.000 claims description 15
- 238000004891 communication Methods 0.000 claims description 4
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 4
- 230000005484 gravity Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/26—Guiding or controlling apparatus, e.g. for attitude control using jets
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/361—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using star sensors
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/363—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using sun sensors
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/365—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using horizon or Earth sensors
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
TEXT OF THE ABSTRACT
A stationkeeping method for a satellite in geostationary orbit comprises the steps of :
- determining at the same time the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, - deducing therefrom a state vector E consis-ting in orbital parameters by the formula:
Z = H.E + C.B
where:
. Z is a measurment vector the components of which are deduced from the angles .alpha.1 and .alpha.2, . H is a measuring matrix, . C is a bias sensitivity matrix, . B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, - determining stationkeeping manoeuvres and applying same by means of thrusters.
A stationkeeping method for a satellite in geostationary orbit comprises the steps of :
- determining at the same time the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, - deducing therefrom a state vector E consis-ting in orbital parameters by the formula:
Z = H.E + C.B
where:
. Z is a measurment vector the components of which are deduced from the angles .alpha.1 and .alpha.2, . H is a measuring matrix, . C is a bias sensitivity matrix, . B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, - determining stationkeeping manoeuvres and applying same by means of thrusters.
Description
- Z1J;~ ?2~ 4 The invention concerns controlling the orbit of space vehicles such as satellites and is more particularly directed to space vehicles in geostationary orbit (in which case the expression "stationkeeping" is S usually used as a synonym for "geostationary orbit control").
The position and speecl of a satellite in orbit can be deduced from a known state vector E formed from six orbital parameters such as, for example:
- the major half-axis of the orbit, usually designated a;
- the eccentricity vector defined in the plane of the orbit by its coordinates:
ex = e.cos . ey = e.sin~
where e is the eccentricity of the orbit (dimensionless parameter) and ~ is the argu~ent of the perigee;
- the inclination vector of the orbit defined by its coordinates:
~ ix = i.cosS2 . iy = i.sinS2 where i is the inclination tin degrees), that is to say the angle of the plane of the orbit relative to the terrestrial equator, and ~ is the right ascension of the ascendant node (which designates the orientation of the line of nodes relative to any predetermined inertial frame of reference);
- the mean longitude lm If the satellite was subjected only to the gravity field of a homogeneous and perfectly spherical' Earth the orbital parameters of the state vector E would remain constant (a = 42.164 km, ex = ey = 0 and ix = iy = 0, lm = parking - or set point longitude) and the satellite would remain strictly geostationary.
~k .
.~ . ~ ,. .
~: .
The position and speecl of a satellite in orbit can be deduced from a known state vector E formed from six orbital parameters such as, for example:
- the major half-axis of the orbit, usually designated a;
- the eccentricity vector defined in the plane of the orbit by its coordinates:
ex = e.cos . ey = e.sin~
where e is the eccentricity of the orbit (dimensionless parameter) and ~ is the argu~ent of the perigee;
- the inclination vector of the orbit defined by its coordinates:
~ ix = i.cosS2 . iy = i.sinS2 where i is the inclination tin degrees), that is to say the angle of the plane of the orbit relative to the terrestrial equator, and ~ is the right ascension of the ascendant node (which designates the orientation of the line of nodes relative to any predetermined inertial frame of reference);
- the mean longitude lm If the satellite was subjected only to the gravity field of a homogeneous and perfectly spherical' Earth the orbital parameters of the state vector E would remain constant (a = 42.164 km, ex = ey = 0 and ix = iy = 0, lm = parking - or set point longitude) and the satellite would remain strictly geostationary.
~k .
.~ . ~ ,. .
~: .
2~ 2~
However, because of disturbances due in particular to the non-spherical shape of the Earth and the non-homegenous nature of terrestrial gravity, the attraction of heavenly bodies such as the Sun and the Moon and solar pressure forces, the orbital parameters change slowly.
The function o~ a geostationary satellite requires in practice that it be held in a narrow window in terms of longitude and latitude (with a width 10typically between 0.05 and 0.1 degree). This requires correction (or "stationkeeping") manoeuvres which are currently computed and trans~itted to the satellite by one or more control stations on the ground and based on measurements by means of antennas on the ground.
15Stationkeeping therefore requires at present a permanently manned infrastructure on the ground (24 hours a day, 365 days a year), resulting in high satellite operating costs. This problem is co~pounded by possible problems with the availability of a suitable location for constructing a control centre and with the need to make it secure.
The object of the invention is to enable autono~ous computation on board the satellite of the orbital parameters and the stationkeeping manoeuvres to be executed so that it is possible to dispense with continuous assistance from the ground and to involve a control centre on the ground only occasionall~.
To this end the invention proposes a station-keeping method for a satellite in geostationary orbit characterised in that:
- there are determined at the same time the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle a2 between the satellite-Pole Star direction and the satellite-Earth direction, .
.
:' ' ~ ' ' - Z~Q2~
- there is deduced therefrom a state vector E consis-ting in orbital parameters by the formula:
Z = H.E + C.B
where:
. Z is a measurement vector the components of which are deduced from the angles al and a . H is a measuring matrix, . C is a bias sensitivity matrix of the form l Xs Ys where Xs is a term corresponding to a O O Yp period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, . B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, - stationkeeping manoeuvres are determined consequently and applied by means of thrusters.
Thus in accordance with the invention the position of the satellite in space is characterised by the angles of the Sun and the Pole Star to the Earth as seen from the satellite.
The angles al and a 2 are preferably determined from measurements effected by at least one terrestrial detector , a plurality of solar detectors and a stellar detector oriented towards the North along the South-North axis of the satellite, the measurements from these detectors being filtered separately so as to obtain for these measurements an overall time-delay (detector intrinsic time-delay plus filter time-delay) that is exactly the same.
The invention also proposes a stationkeeping ~ . ' :, ..
, . ' '' ~"'', .
. ~
2~
system comprising:
- at least one terrestrial detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay plus filter time-delay) has a predetermined value, - a plurality of solar detectors provided with a filter . adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, - at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, - a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle a2 between the satellite-Pole Star direction and the satellite-Earth direction, - an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
- a bias computation auxiliary unit connected to the output of the position computation unit , to a telemetry unit in communication with the ground an~
the output of which is connected to the orbital parameter computation unit, adapted to determine and memorise the bias vector B by comparing the vector E
with corresponding measurements from the ground, and - a control unit connected to stationkeeping ,.~. .,. ^ - :
,; , .
.:: ; , . .
.. , . . ~ . .. . .
, , ~ , . .
thrusters.
In one aspect, the present invention provides a stationkeeping method for a satellite in geostationary orbit comprising the steps of: determining at the same time the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle 2 between the satellite-Pole Star direction and the satellite-Earth direction, deducing therefrom a state vector E consisting in orbital parameters by the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles al and 2~ H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form '~Xs YS ~
O O Y
`
where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, determining stationkeeping manoeuvers from vector E and applying same by means of thrusters.
In another aspect, the present invention provides a stationkeeping system comprising: at least one detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay 5a 200021 4 s plus filter time-delay) has a predetermined value, a plurality of solar detectors provided wi~h a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle 1 between the satellite-Sun direction and the satellite-Earth direction and the angle 2 between the satellite-Pole Star direction and the satellite-Earth direction, an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles ~1 and ~2~ H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form ~Xs Ys ~ 0 Y~ .
where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a pericd substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias 3s vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, a bias - .
~ -20002 1 ~
computation auxiliary unit connected to the output of the position computation unit and to a telemetry unit in communication with the ground, the output of said bias computation auxiliary unit being connected to the orbital parameter computation unit, adapted to determine and memorize the bias vector B by comparing the vector E with corresponding measurements from the ground, and a control unit connected to stationkeeping thrusters.
According to preferred features:
- an attitude determination unit is connected to the outputs of the terrestrial, solar and stellar detectors and to the output of the orbital parameter computation unit, - the orbital parameter computation unit is a KALMAN filter, - the bias computation auxiliary unit is a FRIEDLAND filter.
Objects, characteristics and advantages of the invention will emerge from the following description given by way of non-limiting example only with reference to the appended drawings in which:
- figure 1 is a schematic view of a satellite equipped with an autonomous orbit control system in accordance with the invention;
- figure 2 is a simplified schematic view of this autonomous geostationary orbit control system; and - figure 3 shows the simulation over one year of the path of a satellite of this kind.
Figure 1 shows schematically a satellite 1 in an orbit 2 about the Earth.
The satellite 1 is conventionally associated with a frame of reference X Y Z in which the X axis is tangential to the orbit 2 and oriented in the direction in which the orbit is travelled (from West to East) and the Z axis is directed towards the Earth (geocentric axis).
' ' .' . :~
- ,~
- ' .~ .
.
5c 20002 1 4 The satellite comprises a platform carrying solar panels 3, reflectors 4 and propulsion thrusters of any appropriate known type.
Also in the known way the platform of this satellite comprises at least one terrestrial sensor oriented towards the Earth and schematically represented .
: , ~. . .
: ~ : . - - .
',,, . . ~:
at Tl together with a plurality of solar detectors distributed in a plane parallel to the plane o the X
and Z axes and adapted to face the Sun successively as the satellite completes its orbit Various configurations are known; to give an example, there are in this instance three solar detectors Sl. S2, S3 with one sensor S2 disposed on the side facing towards the Earth and the other two sensors Sl and S3 n the edges opposite this side.
In accordance with the invention the platform of the satellite is also provided, in this instance on its North side (opposite the Y axis), with a stellar detec-tor P of any appropriate known type oriented towards the North along the South-North axis. This is a detector chosen from the SODERN or GALILEO range, for example.
As is kno~n, there are currently important differences between the aforementioned three types of detectors, even if they are all optical detectors.
Firstly, the solar detectors Sl through S3l terrestrial detector Tl and stellar detector P are classified in this order by decreasing incident radiated power.
Also, these known detectors are sensitive to different radiation spectra; the terrestrial detectors are sensitive to infra-red radiation from-the Earth, the solar detectors are formed of photo-electric cells and the stellar detectors are based on rows or arrays of charge-coupled devices (CCD).
Finally, these detectors conventionally comprise integrated filters introducing different time-delays.
The invention exploits the fact that the Pole Star is the only star of its magnitude to be always in the field of view of a stellar detector mounted on the North side of a satellite so that is recognition does not require sophisticated software.
~,. . ~ . : .
,. . :
..
. ' . ' '~. " ' z~
The terrestrial detector Tl~ solar detectors Sl through S3 and stellar detector P make it possible to determine at any time the angular orientation of the Earth, the Sun and the Pole Star in the frame of reference related to the satellite.
Figure 2 shows the structure of the autonomous geostationary orbit control system.
Associated with the detectors Tlr S1 through S3 and P are parallel ~ilter units 10, 11 and 12 in a filter and synchronisation device the outputs of which are connected to a position computation unit 13 followed by an orbital parameter computation unit 14 with which is associated a calibration auxiliary unit 15 also connected to the position computation unit 13 and connected to a telemetry unit 16 in communication with the ground. The output of the orbital parameter computation unit is connected to a computation and command unit 17 controlling the stationkeeping thrusters 18. An attitude determination unit 19 is connected to the unfiltered outputs of the detectors Tl~ S
through S3 and P and to the output of the orbital parameters computation unit 14. In practice a unit 20 is associated with the computation unit 14 to store the orbital parameters computed by the unit 14.
In practice the computation units 13, 14, 15, 17, 18 may be integrated into the onboard computer of the satellite.
In accordance with the invention the position of the satellite in space is characterised by measuring at the same time the respective angular offsets al and a2 between the Sun and the Earth and between the Pole Star and the Earth as seen from the satellite.
The method of computing the angles al and a2, which computation is performed by the computation unit 13, is within the normal competence of those .. .. _ .. . . , . . , ~ . .
: ~ . . . . . .
. , , ,, , .. , . . ,. . ` .
.,, : ., ,:: . :;
-` Z~ Z~.9t skilled in the art.
In outline, these computations may be summarised as follows, where RT~ Rs and Rp denote frames of reference respectively related to the terrestrial detector, to the solar detectors and to the stellar detector. It is known how to define ~atrices Ps and Pp for passing from the frames of reference Rs and Rp to the frame of reference RT. The measurement from the solar and stellar detectors consists of a unit vector Xs or Xp directed from the satellite towards the Sun or the Pole Star.
Each of these vectors may be written in the frame of reference RT:
XS = PS 1 . Xs Xp = p -1 ~
and the navigation angles ~1 and ~2 deduced from:
cos al = ~'XT
cos ~2 = XP'XT
In practice, to obtain the measured angles ~1 and a2 at exactly the same time the invention proposes during a first stage (at 10, 11 and 12) to filter differently the raw measurements from the - various detectors so as not only to obtain an appropriate predetermined residual noise level but also to add a time-delay differing from one detector to another and such that the overall time-delay associated with each detector (or group o~ detectors of the same kind), that is to say the sum of the time-delay inherent to each detector and the additional time-del-ay introduced by the associated filter 10, 11 or 12, is the same for all the measurements. This synchronises the measurements which eliminates their influence on the attitude, which is subject to variations that can be very fast.
The overall time-delay is chosen as equal to ten seconds for example.
- 2~Q~
It is therefore possible to process at 13 a pair of navigation angles with only a moderate period, typically every ten minutes, which does not represent a significant increase in load for the onboard computer.
Eet Z be a measurement vector constituted ~rom the navigation angles ~1 and a~2-The computation unit 14 is adapted to deduce the state vector E combining the orb~ital parameters from a formula of the type:
Z = H.E + C.B
where H is a measurement matrix, B is a bias vector, and C is a bias sensitivity matrix of the type:
/XS Ys \
- C= ..
O O Yp/
where Xs is a term corresponding to a period approximating one day, YS iS a term representing a period approximating one year and Yp is a term representing a period approximating the apparent movement of the Pole Star.
Determining the components of the measurement 25 matrix H is within the normal competence of those skilled in the art, given the definitions selected for the vectors Z and E.
The same goes for the components of matrix C.
In the following example there have been chosen for the vectors 2 and E definitions slightly different to those given hereinabove:
_ Z = (COS Xl - COs Xlsyn, ' X~ X2syn) where the subscript "syn" is associated with the value of the parameter xl r x2 for the satellite assumed to be in an ideal, non-disturbed geosynchronous orbit : :
: -- ~
and where the superscript "T" indicates a vector or matrix transposition;
- E = (~, ~, ex, ey, ix, iy)T
where ~ = a-aSyn ~ = ~ + b~ + M-lm ex = e.cos (n+~) ey - e.sin (a+~) ix = i.cosn iy = i.sinn with M denoting the mean anomaly, e denoting the eccentricity of the orbit, ~ denoting the argument of the perigee; and a denoting the major half-axis of the orbit;
- the matrix H is then writtenrH~
~H2J
Hl = 1 0,xl.sin 1 - yl.cos 1, (l-cos 21).xl - sin 21.yl, -xl.sin 21 + (1 ~ cos 21).yl, -~ .sin 1, z .cos 1]
and H2 = 1 [O,x2.sin 1 - y2.cos l,(l-cos 21).x -sin 21-y , sinx 2syn -x2.sin 21 + (l+cos 21).y2,-z2.sin 1,z2.cosl]
where 1 = lm (parking longitude) + ~ (sidereal time), (xl, Yl~ Zl) is the unit vector (Earth-Sun centre) and (x2, Y2~ Z2) is the unit vector (Earth-Pole Star centre).
The parameters Xs~ Ys and Yp of matrix C
correspond to the coordinates of the sa~.e unit~ vectors if, changing the previous notation, these coordinates are written (Xs~ Ys~ Zs) for the Sun and (Xp, Yp, Zp) for the Pole Star.
As already explained, the form of H depends on the precise form chosen for Z and E.
- The object of the bias vector B is to take into account internal errors and detector alignment errors.
35Even if these errors were carefully calibrated on the ground, the environment and the launch process would be likely to alter them.
Also, it is currently recognised that it is impossible to calibrate them individually in flight because these biases or offsets cannot be observed.
,., : ~ - '' ' ' ' . , .
' ,, - ':
.
: .
. .
2~ 2~
Nevertheless, in the context of the invention it has been noted that it is possible to evaluate the overall effect of all these biases on the two navigation angles by means of a bias sensitivity matrix of the aforementioned type.
According to the invention, the procedure for calibrating the biases is as follows: the position of the satellite and its orbital parameters are determined regularly on the ground (for example, after the satellite is placed in orbit and then once a year) and these orbit parameters and/or the associated navigation angles are then uploaded to the satellite by the telemetry unit 16. The parameters are then processed together with the value Z deduced from the measurements from the detectors by the auxiliary computation unit 15 which deduces from them the new value of the bias vector B. The bias vector, once estimated at the end of a period of around ten days, is stored in the unit 15 until the next calibration. The computation function of the unit 15 may then be deactivated so as not to load unnecessaril~ the onboard computer.
The computation unit 14 is in practice in orbit estimating filter, a KALMAN filter, for example, the equations for which are well known to those skilled in the art.
The computation auxiliary unit 15 is in practice a filter, for example a filter of the FRIEDLAND type the equations of which are also well known. They make it possibe to evaluate B from several measurements.
Simulations carried out over one year have made it possible to verify that it is possible to maintain in this way a satellite within a window of 0.05 degree in longitude and latitude (see figure 3) with the following accuracy (three times the mean standard deviation) for the orbital parameters:
- : - , : :
. . - : . : .
- ,:
- : .
~ 12 200021 4 at a = 50 m at ex = ~3.5)10-5 at ey = (3.5)10-5 at ix and iy = 0.0001 dlegrees at longitude = 0.005 de~rees The raw measurements from the detectors (before filtering, and therefore corresponding to slightly different times) can be used by the unit 18 to determine the three attitude angles of the satellite (yaw, roll, pitch) from the orbital parameters supplied by the unit 14 and therefore autonomously (the time constants of the filters integrated into the detectors are usually around a few tenths of a second for the solar detectors, 0.5 seconds for the terrestrial detector and 0.1 to 0.5 seconds for the stellar detector).
The figure 2 system therefore constitutes an entirely autonomous system for determining the orbit and the attitude.
It goes without saying that the foregoing 2~ description has been given by way of non-limiting example only and that numerous variations thereon may ~e proposed by those skilled in the art without departing from the scope of the invention, in particular with regard to the arrangement of the solar detectors.
A
,.; ~. . . .
. .` . . ` .
. ..
. . .
.. ` :
.
... . .
`. .
However, because of disturbances due in particular to the non-spherical shape of the Earth and the non-homegenous nature of terrestrial gravity, the attraction of heavenly bodies such as the Sun and the Moon and solar pressure forces, the orbital parameters change slowly.
The function o~ a geostationary satellite requires in practice that it be held in a narrow window in terms of longitude and latitude (with a width 10typically between 0.05 and 0.1 degree). This requires correction (or "stationkeeping") manoeuvres which are currently computed and trans~itted to the satellite by one or more control stations on the ground and based on measurements by means of antennas on the ground.
15Stationkeeping therefore requires at present a permanently manned infrastructure on the ground (24 hours a day, 365 days a year), resulting in high satellite operating costs. This problem is co~pounded by possible problems with the availability of a suitable location for constructing a control centre and with the need to make it secure.
The object of the invention is to enable autono~ous computation on board the satellite of the orbital parameters and the stationkeeping manoeuvres to be executed so that it is possible to dispense with continuous assistance from the ground and to involve a control centre on the ground only occasionall~.
To this end the invention proposes a station-keeping method for a satellite in geostationary orbit characterised in that:
- there are determined at the same time the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle a2 between the satellite-Pole Star direction and the satellite-Earth direction, .
.
:' ' ~ ' ' - Z~Q2~
- there is deduced therefrom a state vector E consis-ting in orbital parameters by the formula:
Z = H.E + C.B
where:
. Z is a measurement vector the components of which are deduced from the angles al and a . H is a measuring matrix, . C is a bias sensitivity matrix of the form l Xs Ys where Xs is a term corresponding to a O O Yp period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, . B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, - stationkeeping manoeuvres are determined consequently and applied by means of thrusters.
Thus in accordance with the invention the position of the satellite in space is characterised by the angles of the Sun and the Pole Star to the Earth as seen from the satellite.
The angles al and a 2 are preferably determined from measurements effected by at least one terrestrial detector , a plurality of solar detectors and a stellar detector oriented towards the North along the South-North axis of the satellite, the measurements from these detectors being filtered separately so as to obtain for these measurements an overall time-delay (detector intrinsic time-delay plus filter time-delay) that is exactly the same.
The invention also proposes a stationkeeping ~ . ' :, ..
, . ' '' ~"'', .
. ~
2~
system comprising:
- at least one terrestrial detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay plus filter time-delay) has a predetermined value, - a plurality of solar detectors provided with a filter . adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, - at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, - a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle a2 between the satellite-Pole Star direction and the satellite-Earth direction, - an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
- a bias computation auxiliary unit connected to the output of the position computation unit , to a telemetry unit in communication with the ground an~
the output of which is connected to the orbital parameter computation unit, adapted to determine and memorise the bias vector B by comparing the vector E
with corresponding measurements from the ground, and - a control unit connected to stationkeeping ,.~. .,. ^ - :
,; , .
.:: ; , . .
.. , . . ~ . .. . .
, , ~ , . .
thrusters.
In one aspect, the present invention provides a stationkeeping method for a satellite in geostationary orbit comprising the steps of: determining at the same time the angle ~1 between the satellite-Sun direction and the satellite-Earth direction and the angle 2 between the satellite-Pole Star direction and the satellite-Earth direction, deducing therefrom a state vector E consisting in orbital parameters by the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles al and 2~ H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form '~Xs YS ~
O O Y
`
where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, determining stationkeeping manoeuvers from vector E and applying same by means of thrusters.
In another aspect, the present invention provides a stationkeeping system comprising: at least one detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay 5a 200021 4 s plus filter time-delay) has a predetermined value, a plurality of solar detectors provided wi~h a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle 1 between the satellite-Sun direction and the satellite-Earth direction and the angle 2 between the satellite-Pole Star direction and the satellite-Earth direction, an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles ~1 and ~2~ H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form ~Xs Ys ~ 0 Y~ .
where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a pericd substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias 3s vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, a bias - .
~ -20002 1 ~
computation auxiliary unit connected to the output of the position computation unit and to a telemetry unit in communication with the ground, the output of said bias computation auxiliary unit being connected to the orbital parameter computation unit, adapted to determine and memorize the bias vector B by comparing the vector E with corresponding measurements from the ground, and a control unit connected to stationkeeping thrusters.
According to preferred features:
- an attitude determination unit is connected to the outputs of the terrestrial, solar and stellar detectors and to the output of the orbital parameter computation unit, - the orbital parameter computation unit is a KALMAN filter, - the bias computation auxiliary unit is a FRIEDLAND filter.
Objects, characteristics and advantages of the invention will emerge from the following description given by way of non-limiting example only with reference to the appended drawings in which:
- figure 1 is a schematic view of a satellite equipped with an autonomous orbit control system in accordance with the invention;
- figure 2 is a simplified schematic view of this autonomous geostationary orbit control system; and - figure 3 shows the simulation over one year of the path of a satellite of this kind.
Figure 1 shows schematically a satellite 1 in an orbit 2 about the Earth.
The satellite 1 is conventionally associated with a frame of reference X Y Z in which the X axis is tangential to the orbit 2 and oriented in the direction in which the orbit is travelled (from West to East) and the Z axis is directed towards the Earth (geocentric axis).
' ' .' . :~
- ,~
- ' .~ .
.
5c 20002 1 4 The satellite comprises a platform carrying solar panels 3, reflectors 4 and propulsion thrusters of any appropriate known type.
Also in the known way the platform of this satellite comprises at least one terrestrial sensor oriented towards the Earth and schematically represented .
: , ~. . .
: ~ : . - - .
',,, . . ~:
at Tl together with a plurality of solar detectors distributed in a plane parallel to the plane o the X
and Z axes and adapted to face the Sun successively as the satellite completes its orbit Various configurations are known; to give an example, there are in this instance three solar detectors Sl. S2, S3 with one sensor S2 disposed on the side facing towards the Earth and the other two sensors Sl and S3 n the edges opposite this side.
In accordance with the invention the platform of the satellite is also provided, in this instance on its North side (opposite the Y axis), with a stellar detec-tor P of any appropriate known type oriented towards the North along the South-North axis. This is a detector chosen from the SODERN or GALILEO range, for example.
As is kno~n, there are currently important differences between the aforementioned three types of detectors, even if they are all optical detectors.
Firstly, the solar detectors Sl through S3l terrestrial detector Tl and stellar detector P are classified in this order by decreasing incident radiated power.
Also, these known detectors are sensitive to different radiation spectra; the terrestrial detectors are sensitive to infra-red radiation from-the Earth, the solar detectors are formed of photo-electric cells and the stellar detectors are based on rows or arrays of charge-coupled devices (CCD).
Finally, these detectors conventionally comprise integrated filters introducing different time-delays.
The invention exploits the fact that the Pole Star is the only star of its magnitude to be always in the field of view of a stellar detector mounted on the North side of a satellite so that is recognition does not require sophisticated software.
~,. . ~ . : .
,. . :
..
. ' . ' '~. " ' z~
The terrestrial detector Tl~ solar detectors Sl through S3 and stellar detector P make it possible to determine at any time the angular orientation of the Earth, the Sun and the Pole Star in the frame of reference related to the satellite.
Figure 2 shows the structure of the autonomous geostationary orbit control system.
Associated with the detectors Tlr S1 through S3 and P are parallel ~ilter units 10, 11 and 12 in a filter and synchronisation device the outputs of which are connected to a position computation unit 13 followed by an orbital parameter computation unit 14 with which is associated a calibration auxiliary unit 15 also connected to the position computation unit 13 and connected to a telemetry unit 16 in communication with the ground. The output of the orbital parameter computation unit is connected to a computation and command unit 17 controlling the stationkeeping thrusters 18. An attitude determination unit 19 is connected to the unfiltered outputs of the detectors Tl~ S
through S3 and P and to the output of the orbital parameters computation unit 14. In practice a unit 20 is associated with the computation unit 14 to store the orbital parameters computed by the unit 14.
In practice the computation units 13, 14, 15, 17, 18 may be integrated into the onboard computer of the satellite.
In accordance with the invention the position of the satellite in space is characterised by measuring at the same time the respective angular offsets al and a2 between the Sun and the Earth and between the Pole Star and the Earth as seen from the satellite.
The method of computing the angles al and a2, which computation is performed by the computation unit 13, is within the normal competence of those .. .. _ .. . . , . . , ~ . .
: ~ . . . . . .
. , , ,, , .. , . . ,. . ` .
.,, : ., ,:: . :;
-` Z~ Z~.9t skilled in the art.
In outline, these computations may be summarised as follows, where RT~ Rs and Rp denote frames of reference respectively related to the terrestrial detector, to the solar detectors and to the stellar detector. It is known how to define ~atrices Ps and Pp for passing from the frames of reference Rs and Rp to the frame of reference RT. The measurement from the solar and stellar detectors consists of a unit vector Xs or Xp directed from the satellite towards the Sun or the Pole Star.
Each of these vectors may be written in the frame of reference RT:
XS = PS 1 . Xs Xp = p -1 ~
and the navigation angles ~1 and ~2 deduced from:
cos al = ~'XT
cos ~2 = XP'XT
In practice, to obtain the measured angles ~1 and a2 at exactly the same time the invention proposes during a first stage (at 10, 11 and 12) to filter differently the raw measurements from the - various detectors so as not only to obtain an appropriate predetermined residual noise level but also to add a time-delay differing from one detector to another and such that the overall time-delay associated with each detector (or group o~ detectors of the same kind), that is to say the sum of the time-delay inherent to each detector and the additional time-del-ay introduced by the associated filter 10, 11 or 12, is the same for all the measurements. This synchronises the measurements which eliminates their influence on the attitude, which is subject to variations that can be very fast.
The overall time-delay is chosen as equal to ten seconds for example.
- 2~Q~
It is therefore possible to process at 13 a pair of navigation angles with only a moderate period, typically every ten minutes, which does not represent a significant increase in load for the onboard computer.
Eet Z be a measurement vector constituted ~rom the navigation angles ~1 and a~2-The computation unit 14 is adapted to deduce the state vector E combining the orb~ital parameters from a formula of the type:
Z = H.E + C.B
where H is a measurement matrix, B is a bias vector, and C is a bias sensitivity matrix of the type:
/XS Ys \
- C= ..
O O Yp/
where Xs is a term corresponding to a period approximating one day, YS iS a term representing a period approximating one year and Yp is a term representing a period approximating the apparent movement of the Pole Star.
Determining the components of the measurement 25 matrix H is within the normal competence of those skilled in the art, given the definitions selected for the vectors Z and E.
The same goes for the components of matrix C.
In the following example there have been chosen for the vectors 2 and E definitions slightly different to those given hereinabove:
_ Z = (COS Xl - COs Xlsyn, ' X~ X2syn) where the subscript "syn" is associated with the value of the parameter xl r x2 for the satellite assumed to be in an ideal, non-disturbed geosynchronous orbit : :
: -- ~
and where the superscript "T" indicates a vector or matrix transposition;
- E = (~, ~, ex, ey, ix, iy)T
where ~ = a-aSyn ~ = ~ + b~ + M-lm ex = e.cos (n+~) ey - e.sin (a+~) ix = i.cosn iy = i.sinn with M denoting the mean anomaly, e denoting the eccentricity of the orbit, ~ denoting the argument of the perigee; and a denoting the major half-axis of the orbit;
- the matrix H is then writtenrH~
~H2J
Hl = 1 0,xl.sin 1 - yl.cos 1, (l-cos 21).xl - sin 21.yl, -xl.sin 21 + (1 ~ cos 21).yl, -~ .sin 1, z .cos 1]
and H2 = 1 [O,x2.sin 1 - y2.cos l,(l-cos 21).x -sin 21-y , sinx 2syn -x2.sin 21 + (l+cos 21).y2,-z2.sin 1,z2.cosl]
where 1 = lm (parking longitude) + ~ (sidereal time), (xl, Yl~ Zl) is the unit vector (Earth-Sun centre) and (x2, Y2~ Z2) is the unit vector (Earth-Pole Star centre).
The parameters Xs~ Ys and Yp of matrix C
correspond to the coordinates of the sa~.e unit~ vectors if, changing the previous notation, these coordinates are written (Xs~ Ys~ Zs) for the Sun and (Xp, Yp, Zp) for the Pole Star.
As already explained, the form of H depends on the precise form chosen for Z and E.
- The object of the bias vector B is to take into account internal errors and detector alignment errors.
35Even if these errors were carefully calibrated on the ground, the environment and the launch process would be likely to alter them.
Also, it is currently recognised that it is impossible to calibrate them individually in flight because these biases or offsets cannot be observed.
,., : ~ - '' ' ' ' . , .
' ,, - ':
.
: .
. .
2~ 2~
Nevertheless, in the context of the invention it has been noted that it is possible to evaluate the overall effect of all these biases on the two navigation angles by means of a bias sensitivity matrix of the aforementioned type.
According to the invention, the procedure for calibrating the biases is as follows: the position of the satellite and its orbital parameters are determined regularly on the ground (for example, after the satellite is placed in orbit and then once a year) and these orbit parameters and/or the associated navigation angles are then uploaded to the satellite by the telemetry unit 16. The parameters are then processed together with the value Z deduced from the measurements from the detectors by the auxiliary computation unit 15 which deduces from them the new value of the bias vector B. The bias vector, once estimated at the end of a period of around ten days, is stored in the unit 15 until the next calibration. The computation function of the unit 15 may then be deactivated so as not to load unnecessaril~ the onboard computer.
The computation unit 14 is in practice in orbit estimating filter, a KALMAN filter, for example, the equations for which are well known to those skilled in the art.
The computation auxiliary unit 15 is in practice a filter, for example a filter of the FRIEDLAND type the equations of which are also well known. They make it possibe to evaluate B from several measurements.
Simulations carried out over one year have made it possible to verify that it is possible to maintain in this way a satellite within a window of 0.05 degree in longitude and latitude (see figure 3) with the following accuracy (three times the mean standard deviation) for the orbital parameters:
- : - , : :
. . - : . : .
- ,:
- : .
~ 12 200021 4 at a = 50 m at ex = ~3.5)10-5 at ey = (3.5)10-5 at ix and iy = 0.0001 dlegrees at longitude = 0.005 de~rees The raw measurements from the detectors (before filtering, and therefore corresponding to slightly different times) can be used by the unit 18 to determine the three attitude angles of the satellite (yaw, roll, pitch) from the orbital parameters supplied by the unit 14 and therefore autonomously (the time constants of the filters integrated into the detectors are usually around a few tenths of a second for the solar detectors, 0.5 seconds for the terrestrial detector and 0.1 to 0.5 seconds for the stellar detector).
The figure 2 system therefore constitutes an entirely autonomous system for determining the orbit and the attitude.
It goes without saying that the foregoing 2~ description has been given by way of non-limiting example only and that numerous variations thereon may ~e proposed by those skilled in the art without departing from the scope of the invention, in particular with regard to the arrangement of the solar detectors.
A
,.; ~. . . .
. .` . . ` .
. ..
. . .
.. ` :
.
... . .
`. .
Claims (12)
1. Stationkeeping method for a satellite in geostationary orbit comprising the steps of:
determining at the same time the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, deducing therefrom a state vector E consisting in orbital parameters by the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles .alpha.1 and .alpha.2, H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, determining stationkeeping manoeuvers from vector E
and applying same by means of thrusters.
determining at the same time the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, deducing therefrom a state vector E consisting in orbital parameters by the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles .alpha.1 and .alpha.2, H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, determining stationkeeping manoeuvers from vector E
and applying same by means of thrusters.
2. Method according to claim 1 wherein said angles .alpha.1 and .alpha.2 are determined from measurements effected by at least one terrestrial detector, a plurality of solar detectors and a stellar detector oriented towards the North along the South-North axis of the satellite, said measurements from these detectors being filtered separately so as to obtain for these measurements an overall time-delay (detector intrinsic time-delay plus filter time-delay) that is exactly the same.
3. The method of claim 1 wherein the orbital parameters in the state vector E are six in number, respectively representative of the major half-axis of the orbit, two coordinates of the eccentricity vector in the plane of the orbit, two coordinates of the inclination vector of the orbit and the mean longitude of the satellite, and the components in said measurement vector are two in number.
4. Stationkeeping system comprising:
at least one detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay plus filter time-delay) has a predetermined value, a plurality of solar detectors provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles .alpha.1 and .alpha.2.
H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, a bias computation auxiliary unit connected to the output of the position computation unit and to a telemetry unit in communication with the ground, the output of said bias computation auxiliary unit being connected to the orbital parameter computation unit, adapted to determine and memorize the bias vector B by comparing the vector E with corresponding measurements from the ground, and a control unit connected to stationkeeping thrusters.
at least one detector provided with a filter adapted to generate a time-delay such that the overall time-delay (detector intrinsic time-delay plus filter time-delay) has a predetermined value, a plurality of solar detectors provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, at least one stellar detector oriented towards the North along the South-North axis of the satellite provided with a filter adapted to generate a time-delay such that the overall time-delay is equal to said predetermined value, a position computation unit connected to these filters and adapted to deduce from the filtered measurements from said detectors the angle .alpha.1 between the satellite-Sun direction and the satellite-Earth direction and the angle .alpha.2 between the satellite-Pole Star direction and the satellite-Earth direction, an orbital parameter computation unit connected to the output of the position computation unit and adapted to determine the state vector E of the satellite made up of orbital parameters according to the formula:
Z = H.E + C.B
where Z is a measurement vector the components of which are deduced from the angles .alpha.1 and .alpha.2.
H is a measuring matrix with a number of lines equal to the number of said components in Z and with a number of columns equal to the number of said orbital parameters in E, C is a bias sensitivity matrix of the form where Xs is a term corresponding to a period substantially equal to one day, Ys is a term corresponding to a period substantially equal to one year and Yp is a term characteristic of the movement of the Pole Star, B is a bias vector determined beforehand by comparison of the measured vector Z and measurements made on the ground, a bias computation auxiliary unit connected to the output of the position computation unit and to a telemetry unit in communication with the ground, the output of said bias computation auxiliary unit being connected to the orbital parameter computation unit, adapted to determine and memorize the bias vector B by comparing the vector E with corresponding measurements from the ground, and a control unit connected to stationkeeping thrusters.
5. System according to claim 4 further comprising an attitude determination unit connected to the outputs of the terrestrial, solar and stellar detectors and to the output of the orbital parameter computation unit.
6. System according to claim 5 wherein said orbital parameter computation unit is a KALMAN filter.
7. System according to claim 6 wherein said bias computation auxiliary unit is a FRIEDLAND filter.
8. System according to claim 5 wherein said bias computation auxiliary unit is a FRIEDLAND filter
9. System according to claim 4 wherein said orbital parameter computation unit is a KALMAN filter.
10. System according to claim 9 wherein said bias computation auxiliary unit is a FRIEDLAND filter.
11. System according to claim 4 wherein said bias computation auxiliary unit is a FRIEDLAND filter.
12. The method of claim 4 wherein the orbital parameters in the state vector E are six in number, respectively representative of the major half-axis of the orbit, two coordinates of the eccentricity vector in the plane of the orbit, two coordinates of the inclination vector of the orbit and the mean longitude of the satellite, and the components in said measurement vector are two in number.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR888813121A FR2637564B1 (en) | 1988-10-06 | 1988-10-06 | METHOD AND SYSTEM FOR AUTONOMOUS ORBIT CONTROL OF A GEOSTATIONARY SATELLITE |
FR8813.121 | 1988-10-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2000214A1 CA2000214A1 (en) | 1990-04-06 |
CA2000214C true CA2000214C (en) | 1993-12-21 |
Family
ID=9370753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002000214A Expired - Fee Related CA2000214C (en) | 1988-10-06 | 1989-10-05 | Autonomous orbit control method and system for a geostationary satellite |
Country Status (7)
Country | Link |
---|---|
US (1) | US5108050A (en) |
EP (1) | EP0363243B1 (en) |
JP (1) | JP2847302B2 (en) |
CA (1) | CA2000214C (en) |
DE (1) | DE68911830T2 (en) |
ES (1) | ES2047695T3 (en) |
FR (1) | FR2637564B1 (en) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19703629A1 (en) * | 1997-01-31 | 1998-08-06 | Daimler Benz Aerospace Ag | Method for autonomously determining the position of a satellite |
US5107434A (en) * | 1990-04-19 | 1992-04-21 | General Electric Company | Three-axis spacecraft attitude control using polar star sensor |
DE69109266T2 (en) * | 1990-08-22 | 1995-08-24 | Microcosm Inc | DEVICE FOR HOLDING A SATELLITE ON ITS ORBIT. |
FR2670886B1 (en) * | 1990-12-21 | 1994-07-01 | Aerospatiale | ATTITUDE REACQUISITION METHOD BY STAR RECOGNITION FOR 3-AXIS STABILIZED SATELLITE. |
FR2670746B1 (en) * | 1990-12-21 | 1993-04-16 | Aerospatiale | ATTITUDE MONITORING SYSTEM FOR 3-AXIS SATELLITE ,; IN PARTICULAR FOR OBSERVATION SATELLITE. |
DE4129630A1 (en) * | 1991-09-06 | 1993-05-06 | Deutsche Aerospace Ag, 8000 Muenchen, De | MEASURING ARRANGEMENT AND CONTROL SYSTEM FOR THE POSITION CONTROL OF A THREE-AXIS-STABILIZED SATELLITE AND RELATED MEASURING AND CONTROL PROCEDURES |
US5562266A (en) * | 1992-10-29 | 1996-10-08 | Aerospatiale Societe Nationale Industrielle | Rate gyro calibration method and apparatus for a three-axis stabilized satellite |
FR2697651B1 (en) * | 1992-10-29 | 1995-02-03 | Aerospatiale | Method and apparatus for calibrating the gyros of a stabilized 3-axis satellite. |
US5452869A (en) * | 1992-12-18 | 1995-09-26 | Hughes Aircraft Company | On-board three-axes attitude determination and control system |
RU2033949C1 (en) * | 1993-02-09 | 1995-04-30 | Севастиян Дмитриевич Гнатюк | Self-contained on-board control system for space vehicle |
US5452077A (en) * | 1993-12-09 | 1995-09-19 | Hughes Aircraft Company | Transient-free method of determining satellite attitude |
DE19518617A1 (en) * | 1995-05-24 | 1996-11-28 | Daimler Benz Aerospace Ag | Three-axis stabilized, earth-oriented satellite and associated sun and earth acquisition method using a magnetometer |
US6098929A (en) * | 1998-01-28 | 2000-08-08 | Falbel; Gerald | Three axis attitude readout system for geosynchronous spacecraft |
US6145790A (en) * | 1998-09-22 | 2000-11-14 | Hughes Electronics Corporation | Attitude determination system and method |
EP1076005B1 (en) * | 1999-08-13 | 2007-01-03 | Hughes Electronics Corporation | Spacecraft orbit control using orbit position feedback |
US6691033B1 (en) * | 2000-07-26 | 2004-02-10 | Hughes Electronics Corporation | System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system |
US7454272B1 (en) * | 2005-08-25 | 2008-11-18 | Raytheon Company | Geostationary stationkeeping method |
US8174433B1 (en) | 2009-02-26 | 2012-05-08 | Raytheon Company | Bias estimation and orbit determination |
US9091552B2 (en) * | 2011-10-25 | 2015-07-28 | The Boeing Company | Combined location and attitude determination system and methods |
FR2983981B1 (en) * | 2011-12-09 | 2014-08-22 | Thales Sa | METHOD AND SYSTEM FOR CONTROLLING A SET OF AT LEAST TWO SATELLITES, DESIGNED TO PROVIDE SERVICE ON A GEOSTATIONARY ORBIT, DELIVERING THIS SERVICE TO A NON-GEOSTATIONARY ORBIT |
US9073648B2 (en) * | 2013-02-15 | 2015-07-07 | The Boeing Company | Star tracker rate estimation with kalman filter enhancement |
CN111238567A (en) * | 2020-01-06 | 2020-06-05 | 航天恒星科技有限公司 | Satellite in-orbit health monitoring device |
CN112052560A (en) * | 2020-07-30 | 2020-12-08 | 上海航天控制技术研究所 | Closed-loop simulation design method for system-level on-board computer maintenance mode |
CN112319850B (en) * | 2020-09-30 | 2022-05-24 | 中国卫通集团股份有限公司 | Method and device for automatically keeping position of synchronous orbit satellite |
CN112526561B (en) * | 2020-11-27 | 2024-04-23 | 中国科学院国家天文台 | Method for prolonging forecasting period of two ephemeris of geostationary orbit communication satellite |
CN112769466B (en) * | 2020-12-22 | 2022-08-12 | 火眼位置数智科技服务有限公司 | Low-orbit satellite constellation configuration keeping method |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3417661A1 (en) * | 1983-05-13 | 1984-11-15 | Mitsubishi Denki K.K., Tokio/Tokyo | System for controlling the orientation of an artificial satellite |
US4617634A (en) * | 1983-06-28 | 1986-10-14 | Mitsubishi Denki Kabushiki Kaisha | Artificial satellite attitude control system |
JPS6171300A (en) * | 1984-09-13 | 1986-04-12 | 三菱電機株式会社 | Computer for angle of attitude of artificial satellite |
FR2583873B1 (en) * | 1985-06-20 | 1987-09-11 | Matra | METHOD AND DEVICE FOR INJECTING SATELLITE ONTO GEOSTATIONARY ORBIT WITH STABILIZATION FOLLOWING THE THREE AXES |
GB8616385D0 (en) * | 1986-07-04 | 1986-08-13 | Marconi Space Systems Ltd | Satellite attitude control |
FR2605427A1 (en) * | 1986-10-16 | 1988-04-22 | Centre Nat Etd Spatiales | POINTAGE OF SPATIAL PROBE ANTENNA TO EARTH |
GB8809247D0 (en) * | 1988-04-20 | 1988-05-25 | British Aerospace | Attitude recovery for spacecraft |
-
1988
- 1988-10-06 FR FR888813121A patent/FR2637564B1/en not_active Expired - Fee Related
-
1989
- 1989-09-19 ES ES89402567T patent/ES2047695T3/en not_active Expired - Lifetime
- 1989-09-19 DE DE89402567T patent/DE68911830T2/en not_active Expired - Lifetime
- 1989-09-19 EP EP89402567A patent/EP0363243B1/en not_active Expired - Lifetime
- 1989-10-03 US US07/416,694 patent/US5108050A/en not_active Expired - Fee Related
- 1989-10-05 JP JP1261230A patent/JP2847302B2/en not_active Expired - Lifetime
- 1989-10-05 CA CA002000214A patent/CA2000214C/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
DE68911830T2 (en) | 1994-04-28 |
EP0363243B1 (en) | 1993-12-29 |
US5108050A (en) | 1992-04-28 |
JP2847302B2 (en) | 1999-01-20 |
FR2637564A1 (en) | 1990-04-13 |
DE68911830D1 (en) | 1994-02-10 |
CA2000214A1 (en) | 1990-04-06 |
FR2637564B1 (en) | 1994-10-14 |
ES2047695T3 (en) | 1994-03-01 |
JPH02156312A (en) | 1990-06-15 |
EP0363243A1 (en) | 1990-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2000214C (en) | Autonomous orbit control method and system for a geostationary satellite | |
US5745869A (en) | Techniques for optimizing an autonomous star tracker | |
US5054719A (en) | Active three-axis attitude control system for a geostationary satellite | |
US6023291A (en) | Satellite camera attitude determination and image navigation by means of earth edge and landmark measurement | |
CA1266110A (en) | Star sightings by satellite for image navigation | |
Van Patten et al. | A possible experiment with two counter-orbiting drag-free satellites to obtain a new test of einstein's general theory of relativity and improved measurements in geodesy | |
US8706322B2 (en) | Method and computer program product for controlling inertial attitude of an artificial satellite by applying gyroscopic precession to maintain the spin axis perpendicular to sun lines | |
US9694917B1 (en) | Deployment and control algorithms for wheel cluster formations of satellites | |
US10427806B2 (en) | Deployment and control algorithms for wheel cluster formations of satellites | |
US7877173B2 (en) | Method and apparatus for determining a satellite attitude using crosslink reference signals | |
Bandikova et al. | Characteristics and accuracies of the GRACE inter-satellite pointing | |
Goswami et al. | Analysis of GRACE Follow-On Laser Ranging Interferometer derived inter-satellite pointing angles | |
Wang | Study on center of mass calibration and K-band ranging system calibration of the GRACE mission | |
Zeitlhöfler | Nominal and observation-based attitude realization for precise orbit determination of the Jason satellites | |
Markley et al. | Attitude control system conceptual design for geostationary operational environmental satellite spacecraft series | |
Fujita et al. | Development and ground evaluation of ground-target tracking control of microsatellite RISESAT | |
Paluszek et al. | Optical navigation system | |
Kachmar et al. | Space navigation applications | |
Wie et al. | Attitude and orbit control systems | |
Gadisa et al. | ETRSS-1 attitude and orbit control subsystem and on-orbit 1 performance analysis | |
Challa et al. | Effects of magnetometer calibration and maneuvers on accuracies of magnetometer-only attitude-and-rate determination | |
Ibrahim | Attitude and orbit control of small satellites for autonomous terrestrial target tracking | |
MARKLEY et al. | Attitude control system conceptual design for the GOES-N spacecraft series | |
Johnston-Lemke | High performance attitude determination and control for nanosatellites missions | |
Mueller et al. | Integrated communications and optical navigation system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |