敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法与流程

本发明涉及一种敏捷合成孔径雷达(以下简称为sar)卫星聚束模式姿态机动需求计算方法，尤其适用于天线相位中心偏离星本体，并采用集中馈电的反射面体制进行敏捷sar卫星聚束模式姿态机动需求计算.
背景技术：
：聚束模式是一种常用的高分辨率星载sar工作模式，通过雷达波束扫描使波束中心始终瞄准地面固定点，从而使成像区域始终处于雷达波束的照射之下，延长了合成孔径时间，实现了方位向高分辨率成像。敏捷sar卫星有效载荷固定在卫星上，依靠姿态控制系统控制卫星整体绕俯仰、横滚、偏航3个轴向摆动，实现sar成像所需的波束扫描。敏捷sar卫星具有较高的灵活性和成像能力，能高效地实现高分辨率聚束模式成像。星载sar与光学成像机理不同，例如，sar需侧视成像、sar天线具备波束扫描能力以及sar对相位误差敏感等，导致光学遥感卫星中的姿态机动指标参数计算方法不适用于sar卫星。上海航天控制技术研究所的汪礼成等人在2013年高分辨率对地观测学术年会上发表了《对地聚束成像sar卫星姿态控制技术研究》一文，对聚束模式下sar卫星的姿态跟踪机动规律进行了研究，但未考虑雷达天线相位中心偏离星本体坐标系原点，和雷达天线波束指向偏离星本体坐标系坐标轴的情况。题为“一种基于sar卫星姿态控制实现滑动聚束模式的方法”的发明专利(北京航空航天大学，陈杰等，申请号为cn103076607)提出了一种基于sar卫星姿态控制实现滑动聚束模式的方法，但这种方法是针对滑动聚束模式而非聚束模式设计的，同时它也不能适应雷达天线相位中心偏离星本体坐标系原点，和雷达天线波束指向偏离星本体坐标系坐标轴的情况，不适用于高分辨率聚束模式敏捷sar卫星。因而，敏捷sar卫星聚束模式分辨率高，对雷达波束指向的精度要求提高，需要在高精度卫星轨道和地球模型的基础上，充分考虑雷达天线相位中心偏离星本体坐标系原点，和雷达天线波束指向偏离星本体坐标系坐标轴的影响，设计精确的敏捷sar卫星聚束模式姿态需求计算方案。技术实现要素：为了解决现有技术中存在的问题以及克服现有技术的不足，本发明提出了一种基于迭代 计算的敏捷sar卫星聚束模式姿态机动指标参数计算方案，能适应sar天线在星体上各种不同的安装位置，并充分考虑可能存在的波束离轴角，得到准确的敏捷sar卫星聚束模式成像所需的姿态指标参数，为基于卫星平台姿态机动实现高分辨率星载sar成像提供支撑.本发明提供了一种敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法，包含以下三个主要步骤：瞄准点和时序规划、需求姿态角粗算和需求姿态角精算。步骤1——瞄准点和时序规划1.1)聚束模式工作时雷达波束中心始终瞄准地面待观测场景的几何中心，选择该几何中心作为地面瞄准点.1.2)根据瞄准点坐标以及星历数据，寻找瞄准点回波多普勒中心为零的星历时刻，将其作为成像时序的中心。根据需求方位向分辨率和多普勒调频率计算完整的成像时间，并平均分布于时序中心两侧.1.3)将地面瞄准点和地心坐标变换到轨道坐标系，为后续姿态角计算提供输入.步骤2——需求姿态角粗算2.1)根据天线相位中心和地面瞄准点在轨道系下的坐标，计算期望的波束中心指向向量；2.2)根据卫星质心和地心在轨道系下的坐标，计算起始z轴指向向量；2.3)由于雷达波束可能存在离轴角，根据旋转变换关系，计算波束起始指向向量；2.4)根据欧拉定理，计算以欧拉轴/角参数式表示的需求姿态；2.5)根据欧拉四元数式与欧拉轴/角参数式的转换关系计算需求的欧拉四元数；2.6)根据欧拉角式与欧拉四元数式的转换关系可以计算得到初始欧拉角。步骤3——需求姿态角精算3.1)利用上一步得到的初始欧拉角，计算星本体坐标系到轨道坐标系的转换矩阵，并将天线相位中心星本体坐标系坐标变换成轨道坐标系坐标；3.2)以此坐标为输入，重复步骤二的计算过程，得到新的欧拉角；3.3)比较新欧拉角与初始欧拉角之间的差值，如大于预先设定的门限值，则用新欧拉角代替初始欧拉角，重复上述过程，直到新欧拉角与初始欧拉角之间的差值小于预先设定的门限值，完成姿态机动需求计算.因此，与现有技术相比，本发明可以实现以下的有益效果：1)充分考虑了卫星轨道和地球表面的弯曲特性影响，以及聚束模式星载sar工作特点，为高分辨率星载sar成像提供了一种经济、高效的实现方式；2)充分考虑了可能存在的波束指向偏离星本体坐标系z轴的离轴角，适用于同时存在电 扫描和整星姿态机动的情况，具备更高的普适性，且为机电联合扫描实现聚束模式星载sar成像提供技术基础；3)能适应sar天线在星体上各种不同的安装位置，将天线相位中心与卫星质心偏差的影响降低到可以忽略的程度.附图说明图1是本发明的敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法的流程图；图2是本发明的敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法的计算原理示意图；图3是本发明的敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法的实施验证流程图；图4示出了本发明实施例所涉及的以卫星轨道系为参考系312转序下的需求欧拉角；图5示出了本发明实施例的stk软件中虚拟sar卫星的起始时刻波束指向；图6示出了本发明实施例的stk软件中虚拟sar卫星的中间时刻波束指向；图7示出了本发明实施例的stk软件中虚拟sar卫星的结束时刻波束指向；图8示出了本发明实施例的地面瞄准误差.具体实施方式下面结合附图1-8及具体实施方式对本发明进行详细说明.本发明通过迭代计算得到敏捷sar卫星聚束模式姿态机动需求指标参数，具体的方法流程图如图1所示，包含以下步骤：步骤一，瞄准点和时序规划1.1确定瞄准点为计算姿态机动需求指标参数，需要明确不同时刻雷达波束需要指向的位置信息。聚束模式工作时雷达波束中心始终瞄准地面待观测场景的几何中心，因此，选择该几何中心作为整个成像时间内的地面瞄准点。该瞄准点的具体坐标由任务规划和系统波位设计结果决定，观测任务确定后可以计算得到。1.2时序规划接下来需要确定成像时间序列，为保证回波信号多普勒中心为零，降低数据处理难度，需要以回波多普勒中心为零的星历时刻作为成像时序中心，再根据分辨率需求计算成像总时间，并将其对称分布于成像时序中心两侧，完成时序规划。为寻找瞄准点回波多普勒中心为零的星历时刻，首先需要计算各个星历时刻的回波多普勒中心，计算方法如下：其中，为地心惯性坐标系下卫星的位置矢量，为地心惯性坐标系下瞄准点的位置矢量，为地心惯性坐标系下卫星的速度矢量，为地心惯性坐标系下瞄准点的速度矢量，λ为载波波长，rst为卫星与瞄准点之间的距离。根据上式可以得到所有星历时刻的瞄准点回波多普勒中心值，进而得到中心值为零的星历时刻.接下来根据需求方位向分辨率ρa和多普勒调频率fa计算成像时间，计算公式如下：其中，vg为零多普勒线扫过地面的速度，kwa为多普勒信号处理加权扩展因子，fa为多普勒调频率，由下式计算得到其中，为地心惯性坐标系下卫星的加速度矢量，为地心惯性坐标系下瞄准点的加速度矢量。将成像时间长度ts对称分布于成像时序中心t0两侧，成像时间范围为完成了时序规划。本发明中涉及的坐标系变换方法可查阅卫星控制经典教材，这里不再累述.1.3计算瞄准点和地心轨道坐标系坐标将地面瞄准点和地心坐标变换到轨道坐标系，为后续姿态角计算提供输入。应了解，这个步骤是为了将相关点的坐标变换到同一坐标系下，以便于后续计算.步骤二，需求姿态角粗算计算需求姿态角时，应根据时序规划得到的成像星历时刻逐点进行计算。计算的基本思路是将起始指向向量和期望指向向量变换到同一坐标系下，再利用最短旋转路径准则，计算实现天线波束期望指向的旋转轴和转动角，进而得到所有需求姿态信息。卫星整星构型确定后，可以准确测得天线相位中心在卫星本体坐标系下的坐标，但因姿态角未知，不能得到天线相位中心在轨道坐标系下的精确坐标。粗算阶段，以其在卫星本体 坐标系下的坐标代替轨道坐标系坐标，由于上述两个坐标系原点同为卫星质心，且星体尺寸一般较小，两坐标系下的相位中心坐标值差一般为米级，可以用于姿态角粗算，这种近似引入的误差将在精算阶段予以消除。2.1计算期望指向向量依据上述姿态角计算思路，首先需要计算期望指向向量。根据天线相位中心op和地面瞄准点t在轨道系下的坐标，可得到期望的波束中心指向向量，以表示。根据卫星质心os和天线相位中心op在轨道系下的坐标，可得到卫星质心到相位中心的向量，以表示.由于姿态机动时围绕卫星质心os转动，应以卫星质心为起点计算期望指向向量，以表示，根据向量关系得到：如图2所示，示出了上述的计算原理。2.2计算z轴指向向量接下来需要计算波束起始指向向量，首先计算星本体坐标系z轴指向向量，随后根据该向量和离轴角，通过旋转变换得到波束起始指向向量。根据卫星质心os和地心oe在轨道系下的坐标，可得到起始z轴指向向量，以表示。2.3计算波束起始指向向量由于雷达波束可能存在离轴角，根据旋转变换关系，可得到如下波束起始指向向量：其中，c(θ)为由波束离轴角决定的旋转转换矩阵，θ为离轴角，是天线雷达波束中心偏离卫星本体坐标系z轴的角度，包含天线安装时的斜装角，天线自身的电扫描角等.2.4计算以欧拉轴/角参数式表示的需求姿态首先，计算以欧拉轴/角参数式表示的需求姿态信息，根据欧拉定理，卫星绕质心的姿态机动可以由绕通过质心的某一固定轴转动一定角度得到。该转轴方向矢量用表示，转动角度用φ表示，这即是欧拉轴/角参数式.根据最短路径旋转准则，和φ可由知通过下两式计算得到：其中，×为向量积，·为标量积，|·|为向量长度。2.5计算以欧拉四元数式表示的需求姿态根据欧拉四元数式与欧拉轴/角参数式的转换关系可以计算得到以欧拉四元数表示的需求姿态信息.2.6计算以欧拉角式表示的需求姿态根据欧拉角式与欧拉四元数式的转换关系可以计算得到初始欧拉角.由于在本阶段的计算过程中，使用天线相位中心卫星本体坐标系坐标代替轨道坐标系坐标，所以初始欧拉角中包含一定的误差，下一步将重点剔除该误差。本步骤所涉及的欧拉轴/角参数式、欧拉四元数式、欧拉角式等姿态表示方法之间的转换关系可查阅卫星控制经典教材，在本说明书中不再累述。步骤三，需求姿态角精算将步骤二中得到的欧拉角作为初始姿态角，利用初始姿态角计算天线相位中心在轨道系下的坐标，重复步骤二，得到新的欧拉角，再将新的欧拉角作为初始姿态角，循环迭代计算，直到新旧欧拉角之间误差满足需求精度。由于计算新欧拉角所用的转换矩阵总是比计算初始欧拉角时的转换矩阵更精确，因此，新欧拉角逐渐向准确的需求姿态角逼近，而新欧拉角与初始欧拉角之差逐渐向零收敛。通过选择合适的门限值，能得到满足精度需求的姿态角。这个迭代计算有效剔除了天线相位中心与质心不重合导致的指向误差，得到满足精度要求的姿态角，本步骤具体包含如下过程.3.1计算轨道系天线相位中心坐标利用上一步得到的初始欧拉角，计算星本体坐标系到轨道坐标系的转换矩阵，并将天线相位中心星本体坐标系坐标变换成轨道坐标系坐标。3.2重复步骤二计算欧拉角式表示的需求姿态然后以此坐标为输入，重复步骤二的计算过程，得到新的欧拉角.3.3判决计算结果是否达到判决门限比较新欧拉角与初始欧拉角之间的差值，如大于预先设定的门限值，则用新欧拉角代替初始欧拉角，重复上述过程，直到新欧拉角与初始欧拉角之间的差值小于预先设定的门限值。可根据指向精度的要求确定门限值，本发明方法收敛速度较快，3～4次迭代后新欧拉角与初 始欧拉角之间的差值将达到10-20弧度量级，故该门限值可以设置为10-20弧度。事实上，由于控制系统的执行机构精度有限，更高的计算精度并没有实际价值.通过matlab和卫星仿真工具包(satellitetoolkit，以下简称为stk)软件联合仿真，可以验证本发明提出的方案的有效性，验证流程如图3所示。首先，根据输入条件，在matlab软件中计算姿态机动指标参数，然后将姿态机动指标参数计算结果导入stk软件，显示姿态机动和波束指向变化，并将地面瞄准点位置信息输回matlab软件，分析stk软件中得到的瞄准点位置信息与期望的瞄准点位置信息之间的误差.由于stk软件的正确性已经经过工程验证，它的分析结果能准确地反映卫星运动规律，故上述仿真实验具有较高的可信度。本次仿真实验的sar卫星聚束模式输入参数如以下表1所示，根据本发明的方案，得到实现聚束模式的姿态机动需求指标参数.图4为以轨道系为参考系，312转序下的需求欧拉角.可见，聚束模式成像过程中俯仰机动范围最大，达到1.8835°，偏航机动次之，达到0.4768°，横滚机动最小，为0.0048°。表1参数数值轨道长半轴(km)7078.14轨道偏心率0.01轨道倾角(°)98.1928近地点幅角(°)0升交点赤经(°)0天线离轴角(°)10星本体系下天线相位中心位置(m)[5，5，5]雷达载波波长(m)0.0311多普勒信号处理加权展宽因子1.1期望方位向分辨率(m)0.6期望场景中心经度坐标(°)169.9期望场景中心纬度坐标(°)12.1将上述姿态角导入stk软件，通过stk软件中虚拟sar卫星的雷达波束指向，验证本发明提出方法的正确性和有效性。stk软件中虚拟卫星同样采用表1中的输入参数，在注入计算得到的姿态角后，得到成像起始、中间和结束时刻波束指向，以及地面瞄准点位置(如图5、6、7所示)。在附图中，叉形区域为期望的成像场景中心(瞄准点)，矩形区域为注入本发明方法计算得到的需求姿态角后，虚拟卫星的雷达波束地面足印，可见雷达波束在成像起始、中间和结束时刻都准确地指向了期望的场景中心，证明了上述姿态机动指标参数计算结果的正确性。再将stk软件中得到的实际瞄准点位置信息导回到matlab软件中进行分析，得到了详细的瞄准误差.如图8所示，瞄准误差在10-8米量级，可以忽略不计。综上所述，本发明的敏捷sar卫星聚束模式姿态机动需求计算方法，考虑了雷达波束离轴角的影响，能够适用于包含任意离轴角的敏捷sar卫星聚束模式姿态机动需求计算，也可应用于机械扫描和电扫描联合实现的sar卫星聚束模式姿态机动需求计算。另外，通过迭代计算剔除雷达天线相位中心偏离星本体坐标系原点的影响，得到满足精度要求的需求姿态，该算法能够适用于天线安装于星体任何位置的敏捷sar卫星聚束模式姿态机动需求计算。本发明中未说明部分属于本领域的公知技术。当前第1页12