CN105866808A - 导航接收机定轨误差对卫星姿态精度的影响的确定方法 - Google Patents
导航接收机定轨误差对卫星姿态精度的影响的确定方法 Download PDFInfo
- Publication number
- CN105866808A CN105866808A CN201610452918.9A CN201610452918A CN105866808A CN 105866808 A CN105866808 A CN 105866808A CN 201610452918 A CN201610452918 A CN 201610452918A CN 105866808 A CN105866808 A CN 105866808A
- Authority
- CN
- China
- Prior art keywords
- rightarrow
- delta
- cos
- centerdot
- error
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Manufacturing & Machinery (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
一种导航接收机定轨误差对卫星姿态精度的影响的确定方法,充分利用导航接收机输出的位置速度与轨道根数之间的关系,提取了位置速度误差对姿态确定精度影响之间的表达式,通过他们之间的转换关系,成功分析解算出位置和速度的误差对姿态角在三个方向分别带来的传递关系,可以给出任意位置导航误差对姿态的影响,对工程设计中导航误差分配与姿态确定误差分析提供了直观方法。
Description
技术领域
本发明涉及一种导航接收机定轨误差对卫星姿态精度的影响的确定方法。
背景技术
轨道确定误差对姿态确定的精度有很大的影响,一般来说,定轨误差直接反映到姿态确定误差上。为了保证利用星上测量单机的量测数据的精确定向,必须精确地确定卫星在轨道上的位置,使定轨精度误差与星敏感器的测量精度误差相适应。由于定轨的误差会影响到姿态角的精度,仅利用位置、速度不能直接解算出误差传递函数,需要一种直观的方法处理卫星姿态与定轨误差之间的影响关系。
发明内容
本发明提供一种导航接收机定轨误差对卫星姿态精度的影响的确定方法,充分利用导航接收机输出的位置速度与轨道根数之间的关系,提取了位置速度误差对姿态确定精度影响之间的表达式,通过他们之间的转换关系,成功分析解算出位置和速度的误差对姿态角在三个方向分别带来的传递关系,可以给出任意位置导航误差对姿态的影响,对工程设计中导航误差分配与姿态确定误差分析提供了直观方法。
为了达到上述目的,本发明提供一种导航接收机定轨误差对卫星姿态精度的影响的确定方法,包含以下步骤:
步骤S1、确定位置速度与三轴姿态之间的关系;
根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi,进行泰勒展开并保留二阶以上小量,得到位置和速度对误差矩阵ΔAoi的影响关系,根据卫星姿态与矩阵Aoi的关系,确定位置和速度与三轴姿态的关系;
步骤S2、分别确定位置和速度各分量对三轴姿态的传递函数;
利用近圆轨道的特点和轨道根数与位置速度之间的关系,对位置速度误差与三轴姿态之间的传递函数进行化简,得到直观的表达式。
所述的步骤S1中,确定位置速度与姿态之间的关系的步骤具体包含以下步骤:
设根据理论的位置和速度矢量确定的理论轨道坐标系OXYZ的转换矩阵为Aoi,而根据实际获得的带有误差的位置和速度矢量确定的轨道坐标系的转换矩阵为由于位置速度误差相对位置速度是小量,可以用下式表示其误差ΔA:
另一方面,可以用“姿态角”的概念把Aoi与联系起来,把轨道坐标系的误差折算到相应的滚动、俯仰、偏航三个姿态角误差θ、ψ,当小偏差时,记姿态误差矩阵为C:
有
于是C=(Aoi+ΔA)Aoi -1=E+ΔAAoi T (4)
考虑精确到一阶小量,得:
把式(5)代入式(4)得到:
再由近似的C阵带入式(6)左端,将式(6)右端按矩阵运算法则展开,让左右两端矩阵中诸元相等,可以得到θ、ψ与位置误差[Δx,Δy,Δz]和速度误差的关系,矩阵展开如下:
让对应元素相等,解出滚动、俯仰、偏航三个姿态角误差θ、ψ:
所述的步骤S2中,分别确定位置和速度各分量对三轴姿态的传递函数的步骤具体包含:
步骤S2.1、分别计算位置误差[Δx,Δy,Δz]和速度误差引起的三轴姿态角误差;
步骤S2.2、将各量合并计算总的三轴姿态角误差θ、ψ。
所述的步骤S2.1具体包含以下步骤:
各方向位置误差[Δx,Δy,Δz]和速度误差引起的姿态角误差计算方法类似,记引起的滚动、俯仰、偏航三个姿态角误差分别为θj、ψj,其中这里重点说明X向位置误差Δx引起的姿态角误差θy、ψz;
为了解整个飞行过程中误差的变化规律,把Aoi用相应的根据轨道要素确定的转换矩阵来表示,另外确定位置矢量和速度矢量 与轨道要素纬度幅角u、升交点赤经Ω、倾角i之间的关系,对于近圆轨道,有:
由卫星轨道坐标系的定义可知,根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi为:
Aoi(1,i)=Aoi(2,j)×Aoi(3,k) (15)
其中i=1,2,3,对应于转换矩阵Aoi中每个行向量的三个分量,Aoi(2,j)、Aoi(3,k)为转换矩阵Aoi中第二行和第三行的行向量;
下面计算由于位置误差Δx引起的姿态角偏差矩阵ΔAx:
从而得到位置误差Δx引起的姿态角误差如下:
把和带入式(17),得到:
由于轨道是近圆轨道,速度和位置矢量之间夹角接近90°,因此有:
因此得到:
同理有:
同理得到其他量对三轴姿态的影响。
所述的步骤S2.2具体包含以下步骤:
将各量合并得到总的姿态角误差θ、ψ:
本发明充分利用导航接收机输出的位置速度与轨道根数之间的关系,提取了位置速度误差对姿态确定精度影响之间的表达式,通过他们之间的转换关系,成功分析解算出位置和速度的误差对姿态角在三个方向分别带来的传递关系,可以给出任意位置导航误差对姿态的影响,对工程设计中导航误差分配与姿态确定误差分析提供了直观方法。
具体实施方式
以下具体说明本发明的较佳实施例。
本发明提供一种导航接收机定轨误差对卫星姿态精度的影响的确定方法,包含以下步骤:
步骤S1、确定位置速度与三轴姿态之间的关系;
根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi,进行泰勒展开并保留二阶以上小量,得到位置和速度对误差矩阵ΔAoi的影响关系,根据卫星姿态与矩阵Aoi的关系,确定位置和速度与三轴姿态的关系;
步骤S2、分别确定位置和速度各分量对三轴姿态的传递函数;
利用近圆轨道的特点和轨道根数与位置速度之间的关系,对位置速度误差与三轴姿态之间的传递函数进行化简,得到直观的表达式。
所述的步骤S1中,确定位置速度与姿态之间的关系的步骤具体包含以下步骤:
设根据理论的位置和速度矢量确定的理论轨道坐标系OXYZ的转换矩阵为Aoi,而根据实际获得的带有误差的位置和速度矢量确定的轨道坐标系的转换矩阵为由于位置速度误差相对位置速度是小量,可以用下式表示其误差ΔA:
另一方面,可以用“姿态角”的概念把Aoi与联系起来,把轨道坐标系的误差折算到相应的滚动、俯仰、偏航三个姿态角误差θ、ψ,当小偏差时,记姿态误差矩阵为C:
有
于是C=(Aoi+ΔA)Aoi -1=E+ΔAAoi T (4)
考虑精确到一阶小量,得:
把式(5)代入式(4)得到:
再由近似的C阵带入式(6)左端,将式(6)右端按矩阵运算法则展开,让左右两端矩阵中诸元相等,可以得到θ、ψ与位置误差[Δx,Δy,Δz]和速度误差的关系,矩阵展开如下:
让对应元素相等,解出滚动、俯仰、偏航三个姿态角误差θ、ψ:
所述的步骤S2中,分别确定位置和速度各分量对三轴姿态的传递函数的步骤具体包含:
步骤S2.1、分别计算位置误差[Δx,Δy,Δz]和速度误差引起的三轴姿态角误差;
步骤S2.2、将各量合并计算总的三轴姿态角误差θ、ψ。
所述的步骤S2.1具体包含以下步骤:
各方向位置误差[Δx,Δy,Δz]和速度误差引起的姿态角误差计算方法类似,记引起的滚动、俯仰、偏航三个姿态角误差分别为θj、ψj,其中这里重点说明X向位置误差Δx引起的姿态角误差θy、ψz;
为了解整个飞行过程中误差的变化规律,把Aoi用相应的根据轨道要素确定的转换矩阵来表示,另外确定位置矢量和速度矢量 与轨道要素纬度幅角u、升交点赤经Ω、倾角i之间的关系,对于近圆轨道,有:
由卫星轨道坐标系的定义可知,根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi为:
Aoi(1,i)=Aoi(2,j)×Aoi(3,k) (15)
其中i=1,2,3,对应于转换矩阵Aoi中每个行向量的三个分量,Aoi(2,j)、Aoi(3,k)为转换矩阵Aoi中第二行和第三行的行向量;
下面计算由于位置误差Δx引起的姿态角偏差矩阵ΔAx:
从而得到位置误差Δx引起的姿态角误差如下:
把和带入式(17),得到:
由于轨道是近圆轨道,速度和位置矢量之间夹角接近90°,因此有:
因此得到:
同理有:
同理得到其他量对三轴姿态的影响。
所述的步骤S2.2具体包含以下步骤:
将各量合并得到总的姿态角误差θ、ψ:
本发明充分利用导航接收机输出的位置速度与轨道根数之间的关系,提取了位置速度误差对姿态确定精度影响之间的表达式,通过他们之间的转换关系,成功分析解算出位置和速度的误差对姿态角在三个方向分别带来的传递关系,可以给出任意位置导航误差对姿态的影响,对工程设计中导航误差分配与姿态确定误差分析提供了直观方法。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (5)
1.一种导航接收机定轨误差对卫星姿态精度的影响的确定方法,其特征在于,包含以下步骤:
步骤S1、确定位置速度与三轴姿态之间的关系;
根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi,进行泰勒展开并保留二阶以上小量,得到位置和速度对误差矩阵ΔAoi的影响关系,根据卫星姿态与矩阵Aoi的关系,确定位置和速度与三轴姿态的关系;
步骤S2、分别确定位置和速度各分量对三轴姿态的传递函数;
利用近圆轨道的特点和轨道根数与位置速度之间的关系,对位置速度误差与三轴姿态之间的传递函数进行化简,得到直观的表达式。
2.如权利要求1所述的导航接收机定轨误差对卫星姿态精度的影响的确定方法,其特征在于,所述的步骤S1中,确定位置速度与姿态之间的关系的步骤具体包含以下步骤:
设根据理论的位置和速度矢量确定的理论轨道坐标系OXYZ的转换矩阵为Aoi,而根据实际获得的带有误差的位置和速度矢量确定的轨道坐标系的转换矩阵为由于位置速度误差相对位置速度是小量,可以用下式表示其误差ΔA:
另一方面,可以用“姿态角”的概念把Aoi与联系起来,把轨道坐标系的误差折算到相应的滚动、俯仰、偏航三个姿态角误差θ、ψ,当小偏差时,记姿态误差矩阵为C:
有
于是
考虑精确到一阶小量,得:
把式(5)代入式(4)得到:
再由近似的C阵带入式(6)左端,将式(6)右端按矩阵运算法则展开,让左右两端矩阵中诸元相等,可以得到θ、ψ与位置误差[Δx,Δy,Δz]和速度误差的关系,矩阵展开如下:
让对应元素相等,解出滚动、俯仰、偏航三个姿态角误差θ、ψ:
3.如权利要求1所述的导航接收机定轨误差对卫星姿态精度的影响的确定方法,其特征在于,所述的步骤S2中,分别确定位置和速度各分量对三轴姿态的传递函数的步骤具体包含:
步骤S2.1、分别计算位置误差[Δx,Δy,Δz]和速度误差引起的三轴姿态角误差;
步骤S2.2、将各量合并计算总的三轴姿态角误差θ、ψ。
4.如权利要求3所述的导航接收机定轨误差对卫星姿态精度的影响的确定方法,其特征在于,所述的步骤S2.1具体包含以下步骤:
各方向位置误差[Δx,Δy,Δz]和速度误差引起的姿态角误差计算方法类似,记引起的滚动、俯仰、偏航三个姿态角误差分别为θj、ψj,其中这里重点说明X向位置误差Δx引起的姿态角误差θy、ψz;
为了解整个飞行过程中误差的变化规律,把Aoi用相应的根据轨道要素确定的转换矩阵来表示,另外确定位置矢量和速度矢量 与轨道要素纬度幅角u、升交点赤经Ω、倾角i之间的关系,对于近圆轨道,有:
由卫星轨道坐标系的定义可知,根据某时刻的位置和速度矢量求得卫星轨道坐标系相对地心惯性坐标系的转换矩阵Aoi为:
Aoi(1,i)=Aoi(2,j)×Aoi(3,k) (15)
其中i=1,2,3,对应于转换矩阵Aoi中每个行向量的三个分量,Aoi(2,j)、Aoi(3,k)为转换矩阵Aoi中第二行和第三行的行向量;
下面计算由于位置误差Δx引起的姿态角偏差矩阵ΔAx:
从而得到位置误差Δx引起的姿态角误差如下:
把和带入式(17),得到:
由于轨道是近圆轨道,速度和位置矢量之间夹角接近90°,因此有:
因此得到:
同理有:
同理得到其他量对三轴姿态的影响。
5.如权利要求3所述的导航接收机定轨误差对卫星姿态精度的影响的确定方法,其特征在于,所述的步骤S2.2具体包含以下步骤:
将各量合并得到总的姿态角误差θ、ψ:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610452918.9A CN105866808B (zh) | 2016-06-21 | 2016-06-21 | 导航接收机定轨误差对卫星姿态精度的影响的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610452918.9A CN105866808B (zh) | 2016-06-21 | 2016-06-21 | 导航接收机定轨误差对卫星姿态精度的影响的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866808A true CN105866808A (zh) | 2016-08-17 |
CN105866808B CN105866808B (zh) | 2018-04-20 |
Family
ID=56649475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610452918.9A Active CN105866808B (zh) | 2016-06-21 | 2016-06-21 | 导航接收机定轨误差对卫星姿态精度的影响的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866808B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525492A (zh) * | 2017-07-20 | 2017-12-29 | 航天东方红卫星有限公司 | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 |
CN110505004A (zh) * | 2019-08-13 | 2019-11-26 | 北京邮电大学 | 一种s波段卫星移动通信的卫星信号获取方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5184790A (en) * | 1991-07-22 | 1993-02-09 | Hughes Aircraft Company | Two-axis attitude correction for orbit inclination |
CN102175260A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 一种自主导航***误差校正方法 |
CN102494686A (zh) * | 2011-10-17 | 2012-06-13 | 北京国科环宇空间技术有限公司 | 一种卫星姿态轨道确定***及方法 |
CN103217982A (zh) * | 2013-02-21 | 2013-07-24 | 北京控制工程研究所 | 一种基于轮控模式的轨道控制方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 |
-
2016
- 2016-06-21 CN CN201610452918.9A patent/CN105866808B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5184790A (en) * | 1991-07-22 | 1993-02-09 | Hughes Aircraft Company | Two-axis attitude correction for orbit inclination |
CN102175260A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 一种自主导航***误差校正方法 |
CN102494686A (zh) * | 2011-10-17 | 2012-06-13 | 北京国科环宇空间技术有限公司 | 一种卫星姿态轨道确定***及方法 |
CN103217982A (zh) * | 2013-02-21 | 2013-07-24 | 北京控制工程研究所 | 一种基于轮控模式的轨道控制方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 |
Non-Patent Citations (7)
Title |
---|
D. MORTARI ET.AL: ""Attitude and Orbit Error in n-Dimensional Spaces"", 《JOURNAL OF THE ASTRONAUTICAL SCIENCES》 * |
TZU-PANG TSENG ET.AL: ""Assessing attitude error of FORMOSAT-3/COSMIC satellites and its impact on orbit determination"", 《ADVANCES IN SPACE RESEARCH》 * |
刘庆博: ""空间武器平台的位姿确定与传递对准"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
曹冲等: ""GPS技术在卫星定轨和姿态确定中的应用"", 《电波与天线》 * |
柳仲贵等: ""卫星轨道误差的相关性"", 《飞行器测控学报》 * |
盛传贞等: ""卫星姿态对JASON-2卫星精密定轨影响分析"", 《吉林大学学报(信息学报版)》 * |
黄晨: ""三轴稳定卫星姿态确定与控制***关键技术研究"", 《哈尔滨工程大学博士学位论文》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525492A (zh) * | 2017-07-20 | 2017-12-29 | 航天东方红卫星有限公司 | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 |
CN107525492B (zh) * | 2017-07-20 | 2020-09-18 | 航天东方红卫星有限公司 | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 |
CN110505004A (zh) * | 2019-08-13 | 2019-11-26 | 北京邮电大学 | 一种s波段卫星移动通信的卫星信号获取方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105866808B (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103256928B (zh) | 一种分布式惯性导航***及其姿态传递对准方法 | |
CN101788296B (zh) | 一种sins/cns深组合导航***及其实现方法 | |
CN104792340B (zh) | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 | |
CN103575299B (zh) | 利用外观测信息的双轴旋转惯导***对准及误差修正方法 | |
CN103676941B (zh) | 基于运动学和动力学模型的卫星控制***故障诊断方法 | |
CN104567930A (zh) | 一种能够估计和补偿机翼挠曲变形的传递对准方法 | |
CN101246012B (zh) | 一种基于鲁棒耗散滤波的组合导航方法 | |
CN103675861A (zh) | 一种基于星载gnss多天线的卫星自主定轨方法 | |
CN104165640A (zh) | 基于星敏感器的近空间弹载捷联惯导***传递对准方法 | |
CN101650569B (zh) | 三维空间中多运动体的寻迹编队控制方法 | |
CN104655152A (zh) | 一种基于联邦滤波的机载分布式pos实时传递对准方法 | |
CN102519485B (zh) | 一种引入陀螺信息的二位置捷联惯性导航***初始对准方法 | |
CN104748761B (zh) | 基于最优姿态匹配的动基座传递对准时延补偿方法 | |
CN104236586B (zh) | 基于量测失准角的动基座传递对准方法 | |
CN104457748A (zh) | 一种嵌入式瞄准吊舱测姿***及其传递对准方法 | |
CN106525034A (zh) | 一种基于对偶四元数的惯导***传递对准建模方法 | |
CN105184002A (zh) | 一种数传天线指向角度的仿真分析方法 | |
CN105115508A (zh) | 基于后数据的旋转制导炮弹快速空中对准方法 | |
CN105180728A (zh) | 基于前数据的旋转制导炮弹快速空中对准方法 | |
CN108181916A (zh) | 小卫星相对姿态的控制方法及装置 | |
CN110304279A (zh) | 一种电推进卫星的质心在轨标定补偿方法 | |
CN103968844B (zh) | 基于低轨平台跟踪测量的大椭圆机动航天器自主导航方法 | |
CN105069311A (zh) | 一种远程火箭发射初态误差传播估计方法 | |
CN109059914A (zh) | 一种基于gps和最小二乘滤波的炮弹滚转角估计方法 | |
CN108592946A (zh) | 一种基于两套旋转惯导冗余配置下的惯性器件漂移在线监控方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |