CN103808330B - A kind of star sensor field trial satellite ephemeris computational methods - Google Patents
A kind of star sensor field trial satellite ephemeris computational methods Download PDFInfo
- Publication number
- CN103808330B CN103808330B CN201210439239.XA CN201210439239A CN103808330B CN 103808330 B CN103808330 B CN 103808330B CN 201210439239 A CN201210439239 A CN 201210439239A CN 103808330 B CN103808330 B CN 103808330B
- Authority
- CN
- China
- Prior art keywords
- test
- moment
- angle
- satellite
- calculate
- 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
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
-
- 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/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Manufacturing & Machinery (AREA)
- Navigation (AREA)
Abstract
The invention discloses a kind of star sensor field trial satellite ephemeris computational methods, it includes the steps of determining that calculating input and output amount parameter, utilizes input quantity, calculates output according to formula such as spherical triangle, astronomical ephemeris computation.Compared with prior art, its advantage and providing the benefit that: solve satellite under ground static operating mode with the problem that flight track ephemeris is inconsistent in-orbit, it is achieved that the simulation of process of the test Satellite Live Flying track.The method empirical tests result is correct, respond well, except can be applicable to star sensor field trial, applies also for the system-level field trial of other optical measurement such as sun sensor.
Description
Technical field
The present invention relates to a kind of star sensor field trial technology, the especially track ephemeris computation method of satellite when carrying out star sensor system level field trial.
Background technology
Field trial currently for star sensor is mainly product-level, is connected with spaceborne computer by star sensor so that it is the system-level test in outfield accessing satellite attitude control system closed loop did not also carry out.But in order to fully verify that in star sensor installation polarity in satellite attitude control system, service behaviour and systems soft ware, star sensor data processes the correctness with attitude algorithm, it is necessary for carrying out the system-level test in star sensor outfield.
System-level test is exactly the field trial place at broad view, select fine night, mounting means on analog satellite, star sensor is arranged on turntable, and be connected with spaceborne computer, allowing star sensor be operated in and operating mode approximate in orbit, output inertial attitude data are to spaceborne computer, and spaceborne computer calculates the track profile of satellite according to inertial attitude data in conjunction with track almanac data.
Test there is the problem that, the ephemeris of satellite and orbit parameter and situation about truly flying in-orbit inconsistent at present, test location longitude and latitude is exactly test the substar longitude and latitude that moment satellite orbital position is corresponding, the known quantity determined that, the semi-major axis of track, eccentricity, orbit inclination angle takes satellite actual value in-orbit, circular orbit satellite is constant value by these several parameters, above-mentioned parameter can be arranged in spaceborne computer software as constant in advance, and track argument, right ascension of ascending node, during test, on the star of spaceborne computer, the parameter such as clock can not directly obtain, therefore situation in-orbit can not truly be simulated.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of star sensor field trial satellite ephemeris computational methods.It can meet star sensor field trial needs.
For solving above-mentioned technical problem, a kind of star sensor field trial satellite ephemeris computational methods of the present invention, are utilize known inclination of satellite orbit i, test location latitudeAnd longitude, calculate local correspondence time on-test in moment of ascending node of orbit and the track argument of spaceborne computer clock, satellite position in track segmental arc and the method reflecting orbital plane and the right ascension of ascending node of equatoriat plane relative position.
Described star sensor field trial satellite ephemeris computational methods, it comprises the steps:
1) calculating input and output amount parameter is determined
Determine test location longitude, latitude, inclination of satellite orbit i, ascending node local time tdAs input quantity, moment on-test t0And the track argument U that the moment is corresponding0, right ascension of ascending nodeFor output;
2) test moment track argument U is calculated0
Utilize orbit inclination angle i first to calculate orbital plane and equatoriat plane angle, utilize its angle, location latitude, calculate test moment track argument U according to spherical triangle formula0;
3) moment on-test t is calculated0
Utilize orbit inclination angle i first to calculate orbital plane and equatoriat plane angle, recycle its angle andAscending node and meridian plane angle is calculated, the t local time in of utilizing its angle and ascending node according to spherical triangle formulad, longitudeCalculate moment on-test t0;
4) right ascension of ascending node is calculated
Utilize the t calculated0Calculate Julian date, and then calculate t0Greenwich sidereal time in moment angle, utilize t0、td、Calculate right ascension of ascending node。
The method that the present invention adopts, compared with prior art, its advantage and providing the benefit that: the present invention proposes a kind of system-level test satellite ephemeris computation method in general star sensor outfield, solve satellite under ground static operating mode with flight track ephemeris is inconsistent in-orbit problem, it is achieved that the simulation of process of the test Satellite Live Flying track.The method empirical tests result is correct, respond well, except can be applicable to star sensor field trial, applies also for the system-level field trial of other optical measurement such as sun sensor.
Accompanying drawing explanation
Below with reference to drawings and Examples, the invention will be further described.
Satellite orbital position and attitude embodiment schematic diagram when Fig. 1 is for test.
Detailed description of the invention
As it is shown in figure 1, satellite orbital position and attitude schematic diagram during for test, in figure, symbol implication is as follows:
For the earth's core;
N, S represent earth north, the South Pole respectively;
C is ascending node;
Face, place is earth equatorial plane;
For test location warp;
For ascending node place warp;
Face, place is satellite orbit plane;
For orbital plane withIntersection point, namely tests location;
For the celestial body coordinate system in zero attitude situation.
Known input quantity is:
Orbit inclination angle;
B point longitude and latitude、;
The t local time in of ascending node of orbitd。
Output to be asked is:
Track argument U0;
Moment on-test t0;
Right ascension of ascending node。
Calculating process is as follows:
From trirectangular spherical triangle
,, then by spherical triangle formula
Can obtain
So obtaining track argument U0:
Again by spherical triangle formula
Can obtain
Longitude corresponding to Beijing time is
Then should be during Beijing during test
t0=td+
If date respectively year, month, day during test, first calculate Julian date
Greenwich sidereal time, angle was:
Right ascension of ascending nodeFor:
Spaceborne computer calculates the parametric results of gained in aforementioned manners and carries out Orbit simulation computing, operation result can keep consistent with satellite Live Flying track when the overtesting location, analog satellite situation in orbit more truly, for the operating condition of test that the checking test offer of star sensor is more true to nature.
Claims (1)
1. star sensor field trial satellite ephemeris computational methods, utilize known inclination of satellite orbit i, test location latitudeAnd longitudeIt is characterized in that: the calculating to local correspondence time on-test in moment of ascending node of orbit and spaceborne computer clock, the track argument of reflection satellite position in track segmental arc and reflection orbital plane with the right ascension of ascending node of equatoriat plane relative position, including: 1) determine calculating input and output amount parameter
Determine test location longitude, latitude, inclination of satellite orbit i, ascending node local time tdAs input quantity, moment on-test t0And the track argument U that the moment is corresponding0, right ascension of ascending nodeFor output;
2) test moment track argument U is calculated0
Utilize orbit inclination angle i first to calculate orbital plane and equatoriat plane angle, utilize its angle, location latitude, calculate test moment track argument U according to spherical triangle formula0;
3) moment on-test t is calculated0
Utilize orbit inclination angle i first to calculate orbital plane and equatoriat plane angle, recycle its angle andAscending node and meridian plane angle is calculated, the t local time in of utilizing its angle and ascending node according to spherical triangle formulad, longitudeCalculate moment on-test t0;
4) right ascension of ascending node is calculated
Utilize the t calculated0Calculate Julian date, and then calculate t0Greenwich sidereal time in moment angle, utilize t0、td、Calculate right ascension of ascending node。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210439239.XA CN103808330B (en) | 2012-11-07 | 2012-11-07 | A kind of star sensor field trial satellite ephemeris computational methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210439239.XA CN103808330B (en) | 2012-11-07 | 2012-11-07 | A kind of star sensor field trial satellite ephemeris computational methods |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103808330A CN103808330A (en) | 2014-05-21 |
CN103808330B true CN103808330B (en) | 2016-07-20 |
Family
ID=50705439
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210439239.XA Active CN103808330B (en) | 2012-11-07 | 2012-11-07 | A kind of star sensor field trial satellite ephemeris computational methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103808330B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965208A (en) * | 2015-04-01 | 2015-10-07 | 北京理工雷科电子信息技术有限公司 | Ephemeris for simulation of satellite navigation signals, and almanac obtaining method |
CN108710379B (en) * | 2018-06-14 | 2021-02-02 | 上海卫星工程研究所 | Method for calculating yaw guide angle in imaging of stationary satellite |
CN113532423A (en) * | 2021-07-22 | 2021-10-22 | 哈尔滨工业大学 | Astronomical platform autonomous positioning method based on tilt angle sensor |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101082494A (en) * | 2007-06-19 | 2007-12-05 | 北京航空航天大学 | Self boundary marking method based on forecast filtering and UPF spacecraft shading device |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100666160B1 (en) * | 2004-12-29 | 2007-01-09 | 한국항공우주연구원 | Satellite Magnetometer Bias Correction Method Using Orbit Geometry |
-
2012
- 2012-11-07 CN CN201210439239.XA patent/CN103808330B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101082494A (en) * | 2007-06-19 | 2007-12-05 | 北京航空航天大学 | Self boundary marking method based on forecast filtering and UPF spacecraft shading device |
Non-Patent Citations (2)
Title |
---|
基于GPS非差观测值的CHAMP卫星约化动力学定轨方法研究;周建;《中国优秀硕士学位论文全文数据库信息科技辑》;20120715(第7期);正文第2,10-12,38-39页 * |
基于星敏感器的卫星瞬时姿态计算方法;郑万波;《吉林大学学报》;20030228;第21卷(第1期);第27-30页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103808330A (en) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101216319B (en) | Low orbit satellite multi-sensor fault tolerance autonomous navigation method based on federal UKF algorithm | |
CN104833352B (en) | High-precision vision/inertia integrated navigation method under multimedium complex environment | |
CN102252673B (en) | Correction method for on-track aberration of star sensor | |
Ning et al. | Autonomous satellite navigation using starlight refraction angle measurements | |
CN103453906B (en) | The Forecasting Methodology of satellite orbit | |
CN103646127B (en) | Satellite orbit gesture visual three-dimensional display packing | |
WO2013070518A1 (en) | Solar timer using gps technology | |
CN103616028A (en) | Star light refraction satellite autonomous navigation method based on single star sensor | |
CN104048664A (en) | Autonomous orbit determination method of navigation satellite constellation | |
CN109032176A (en) | A kind of geostationary orbit based on differential algebra is determining and parameter determination method | |
CN101893712B (en) | Weight selection fitting method for precise orbit determination of geostationary satellite | |
CN102426025A (en) | Simulation analysis method for drift correction angle during remote sensing satellite attitude maneuver | |
CN113359510B (en) | Beidou satellite navigation system signal simulator data real-time simulation system | |
CN103808330B (en) | A kind of star sensor field trial satellite ephemeris computational methods | |
US9806528B2 (en) | Systems and methods for estimating net solar energy production for airborne photovoltaic systems | |
Ma | Very long baseline interferometry applied to polar motion, relativity and geodesy | |
Hugentobler et al. | Satellite orbits and attitude | |
CN103852812B (en) | A kind of forward simulation occultation method | |
CN108072387B (en) | Ground correction method and system for on-orbit deviation of low-precision sensor | |
CN102506894A (en) | Stationary base platform absolute space attitude reference establishing method based on precise timing | |
CN103743488B (en) | Infrared imaging simulation method for globe limb background characteristics of remote sensing satellite | |
US9182584B2 (en) | Method and system for evaluating stare-time by a pointing system | |
Celik et al. | A generic three-dimensional model for solar energy reflected from mirrors in circular orbits | |
Jenness et al. | Pal: a positional astronomy library | |
Leonard et al. | Liaison-supplemented navigation for geosynchronous and lunar l1 orbiters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |