CN110553667B - Method for carrying out thermal deformation compensation on star sensor - Google Patents
Method for carrying out thermal deformation compensation on star sensor Download PDFInfo
- Publication number
- CN110553667B CN110553667B CN201910857315.0A CN201910857315A CN110553667B CN 110553667 B CN110553667 B CN 110553667B CN 201910857315 A CN201910857315 A CN 201910857315A CN 110553667 B CN110553667 B CN 110553667B
- Authority
- CN
- China
- Prior art keywords
- star sensor
- deformation
- temperature
- mounting surface
- thermal deformation
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Navigation (AREA)
Abstract
The invention provides a method for carrying out thermal deformation compensation on a star sensor, which comprises the following steps: s1, establishing a first rectangular coordinate system on the star sensor mounting surface without thermal deformation; s2, controlling the temperature of the star sensor mounting surface, and determining the reference temperature at which the star sensor mounting surface is not subjected to thermal deformation; s3, taking the reference temperature as a starting point, carrying out temperature rise operation on the star sensor mounting surface, selecting a measuring point on the star sensor mounting surface, and measuring the deformation degree of the measuring point relative to the first rectangular coordinate system at different temperatures to obtain a plurality of groups of deformation degree measuring vectors; s4, fitting the multiple groups of deformation measurement vectors through a polynomial to obtain a temperature-deformation fitting formula of the deformation measurement vectors; s5, calculating a thermal deformation correction quaternion of the star sensor according to a temperature-deformation fitting formula; and S6, correcting the quaternion of the inertial coefficient measured by the star sensor according to the thermal deformation correction quaternion. The method is simple and reliable, and improves the measurement precision of the on-orbit attitude of the satellite.
Description
Technical Field
The invention relates to the field of spacecraft control, in particular to a method for compensating thermal deformation of a star sensor, which is used for compensating attitude measurement errors of the star sensor caused by thermal deformation of a star sensor mounting surface.
Background
With the improvement of satellite attitude control technology, the realization of high-precision attitude determination becomes a basic requirement for a satellite attitude determination system, and the accuracy of a star sensor which is taken as the most main measurement mechanism on the current satellite directly influences the accuracy of the whole attitude determination system.
The attitude determination system using the star sensor as a main measuring mechanism mainly has the following error sources: 1) the sensor has measurement errors such as noise, constant drift, various high and low frequency errors and the like; 2) the correction of the fusion algorithm may cause satellite attitude measurement errors, for example, the attitude determination accuracy finally obtained by selecting different filtering methods or fusion methods is different; 3) the star sensor and the star body have installation deviation, thermal deformation deviation and the like.
The prior art provides different solutions for the error sources, for example, the accuracy of fusion of the star sensor and other sensors is improved by optimizing a filtering algorithm, the measurement error is compensated by identifying the noise of the star sensor, the structural thermal deformation is reduced by a new material and new technology, and the like. The method for solving the measurement error of the star sensor in the prior art needs a large amount of theoretical research and engineering practice to obtain a certain effect, and has high calculation complexity and complex process.
Disclosure of Invention
The invention aims to provide a method for carrying out thermal deformation compensation on a star sensor, which is used for calibrating an inertia coefficient quaternion measured by the star sensor. According to the invention, a measuring point is selected on the star sensor mounting surface, and a temperature-deformation fitting formula of the measuring point is obtained through a thermal test. And calculating to obtain a thermal deformation correction quaternion according to the temperature-deformation fitting formula and an attitude conversion matrix from an actual installation coordinate system to a theoretical installation coordinate system of the star sensor after the star sensor is thermally deformed, monitoring the temperature of the measuring point in an on-orbit manner, and performing thermal deformation compensation on the quaternion of the inertial system measured by the star sensor through the thermal deformation correction quaternion to improve the satellite attitude determination precision.
In order to achieve the above object, the present invention provides a method for compensating thermal deformation of a star sensor, which is used for correcting an inertia system quaternion measured by the star sensor, and comprises the following steps:
s1, establishing a first rectangular coordinate system on the star sensor mounting surface without thermal deformation by taking the star sensor theoretical mounting coordinate system as a reference;
s2, performing a thermal test on the ground, performing temperature control on the star sensor mounting surface, and determining a reference temperature at which thermal deformation does not occur on the star sensor mounting surface; determining the reference temperature T of the star sensor mounting surface without thermal deformation by measuring whether the star sensor mounting surface deforms or not0；
S3, taking the reference temperature as a starting point, and carrying out temperature rise operation on the star sensor mounting surface; selecting a measuring point on a star sensor mounting surface, and monitoring the deformation degree of the measuring point relative to the first rectangular coordinate system under different heating temperatures to obtain a plurality of groups of deformation degree measuring vectors; the heating temperature is higher than the reference temperature T0；
S4, fitting the multiple groups of deformation measurement vectors through a polynomial to obtain a temperature-deformation fitting formula of the deformation measurement vectors;
s5, obtaining an attitude transformation matrix from the actual installation coordinate system of the thermally deformed star sensor to the theoretical installation coordinate system of the star sensor according to the temperature-deformation fitting formula; calculating and generating a thermal deformation correction quaternion according to the attitude transformation matrix;
s6, monitoring the temperature information of the measuring point when the satellite runs in orbit; when the temperature information exceeds the reference temperature T0And correcting the quaternion of the inertial system measured by the star sensor according to the thermal deformation correction quaternion.
In step S1, the first rectangular coordinate system is an orthogonal rectangular coordinate system and includes mutually perpendicular XrAxis, YrAxis and ZrAxis, said XrAxis, YrAxis, ZrThe axes intersecting at a point Or: wherein XrAxis, YrAxis falling on star sensor mounting surface without thermal deformation, ZrThe axis is perpendicular to the star sensor mounting surface without thermal deformation.
The distortion measurement vector in step S3 is (T)i,αi,βi) (ii) a Where i is the number of measurements, TiFor the temperature measured i, alphaiIs a temperature TiThe measurement points are relative to YrAngle of deformation of the shaft, betaiIs a temperature TiThe measurement points are relative to XrThe deformation angle of the shaft.
In step S4, the sets of deformation measurement vectors are specifically fitted by a first-order polynomial, and the temperature-deformation fitting formula is:
wherein T is the current temperature, and alpha and beta are the relative position of the measuring point to Y at the current temperature TrAxis, XrThe deformation angle of the shaft; kTα、KTβRespectively, measured point relative to YrAxis, XrTemperature coefficient of the shaft.
In step S5, the posture conversion matrix is:
converting A into quaternion formula according to direction cosine matrixs→rConverted to a thermal distortion corrected quaternion qs→r。
Inertial system quaternion q measured by the modified star sensor in step S6i→aComprises the following steps:
wherein q isicThe quaternion of the inertial system is obtained by the measurement of the star sensor; q. q.serrInstalling deviation for the star sensor; q. q.ss→bThe star sensor is a theoretical installation matrix from a theoretical installation coordinate system to the main system.
Compared with the prior art, the method carries out thermal deformation compensation on the star sensor by a method of ground test calibration and on-orbit temperature compensation. The deformation quantity of a measuring point of the star sensor mounting surface is measured for multiple times at different temperatures through a ground thermal test, so that the reference temperature of the star sensor mounting surface without deformation and a plurality of groups of deformation measurement vectors during deformation are obtained, and a degree-deformation formula of the deformation measurement vector temperature is obtained by utilizing a polynomial fitting method. When the satellite runs in orbit, the thermal deformation degree of the mounting surface of the star sensor can be calculated only according to the monitored temperature information of the measuring point, and then the quaternion of the inertia coefficient measured by the star sensor is compensated, so that a more accurate satellite attitude is obtained, and the satellite attitude determination precision is improved. According to the invention, the attitude measurement error caused by the thermal deformation of the mounting surface of the on-orbit can be effectively compensated by only adding one temperature measurement point to the mounting point of the star sensor, the compensation method is simple and reliable, and the determination precision of the on-orbit attitude can be effectively improved.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings used in the description will be briefly introduced, and it is obvious that the drawings in the following description are an embodiment of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts according to the drawings:
FIG. 1 is a flow chart of a method for compensating for thermal deformation of a star sensor according to the present invention;
FIG. 2 is a schematic diagram of a first rectangular coordinate system when the star sensor mounting surface is not thermally deformed according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating deformation angles of measurement points when thermal deformation occurs on a star sensor mounting surface according to an embodiment of the present invention;
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention provides a method for compensating thermal deformation of a star sensor, which is used for correcting an inertia system quaternion measured by the star sensor, and comprises the following steps:
s1, establishing a first rectangular coordinate system on the star sensor mounting surface without thermal deformation by taking the star sensor theoretical mounting coordinate system as a reference; as shown in FIG. 1, the first rectangular coordinate system comprises X's perpendicular to each otherrAxis, YrAxis and ZrAxis, said XrAxis, YrAxis, ZrThe axes intersecting at a point Or: wherein XrAxis, YrAxis falling on star sensor mounting surface without thermal deformation, ZrThe axis is perpendicular to the star sensor mounting surface without thermal deformation. Recording the first rectangular coordinate system as OrXrYrZr；
S2, performing a thermal test on the ground, heating the mounting surface of the star sensor, and determining the reference temperature T at which the mounting surface of the star sensor is not subjected to thermal deformation by measuring whether the mounting surface of the star sensor is deformed or not0；
S3, taking the reference temperature as a starting point, and carrying out temperature rise operation on the star sensor mounting surface; selecting a measuring point on a star sensor mounting surface, and measuring the deformation degree of the measuring point relative to the first rectangular coordinate system under different heating temperatures to obtain a plurality of groups of deformation degree measuring vectors; the deformation degree measurement vector is (T)i,αi,βi) (ii) a Where i is the number of measurements, TiFor the temperature measured i, alphaiIs a temperature TiThe measurement points are relative to YrAngle of deformation of the shaft, betaiIs a temperature TiThe measurement points are relative to XrThe deformation angle of the shaft. The heating temperature is higher than the reference temperature T0。
As shown in fig. 2, in the embodiment of the present application, the first rectangular coordinate system is established with the measurement point as the intersection point of the first rectangular coordinate system. When the heating temperature of the star sensor mounting surface exceeds the reference temperature T0Deformation occurs, temperature TiLower deformationThe section of the latter measurement point is denoted S. As shown in FIG. 3, the first rectangular coordinate system is around XrAxis of rotation alphaiAfter angling, rewinding YrShaft rotation betaiThe angle forms a second rectangular coordinate system OrXsYsZs. The second rectangular coordinate system is also an orthogonal rectangular coordinate system, including corresponding to X respectivelyrAxis, YrAxis, ZrX of axissAxis, YsAxis, ZsA shaft. XsAxis, YsPlane O formed by the shaftrXsYsFalls on the tangent plane S of the measuring point. Angle alphaiI.e. the temperature TiThe measurement points are relative to YrAngle of deformation of the shaft, betaiIs a temperature TiThe measurement points are relative to XrThe deformation angle of the shaft.
S4, fitting the multiple groups of deformation measurement vectors through a polynomial to obtain a temperature-deformation fitting formula of the deformation measurement vectors; in an application embodiment of the present invention, the sets of deformation measurement vectors are fitted by a first-order polynomial, and the temperature-deformation fitting formula is:
wherein T is the current temperature, and alpha and beta are the relative position of the measuring point to Y at the current temperature TrAxis, XrThe deformation angle of the shaft; kTα、KTβRespectively, measured point relative to YrAxis, XrTemperature coefficient of the shaft.
S5, calculating an attitude transformation matrix A from the actual installation coordinate system of the star sensor to the theoretical installation coordinate system of the star sensor after thermal deformation according to the temperature-deformation fitting formulas→rThis is prior art;
according to the direction cosine matrix to the quaternionAccording to formula As→rCalculating and generating thermal deformation correction quaternion qs→rThis is the prior art.
S6, monitoring the temperature information of the measuring point when the satellite operates in orbit, and when the temperature information exceeds the reference temperature T0And correcting the quaternion of the inertial system measured by the star sensor according to the thermal deformation correction quaternion.
Inertial system quaternion q measured by the modified star sensor in step S6i→aComprises the following steps:
wherein q isicThe quaternion of the inertial system is obtained by the measurement of the star sensor; q. q.serrInstalling deviation for the star sensor; q. q.ss→bThe star sensor is a theoretical installation matrix from a theoretical installation coordinate system to the main system. The system is a spacecraft self-fixed coordinate system, and a theoretical installation matrix q from a theoretical installation coordinate system to the system of the star sensor is obtaineds→bAre known in the art.
Compared with the prior art, the method carries out thermal deformation compensation on the star sensor by a method of ground test calibration and on-orbit temperature compensation. The deformation quantity of a measuring point of the star sensor mounting surface is measured for multiple times at different temperatures through a ground thermal test, so that the reference temperature of the star sensor mounting surface without deformation and a plurality of groups of deformation measurement vectors during deformation are obtained, and a degree-deformation formula of the deformation measurement vector temperature is obtained by utilizing a polynomial fitting method. When the satellite runs in orbit, the thermal deformation degree of the mounting surface of the star sensor can be calculated only according to the monitored temperature information of the measuring point, and then the quaternion of the inertia coefficient measured by the star sensor is compensated, so that a more accurate satellite attitude is obtained, and the satellite attitude determination precision is improved. According to the invention, the attitude measurement error caused by the thermal deformation of the mounting surface of the on-orbit can be effectively compensated by only adding one temperature measurement point to the mounting point of the star sensor, the compensation method is simple and reliable, and the determination precision of the on-orbit attitude can be effectively improved.
While the invention has been described with reference to specific embodiments, the invention is not limited thereto, and various equivalent modifications and substitutions can be easily made by those skilled in the art within the technical scope of the invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (5)
1. A method for carrying out thermal deformation compensation on a star sensor is used for correcting an inertia system quaternion measured by the star sensor and is characterized by comprising the following steps:
s1, establishing a first rectangular coordinate system on the star sensor mounting surface without thermal deformation by taking the star sensor theoretical mounting coordinate system as a reference;
s2, performing a thermal test on the ground, performing temperature control on the star sensor mounting surface, and determining a reference temperature at which thermal deformation does not occur on the star sensor mounting surface; determining the reference temperature T of the star sensor mounting surface without thermal deformation by measuring whether the star sensor mounting surface deforms or not0；
S3, taking the reference temperature as a starting point, and carrying out temperature rise operation on the star sensor mounting surface; selecting a measuring point on a star sensor mounting surface, and monitoring the deformation degree of the measuring point relative to the first rectangular coordinate system under different heating temperatures to obtain a plurality of groups of deformation degree measuring vectors; the heating temperature is higher than the reference temperature T0；
S4, fitting the multiple groups of deformation measurement vectors through a polynomial to obtain a temperature-deformation fitting formula of the deformation measurement vectors;
in step S4, the sets of deformation measurement vectors are fitted by a first-order polynomial, and the temperature-deformation fitting formula is:
wherein T is the current temperature, and alpha and beta are respectivelyThe measurement point is relative to Y at the current temperature TrAxis, XrThe deformation angle of the shaft; kTα、KTβRespectively, measured point relative to YrAxis, XrTemperature coefficient of the shaft;
s5, obtaining an attitude transformation matrix from the actual installation coordinate system of the thermally deformed star sensor to the theoretical installation coordinate system of the star sensor according to the temperature-deformation fitting formula; calculating and generating a thermal deformation correction quaternion according to the attitude transformation matrix;
s6, monitoring the temperature information of the measuring point when the satellite runs in orbit; when the temperature information exceeds the reference temperature T0And correcting the quaternion of the inertial system measured by the star sensor according to the thermal deformation correction quaternion.
2. The method for compensating for thermal deformation of a star sensor of claim 1, wherein the first rectangular coordinate system is an orthogonal rectangular coordinate system in step S1, and comprises mutually perpendicular X' SrAxis, YrAxis and ZrAxis, said XrAxis, YrAxis, ZrThe axes intersecting at a point Or: wherein XrAxis, YrAxis falling on star sensor mounting surface without thermal deformation, ZrThe axis is perpendicular to the star sensor mounting surface without thermal deformation.
3. The method for compensating for thermal deformation of a star sensor as claimed in claim 1, wherein the distortion measure vector in step S3 is (T)i,αi,βi) (ii) a Where i is the number of measurements, TiFor the temperature measured i, alphaiIs a temperature TiThe measurement points are relative to YrAngle of deformation of the shaft, betaiIs a temperature TiThe measurement points are relative to XrThe deformation angle of the shaft.
4. The method for thermally deforming and compensating the star sensor as claimed in claim 1, wherein the attitude transformation matrix in step S5 is:
converting A into quaternion formula according to direction cosine matrixs→rConverted to a thermal distortion corrected quaternion qs→r。
5. The method for compensating for thermal deformation of a star sensor as claimed in claim 4, wherein the modified quaternion q of the inertia coefficient measured by the star sensor in step S6i→aComprises the following steps:
wherein q isicThe quaternion of the inertial system is obtained by the measurement of the star sensor; q. q.serrInstalling deviation for the star sensor; q. q.ss→bThe star sensor is a theoretical installation matrix from a theoretical installation coordinate system to the main system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910857315.0A CN110553667B (en) | 2019-09-11 | 2019-09-11 | Method for carrying out thermal deformation compensation on star sensor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910857315.0A CN110553667B (en) | 2019-09-11 | 2019-09-11 | Method for carrying out thermal deformation compensation on star sensor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110553667A CN110553667A (en) | 2019-12-10 |
| CN110553667B true CN110553667B (en) | 2021-06-29 |
Family
ID=68739935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910857315.0A Active CN110553667B (en) | 2019-09-11 | 2019-09-11 | Method for carrying out thermal deformation compensation on star sensor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110553667B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111380567A (en) * | 2020-03-24 | 2020-07-07 | 上海卫星工程研究所 | Satellite star sensor pointing in-orbit thermal deformation correction system |
| CN114234962B (en) * | 2021-11-10 | 2023-09-12 | 上海航天控制技术研究所 | Multi-star sensor on-orbit thermal deformation correction method, storage medium and electronic equipment |
| CN114295092B (en) * | 2021-12-29 | 2023-08-01 | 航天科工智能运筹与信息安全研究院(武汉)有限公司 | Satellite radiometer thermal deformation error compensation method based on quaternion scanning imaging model |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04237038A (en) * | 1991-01-22 | 1992-08-25 | Nec Corp | Image pickup action correction system |
| JP2017024108A (en) * | 2015-07-21 | 2017-02-02 | ファナック株式会社 | Machine tool thermal displacement correction apparatus |
| CN106525078B (en) * | 2016-11-08 | 2019-03-15 | 上海航天控制技术研究所 | Rotating Platform for High Precision Star Sensor vacuum demarcates auxiliary device and thermal distortion correction method |
| CN106525001B (en) * | 2016-11-16 | 2019-05-03 | 上海卫星工程研究所 | It is directed toward calculation method in geostationary orbit remote sensing satellite camera optical axis space |
| CN109696179B (en) * | 2018-11-15 | 2022-10-18 | 上海航天控制技术研究所 | Thermoelastic error estimation method for remote sensing satellite star sensor |
-
2019
- 2019-09-11 CN CN201910857315.0A patent/CN110553667B/en active Active
Non-Patent Citations (1)
| Title |
|---|
| "相机与星敏组合体结构件常压热变形复合测量试验方法";张景川 等;《纳米技术与精密工程》;20140530;第12卷(第3期);正文第168-174页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110553667A (en) | 2019-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110553667B (en) | Method for carrying out thermal deformation compensation on star sensor | |
| CN111307293B (en) | Infrared human body temperature measurement system and temperature measurement method thereof | |
| CN112461261B (en) | Device and method for correcting temperature drift of MEMS gyroscope | |
| US20150310276A1 (en) | Method for the automatic correction of alignment errors in star tracker systems | |
| CN107246883B (en) | High-precision star sensor mounting matrix on-orbit real-time calibration method | |
| CN106767673B (en) | A kind of direction measurement method of satellite high-precision optical sensitive load | |
| CN112082574B (en) | Star sensor correction method and system | |
| CN107564057B (en) | High-orbit planar array optical satellite in-orbit geometric calibration method considering atmospheric refraction correction | |
| JP2008216062A (en) | Mobile-unit attitude measuring apparatus | |
| CN110411477A (en) | Based on the quick installation error on-orbit calibration method of the motor-driven star of sequence | |
| CN106679695B (en) | Focal plane thermal deformation testing method based on star sensor | |
| CN114754798A (en) | On-orbit identification and calibration method for gyro error characteristic parameters | |
| CN110674888B (en) | Head posture recognition method based on data fusion | |
| CN109655080B (en) | On-orbit calibration method for digital sun sensor | |
| CN110411481B (en) | Method and system for calibrating non-orthogonal error of gyroscope | |
| CN107228672B (en) | Star sensor and gyroscope data fusion method suitable for attitude maneuver working condition | |
| CN107228683B (en) | Slow-variation error real-time on-orbit correction method among multiple star sensors | |
| CN116256004A (en) | Star sensor on-orbit calibration and attitude resolving method based on improved particle swarm algorithm | |
| CN113252075B (en) | Star sensor precision correction method based on orbit periodic characteristics | |
| CN109916429A (en) | A kind of micro-electro-mechanical gyroscope calibrating installation, method and UAV Navigation System | |
| CN111220179B (en) | Inertial reference space-time accurate alignment method of optical navigation sensor | |
| CN113063444B (en) | Sub-angle second precision star sensor optical axis measurement reference deviation calibration method and system | |
| CN114838739B (en) | Satellite image geometric calibration method considering complete regression period | |
| CN106802225B (en) | Six-degree-of-freedom non-decoupling mechanism end pose error compensation method in time-varying temperature field | |
| CN110715795A (en) | Calibration and measurement method for fast reflector in photoelectric tracking system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 |