CN101214860A - Method for self-determination choosing attitude determination mode during rail controlling course - Google Patents

Method for self-determination choosing attitude determination mode during rail controlling course Download PDF

Info

Publication number
CN101214860A
CN101214860A CNA2007103015901A CN200710301590A CN101214860A CN 101214860 A CN101214860 A CN 101214860A CN A2007103015901 A CNA2007103015901 A CN A2007103015901A CN 200710301590 A CN200710301590 A CN 200710301590A CN 101214860 A CN101214860 A CN 101214860A
Authority
CN
China
Prior art keywords
attitude
satellite
gyro
star sensor
determination
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CNA2007103015901A
Other languages
Chinese (zh)
Inventor
宗红
陈义庆
王淑一
黄江川
李铁寿
太萍
程莉
王寨
韩冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CNA2007103015901A priority Critical patent/CN101214860A/en
Publication of CN101214860A publication Critical patent/CN101214860A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

A method for autonomously selecting an attitude-determination mode during the orbital control process includes: (1) predicting the satellite inertia attitude according to gyro measurement data; (2) judging that whether the star sensor attitude correction is required to be induced: judging that whether the attitude angular velocity of the satellite exceeds the measurement scope of a gyro according to the measured information of the gyro; if not, continuously using the gyro to predict the inertia attitude of the satellite, transferring to the step (1) and predicting the inertia attitude of the satellite in the next period; if the attitude angular velocity of the satellite exceeds the threshold which is set according to the gyro measurement scope, a gyro transfinite mark is arranged; then judging that whether the attitude angular velocity of three axils of the satellite satisfies the working condition requirements of the star sensor, and if yes, transferring to the step (3) and if no, transferring to the step (1); (3)introducing the star sensor for the attitude correction. The method of the invention reduces the attitude error during the orbital control process, improves the orbital control precision and effectively ensures that the orbital control is accurately finished in time.

Description

The method of self-determination choosing attitude determination mode in the rail control process
Technical field
The present invention relates to the method that attitude is determined in a kind of spacecraft rail control process.Particularly bigger at rail control engine disturbance torque, cause measuring satellite angular velocities bigger, exceed under the gyro to measure scope situation, satellite can independently be introduced the method that the star sensor metrical information is accurately revised satellite attitude.
Background technology
Gyro attitude prediction method is comparatively commonly used in satellite attitude control process, but because the influence of its gyro wander makes it have only application within a short period of time could guarantee accuracy of attitude determination.And the scope of testing the speed of general high accuracy gyroscope is more limited, change greatly in measuring satellite angular velocities, exceed under the gyro to measure scope situation, only estimate that by gyro data satellite attitude can make the true attitude of estimation pose deviation satellite, influences the satellite control accuracy.The method of using gyro to estimate in the Satellite Orbit Maneuver process has been avoided comparatively complicated filtering algorithm, simplified on the star and realized, help improving the reliability of rail control, but the attitude misalignment that brings that is limited in scope owing to test the speed will make the firing attitude of satellite depart from needed nominal firing attitude, directly influence the rail control effect of satellite.For the satellites that control is had relatively high expectations concerning rail such as moon exploration and even survey of deep space, rail control error may cause the failure of whole task.
The method of gyro and star sensor integrated attitude determination is also comparatively commonly used in the three axis stabilization control of high precision satellite.But it is mainly used in the Steady-State Control process of satellite.When the attitude angular velocity that has when satellite is big, can produce adverse influence to the importance in star map recognition ability of star sensor, star sensor is output data normally, estimates attitude thereby can not revise satellite in time, still can make satellite estimate the actual attitude of pose deviation.Can influence rail control effect equally if in rail control process, use, cause rail control mission failure when serious.
When most geostationary orbits or middle low-orbit satellite transfer orbital control, attitude is determined only to utilize the metrical information of gyro to carry out.When the orbit maneuver engine disturbance torque is big, might cause satellite bigger attitude angular velocity to occur, when this cireular frequency surpasses the measurement range of gyro, only carrying out attitude prediction by gyro can cause satellite to determine that attitude and the actual attitude of satellite are not inconsistent, make the actual attitude of satellite constantly depart from firing attitude, influence rail control precision, may cause rail control failure when serious.
Summary of the invention
Technology of the present invention is dealt with problems: overcome the deficiencies in the prior art, attitude angular velocity is big during at the rail control that may occur, the situation that surpasses the gyro to measure scope, when having proposed that attitude angular velocity is in the gyro measurable range during a kind of rail control, utilize gyro data to estimate satellite attitude, occur exceeding the gyro to measure scope than great-attitude angle speed the time, the autonomous method of introducing the correction of star sensor attitude, promptly propose during a kind of rail control according to the measuring satellite angular velocities observed reading, the method of self-determination choosing attitude determination mode, reduce the attitude error in the rail control process, improve rail control precision, guarantee that effectively the rail control is timely, finish exactly.
The method of self-determination choosing attitude determination mode is meant in the rail control process celestial body cireular frequency hour in the rail control process, only selects gyro to estimate to carry out attitude to determine; When measuring satellite angular velocities is bigger, when surpassing the gyro to measure scope, independently select gyro to estimate to carry out attitude to determine with the method for star sensor attitude correction.
Technical solution of the present invention: the method for self-determination choosing attitude determination mode in the rail control process is characterized in that comprising:
(1) estimate satellite inertia attitude according to the gyro to measure data:
At first calculate the attitude angular velocity of satellite, calculate three absolute angle speed then,, estimate the inertia attitude of satellite at last according to absolute angle speed according to the metrical information of gyro;
(2) judge whether to need to introduce the correction of star sensor attitude:
Metrical information according to gyro, judge whether the attitude angular velocity of satellite has surpassed the measurement range of gyro, if do not exceed measurement range, then continue to utilize gyro to estimate satellite inertia attitude, change step (1) over to, carry out the satellite inertia attitude prediction of following one-period, if surpass the threshold value of setting according to the gyro to measure scope, then set the sign that gyro transfinites, judge then whether the attitude angular velocity of three on satellite satisfies the service conditions requirement of star sensor, when the service conditions that satisfies star sensor requires, then change step (3) over to, otherwise continue to adopt the gyro to measure data to estimate satellite inertia attitude, change (1) over to.
(3) introduce star sensor and carry out the attitude correction.
The present invention's advantage compared with prior art is:
(1) the present invention can independently select different method for determining posture for use according to actual conditions in rail control process, does not end the rail control easily, can guarantee that the rail control carries out smoothly.
(2) the present invention has guaranteed effectively that by the introducing of star sensor information high thrust rail control engine becomes the precision that attitude is determined in the rail process.Thereby in rail control process, carry out attitude control accurately and timely.
(3) this method is applicable to the high thrust transfer orbital control of spacecraft.Be particularly suitable in follow-up survey of deep space series satellite, using.Have good, the portable good characteristics of inheritance.
Description of drawings
Fig. 1 is the diagram of circuit of the inventive method;
Fig. 2 is rail control process attitude determination mode flow path switch figure;
Fig. 3 is the simulation curve of the inventive method;
Fig. 4 is the simulation curve of conventional approach;
Wherein, in the simulation curve of Fig. 3, Fig. 4, on behalf of qw (0), curve 2, curve 1 represent qw (1), curve 3 to represent qw (2), is respectively the vector part of the actual attitude quaternion of satellite; On behalf of q (0), curve 5, curve 4 represent q (1), curve 6 to represent q (2), is respectively the star upper estimate of satellite attitude quaternion vector part.
The specific embodiment
The measurement range that is installed in the gyro that tests the speed on the satellite is limited.And in the Satellite Orbit Maneuver process, because rocking the interference that brings etc., interference that rail control engine start brings and liquid fuel all may cause bigger measuring satellite angular velocities, when this cireular frequency surpasses the gyro to measure scope, gyro can't provide correct angular velocity information, thereby satellite can't obtain attitude information accurately, can not guarantee attitude control accuracy, thereby influence rail control precision.
The present invention takes to estimate in the attitude that rail control process at first adopts the take off data of gyro to carry out satellite, the measurement of real-time judge gyro output simultaneously, near saturation value the time, provide the attitude angular velocity sign that transfinites, after treating that measuring satellite angular velocities steadily (satisfies star sensor and measures requirement), introduce the method that the star sensor data are revised the attitude of gyro estimation, to reduce the saturated evaluated error of bringing of gyro.
As shown in Figure 1, concrete implementation step of the present invention is as follows:
(1) gyro is estimated the inertia attitude: estimate satellite inertia attitude according to the gyro to measure data.
A, gyro data switching rate integration
Utilize the metrical information of three gyros can obtain satellite three-axis attitude angular velocity information.With the rate is example, and according to the principle of work of rate, each sampling can obtain the angle step of gyro in this sampling period, through calculating the attitude angular velocity of three on satellite after suitably changing.
(k) the installation matrix B in the satellite body system of axes is calculated gyro output transition matrix A=B for numbering i, j according to three gyros that participate in deciding appearance -1, in conjunction with the measurement output △ g of three gyros i, △ g j, △ g k, can obtain satellite three-axis attitude angle increment information △ g x, △ g y, △ g zCalculate by following formula:
Δ g x Δ g y Δ g z = A Δ g i Δ g j Δ g k
B, three absolute angle speed calculation
According to the result of calculation of gyro data conversion, calculate the mean angular velocity in the sampling period: ω ^ = [ ω ^ x ω ^ y ω ^ z ] T . Note cutting gyroscope constant value drift amount =[ in the calculating through demarcating xyz] T, unit: radian/hour.The demarcation of gyroscope constant value drift can take ground demarcate in advance or star on method by other sensor information real-time calibrations such as star sensors obtain.
ω ^ x = Δ g x / Δt - b ^ x / 3600
ω ^ y = Δ g y / Δt - b ^ y / 3600
ω ^ z = Δ g z / Δt - b ^ z / 3600
C, inertia attitude prediction
According to three absolute angle speed that calculate, calculate the three axis angular rate increments of satellite in this cycle, adopt following formula to estimate the inertia attitude of satellite, with the quaternion form q - ^ = q ^ 1 q ^ 2 q ^ 3 q ^ 4 Provide.
Δ g ^ = ω ^ · Δt
q - ^ = q - ^ + 1 2 Eq ( q - ^ ) Δ g ^
If q ^ 4 < 0 , Then q &OverBar; ^ = - q &OverBar; ^
q - ^ = q - ^ / Norm ( q - ^ )
Δ t is the sampling period in the above-mentioned formula, and function Eq utilizes measuring satellite angular velocities information to estimate the computing formula of satellite attitude quaternion.
(2) judge whether to need to introduce the correction of star sensor attitude
Metrical information according to gyro, judge whether the attitude angular velocity of satellite has surpassed the measurement range of gyro, if do not exceed measurement range, then continue to utilize gyro to estimate satellite inertia attitude, change step (1) over to, carry out the satellite inertia attitude prediction of following one-period, if surpass the threshold value of setting according to the gyro to measure scope, then set the sign that gyro transfinites, judge then whether the attitude angular velocity of three on satellite satisfies the service conditions requirement of star sensor, when the service conditions that satisfies star sensor requires, then change step (3) over to, otherwise continue to adopt the gyro to measure data to estimate satellite inertia attitude, change (1) over to.
(3) star sensor attitude modification method
After introducing the star sensor correction, adopt conventional star sensor filtering algorithm to carry out real-time attitude correction.The filtering algorithm design gets final product according to the Kalman filtering principle of classics.
Provide the concrete application of this method in lunar exploration satellite below:
As shown in Figure 2, lunar spacecraft adopts 6 three floating gyros, and wherein any three gyros all can be determined satellite inertia attitude angular velocity.The measurement range of each gyro is-0.9 °/s~0.9 °/s, considers certain design margin in engineering is used, and choosing the threshold value that attitude angular velocity surpasses the gyro to measure scope is 0.8 °/s.
Using 3 medium accuracy star sensors to carry out the satellite three-axis attitude in the lunar spacecraft determines.Star sensor has under the 0.5 °/s cireular frequency situation and can normally discern star chart at celestial body, provides attitude data.In engineering is used, consider certain design margin, choose and introduce the star sensor data to carry out the threshold values of attitude correction be 0.35 °/s.
According to above-mentioned implementation step, carried out l-G simulation test at certain lunar exploration satellite, simulation curve is as shown in Figure 3.Simulation process is: the rail control begins, measuring satellite angular velocities big (1.5 °/s), gyro output is saturated, the actual attitude of Satellite Attitude Estimation value and satellite engenders deviation, causes attitude deviation target ignition attitude.Behind the simulation time 20s, introduce the quick correction of star, as can be seen, very rapid convergence is in true attitude in the satellite attitude valuation from curve, and attitude error is less than 0.001 (rad) behind the 100s.
In order to compare the effect that the inventive method reaches, Fig. 4 has provided the simulation curve of not introducing the star sensor correction in rail control process.The simulation process condition is consistent with said process, and satellite only carries out attitude prediction by gyro data.From simulation curve, owing to do not introduce star sensor filtering correction, the satellite attitude valuation departs from the true attitude of satellite all the time, and attitude error is about 0.08 (rad) behind the 100s, and this error will influence rail control precision.
System described above is a kind of situation of the present invention, and those skilled in the art can carry out under the situation of the present invention variously augmenting, improving and change not departing from according to different requirements and design parameters, and therefore, the present invention is widely.

