CN107246875A - Relative attitude determines method between star under a kind of accurate formation task - Google Patents
Relative attitude determines method between star under a kind of accurate formation task Download PDFInfo
- Publication number
- CN107246875A CN107246875A CN201710532819.6A CN201710532819A CN107246875A CN 107246875 A CN107246875 A CN 107246875A CN 201710532819 A CN201710532819 A CN 201710532819A CN 107246875 A CN107246875 A CN 107246875A
- Authority
- CN
- China
- Prior art keywords
- star
- relative
- relative attitude
- primary
- attitude
- 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
-
- 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/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
Method is determined the invention discloses relative attitude between star under a kind of accurate formation task, is comprised the following steps:S1, transmit/receive antenna layout configuration design is completed based on measurement dimension constraint;S2, according to Cross-Link measurement information, sets up between star relative attitude relation and independent of relative attitude relation between the star of centroid of satellite;S3, the real time information of relative status between star is calculated by filtering algorithm.Relative attitude is determined between the present invention can realize the in-orbit real-time star of complex condition.
Description
Technical field
Method is determined the present invention is more particularly directed to relative attitude between star under a kind of accurate formation task.
Background technology
With continuing to develop for satellite technology, increasing country and tissue all increasingly tend to utilize multi-satellite
Satellite Formation Flying system is constituted to realize space tasks target.Therefore Satellite Formation Flying system control tasks become increasingly complex, while
It is faced with pressure in terms of high control precision and autonomy-oriented.
Relative attitude determines that technology is to meet posture Collaborative Control demand between accurate formation star between star, it is ensured that formation control essence
The support technology of degree.At present, relative attitude determines that method mainly has 3 kinds between the star proposed both at home and abroad:
1st kind of calculus of finite differences, i.e., realize that single star high-precision attitude is determined by the self-contained attitude sensor of satellite, utilize star
Between communication link carry out information exchange, pass through simple mathematical difference calculate obtain star between relative attitude;2nd species GPS methods, i.e.,
Cross-Link measurement is carried out using class GPS sensors, realizes that relative attitude is determined by filtering algorithm;3rd kind of direct method, when formation is defended
Relative distance between star relatively hundred meters of magnitudes when, relative attitude between star is determined come direct measurement using visual acuity device.
Above-mentioned 3 kinds of methods, calculus of finite differences is limited by single star sensor attitude determination accuracy in itself, while the processing of data difference
Mode meeting measurement by magnification error, attitude determination accuracy is limited between realizing star;Existing class GPS method is coupled with centroid of satellite,
Relative attitude measurement accuracy is limited by the in-orbit estimate accuracy of centroid of satellite between star, and the in-orbit of special centroid of satellite estimates non-reality
Posture determines requirement of real-time between when property influence star.Direct method is limited by interstellar distance and the visual field of visual acuity device is required,
Attitude determination accuracy and engineering practice are limited between star.
The content of the invention
Method is determined it is an object of the invention to provide relative attitude between star under a kind of accurate formation task, complex conditions are realized
Under between in-orbit real-time star relative attitude determine.
In order to realize the above object the present invention is achieved by the following technical solutions:
Relative attitude determines method between star under a kind of accurate formation task, is characterized in, comprises the following steps:
S1, transmit/receive antenna layout configuration design is completed based on measurement dimension constraint;
S2, according to Cross-Link measurement information, sets up between star relative attitude relation and independent of relative between the star of centroid of satellite
Posture relation;
S3, the real time information of relative status between star is calculated by filtering algorithm.
Transmitting antenna and reception antenna quantity are 2 in described step S1.
Described step S2 is specifically included:
S2.1, sets up independent of relative attitude state relation formula between the star of centroid of satellite:
ρij=R+Mr2j-r1i (1)
Wherein i=1,2, j=1, star known datum point is relative to the position vector of primary known datum point, ρ supplemented by 2, Rij
For the distance of j-th of reception antenna of i-th of transmitting antenna of primary to auxiliary star, r1iIt is the relative primary of i-th of transmitting antenna of primary
Position vector of the known datum point in primary body coordinate system, r2jSupplemented by star the relatively auxiliary star of j-th of reception antenna known to base
Position vector in auxiliary star body coordinate system on schedule, posture of the star body coordinate system relative to primary body coordinate system supplemented by M
Matrix;
S2.2, if ψ, θ,Supplemented by star with respect to relative attitude angle between the star of primary, then relative attitude state relation formula can be represented
For:
Wherein, ω1, ω2The respectively attitude angular velocity of primary and auxiliary star, celestial body coordinate is tied to primary body coordinate system supplemented by M
Transformation matrix of coordinates:
Wherein c represents that cos, s represent sin.
The real time information of relative status between star is calculated in described step S3 using EKF wave filters or UKF wave filters.
The present invention compared with prior art, with advantages below:
1st, engineering realizability is strong, attitude determination method between the high existing star at present of precision, or uses direct differential mode
Expand noise, either rely on centroid estimation precision or rely on visual field and interstellar distance, relative attitude determines precision and applies model
Limit is with, this method is laid out by optimizing transmit/receive antenna, sets up the measurement equation independent of centroid estimation, and constraint is few, no
Introduce attitude determination accuracy between extra error, star high.
2nd, relative attitude determines that overall process is clear between the star of method applicability and highly reliable Satellite Formation Flying, from correlation transmitting-receiving
The layout of antenna and datum mark and it is optimized to relative attitude and determines Real-Time Filtering algorithm, clearly, no visual field, barycenter is in-orbit estimates for meaning
Meter, interstellar distance etc. constrain, can in-orbit autonomous realization, ensure that the quality of accurate formation task run, be appearance between star
State Collaborative Control provides input and basis.
3rd, the strong this method of economy compares existing class GPS method, and transmitting antenna and reception antenna quantity are optimal, are not required to simultaneously
Will attitude measurement sensor between the star of such as optical camera, financial cost is low.
Brief description of the drawings
Fig. 1 determines the flow chart of method for relative attitude between star under a kind of accurate formation task of the present invention;
Fig. 2 is inventive antenna schematic layout pattern.
Embodiment
Below in conjunction with accompanying drawing, by describing a preferably specific embodiment in detail, the present invention is further elaborated.
Posture difference directly affects the precision of payload work between accurate formation task culminant star, with posture between clock star to compiling
There is coupling influence in team's control accuracy, therefore precision formation must enter posture determination between planet.For this engineer applied problem,
The present invention proposes under a kind of accurate formation task condition high-precision relative attitude between star and determines method, realizes that complex condition exists
Relative attitude is determined between the real-time star of rail.
As shown in figure 1, relative attitude determines method between star under a kind of accurate formation task, comprise the following steps:
S1, based on measurement dimension constraint complete transmit/receive antenna layout configuration design, design requirement be transmitting antenna and
The number of reception antenna is minimum;
S2, according to Cross-Link measurement information, sets up between star relative attitude relation and independent of relative between the star of centroid of satellite
Posture relation;
S3, the real time information of relative status between star is calculated by filtering algorithm.
Transmitting antenna and reception antenna quantity are 2 in above-mentioned step S1.
Above-mentioned step S2 is specifically included:
It is inventive antenna schematic layout pattern referring to Fig. 2, wherein triangle is transmitting antenna, circle is reception antenna, five
Angle star is the known datum point (being usually GNSS receiver antenna phase center) that satellite body is selected.
S2.1, sets up independent of relative attitude state relation formula between the star of centroid of satellite:
ρij=R+Mr2j-r1i (1)
Wherein i=1,2, j=1, star known datum point is relative to the position vector of primary known datum point, ρ supplemented by 2, Rij
For the distance of j-th of reception antenna of i-th of transmitting antenna of primary to auxiliary star, r1iIt is the relative primary of i-th of transmitting antenna of primary
Position vector of the known datum point in primary body coordinate system, r2jSupplemented by star the relatively auxiliary star of j-th of reception antenna known to base
Position vector in auxiliary star body coordinate system on schedule, posture of the star body coordinate system relative to primary body coordinate system supplemented by M
Matrix;
S2.2, if ψ, θ,Supplemented by star with respect to relative attitude angle between the star of primary, then relative attitude state relation formula can be represented
For:
Wherein, ω1, ω2The respectively attitude angular velocity of primary and auxiliary star, celestial body coordinate is tied to primary body coordinate system supplemented by M
Transformation matrix of coordinates:
Wherein c represents that cos, s represent sin.
Using EKF filtering (EKF) or UKF filtering (Unscented kalman filtering) meter in described step S3
The real time information of relative status between calculation star.
In summary, relative attitude determines method between star under a kind of accurate formation task of the present invention, realizes complex condition
Relative attitude is determined between in-orbit real-time star.
Although present disclosure is discussed in detail by above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read the above, for the present invention's
A variety of modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (4)
1. relative attitude determines method between star under a kind of accurate formation task, it is characterised in that comprise the following steps:
S1, transmit/receive antenna layout configuration design is completed based on measurement dimension constraint;
S2, according to Cross-Link measurement information, sets up between star relative attitude relation and independent of relative attitude between the star of centroid of satellite
Relation;
S3, the real time information of relative status between star is calculated by filtering algorithm.
2. relative attitude determines method between star under precision formation task as claimed in claim 1, it is characterised in that described step
Transmitting antenna and reception antenna quantity are 2 in rapid S1.
3. relative attitude determines method between star under precision formation task as claimed in claim 1, it is characterised in that described step
Rapid S2 is specifically included:
S2.1, sets up independent of relative attitude state relation formula between the star of centroid of satellite:
ρij=R+Mr2j-r1i (1)
Wherein i=1,2, j=1, star known datum point is relative to the position vector of primary known datum point, ρ supplemented by 2, RijBased on
I-th of transmitting antenna of star is to the distance of j-th of reception antenna of auxiliary star, r1iIt is known to the relative primary of i-th of transmitting antenna of primary
Position vector of the datum mark in primary body coordinate system, r2jSupplemented by star the relatively auxiliary star known datum point of j-th of reception antenna
Position vector in auxiliary star body coordinate system, attitude matrix of the star body coordinate system relative to primary body coordinate system supplemented by M;
S2.2, if ψ, θ,Supplemented by star with respect to relative attitude angle between the star of primary, then relative attitude state relation formula is represented by:
Wherein, ω1, ω2The respectively attitude angular velocity of primary and auxiliary star, celestial body coordinate is tied to the seat of primary body coordinate system supplemented by M
Mark transformation matrix:
Wherein c represents that cos, s represent sin.
4. relative attitude determines method between star under precision formation task as claimed in claim 1, it is characterised in that described step
The real time information of relative status between star is calculated in rapid S3 using EKF wave filters or UKF wave filters.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710532819.6A CN107246875B (en) | 2017-07-03 | 2017-07-03 | Method for determining inter-satellite relative attitude under precise formation task |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710532819.6A CN107246875B (en) | 2017-07-03 | 2017-07-03 | Method for determining inter-satellite relative attitude under precise formation task |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107246875A true CN107246875A (en) | 2017-10-13 |
CN107246875B CN107246875B (en) | 2020-07-31 |
Family
ID=60015021
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710532819.6A Active CN107246875B (en) | 2017-07-03 | 2017-07-03 | Method for determining inter-satellite relative attitude under precise formation task |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107246875B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110632935A (en) * | 2019-08-22 | 2019-12-31 | 上海航天控制技术研究所 | Autonomous control method for formation satellite around flight |
CN112208798A (en) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | Flight-around formation high code rate inter-satellite link switching method and system |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020004691A1 (en) * | 2000-03-10 | 2002-01-10 | Yasuhiro Kinashi | Attitude determination and alignment using electro-optical sensors and global navigation satellites |
US20060064244A1 (en) * | 1994-01-03 | 2006-03-23 | Robbins James E | Differential GPS corrections using virtual stations |
CN102262248A (en) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | Satellite gravity inversion method based on double-satellite spatial three-dimensional interpolation principle |
CN103546211A (en) * | 2013-10-31 | 2014-01-29 | 中国人民解放军国防科学技术大学 | Space division and time division intersatellite link rapid building method based on space-time prior link building information |
CN103852082A (en) * | 2012-11-30 | 2014-06-11 | 上海航天控制工程研究所 | Inter-satellite measurement and gyro attitude orbit integrated smoothing estimation method |
CN106154298A (en) * | 2016-08-12 | 2016-11-23 | 浙江大学 | The satellites formation GPS relative navigation system strengthened based on H_2O maser and method |
CN106443577A (en) * | 2016-09-05 | 2017-02-22 | 北京航空航天大学 | Multi-path error detection and elimination method in allusion to inter-satellite radio frequency relative measurement |
CN106679653A (en) * | 2016-12-12 | 2017-05-17 | 东南大学 | Relative measurement method of HEO (High Elliptical Orbit) satellite group based on satellite sensor and inter-satellite link |
-
2017
- 2017-07-03 CN CN201710532819.6A patent/CN107246875B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060064244A1 (en) * | 1994-01-03 | 2006-03-23 | Robbins James E | Differential GPS corrections using virtual stations |
US20020004691A1 (en) * | 2000-03-10 | 2002-01-10 | Yasuhiro Kinashi | Attitude determination and alignment using electro-optical sensors and global navigation satellites |
CN102262248A (en) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | Satellite gravity inversion method based on double-satellite spatial three-dimensional interpolation principle |
CN103852082A (en) * | 2012-11-30 | 2014-06-11 | 上海航天控制工程研究所 | Inter-satellite measurement and gyro attitude orbit integrated smoothing estimation method |
CN103546211A (en) * | 2013-10-31 | 2014-01-29 | 中国人民解放军国防科学技术大学 | Space division and time division intersatellite link rapid building method based on space-time prior link building information |
CN106154298A (en) * | 2016-08-12 | 2016-11-23 | 浙江大学 | The satellites formation GPS relative navigation system strengthened based on H_2O maser and method |
CN106443577A (en) * | 2016-09-05 | 2017-02-22 | 北京航空航天大学 | Multi-path error detection and elimination method in allusion to inter-satellite radio frequency relative measurement |
CN106679653A (en) * | 2016-12-12 | 2017-05-17 | 东南大学 | Relative measurement method of HEO (High Elliptical Orbit) satellite group based on satellite sensor and inter-satellite link |
Non-Patent Citations (1)
Title |
---|
陶久亮: ""卫星编队星间相对测量技术研究"", 《中国优秀硕士学位论文全文数据库-工程科技II辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110632935A (en) * | 2019-08-22 | 2019-12-31 | 上海航天控制技术研究所 | Autonomous control method for formation satellite around flight |
CN110632935B (en) * | 2019-08-22 | 2022-09-13 | 上海航天控制技术研究所 | Autonomous control method for formation satellite flying around |
CN112208798A (en) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | Flight-around formation high code rate inter-satellite link switching method and system |
CN112208798B (en) * | 2020-10-22 | 2022-03-22 | 上海卫星工程研究所 | Flight-around formation high code rate inter-satellite link switching method and system |
Also Published As
Publication number | Publication date |
---|---|
CN107246875B (en) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU775676B2 (en) | System for determining the heading and/or attitude of a body | |
US10763579B2 (en) | Mobile terminal antenna alignment using arbitrary orientation attitude | |
CN102506857B (en) | Relative attitude measurement real-time dynamic filter method based on dual-inertial measurement unit/differential global positioning system (IMU/DGPS) combination | |
CN105891863B (en) | It is a kind of based on highly constrained EKF localization method | |
CN100476360C (en) | Integrated navigation method based on star sensor calibration | |
CN106767793A (en) | A kind of AUV underwater navigation localization methods based on SINS/USBL tight integrations | |
CN104655152A (en) | Onboard distributed type POS real-time transmission alignment method based on federal filtering | |
CN103852082B (en) | Inter-satellite measurement and gyro attitude orbit integrated smoothing estimation method | |
CN101666868A (en) | Satellite signal vector tracking method based on SINS/GPS deep integration data fusion | |
CN106292677B (en) | Attitude control method and system based on sidereal hour angle | |
CN104422948A (en) | Embedded type combined navigation system and method thereof | |
CN105068102A (en) | DSP+FPGA-based ultra-tight combined navigation method | |
CN104567930A (en) | Transfer alignment method capable of estimating and compensating wing deflection deformation | |
CN103675861A (en) | Satellite autonomous orbit determination method based on satellite-borne GNSS multiple antennas | |
EP2585852A1 (en) | Moving platform ins range corrector (mpirc) | |
CN105547290B (en) | It is a kind of based on ultra short baseline locating system from latent device air navigation aid | |
CN101261130A (en) | On-board optical fibre SINS transferring and aligning accuracy evaluation method | |
CN109813306A (en) | A kind of unmanned vehicle planned trajectory satellite location data confidence level calculation method | |
CN103454665A (en) | Method for measuring double-difference GPS/SINS integrated navigation attitude | |
CN103279127A (en) | Autonomous control method for GEO (Geosynchronous) orbit satellite by using angle information only | |
CN106707322A (en) | RTK/SINS-based high dynamic positioning attitude-determining system and method | |
CN107246875A (en) | Relative attitude determines method between star under a kind of accurate formation task | |
CN102607558A (en) | Relative deformation measuring method based on inertial measurement units | |
CN109084760A (en) | Navigation system between a kind of building | |
Romanovas et al. | A method for IMU/GNSS/Doppler Velocity Log integration in marine applications |
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 |