CN106197434A - 基于地面目标点位置的卫星姿态角的计算方法 - Google Patents
基于地面目标点位置的卫星姿态角的计算方法 Download PDFInfo
- Publication number
- CN106197434A CN106197434A CN201610502752.7A CN201610502752A CN106197434A CN 106197434 A CN106197434 A CN 106197434A CN 201610502752 A CN201610502752 A CN 201610502752A CN 106197434 A CN106197434 A CN 106197434A
- Authority
- CN
- China
- Prior art keywords
- satellite
- target point
- coordinate system
- ground target
- axis
- 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.)
- Pending
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种基于地面目标点位置的卫星姿态角的计算方法,其方法过程如下:已知地面目标点位置(大地坐标系下经度B、纬度L、高程H),),将地面目标点在大地坐标系中的经度、纬度和高程转换为地固坐标系中的直角坐标;将地面目标点在地固坐标系中的直角坐标转换为地心J2000惯性坐标系中的位置矢量;根据卫星在地心J2000惯性坐标系中的位置速度矢量和速度矢量求出地心J2000惯性坐标系到卫星轨道坐标系的转换矩阵C;用转换矩阵C将地面目标点在地心J2000惯性坐标系中的位置矢量转换成在卫星轨道坐标系中的位置矢量;根据地面目标点在卫星轨道坐标系中的位置矢量求得卫星的姿态角。
Description
技术领域
本发明属于航天测量与控制应用领域,涉及航天器对地面目标观测的姿态机动计算方法,适用于地球资源航天器及遥感航天器在飞行过程中对地面目标观测的姿态机动计算。
背景技术
随着空间对地观测技术的迅猛发展,遥感航天器越来越多,对地面目标的访问模式也由以往的固定星下点模式、单侧摆模式发展可多自由度姿态调节传感器模式。以往的固定星下点模式、单侧摆模式无法对同时能够两个方向、多自由度姿态调节传感器模式进行建模计算,且计算过程冗长复杂。随着新型航天器和传感器的出现,需要研究新的可多自由度姿态调节的对地面目标观测姿态机动计算方法。
发明内容
本发明所要解决的技术问题是航天器对地面目标观测的姿态计算方法,建立了地面目标点和卫航天器坐标转换关系,依据空间几何关系给出了航天器对地面目标观测的方程,并得到了该方程的分析解,计算航天器观测该目标所需姿态。该结果可用于多自由度姿态调节传感器模式的航天器对地面目标观测的姿态计算中。
本发明采用的技术方案为:
基于地面目标点位置的卫星姿态角的计算方法,包括以下步骤:
(1)将地面目标点在大地坐标系中的经度、纬度和高程转换为地固坐标系中的直角坐标;
(2)将地面目标点在地固坐标系中的直角坐标转换为地心J2000惯性坐标系中的位置矢量;
(3)根据卫星在地心J2000惯性坐标系中的位置速度矢量和速度矢量求出地心J2000惯性坐标系到卫星轨道坐标系的转换矩阵C;
(4)用转换矩阵C将地面目标点在地心J2000惯性坐标系中的位置矢量转换成卫星在轨道坐标系中的位置矢量;
(5)根据地面目标点在卫星轨道坐标系中的位置矢量求得卫星的姿态角;
完成基于地面目标点位置的卫星姿态角的计算。
其中,步骤(5)具体为:卫星绕X轴旋转时的姿态旋转角为横滚角绕Y轴旋转时的姿态旋转角为俯仰角θ,绕Z轴旋转时的姿态旋转角为偏航角ψ;
卫星绕两个轴旋转,有以下四种方法可得到卫星的姿态角;
方法1:卫星先绕Y轴旋转,再绕X轴旋转;
由:
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法2:卫星先绕X轴旋转,再绕Y轴旋转:
由:
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法3:卫星先绕Z轴旋转,再绕Y轴旋转:
由
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法4:卫星先绕Z轴旋转,再绕X轴旋转:
由
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值。
与背景技术相比本发明的有益效果为:
(1)基于空间几何关系给出航天器对地观测方程。该方法从几何观测出发,适应航天器两个姿态角可调节对地观测的需求。
(2)航天器对地实时观测计算效率及精度高。
(3)增长卫星的对地观测时间,实现卫星对地面目标的凝视效果。
附图说明
图1是卫星、地面目标和地球空间几何关系。
图2是已知地面目标点位置卫星J2000系位置速度求卫星姿态角度步骤图。
具体实施方式
1 用到的坐标系统和转换
1.1 坐标系的定义
表1 坐标系的定义
1.2 定义旋转矩阵
在右旋直角坐标系中,以坐标轴的旋转而实现的坐标变换,可以通过旋转矩阵来表示新老坐标系之间的变换关系。定义如下三个旋转矩阵:
其中RK(θ)表示绕K轴正向旋转θ角(逆时针方向)。例,设新坐标系由原坐标系绕Z轴逆时针旋转θ角而得,则目标在新坐标系中的坐标(X’,Y’,Z’)与原坐标系中的坐标(X,Y,Z)存在如下关系:
1.3 J2000.0惯性坐标系和地固坐标系之间的转换
这两个坐标系之间的差别是极移、地球自转以及岁差章动。对于J2000.0惯性坐标系和地固系中的位置矢量rJ2000和rEG,有:
rEG=(EP)(ER)(NR)(PR)rJ2000
(PR)为岁差矩阵,(NR)为章动矩阵,(ER)为地球自转矩阵,(EP)为极移矩阵。它们分别由下列各式表达:
(EP)=Ry(-xp)Rx(-yp)
(ER)=Rz(SG)
(NR)=Rx(-Δε)Ry(Δθ)Rz(-Δμ)
(PR)=Rz(-zA)Ry(θA)Rz(-ζA)
各式中的xp,yp为极移分量,SG为格林尼治恒星时,ζA,θA,zA为岁差量,Δε,Δθ,Δμ为章动量。
1.4 大地坐标系与地固坐标系直角坐标的转换
(X,Y,Z)为地固坐标系直角坐标;(B,L,H)为大地纬度、大地经度、大地高程;Re=6378137m为地球赤道半径,ec为地球子午圈的偏心率。则
1.5 J2000.0惯性坐标系和卫星轨道坐标系之间的转换
记位置矢量从J2000.0惯性坐标系到卫星轨道坐标系的转换矩阵为C,由卫星轨道坐标系的定义可得转换矩阵C的各元素为:
C(1,i)=C(2,j)×C(3,k)
其中:i=1,2,3对应着转换矩阵C中每个行向量的三个分量。
C(1,i),C(2,j),C(3,k)分别为矩阵C中第一行、第二行和第三行的行向量。
1.6 卫星轨道坐标系和卫星本体坐标系之间的转换
在卫星轨道系中,定义绕X轴旋转为滚动角φ(也可称为侧摆),定义绕Y轴旋转俯仰角θ,绕Z轴旋转偏航角:ψ,经坐标旋转到卫星本体系。根据旋转次序不同有个定义。
1.7 卫星本体坐标系和传感器基准坐标系之间的转换
可认为卫星本体坐标系和传感器基准坐标系的坐标原点相同,转换关系可根据安装位置定义的通过坐标旋转得到。传感器安装于卫星本体系Z轴。
2、已知地面目标点位置航天器J2000系位置速度求卫星姿态角度
航天器、地面目标和地球空间几何关系如图1所示,已知地面目标点位置(大地坐标系下经度B、纬度L、高程H),J2000坐标系下航天器某一时刻的位置速度矢量rJ2000和速度矢量
步骤如下:
1)将地面目标点在大地坐标系中的经度B、纬度L和高程H转换为地固坐标系中的直角坐标;
2)将地面目标点在地固坐标系中的直角坐标转换为地心J2000惯性坐标系中的位置矢量robj;
3)根据卫星在地心J2000坐标系中的位置速度矢量rJ2000和速度矢量求出地心J2000惯性坐标系到卫星轨道坐标系的转换矩阵;
4)用转换矩阵C将地面目标点在地心J2000惯性坐标系中的位置矢量robj转换成卫星在轨道坐标系中的位置矢量r(x,y,z);
5)根据地面目标点在卫星轨道坐标系中的位置矢量r(x,y,z)求得卫星的姿态角,定义卫星绕X轴旋转时的姿态旋转角为横滚角绕Y轴旋转时的姿态旋转角为俯仰角θ,绕Z轴旋转时的姿态旋转角为偏航角ψ;
卫星绕两个轴旋转,有以下四种方法可得到卫星的姿态角;
方法1:卫星先绕Y轴旋转,再绕X轴旋转;
由:
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法2:卫星先绕X轴旋转,再绕Y轴旋转:
由:
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法3:卫星先绕Z轴旋转,再绕Y轴旋转:
由
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值;
方法4:卫星先绕Z轴旋转,再绕X轴旋转:
由
得
式中,(x,y,z)为地面目标点在卫星轨道坐标系中的位置矢量,r为位置矢量(x,y,z)的绝对值。
完成基于地面目标点位置的卫星姿态角度的计算。
现结合实施例、附图对本发明作进一步描述:
已知地面目标点位置(大地坐标系下经度B、纬度L、高程H),J2000坐标系下航天器某一时刻的位置速度矢量rJ2000和速度矢量按照图2的步骤计算出航天器的姿态角(φ,θ,ψ)。
例:北京时2016年1月1日21:30航天器位置速度为/4406328.0,5117483.0,1311255.0,1699.0,448.0,-7406.0/若已知地面目标点坐标为(B=10.0度L=120.0度H=0.0米)则计算出的四种转序的姿态角见表2。
表2 由地面点求出的四种转序的姿态角
Claims (2)
1.基于地面目标点位置的卫星姿态角的计算方法,其特征在于,包括以下步骤:
(1)将地面目标点在大地坐标系中的经度、纬度和高程转换为地固坐标系中的直角坐标;
(2)将地面目标点在地固坐标系中的直角坐标转换为地心J2000惯性坐标系中的位置矢量;
(3)根据卫星在地心J2000惯性坐标系中的位置速度矢量和速度矢量求出地心J2000惯性坐标系到卫星轨道坐标系的转换矩阵C;
(4)用转换矩阵C将地面目标点在地心J2000惯性坐标系中的位置矢量转换成在卫星轨道坐标系中的位置矢量;
(5)根据地面目标点在卫星轨道坐标系中的位置矢量求得卫星的姿态角;
完成基于地面目标点位置的卫星姿态角的计算。
2.根据权利要求1所述的基于地面目标点位置的卫星姿态角的计算方法,其特征在于,步骤(5)具体为:卫星绕X轴旋转时的姿态旋转角为横滚角绕Y轴旋转时的姿态旋转角为俯仰角θ,绕Z轴旋转时的姿态旋转角为偏航角ψ;
卫星绕两个轴旋转,根据以下方法1、方法2、方法3或方法四可得到卫星的姿态角;
方法1:卫星先绕Y轴旋转,再绕X轴旋转;
由:
得
式中,x为地面目标点在卫星轨道坐标系中x轴方向的位置矢量,y为地面目标点在卫星轨道坐标系中y轴方向的位置矢量,z为地面目标点在卫星轨道坐标系中z轴方向的位置矢量;r为位置矢量的绝对值;
方法2:卫星先绕X轴旋转,再绕Y轴旋转:
由:
得
式中,x为地面目标点在卫星轨道坐标系中x轴方向的位置矢量,y为地面目标点在卫星轨道坐标系中y轴方向的位置矢量,z为地面目标点在卫星轨道坐标系中z轴方向的位置矢量;r为位置矢量的绝对值;
方法3:卫星先绕Z轴旋转,再绕Y轴旋转:
由
得
式中,x为地面目标点在卫星轨道坐标系中x轴方向的位置矢量,y为地面目标点在卫星轨道坐标系中y轴方向的位置矢量,z为地面目标点在卫星轨道坐标系中z轴方向的位置矢量;r为位置矢量的绝对值;
方法4:卫星先绕Z轴旋转,再绕X轴旋转:
由
得
式中,x为地面目标点在卫星轨道坐标系中x轴方向的位置矢量,y为地面目标点在卫星轨道坐标系中y轴方向的位置矢量,z为地面目标点在卫星轨道坐标系中z轴方向的位置矢量;r为位置矢量的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610502752.7A CN106197434A (zh) | 2016-06-30 | 2016-06-30 | 基于地面目标点位置的卫星姿态角的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610502752.7A CN106197434A (zh) | 2016-06-30 | 2016-06-30 | 基于地面目标点位置的卫星姿态角的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106197434A true CN106197434A (zh) | 2016-12-07 |
Family
ID=57463243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610502752.7A Pending CN106197434A (zh) | 2016-06-30 | 2016-06-30 | 基于地面目标点位置的卫星姿态角的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106197434A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287352A (zh) * | 2018-02-06 | 2018-07-17 | 西安四方星途测控技术有限公司 | 一种卫星圆锥传感器作用范围的确认方法及系统 |
CN108491633A (zh) * | 2018-03-24 | 2018-09-04 | 中国人民解放军海军航空大学 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
CN109001762A (zh) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | 一种同步带目标光学观测卫星的地气光抑制方法及系统 |
CN110162069A (zh) * | 2019-05-10 | 2019-08-23 | 北京航空航天大学 | 一种近地轨道航天器阳光反射凝视期望姿态解析求解方法 |
CN111552272A (zh) * | 2020-05-07 | 2020-08-18 | 上海航天电子有限公司 | 无源相控阵天线波束控制算法 |
CN111751809A (zh) * | 2020-06-09 | 2020-10-09 | 军事科学院系统工程研究院后勤科学与技术研究所 | 一种点源靶标反射镜的调整角度计算方法 |
CN111966953A (zh) * | 2020-08-07 | 2020-11-20 | 四川泛华航空仪表电器有限公司 | 一种基于坐标转换的飞机姿态角计算方法 |
CN114282166A (zh) * | 2021-12-22 | 2022-04-05 | 中国西安卫星测控中心 | 基于环境函数矩阵的卫星经度计算方法、变轨判别方法 |
CN116975504A (zh) * | 2023-09-22 | 2023-10-31 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103809178A (zh) * | 2014-01-17 | 2014-05-21 | 西安空间无线电技术研究所 | 一种地球同步轨道合成孔径雷达实现覆盖区连续观测方法 |
CN103954289A (zh) * | 2014-05-20 | 2014-07-30 | 哈尔滨工业大学 | 一种光学成像卫星敏捷机动姿态确定方法 |
CN105116910A (zh) * | 2015-09-21 | 2015-12-02 | 中国人民解放军国防科学技术大学 | 一种对地面点凝视成像的卫星姿态控制方法 |
CN105115477A (zh) * | 2015-07-27 | 2015-12-02 | 上海卫星工程研究所 | 对地面点目标推扫成像任务参数的星载求解方法 |
-
2016
- 2016-06-30 CN CN201610502752.7A patent/CN106197434A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103809178A (zh) * | 2014-01-17 | 2014-05-21 | 西安空间无线电技术研究所 | 一种地球同步轨道合成孔径雷达实现覆盖区连续观测方法 |
CN103954289A (zh) * | 2014-05-20 | 2014-07-30 | 哈尔滨工业大学 | 一种光学成像卫星敏捷机动姿态确定方法 |
CN105115477A (zh) * | 2015-07-27 | 2015-12-02 | 上海卫星工程研究所 | 对地面点目标推扫成像任务参数的星载求解方法 |
CN105116910A (zh) * | 2015-09-21 | 2015-12-02 | 中国人民解放军国防科学技术大学 | 一种对地面点凝视成像的卫星姿态控制方法 |
Non-Patent Citations (1)
Title |
---|
胡雅斯等: "适用多传感器的卫星对地覆盖计算模型", 《遥感技术与应用》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287352A (zh) * | 2018-02-06 | 2018-07-17 | 西安四方星途测控技术有限公司 | 一种卫星圆锥传感器作用范围的确认方法及系统 |
CN108287352B (zh) * | 2018-02-06 | 2021-11-02 | 中科星图(西安)测控技术有限公司 | 一种卫星圆锥传感器作用范围的确认方法及系统 |
CN108491633A (zh) * | 2018-03-24 | 2018-09-04 | 中国人民解放军海军航空大学 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
CN108491633B (zh) * | 2018-03-24 | 2021-10-01 | 中国人民解放军海军航空大学 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
CN109001762A (zh) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | 一种同步带目标光学观测卫星的地气光抑制方法及系统 |
CN110162069A (zh) * | 2019-05-10 | 2019-08-23 | 北京航空航天大学 | 一种近地轨道航天器阳光反射凝视期望姿态解析求解方法 |
CN111552272A (zh) * | 2020-05-07 | 2020-08-18 | 上海航天电子有限公司 | 无源相控阵天线波束控制算法 |
CN111751809B (zh) * | 2020-06-09 | 2023-11-14 | 军事科学院系统工程研究院后勤科学与技术研究所 | 一种点源靶标反射镜的调整角度计算方法 |
CN111751809A (zh) * | 2020-06-09 | 2020-10-09 | 军事科学院系统工程研究院后勤科学与技术研究所 | 一种点源靶标反射镜的调整角度计算方法 |
CN111966953A (zh) * | 2020-08-07 | 2020-11-20 | 四川泛华航空仪表电器有限公司 | 一种基于坐标转换的飞机姿态角计算方法 |
CN114282166A (zh) * | 2021-12-22 | 2022-04-05 | 中国西安卫星测控中心 | 基于环境函数矩阵的卫星经度计算方法、变轨判别方法 |
CN114282166B (zh) * | 2021-12-22 | 2024-04-19 | 中国西安卫星测控中心 | 基于环境函数矩阵的卫星经度计算方法、变轨判别方法 |
CN116975504A (zh) * | 2023-09-22 | 2023-10-31 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
CN116975504B (zh) * | 2023-09-22 | 2023-12-15 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106197434A (zh) | 基于地面目标点位置的卫星姿态角的计算方法 | |
CN106197425A (zh) | 基于卫星姿态角的地面目标点位置的计算方法 | |
CN104848860B (zh) | 一种敏捷卫星成像过程姿态机动规划方法 | |
CN106124170B (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
CN106155074B (zh) | 一种保证星地链路的卫星三轴对日定向控制方法 | |
CN104332707B (zh) | 一种用于低轨星载天线跟踪地面站的方法 | |
CN104570742B (zh) | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 | |
Strange et al. | Mapping the V-infinity Globe | |
CN104729457B (zh) | 太阳相对近地轨道微小卫星位置的确定方法 | |
CN103983254A (zh) | 一种新型敏捷卫星机动中成像方法 | |
JPH06510502A (ja) | 3軸安定衛星の姿勢制御で使用する測定装置と付属する評価方法、制御系および制御方法 | |
CN107607947A (zh) | 基于卡尔曼滤波的星载雷达成像参数在线估计方法 | |
CN101109959A (zh) | 一种适用于任意运动微小型系统的定姿系统 | |
CN106483466B (zh) | 一种卫星入轨阶段太阳电池阵输出电流的估算方法 | |
CN108710379A (zh) | 静止卫星成像偏航导引角计算方法 | |
CN106096204B (zh) | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 | |
CN107450582A (zh) | 一种基于星上实时规划的相控阵数传引导控制方法 | |
CN108181916A (zh) | 小卫星相对姿态的控制方法及装置 | |
CN106125745A (zh) | 一种对空间合作目标跟踪成像的卫星姿态控制方法 | |
CN107525492A (zh) | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 | |
CN103675760A (zh) | 一种星载地球同步轨道合成孔径雷达姿态导引方法 | |
Bangert et al. | Performance characteristics of the UWE-3 miniature attitude determination and control system | |
CN103487032A (zh) | 低轨空间相机自由指向像移矢量计算方法 | |
CN104567881A (zh) | 一种基于太阳、地心及月表重力矢量的巡视器位置姿态确定方法 | |
CN104714243B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161207 |
|
RJ01 | Rejection of invention patent application after publication |