Claims (4)

1. the method for self-determination choosing attitude determination mode in the rail control process is characterized in that comprising:
(1) estimate satellite inertia attitude according to the gyro to measure data:
At first calculate the attitude angular velocity of satellite, calculate three absolute angle speed then,, estimate the inertia attitude of satellite at last according to absolute angle speed according to the metrical information of gyro;
(2) judge whether to need to introduce the correction of star sensor attitude:
Metrical information according to gyro, judge whether the attitude angular velocity of satellite has surpassed the measurement range of gyro, if do not exceed measurement range, then continue to utilize gyro to estimate satellite inertia attitude, change step (1) over to, carry out the satellite inertia attitude prediction of following one-period, if surpass the threshold value of setting according to the gyro to measure scope, then set the sign that gyro transfinites, judge then whether the attitude angular velocity of three on satellite satisfies the service conditions requirement of star sensor, when the service conditions that satisfies star sensor requires, then change step (3) over to, otherwise continue to adopt the gyro to measure data to estimate satellite inertia attitude, change step (1) over to.
(3) introduce star sensor and carry out the attitude correction.
2. the method for self-determination choosing attitude determination mode in the rail control process according to claim 1, it is characterized in that: the gyro in the described step (1) is selected rate for use.
3. the method for self-determination choosing attitude determination mode in the rail control process according to claim 1, it is characterized in that: the inertia attitude in the described step (1) provides with the quaternion form.
4. the method for self-determination choosing attitude determination mode in the rail control process according to claim 1 is characterized in that: the star sensor in the described step (2) is a CCD visual star sensor, and it is that attitude angular velocity is smaller or equal to 0.35 degree/second that its service conditions requires.
CNA2007103015901A 2007-12-26 2007-12-26 Method for self-determination choosing attitude determination mode during rail controlling course Pending CN101214860A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2007103015901A CN101214860A (en) 2007-12-26 2007-12-26 Method for self-determination choosing attitude determination mode during rail controlling course

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2007103015901A CN101214860A (en) 2007-12-26 2007-12-26 Method for self-determination choosing attitude determination mode during rail controlling course

Publications (1)

Publication Number Publication Date
CN101214860A true CN101214860A (en) 2008-07-09

Family

ID=39621388

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2007103015901A Pending CN101214860A (en) 2007-12-26 2007-12-26 Method for self-determination choosing attitude determination mode during rail controlling course

Country Status (1)

Country Link
CN (1) CN101214860A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102424116A (en) * 2011-12-08 2012-04-25 中国空间技术研究院 Method for optimizing orbital transfer strategy of geostationary orbit satellite
CN101696885B (en) * 2009-11-05 2012-05-23 中国人民解放军国防科学技术大学 Method for improving data processing precision of star sensors
CN101696884B (en) * 2009-11-05 2012-08-08 中国人民解放军国防科学技术大学 Method for determining spatial attitude accuracy of satellite
CN101758933B (en) * 2009-12-30 2012-08-22 北京控制工程研究所 Attitude and orbit control method based on fore and after arrangement of engine
CN104085539A (en) * 2014-06-26 2014-10-08 北京控制工程研究所 Method for imaging calibration attitude control
CN106249590A (en) * 2016-08-09 2016-12-21 中国科学院软件研究所 The method that integrated self-adaptive Nano satellite attitude determines
CN109489661A (en) * 2018-11-02 2019-03-19 上海航天控制技术研究所 Gyro constant value drift estimation method when a kind of satellite is initially entered the orbit
CN110209185A (en) * 2019-06-26 2019-09-06 北京控制工程研究所 A kind of Spacecraft During Attitude Maneuver antihunt means using the quick information posture resetting of star
CN111610795A (en) * 2020-05-12 2020-09-01 北京控制工程研究所 Pseudo-inverse solvable minimum configuration attitude control thruster instruction distribution method
CN113720355A (en) * 2021-09-10 2021-11-30 北京控制工程研究所 Gyro output saturation autonomous diagnosis method and system
CN116374208A (en) * 2022-12-30 2023-07-04 中国科学院空间应用工程与技术中心 On-orbit inflation and deflation method and system for space station

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101696885B (en) * 2009-11-05 2012-05-23 中国人民解放军国防科学技术大学 Method for improving data processing precision of star sensors
CN101696884B (en) * 2009-11-05 2012-08-08 中国人民解放军国防科学技术大学 Method for determining spatial attitude accuracy of satellite
CN101758933B (en) * 2009-12-30 2012-08-22 北京控制工程研究所 Attitude and orbit control method based on fore and after arrangement of engine
CN102424116A (en) * 2011-12-08 2012-04-25 中国空间技术研究院 Method for optimizing orbital transfer strategy of geostationary orbit satellite
CN104085539A (en) * 2014-06-26 2014-10-08 北京控制工程研究所 Method for imaging calibration attitude control
CN104085539B (en) * 2014-06-26 2015-12-30 北京控制工程研究所 The attitude control method of imaging calibration
CN106249590B (en) * 2016-08-09 2019-03-19 中国科学院软件研究所 The method that integrated self-adaptive Nano satellite posture determines
CN106249590A (en) * 2016-08-09 2016-12-21 中国科学院软件研究所 The method that integrated self-adaptive Nano satellite attitude determines
CN109489661A (en) * 2018-11-02 2019-03-19 上海航天控制技术研究所 Gyro constant value drift estimation method when a kind of satellite is initially entered the orbit
CN109489661B (en) * 2018-11-02 2020-06-09 上海航天控制技术研究所 Gyro combination constant drift estimation method during initial orbit entering of satellite
CN110209185A (en) * 2019-06-26 2019-09-06 北京控制工程研究所 A kind of Spacecraft During Attitude Maneuver antihunt means using the quick information posture resetting of star
CN110209185B (en) * 2019-06-26 2022-02-01 北京控制工程研究所 Spacecraft attitude maneuver stabilization method by means of attitude reset of satellite sensitive information
CN111610795A (en) * 2020-05-12 2020-09-01 北京控制工程研究所 Pseudo-inverse solvable minimum configuration attitude control thruster instruction distribution method
CN111610795B (en) * 2020-05-12 2023-04-14 北京控制工程研究所 Pseudo-inverse solvable minimum configuration attitude control thruster instruction distribution method
CN113720355A (en) * 2021-09-10 2021-11-30 北京控制工程研究所 Gyro output saturation autonomous diagnosis method and system
CN113720355B (en) * 2021-09-10 2023-11-10 北京控制工程研究所 Gyroscope output saturation autonomous diagnosis method and system
CN116374208A (en) * 2022-12-30 2023-07-04 中国科学院空间应用工程与技术中心 On-orbit inflation and deflation method and system for space station
CN116374208B (en) * 2022-12-30 2023-11-07 中国科学院空间应用工程与技术中心 On-orbit inflation and deflation method and system for space station

