CN103438886A - Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model - Google Patents

Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model Download PDF

Info

Publication number
CN103438886A
CN103438886A CN2013103347599A CN201310334759A CN103438886A CN 103438886 A CN103438886 A CN 103438886A CN 2013103347599 A CN2013103347599 A CN 2013103347599A CN 201310334759 A CN201310334759 A CN 201310334759A CN 103438886 A CN103438886 A CN 103438886A
Authority
CN
China
Prior art keywords
attitude
satellite
thick
value
data
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
CN2013103347599A
Other languages
Chinese (zh)
Other versions
CN103438886B (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.)
STATE SATELLITE METEROLOGICAL CENTER
National Satellite Meteorological Center
Original Assignee
STATE SATELLITE METEROLOGICAL CENTER
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 STATE SATELLITE METEROLOGICAL CENTER filed Critical STATE SATELLITE METEROLOGICAL CENTER
Priority to CN201310334759.9A priority Critical patent/CN103438886B/en
Publication of CN103438886A publication Critical patent/CN103438886A/en
Application granted granted Critical
Publication of CN103438886B publication Critical patent/CN103438886B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

The invention relates to a determination method for attitudes of a spinning stabilized meteorological satellite based on a coarse-fine attitude relation model. The method is characterized by comprising the following steps that 1: the coarse attitude of the satellite is calculated; 2: a random error in the coarse attitude of the satellite is eliminated; and 3: the coarse attitude is modified by utilizing the coarse-fine attitude relation model. The method has the benefits that 1, the precision is high, and a calculated fine attitude result is basically consistent with the practical fine attitude of the satellite; 2, the applicability is high, and the fine attitude of the satellite can be calculated by the method to achieve accurate positioning of an observation cloud picture when the fine attitude cannot be obtained within 24h after attitude orbit control of the satellite or in an area observation mode; and 3, a processing procedure is distinct, the operand is not great, and implementation is facilitated.

Description

Spin stabilization weather satellite attitude based on thick smart attitude relational model is determined method
Technical field
The invention belongs to weather satellite observed image positioning field, be specifically related to a kind of spin stabilization weather satellite attitude based on thick smart attitude relational model and determine method.
Background technology
Accurately attitude determines it is that the weather satellite observed image can pinpoint basis, is also the basis that the weather satellite observed image can be applied.At present, China's spin stabilization weather satellite depends on pie chart and looks like to derive the accurate attitude parameter of satellite and carry out the framing to following 24 hours.This method accurate positioning is widely applied in business.But this method can't be derived the accurate attitude of satellite, the accurate location that can not realize satellite cloud picture when the pie chart data accumulation of not receiving figure or collecting is not enough.When satellite carries out region-wide figure observation or carried out after the rail control, in 48 hours, now there is no pie chart picture or disk image document very little at satellite, the accurate attitude of satellite is difficult to determine.
Utilize the satellite telemetry value can calculate the attitude of satellite, find that in practical business the attitude of satellite and the actual attitude of satellite that this remote measurement value is via satellite calculated have certain error, can produce the image lattice skew while utilizing this attitude to framing, therefore we call thick attitude to this attitude of the satellite that the satellite telemetry value calculates of utilizing, and are called smart attitude by the anti-attitude of satellite of releasing of image.Because spin stabilization weather satellite remote measurement value is unbroken, at any time telemometer is calculated the thick attitude of satellite via satellite.
In order to improve the computational accuracy of the thick attitude of satellite, Li Yuheng etc. are a lot of, and the expert has done deep research in this respect.Point out in correlative study that attitude determines that result is often inconsistent owing to being subject to how much observation condition restrictions and measuring error different with computing method.In order to eliminate this error, often all pass through the wide Changing Pattern of infrared earth sensor earth observation string on the research star, determine that spin axis and orbital method are to angle and the wide relation of earth measurement string, derive telemetry difference value calculating spin axle and the orbital method computing formula to angle, due to this method, can to eliminate the measurement data system poor, precision is improved, and can well meet the Satellite TT demand at present.But the attitude that this method is calculated still can not be directly used in framing, framing need to be determined the exact position of each pixel.The practical business operation shows, there is certain relation between thick attitude and smart attitude, if can set up the relational model of a kind of thick attitude and smart attitude, can in the time cannot calculating smart attitude, according to thick attitude, accurately estimate satellite essence attitude, thus the location of realizing observed image.
Summary of the invention
The purpose of this invention is to provide a kind of spin stabilization weather satellite attitude based on thick smart attitude relational model and determine method, set up the relational model of the thick attitude of satellite and smart attitude, when satellite does not have smart attitude, utilize relational model to derive the accurate attitude of satellite, realize that cloud atlas accurately locates.
The objective of the invention is to be achieved through the following technical solutions:
A kind of spin stabilization weather satellite attitude based on thick smart attitude relational model is determined method, comprises the following steps:
Step 1: calculate the thick attitude of satellite:
Decoding remote measurement bag, extract the information such as horizon instrument observation in satellite telemetering data, realize the thick Attitude Calculation algorithm of satellite, and export thick attitude result, wherein, the decoding remote measurement is to realize that remote measurement obtains, telemetry format decoding and quality inspection function, the method that its adopts is to judge the quality of every frame remote measurement according to the frame head synchronous code of telemetry frame and every frame end part with proof test value, for the process quality inspection, satisfactory telemetry frame, according to telemetry format, definition is decoded, then extract a plurality of T values in remote measurement, calculate corresponding physical quantity according to conversion formula,
The method that thick attitude solves employing is to utilize the sun in satellite telemetry-earth metrical information to solve the attitude of satellite, calculate the satellite rotating speed by the T value, the angle of satellite spin axle and solar direction, the angle between satellite spin axle and the earth's core and dihedral, then solve attitude parameter;
Step 2: eliminate the stochastic error in the thick attitude of satellite:
The thick attitude of the satellite calculated according to the remote measurement value can be subject to the impact of various factors to have stochastic error, the average complex filter of employing based on the slip intermediate value revised the thick attitude of satellite, by the weighted mean complex filter based on the slip intermediate value, the thick attitude data of abnormal satellite is differentiated, adopted the attitude moving window of a m length in algorithm
Figure 2013103347599100002DEST_PATH_IMAGE002
, utilize m continuous attitude value to determine the validity of the thick attitude data of satellite, if wave filter is judged this satellite, thick attitude data is effective, output; Otherwise, judge that the thick attitude data of this satellite is abnormal data, with weighted mean, revise;
Step 3: utilize thick attitude and the thick attitude of smart attitude relational model correction:
Thick attitude forecast attitude of satellite part is mainly the relation of the thick attitude of research and smart attitude, and calculates and forecast the attitude of satellite by thick attitude; Usually the accurate attitude of satellite comprises the mismatch parameter of satellite spin vector scanning radiometer, and the thick attitude of satellite is the attitude of satellite platform, does not contain the mismatch parameter of satellite scanning radiometer.Utilize smart attitude and the slightly comparison of attitude in historical data, the computing system error, used polynomial least square curve method matching Changing Pattern, and so far, spin stabilization weather satellite essence Attitude Calculation completes.
Beneficial effect of the present invention is:
1, precision is high, and the smart attitude result of calculating and the actual smart attitude of satellite are basically identical;
2, applicability is strong, at satellite, has carried out after the rail control in 24 hours or under the area observation pattern, in the time of can't obtaining smart attitude, can calculate satellite essence attitude by this method and realize observing cloud atlas accurately to locate;
3, processing procedure is clear, and operand is little, is convenient to implement.
The accompanying drawing explanation
Below with reference to the accompanying drawings the present invention is described in further detail.
Fig. 1 is the process flow diagram that the described spin stabilization weather satellite attitude based on thick smart attitude relational model of the embodiment of the present invention is determined method.
Embodiment
As shown in Figure 1, the described a kind of spin stabilization weather satellite attitude based on thick smart attitude relational model of the embodiment of the present invention is determined method, comprises the following steps:
Step 1: calculate the thick attitude of satellite:
Step 1.1): the satellite telemetering data obtained be take to frame and carry out quality inspection as unit, judge the quality of every frame satellite telemetering data according to frame head synchronous code and postamble proof test value, frame head and postamble be undesirable to be considered as off qualityly, will be rejected;
Step 1.2): for quality inspection, the satisfactory satellite telemetering data frame through step 1.1, according to the formal definition of satellite telemetering data, decoded, each radio frequency channel value is deposited in and specifies array I;
Step 1.3): from array I, extract T 1to T 11the former code value of scale-of-two, then convert to and calculate the needed physical quantity P of the thick attitude of satellite 1to P 11;
Step 1.4): use P 1to P 11the solar vector S of the unit of deriving, the nadir vector E of earth center position unit, the axial unit vector Z of satellite spin, spin axis and sun angle theta simultaneously that according to sun sensor, record s, the earth sensor spin axis and the nadir vector angle θ that record e, under the line in inertial coordinates system, by following attitude equation:
Figure 2013103347599100002DEST_PATH_IMAGE004
Can calculate the sun-spin axis plane and the interplanar included angle X of the earth's core-spin axis se;
Step 1.5): according to known position of sun (α s, δ s), by following formula:
Figure 2013103347599100002DEST_PATH_IMAGE006
Can try to achieve the thick attitude right ascension of satellite α ewith declination δ e, wherein:
Figure 2013103347599100002DEST_PATH_IMAGE008
=
Figure 2013103347599100002DEST_PATH_IMAGE010
=
Figure 2013103347599100002DEST_PATH_IMAGE014
=
Step 1.6): according to known position of sun (α s, δ s), inclination of satellite orbit i, true perigee angle f and right ascension of ascending node Ω, S and E can be expressed as respectively:
S=
Figure 2013103347599100002DEST_PATH_IMAGE020
E=
Figure 2013103347599100002DEST_PATH_IMAGE022
Therefore, the position of satellite is:
Step 2: eliminate the stochastic error in the thick attitude of satellite:
Described step 2 is differentiated the thick attitude data of abnormal satellite by the weighted mean complex filter based on the slip intermediate value, adopts the attitude moving window of a m length in algorithm , utilize m continuous attitude value to determine the validity of the thick attitude data of satellite, if wave filter is judged this satellite, thick attitude data is effective, output; Otherwise, judge that the thick attitude data of this satellite is abnormal data, with weighted mean, to revise, it further comprises:
Step 2.1): suppose the thick attitude data sequence of initial satellite
Figure 2013103347599100002DEST_PATH_IMAGE026
intermediate value be Z, by the thick attitude intermediate value of satellite absolute deviation tectonic sequence
Figure 2013103347599100002DEST_PATH_IMAGE028
=
Figure 2013103347599100002DEST_PATH_IMAGE030
, hypothetical sequence intermediate value be D, wherein, the thick attitude intermediate value of satellite absolute deviation Q=1.4826 * D, Q can replace standard deviation;
Step 2.2): if the thick attitude intermediate value of satellite absolute deviation sequence the value of the middle k of existence all is greater than
Figure 2013103347599100002DEST_PATH_IMAGE034
, , basis
Figure 2013103347599100002DEST_PATH_IMAGE038
value sequence
Figure 2013103347599100002DEST_PATH_IMAGE026A
split into
Figure 2013103347599100002DEST_PATH_IMAGE040
with
Figure 2013103347599100002DEST_PATH_IMAGE042
two sequences, wherein,
Figure 2013103347599100002DEST_PATH_IMAGE044
, ;
Step 2.3): the sequence obtained by step 2.2
Figure 2013103347599100002DEST_PATH_IMAGE040A
in value be normal value, this type of value need not be revised; Sequence
Figure 2013103347599100002DEST_PATH_IMAGE042A
in value be exceptional value, need to be revised the thick attitude data sequence of satellite
Figure 2013103347599100002DEST_PATH_IMAGE026AA
final revised result is
Figure 2013103347599100002DEST_PATH_IMAGE049
, and meet following condition:
If
Figure 2013103347599100002DEST_PATH_IMAGE051
,
Figure 2013103347599100002DEST_PATH_IMAGE053
;
If
Figure 2013103347599100002DEST_PATH_IMAGE055
Figure 2013103347599100002DEST_PATH_IMAGE042AA
, formula 1 is arranged, as follows:
Figure 2013103347599100002DEST_PATH_IMAGE057
Step 2.4): at first accumulate the thick attitude sequence of initial satellite in algorithm, the x, y and z axes of the thick attitude of satellite are accumulated respectively to sequence
Figure 2013103347599100002DEST_PATH_IMAGE026AAA
,
Figure 2013103347599100002DEST_PATH_IMAGE059
with
Figure 2013103347599100002DEST_PATH_IMAGE061
, when sequence runs up to length and is m, according to formula 1, once revise;
Step 2.5): fixed sequence program length is m, slippery sequence when having new thick attitude value to calculate
Figure DEST_PATH_IMAGE026AAAA
,
Figure DEST_PATH_IMAGE059A
with
Figure DEST_PATH_IMAGE061A
, new calculated value is put into to tail of the queue, and removes data of original head of the queue, and then carry out filter correction according to m value in 1 pair of sequence of formula.
The characteristic of the described weighted mean complex filter based on the slip intermediate value can and be revised these three parameters of weights C and be adjusted with window width m, thresholding L.Wherein: m affects total consistance of wave filter; Threshold parameter L directly determines the wave filter degree of initiatively keeping forging ahead, and the L value increases, and will be judged to be singular data and the possibility that replaces by intermediate value reduces, and when L=0, wave filter is determined all the time; Revise weights C and directly determine wave filter correction amplitude, the C value increases, and wave filter correction amplitude is larger, when C=0, uses the abnormal attitude of average correction.
Step 3: utilize thick attitude and the thick attitude of smart attitude relational model correction:
Thick attitude and smart attitude that described step 3 will have been removed stochastic error compare, and the computing system error is used the Changing Pattern of two kinds of attitudes of polynomial least square curve method matching, and it further comprises:
Step 3.1): smart attitude and the relation table of the x, y and z axes of thick attitude are shown to three groups of discrete data points
Figure 2013103347599100002DEST_PATH_IMAGE063
,
Figure 2013103347599100002DEST_PATH_IMAGE065
with , respectively by three relation functions
Figure 2013103347599100002DEST_PATH_IMAGE069
,
Figure 2013103347599100002DEST_PATH_IMAGE071
with
Figure 2013103347599100002DEST_PATH_IMAGE073
with three groups of data fittings of give;
Step 3.2): to attitude data
Figure 2013103347599100002DEST_PATH_IMAGE075
,
Figure 685538DEST_PATH_IMAGE065
with
Figure 678902DEST_PATH_IMAGE067
, at all number of times, be no more than
Figure 2013103347599100002DEST_PATH_IMAGE077
the function class that forms of polynomial expression in, ask respectively
Figure 2013103347599100002DEST_PATH_IMAGE081
,
Figure 2013103347599100002DEST_PATH_IMAGE083
with
Figure 2013103347599100002DEST_PATH_IMAGE085
, make error ,
Figure 2013103347599100002DEST_PATH_IMAGE089
with
Figure 2013103347599100002DEST_PATH_IMAGE091
the quadratic sum minimum,
Figure 2013103347599100002DEST_PATH_IMAGE093
Figure 2013103347599100002DEST_PATH_IMAGE095
Figure 2013103347599100002DEST_PATH_IMAGE097
Wherein,
Figure 2013103347599100002DEST_PATH_IMAGE099
,
Figure 2013103347599100002DEST_PATH_IMAGE101
,
Figure 2013103347599100002DEST_PATH_IMAGE103
(i=0,1 ..., m);
Figure 11661DEST_PATH_IMAGE081
,
Figure 935624DEST_PATH_IMAGE083
with
Figure 396692DEST_PATH_IMAGE085
,
Make error
Figure 13487DEST_PATH_IMAGE087
, with minimum,
Figure 2013103347599100002DEST_PATH_IMAGE105
,
Figure 2013103347599100002DEST_PATH_IMAGE107
,
Figure 2013103347599100002DEST_PATH_IMAGE109
it is required polynomial expression.Owing between X, Y and Z three, having relation
Figure 2013103347599100002DEST_PATH_IMAGE111
, the X after matching, Y and Z value generally are difficult to meet this relation.So with in only need to select in X, Y and Z two to calculate, remain one and pass through
Figure 428834DEST_PATH_IMAGE111
this relation is asked calculation.Actual with in general selection X, Y calculate, and then calculate Z by X, Y.
The present invention is not limited to above-mentioned preferred forms; anyone can draw other various forms of products under enlightenment of the present invention; no matter but do any variation on its shape or structure; every have identical with a application or akin technical scheme, within all dropping on protection scope of the present invention.

Claims (6)

1. the spin stabilization weather satellite attitude based on thick smart attitude relational model is determined method, it is characterized in that, comprises the following steps:
Step 1: calculate the thick attitude of satellite;
Step 2: eliminate the stochastic error in the thick attitude of satellite; And
Step 3: utilize thick attitude and the thick attitude of smart attitude relational model correction.
2. the spin stabilization weather satellite attitude based on thick smart attitude relational model according to claim 1 is determined method, it is characterized in that, described step 1 further comprises:
Step 1.1): the satellite telemetering data obtained be take to frame and carry out quality inspection as unit, judge the quality of every frame satellite telemetering data according to frame head synchronous code and postamble proof test value, be considered as off qualityly for undesirable satellite telemetering data value, it is given up;
Step 1.2): for quality inspection, the satisfactory satellite telemetering data frame through step 1.1, according to the formal definition of satellite telemetering data, decoded, each radio frequency channel value is deposited in and specifies array I;
Step 1.3): from array I, extract T 1to T 11the former code value of scale-of-two, then convert to and calculate the needed physical quantity P of the thick attitude of satellite 1to P 11;
Step 1.4): according to ripe Attitude Calculation theory, use P 1to P 11the solar vector S of the unit of deriving, the nadir vector E of earth center position unit, the axial unit vector Z of satellite spin, spin axis and sun angle theta simultaneously that according to sun sensor, record s, the earth sensor spin axis and the nadir vector angle θ that record e, under the line in inertial coordinates system, by following attitude equation:
Can calculate the sun-spin axis plane and the interplanar included angle X of the earth's core-spin axis se;
Step 1.5): according to known position of sun (α s, δ s), by following formula:
Figure 2013103347599100001DEST_PATH_IMAGE004
Can try to achieve the thick attitude right ascension of satellite α ewith declination δ e, wherein:
Figure 2013103347599100001DEST_PATH_IMAGE006
=
Figure 2013103347599100001DEST_PATH_IMAGE008
Figure 2013103347599100001DEST_PATH_IMAGE010
=
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE014
=
Figure DEST_PATH_IMAGE016
Step 1.6): according to known position of sun (α s, δ s), inclination of satellite orbit i, true perigee angle f and right ascension of ascending node Ω, S and E can be expressed as respectively:
S=
E=
Figure DEST_PATH_IMAGE020
Therefore, the position of satellite is:
Figure DEST_PATH_IMAGE022
Spin stabilization weather satellite attitude based on thick smart attitude relational model according to claim 2 is determined method, it is characterized in that: in described step 2, by the weighted mean complex filter based on the slip intermediate value, the thick attitude data of abnormal satellite is differentiated, adopted the attitude moving window of a m length in algorithm , utilize m continuous attitude value to determine the validity of the thick attitude data of satellite, if wave filter is judged this satellite, thick attitude data is effective, output; Otherwise, judge that the thick attitude data of this satellite is abnormal data, with weighted mean, revise.
3. the spin stabilization weather satellite attitude based on thick smart attitude relational model according to claim 2 is determined method, it is characterized in that, described step 2 comprises:
Step 2.1): suppose the thick attitude data sequence of initial satellite
Figure DEST_PATH_IMAGE026
intermediate value be Z, by the thick attitude intermediate value of satellite absolute deviation tectonic sequence
Figure DEST_PATH_IMAGE028
=
Figure DEST_PATH_IMAGE030
, hypothetical sequence
Figure DEST_PATH_IMAGE028A
intermediate value be D, wherein, the thick attitude intermediate value of satellite absolute deviation Q=1.4826 * D, Q can replace standard deviation;
Step 2.2): if the thick attitude intermediate value of satellite absolute deviation sequence
Figure DEST_PATH_IMAGE032
the value of the middle k of existence all is greater than
Figure DEST_PATH_IMAGE034
, , basis value sequence split into
Figure DEST_PATH_IMAGE040
with
Figure DEST_PATH_IMAGE042
two sequences, wherein, ,
Figure DEST_PATH_IMAGE046
;
Step 2.3): the sequence obtained by step 2.2
Figure DEST_PATH_IMAGE040A
in value be normal value, this type of value need not be revised; Sequence
Figure DEST_PATH_IMAGE042A
in value be exceptional value, need to be revised the thick attitude data sequence of satellite
Figure DEST_PATH_IMAGE026AA
final revised result is , and meet following condition:
If
Figure DEST_PATH_IMAGE051
,
Figure DEST_PATH_IMAGE053
;
If
Figure DEST_PATH_IMAGE055
Figure DEST_PATH_IMAGE042AA
, formula 1 is arranged, as follows:
Figure DEST_PATH_IMAGE057
Step 2.4): at first accumulate the thick attitude sequence of initial satellite in algorithm, the x, y and z axes of the thick attitude of satellite are accumulated respectively to sequence
Figure DEST_PATH_IMAGE059
, with
Figure DEST_PATH_IMAGE063
, when sequence runs up to length and is m, according to formula 1, once revise;
Step 2.5): fixed sequence program length is m, slippery sequence when having new thick attitude value to calculate
Figure DEST_PATH_IMAGE026AAA
, with
Figure DEST_PATH_IMAGE067
, new calculated value is put into to tail of the queue, and removes data of original head of the queue, and then carry out filter correction according to m value in 1 pair of sequence of formula.
4. the spin stabilization weather satellite attitude based on thick smart attitude relational model according to claim 4 is determined method, it is characterized in that: window width m, the thresholding L for characteristic of the described weighted mean complex filter based on the slip intermediate value and these three parameters of correction weights C are adjusted, and wherein: m affects total consistance of wave filter; Threshold parameter L directly determines the wave filter degree of initiatively keeping forging ahead, and the L value increases, and will be judged to be singular data and the possibility that replaces by intermediate value reduces, and when L=0, wave filter is determined all the time; Revise weights C and directly determine wave filter correction amplitude, the C value increases, and wave filter correction amplitude is larger, when C=0, uses the abnormal attitude of average correction.
5. the spin stabilization weather satellite attitude based on thick smart attitude relational model according to claim 5 is determined method, it is characterized in that, in described step 3, thick attitude and the smart attitude of having removed stochastic error are compared, the computing system error, used the Changing Pattern of two kinds of attitudes of polynomial least square curve method matching.
6. determine method according to the described spin stabilization weather satellite attitude based on thick smart attitude relational model of claim 5 or 6, it is characterized in that, described step 3 comprises:
Step 3.1): smart attitude and the relation table of the x, y and z axes of thick attitude are shown to three groups of discrete data points ,
Figure DEST_PATH_IMAGE071
with
Figure DEST_PATH_IMAGE073
, respectively by three relation functions
Figure DEST_PATH_IMAGE075
, with
Figure DEST_PATH_IMAGE079
with three groups of data fittings of give;
Step 3.2): to attitude data
Figure DEST_PATH_IMAGE081
,
Figure 979260DEST_PATH_IMAGE071
with
Figure 896400DEST_PATH_IMAGE073
, at all number of times, be no more than
Figure DEST_PATH_IMAGE083
the function class that forms of polynomial expression in, ask respectively
Figure DEST_PATH_IMAGE087
,
Figure DEST_PATH_IMAGE089
with
Figure DEST_PATH_IMAGE091
, make error
Figure DEST_PATH_IMAGE093
,
Figure DEST_PATH_IMAGE095
with
Figure DEST_PATH_IMAGE097
the quadratic sum minimum,
Figure DEST_PATH_IMAGE099
Figure DEST_PATH_IMAGE101
Figure DEST_PATH_IMAGE103
Wherein,
Figure DEST_PATH_IMAGE105
,
Figure DEST_PATH_IMAGE107
,
Figure DEST_PATH_IMAGE109
(i=0,1 ..., m);
Figure 569082DEST_PATH_IMAGE087
,
Figure 151242DEST_PATH_IMAGE089
with
Figure 965614DEST_PATH_IMAGE091
,
Make error , with
Figure 384503DEST_PATH_IMAGE097
minimum,
Figure DEST_PATH_IMAGE111
,
Figure DEST_PATH_IMAGE113
,
Figure DEST_PATH_IMAGE115
it is required polynomial expression.
CN201310334759.9A 2013-08-02 2013-08-02 Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model Expired - Fee Related CN103438886B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310334759.9A CN103438886B (en) 2013-08-02 2013-08-02 Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310334759.9A CN103438886B (en) 2013-08-02 2013-08-02 Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model

Publications (2)

Publication Number Publication Date
CN103438886A true CN103438886A (en) 2013-12-11
CN103438886B CN103438886B (en) 2017-04-19

Family

ID=49692582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310334759.9A Expired - Fee Related CN103438886B (en) 2013-08-02 2013-08-02 Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model

Country Status (1)

Country Link
CN (1) CN103438886B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106628258A (en) * 2016-10-10 2017-05-10 北京控制工程研究所 Satellite spin attitude determination method based on sun vector information
CN107357758A (en) * 2017-06-29 2017-11-17 中国人民解放军63796部队 The multinomial least square regression Memorability of location information seeks fast method
CN110203424A (en) * 2019-05-05 2019-09-06 中国人民解放军63921部队 Utilize the method and apparatus of measurement data estimation spacecraft spin motion
CN110297826A (en) * 2019-05-31 2019-10-01 南京理工大学 Method based on json dynamic analysis satellite telemetering data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0840094A2 (en) * 1996-10-31 1998-05-06 Samsung Electronics Co., Ltd. Navigation system
CN102538792A (en) * 2012-02-08 2012-07-04 北京航空航天大学 Filtering method for position attitude system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0840094A2 (en) * 1996-10-31 1998-05-06 Samsung Electronics Co., Ltd. Navigation system
CN102538792A (en) * 2012-02-08 2012-07-04 北京航空航天大学 Filtering method for position attitude system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
L.KAFFER ETC: "ROSAT卫星轨在姿态测量功能的恢复", 《控制工程》 *
李于衡等: "FY-2气象卫星在轨管理工程测控关键技术(上)", 《上海航天》 *
李于衡等: "地球静止轨道自旋卫星姿态确定及控制策略", 《飞行器测控学报》 *
桂志先等: "高精度速度数据体系统的建立与应用", 《石油工业计算机应用》 *
赵桂玲等: "基于级联FIR 滤波器的系泊精对准算法", 《中国惯性技术学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106628258A (en) * 2016-10-10 2017-05-10 北京控制工程研究所 Satellite spin attitude determination method based on sun vector information
CN106628258B (en) * 2016-10-10 2019-05-24 北京控制工程研究所 A kind of satellite spin attitude determination method based on solar vector information
CN107357758A (en) * 2017-06-29 2017-11-17 中国人民解放军63796部队 The multinomial least square regression Memorability of location information seeks fast method
CN107357758B (en) * 2017-06-29 2021-04-13 中国人民解放军63796部队 Polynomial least square regression memorability speed calculation method for positioning information
CN110203424A (en) * 2019-05-05 2019-09-06 中国人民解放军63921部队 Utilize the method and apparatus of measurement data estimation spacecraft spin motion
CN110297826A (en) * 2019-05-31 2019-10-01 南京理工大学 Method based on json dynamic analysis satellite telemetering data
CN110297826B (en) * 2019-05-31 2020-12-11 南京理工大学 Method for dynamically analyzing satellite telemetry data based on json

Also Published As

Publication number Publication date
CN103438886B (en) 2017-04-19

Similar Documents

Publication Publication Date Title
CN108332649B (en) Landslide deformation comprehensive early warning method and system
CN105629263B (en) A kind of troposphere atmosphere delay estimation error correcting method and correction system
CN105486312B (en) A kind of star sensor and high frequency angular displacement sensor integrated attitude determination method and system
CN103969643B (en) One carries out X-band pathfinder inverting ocean wave parameter method based on novel wave dispersion relation band filter
CN103631681B (en) A kind of method of online reparation abnormal data of wind power plant
CN103674021A (en) Integrated navigation system and method based on SINS (Strapdown Inertial Navigation System) and star sensor
CN103106632B (en) A kind of fusion method of the different accuracy three dimensional point cloud based on average drifting
CN103438886A (en) Determination method for attitudes of spinning stabilized meteorological satellite based on coarse-fine attitude relation model
CN102506876B (en) Self-contained navigation method for measurement of earth ultraviolet sensor
CN103399300A (en) Wave packet superposition microseism ground location method
CN107656278B (en) Quantitative precipitation estimation method based on dense rainfall station
CN103954953A (en) Method for performing blind source error compensation on airborne laser radar based on data driving
CN102707336A (en) Novel method of using A-Train series satellite data for synergetic inversion of cloud phase states and cloud parameters
CN103617336A (en) Method for drawing aircraft noise contour map
CN104112062A (en) Method for obtaining wind resource distribution based on interpolation method
CN105068032A (en) Photovoltaic convergence box current acquisition channel temperature drift coefficient calibration method
CN105547287A (en) Irregular coelonavigation sight line information extracting method
CN102168979B (en) Isoline matching method for passive navigation based on triangular constraint model
CN103197303A (en) Earth surface two-direction reflection characteristic retrieval method and earth surface two-direction reflection characteristic retrieval system based on multiple sensors
CN103338024A (en) Complementation Kalman filtering device and method of time delay in antenna array
CN104933235B (en) A kind of method for merging coastal waters multi-satellite sea level height abnormal data
CN104133993A (en) Small-field fixed star matching method based on error vector matching
CN103364842A (en) Error separation method of strapdown airborne gravitometer
CN102445712B (en) Character window weighting related spectrum matching method facing rocks and minerals
CN107328409B (en) A kind of X-ray pulsar navigation method based on dynamic pulse accumulation window

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170419

Termination date: 20170802