CN101196398A - Spacecraft posture confirming method based on Euler-q algorithm and DD2 filtering - Google Patents
Spacecraft posture confirming method based on Euler-q algorithm and DD2 filtering Download PDFInfo
- Publication number
- CN101196398A CN101196398A CNA2007103012797A CN200710301279A CN101196398A CN 101196398 A CN101196398 A CN 101196398A CN A2007103012797 A CNA2007103012797 A CN A2007103012797A CN 200710301279 A CN200710301279 A CN 200710301279A CN 101196398 A CN101196398 A CN 101196398A
- Authority
- CN
- China
- Prior art keywords
- attitude
- delta
- spacecraft
- rightarrow
- euler
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
A spacecraft attitude fixing method of DD2 (Divided Difierence 2) filter and a counting method based on Euler-q (Euler-quaternion) are provided, which relates to a state estimation and attitude determination method based on vector observation. The invention is characterized in that: according to the reference ephemeris and the reference vector and the observation vector observed by a star sensor, the spacecraft attitude quaternion can be gotten through Euler-q counting method; through DD2 filter, on one hand, feedback compensation the attitude matrix obtained through resolving by gyro; on the other hand, revise the palstance output of three-axle gyro and conduct compensation to the constant drift of gyroscope. As for that, the current spacecraft attitude information when getting high precision can be achieved, the constant drift of gyroscope can be scaled on line. The invention is an unattended attitude determining method, which is characterized in high precision and good timeliness, being convenient and feasible, and can be used in the attitude fixing system of various spacecrafts.
Description
Technical field
The attitude that the present invention relates to a kind of spacecraft is determined method, particularly a kind of spacecraft attitude based on Euler-q (Euler-quaternion) algorithm and DD2 (Divided Difference 2) filtering is determined method, and the attitude that is used for various high-precision inertia/celestial combined navigation system is determined.
Background technology
It is to utilize spaceborne attitude sensor to measure resulting information that the attitude of spacecraft sets the tasks, through suitable processing, try to achieve the information that is fixed on a certain reference frame of spacecraft body coordinate system space, spacecraft is because the requirement of various tasks, need high-precision attitude information, this is one of gordian technique of spacecraft development and widespread use.
Inertia/astronomy, be that INS (Inertial Navigation System)/CNS (CelestialNavigation System) attitude and heading reference system is that a kind of autonomous attitude is fully determined navigational system, it utilizes the carrier angular velocity information of gyro to measure and the starlight information in the star sensor measurement visual field, on the basis of initial information, carry out attitude algorithm and filtering, attitude can be provided continuously, in real time, have that independence is strong, good concealment, be not subjected to advantage such as weather condition restriction, thereby be widely used in field such as space flight.The spacecraft attitude that uses is determined model at present, be state mostly with the attitude quaternion, because the degree of freedom of three-axis attitude is 3, therefore attitude quaternion exists redundant, if and the vector of getting attitude quaternion merely is the state of estimating, then must simplify, cause the precision of model to reduce the estimation model linearization.Be applied to spacecraft attitude and determine in the method, use determinacy vector observation merely, for example Davenport q method is owing to need computation of characteristic values and eigenmatrix, online dyscalculia, thereby be difficult to be applied.And on Davenport q method improved QuEsT (Quaternion Estimator) method, though on algorithm, obtained great improvement, still existed too to rely on the star sensor precision and can not revise gyroscopic drift and export the shortcoming of real-time attitude.SVD (the Singular ValueDecomposition) robustness is very strong but need svd, because svd is a complex work that operand is very big, is not very effective therefore.FOAM (A Fast Optimal MatrixAlgurithm) algorithm can not estimate that gyroscopic drift makes gyro output obtain on-line proving; And use other state estimation method, for example, EKF EKF (Extend Kalman Filter), because using the first order nonlinear of Taylor expansion handles, though can estimate nonlinear model, but it is not high to the nonlinear model estimated accuracy, UKF (Unsented Kalman Filter) then utilizes the sampled point of a series of approximate Gaussian distribution, carry out the recursion and the renewal of state and error covariance by the Unscented conversion, though precision increases, but calculate loaded down with trivial detailsly, real-time is not high, thereby is difficult to satisfy the requirement that attitude is determined real-time.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, provide a kind of spacecraft attitude to determine method based on Euler-q algorithm and DD2 filtering, the spacecraft attitude that this method is obtained by the observation of Euler-q algorithm computation vector, while is by the constant value drift of DD2 Filtering Estimation attitude error and gyro, high-accuracy posture information not only is provided, and has proofreaied and correct the output of gyro angular speed.
Technical solution of the present invention is: a kind of spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method, it is characterized in that using Euler-q algorithm and DD2 Filtering Estimation model, utilize star sensor and gyroscope, spacecraft attitude in the time of both can obtaining high-precision real, can estimate gyrostatic constant value drift again, its concrete steps are as follows:
(1) by the responsive spacecraft attitude of gyroscope, the output tri-axis angular rate, and use the angle increment algorithm to upgrade the attitude battle array, obtain the real-time attitude of spacecraft under geocentric inertial coordinate system: course angle , pitching angle theta, roll angle γ;
(2) utilize the responsive starlight of star sensor, pre-service obtains observing the coordinate of starlight vector s under the spacecraft body coordinate system through star chart; And through obtaining the coordinate of corresponding reference starlight vector v under geocentric inertial coordinate system with it after the identification of star chart coupling;
(3) use the Euler-q algorithm, utilize observation starlight vector and corresponding reference starlight vector with it, acquisition spacecraft attitude hypercomplex number;
(4) difference of the spacecraft attitude that resolves of the spacecraft attitude hypercomplex number that obtains with starlight observation and the middle gyroscope of step (1) is as observed quantity, and DD2 filtering is carried out in gyroscopic drift to the spacecraft attitude sum of errors, and quantity of state is: (Δ q
1Δ q
2Δ q
3b
1b
2b
3);
(5) the spacecraft attitude sum of errors gyroscope constant value drift that obtains by DD2 filtering, feedback compensation gyroscope output attitude and gyroscope output angle speed obtain the output of high-precision spacecraft attitude and high-precision gyro;
(6) behind the feedback compensation, repeat (1)-(5) step.
Principle of the present invention is: use gyroscope output angle speed to resolve attitude as the short-term reference, real-time attitude is provided, higher precision was arranged in the short time, but attitude information is constantly dispersed, the starlight vector that is obtained by star sensor provides high-precision long term reference spacecraft attitude by the Euler-q algorithm with certain frequency, with this long term reference attitude as observed quantity, by the DD2 wave filter, revise gyro output attitude, information is dispersed in correction, estimates gyroscope constant value drift simultaneously, the output of feedback compensation gyro, gyro is carried out on-line proving, keep high-accuracy posture to export in real time.
The present invention's advantage compared with prior art is: the invention provides the Euler q method of being calculated the attitude battle array by the revolving property at Euler axle and Euler angle, computing velocity is faster, adopting with attitude quaternion vector and gyro error simultaneously is the model of state, improved the precision of model, the singularity in having avoided calculating; And in the integrated attitude determination filtering, adopt DD2 Filtering Estimation hypercomplex number vector error and gyroscopic drift, the new wave filter that approaches with polynomial expression has been considered the uncertainty of new model state estimation, easier realization and do not need differentiate, precision is better than EKF and is not less than UKF, the computing velocity that when guaranteeing accuracy of attitude determination, improves greatly, estimate the constant value drift of gyro simultaneously, greatly facilitate the engineering practical application.
Description of drawings
Fig. 1 determines method flow diagram for the spacecraft attitude of use Euler-q algorithm of the present invention and DD2 filtering;
Fig. 2 is the relation between spacecraft body coordinate system and the geocentric inertial coordinate system.
Embodiment
As shown in Figure 1, 2, specific implementation method of the present invention is as follows:
1, by the responsive spacecraft attitude of gyroscope, the output tri-axis angular rate, and use the angle increment algorithm to upgrade the attitude battle array, and obtain spacecraft at geocentric inertial coordinate system, promptly earth the earth's core is an initial point, the sensing zenith is z
i, sensing first point of Aries is x
i, y
iWith x
i, z
iBecome the real-time attitude under the right-handed helix: course angle , pitching angle theta, roll angle γ.Gyro output attitude frequency is determined that by actual gyro sampling rate adopting output frequency in this example is 100HZ;
Calculate initial attitude hypercomplex number q battle array [q by initial attitude
1q
2q
3q
4], subscript 0 is an initial value:
The renewal matrix is:
Δ θ=[Δ θ
1Δ θ
2Δ θ
3] be the gyro output angle increment.Δ θ is the skew symmetry battle array of Δ θ:
Calculating attitude cosine battle array is:
Calculate real-time attitude angle under the geocentric inertial coordinate system by the direction cosine battle array.
The attitude angle solution formula is as follows:
The angle of pitch is:
θ=sin
-1(C
23) (8)
Crab angle such as following table:
C 21Value is judged | C 22Value is judged | Crab angle |
=0 | >0 | 0 |
<0 | >0 | atan -1(-C 21/C 22) |
<0 | =0 | π/2 |
<0 | <0 | atan -1(-C 21/C 22)+π |
=0 | <0 | π |
>0 | <0 | atan -1(-C 21/C 22)+π |
>0 | =0 | 3π/2 |
>0 | >0 | atan -1(-C 21/C 22)+2π |
Roll angle value such as following table:
C 13Value is judged | C 33Value is judged | Roll angle |
=0 | <0 | -π |
>0 | <0 | atan -1(-C 13/C 33)-π |
>0 | =0 | -π/2 |
Arbitrary value | >0 | atan -1(-C 13/C 33) |
<0 | =0 | π/2 |
<0 | <0 | atan -1(-C 13/C 33)+π |
2, utilize the responsive starlight of star sensor, pre-service obtains observing the coordinate of starlight vector s under the spacecraft body coordinate system through star chart, and promptly the spacecraft barycenter is an initial point, is reversed the zb axle with the earth's core line, with the tangent sensing of sky latitude east be x
b, y
bWith z
b, x
bBecome right-handed helix; And through obtaining the coordinate of corresponding reference starlight vector v under geocentric inertial coordinate system with it after the identification of star chart coupling;
3, use the Euler-q algorithm, utilize observation starlight vector and corresponding reference starlight vector with it, acquisition spacecraft attitude hypercomplex number;
Star sensor phase closing precision:
Calculate the B matrix:
Compute vector Z:Z={b
23-b
32b
31-b
13b
12-b
21}
T(11)
Calculate associated weight:
Wherein subscript i is an i starlight, β
iPrecision when observing i starlight for star sensor, n is the starlight number, is as the criterion with the number of importance in star map recognition starlight.
Introduce vector of unit length:
Calculate symmetric matrix:
The Euler axle is the proper vector of the minimal eigenvalue of H battle array.The computing formula of the proper vector of the minimal eigenvalue of H battle array is as follows:
Can directly from H battle array 3 rank feature equatioies, directly resolve out:
λ
3+aλ
2+bλ+c=0 (16)
Here a, b, the c coefficient can usually be represented according to the unit of H battle array:
a=-tr[H]=-h
11-h
22-h
33
Setup parameter is:
p
2=(a/3)
2-(b/3)
q=[(b/3)-(a/3)
2](a/3)-(c/3) (18)
Three real roots separate into:
λ
1=-p(sinw+cosw)-a/3
λ
2=q(sinw-cosw)-a/3 (19)
λ
3=2pcosw-a/3
Because 0≤w≤π/3 draw λ
iSatisfy condition:
0≤λ
1≤λ
2≤λ
3 (20)
λ
MinBe minimum λ value.Calculate λ once
MinProper vector, i.e. Euler q axle:
Provide three selectable separating:
All
(k=1,2,3) all are parallel, and the mould that calculates each vector is:
Select p near unit value
k(k=1,2,3), then by
Draw optimum Euler-q axle.
Calculate relevant Euler-q angle Φ
Opt:
Wherein:
The spacecraft attitude matrix A can be represented by Euler-q axle and Euler-q angle:
4, the spacecraft attitude that resolves of spacecraft attitude hypercomplex number that obtains with starlight observation and gyroscope is poor, as observed quantity, DD2 filtering is carried out in gyroscopic drift to the spacecraft attitude sum of errors, and the filtering cycle is by importance in star map recognition decision interval time, this routine medium frequency is 1HZ, and quantity of state is: (δ q
1δ q
2δ q
3b
1b
2b
3)
State equation:
Observation equation:
y=[I
3 0]X+w (27)
Can get after the discretize:
x
k+1=f(x
k,v
k) (28)
y
k=g(x
k,w
k) (29)
Below various in, subscript-for the amount average, subscript ^ for the amount valuation, subscript ~ for the amount error.
Introducing four cholesky factors decomposes:
Wherein,
Q is a state-noise variance battle array, and R is the measuring noise square difference battle array, Y
K-1=[y
0y
1Y
K-1], be the measurement matrix in past.
Use S
X, jExpression S
xThe j column element of matrix.
Can obtain following each amount:
Suppose that evaluated error is Gauss and does not have partially that establishing gap length is h
2=3.
One step upgraded:
From compound matrice, obtain:
Gain matrix is:
K
k=P
xy(k)[S
y(k)S
y(k)
T]
-1 (38)
Covariance matrix:
5, the spacecraft attitude error that obtains by DD2 filtering, feedback compensation gyroscope output attitude obtains high-precision spacecraft attitude; By the gyroscope constant value drift that filtering obtains, feedback compensation gyroscope output angle speed, the gyro output after obtaining proofreading and correct:
[δ q in the state that draws by filtering
1δ q
2δ q
3], obtain error quaternion:
Calculate the attitude behind the feedback compensation thus.
In must doing well by filtering
The output of correction gyro:
In the formula, w
OutBe the output of the gyro before revising.
6, behind the feedback compensation, repeat the 1-5 step.
Claims (3)
1. the spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method, it is characterized in that may further comprise the steps:
(1) by the responsive spacecraft attitude of gyroscope, the output tri-axis angular rate obtains the real-time attitude of spacecraft under geocentric inertial coordinate system;
(2) utilize the responsive starlight of star sensor, pre-service obtains observing the coordinate of starlight vector s under the spacecraft body coordinate system through star chart; And through obtaining the coordinate of corresponding reference starlight vector v under geocentric inertial coordinate system with it after the identification of star chart coupling;
(3) use the Euler-q algorithm, utilize observation starlight vector and corresponding reference starlight vector with it, acquisition spacecraft attitude hypercomplex number;
(4) spacecraft attitude that resolves of spacecraft attitude hypercomplex number that obtains with starlight observation and gyroscope is poor, and as observed quantity, DD2 filtering is carried out in gyroscopic drift to the spacecraft attitude sum of errors;
(5) the spacecraft attitude sum of errors gyroscope constant value drift that obtains by DD2 filtering, feedback compensation gyroscope output attitude and gyroscope output angle speed obtain attitude and the gyro output of spacecraft body coordinate system with respect to geocentric inertial coordinate system;
(6) behind the feedback compensation, repeat (1)-(5) step.
2. the spacecraft attitude based on Euler-q algorithm and DD2 filtering according to claim 1 is determined method, it is characterized in that: in the described step (4), estimation to the gyroscopic drift of spacecraft attitude sum of errors realizes by the DD2 wave filter, and attitude determines that model is as follows:
Quantity of state is: (δ q
1δ q
2δ q
3b
1b
2b
3), δ q wherein
1δ q
2δ q
3Be the vector part of attitude error hypercomplex number, b
1b
2b
3It is respectively the constant value drift of three outputs of gyro;
State equation, observation equation is respectively:
y=[I
3 0]X+W (2)
Wherein, y is observed quantity, ε
1, ε
2Be white Gaussian noise, w is a gyro output angle speed, and W is a measurement noise, I
3Be 3 * 3 unit vectors.
3. the spacecraft attitude based on Euler-q algorithm and DD2 filtering according to claim 1 is determined method, it is characterized in that: in the described step (5), the feedback compensation of spacecraft attitude sum of errors gyroscopic drift is realized by following steps:
[δ q in the state that draws by filtering
1δ q
2δ q
3], obtain error quaternion:
In must doing well by filtering
The output of correction gyro:
In the formula, w
OutBe the output of the gyro before revising.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007103012797A CN100559125C (en) | 2007-05-25 | 2007-12-19 | A kind of spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710099610.1 | 2007-05-25 | ||
CN200710099610 | 2007-05-25 | ||
CNB2007103012797A CN100559125C (en) | 2007-05-25 | 2007-12-19 | A kind of spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101196398A true CN101196398A (en) | 2008-06-11 |
CN100559125C CN100559125C (en) | 2009-11-11 |
Family
ID=39546948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007103012797A Expired - Fee Related CN100559125C (en) | 2007-05-25 | 2007-12-19 | A kind of spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100559125C (en) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846510A (en) * | 2010-05-28 | 2010-09-29 | 北京航空航天大学 | High-precision satellite attitude determination method based on star sensor and gyroscope |
CN101893440A (en) * | 2010-05-19 | 2010-11-24 | 哈尔滨工业大学 | Celestial autonomous navigation method based on star sensors |
CN101943585A (en) * | 2010-07-02 | 2011-01-12 | 哈尔滨工程大学 | Calibration method based on CCD star sensor |
CN101957203A (en) * | 2010-06-07 | 2011-01-26 | 哈尔滨工业大学 | High-accuracy star tracking method of star sensor |
CN102252673A (en) * | 2011-06-03 | 2011-11-23 | 哈尔滨工业大学 | Correction method for on-track aberration of star sensor |
CN102495831A (en) * | 2011-11-17 | 2012-06-13 | 西北工业大学 | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight |
CN102506865A (en) * | 2011-11-17 | 2012-06-20 | 西北工业大学 | Four-ary number polynomial approximate output method during extreme aerobat flight based on angular velocity |
CN102937450A (en) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | Relative attitude determining method based on gyroscope metrical information |
CN103727937A (en) * | 2013-11-20 | 2014-04-16 | 中国人民解放军海军大连舰艇学院 | Star sensor based naval ship attitude determination method |
CN103954289A (en) * | 2014-05-20 | 2014-07-30 | 哈尔滨工业大学 | Method for determining agile motor gesture of optical imaging satellite |
CN104034334A (en) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | Single-star and double-star attitude determination method of small-field star sensor |
KR101491891B1 (en) | 2010-11-08 | 2015-02-11 | 알버트 이바노비치 베구노프 | Method for producing aluminium by metallothermic reduction of trichloride with magnesium and apparatus for carrying out said method |
CN104833359A (en) * | 2015-05-27 | 2015-08-12 | 北京航空航天大学 | Discrete hidden Markov model characteristic sequence model based star chart pattern recognition method |
CN104833358A (en) * | 2015-05-13 | 2015-08-12 | 上海交通大学 | Method for determining geometric linear gesture of star sensor based on Rodrigues parameters |
CN105180946A (en) * | 2015-09-02 | 2015-12-23 | 上海新跃仪表厂 | Wideband measurement-based satellite high-precision attitude determination method and system thereof |
CN105424047A (en) * | 2015-10-30 | 2016-03-23 | 上海新跃仪表厂 | Spacecraft attitude pointing error identification method based on road sign information |
CN105957058A (en) * | 2016-04-21 | 2016-09-21 | 华中科技大学 | Preprocessing method of star map |
CN108344410A (en) * | 2018-01-23 | 2018-07-31 | 东南大学 | A method of the raising star sensor output frequency based on gyroscope auxiliary |
WO2018214227A1 (en) * | 2017-05-22 | 2018-11-29 | 深圳市靖洲科技有限公司 | Unmanned vehicle real-time posture measurement method |
CN109489661A (en) * | 2018-11-02 | 2019-03-19 | 上海航天控制技术研究所 | Gyro constant value drift estimation method when a kind of satellite is initially entered the orbit |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101943584B (en) * | 2010-07-02 | 2012-05-09 | 哈尔滨工程大学 | Alignment method based on CCD (Charge Coupled Device) star sensor |
-
2007
- 2007-12-19 CN CNB2007103012797A patent/CN100559125C/en not_active Expired - Fee Related
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893440A (en) * | 2010-05-19 | 2010-11-24 | 哈尔滨工业大学 | Celestial autonomous navigation method based on star sensors |
CN101893440B (en) * | 2010-05-19 | 2011-12-14 | 哈尔滨工业大学 | Celestial autonomous navigation method based on star sensors |
CN101846510B (en) * | 2010-05-28 | 2013-03-27 | 北京航空航天大学 | High-precision satellite attitude determination method based on star sensor and gyroscope |
CN101846510A (en) * | 2010-05-28 | 2010-09-29 | 北京航空航天大学 | High-precision satellite attitude determination method based on star sensor and gyroscope |
CN101957203A (en) * | 2010-06-07 | 2011-01-26 | 哈尔滨工业大学 | High-accuracy star tracking method of star sensor |
CN101957203B (en) * | 2010-06-07 | 2012-05-23 | 哈尔滨工业大学 | High-accuracy star tracking method of star sensor |
CN101943585A (en) * | 2010-07-02 | 2011-01-12 | 哈尔滨工程大学 | Calibration method based on CCD star sensor |
CN101943585B (en) * | 2010-07-02 | 2013-07-31 | 哈尔滨工程大学 | Calibration method based on CCD star sensor |
KR101491891B1 (en) | 2010-11-08 | 2015-02-11 | 알버트 이바노비치 베구노프 | Method for producing aluminium by metallothermic reduction of trichloride with magnesium and apparatus for carrying out said method |
CN102252673A (en) * | 2011-06-03 | 2011-11-23 | 哈尔滨工业大学 | Correction method for on-track aberration of star sensor |
CN102252673B (en) * | 2011-06-03 | 2012-10-24 | 哈尔滨工业大学 | Correction method for on-track aberration of star sensor |
CN102506865A (en) * | 2011-11-17 | 2012-06-20 | 西北工业大学 | Four-ary number polynomial approximate output method during extreme aerobat flight based on angular velocity |
CN102495831A (en) * | 2011-11-17 | 2012-06-13 | 西北工业大学 | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight |
CN102495831B (en) * | 2011-11-17 | 2014-08-06 | 西北工业大学 | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight |
CN102937450A (en) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | Relative attitude determining method based on gyroscope metrical information |
CN102937450B (en) * | 2012-10-31 | 2015-11-25 | 北京控制工程研究所 | A kind of relative attitude defining method based on gyro to measure information |
CN103727937A (en) * | 2013-11-20 | 2014-04-16 | 中国人民解放军海军大连舰艇学院 | Star sensor based naval ship attitude determination method |
CN103954289B (en) * | 2014-05-20 | 2016-06-22 | 哈尔滨工业大学 | The quick motor-driven attitude determination method of a kind of Optical Imaging Satellite |
CN103954289A (en) * | 2014-05-20 | 2014-07-30 | 哈尔滨工业大学 | Method for determining agile motor gesture of optical imaging satellite |
CN104034334A (en) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | Single-star and double-star attitude determination method of small-field star sensor |
CN104034334B (en) * | 2014-06-05 | 2016-09-14 | 哈尔滨工程大学 | Single star of a kind of small field of view star sensor and double star method for determining posture |
CN104833358A (en) * | 2015-05-13 | 2015-08-12 | 上海交通大学 | Method for determining geometric linear gesture of star sensor based on Rodrigues parameters |
CN104833359A (en) * | 2015-05-27 | 2015-08-12 | 北京航空航天大学 | Discrete hidden Markov model characteristic sequence model based star chart pattern recognition method |
CN104833359B (en) * | 2015-05-27 | 2016-06-15 | 北京航空航天大学 | A kind of star chart mode identification method based on discrete Markov characteristic sequence model |
CN105180946B (en) * | 2015-09-02 | 2019-01-01 | 上海新跃仪表厂 | Satellite high-precision attitude determination method and system based on wideband measurement |
CN105180946A (en) * | 2015-09-02 | 2015-12-23 | 上海新跃仪表厂 | Wideband measurement-based satellite high-precision attitude determination method and system thereof |
CN105424047A (en) * | 2015-10-30 | 2016-03-23 | 上海新跃仪表厂 | Spacecraft attitude pointing error identification method based on road sign information |
CN105424047B (en) * | 2015-10-30 | 2018-10-30 | 上海新跃仪表厂 | Spacecraft attitude error in pointing discrimination method based on mark information |
CN105957058A (en) * | 2016-04-21 | 2016-09-21 | 华中科技大学 | Preprocessing method of star map |
CN105957058B (en) * | 2016-04-21 | 2019-01-04 | 华中科技大学 | A kind of preprocess method of star chart |
WO2018214227A1 (en) * | 2017-05-22 | 2018-11-29 | 深圳市靖洲科技有限公司 | Unmanned vehicle real-time posture measurement method |
CN108344410A (en) * | 2018-01-23 | 2018-07-31 | 东南大学 | A method of the raising star sensor output frequency based on gyroscope auxiliary |
CN108344410B (en) * | 2018-01-23 | 2022-02-11 | 东南大学 | Gyroscope-assisted method for improving output frequency of star sensor |
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 |
Also Published As
Publication number | Publication date |
---|---|
CN100559125C (en) | 2009-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100559125C (en) | A kind of spacecraft attitude based on Euler-q algorithm and DD2 filtering is determined method | |
CN101949703B (en) | Strapdown inertial/satellite combined navigation filtering method | |
CN109324330B (en) | USBL/SINS tight combination navigation positioning method based on mixed derivative-free extended Kalman filtering | |
CN100487618C (en) | High precision combined posture-determining method based on optimally genetic REQUEST and GUPF | |
Wu et al. | Velocity/position integration formula part I: Application to in-flight coarse alignment | |
CN106291645A (en) | Be suitable to the volume kalman filter method that higher-dimension GNSS/INS couples deeply | |
CN102692225B (en) | Attitude heading reference system for low-cost small unmanned aerial vehicle | |
CN102353378B (en) | Adaptive federal filtering method of vector-form information distribution coefficients | |
Stančić et al. | The integration of strap-down INS and GPS based on adaptive error damping | |
CN101706284B (en) | Method for increasing position precision of optical fiber gyro strap-down inertial navigation system used by ship | |
CN102853837B (en) | MIMU and GNSS information fusion method | |
CN103630137A (en) | Correction method used for attitude and course angles of navigation system | |
CN103364817B (en) | POS system double-strapdown calculating post-processing method based on R-T-S smoothness | |
CN102175260A (en) | Error correction method of autonomous navigation system | |
CN110567455B (en) | Tightly-combined navigation method for quadrature updating volume Kalman filtering | |
KR100443550B1 (en) | IMU-GPS Integrated System including error correction system, Method for reducing search space of integer ambiguity, Method for detecting Cycle slip, and position, velocity, attitude determination Method using the same | |
CN111623771B (en) | Polarization inertial navigation tight combination navigation method based on light intensity | |
CN102654406A (en) | Initial alignment method for moving bases based on combination of nonlinear prediction filtering and cubature Kalman filtering | |
CN113375659B (en) | Pulsar navigation method based on starlight angular distance measurement information | |
CN102116634A (en) | Autonomous dimensionality reduction navigation method for deep sky object (DSO) landing detector | |
CN102980580A (en) | No-gyro satellite gesture determination method based on tensor product multi-cell robust heavy hydrogen (H2) filtering | |
CN103454665A (en) | Method for measuring double-difference GPS/SINS integrated navigation attitude | |
CN110954102A (en) | Magnetometer-assisted inertial navigation system and method for robot positioning | |
Liu et al. | Modified sage-husa adaptive Kalman filter-based SINS/DVL integrated navigation system for AUV | |
CN105988129A (en) | Scalar-estimation-algorithm-based INS/GNSS combined navigation 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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091111 Termination date: 20211219 |