具体实施方式

下面结合实施例对本发明作进一步阐述。

火星上层大气所产生的阻力影响火星探测器轨道，使得环绕火星飞行的探测器轨道半长轴减小，影响卫星寿命。同时，大气阻力是轨道计算过程中的一项重要摄动力，定量描述探测器受到的大气阻力有利于提高火星探测器轨道计算和轨道预报的精度。我国目前的火星探测计划中，由于探测器绕火星近距离飞行(500km高度以下)的时间较短，所以轨道计算和预报中没有考虑火星大气阻力的影响。随着后续对火星的深入探测，必将开展近距离、长时间环绕火星飞行探测任务，届时，火星大气阻力将成为不可忽略的影响因素。

本发明提供一种计算火星探测器大气阻力的方法，能够有效提高火星探测器轨道计算和预报的精度，能够解决火星环绕探测器寿命预报问题。本方法操作方便，准确性高，可用于实时和准实时轨道计算。

火星探测器大气阻力的计算流程如图1所示，具体包括如下步骤：

步骤一、计算火星大气的阻力系数C_D；大气阻力系数CD是表征卫星表面与碰撞上来的粒子之间动量交换的参数，它的确定基于两项假设，一是大气自由分子流假设。火星大气稀薄，对于特征长度在10米以内的探测器，100km以上区域均可视为自由分子流；二是大气粒子的散射假设。试验表明，粒子反射性质与物体表面的温度有关，如果物体表面温度接近室温为散射，而如果温度达到1600K的高温时，则几乎是镜面反射。火星大气温度在120～220K范围内，所以大气粒子的散射假设是成立的。因此火星大气的阻力系数C_D的计算方法为：

式中，T_si为探测器各表面温度；

T_a为大气的温度；

β_i为大气的运动方向与探测器各表面的夹角，通过卫星探测器姿态监视系统获得；

μ_i为入射气体分子的平均质量与卫星探测器散射表面分子质量比。

步骤二、根据探测器所处的位置，综合考虑太阳活动，火星季节，火星地方时和火星尘暴情况，通过火星气候数据库MCD_V5.3获得火星大气密度ρ；火星气候数据库MCD_V5.3是从火星大气循环模型数值模拟获得的气象领域数据库，并使用观测数据完成验证；火星气候数据库MCD_V5.3提供了0～260公里范围内的温度、压强、风场和大气密度，并根据太阳活动的程度设置了平静、中等和活跃三个场景。

通过火星气候数据库MCD_V5.3获得火星大气密度ρ的具体方法为：

调用火星气候数据库MCD_V5.3的子程序call_mcd，输入参数为预设纵坐标类型zkey、预设火星地面点纵坐标xz、预设火星地面点经度xlon、预设火星地面点纬度xlat、预设数据分辨率hireskey、预设历元表示方式datekey、预设历元时刻xdate、火星地方时localtime、数据库文件存放路径dset和预设的火星尘暴和太阳EUV辐射的场景scena，输出结果即为火星大气密度ρ。

步骤三、计算卫星探测器相对火星大气的速度v_rel；具体为：

S31、建立火星固联坐标系oxyz；其中，原点o为火星质心；选取火星平赤道面作为第一参考平面；x轴指向本初子午线与第一参考平面的焦点；z轴指向火星自转轴；y轴由右手定则确定。

S32、建立火星天球坐标系o1x1y1z1；火星天球坐标系是将地球天球坐标系平移至火星，原点o1为火星质心；第二参考平面为地球历元平赤道；x1轴指向历元平春分点方向；z轴指向地球自转轴；y轴由右手定则确定。

S33、计算火星固联坐标系oxyz到火星天球坐标系o1x1y1z1的各轴的单位旋转矩阵R_i(α)＝[R_x(α)、R_y(α)、R_z(α)]；其中，设定x轴到x1轴的旋转角度、y轴到y1轴的旋转角度、z轴到z1轴的旋转角度均为α；旋转矩阵R_i(α)的计算方法为：

S34、火星自转周期为T＝24小时37分22.6秒，由此可计算火星固联坐标系下火星大气的自转角速度为由于探测器的速度一般在火星天球系中描述，还需要将Ω转换到天球系，计算t时刻火星固联坐标系oxyz到火星天球坐标系o1x1y1z1的坐标转换矩阵RT2C；坐标转换矩阵RT2C的计算方法为：

式中，N为指向参数，N＝3.37919183°；

J为指向参数，J＝24.67682669°；

为火星章动；计算公式为

I为火星轨道倾角；

ψ₀和I₀为J2000.0历元时刻的常量，为火星章动变率，为火星轨道倾角的长期变化率，在表1中给出；ψ_nut和I_nut分别为经度和倾角章动修正，由下式计算：

式中，α_m和θ_m的计算方法为：

其中，n′表示火星的平均公转角速度，l′₀为J2000.0火星平近点角，取值为19.3871°，l′为火星平近点角,q＝2ω，ω为近星点幅角，q也可表示为q＝142°.0+1°.3t，t表示J2000.0起算的儒略世纪数，表1～2给出了相关系数。

表1

表2

Φ为火星自转角；

X_p和Y_p为火星自转轴的极移；

S35、计算火星天球坐标系下的火星大气速度Ω_in、卫星探测器的位置矢量r_in、卫星探测器的速度矢量v_in；计算卫星探测器相对火星大气的速度v_rel。星大气速度Ω_in的计算方法为：

卫星探测器在火星天球坐标系下的位置矢量r_in测量获得；

卫星探测器在火星天球坐标系下的速度矢量v_in测量获得；

卫星探测器相对火星大气的速度v_rel的计算方法为：

步骤四、计算卫星探测器的迎风面积A；卫星探测器的迎风面积A的计算方法为：

式中，A_s为卫星探测器本体迎风面积，在地测量获得；

A_PN为帆板迎风面积；

A_P为帆板总面积；

N为帆板法向量，根据太阳和探测器位置测量获得。

步骤五、计算卫星探测器受到的大气阻力f：

根据本发明设计的火星探测器大气阻力的方法，能够有效提高火星探测器轨道计算和预报的精度，能够解决火星环绕探测器寿命预报问题。本方法操作方便，准确性高，可用于实时和准实时轨道计算。

本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。