一种可变步长掩星预报方法技术

一种可变步长掩星预报方法，先定义“卫星-目标源”坐标系并建立其相对于WGS-84坐标系的转换矩阵；然后建立卫星在“卫星-目标源”坐标系中掩星事件的判据，若满足判据，进一步判断卫星与目标源连线的矢量是否在卫星掩星天线可观测范围内，若在，则认为发生一次掩星事件，进而对掩星事件进行有效预报，本发明专利技术采用坐标变换的思想，将掩星发生过程变换为“卫星-目标源”平面内定点转动射线切割圆周的几何问题，从本质上揭示了掩星事件发生的物理机理，是一种可变步长的掩星预报算法，该方法可有效地提高计算效率，且与卫星掩星天线安装相结合，适合实际工程应用，尤其适用于长期地对太阳或GPS掩星事件进行预报。

【技术实现步骤摘要】
一种可变步长掩星预报方法
本专利技术涉及一种掩星预报方法，特别是一种可变步长掩星预报方法，该方法基于坐标变换的思想获得掩星事件的判据，且结合卫星掩星天线的相关信息，适用于长期地对太阳掩星或GPS掩星事件进行预报。
技术介绍
目前，气候变化问题已成为全世界共同发展的焦点问题。大气污染是世界各国面临的最大挑战之一，各国政府非常重视，尤其是温室气体的排放。因此，开展大气探测技术有助于进一步了解人类在温室气体排放和导致全球气候变化方面所扮演的角色。全球变暖危及人类的可持续发展，需要大力减少温室气体排放，减排涉及能源安全、经济发展和生活水平，具有现实代价。对大气成分的观测模式主要包括：天底模式、闪耀模式、目标模式、临边模式和掩星观测模式。根据卫星观测目标源的不同，掩星探测可分为太阳掩星(目标源为太阳)和无线电掩星(目标源为GPS或LEO卫星)。基于掩星观测模式的太阳掩星，可以有效地利用卫星获取的对地观测数据，获取全球温室气体分布与变化信息，实时监测大气温室气体变化，估算全球碳源汇区域格局，为减缓和适应全球变化提供数据和理论依据。GPS掩星大气探测技术是数值天气预报、大气空间环境探测的重要手段，可用于获取常规技术和手段不能达到的敏感地区的气象情报资料，为现代主战武器系统、侦察系统和通信系统提供高时间分辨率、高时效性、高精度的气象水文保障。发展掩星大气探测技术并开展其仿真研究，模拟计算电波在大气中的传播路径，有利于我们了解掩星观测的物理过程、分析大气中的物理过程对掩星的影响、掌握掩星观测量的特征和变化规律等，而对掩星事件的有效预报直接影响着大气成份探测的结果。现有针对掩星事件进行预报的方法通常以距离或夹角等为判据，通过采取定步长轨道递推的方式对掩星事件(发生时间和持续长度等)进行预报，这种方法类似于权举计算，所涉及的计算量较大，计算效率较低，尤其不适合针对掩星事件进行长期的预报。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提供了一种可变步长掩星预报方法，采用坐标变换的方法，将掩星的发生过程变换为“卫星-地球-目标源”平面内定点转动射线切割圆周的几何问题，进而从本质上揭示了掩星事件发生的物理机理，基于该思路得到了可变步长的掩星事件预报算法，可有效地提高计算效率，且该掩星预报方法与卫星掩星天线相关信息结合，适合实际工程应用。本专利技术的技术解决方案是：一种可变步长掩星预报方法，步骤如下：(1)建立“卫星-目标源”坐标系，即Sst坐标系，所述Sst坐标系具体为：以地球中心为坐标原点，X轴为从地心指向卫星的矢量RSat，卫星至目标源连线的矢量在XY平面内，Z轴由RSat与RSat-Taget的叉乘确定，即Z轴的单位矢量为Y轴方向由X轴方向和Z轴方向的单位矢量叉乘确定，即Y轴的单位矢量为(2)计算步骤(1)中建立的Sst坐标系与WGS-84坐标系之间的转换矩阵；并求得卫星至目标源连线在“卫星-目标源”坐标系中的单位向量；所述转换矩阵由公式：给出，其中R为转换矩阵，其中由公式：给出，式中，λ和分别为卫星在WGS-84坐标系中的地理经度和地理纬度；由公式：给出，其中，由公式：给出，其中，λT和分别为WGS-84坐标系中的目标源的经度和纬度；由公式：给出；卫星至目标源连线在“卫星-目标源”坐标系中的单位向量由公式：给出；(3)计算步骤(2)中得到的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角，并利用卫星在Sst坐标系中掩星事件判据，判断是否发生掩星事件；若满足掩星事件判据，则计算满足掩星事件判据的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量，进入步骤(4)；并利用中心切点判据求得发生掩星事件的中心切点；(4)利用步骤(2)得到的Sst坐标系与WGS-84坐标系之间的转换矩阵将步骤(3)中满足掩星事件判据的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量转换到WGS-84坐标系中，然后利用WGS-84坐标系与卫星本体系之间的转换矩阵将卫星至目标源连线的单位向量转换到卫星本体系中，并判断是否满足掩星事件有效判据，若满足，则发生了一次掩星事件，否则没有发生掩星事件。所述步骤(3)中的卫星在Sst坐标系中掩星事件判据具体为：asin((hmin+Re)/a)≤π-α≤asin((hmax+Re)/a)式中，α为卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角，a为卫星的轨道半长轴，hmin、hmax和hmean分别为发生掩星事件过程中对应的大气层的最小高度、最大高度和中心切点的大气高度，hmin、hmax和hmean均为预先给定的固定值；Re为地球半径。所述步骤(3)中的中心切点判据为：π-α′＝asin((hmean+Re)/a)由上述判据计算得到求得中心切点时对应的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角α′，所述中心切点为Sst坐标系XY平面内卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角为α′时，通过地心且与卫星至目标源连线垂直的直线与卫星至目标源连线的交点。所述步骤(4)中的Sst坐标系与WGS-84坐标系之间的转换矩阵，具体由公式：Reo＝R0·Rx(-π/2)·Rz(u+π/2)·Rx(i)·Rz(Ω-αG)给出，其中R0为卫星在1-2-3转序下三轴姿态机动角度分别为θ和φ情况下的姿态矩阵，Rx、Ry和Rz分别为卫星本体系中绕X轴、Y轴和Z轴的基元变换矩阵；i、Ω和u分别为卫星运行轨道的轨道倾角、升交点赤经和纬度幅角；αG为格林威治赤经。所述步骤(4)中的掩星事件有效判据具体为：判断矢量和Rantenna的夹角是否小于掩星天线半张角β，若满足掩星事件有效判据，则发生一次有效掩星事件，进而对该次掩星事件进行有效预报；否则没有发生有效掩星事件，不进行掩星事件有效预报，所述矢量为卫星至目标源连线的在卫星本体系中的单位矢量，Rantenna为卫星掩星天线在卫星本体系中的安装方位。本专利技术与现有技术相比的有益效果是：(1)本专利技术将掩星事件发生过程变换为“卫星-目标源”坐标系中定点转动射线切割圆周的几何问题，从本质上揭示了掩星事件发生的物理机理，原理更加直观，且只需对一个夹角进行掩星事件判断，计算更加简单；(2)本专利技术中的掩星事件判据规定了卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角的有效范围，在实际掩星事件预报中，只需对处于有效范围内的夹角进行判断，在有效范围外采用可变步长的方法进行轨道递推，简化了掩星事件预报的计算，提高了掩星事件预报的效率；(3)本专利技术在卫星坐标系中进行掩星事件有效判据的判断，综合考虑了卫星实际的在轨姿态和掩星天线的安装方位，有效判据所需信息获取简单，贴近实际工程应用。附图说明图1为“卫星-目标源”坐标系中掩星事件几何示意图；图2为本专利技术流程图图2；图3为低轨卫星针对各颗GPS卫星的发生电离层掩星事件次数分布示意图；图4为低轨卫星电离层掩星发生时刻随时间的分布示意图；图5为中心切点在WGS-84坐标系下经纬度分布示意图；图6为某颗GPS卫星一天内每次掩星的星下点地理经纬度分布示意图；图7为某颗GPS卫星掩星时长分本文档来自技高网...

【技术保护点】
一种可变步长掩星预报方法，其特征在于步骤如下：(1)建立“卫星‑目标源”坐标系，即Sst坐标系，所述Sst坐标系具体为：以地球中心为坐标原点，X轴为从地心指向卫星的矢量RSat，卫星至目标源连线的矢量在XY平面内，Z轴由RSat与RSat‑Taget的叉乘确定，即Z轴的单位矢量为Y轴方向由X轴方向和Z轴方向的单位矢量叉乘确定，即Y轴的单位矢量为(2)计算步骤(1)中建立的Sst坐标系与WGS‑84坐标系之间的转换矩阵；并求得卫星至目标源连线在“卫星‑目标源”坐标系中的单位向量；所述转换矩阵由公式：R=[x→,y→,z→]T]]>给出，其中R为转换矩阵，其中由公式：给出，式中，λ和分别为卫星在WGS‑84坐标系中的地理经度和地理纬度；由公式：z→=x→×R→Sat-Taget]]>给出，其中，由公式：给出，其中，λT和分别为WGS‑84坐标系中的目标源的经度和纬度；由公式：y→=z→×x→]]>给出；卫星至目标源连线在“卫星‑目标源”坐标系中的单位向量由公式：R→ss_new=R·R→Sat-Taget]]>给出；(3)计算步骤(2)中得到的卫星至目标源连线在“卫星‑目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角，并利用卫星在Sst坐标系中掩星事件判据，判断是否发生掩星事件；若满足掩星事件判据，则计算满足掩星事件判据的卫星至目标源连线在“卫星‑目标源”坐标系中的单位向量，进入步骤(4)；并利用中心切点判据求得发生掩星事件的中心切点；(4)利用步骤(2)得到的Sst坐标系与WGS‑84坐标系之间的转换矩阵将步骤(3)中满足掩星事件判据的卫星至目标源连线在“卫星‑目标源”坐标系中的单位向量转换到WGS‑84坐标系中，然后利用WGS‑84坐标系与卫星本体系之间的转换矩阵将卫星至目标源连线的单位向量转换到卫星本体系中，并判断是否满足掩星事件有效判据，若满足，则发生了一次掩星事件，否则没有发生掩星事件。...

【技术特征摘要】
1.一种可变步长掩星预报方法，其特征在于步骤如下：(1)建立“卫星-目标源”坐标系，即Sst坐标系，所述Sst坐标系具体为：以地球中心为坐标原点，X轴为从地心指向卫星的矢量RSat，卫星至目标源连线的矢量在XY平面内，Z轴由RSat与RSat-Taget的叉乘确定，即Z轴的单位矢量为Y轴方向由X轴方向和Z轴方向的单位矢量叉乘确定，即Y轴的单位矢量为(2)计算步骤(1)中建立的Sst坐标系与WGS-84坐标系之间的转换矩阵；并求得卫星至目标源连线在“卫星-目标源”坐标系中的单位向量；所述转换矩阵由公式：给出，其中R为转换矩阵，其中由公式：给出，式中，λ和分别为卫星在WGS-84坐标系中的地理经度和地理纬度；由公式：给出，其中，由公式：给出，其中，λT和分别为WGS-84坐标系中的目标源的经度和纬度；由公式：给出；卫星至目标源连线在“卫星-目标源”坐标系中的单位向量由公式：给出；(3)计算步骤(2)中得到的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量与Sst坐标系X轴之间的夹角，并利用卫星在Sst坐标系中掩星事件判据，判断是否发生掩星事件；若满足掩星事件判据，则计算满足掩星事件判据的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量，进入步骤(4)；并利用中心切点判据求得发生掩星事件的中心切点；(4)利用步骤(2)得到的Sst坐标系与WGS-84坐标系之间的转换矩阵将步骤(3)中满足掩星事件判据的卫星至目标源连线在“卫星-目标源”坐标系中的单位向量转换到WGS-84坐标系中，然后利用WGS-84坐标系与卫星本体系之间的转换矩阵将卫星至目标源连线的单位向量转换到...

【专利技术属性】
技术研发人员：吕秋杰，徐明，谭田，蒙薇，张燕，李志武，陈罗婧，
申请(专利权)人：航天东方红卫星有限公司，
类型：发明
国别省市：北京;11