Similar Documents

Publication Publication Date Title
CN101214860A (en) Method for self-determination choosing attitude determination mode during rail controlling course
CN100529667C (en) Star sensor attitude determination method at self-determination retrieve rail controlling fault
CN103134491B (en) GEO orbit transfer vehicle SINS/CNS/GNSS integrated navigation system
CN101788296B (en) SINS/CNS deep integrated navigation system and realization method thereof
CN101949703B (en) Strapdown inertial/satellite combined navigation filtering method
CN104344837B (en) Speed observation-based redundant inertial navigation system accelerometer system level calibration method
CN112414413B (en) Relative angular momentum-based angle-only maneuvering detection and tracking method
CN103017760B (en) A kind of highly elliptic orbit Mars probes are independently to fiery orientation method
EP2578995B1 (en) Modified Kalman filter for generation of attitude error corrections
CN101393025A (en) AUV combined navigation system non-tracing switch method
CN103792561B (en) A kind of tight integration reduced-dimensions filtering method based on GNSS passage difference
CN105891863A (en) High-constraint based extended Kalman filter (EKF) positioning method
CN103676941A (en) Satellite control system fault diagnosis method based on kinematics and dynamics model
CN103198187A (en) Track design method of deep space probe and based on differential modification
CN105241456A (en) Loitering munition high-precision combination navigation method
Yang et al. Two antennas GPS-aided INS for attitude determination
CN109489661A (en) Gyro constant value drift estimation method when a kind of satellite is initially entered the orbit
CN113447043B (en) GNSS-based satellite astronomical navigation system error autonomous calibration method and system
Li et al. Geostationary station-keeping with electric propulsion in full and failure modes
CN100588907C (en) Ship optical fiber gyroscope strapping system initial posture determination method
CN109974695A (en) The robust adaptive filtering method of surface ship navigation system based on the space Krein
CN105300407A (en) Marine dynamic starting method for uniaxial modulation laser gyro inertial navigation system
CN110667892B (en) Satellite despinning control method based on geomagnetic measurement
CN104634348A (en) Attitude angle calculation method in integrated navigation
CN112306075A (en) High-precision off-orbit reverse iterative guidance method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Open date: 20080709