CN105424047A - Spacecraft attitude pointing error identification method based on road sign information - Google Patents
Spacecraft attitude pointing error identification method based on road sign information Download PDFInfo
- Publication number
- CN105424047A CN105424047A CN201510729265.XA CN201510729265A CN105424047A CN 105424047 A CN105424047 A CN 105424047A CN 201510729265 A CN201510729265 A CN 201510729265A CN 105424047 A CN105424047 A CN 105424047A
- Authority
- CN
- China
- Prior art keywords
- sensor
- vector
- error
- attitude
- matrix
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
The invention provides a spacecraft attitude pointing error identification method based on road sign information. The method comprises the steps that 1, a mathematical statistic model for sensor vector measurement and installation error angles is built; 2, an installation error angle observation equation between a full attitude sensor and a single vector sensor is constructed; 3, an optimal estimation method is designed to estimate an installation error angle between the full attitude sensor and the single vector sensor within some interval; 4, data fitting is performed on error angle data within different intervals to obtain a basic time varying installation error angle equation. The pointing error identification method is simple and accurate, the attitude angle and the angular speed of a spacecraft do not need to be calculated, attitude data of the spacecraft do not need to be acquired in real time, insensitivity on long-term data loss is achieved, and the practicability is high.
Description
Technical field
The present invention relates to spacecraft attitude error determination technology in-orbit, particularly, relate to a kind of spacecraft based on mark information attitude error method of estimation in-orbit.
Background technology
Remote sensing satellite points to in-orbit determines precision index, is directly connected to target geographic position determination precision level.The size of spacecraft error in pointing can comprehensively weigh out spacecraft in-orbit attitude determine and gesture stability situation, along with the continuous lifting of space mission demand, need badly and accurately provide spacecraft error in pointing and effectively revise.
Current domestic correlative study work mainly concentrates on spacecraft body attitude determines, points to the research determined less for spacecraft.Experience intense impact effect in Spacecraft Launch process, in orbit, be subject to the impact of the external environments such as different illumination conditions, cause the attitude of spacecraft itself and load to be pointed to and there is larger difference.Therefore, for the correlated results that spacecraft body attitude is determined, spacecraft can not be solved and point to the problem determined.In addition, on a small quantity for the result of study that spacecraft error in pointing is determined, need the attitude information utilizing spacecraft itself, need to carry out identification by the data of several orbital period, the requirement for data acquisition is higher.
The invention provides a kind of spacecraft attitude error in pointing discrimination method based on mark information, according to road sign point benchmark and observation information, spacecraft attitude information and auxiliary data thereof, obtain spacecraft real directional information in-orbit.Compared with the prior art, this discrimination method simple structure, identification result are accurate; The method, without the need to calculating attitude angle and the angular velocity of spacecraft, without the need to Real-time Obtaining spacecraft attitude data, is lost insensitive, practical to long term data.
Summary of the invention
For defect of the prior art, the object of this invention is to provide a kind of spacecraft attitude error in pointing discrimination method based on mark information.
According to the spacecraft attitude error in pointing discrimination method based on mark information provided by the invention, comprise the steps:
Step 1: the mathematical statistical model setting up sensor measurement and fix error angle;
Set up the mathematical statistical model of sensor vector measurement, alignment error, provide the estimated value of the initial value of single vector sensor observation vector under body coordinate system, full attitude sensor observation vector;
Step 2: build error in pointing angle observation equation according to the initial value of the observation vector of single vector sensor under body coordinate system, the estimated value of the full observation vector of attitude sensor under body coordinate system;
Step 3: according to the statistical property at error in pointing angle observation equation and error in pointing angle, statistical estimate is carried out to attitude error angle, obtains attitude error angular estimation value;
Step I: repeat step 1 to step 3, obtain different time sections attitude error angular estimation value in same track, carries out curve fitting to described different time sections attitude error angular estimation value and obtains long period attitude error fundamental equation.
Preferably, described step 1 comprises:
Step 1.1: the mathematical statistical model setting up sensor vector measurement, alignment error, specifically comprises following sub-step:
Step 1.1.1, sets up sensor reference measurement vector as follows:
In formula,
the direction vector measured value of sensor j,
direction of measurement true value,
be the measurement noises of sensor j, subscript k is the numbering of time, k=1,2 ..., N;
Step 1.1.2, determines that vector and installation relation is measured by satellite body system as follows:
In formula, S
jfor sensor j installs matrix,
for the direction vector under spacecraft body series,
for the direction vector true value under spacecraft body series,
for the direction vector noise under spacecraft body series;
Step S1.1.3, alignment error define method is as follows:
Wherein, S
jthe installation matrix of sensor j,
demarcate for ground and obtain matrix initial value being installed, M
jfor alignment error matrix, M
jdefinition as shown in formula (4):
I is unit matrix, and O () represents higher-order shear deformation, θ
jfor alignment error matrix corner vector,
for vectorial θ
jnegative symmetric matrix, as shown in formula (5).
Step 1.2: utilize the installation matrix initial value of the observed reading of a certain single vector sensor, this single vector sensor to obtain at body coordinate system observation vector initial value, computing formula is as follows:
In formula:
represent that single vector sensor j is at t
kthe initial value of the observation vector during moment under body coordinate system,
represent the initial installation matrix of single vector sensor j,
represent that single vector sensor j is at t
kthe measured value in moment,
represent the benchmark terrestrial reference vector of a certain single vector sensor j, A
krepresent attitude matrix, θ
jrepresent the fix error angle of a certain single vector sensor j,
represent that a certain single vector sensor j is at t
kobservation vector error during the moment under body coordinate system.
Step 1.3: according to a certain full attitude sensor at moment t
kthe base vector of the attitude of satellite obtained, a certain single vector sensor obtains the estimated value of this full attitude sensor at body coordinate system observation vector, and computing formula is as follows:
In formula:
represent that a certain full attitude sensor i is at moment t
kthe estimated value of the observation vector under body coordinate system that corresponding a certain single vector sensor j obtains, ξ
i,krepresent that a certain full attitude sensor i is at t
kmoment measuring error, θ
irepresent the fix error angle of a certain full attitude sensor i,
represent that a certain full attitude sensor i is at t
kthe body attitude matrix estimated value that moment obtains, comprises full attitude sensor i measuring error ξ
i,kwith fix error angle θ
ithe evaluated error introduced.
Preferably, described step 2 comprises:
Build error in pointing angle observation equation according to the estimated value of the initial value of the observation vector of a certain single vector sensor under body coordinate system, the observation vector of a certain full attitude sensor under body coordinate system, the equation of structure is as follows:
In formula: (θ
j-θ
i) be the error angle between a certain full attitude sensor i and a certain single vector sensor j, namely need the error in pointing solved,
represent and vector
vertical matrix, △ z
ij, kbetween a certain full attitude sensor i and a certain single vector sensor j in error angle observation equation at moment t
kobservation noise.
Preferably, described step 3 comprises:
According to error in pointing angle observation equation and statistical property, application least square method or Maximum Likelihood Estimation carry out statistical estimate to attitude error angle;
-application least square method obtains following accounting equation:
In formula
represent state variable to be estimated, H represents observing matrix, and Z represents observed quantity, ()
t()
-1representing matrix transposition, matrix inversion operation respectively;
-application maximum-likelihood method obtains following accounting equation:
In formula, R represents random vector △ z
ij, kcovariance matrix.
Preferably, described step I comprises:
According to benchmark mark information in different time sections in same track
namely on ground certain road sign point to the direction vector R of spacecraft centroid
ts, full attitude sensor i is at t
kthe body attitude matrix that moment obtains
single vector sensor j is at t
kthe measured value in moment
and the initial installation matrix of sensor j
repeat step 2 to step 3, obtain different time sections attitude error angle (θ in same track
j-θ
i) estimated value; The form of Fourier series is adopted to represent fundamental equation F (u) at attitude error angle in the whole orbital period, as shown in formula (11),
F(u)=x
0+x
1sin(u)+x
2cos(u)+x
3sin(u)+x
4cos(u)+...(11)
U is satellite ascending node argument.Utilize nonlinear least square fitting method, find the coefficient x meeting following equation,
Wherein ∑ is for adding and sign of operation, and min, for getting minimum operation symbol, obtains the Fourier coefficient (x in fundamental equation F (u)
0, x
1, x
2, x
3, x
4...) estimated value, x represents Fourier coefficient, and xdata represents the input data of fundamental equation F (u), xdata
irepresent i-th group of data of xdata.Ydata represents the observation data of fundamental equation F (u), ydata
irepresent i-th group of data of ydata.
Compared with prior art, the present invention has following beneficial effect:
1, the spacecraft attitude error in pointing discrimination method based on mark information provided by the invention, according to road sign point benchmark and observation information, spacecraft attitude information and auxiliary data thereof, obtain spacecraft real directional information in-orbit, compared with the prior art, this discrimination method simple structure, identification result are accurate.
2, the spacecraft attitude error in pointing discrimination method based on mark information provided by the invention, the method, without the need to calculating attitude angle and the angular velocity of spacecraft, without the need to Real-time Obtaining spacecraft attitude data, is lost insensitive, practical to long term data.
Accompanying drawing explanation
By reading the detailed description done non-limiting example with reference to the following drawings, other features, objects and advantages of the present invention will become more obvious:
Fig. 1 is the process flow diagram of the spacecraft attitude error in pointing discrimination method based on mark information provided by the invention.
Fig. 2 is road sign base vector graph of a relation provided by the invention.
Embodiment
Below in conjunction with concrete enforcement, the present invention is described in detail.Below implement to contribute to those skilled in the art and understand the present invention further, but do not limit the present invention in any form.It should be pointed out that to those skilled in the art, without departing from the inventive concept of the premise, some distortion and improvement can also be made.These all belong to protection scope of the present invention.
According to the spacecraft attitude error in pointing discrimination method based on mark information provided by the invention, provide concrete implementation step in practical application.
1, identification desired data and acquisition methods thereof
Star responsive relatively and between load the estimation of fix error angle need J2000.0 geocentric equatorial polar coordinate lower platform attitude data
with road sign base vector
road sign measurement vector under camera coordinates system
and matrix initial value installed by load camera
further, the observation vector estimated value under body series is converted it into
and observation vector initial value under body series
1.1 platform stance data
The attitude matrix A of platform relative orbit system can be obtained by image auxiliary data
boand track system relative inertness Conversion Matrix of Coordinate A
oi, through type (1) can obtain
1.2 road sign base vectors
Road sign base vector
for road sign point is to the direction vector R of spacecraft centroid
ts, R
etfor I system initial point O
ethe direction vector of direction disk point T, R
esfor I system initial point O
epoint to spacecraft centroid O
sdirection vector.Base vector
with vectorial R
et, R
esrelation such as formula (2) ~ (3), || || represent the computing of delivery value.
R
ts=R
es-R
et(2)
(1) GoogleEarth first, is utilized to find the longitude α of certain road sign point under the earth's core sphere is connected coordinate system (WGS84)
ewith latitude δ
e, utilize formula (4) that the direction vector R ' of road sign point under WGS84 coordinate system can be obtained
et.According to the temporal information in image auxiliary data, obtain WGS84 coordinate and be tied to J2000 Conversion Matrix of Coordinate A
wGS84ToJ2000, vector R can be obtained according to formula (5)
et, r in formula
efor the radius of the earth;
R
et=r
eA
WGS84ToJ2000R′
et(5)
(2) secondly, according to the position x of spacecraft in inertial space, y, z can obtain vector R
es, as formula (6);
R
es=[xyz]
T(6)
(3) by vector R
etand R
esroad sign base vector can be tried to achieve in substitution formula (2) ~ (3)
1.3 road sign measurement vectors
This coordinate (x under camera coordinates is obtained according to certain road sign point position in load diagram picture
c, y
c), combining camera focal distance f can obtain measurement vector
shown in (7),
1.4 full attitude sensor i are at t
kthe vector that moment attitude obtains
observation vector estimated value under body series
1.5 single vector sensor j are at t
kmoment observed reading
and initially matrix is installed
observation vector initial value under the body series obtained
2, constant error angular estimation in certain time interval
Image and auxiliary data thereof are converted to the data of above-mentioned needs, and long data are carried out segmentation.Recycle formula (10), (11) or (12) to obtain in different time interval corresponding constant error angle.
represent and vector
vertical matrix, △ z
ij, kbetween a certain full attitude sensor i and a certain single vector sensor j in error angle observation equation at moment t
kobservation noise.
represent state variable to be estimated, H represents observing matrix, and Z represents observed quantity, ()
t()
-1representing matrix transposition and matrix inversion operation respectively, in formula, R represents random vector △ z
ij, kcovariance matrix.
3, constant error angular estimation in certain time interval
Application of formula (13) ~ (14), become the coefficient (x of Fourier series in alignment error fundamental equation (13) when obtaining
0, x
1, x
2, x
3, x
4...).And then become fix error angle contour curve when can obtain.
F(u)=x
0+x
1sin(u)+x
2cos(u)+x
3sin(u)+x
4cos(u)+...(13)
Wherein, u is satellite ascending node argument, and ∑ is for adding and sign of operation, and min is for getting minimum operation symbol.
Above specific embodiments of the invention are described.It is to be appreciated that the present invention is not limited to above-mentioned particular implementation, those skilled in the art can make various distortion or amendment within the scope of the claims, and this does not affect flesh and blood of the present invention.
Claims (5)
1., based on a spacecraft attitude error in pointing discrimination method for mark information, it is characterized in that, comprise the steps:
Step 1: the mathematical statistical model setting up sensor measurement and fix error angle;
Set up the mathematical statistical model of sensor vector measurement, alignment error, provide the estimated value of the initial value of single vector sensor observation vector under body coordinate system, full attitude sensor observation vector;
Step 2: build error in pointing angle observation equation according to the initial value of the observation vector of single vector sensor under body coordinate system, the estimated value of the full observation vector of attitude sensor under body coordinate system;
Step 3: according to the statistical property at error in pointing angle observation equation and error in pointing angle, statistical estimate is carried out to attitude error angle, obtains attitude error angular estimation value;
Step I: repeat step 1 to step 3, obtain different time sections attitude error angular estimation value in same track, carries out curve fitting to described different time sections attitude error angular estimation value and obtains long period attitude error fundamental equation.
2. the spacecraft attitude error in pointing discrimination method based on mark information according to claim 1, it is characterized in that, described step 1 comprises:
Step 1.1: the mathematical statistical model setting up sensor vector measurement, alignment error, specifically comprises following sub-step:
Step 1.1.1, sets up sensor reference measurement vector as follows:
In formula,
the direction vector measured value of sensor j,
direction of measurement true value,
be the measurement noises of sensor j, subscript k is the numbering of time, k=1,2 ..., N;
Step 1.1.2, determines that vector and installation relation is measured by satellite body system as follows:
In formula, S
jfor sensor j installs matrix,
for the direction vector under spacecraft body series,
for the direction vector true value under spacecraft body series,
for the direction vector noise under spacecraft body series;
Step S1.1.3, alignment error define method is as follows:
Wherein, S
jthe installation matrix of sensor j,
demarcate for ground and obtain matrix initial value being installed, M
jfor alignment error matrix, M
jdefinition as shown in formula (4):
I is unit matrix, and O () represents higher-order shear deformation, θ
jfor alignment error matrix corner vector,
for vectorial θ
jnegative symmetric matrix, as shown in formula (5);
Step 1.2: utilize the installation matrix initial value of the observed reading of a certain single vector sensor, this single vector sensor to obtain at body coordinate system observation vector initial value, computing formula is as follows:
In formula:
represent that single vector sensor j is at t
kthe initial value of the observation vector during moment under body coordinate system,
represent the initial installation matrix of single vector sensor j,
represent that single vector sensor j is at t
kthe measured value in moment,
represent the benchmark terrestrial reference vector of a certain single vector sensor j, A
krepresent attitude matrix, θ
jrepresent the fix error angle of a certain single vector sensor j,
represent that a certain single vector sensor j is at t
kobservation vector error during the moment under body coordinate system;
Step 1.3: according to a certain full attitude sensor at moment t
kthe base vector of the attitude of satellite obtained, a certain single vector sensor obtains the estimated value of this full attitude sensor at body coordinate system observation vector, and computing formula is as follows:
In formula:
represent that a certain full attitude sensor i is at moment t
kthe estimated value of the observation vector under body coordinate system that corresponding a certain single vector sensor j obtains, ξ
i,krepresent that a certain full attitude sensor i is at t
kmoment measuring error, θ
irepresent the fix error angle of a certain full attitude sensor i,
represent that a certain full attitude sensor i is at t
kthe body attitude matrix estimated value that moment obtains, comprises full attitude sensor i measuring error ξ
i,kwith fix error angle θ
ithe evaluated error introduced.
3. the spacecraft attitude error in pointing discrimination method based on mark information according to claim 1, it is characterized in that, described step 2 comprises:
Build error in pointing angle observation equation according to the estimated value of the initial value of the observation vector of a certain single vector sensor under body coordinate system, the observation vector of a certain full attitude sensor under body coordinate system, the equation of structure is as follows:
In formula: (θ
j-θ
i) be the error angle between a certain full attitude sensor i and a certain single vector sensor j, namely need the error in pointing solved,
represent and vector
vertical matrix, Δ z
ij, kbetween a certain full attitude sensor i and a certain single vector sensor j in error angle observation equation at moment t
kobservation noise.
4. the spacecraft attitude error in pointing discrimination method based on mark information according to claim 1, it is characterized in that, described step 3 comprises:
According to error in pointing angle observation equation and statistical property, application least square method or Maximum Likelihood Estimation carry out statistical estimate to attitude error angle;
-application least square method obtains following accounting equation:
In formula
represent state variable to be estimated, H represents observing matrix, and Z represents observed quantity, ()
t, ()
-1representing matrix transposition, matrix inversion operation respectively;
-application maximum-likelihood method obtains following accounting equation:
In formula, R represents random vector Δ z
ij, kcovariance matrix.
5. the spacecraft attitude error in pointing discrimination method based on mark information according to claim 1, it is characterized in that, described step I comprises:
According to benchmark mark information in different time sections in same track
namely on ground certain road sign point to the direction vector R of spacecraft centroid
ts, full attitude sensor i is at t
kthe body attitude matrix that moment obtains
single vector sensor j is at t
kthe measured value in moment
and the initial installation matrix of sensor j
repeat step 2 to step 3, obtain different time sections attitude error angle (θ in same track
j-θ
i) estimated value; The form of Fourier series is adopted to represent fundamental equation F (u) at attitude error angle in the whole orbital period, as shown in formula (11),
F(u)=x
0+x
1sin(u)+x
2cos(u)+x
3sin(u)+x
4cos(u)+...(11)
U is satellite ascending node argument; Utilize nonlinear least square fitting method, find the coefficient x meeting following equation,
Wherein Σ is for adding and sign of operation, and min, for getting minimum operation symbol, obtains the Fourier coefficient (x in fundamental equation F (u)
0, x
1, x
2, x
3, x
4...) estimated value, x represents Fourier coefficient, and xdata represents the input data of fundamental equation F (u), xdata
irepresent i-th group of data of xdata.Ydata represents the observation data of fundamental equation F (u), ydata
irepresent i-th group of data of ydata.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510729265.XA CN105424047B (en) | 2015-10-30 | 2015-10-30 | Spacecraft attitude error in pointing discrimination method based on mark information |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510729265.XA CN105424047B (en) | 2015-10-30 | 2015-10-30 | Spacecraft attitude error in pointing discrimination method based on mark information |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105424047A true CN105424047A (en) | 2016-03-23 |
CN105424047B CN105424047B (en) | 2018-10-30 |
Family
ID=55502420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510729265.XA Active CN105424047B (en) | 2015-10-30 | 2015-10-30 | Spacecraft attitude error in pointing discrimination method based on mark information |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105424047B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973240A (en) * | 2016-07-15 | 2016-09-28 | 哈尔滨工大服务机器人有限公司 | Conversion method of navigation module coordinate system and robot coordinate system |
CN109459042A (en) * | 2018-12-07 | 2019-03-12 | 上海航天控制技术研究所 | A kind of spacecraft multi-mode autonomous navigation system based on world image |
CN114329943A (en) * | 2021-12-23 | 2022-04-12 | 哈尔滨工业大学(深圳) | Control performance boundary design method, device and medium based on attitude rotation matrix |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0683098A1 (en) * | 1994-05-16 | 1995-11-22 | Hughes Aircraft Company | Spacecraft attitude determination using sun sensor, earth sensor and space-to-ground link |
US6142423A (en) * | 1999-06-29 | 2000-11-07 | Trw Inc. | Ephemeris/attitude reference determination using on-board optics and other satellite ephemeris |
US6336062B1 (en) * | 1999-12-10 | 2002-01-01 | Nec Corporation | Attitude angle sensor correcting apparatus for an artificial satellite |
CN101196398A (en) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | Spacecraft posture confirming method based on Euler-q algorithm and DD2 filtering |
CN102435763A (en) * | 2011-09-16 | 2012-05-02 | 中国人民解放军国防科学技术大学 | Measuring method for attitude angular velocity of spacecraft based on star sensor |
CN102865866A (en) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | Satellite attitude determination method and attitude determination error analytical method based on two star sensors |
GB2496042A (en) * | 2011-10-25 | 2013-05-01 | Boeing Co | Spacecraft attitude and position determination system |
CN103674031A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Method for measuring attitude of spacecraft by using pulsar radiation vector and linear polarization information |
CN104165642A (en) * | 2014-08-29 | 2014-11-26 | 东南大学 | Method for directly correcting and compensating course angle of navigation system |
CN104296752A (en) * | 2014-09-24 | 2015-01-21 | 上海卫星工程研究所 | Autonomous spacecraft navigation system with combination of astronomical angle measurement and speed measurement, and navigation method of autonomous spacecraft navigation system |
CN104792340A (en) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | Star sensor installation error matrix and navigation system star-earth combined calibration and correction method |
-
2015
- 2015-10-30 CN CN201510729265.XA patent/CN105424047B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0683098A1 (en) * | 1994-05-16 | 1995-11-22 | Hughes Aircraft Company | Spacecraft attitude determination using sun sensor, earth sensor and space-to-ground link |
US6142423A (en) * | 1999-06-29 | 2000-11-07 | Trw Inc. | Ephemeris/attitude reference determination using on-board optics and other satellite ephemeris |
US6336062B1 (en) * | 1999-12-10 | 2002-01-01 | Nec Corporation | Attitude angle sensor correcting apparatus for an artificial satellite |
CN101196398A (en) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | Spacecraft posture confirming method based on Euler-q algorithm and DD2 filtering |
CN102435763A (en) * | 2011-09-16 | 2012-05-02 | 中国人民解放军国防科学技术大学 | Measuring method for attitude angular velocity of spacecraft based on star sensor |
GB2496042A (en) * | 2011-10-25 | 2013-05-01 | Boeing Co | Spacecraft attitude and position determination system |
CN103674031A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Method for measuring attitude of spacecraft by using pulsar radiation vector and linear polarization information |
CN102865866A (en) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | Satellite attitude determination method and attitude determination error analytical method based on two star sensors |
CN104165642A (en) * | 2014-08-29 | 2014-11-26 | 东南大学 | Method for directly correcting and compensating course angle of navigation system |
CN104296752A (en) * | 2014-09-24 | 2015-01-21 | 上海卫星工程研究所 | Autonomous spacecraft navigation system with combination of astronomical angle measurement and speed measurement, and navigation method of autonomous spacecraft navigation system |
CN104792340A (en) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | Star sensor installation error matrix and navigation system star-earth combined calibration and correction method |
Non-Patent Citations (1)
Title |
---|
马广富等: "基于反步法的主从航天器相对姿态控制", 《控制理论与应用》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973240A (en) * | 2016-07-15 | 2016-09-28 | 哈尔滨工大服务机器人有限公司 | Conversion method of navigation module coordinate system and robot coordinate system |
CN105973240B (en) * | 2016-07-15 | 2018-11-23 | 哈尔滨工大服务机器人有限公司 | A kind of conversion method of navigation module coordinate system and robot coordinate system |
CN109459042A (en) * | 2018-12-07 | 2019-03-12 | 上海航天控制技术研究所 | A kind of spacecraft multi-mode autonomous navigation system based on world image |
CN114329943A (en) * | 2021-12-23 | 2022-04-12 | 哈尔滨工业大学(深圳) | Control performance boundary design method, device and medium based on attitude rotation matrix |
CN114329943B (en) * | 2021-12-23 | 2023-01-24 | 哈尔滨工业大学(深圳) | Control performance boundary design method, device and medium based on attitude rotation matrix |
Also Published As
Publication number | Publication date |
---|---|
CN105424047B (en) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | A new direct filtering approach to INS/GNSS integration | |
CN103076017B (en) | Method for designing Mars entry phase autonomous navigation scheme based on observability degree analysis | |
CN105371870B (en) | A kind of in-orbit accuracy measurement method of star sensor based on star chart data | |
US8188919B2 (en) | Globally-convergent geo-location algorithm | |
CN102435763B (en) | Measuring method for attitude angular velocity of spacecraft based on star sensor | |
CN105606094A (en) | Information condition matched-filtering estimation method based on MEMS/GPS combined system | |
CN107132542A (en) | A kind of small feature loss soft landing autonomic air navigation aid based on optics and Doppler radar | |
Wang et al. | Technology of gravity aided inertial navigation system and its trial in South China Sea | |
CN103630139A (en) | Underwater vehicle all-attitude determination method based on magnetic gradient tensor measurement | |
CN103697883B (en) | A kind of aircraft horizontal attitude defining method based on skyline imaging | |
CN104266650B (en) | A kind of Mars landing device air approach section air navigation aid based on sampled point inheritance strategy | |
CN104165642A (en) | Method for directly correcting and compensating course angle of navigation system | |
CN105509750A (en) | Astronomical velocity measurement and ground radio combined Mars acquisition phase navigation method | |
CN102116626A (en) | Prediction and correction method of node of star point track image | |
CN104655135A (en) | Landmark-recognition-based aircraft visual navigation method | |
CN104034334B (en) | Single star of a kind of small field of view star sensor and double star method for determining posture | |
CN108364279A (en) | Determine the method that stationary orbit remote sensing satellite is directed toward deviation | |
Fang et al. | UKF for integrated vision and inertial sensors based on three-view geometry | |
CN105424047A (en) | Spacecraft attitude pointing error identification method based on road sign information | |
CN111428369A (en) | Method for calculating confidence of space target collision early warning result | |
CN104931071A (en) | Star sensor on-orbit geometric calibration method and system based on selecting weight iteration | |
CN105737850A (en) | Variable-scale unidirectional gravity sampling vector matching and positioning method based on particle filtering | |
CN103616024B (en) | A kind of planetary detection approach section autonomous navigation system observability degree is determined method | |
Chu et al. | Performance comparison of tight and loose INS-Camera integration | |
CN105389819A (en) | Robust semi-calibrating down-looking image epipolar rectification method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |