CN106843250A - A kind of limited plane of illumination determination methods in 01 formula sun sensor visual field - Google Patents

A kind of limited plane of illumination determination methods in 01 formula sun sensor visual field Download PDF

Info

Publication number
CN106843250A
CN106843250A CN201710081265.2A CN201710081265A CN106843250A CN 106843250 A CN106843250 A CN 106843250A CN 201710081265 A CN201710081265 A CN 201710081265A CN 106843250 A CN106843250 A CN 106843250A
Authority
CN
China
Prior art keywords
satellite
illumination
cos
sun
blocked
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
Application number
CN201710081265.2A
Other languages
Chinese (zh)
Other versions
CN106843250B (en
Inventor
张艳召
林荣峰
钟金凤
陈撼
裴甲瑞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN201710081265.2A priority Critical patent/CN106843250B/en
Publication of CN106843250A publication Critical patent/CN106843250A/en
Application granted granted Critical
Publication of CN106843250B publication Critical patent/CN106843250B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Navigation (AREA)

Abstract

The invention provides a kind of limited plane of illumination determination methods in 01 formula sun sensor visual field, wherein being projected in the position of satellite orbit face projection vector according to the sun ephemeris computation sun;Calculate the position interval of the ground shadow on track;Whether interval in illumination calculate current time satellite;Whether whether the 01 formula sun sensor that judgement is not blocked is shone, and judges there is illumination by other five quadrants blocked beyond satellite quadrant;Judgement is blocked whether satellite quadrant is shone.01 formula sensor and in-orbit recursion shade illumination ephemeris that i.e. the present invention can not blocked by judgement, joint judge solar azimuth, solve to block sun sensor because of load, the sensitive solar irradiation loss of learning problem for causing.

Description

A kind of limited plane of illumination determination methods in 01 formula sun sensor visual field
Technical field
The present invention relates to satellite guidance, navigation and control technology, a kind of 01 formula sun sensor visual field is related in particular to Limited plane of illumination determination methods.
Background technology
After attitude of satellite unstability, security control pattern can be entered, attitude control system need to be according to the information of 01 formula sensor Judge the orientation of sunshine, the sun can be directed at by attitude maneuver solar battery array, to ensure that whole energy source of star is supplied.According to this It is required that, the visual field of 01 formula sensor need to cover 6 quadrants of celestial body.At present, 01 formula sensor generally uses one neutralizing Whether certainly scheme, a 01 formula sun sensor can measure three quadrants by Lighting information, and it is sensitive that Fig. 1 show the 01 formula sun The configuration picture of device.
In recent years, the resolution ratio of load over the ground of earth observation satellite is improved constantly, and its physical dimension is continuously increased, can be to defending The 01 formula sun sensor of+Zb direction of principal axis of star causes to block.Fig. 2 show the whole star of 01 formula sun sensor and installs polarity letter Figure, it is known that satellite is furnished with two 01 formula sensors, ideally can just calculate sunny side according to two 01 formula sensors Position, but due to the interference of load over the ground ,+Zb direction of principal axis can not sensitivity go out whether solar irradiation information, so can not be according to biography System method is processed.
The content of the invention
In view of the shortcomings of the prior art, the present invention provides a kind of limited plane of illumination in 01 formula sun sensor visual field and judges Method, solves to block sun sensor because of load, the sensitive solar irradiation loss of learning problem for causing.
In order to achieve the above object, one technical scheme of the present invention is to provide a kind of limited light in 01 formula sun sensor visual field Show up determination methods, it includes procedure below:
Control software calculates the position that the sun is projected in satellite orbit face projection vector according to sun ephemeris on step one, star Put;
Control software calculates the position interval of the ground shadow on track on step 2, star;
Whether control software calculates current time satellite in illumination interval on step 3, star;
Control software judges whether the 01 formula sun sensor not blocked is shone on step 4, star;
Control software judges to be blocked whether satellite quadrant is shone on step 5, star.
In the present invention, whether satellite is in-orbit is in area of illumination, can be carried out according to satellite orbit parameter, sun ephemeris information The judgement in geometry domain, then the information according to other formula sensors of 5 quadrants 01 whether judge the quadrant that is blocked in light According to so as to realize the determination to solar azimuth.
The method that the present invention is used, compared with prior art, its advantage and beneficial effect are:The present invention can be in 01 sensor When being blocked by load, the orientation of the sun is judged, so as to smoothly realize to sun-orientation control, it is ensured that whole energy source of star safety.
Brief description of the drawings
Fig. 1 is the configuration picture of 01 formula sun sensor multiple directions.
Fig. 2 is that the whole star of 01 formula sun sensor installs polarity sketch.
Fig. 3 is the flow chart of the limited plane of illumination determination methods in 01 formula sun sensor visual field of the present invention.
Fig. 4 is the ground sketch that calculates of shadow in the present invention.
Specific embodiment
The present invention can judge what is blocked according to satellite orbit, sun ephemeris, 01 sensor not blocked combination Whether satellite quadrant is shone.
As shown in figure 3, the present invention provides a kind of 01 formula sun sensor visual field limited plane of illumination determination methods, it include with Lower step:
Control software is projected in the position of satellite orbit face projection vector according to the sun ephemeris computation sun on step one, star Put;
Control software calculates the position interval of the ground shadow on track on step 2, star;
Whether control software calculates current time satellite in illumination interval on step 3, star;
Control software calculates whether 01 sensor not blocked is shone on step 4, star;
Control software judges to be blocked whether satellite quadrant is shone on step 5, star.
Below with reference to drawings and Examples, each step to the inventive method is described further.
1) in step one, calculated comprising the projected angle to the sun in satellite orbit face:
Sun Dec, declination are calculated
δs=arcsin (sin Λ sin ε)
Wherein, Λ, ε represent solar celestial longitude and yellow red angle respectively.
The corresponding latitude argument of projection of the solar vector in orbital plane
Wherein, i, Ω represent the orbit inclination angle and right ascension of ascending node of satellite respectively.
2) in step 2, calculated comprising the ground shadow on track:
Coordinate of the satellite orbit face normal direction under J2000 can be expressed as
Ssat=-[sin Ω sini-cos Ω sini cosi]T (3)
According to Sun Dec αsWith declination δs, there is shown coordinate of the solar vector under inertial system J2000
Ssun=[cos αscosδs sinαscosδs sinδs]T (4)
In can be to calculate the angle between satellite orbit face normal negative sense and solar vector
The earth regards a circle, its radius length R aseRepresent, then the earth is one on the vertical face of solar vector Circle, the roundlet in corresponding diagram 4;The track of satellite is near-circular orbit, its on the vertical face of solar vector be projected as one it is ellipse Circle, the oval semi-major axis is satellite orbit semi-major axis a, semi-minor axis is a | cos β |.
Can be obtained by the relation in Fig. 4, the half-shade angle α of ground shadow is
3) whether in step 3, interval in illumination current time satellite is calculated:
So in t, if the latitude argument of satellite meets following relation, in area of illumination.
ut-u0∈[0,π-α]∪[π+α,2π] (7)
Wherein, utIt is the latitude argument of satellite.
4) in step 4, whether 01 sensor that judgement is not blocked is shone
AD sample rate current results according to 01 sensor not blocked, On board computer can judge to be blocked satellite as Whether other five quadrants beyond limit have the sun.
5) in step 5, judgement is blocked whether satellite quadrant is shone
According to the result of step 4, if in area of illumination, the 01 sensor no signal do not blocked, then presumption is blocked and defended Astrology limit has sun signal.
In sum, after satellite unstability, attitude control system need to be according to where 01 formula sun sensor judges sunshine The quadrant of satellite, the motor-driven guarantee windsurfing of the attitude of satellite is carried out to day state according to solar direction, it is ensured that whole energy source of star supply.Due to Satellite health external applied load physical size is larger, a certain quadrant can be caused to block 01 formula sun sensor visual field, and the present invention can With the 01 formula sensor and in-orbit recursion shade illumination ephemeris that are not blocked by judgement, joint judges solar azimuth.
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 of the invention Various modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.

Claims (5)

1. a kind of 01 formula sun sensor visual field is limited plane of illumination determination methods, it is characterised in that comprising procedure below:
Control software calculates the position that the sun is projected in satellite orbit face projection vector according to sun ephemeris on step one, star;
Control software calculates the position interval of the ground shadow on track on step 2, star;
Whether control software calculates current time satellite in illumination interval on step 3, star;
Control software judges whether the 01 formula sun sensor that is not blocked is shone on step 4, star, judge to be blocked satellite as Whether other five quadrants beyond limit have illumination;
Control software judges to be blocked whether satellite quadrant is shone on step 5, star:If do not blocked during in area of illumination 01 formula sun sensor no signal, then judge that being blocked satellite quadrant has sun signal.
2. 01 formula sun sensor visual field described in claim 1 is limited plane of illumination determination methods, it is characterised in that in step one In, ask for the corresponding latitude argument of projection of the solar vector on satellite orbit face
u 0 = a r c t a n ( - s i n ( Ω - α s ) cos i cosδ s + sin i sinδ s | cos ( Ω - α s ) cosδ s | )
Wherein, i, Ω represent the orbit inclination angle and right ascension of ascending node of satellite respectively;
Sun Dec αs, declination δsFor
&alpha; s = a r c t g ( tan &Lambda; c o s &epsiv; ) cos &Lambda; &GreaterEqual; 0 &pi; + a r c t g ( tan &Lambda; c o s &epsiv; ) cos &Lambda; < 0
δs=arcsin (sin Λ sin ε)
Wherein, Λ, ε represent solar celestial longitude and yellow red angle respectively.
3. 01 formula sun sensor visual field described in claim 2 is limited plane of illumination determination methods, it is characterised in that in step 2 In, the half-shade angle α for asking for ground shadow is
&alpha; = a r c s i n ( R e 2 - a 2 cos 2 &beta; a 2 sin 2 &beta; )
Wherein, ReIt is earth radius length;Satellite orbit projects formation ellipse, the half of the ellipse on the vertical face of solar vector Major axis is satellite orbit semi-major axis a, semi-minor axis is a | cos β |.
Angle between satellite orbit face normal negative sense and solar vector
β=arccos (Ssat·Ssun)
=arccos (- cos δssin i sinΔα-cos i sinδs)
Coordinate of the satellite orbit face normal direction under inertial system J2000
Ssat=-[sin Ω sin i-cos Ω sin i cos i]T
Coordinate of the solar vector under inertial system J2000
Ssun=[cos αscosδs sinαscosδs sinδs]T
4. 01 formula sun sensor visual field described in claim 3 is limited plane of illumination determination methods, it is characterised in that
In step 3, the latitude argument of current time t satellite meets following relation, judges that current time satellite is in area of illumination Between:
ut-u0∈[0,π-α]∪[π+α,2π]
Wherein, utIt is the latitude argument of satellite.
5. 01 formula sun sensor visual field described in claim 1 or 4 is limited plane of illumination determination methods, it is characterised in that
In step 4, On board computer judges to be hidden according to the AD sample rate current results of the 01 formula sun sensor not blocked Whether other five quadrants beyond gear satellite quadrant have illumination.
CN201710081265.2A 2017-02-15 2017-02-15 method for judging limited illumination surface of 01-type sun sensor view field Active CN106843250B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710081265.2A CN106843250B (en) 2017-02-15 2017-02-15 method for judging limited illumination surface of 01-type sun sensor view field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710081265.2A CN106843250B (en) 2017-02-15 2017-02-15 method for judging limited illumination surface of 01-type sun sensor view field

Publications (2)

Publication Number Publication Date
CN106843250A true CN106843250A (en) 2017-06-13
CN106843250B CN106843250B (en) 2019-12-10

Family

ID=59127377

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710081265.2A Active CN106843250B (en) 2017-02-15 2017-02-15 method for judging limited illumination surface of 01-type sun sensor view field

Country Status (1)

Country Link
CN (1) CN106843250B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101231170A (en) * 2008-01-31 2008-07-30 北京控制工程研究所 Method for processing information of APS sun sensor
CN101726299A (en) * 2009-12-09 2010-06-09 哈尔滨工业大学 Sun sensor simulator for satellite closed cycle simulation test
US20130179073A1 (en) * 2012-01-11 2013-07-11 Jena-Optronik Gmbh Method for increasing the reliability of sensor systems
CN104097793A (en) * 2014-06-24 2014-10-15 上海微小卫星工程中心 Zero momentum magnetic control sun capture device and method of satellite
CN105444760A (en) * 2015-11-30 2016-03-30 北京控制工程研究所 Sun sensor based on collection of photovoltaic micro energy
CN105737822A (en) * 2016-02-06 2016-07-06 上海微小卫星工程中心 All-day ball view field sun sensor and determination method for sun vector thereof
CN105806369A (en) * 2016-05-20 2016-07-27 上海航天控制技术研究所 Star sensor in-orbit aberration correction method
CN106289237A (en) * 2016-07-27 2017-01-04 上海航天控制技术研究所 APS digital sun sensor solar angle computational methods and system

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101231170A (en) * 2008-01-31 2008-07-30 北京控制工程研究所 Method for processing information of APS sun sensor
CN101726299A (en) * 2009-12-09 2010-06-09 哈尔滨工业大学 Sun sensor simulator for satellite closed cycle simulation test
US20130179073A1 (en) * 2012-01-11 2013-07-11 Jena-Optronik Gmbh Method for increasing the reliability of sensor systems
CN104097793A (en) * 2014-06-24 2014-10-15 上海微小卫星工程中心 Zero momentum magnetic control sun capture device and method of satellite
CN105444760A (en) * 2015-11-30 2016-03-30 北京控制工程研究所 Sun sensor based on collection of photovoltaic micro energy
CN105737822A (en) * 2016-02-06 2016-07-06 上海微小卫星工程中心 All-day ball view field sun sensor and determination method for sun vector thereof
CN105806369A (en) * 2016-05-20 2016-07-27 上海航天控制技术研究所 Star sensor in-orbit aberration correction method
CN106289237A (en) * 2016-07-27 2017-01-04 上海航天控制技术研究所 APS digital sun sensor solar angle computational methods and system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王俊等: "四象限差动式模拟太阳敏感器设计", 《传感技术学报》 *
肖红光等: "基于光线跟踪的星敏感视场仿真分析", 《***仿真学报》 *

Also Published As

Publication number Publication date
CN106843250B (en) 2019-12-10

Similar Documents

Publication Publication Date Title
US10220965B2 (en) Attitude determination using earth horizon sensors
CN105905317A (en) Sun-pointing control system for satellite and control method of sun-pointing control system
CN106124170B (en) A kind of camera optical axis direction computational methods based on high-precision attitude information
CN103389099B (en) Spacecraft attitude based on X-ray pulsar and position measuring system and method
CN104880191B (en) A kind of polarization auxiliary navigation method based on solar vector
Gounley et al. Autonomous satellite navigation by stellar refraction
CN106155074A (en) A kind of satellite three axle Direct to the sun control method ensureing satellite-ground link
CN103900571B (en) A kind of carrier posture measuring method based on the rotary-type SINS of inertial coodinate system
CN108871326B (en) A kind of single-shaft-rotation modulation inertia-astronomy deep integrated navigation method
CN106767705A (en) Imaging posture modeling method during a kind of Optical remote satellite point target observation mission is dynamic
US11396354B2 (en) Covert underwater navigation via polarimetry
CN106679645A (en) Multi-directional polarized light-based real-time navigation device
CN106197434A (en) The computational methods at attitude of satellite angle based on ground target point position
CN106197425A (en) The computational methods of ground target point position based on attitude of satellite angle
CN107085634A (en) It is quick to calculate sunshine and the method for sun synchronous satellite star sensor minimum angle
WO1994018073A1 (en) Autonomous on-board satellite control system
CN102589544A (en) Three-dimensional attitude acquisition method based on space characteristics of atmospheric polarization mode
CN101858747A (en) Method for resolving and determining counterglow oriented object posture of satellite sailboard by effectively using earth irradiation energy
CN104729457A (en) Method for determining position of sun relative to near-earth microsatellites
CN109813306A (en) A kind of unmanned vehicle planned trajectory satellite location data confidence level calculation method
CN109269510A (en) HEO satellite formation flying autonomous navigation method based on star sensor and inter-satellite link
Zhang et al. A self-contained interactive iteration positioning and orientation coupled navigation method based on skylight polarization
CN104123461A (en) Lighting visual relation calculating method for photometric analysis of space object
CN103913167A (en) Method for determining spatial attitude of aircraft in atmosphere by utilizing natural light polarization mode
CN106918317B (en) Side-sway observation covering band determines method and device to a kind of scan-type satellite over the ground

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