CN105737848B - System-level star sensor star viewing system and star viewing method - Google Patents
System-level star sensor star viewing system and star viewing method Download PDFInfo
- Publication number
- CN105737848B CN105737848B CN201410748618.6A CN201410748618A CN105737848B CN 105737848 B CN105737848 B CN 105737848B CN 201410748618 A CN201410748618 A CN 201410748618A CN 105737848 B CN105737848 B CN 105737848B
- Authority
- CN
- China
- Prior art keywords
- star sensor
- theodolite
- star
- axes
- coordinate system
- 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.)
- Active
Links
Images
Landscapes
- Navigation (AREA)
Abstract
The invention relates to a star sensor system-level star observation test technology, and provides a star sensor system-level outfield star observation test method for verifying the data correctness of a star sensor in real starry sky and the attitude determination precision of the star sensor driven by the rotation speed of the earth, which can test the correctness of the output data of the star sensor. Compared with the prior art, the method has the main advantages that the star sensor and the inertia measurement combination are installed on the rotary table together, the tool ensures that the three-axis directions of the inertia measurement combination and the star sensor are consistent, the Kalman filtering method is initially utilized in the test according to the data of the inertia measurement combination, the three axes of the star sensor are aligned with the geographic coordinate system under the condition of a static base, the star sensor installation matrix of a certain test time is determined and known on the basis of alignment, the correctness of the output data of the star sensor can be further judged, and the measurement precision of the star sensor under the driving of the earth rotating speed can be calculated by comparing coaxial gyro data.
Description
Technical Field
The invention relates to a star sensor system-level star observation and static base initial alignment test technology, in particular to a spacecraft ground physical simulation test method.
Background
The system-level outfield star observation test is an important test for verifying the performance of a star sensor and the function of an aircraft attitude determination module, and the existing star observation tests generally comprise two types: the first is that the star sensor independently views the star, and the method can simply and conveniently measure the polarity of the star sensor and the relative accuracy of the star sensor at the earth rotating speed. The second is that the star sensor and the aircraft control system view the star together, the method can conveniently examine the polarity of the star sensor and the control system together, but the three-axis direction of the star sensor is only approximately aligned with the geographic coordinate system at the initial test, and the absolute attitude determination precision of the star sensor is not easy to examine.
Based on the defects of the two field star viewing methods, the invention provides a novel field star viewing method, namely, a star sensor and an inertia measurement combination are simultaneously installed and tested on a rotary table, so that the three shafts of the star sensor and the inertia measurement combination are respectively parallel to each other for facilitating data processing, and initial alignment is carried out before an experiment under the condition of a static base, so that the system-level absolute measurement of the star sensor and the verification of the absolute precision under the drive of the earth rotating speed are realized.
Disclosure of Invention
The problem to be solved by the invention is that the absolute measurement precision of the star sensor cannot be measured in the prior art; in order to solve the problems, the invention provides a system-level star sensor star viewing system and a star viewing method.
The system-level star sensor star viewing system provided by the invention comprises: the star sensor and inertia measurement combination is arranged on the test turntable; under the condition of a static base, determining a three-axis and three-axis transformation matrix of a geographic coordinate system of the star sensor, and representing q as q by using a quaternionsl。
Further, still include: the theodolite comprises a first theodolite, a second theodolite, a third theodolite and a fourth theodolite; the first theodolite and the second theodolite are respectively in auto-collimation with the star sensor; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; and the third theodolite and the fourth theodolite are mutually aimed.
Further, the inertia measurement assembly includes a gyroscope to measure angular velocity and an accelerometer to measure acceleration.
The invention also provides a star viewing method of the system-level star sensor star viewing system, which comprises the following steps:
firstly, determining a conversion matrix of three axes of the inertial measurement combination and three axes of a geographic coordinate system under the condition of a static base;
secondly, determining a conversion matrix of the three axes of the inertia measurement combination and the three axes of the star sensor under the condition of a static base;
step three, determining a transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate system according to the step one and the step two;
fourthly, calculating to obtain an output quaternion of the star sensor under the geographic coordinate system according to a conversion matrix of the three axes of the star sensor and the three axes of the geographic coordinate system at a certain moment; and comparing with the actual output of the star sensor to obtain the measurement precision of the star sensor.
Further, under the condition that the gravity acceleration and the latitude of the test point are known, the conversion matrix of three axes of the inertial measurement combination and three axes of the geographic coordinate system is determined through a Kalman filtering algorithm by utilizing gyroscope data and accelerometer data in the inertial measurement combination.
Further, the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the first theodolite and the second theodolite are respectively in auto-collimation with the star sensor; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the third theodolite and the fourth theodolite are mutually aimed; and the first theodolite and the third theodolite are mutually aimed to obtain a conversion matrix of the three axes of the inertial measurement combination and the three axes of the star sensor.
Further, still include: obtaining the correctness of the attitude determination algorithm on the satellite, comprising the following steps:
step five, calculating a conversion quaternion q of the star sensor on the earth relative to a stationary orbitbs;
Sixthly, calculating the orbit quaternion q of the star sensoroi;
Step seven, calculating the quaternion of the attitude on the satelliteqliOutputting quaternion of the star sensor under a geographic coordinate system; q. q.sslIs the quaternion representation of the transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate systemAnd meanwhile, the on-satellite attitude determination algorithm is correct.
Further, when the three axes of the star sensor and the three axes of the geographic coordinate system convert the matrix into a unit matrix, the conversion quaternion of the star sensor on the earth relative to the stationary orbit is as follows:
The advantages of the invention include:
1) the measuring precision of the star sensor and the correctness of the on-satellite algorithm are accurately verified;
2) the absolute accuracy of the star sensor under the condition of the earth rotation speed can be obtained;
3) the star sensor data can be compared with the gyro data, and abnormal data caused by vibration interference in a star observation environment can be eliminated.
Drawings
FIG. 1 is a system level star sensor star system view; FIG. 2 is a schematic diagram of a star sensor mounting matrix determined by a theodolite stationing; FIG. 3 is a diagram of relative positions of three axes of the star sensor and a track coordinate system; FIG. 4 is a schematic diagram of the transformation between the three-axis coordinate system and the orbit coordinate system of the star sensor.
Detailed Description
The invention is further illustrated below with reference to the figures and examples.
The invention also provides a star viewing method of the system-level star sensor star viewing system, which comprises the following steps:
firstly, determining a conversion matrix of three axes of the inertial measurement combination and three axes of a geographic coordinate system under the condition of a static base;
the method for accurately calculating the transformation matrix between the inertial measurement unit and the geographic coordinate system by using kalman filtering based on the data of the gyroscope and the acceleration is well known to those skilled in the art and will not be described in detail herein.
In the embodiment, for convenience of data processing, the precise turntable is adjusted for multiple times by combining the conversion matrix between the inertial measurement combination and the geographic coordinate system for multiple times, so that the conversion matrix between the inertial measurement combination and the geographic coordinate system is similar to a unit matrix, and the nonzero element is within the allowable range of the test error.
Secondly, determining a conversion matrix of the three axes of the inertia measurement combination and the three axes of the star sensor under the condition of a static base;
three axes of the star sensor are respectively Xs、Ys、ZsRespectively pointing to south, sky and west, and an orbital coordinate system OO-XOYOZOThe relative positions are shown in figure 3.
The consistency of the star sensor and the inertia measurement combination measuring system can be ensured from hardware through finish machining, the star sensor is a high-precision attitude determination instrument, slight deviation of an installation matrix can have great influence on a test result, the installation matrix can be obtained through a transit station distribution method, the star sensor 03 and the inertia measurement combination 04 are installed on a test platform 06 by combining with a reference figure 1 and a reference figure 2, and the test platform 06 and the star sensor 03 are connected with a space GNC computer 02; the space power supply supplies power to the 01 star sensor 03, the inertia measurement assembly 04 and the space GNC computer 02, and the space GNC computer 02 is further connected with the test computer 05. After the star sensor 03 and the inertia measurement assembly 04 are installed on the test platform 06, a first theodolite T1 and a second theodolite cross sight T2 are erected; a third theodolite T3, a fourth theodolite T4; the first theodolite and the second theodolite are respectively in auto-collimation with the star sensor 03; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the third theodolite and the fourth theodolite are mutually aimed; and the first theodolite and the third theodolite are mutually aimed to obtain a conversion matrix of the three axes of the inertial measurement combination and the three axes of the star sensor.
The method for determining the mounting matrix of the star sensor and further determining the transformation matrix of the three axes of the inertia measurement combination and the star sensor through the transit station is well known to those skilled in the art and will not be described in detail herein.
Step three, determining a transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate system according to the step one and the step two;
fourthly, calculating to obtain an output quaternion of the star sensor under the geographic coordinate system according to a conversion matrix of the three axes of the star sensor and the three axes of the geographic coordinate system at a certain moment; and comparing with the actual output of the star sensor to obtain the measurement precision of the star sensor.
The star sensor outputs attitude quaternion under a J2000 inertial system, under the condition of a static base, the three-axis and geographic coordinate system conversion matrix of the star sensor is accurately determined through a bridge of an inertial measurement combination, the theoretical output of the star sensor under the geographic coordinate system at a certain moment can be obtained through calculation of a known formula according to the known calculation formula, and the theoretical output is compared with the actual output quaternion of the star sensor, so that the system-level measurement accuracy of the star sensor can be calculated.
The measurement accuracy of the star sensor can be calculated through the steps from the first step to the fourth step. When the star sensor is placed on the earth, the actual operation orbit of the star sensor is considered to be the same as the earth static orbit only in consideration of the aspect of attitude determination, orbit data of the star at the initial moment is observed, and the attitude angle on the star obtained after conversion of the orbit quaternion and the star sensor quaternion is 0 degree theoretically, so that the attitude determination algorithm on the star is verified. The step of verifying the correctness of the attitude determination algorithm on the satellite comprises the following steps:
step five, calculating a conversion quaternion q of the star sensor on the earth relative to a stationary orbitbs(ii) a When the three-axis of the star sensor and the three-axis conversion matrix of the geographic coordinate system are unit arrays, the conversion quaternion of the star sensor on the earth relative to a static orbit is as follows:
Sixthly, calculating the orbit quaternion q of the star sensoroi,
step seven, calculating the quaternion of the attitude on the satelliteqliOutputting quaternion of the star sensor under a geographic coordinate system; q. q.sslIs the quaternion representation of the transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate systemAnd meanwhile, the on-satellite attitude determination algorithm is correct.
In one embodiment of the present invention, the star sensor mounting quaternion noted above the star as shown in fig. 3 and 4 can be obtained through the following rotation sequence: the star sensor rotates around the axis Zs for 90+ Lat degrees in a negative way, Lat refers to the latitude of a star observation place, 31.1731 degrees is taken, and then the star sensor rotates around the axis Ys for 90 degrees in a positive way, so that the installation quaternion of the star sensor is obtained as follows:
obtaining an orbit quaternion at the initial satellite observation time from six static orbits, wherein if the initial satellite observation time is 20 o' clock 0 min 0 s in 4/20/2013, the ascent intersection point Ω is 132.40307 °, the orbit inclination angle i is 0 °, the latitude argument u is 120 °, and the initial orbit quaternion can be obtained:
qoi=[-0.5906 0 0 0.8069]Tthe initial quaternion of the star sensor is as follows:
although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (3)
1. A star viewing method of a system-level star sensor star viewing system is characterized in that the system-level star sensor star viewing system comprises the following steps: the star sensor and inertia measurement combination is arranged on the test turntable; under the condition of a static base, determining a three-axis and three-axis transformation matrix of a geographic coordinate system of the star sensor, and representing q as q by using a quaternionsl(ii) a Further comprising: the theodolite comprises a first theodolite, a second theodolite, a third theodolite and a fourth theodolite; the first theodolite and the second theodolite are respectively in auto-collimation with the star sensor; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the third theodolite and the fourth theodolite are mutually aimed; the inertia measurement combination comprises a gyroscope for measuring angular velocity and an accelerometer for measuring acceleration; the star viewing method comprises the following steps:
firstly, determining a conversion matrix of three axes of the inertial measurement combination and three axes of a geographic coordinate system under the condition of a static base; under the condition that the gravity acceleration and the latitude of the test point are known, determining a conversion matrix of three axes of the inertial measurement combination and three axes of a geographic coordinate system by using gyroscope data and accelerometer data in the inertial measurement combination through a Kalman filtering algorithm;
secondly, determining a conversion matrix of the three axes of the inertia measurement combination and the three axes of the star sensor under the condition of a static base; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the first theodolite and the second theodolite are respectively in auto-collimation with the star sensor; the first theodolite and the second theodolite are mutually aimed; the third theodolite and the fourth theodolite are respectively and automatically collimated with the inertial measurement combination; the third theodolite and the fourth theodolite are mutually aimed; the first theodolite and the third theodolite are mutually aimed to obtain a conversion matrix of the three axes of the inertial measurement combination and the three axes of the star sensor;
step three, determining a transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate system according to the step one and the step two;
fourthly, calculating to obtain an output quaternion of the star sensor under the geographic coordinate system according to a conversion matrix of the three axes of the star sensor and the three axes of the geographic coordinate system at a certain moment; comparing with the actual output of the star sensor to obtain the measurement precision of the star sensor;
the star viewing method further comprises the following steps: verifying the correctness of the attitude determination algorithm on the satellite, comprising the following steps:
step five, calculating a conversion quaternion q of the star sensor on the earth relative to a stationary orbitbs;
Sixthly, calculating the orbit quaternion q of the star sensoroi;
Step seven, calculating the quaternion of the attitude on the satelliteqliOutputting quaternion of the star sensor under a geographic coordinate system; q. q.sslIs the quaternion representation of the transformation matrix of the three axes of the star sensor and the three axes of the geographic coordinate systemAnd meanwhile, the on-satellite attitude determination algorithm is correct.
3. The system-level star sensor star viewing system observation method according to claim 2, wherein when the three axes of the star sensor and the three axes of the geographic coordinate system are transformed into a unit matrix, the transformed quaternion of the star sensor on the earth relative to the stationary orbit is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410748618.6A CN105737848B (en) | 2014-12-10 | 2014-12-10 | System-level star sensor star viewing system and star viewing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410748618.6A CN105737848B (en) | 2014-12-10 | 2014-12-10 | System-level star sensor star viewing system and star viewing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105737848A CN105737848A (en) | 2016-07-06 |
CN105737848B true CN105737848B (en) | 2021-03-09 |
Family
ID=56239545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410748618.6A Active CN105737848B (en) | 2014-12-10 | 2014-12-10 | System-level star sensor star viewing system and star viewing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105737848B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106908085B (en) * | 2017-04-10 | 2019-05-24 | 北京航天自动控制研究所 | A kind of starlight measurement direction control method based on earth rotation compensation |
CN108387245B (en) * | 2017-12-27 | 2020-09-22 | 北京航天时代光电科技有限公司 | Installation polarity measurement method for rate gyroscope and meter adding device for carrying |
CN109682395B (en) * | 2018-12-13 | 2021-02-02 | 上海航天控制技术研究所 | Star sensor dynamic noise equivalent angle evaluation method and system |
CN113607188B (en) * | 2021-08-02 | 2022-07-05 | 北京航空航天大学 | Theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method |
CN114001756B (en) * | 2021-11-03 | 2023-09-19 | 中国科学院光电技术研究所 | Small-field-of-view star sensor outfield ground star finding method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6131068A (en) * | 1999-08-30 | 2000-10-10 | Honeywell International Inc. | Accuracy of an inertial measurement unit |
CN102032918B (en) * | 2010-10-20 | 2012-11-07 | 郑州辰维科技股份有限公司 | Method for calibrating direction of three-probe start sensor |
CN102081360B (en) * | 2011-02-25 | 2012-12-05 | 哈尔滨工业大学 | Inertial astronomical combined navigation semi-physical experimentt system |
CN102288199B (en) * | 2011-06-22 | 2012-11-28 | 哈尔滨工业大学 | Ground testing method of star sensor |
CN103245364B (en) * | 2013-05-14 | 2015-06-03 | 哈尔滨工业大学 | Method for testing dynamic performance of star sensor |
-
2014
- 2014-12-10 CN CN201410748618.6A patent/CN105737848B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105737848A (en) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6613236B2 (en) | High-precision optical fiber gyro inertial surveying method | |
CN102252673B (en) | Correction method for on-track aberration of star sensor | |
CN105737848B (en) | System-level star sensor star viewing system and star viewing method | |
CN102927994B (en) | A kind of quick calibrating method of oblique redundant strapdown inertial navigation system | |
CN103575299B (en) | Utilize dual-axis rotation inertial navigation system alignment and the error correcting method of External Observation information | |
CN105371844B (en) | A kind of inertial navigation system initial method based on inertia/astronomical mutual assistance | |
CN102879014B (en) | Optical imaging autonomous navigation semi-physical simulation testing system for deep space exploration proximity process | |
CN104501835B (en) | The ground system test and method that a kind of space-oriented application heterogeneity IMU is initially aligned | |
CN104374388B (en) | Flight attitude determining method based on polarized light sensor | |
CN102288199B (en) | Ground testing method of star sensor | |
CN103090866B (en) | Method for restraining speed errors of single-shaft rotation optical fiber gyro strapdown inertial navigation system | |
CN103900608B (en) | A kind of low precision inertial alignment method based on quaternary number CKF | |
CN105160125B (en) | A kind of simulating analysis of star sensor quaternary number | |
CN106017507A (en) | Method for fast calibration of medium-and-low-precision optical fiber inertia units | |
CN103245364B (en) | Method for testing dynamic performance of star sensor | |
CN103759729B (en) | Adopt the soft lunar landing ground experiment initial attitude acquisition methods of inertial navigation | |
CN104296908A (en) | Three-degree-of-freedom air bearing table disturbance torque composition measuring device | |
CN103344872A (en) | Test method of installation polarity of star sensor | |
CN102707080A (en) | Method for simulating strapdown inertial navigation gyroscope by using star sensor | |
Xiaolin et al. | A tightly coupled rotational SINS/CNS integrated navigation method for aircraft | |
CN111879335A (en) | Calibration method for drift coefficient of multi-position gyroscope based on centrifugal machine | |
CN102519454B (en) | Selenocentric direction correction method for sun-earth-moon navigation | |
CN102183263A (en) | Method for calibrating fiber optic gyroscope constant drift | |
CN103913169A (en) | Strap-down inertial/starlight refraction combined navigation method of aircrafts | |
CN113465570B (en) | Air bearing table initial alignment method and system based on high-precision IMU |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |