[发明专利]卫星对目标观测下视角计算方法及系统有效
申请号: | 202110154808.5 | 申请日: | 2021-02-04 |
公开(公告)号: | CN112966211B | 公开(公告)日: | 2022-03-18 |
发明(设计)人: | 俞航;吕旺;范季夏;黄金生;黄欣;王田野;袁伟 | 申请(专利权)人: | 上海卫星工程研究所 |
主分类号: | G06F17/15 | 分类号: | G06F17/15;G06F17/16;G06F16/29 |
代理公司: | 上海段和段律师事务所 31334 | 代理人: | 李佳俊;郭国中 |
地址: | 200240 *** | 国省代码: | 上海;31 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 卫星 目标 观测 视角 计算方法 系统 | ||
1.一种卫星对目标观测下视角计算方法,其特征在于,包括:
步骤S102:计算卫星与目标点的距离矢量在地固坐标系下的分量A;
步骤S103:计算卫星与目标点的距离矢量在卫星天线坐标系下的分量B;
步骤S104:结合所述分量A与所述分量B,计算目标下视角角度;
所述步骤S102中:
在计算目标点的地理位置时,考虑地球椭球面模型,地球椭球面模型采用如下公式:
获得目标在地固系下的位置矢量,其中,r为目标点在地固系下的位置矢量,rx,ry,rz分别表示在地固系下的X方向、Y方向、Z方向的位置矢量,I为目标的地理经度,μ为目标的地理纬度,h为目标的水平地高,λs和rs由以下公式表示:
λs=arc tan((1-f)2tanμ)
其中,f为地球椭率,数值为R为地球赤道半径,数值为6378137米;
根据GPS数据读取T时刻后目标点在地固系下的位置rT:
对卫星和目标点进行矢量作差,求出卫星至目标点在地固坐标系下的距离矢量drd:
所述步骤S103中:
根据卫星轨道参数,时刻信息,计算地固系至卫星天线坐标系的转化矩阵CEF:
CEF=Ca,O*CO,I*INV(C′I,EF*CI,EF)*C′I,EF
其中CI,EF为J2000地心惯性坐标系到地固坐标系的变换矩阵;C′I,EF为CI,EF的逆矩阵;INV为矩阵的伪逆函数;CO,I为J2000地心惯性坐标系到轨道坐标系的变换矩阵;Ca,O为轨道坐标系到卫星天线坐标系的变换矩阵;
计算在卫星天线坐标系下的距离矢量drb:
所述步骤S104中:
计算卫星对目标观测的下视角θb,t0:
其中,ryb,t0、rzb,t0分别为t0秒的距离矢量在卫星天线坐标系下的Y、Z方向分量。
2.一种卫星对目标观测下视角计算系统,其特征在于,包括:
模块S102:计算卫星与目标点的距离矢量在地固坐标系下的分量A;
模块S103:计算卫星与目标点的距离矢量在卫星天线坐标系下的分量B;
模块S104:结合所述分量A与所述分量B,计算目标下视角角度;
所述模块S102中:
在计算目标点的地理位置时,考虑地球椭球面模型,地球椭球面模型采用如下公式:
获得目标在地固系下的位置矢量,其中,r为目标点在地固系下的位置矢量,rx,ry,rz分别表示在地固系下的X方向、Y方向、Z方向的位置矢量,I为目标的地理经度,μ为目标的地理纬度,h为目标的水平地高,λs和rs由以下公式表示:
λs=arc tan((1-f)2tanμ)
其中,f为地球椭率,数值为R为地球赤道半径,数值为6378137米;
根据GPS数据读取T时刻后目标点在地固系下的位置rT:
对卫星和目标点进行矢量作差,求出卫星至目标点在地固坐标系下的距离矢量drd:
所述模块S103中:
根据卫星轨道参数,时刻信息,计算地固系至卫星天线坐标系的转化矩阵CEF:
CEF=Ca,O*CO,I*INV(C′I,EF*CI,EF)*C′I,EF
其中CI,EF为J2000地心惯性坐标系到地固坐标系的变换矩阵;C′I,EF为CI,EF的逆矩阵;INV为矩阵的伪逆函数;CO,I为J2000地心惯性坐标系到轨道坐标系的变换矩阵;Ca,O为轨道坐标系到卫星天线坐标系的变换矩阵;
计算在卫星天线坐标系下的距离矢量drb:
所述模块S104中:
计算卫星对目标观测的下视角θb,t0:
其中,ryb,t0、rzb,t0分别为t0秒的距离矢量在卫星天线坐标系下的Y、Z方向分量。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于上海卫星工程研究所,未经上海卫星工程研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110154808.5/1.html,转载请声明来源钻瓜专利网。