CN105775170A - 太阳同步轨道卫星双视场星敏感器安装指向设计方法 - Google Patents

Abstract

本发明公开了一种太阳同步轨道卫星双视场星敏感器安装指向设计方法，基于卫星在轨运行期间杂散光与卫星本体坐标系间的矢量模型，分析出杂散光与卫星本体坐标系的矢量关系，通过建立杂散光边界曲线数学方程，解算出不受杂散光影响区域，通过边界曲线夹角关系，解算出星敏感器指向方程，并得出为抑制杂散光星敏感器遮光罩设计最大太阳光、地气光遮蔽角之间关系；同时提出不受杂光影响区域不满足双星敏安装时，在单一工况下至少有一个星敏可工作方案。本发明具有如下优点：有效解决太阳同步轨道航天器所安装双视场星敏感器为避免杂散光影响安装指向问题，同时为不同安装条件下双视场星敏感器遮光罩最大太阳光、地气光遮蔽角设计提供理论依据。

Description

太阳同步轨道卫星双视场星敏感器安装指向设计方法

技术领域

本发明涉及航天器姿态测量领域，具体涉及一种太阳同步轨道卫星双视场星敏感器安装指向设计方法。

背景技术

星敏感器作为航天器姿态测量器件，以恒星作为姿态测量参考基准，与其它卫星姿态敏感器件(太阳敏感器、陀螺、磁强计等)相比具有精度高、无漂移、功耗低，并且输出绝对姿态信息的优点，是当前应用最广泛的姿态敏感器。在卫星运行过程中，星敏感器有可能受到太阳光或地气光的干扰，使星敏感器的像面背景噪声提高，影响星点提取的准确性，甚至星点会被淹没在背景杂光中，导致星敏感器无法正常工作。因此，星敏感器的抗杂散光问题是目前星敏感器研究的一项重要内容。

在星敏感器杂散光抑制方面，一方面可以通过设计星敏感器光学镜头、遮光罩结构或改变遮光罩内壁涂层材料提高星敏感器自身的抗杂散光性能，另一方面可以通过分析和优化星敏感器在卫星上的安装位置来尽量避免杂散光对星敏感器的影响。

相比于单视场星敏感器，双视场星敏感器同时工作具有可测视场角大、精度高等优点，为保证双视场星敏同时正常工作，其在航天器上安装条件更为苛刻，同时当各视轴相互正交时，达到最高精度；即使在光照条件极为复杂的情况下，当其中一个星敏感器受杂光影响不能正常工作时，另一个星敏感器也可作为备份星敏感器进行工作。因此，双视场星敏感器在测量精度及可靠性方面具有明显优势，是星敏感器技术发展的重要方向。

在星敏感器安装指向方法研究方面，有许多研究机构针对实际应用提出了设计方法，其中，针对卫星多种工况姿态条件，2010年航天东方红公司提出了基于任务规划的星敏感器安装角度确定方法可以解决敏捷机动工况下给定星敏感器安装指向问题，但该方法通过不断调整星敏感器的安装角度和工况条件确定出合适的星敏感器安装角度，设计方法十分复杂，同时不能给出在不同安装指向条件下星敏感器遮光罩最大太阳光、地气光遮蔽角。

发明内容

本发明旨在至少解决上述技术问题之一。

为此，本发明的一个目的在于提出一种太阳同步轨道卫星双视场星敏感器安装指向设计方法。

为了实现上述目的，本发明的实施例公开了一种太阳同步轨道卫星双视场星敏感器安装指向设计方法，包括以下步骤：S1：设置太阳同步轨道参数、卫星参数、STK软件中HPOP阻力模型参数；S2：在STK软件中进行HPOP仿真，确定卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，其中，对于卫星本体坐标系，Z轴指向地心，X轴为卫星前进方向，Y轴由右手定则确定，α为太阳光与X轴夹角，β为太阳光与Y轴夹角，γ为太阳光与Z轴夹角；S3:根据地球半径、大气层高度计算确定地气光与卫星本体边界夹角；S4：根据所述卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，确定所述卫星在无机动状态下，太阳光相对卫星本体坐标系边界数学方程；根据所述地气光与卫星本体边界夹角确定所述卫星在无机动状态下，地气光相对卫星本体坐标系边界数学方程；S5：设置卫星机动状态下最大侧摆角和俯仰角，根据坐标系旋转矩阵确定所述卫星在不同机动状态下太阳光、地气光相对卫星本体坐标系边界曲线数学方程，以及各种工况条件下卫星本体坐标系不受杂光影响区域；S6：判断各种工况条件下不受杂光影响区域能否满足双星敏安装，如满足双星敏安装条件，根据方程所述卫星在各种工况机动状态下，太阳光相对卫星本体坐标系边界数学方程、地气光相对卫星本体坐标系边界数学方程，分析杂散光边界夹角关系，确定双视场星敏感器安装条件下为避免杂散光影响星敏感器安装指向范围数学方程，其中ω为双视场星敏感器单视场指向偏离YZ平面角度；S7：根据双视场星敏感器安装平面与卫星本体坐标系XY平面夹角与进动角ω关系方程确定双视场星敏感器安装平面与卫星本体坐标系XY平面夹角范围及其与进动角ω关系；S8：根据关系式计算太阳光最大遮蔽角；根据关系式极值点z<0和极值点z>0计算地气光最大遮蔽角，其中，S_max为太阳光最大遮蔽角，E_max为地气光最大遮蔽角，为双视场星敏感器安装平面与卫星本体坐标系XY平面夹角，θ为极值点所在指向Z平面投影线与Z轴夹角，同时，极值点为星敏不同指向太阳光/地气光边界曲线在其轴线上最高投影点；确定双视场星敏感器遮光罩最大太阳光遮蔽角与地气光遮蔽角关系，根据遮光罩设计原则判断安装区间是否满足设计要求；若满足设计要求，设计结束；若不满足设计要求，如可进行工况调整，执行步骤S5，重新设置卫星机动状态下最大侧摆角、俯仰角，如不可进行工况调整，执行步骤S9；以及S9：随卫星侧摆角增大，在各种工况条件下不受杂光影响区域如不能满足双星敏安装，双星敏不能在各种工况条件下不受杂散光影响而进行正常工作时，为保证在不同工况条件下至少有一个星敏可正常工作，需保证：Ⅰ在卫星左侧摆及正常工况情况下有一个星敏可正常工作，Ⅱ在卫星右侧摆及正常工况有一个星敏可正常工作。

根据本发明实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法，可以有效解决太阳同步轨道航天器所安装双视场星敏感器为避免杂散光影响安装指向问题，同时为不同安装指向条件下双视场星敏感器遮光罩最大太阳光、地气光遮蔽角设计提供理论依据。

另外，根据本发明上述实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法，还可以具有如下附加的技术特征：

进一步地，在步骤S1中，所述太阳同步轨道参数包括高度h、降交点地方时T、在轨运行周期t；所述卫星参数包括卫星质量M、卫星体积V；所述STK软件中HPOP阻力模型参数包括大气阻力C_d、太阳光压C_r、卫星面质比A/M和地磁力K_p，其中，所述卫星面质比A/M为1/4卫星表面积S与质量M的比值。

进一步地，步骤S3进一步包括：设置地球半径R_e，大气层高度d，根据γ_atm＝asin((R_e+d)/h)，计算确定地气光与卫星本体边界夹角γ_atm。

进一步地，在步骤S5中，所述卫星的机动状态包括左侧摆、右侧摆、前仰和后仰。

进一步地，步骤S6进一步包括：S601：所述双视场星敏感器相互垂直，联立方程与卫星前侧摆地气光边界曲线方程、左侧摆地气光边界曲线方程可求解绕Z轴不同偏转角度平面与地气光边界曲线交角；S602：联立方程与卫星右侧摆太阳光边界曲线方程可求解绕Z轴不同偏转角度平面与太阳光边界曲线交角；S603：所述双视场星敏感器的安装指向与YZ平面成η角，η角星敏安装线与Z轴夹角范围为：sinδ～sinε，其中，ε为星敏感器安装线与相应地气光边界曲线交点与卫星本体坐标系原点的连线与Z轴的夹角，δ为星敏感器安装线与相应太阳光边界曲线交点与卫星本体坐标系原点的连线与Z轴的夹角，安装指向范围曲线方程分别为：

进一步地，步骤S9进一步包括：在卫星左侧摆及正常工况下，星敏可安装指向区域包含：①卫星左侧摆地气光边界曲线到卫星正常工况太阳光左侧边界曲线之间，②卫星正常工况地气光边界曲线到卫星左侧摆太阳光边界曲线之间；

在卫星右侧摆及正常工况下，星敏可安装指向区域包含：③卫星正常工况地气光边界曲线到卫星右侧摆太阳光边界曲线之间，④卫星正常工况太阳光右侧边界曲线到卫星右侧摆地气光边界曲线之间；

以上两种星敏指向在卫星本体坐标系YZ平面内安装其遮光罩最大太阳光遮蔽角、最大地气光遮蔽角设计可达到最优；为使双星敏垂直安装，同时避免卫星前后摆动地气光对星敏的影响，其最优安装区域为①和④；同时随侧摆角的不断增大，各区域安装区域逐渐减小，当区域④不能满足星敏安装时，双星敏最优指向区域为①和③。

本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1是本发明一个实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法的流程图；

图2是本发明一个实施例的太阳光相对卫星本体坐标系示意图；

图3是本发明一个实施例的三年在轨运行期间太阳光与卫星本体坐标系关系模型；

图4是本发明一个实施例的卫星无机动状态下太阳光边界曲线与卫星本体坐标系关系图；

图5是本发明一个实施例的卫星无机动状态下地气光边界曲线与卫星本体坐标系关系图；

图6是本发明一个实施例的卫星不同工况下太阳光、地气光边界曲线在YZ平面投影；

图7是本发明一个实施例的卫星不同工况下太阳光、地气光边界曲线在XZ平面投影；

图8是本发明一个实施例的不同安装偏角情况下太阳光、地气光与Z轴夹角；

图9是本发明一个实施例的双视场星敏感器安装指向线与Z轴夹角；

图10是本发明一个实施例的双视场安装平面与XY平面的夹角与进动角关系；

图11是本发明一个实施例的最大太阳光、地气光遮蔽角计算示意图；

图12是本发明一个实施例的星敏感器指向与太阳光最大遮蔽角关系图；

图13是本发明一个实施例的双视场星敏感器遮光罩最大太阳光遮蔽角计算流程图；

图14是本发明一个实施例的极值点Z坐标小于0时星敏感器指向与地气光关系图；

图15是本发明一个实施例的极值点Z坐标大于0时星敏感器指向与地气光关系图；

图16是本发明一个实施例的双视场星敏感器遮光罩最大地气光遮蔽角计算流程图；

图17是本发明一个实施例的双视场星敏感器遮光罩最大太阳遮蔽角与最大地气光遮蔽角的关系图。

图18是本发明一个实施例太阳同步轨道高度535km、10:30轨道卫星在左右摆动及前后摆动均为45°情况下，卫星本体坐标系受杂散光影响及星敏安装指向示意图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。

在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。

在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

参照下面的描述和附图，将清楚本发明的实施例的这些和其他方面。在这些描述和附图中，具体公开了本发明的实施例中的一些特定实施方式，来表示实施本发明的实施例的原理的一些方式，但是应当理解，本发明的实施例的范围不受此限制。相反，本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。

以下结合附图描述根据本发明实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法。

请参考图1，一种太阳同步轨道卫星双视场星敏感器安装指向设计方法，包括以下步骤：

S1：设置太阳同步轨道参数：包括高度h、降交点地方时T、在轨运行周期t；卫星参数：包括卫星质量M、卫星体积V；STK(SatelliteToolKit)软件中HPOP(高精度轨道预报)阻力模型参数：包括大气阻力C_d、太阳光压C_r、卫星面质比A/M、地磁力K_p，其中卫星面质比A/M为1/4卫星表面积S与质量M的比值。

具体地，考虑卫星对地成像、所受摄动及阻力因素，本实施例设置卫星轨道高度h为700km；同时考虑对地观测成像，本实施例设置降交点地方时T为受光照较为复杂的10:30轨道；本实施例设置卫星运行周期t为3年。

本实施例设置卫星质量M为20.5kg，卫星体积V为324*360*400mm。

本实施例设置大气阻力C_d为2.2，太阳光压C_r为1.0，地磁力K_p为3.0，根据卫星面质比A/M为1/4卫星表面积S与质量M的比值，计算可得卫星面质比A/M≈0.01。

S2：在STK(SatelliteToolKit)中进行HPOP(高精度轨道预报)仿真，确定卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，其中，对于卫星本体坐标系，Z轴指向地心，X轴为卫星前进方向，Y轴由右手定则确定，α为太阳光与X轴夹角，β为太阳光与Y轴夹角，γ为太阳光与Z轴夹角。

具体地，本实施例中，对于卫星本体坐标系，Z轴指向地心，X轴为卫星前进方向，Y轴由右手定则确定，α为太阳光与X轴夹角，β为太阳光与Y轴夹角，γ为太阳光与Z轴夹角，如图2所示。根据步骤S1参数设置，本实施例通过STK软件得出卫星在轨运行3年间太阳光与卫星本体坐标系三轴矢量夹角，其中，太阳光与X轴向矢量夹角余弦最大值α_max为0.9946，最小值α_min为-0.9946，因此太阳光与X轴夹角范围为5.9620°～174.0470°；太阳光与Y轴向矢量夹角余弦最大值β_max为-0.0841，最小值β_min为-0.4537，因此太阳光与Y轴夹角范围为94.8270°～116.9830°；太阳光与Z轴向矢量夹角余弦最大值γ_max为0.4376，最小值γ_min为-0.9914，因此太阳光与Y轴夹角范围为64.0470°～172.4730°。将各时刻太阳光矢量方向在卫星本体坐标系中描述出来，如图3所示，光点表示太阳光相对于卫星本体坐标系原点的来光矢量方向。

S3:设置地球半径R_e，大气层高度d，根据γ_atm＝asin((R_e+d)/h)，计算确定地气光与卫星本体边界夹角γ_atm，即与Z轴正方向夹角。

具体地，由于地心引力作用，几乎全部气体集中在离地面100公里的高度范围内，其中75％的大气又集中在地面至10公里高度的对流层范围内。本实施例设定大气层高度d＝100km，地球半径R_e＝6378km。根据关系式γ_atm＝asin((R_e+d)/h)，将上述参数代入上式计算可得地气光与卫星本体边界夹角(即与Z轴正方向夹角)γ_atm为66.24°，考虑半长轴变化，本实施例采用67°(对应659.44km)冗余设计。

S4.根据步骤S2中卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，确定卫星无机动状态下太阳光相对卫星本体坐标系边界数学方程；根据步骤S3确定卫星无机动状态下地气光相对卫星本体坐标系边界数学方程。

具体地，卫星无机动工作状态下，太阳光边界方程通过步骤S2中太阳光与卫星本体坐标系各轴夹角极值确定，根据步骤S2分析，太阳光在卫星本体坐标系边界方程由两条曲线构成，分别如下公式，其中z<0.4376，θ₁＝94.8270°，θ₂＝116.9830°，如图4。

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ y={\mathrm{sin\&theta;}}_1\end{array}\right.---\left(1\right)$

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ y={\mathrm{sin\&theta;}}_2\end{array}\right.---\left(2\right)$

由于卫星本体坐标系Z轴始终指向地心，因此，卫星无机动状态下地气光边界可通过不同轨道高度地气光半锥角θ直接确定，地气光在卫星本体坐标系边界方程如下公式，其中θ＝67°，如图5。

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ z=cos\mathrm{\&theta;}\end{array}\right.---\left(3\right)$

S5：设置卫星机动状态下最大侧摆角、俯仰角，根据坐标系旋转矩阵确定卫星在不同机动状态下(左侧摆、右侧摆、前仰、后仰)太阳光、地气光相对卫星本体坐标系边界曲线数学方程，以及各种工况条件下卫星本体坐标系不受杂光影响区域。

具体地，本实施例定义卫星前摆为：卫星绕本体坐标系绕Y轴由+X向-Z轴向偏转，最大偏转角为-30°；卫星后摆为：卫星绕本体坐标系绕Y轴由+X向+Z轴向偏转，最大偏转角为30°。本实施例定义卫星左侧摆：卫星绕本体坐标系绕X轴由-Y向-Z轴向偏转，最大偏转角为30°；右侧摆：卫星绕本体坐标系绕X轴由-Y向+Z轴向偏转，最大偏转角为30°。

卫星左侧摆30°地气光边界曲线方程等价于步骤S4中地气光边界方程(3)绕X轴旋转-30°，经旋转矩阵变换后，如下式：

卫星左侧摆30°太阳光边界曲线方程等价于步骤S4中太阳光边界方程(1)、(2)绕X轴旋转-30°，只需分析其左侧边界曲线(1)的影响，经旋转矩阵变换后，如下式：

卫星右侧摆30°地气光边界曲线方程等价于步骤S4中地气光边界方程(3)绕X轴旋转30°，经旋转矩阵变换后，如下式：

卫星右侧摆30°太阳光边界曲线方程等价于步骤S4中太阳光边界方程(1)、(2)绕X轴旋转30°，只需分析其右侧边界曲线(2)的影响，经旋转矩阵变换后，如下式：

卫星前摆30°地气光边界曲线方程等价于步骤S4中地气光边界方程(3)绕Y轴旋转-30°，经旋转矩阵变换后，如下式：

卫星前摆30°太阳光边界曲线方程等价于步骤S4中太阳光边界方程(1)、(2)绕Y轴旋转-30°，太阳光相对卫星本体坐标系的影响区域在地气光覆盖范围内，不需分析。

卫星后摆30°地气光边界曲线方程等价于步骤S4中地气光边界方程(3)绕Y轴旋转30°，经旋转矩阵变换后，如下式：

卫星后摆30°太阳光边界曲线方程等价于步骤S4中太阳光边界方程(1)、(2)绕Y轴旋转-30°，太阳光相对卫星本体坐标系的影响区域在地气光覆盖范围内，不需分析。

根据步骤S4和S5上述部分对太阳光、地气光在各种工况下的边界方程的提取及变换，可得出在各种工况下卫星本体坐标系不受太阳光、地气光干扰的区域，边界曲线在卫星本体YZ平面投影如图6，可知不受干扰光影响区域共两部分，夹角范围分别为57.8°和26°。本实施例设置星敏感器视场角为15°，因此，57.8°范围区域满足双星敏安装条件，26°范围区域不满足星敏感器安装最小条件。边界曲线在XZ平面投影如图7，可知不受干扰区域由卫星左侧摆地气光边界曲线(6)、卫星右侧摆太阳光边界曲线(5)、卫星后侧摆地气光边界曲线(9)和卫星前侧摆地气光边界曲线(8)构成。

S6.根据方程与步骤S4、S5确定的边界方程，分析杂散光边界夹角关系，确定双视场星敏感器安装条件下为避免杂散光影响星敏感器安装指向范围数学方程，其中ω为双视场星敏感器单视场指向偏离YZ平面角度。

具体地，本实施例双视场星敏感器相互垂直，为使双视场星敏感器在卫星各种工况下不受杂散光影响可同时正常工作，双视场星敏感器单视场指向偏离YZ平面ω与球面交线方程为：

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ x=sin\mathrm{\&omega;}\end{array}\right.---\left(10\right)$

联立方程(10)与卫星前侧摆地气光边界曲线方程(8)、左侧摆地气光边界曲线方程(6)可求解绕Z轴不同偏转角度平面与地气光边界曲线交角；同时，联立方程(10)与卫星右侧摆太阳光边界曲线方程(5)可求解绕Z轴不同偏转角度平面与太阳光边界曲线交角，如图8。随着星敏感器指向绕Z轴偏离卫星本体坐标系YZ平面角度的增大，太阳光边界曲线与地气光边界曲线夹角逐渐减小。因此，星敏感器安装指向与YZ平面夹角η为45°，45°星敏安装线与Z轴夹角范围sinδ～sinε为：sin(-41.01°)～sin(2.696°)，如图9。45°双星敏安装指向范围曲线方程分别为：

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ x=\sqrt{2}/2\\ -0.656<z<0.047,y>0\end{array}\right.---\left(11\right)$

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ x=-\sqrt{2}/2\\ -0.656<z<0.047,y>0\end{array}\right.---\left(12\right)$

S7.根据双视场星敏感器安装平面与卫星本体坐标系XY平面夹角与进动角ω关系方程确定双视场星敏感器安装平面与卫星本体坐标系XY平面夹角范围及其与进动角ω关系。

具体地，根据双视场星敏感器安装平面与XY平面夹角与进动角ω关系方程：

计算可得：①双视场星敏感器安装面与XY平面夹角范围为-65°～3°，其中ω为星敏感器安装线与YZ平面夹角，本实施例ω为45°。②双视场星敏感器安装平面与XY平面夹角为0°时，星敏感器指向所在绕Z轴旋转平面与YZ平面夹角(进动角)为45°，随着安装平面角度的增大，与YZ平面夹角(进动角)由45°逐渐增大到67°，如图10。因此，随着安装平面角度的变化，星敏感器指向所在绕Z轴旋转平面YZ平面夹角不断变化。

S8：根据关系式计算太阳光最大遮蔽角；根据关系式(极值点z<0)和(极值点z>0)计算地气光最大遮蔽角，其中，S_max为太阳光最大遮蔽角，E_max为地气光最大遮蔽角，为双视场星敏感器安装平面与卫星本体坐标系XY平面夹角，θ为极值点所在指向Z平面投影线与Z轴夹角，同时，极值点为星敏不同指向太阳光/地气光边界曲线在其轴线上最高投影点。确定双视场星敏感器遮光罩最大太阳光遮蔽角与地气光遮蔽角关系，根据遮光罩设计原则判断安装区间是否满足设计要求；若满足设计要求，设计结束；若不满足设计要求，执行步骤S5，重新设置卫星机动状态下最大侧摆角、俯仰角。

具体地，为计算星敏感器在安装线内不同指向遮光罩设计最大太阳光、地气光遮蔽角，需分析太阳光、地气光边界曲线在星敏感器指向所在绕Z轴旋转平面的投影问题，星敏感器指向不同，其投影平面不同，同时太阳光、地气光边界曲线在指向轴平面上的投影高度不同，而最大太阳光、地气光遮蔽角可通过安装指向线与相应最高点与原点连线的角度求得，如图11。

最大太阳遮蔽角可由双星敏安装平面与XY平面夹角太阳光极值点投影线与Z轴夹角θ求得，太阳最大遮蔽角其关系如图12。双视场星敏感器遮光罩最大太阳遮蔽角计算流程如图13。

与太阳光最大遮蔽角计算类似，地气光最大遮蔽角也由上述关系解算得出，但地气光最大遮蔽角的计算需对极值点Z坐标数据做出判别。当极值点Z坐标小于0时，地气光最大遮蔽角如图14；当极值点Z坐标大于0时，地气光最大遮蔽角如图15。双视场星敏感器遮光罩最大地气光遮蔽角计算流程如图16。

由上计算可得在双视场星敏感器指向平面与XY平面不同夹角条件下，遮光罩最大太阳遮蔽角与最大地气光遮蔽角的关系，如图17。

由于太阳光强远高于地气光，太阳遮蔽角需大于地气光遮蔽角，因此双星敏安装指向平面与XY平面夹角应大于-32°，同时，本实施例星敏感器视场角为15°，设计地气光遮蔽角不小于10°，由图17可知，双星敏安装指向平面与XY平面夹角应小于-11°。具体星敏感器安装平面指向需与星敏感器遮光罩设计要求相结合。若满足遮光罩尺寸及遮光性能设计要求，设计结束；若不满足遮光罩尺寸及遮光性能设计要求，执行步骤S5，适当减小卫星机动状态下最大侧摆角、俯仰角。

根据本发明实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法，可以有效解决太阳同步轨道航天器所安装双视场星敏感器为避免杂散光(太阳光、地气光)影响安装指向问题，同时为不同安装指向条件下双视场星敏感器遮光罩最大太阳光、地气光遮蔽角设计提供理论依据。为星敏感器避免杂散光影响进行安装指向及设计提供了一种有实际应用价值的设计方案。

S9：随卫星侧摆角增大，卫星各种工况条件下不受杂散光影响区域逐渐减小，当该区域不满足双星敏安装条件时，双星敏不能在各种工况条件下不受杂散光影响而进行正常工作。为保证在不同工况条件下至少有一个星敏可正常工作，需对星敏指向进行优化。

卫星正常工况受杂散光影响区域处于卫星左侧摆及右侧摆之间，当不满足任何工况双星敏可同时工作条件时，为使星敏使用效率最大化，需保证：Ⅰ在卫星左侧摆及正常工况情况下有一个星敏可正常工作，Ⅱ在卫星右侧摆及正常工况有一个星敏可正常工作。

在卫星左侧摆及正常工况下，星敏可安装指向区域包含：①卫星左侧摆地气光边界曲线到卫星正常工况太阳光左侧边界曲线之间，②卫星正常工况地气光边界曲线到卫星左侧摆太阳光边界曲线之间。

在卫星右侧摆及正常工况下，星敏可安装指向区域包含：③卫星正常工况地气光边界曲线到卫星右侧摆太阳光边界曲线之间，④卫星正常工况太阳光右侧边界曲线到卫星右侧摆地气光边界曲线之间。

以上两种星敏指向在卫星本体坐标系YZ平面内安装其遮光罩最大太阳光遮蔽角、最大地气光遮蔽角设计可达到最优。为使双星敏垂直安装，同时避免卫星前后摆动地气光对星敏的影响，其最优安装区域为①和④。同时随侧摆角的不断增大，各区域安装区域逐渐减小，当区域④不能满足星敏安装时，双星敏最优指向区域为①和③。

图18为太阳同步轨道高度535km、10:30轨道卫星在左右摆动及前后摆动均为45°情况下，卫星本体坐标系受杂散光影响示意图。在卫星左侧摆及正常工况下，星敏安装指向区域包含在YZ平面内区域①及区域②；在卫星右侧摆及正常工况下，星敏安装指向区域包含在YZ平面内区域③及区域④。为使双星敏垂直安装，同时避免卫星前后摆动地气光对星敏的影响，其最优安装区域为①和④。同时随侧摆角的不断增大，各区域安装区域逐渐减小，当区域④不能满足星敏安装时，双星敏最优指向区域为①和③。

另外，本发明实施例的太阳同步轨道卫星双视场星敏感器安装指向设计方法的其它构成以及作用对于本领域的技术人员而言都是已知的，为了减少冗余，不做赘述。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同限定。

Claims (6)

1.一种太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，包括以下步骤：

S1：设置太阳同步轨道参数、卫星参数、STK软件中HPOP阻力模型参数；

S2：在STK软件中进行HPOP仿真，确定卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，其中，对于卫星本体坐标系，Z轴指向地心，X轴为卫星前进方向，Y轴由右手定则确定，α为太阳光与X轴夹角，β为太阳光与Y轴夹角，γ为太阳光与Z轴夹角；

S3:根据地球半径、大气层高度计算确定地气光与卫星本体边界夹角；

S4：根据所述卫星在轨运行期间太阳光与卫星本体坐标系三轴矢量夹角，确定所述卫星在无机动状态下，太阳光相对卫星本体坐标系边界数学方程；根据所述地气光与卫星本体边界夹角确定所述卫星在无机动状态下，地气光相对卫星本体坐标系边界数学方程；

S5：设置卫星机动状态下最大侧摆角和俯仰角，根据坐标系旋转矩阵确定所述卫星在不同机动状态下太阳光、地气光相对卫星本体坐标系边界曲线数学方程，以及各种工况条件下卫星本体坐标系不受杂光影响区域；

S6：判断各种工况条件下不受杂光影响区域能否满足双星敏安装，如满足双星敏安装条件，根据方程所述卫星在各种工况机动状态下，太阳光相对卫星本体坐标系边界数学方程、地气光相对卫星本体坐标系边界数学方程，分析杂散光边界夹角关系，确定双视场星敏感器安装条件下为避免杂散光影响星敏感器安装指向范围数学方程，其中ω为双视场星敏感器单视场指向偏离YZ平面角度；

S7：根据双视场星敏感器安装平面与卫星本体坐标系XY平面夹角与进动角ω关系方程确定双视场星敏感器安装平面与卫星本体坐标系XY平面夹角范围及其与进动角ω关系；

S8：根据关系式计算太阳光最大遮蔽角；根据关系式极值点z<0和极值点z>0计算地气光最大遮蔽角，其中，S_max为太阳光最大遮蔽角，E_max为地气光最大遮蔽角，为双视场星敏感器安装平面与卫星本体坐标系XY平面夹角，θ为极值点所在指向Z平面投影线与Z轴夹角，同时，极值点为星敏不同指向太阳光/地气光边界曲线在其轴线上最高投影点；确定双视场星敏感器遮光罩最大太阳光遮蔽角与地气光遮蔽角关系，根据遮光罩设计原则判断安装区间是否满足设计要求；若满足设计要求，设计结束；若不满足设计要求，如可进行工况调整，执行步骤S5，重新设置卫星机动状态下最大侧摆角、俯仰角；如不可进行工况调整，执行步骤S9；以及

S9：随卫星侧摆角增大，在各种工况条件下不受杂光影响区域如不能满足双星敏安装，双星敏不能在各种工况条件下不受杂散光影响而进行正常工作时，为保证在不同工况条件下至少有一个星敏可正常工作，需保证：Ⅰ在卫星左侧摆及正常工况情况下有一个星敏可正常工作，Ⅱ在卫星右侧摆及正常工况有一个星敏可正常工作。

2.根据权利要求1所述的太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，在步骤S1中，所述太阳同步轨道参数包括高度h、降交点地方时T、在轨运行周期t；

所述卫星参数包括卫星质量M、卫星体积V；

所述STK软件中HPOP阻力模型参数包括大气阻力C_d、太阳光压C_r、卫星面质比A/M和地磁力K_p，其中，所述卫星面质比A/M为1/4卫星表面积S与质量M的比值。

3.根据权利要求2所述的太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，步骤S3进一步包括：

设置地球半径R_e，大气层高度d，根据γ_atm＝asin((R_e+d)/h)，计算确定地气光与卫星本体边界夹角γ_atm。

4.根据权利要求3所述的太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，在步骤S5中，所述卫星的机动状态包括左侧摆、右侧摆、前仰和后仰。

5.根据权利要求4所述的太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，步骤S6进一步包括：

S601：所述双视场星敏感器相互垂直，联立方程与卫星前侧摆地气光边界曲线方程、左侧摆地气光边界曲线方程可求解绕Z轴不同偏转角度平面与地气光边界曲线交角；

S602：联立方程与卫星右侧摆太阳光边界曲线方程可求解绕Z轴不同偏转角度平面与太阳光边界曲线交角；

S603：所述双视场星敏感器的安装指向与YZ平面成η角，η角星敏安装线与Z轴夹角范围为：sinδ～sinε，其中，ε为星敏感器安装线与相应地气光边界曲线交点与卫星本体坐标系原点的连线与Z轴的夹角，δ为星敏感器安装线与相应太阳光边界曲线交点与卫星本体坐标系原点的连线与Z轴的夹角，安装指向范围曲线方程分别为：

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ x=sin\mathrm{\&eta;}\\ sin\mathrm{\&delta;}<z<sin\mathrm{\&epsiv;},y>0\end{array}\right.$

$\left\{\begin{array}{c}{x}^2+y^2+z^2=1\\ x=-sin\mathrm{\&eta;}\\ sin\mathrm{\&delta;}<z<sin\mathrm{\&epsiv;},y>0\end{array}\right..$

6.根据权利要求5所述的太阳同步轨道卫星双视场星敏感器安装指向设计方法，其特征在于，步骤S9进一步包括：

在卫星左侧摆及正常工况下，星敏可安装指向区域包含：①卫星左侧摆地气光边界曲线到卫星正常工况太阳光左侧边界曲线之间，②卫星正常工况地气光边界曲线到卫星左侧摆太阳光边界曲线之间；

在卫星右侧摆及正常工况下，星敏可安装指向区域包含：③卫星正常工况地气光边界曲线到卫星右侧摆太阳光边界曲线之间，④卫星正常工况太阳光右侧边界曲线到卫星右侧摆地气光边界曲线之间；

以上两种星敏指向在卫星本体坐标系YZ平面内安装其遮光罩最大太阳光遮蔽角、最大地气光遮蔽角设计可达到最优；为使双星敏垂直安装，同时避免卫星前后摆动地气光对星敏的影响，其最优安装区域为①和④；同时随侧摆角的不断增大，各区域安装区域逐渐减小，当区域④不能满足星敏安装时，双星敏最优指向区域为①和③。