基于蒙特卡罗法的星载太阳辐照度光谱仪对日指向误差分析

基于蒙特卡罗法的星载太阳辐照度光谱仪对日指向误差分析 Analysis of sun pointing error of spaceborne solar spectroradiometer based on Monte Carlo method 1 2 corresp first-author E-mail： guoxuhuhu@163.com E-mail： guoxuhuhu@163.com 郭　旭（1991-），女，吉林长春人，硕士，助理研究员，2015年于北京航空航天大学获得硕士学位，主要从事航天遥感电子学设计、伺服系统设计等研究工作。E-mail：guoxuhuhu@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226751&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226756&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226753&type= 22.01333618 32.00407028 郭　旭 （1991-），女，吉林长春人，硕士，助理研究员，2015年于北京航空航天大学获得硕士学位，主要从事航天遥感电子学设计、伺服系统设计等研究工作。E-mail： guoxuhuhu@163.com guoxuhuhu@163.com 1 1 颜昌翔（1973-），男，湖北洪湖人，博士，研究员，2001年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事空间光学遥感仪器的光机电一体化技术，多光谱、超光谱空间遥感成像技术、偏振探测技术等方面的研究。E-mail：yancx0128@126.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226758&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226762&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226760&type= 22.01333618 32.00402832 颜昌翔 （1973-），男，湖北洪湖人，博士，研究员，2001年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事空间光学遥感仪器的光机电一体化技术，多光谱、超光谱空间遥感成像技术、偏振探测技术等方面的研究。E-mail： yancx0128@126.com yancx0128@126.com 1 1 1 中国科学院 长春光学精密机械与物理研究所，长春 吉林 130033 中国科学院 长春光学精密机械与物理研究所，长春 吉林 130033 Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 中国科学院大学，北京 100049 中国科学院大学，北京 100049 University of Chinese Academy of Sciences， Beijing 100049， China University of Chinese Academy of Sciences， Beijing 100049， China 为保证光谱仪能够获取到有效太阳光谱数据，要求搭载光谱仪的二维转台能够长时间高精度跟踪太阳，而跟踪太阳的前提是转台实现精准指向，使太阳进入导行镜的有效视场。根据从卫星接收的太阳矢量和姿态数据对安装误差进行修正，计算修正后转台指向的参考值，并确保指向误差小于1°。在转台各个基准立方镜上建立坐标系，采用坐标变换法建立了对日指向的数学模型，通过运动学反解给出转台调整角度与轨道坐标系下太阳矢量的解析关系。然后，结合光谱仪的在轨运动形式，对其工作过程中的主要误差来源进行分析，并在MatLab平台上搭建基于蒙特卡罗法的指向误差模型。仿真结果表明指向误差优于0.35°。最后，在地面开展太阳指向模拟实验，测得春分轨道下指向误差低于0.16°，满足导行镜捕获视场需求。实验结果验证了数学模型的正确性及有效性，同时为光谱仪在轨工作提供了设计参考。 To guarantee the efficient operation， that is， capture of solar spectral data， of a solar spectroradiometer， the two-dimensional pointing turntable that carries the spectrometer should precisely track the sun over long periods. Highly precise sun pointing ensures that the sun stays in the effective field of vision of the guide telescope. Specifically， according to the solar vector and positioning data received from the satellite， the assembly error is corrected and the reference data for pointing are calculated after correction， ensuring that the pointing error is less than 1°. To accomplish this， coordinate systems were built into all the hexahedral reference prisms on the turntable. and the mathematical model of sun pointing was implemented using the transformation matrices of the coordinate systems. Further， the relationship between the turntable's angle of adjustment and the solar vector on the orbit coordinate system was derived using inverse kinematics. then the main source of error during operation of the spectrometer was analyzed along with its mode of operation in orbit； from these analyses， the pointing error model was built based on the Monte Carlo method on MATLAB. The model was employed to perform a simulation and analysis of sun pointing. The simulation results demonstrate that the pointing error is less than 0.35°. Finally， the simulationwas performed to obtain the result that the pointing error under the vernal equinox orbit is less than 0.16°， meeting the requirement for the field of view of the guide telescope. The correctness and effectiveness of the mathematical model presented in this paper are validated based on the test results， which simultaneously offer a design reference for the spectrometer. 星载光谱仪 坐标变换 运动学反解 蒙特卡洛 指向误差 spaceborne spectrometer transformation of coordinate system inverse kinematics Monte Carlo pointing error 1 引　言 1 太阳辐照度光谱仪作为太阳光谱监测的主要光学载荷，为气象领域提供了大量的数据支撑。光谱仪载荷主要包括不同谱段的谱仪设备、用于对日指向跟踪的导行镜、以及用于提供指向功能的二维转台等部件。目前，国内外对于太阳辐照度光谱仪的研究均取得了一定的进展，研究内容主要围绕光谱仪载荷本身的设计以及光谱仪对日指向精度等方面。 在光谱仪载荷设计方面，李占峰 ［ 1 1 ］ 等研究了多通道光谱仪中光谱仪输出波长与光谱仪内部调节丝杠的位移非线性问题，并给出了相应的理论公式。曹佃生 ［ 2 2 ］ 等根据太阳光谱仪中对于波长重复性的指标要求，对光谱仪内部的波长扫描机构进行了设计，并对扫描机构的精度进行了分析。李寒霜 ［ 3 3 ］ 研究了太阳辐照度光谱仪中紫外-真空紫外波段的光谱定标以及辐射定标的问题，并构建了一套可以实现150~300 nm波段测量的定标系统。孙德贝 ［ 4 4 ］ 等根据光的叠加原理研制了一台太阳光谱仪探测系统线性度的测试装置，该装置的工作波段为200~2 400 nm，可模拟多谱段的太阳光谱特性。李福田 ［ 5 5 ］ 等针对高精度空间太阳光谱仪的辐射定标问题，研制了一种采用数字微镜器件的光谱辐亮度标准光源，为太阳光谱仪的辐射定标提供硬件基准。李占峰 ［ 6 6 ］ 等介绍了一种由卫星平台当前广播事件和轨道瞬根来推导预报短时太阳角度的方法，该方法预报的最大角度误差为0.5°，预热时间的最大偏差为20″。孙立微 ［ 7 7 ］ 针对太阳高光谱的辐射定标问题，研制了由激光二极管配合钨灯的定标光源，为太阳辐照度光谱仪的辐射定标提供了良好的光源基准。窦晨浩 ［ 8 8 ］ 针对太阳光谱仪在轨工作过程中获取到的偏振测量数据、修正数据以及大气探测数据等，采用定量化分析法进行了分析与研究。 在光谱仪对日指向精度方面，其研究主要为二维指向机构、导行镜及光谱仪在各环节中的加工、装调误差最终合成后对于太阳指向性能的影响。王金元 ［ 9 9 ］ 等分析了一种小口径的相控阵天线的指向误差，并提出了降低各项误差的方法，可有效提高指向精度。胡晓炜 ［ 10 10 ］ 等以地平式望远镜为模型，对其指向误差进行建模及分析，并介绍了指向误差数据的获取方法。通过修正系数可有效降低设备的指向误差，提升指向精度。康跃然 ［ 11 11 ］ 等分析了三轴转台姿态变化时各误差分量对于最终的系统指向误差的影响，同时通过仿真获取各个误差分量的变化规律，为后续的误差分配及补偿提供了理论依据。吴伟平 ［ 12 12 ］ 等建立了在轨视轴临边指向的误差补偿方法，消除了探测仪器由于外界客观因素所导致的指向误差。该方法可以将指向误差控制在±3.08″的范围内。 为实现光谱仪对日指向，将太阳引入导行镜视场，需要建立太阳矢量、卫星姿态和安装误差等因素与转台转角的数学关系。本文首先采用坐标变换法建立光谱仪载荷系统的坐标系及相应的变换矩阵的数学模型；然后通过蒙特卡罗法分析载荷整机的各项误差，仿真得出对日指向误差；最后开展对日指向模拟实验，验证了理论分析及仿真结果的准确性及有效性。 2 光谱仪系统模型 1 2.1　系统组成 2 根据太阳辐照度光谱仪工作的轨道参数，一年内太阳光在轨道面的入射角在约47°的范围内变化，远大于光谱仪的有效视场。根据光谱仪的工作特点和任务需求，需设计二维转动机构，对卫星偏航和俯仰方向的轨道运动进行补偿。 星载太阳辐照度光谱仪载荷整机是由两台太阳某谱段的光谱仪、一台导行镜、以及搭载上述光学有效载荷的二维转台组成。二维转台作为该载荷的核心，主要由基座、方位轴系（U形架）和俯仰轴系（O形架）等部组件构成。方位U形架为外框架，安装在转台基座上，俯仰O形架为内框架，安装在U形架上，导行镜和光谱仪安装在O形架上。光谱仪二维转台的三维模型如 图1 图1 所示。 2.2　工作原理 2 载荷在轨工作时，首先调整二维转台方位、俯仰角位置实现对日指向功能，使太阳进入导行镜有效视场；然后依据导行镜偏移量实时调整两轴角位置实现对日跟踪；只有当跟踪误差低于一定范围时，即太阳小角度偏离谱仪视场中心时，两台光谱仪方能获取有效的光谱数据。因此，转台对日指向是载荷能够有效工作的前提。 由于卫星姿态的实时调整、载荷整机相对卫星平台存在安装误差以及载荷内部组件间存在安装误差等因素，需要实时计算二维转台的工作角度。具体地，根据卫星平台的广播数据，获取当前轨道坐标系下的太阳矢量和卫星姿态角测量值，并对安装误差进行修正，进而计算出转台方位轴和俯仰轴的转动角度，最终保证太阳位于导行镜的有效视场内，并尽量靠近视场中心。 3 对日指向数学模型 1 3.1　系统坐标系建立 2 载荷整机的坐标系建立方式由它在卫星平台上的安装方式决定。对于卫星平台， X 方向为飞行方向， Z 方向指向地心， Y 方向由右手法则确定。载荷在卫星平台上的安装方式如 图2 图2 所示。 首先建立轨道坐标系 O （ O O -X O Y O Z O ），坐标原点 O O 位于卫星质心， Z O 指向地心， X O 轴位于轨道平面内并指向卫星飞行方向， Y O 通过右手坐标系确定。卫星坐标系 S （ O S -X S Y S Z S ），坐标原点 O S 与 O O 重合，卫星无姿态运动时， O 系与 S 系重合，卫星的三轴姿态角 φ ， θ ， ψ 即指 S 系在 O 系内的三轴姿态描述。确定卫星坐标系后，如 图1 图1 所示，在转台基座上的精测棱镜1上建立基座坐标系 C （ O C -X C Y C Z C ），若忽略精测棱镜的安装误差、角度加工误差， X C ， Y C ， Z C 分别与精测棱镜1的3个表面垂直，正方向与卫星坐标系相同，后续棱镜坐标系的定义与之相同，不再赘述。在方位轴U形架的精测棱镜3上建立方位轴坐标系 A （ O A -X A Y A Z A ），在俯仰轴O形架的精测棱镜2上建立俯仰轴坐标系 E （ O E -X E Y E Z E ），在忽略U形架在载荷基座的安装误差时，方位零位时 A 系与 C 系三轴平行。类似地，若忽略O形架在U形架的安装误差时，俯仰零位时 E 系与 A 系三轴平行。在导行镜精测棱镜4上建立导行镜坐标系 G （ O G -X G Y G Z G ）， -Y 方向代表视轴。 3.2　坐标变换矩阵推导 2 已知各个组件的局部坐标系后，基于坐标变换的原理，太阳矢量从轨道坐标系下的矢量表达至导行镜坐标系下的矢量表达式为： <math display="block" id="M1"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>G</mi></mrow></msub><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>E</mi><mi>G</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>E</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>A</mi><mi>E</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>A</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>C</mi><mi>A</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>S</mi><mi>C</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>O</mi><mi>S</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226777&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226775&type= 87.29133606 3.81000018 ， 其中： S G 为导行镜坐标系下的太阳矢量，规定其坐标为［0　-1　0］ T ； S O 为轨道坐标系下的太阳矢量，假设其坐标为［ X sun 　 Y sun 　 Z sun ］ T ，该太阳矢量数据由平台广播数据提供； M OS 是由卫星相对于轨道姿态变化引起的变换矩阵； M SC 为卫星坐标系至转台基座坐标系的变换矩阵，通过标定精测棱镜1的姿态可以获得； M CA 为基座坐标系至方位轴系坐标系的变换矩阵，通过标定精测棱镜1和3的姿态可以获得； M AE 为方位轴系坐标系至俯仰坐标系的变换矩阵，通过标定精测棱镜2和3的姿态可以获得； M EG 为俯仰轴坐标系至导行镜坐标系的变换矩阵，通过标定精测棱镜2和4的姿态可以获得； M A ， M E 分别为方位轴系、俯仰轴系转动前后局部坐标系的变换矩阵。 在明确上述变换环节中各个矩阵的意义后，对各变换矩阵的形式进行推导 ［ 13 13 - 15 15 13-15 ］ 。假设卫星相对轨道坐标系的滚动角（绕 X O 轴）、俯仰角（绕 Y O 轴）和偏航角（绕 Z O 轴）分别为 φ ， θ ， ψ ，则 M OS 为： <math display="block" id="M2"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>O</mi><mi>S</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ψ</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ψ</mi><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ϕ</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ϕ</mi></mtd></mtr><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ψ</mi><mo>+</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ψ</mi><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ϕ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ψ</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>ϕ</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226782&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226780&type= 126.74599457 16.08666611 . M SC 可根据卫星平台精测棱镜及转台精测棱镜1的经纬仪标定结果给出，其值如下： <math display="block" id="M3"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>S</mi><mi>C</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="right"><mtr><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">28</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">20</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">30</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">08</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">19</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">08</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226783&type= 67.98733521 13.54666710 . 同理， M CA ， M AE 及 M EG 也可根据经纬仪对相关精测棱镜的标定结果给出，模型中4个精测棱镜的标定结果详见 表1 表1 。其中，棱镜 i - i 的值代表自准值；棱镜 i - j 的值代表互瞄值， i ， j 的取值均为1~4。 棱镜名称 棱镜1 棱镜2 棱镜3 棱镜4 2 棱镜1 0 — — 174°43′28″ 89°59′54″ — — 90°11′40″ 2 棱镜2 — 0 — 175°27′19″ — 90°1′48″ — 91°49′20″ 2 棱镜3 — — 0 172°22′51″ — — 89°59′50″ 90°02′12″ 2 棱镜4 5°18′25″ 4°35′24″ 352°23′41″ 0 89°48′12″ 88°10′35″ 89°57′42″ 90°0′0″ 根据 表1 表1 的数据求得 M CA ， M AE 及 M EG 的表达式，分别为： <math display="block" id="M4"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>C</mi><mi>A</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="right"><mtr><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">31</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">73</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">31</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">08</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">73</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">08</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226793&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226791&type= 68.32600403 13.54666710 ， <math display="block" id="M5"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>A</mi><mi>E</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="right"><mtr><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">55</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">51</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">55</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">48</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">51</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">47</mn></mtd><mtd><mtext> </mtext><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226797&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226795&type= 68.32600403 13.54666710 ， <math display="block" id="M6"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>E</mi><mi>G</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="right"><mtr><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">79</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">02</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">79</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd><mtd><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">52</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">02</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.000</mn><mtext> </mtext><mn mathvariant="normal">52</mn></mtd><mtd><mn mathvariant="normal">1.000</mn><mtext> </mtext><mn mathvariant="normal">00</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226801&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226799&type= 68.32600403 13.54666710 . M A 为方位轴转动前后的变换矩阵，其值为： <math display="block" id="M7"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>A</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mtext>    </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>A</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>A</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226805&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226803&type= 42.58733368 13.54666710 . 同理， M E 的形式为： <math display="block" id="M8"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>E</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">1</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mtext>    </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>E</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>E</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>E</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>E</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226809&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226807&type= 41.90999985 13.54666710 . 3.3　运动学反解结果 2 当忽略卫星姿态角度、载荷整机在卫星平台上的安装误差以及方位轴系、俯仰轴系、导行镜等各组件在转台上的安装误差时，各个坐标系的坐标轴平行，变换矩阵简化为单位阵。则 式（1） 式（1） 化简为： <math display="block" id="M9"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>G</mi></mrow></msub><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">Μ</mi></mrow><mrow><mi>E</mi></mrow></msub><mo>⋅</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>A</mi></mrow></msub><mo>⋅</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226813&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226811&type= 28.19400024 3.81000018 . 将 式（7） 式（7） 和 式（8） 式（8） 代入 式（9） 式（9） 后，可解算出转台方位轴及俯仰轴的转动角度为： <math display="block" id="M10"><mi>A</mi><mo>=</mo><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>X</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226817&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226815&type= 27.51666832 8.46666718 ， <math display="block" id="M11"><mi>E</mi><mo>=</mo><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>Z</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226821&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226819&type= 25.73866653 3.72533321 . 上述求解出的转台运动角度为忽略平台姿态调整和安装误差等因素得到的结果，即理论真值。但实际工程中会存在多个环节的误差，为兼顾工程中实际误差的必然性以及解析计算的便利性，鉴于 式（5） 式（5） 俯仰轴系与方位轴系的安装误差相对较小，考虑软件的在轨计算效率，将它忽略不计。其余误差变换矩阵纳入计算环节，通过坐标变换矩阵反解出转台方位轴及俯仰轴的工作转角。 首先，计算太阳矢量在转台基座坐标系下的矢量表达，记为： <math display="block" id="M12"><mi mathvariant="bold-italic">B</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>C</mi><mi>A</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>S</mi><mi>C</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>B</mi><mi>S</mi></mrow></msub><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226825&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226823&type= 55.54133606 15.40933418 . 其次，得到太阳矢量在俯仰轴坐标系下的矢量表达，记为： <math display="block" id="M13"><mi mathvariant="bold-italic">A</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><msubsup><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi>E</mi><mi>G</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>G</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226829&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226827&type= 33.95133591 15.40933418 . 将上述表达式代入 式（1） 式（1） ，即： <math display="block" id="M14"><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">1</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mtext>   </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>E</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>E</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>E</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>E</mi></mtd></mtr></mtable></mrow></mfenced><mo>×</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mtext>    </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>A</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>A</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mrow></mfenced><mo>×</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226833&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226831&type= 90.50866699 15.40933418 . 通过 式（14） 式（14） 可以反解出方位轴及俯仰轴的转动角度，即： <math display="block" id="M15"><mtable columnalign="left"><mtr><mtd><mi>A</mi><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mroot><mrow><msubsup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow/></mroot></mrow></mfrac><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mi>E</mi><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mroot><mrow><msubsup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow/></mroot></mrow></mfrac><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226837&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226835&type= 51.13866806 20.82799911 . 4 指向误差分析 1 由第三部分的计算结果，可以得到任一组太阳矢量及相应安装误差下的转台工作转角。由于各项误差均存在一定的分布区间，需要保证在统计学角度的极端情况下，太阳矢量依然位于导行镜视场中心。 4.1　误差仿真 2 误差仿真流程如下： （1）指定某种工况并设定仿真模型的参数，将上述参数作为仿真场景的真值，参数包括太阳矢量、卫星姿态和安装误差等； （2）根据第一步设定的仿真理论参数，计算转台方位、俯仰调整角的真值； （3）根据平台姿态和太阳矢量的数据误差，以及转台内部各部组件的安装测量误差，将第一步的无误差数据误差化，得到实际的转台方位、俯仰目标调整角； （4）根据第三步的计算结果控制转台转动，转台定位误差为编码器的测量误差，可得实际控制输出的方位、俯仰调整角，将实际值与真值比较，得出转台对日指向的最终实际误差。 4.2　蒙特卡罗分析 2 根据前文内容可知，系统工作过程中有多个环节存在误差，从误差性质上可分为系统误差和随机误差，从误差种类上分为运动误差、安装误差、加工误差及测量误差等。分析多个误差源对最终结果的影响，需要引入基于蒙特卡罗法的误差模型。 蒙特卡罗法是基于变量的统计学特性，使用随机变量代替常量，同时保证随机变量满足一定的概率分布，使最终的计算结果更加逼近使用常量所得到的结果 ［ 16 16 ］ 。鉴于本系统中误差项较多，在后续的计算中将每一项误差均视为一个随机变量，并且它满足正态分布。从分布区间中取出若干组数据代入公式进行计算，从而得到最终结果的分布情况，即可确认载荷的整体误差是否满足工况需求。 根据二维转台内各个精测棱镜的测量结果，以及外部提供的设计输入，可以获得下列仿真参数： （1）太阳矢量 <math id="M16"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi>O</mi></mrow></msub><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>X</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></mtd><mtd><msub><mrow><mi>Y</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></mtd><mtd><msub><mrow><mi>Z</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0.085</mn></mtd><mtd><mo>-</mo><mn mathvariant="normal">0.981</mn></mtd><mtd><mn mathvariant="normal">0.176</mn></mtd></mtr></mtable></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226842&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226839&type= 83.48133087 4.74133348 ； （2）卫星姿态［ <math id="M17"><mtable><mtr><mtd><mi>φ</mi></mtd><mtd><mi>θ</mi></mtd><mtd><mi>ψ</mi></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226847&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226845&type= 13.54666710 3.64066648 ］ <math id="M18"><mo>=</mo><mo stretchy="false">[</mo><mn mathvariant="normal">0.060</mn><mtext> </mtext><mn mathvariant="normal">3</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226851&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226849&type= 16.84866714 3.55599999 <math id="M19"><mo>-</mo><mn mathvariant="normal">0.051</mn><mtext> </mtext><mn mathvariant="normal">8</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226855&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226853&type= 15.66333389 2.79399991 <math id="M20"><mn mathvariant="normal">0.028</mn><mtext> </mtext><mn mathvariant="normal">3</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226858&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226856&type= 11.93799973 2.79399991 ］； （3）载荷基座相对于平台的安装误差［ <math id="M21"><mtable><mtr><mtd><msub><mrow><mi>φ</mi></mrow><mrow><mi>C</mi></mrow></msub></mtd><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mi>C</mi></mrow></msub></mtd><mtd><msub><mrow><mi>ψ</mi></mrow><mrow><mi>C</mi></mrow></msub></mtd></mtr></mtable><mo stretchy="false">]</mo><mo>=</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226862&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226860&type= 24.21466637 3.72533321 ［ <math id="M22"><mtable><mtr><mtd><mn mathvariant="normal">0.013</mn><mtext> </mtext><msup><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mn mathvariant="normal">0.012</mn><mtext> </mtext><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mn mathvariant="normal">0.016</mn><mtext> </mtext><msup><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226866&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226864&type= 43.09533310 2.79399991 ］； （4）转台轴系相对于基座的安装误差 <math id="M23"><mo stretchy="false">[</mo><mtable><mtr><mtd><msub><mrow><mi>φ</mi></mrow><mrow><mi>A</mi><mi>E</mi></mrow></msub></mtd><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mi>A</mi><mi>E</mi></mrow></msub></mtd><mtd><msub><mrow><mi>ψ</mi></mrow><mrow><mi>A</mi><mi>E</mi></mrow></msub></mtd></mtr></mtable><mo stretchy="false">]</mo><mo>=</mo><mo stretchy="false">[</mo><mtable><mtr><mtd><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">031</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mo>-</mo><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">073</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">014</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd></mtr></mtable><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226870&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226869&type= 75.09933472 3.72533321 ； （5）导行镜相对于转台轴系的安装误差 <math id="M24"><mo stretchy="false">[</mo><mtable><mtr><mtd><msub><mrow><mi>φ</mi></mrow><mrow><mi>G</mi></mrow></msub></mtd><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mi>G</mi></mrow></msub></mtd><mtd><msub><mrow><mi>ψ</mi></mrow><mrow><mi>G</mi></mrow></msub></mtd></mtr></mtable><mo stretchy="false">]</mo><mo>=</mo><mo stretchy="false">[</mo><mtable><mtr><mtd><mo>-</mo><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">03</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">083</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd><mtd><mo>-</mo><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">045</mn></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msup></mtd></mtr></mtable><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226875&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226873&type= 73.57533264 3.72533321 。 在MatLab中建立变换矩阵的数学模型，并代入上述仿真参数，可以求得当方位轴转角 A 0 =3.951°，俯仰轴转角 E 0 = - 13.316 9°时，太阳处于导行镜视场中心。 上述计算是基于各项误差为常数，下面引入基于蒙特卡罗的误差模型。采用正态分布函数，并适当选取期望及标准差。其中，太阳矢量误差和姿态测量误差通过分析卫星运动的仿真数据和在轨实测数据预估；各项安装误差修正后为经纬仪测量误差，依据工程经验预估；时延误差针对在轨软件的具体实现进行估算；转台定位误差取决于编码器的测量误差。各项误差的估计值详见 表2 表2 。 误差源 期望 μ 标准差(3 σ )/(°) 预估依据 太阳矢量误差 0 0.238 0 仿真最大值 姿态测量误差 0 0.006 0 仿真最大值 仪器安装误差 0 0.002 8 3 测量误差， 与测试工具、测试方法有关， 参考经验 转台安装误差 0 0.002 8 跟踪探测器安装误差 0 0.002 8 时延误差 0 0.017 5 计算最大值 转台定位误差 0 0.014 0 编码器误差 4.3　结果与讨论 2 根据 表2 表2 中的各项误差值，在各个误差维度下分别生成10 000组随机数进行仿真。经过仿真，转台方位转角及俯仰转角数据的直方图如 图 3 图 3 所示。从 图3 图3 可以看出，转台转角的仿真值大部分分布在理论真值附近，在方位、俯仰两个维度上呈正态分布规律。转台对日指向误差的仿真结果如 图4 图4 所示。 根据指向误差的仿真结果，方位轴的最大误差为0.338 7°，俯仰轴的最大误差为0.294 5°。从 图4 图4 可以看出，指向误差分布在1°的圆形视场中，且全部分布在0.35°以内，意味着通过调整转台的方位角和俯仰角，能够保证太阳进入导行镜的有效视场，仿真结果满足指向任务需求。 5 模拟实验 1 5.1　实验环境 2 为了验证上述理论分析及仿真结果的准确性及有效性，在实验室环境下，搭建光谱仪对日指向模拟实验，用于模拟转台对日指向过程。 实验环境如下：在导行镜前端放置一台平行光管，并在焦面处放置一个强光光源，以模拟太阳光，并且该光源在导行镜视场内的光斑大小与太阳近似；将二维转台置于高精度六自由度平台上，该平台的转动范围大于转台转动角范围，控制精度优于卫星平台的在轨姿态控制精度，能够模拟卫星姿态以及卫星飞行时太阳位置的变化。 初始状态为光源光线沿+ Y 轴方向射入导行镜视场。实验现场示意图如 图5 图5 所示。 5.2　指向实验 2 根据光谱仪载荷在轨的实际工作状态，选取有代表性的工况进行模拟实验。实验分两步：第一步是选取4种典型工况，令转台方位轴和俯仰轴自转台零位向不同方向转动，验证对日指向的数学模型；第二步是选取典型工作轨道，令转台在一轨工作中长时间保持对日指向，测量指向误差低于1°。具体实验流程如下： （1）测量光谱仪转台基座在六自由度平台上的安装误差，记为［ φ c 　 θ c 　 ψ c ］； （2）选取4组太阳矢量［ X sun 　 Y sun 　 Z sun ］、卫星姿态［ φ 　 θ 　 ψ ］，将参试设备安装误差的测量结果、选取的太阳矢量和平台姿态数据带入对日指向数学模型中，反解转台方位和俯仰目标角； （3）根据选取的数据分别驱动六自由度平台和二维转台转动； （4）记录导行镜测得的角偏移量，该参数为指向误差与导行镜测量误差的合成结果，其中导行镜的测量误差为3″（3 σ ）； （5）选取春分轨道的太阳矢量和卫星姿态数据，重复步骤（1）~（4），对转台长时间工作的指向误差进行验证和分析。 5.3　结果与讨论 2 将实验环境搭建完成后，首先测量光谱仪载荷基座相对于平台的安装误差为： <math display="block" id="M25"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></mtd><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></mtd><mtd><msub><mrow><mi>ψ</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0.011</mn><mo>°</mo></mtd><mtd><mn mathvariant="normal">0.005</mn><mo>°</mo></mtd><mtd><mn mathvariant="normal">0.007</mn><mo>°</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226917&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10226910&type= 63.83866501 4.74133348 . 然后依次测量转台方位轴和俯仰轴正、反转4种工况下的指向精度。实验参数（太阳矢量、卫星姿态参数和转台转角）和实验结果详见 表3 表3 。 2 序号 2 太阳矢量 [ X sun 　 Y sun 　 Z sun ] 2 卫星姿态/(rad) [ φ 　 θ 　 ψ ] 2 转台转角/(°) 2 角偏移量/(°) 方位轴 俯仰轴 方位轴 俯仰轴 1 [-0.152 1 -0.974 5 -0.164 7] [0.004 0.001 0.000 88] -8.891 6 9.240 2 0.012 3 -0.024 7 2 [-0.182 9 -0.974 5 0.129 9] [0.002 4 0.003 2 0.000 1] -10.630 4 -7.588 9 0.094 8 -0.050 3 3 [0.144 2 -0.974 4 0.172 0] [0.001 4 0.000 94 0.001 8] 8.345 2 -9.993 0 -0.112 0.100 8 4 [0.160 3 -0.974 4 -0.157 1] [0.000 5 0.000 64 0.000 73] 9.318 1 9.004 0 -0.074 9 -0.125 7 由实验结果中测得的导行镜偏移量数据可知，两轴最大指向误差为0.125 7°，低于1°，因此采用本文提出的对日指向数学模型，转台能够将太阳引入导行镜的有效视场。 重复上述实验过程，选取一轨春分轨道的太阳矢量和卫星姿态数据，令六自由度平台依据当前数据运动，实时解算转台转角驱动转台转动，测量600组导行镜偏移量，指向误差如 图6 图6 所示。由图可知，方位轴的最大指向误差为0.126 5°，俯仰轴的最大指向误差为0.154 2°，实验测得指向误差处于仿真分析的误差范围内。实验中六自由度平台未模拟太阳矢量的误差数据，因此实验测得的误差最大值稍低于仿真分析结果。 对日指向模拟实验结果一方面验证了本文提出的转台对日指向数学模型的准确性；另一方面，测得的指向误差验证了误差仿真分析结果，说明转台对日指向误差在导行镜的有效视场内，能够满足载荷在轨工作的要求。 6 结　论 1 光谱仪的对日高精度指向对卫星在轨工作至关重要。本文采用坐标变换矩阵法，对太阳辐照度光谱仪的系统坐标系及矩阵变换关系进行了分析并建立了相应的数学模型，得到了二维转台工作转角的解析解；利用基于蒙特卡罗法的误差模型分析光谱仪载荷系统的整体误差，并进行了对日指向误差仿真实验，实验结果表明指向误差优于0.35°；最后，通过转台对日指向模拟实验对所建立的数学模型和误差仿真分析结果进行了验证，测得不同工况下转台的指向误差低于0.16°。通过转台的高精度指向，太阳能够进入导行镜视场内，满足光谱仪系统的设计指标要求及在轨工作要求。 参考文献 李占峰 ， 黄煜 ， 王淑荣 ， 等 . 多通道共光栅转轴太阳光谱仪波长非线性研究 ［J］. 中国激光 ， 2019 ， 46 （ 9 ）： 0911004 . 多通道共光栅转轴太阳光谱仪波长非线性研究 LI ZH F ， HUANG Y ， WANG SH R ， et al . . Wavelength nonlinearity of multi-channel solar spectrometer with common grating rotating axis ［J］. Chinese Journal of Lasers ， 2019 ， 46 （ 9 ）： 0911004 . （in Chinese） . Wavelength nonlinearity of multi-channel solar spectrometer with common grating rotating axis 曹佃生 ， 林冠宇 ， 闻宝朋 ， 等 . 红外双光栅光谱仪波长扫描机构设计与精度分析 ［J］. 中国激光 ， 2017 ， 44 （ 5 ）： 0511002 . 红外双光栅光谱仪波长扫描机构设计与精度分析 CAO D SH ， LIN G Y ， WEN B P ， et al . . Wavelength scanning mechanism design and accuracy analysis of infrared double grating spectrometer ［J］. Chinese Journal of Lasers ， 2017 ， 44 （ 5 ）： 0511002 . （in Chinese） . Wavelength scanning mechanism design and accuracy analysis of infrared double grating spectrometer 李寒霜 . 紫外-真空紫外太阳光谱仪及光谱/辐射定标研究 ［D］. 长春 ： 中国科学院长春光学精密机械与物理研究所 ， 2019 . 紫外-真空紫外太阳光谱仪及光谱/辐射定标研究 LI H SH . Study on UV-VUV Solar Spectrometer and Spectral/Radiation Calibration ［D］. Changchun ： Changchun Institute of Optics ， Fine Mechanics and Physics ， Chinese Academy of Sciences ， 2019 . （in Chinese） Study on UV-VUV Solar Spectrometer and Spectral/Radiation Calibration 孙德贝 ， 李志刚 ， 李福田 . 用于太阳光谱仪的光电探测系统线性度测试装置 ［J］. 中国光学 ， 2019 ， 12 （ 2 ）： 294 - 301 . 用于太阳光谱仪的光电探测系统线性度测试装置 SUN D B ， LI ZH G ， LI F T . Linearity testing device for the photoelectric detecting system of solar spectrometers ［J］. Chinese Optics ， 2019 ， 12 （ 2 ）： 294 - 301 . （in Chinese） Linearity testing device for the photoelectric detecting system of solar spectrometers 李福田 ， 李志刚 ， 王晓旭 ， 等 . 自校准型光谱辐亮度标准光源 ［J］. 光学 精密工程 ， 2017 ， 25 （ 8 ）： 2004 - 2010 . 自校准型光谱辐亮度标准光源 LI F T ， LI ZH G ， WANG X X ， et al . . Self-calibrated spectral radiance standard source ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 8 ）： 2004 - 2010 . （in Chinese） . Self-calibrated spectral radiance standard source 李占峰 ， 王淑荣 ， 黄煜 ， 等 . 星载太阳辐照度光谱仪短时高精度太阳预报 ［J］. 光学学报 ， 2019 ， 39 （ 7 ）： 167 - 171 . 星载太阳辐照度光谱仪短时高精度太阳预报 LI ZH F ， WANG SH R ， HUANG Y ， et al . . High-precision and short-time solar forecasting by spaceborne solar irradiance spectrometer ［J］. Acta Optica Sinica ， 2019 ， 39 （ 7 ）： 167 - 171 . （in Chinese） . High-precision and short-time solar forecasting by spaceborne solar irradiance spectrometer 孙立微 . 空间高光谱遥感仪器辐射定标技术研究 ［D］. 长春 ： 中国科学院长春光学精密机械与物理研究所 ， 2018 . 空间高光谱遥感仪器辐射定标技术研究 SUN L W . Research on Radiometric Calibration for Space Hyperspectral Remote Sensor ［D］. Changchun ： Changchun Institute of Optics ， Fine Mechanics and Physics ， Chinese Academy of Sciences ， 2018 . （in Chinese） Research on Radiometric Calibration for Space Hyperspectral Remote Sensor 窦晨浩 . 空间大气/太阳光谱仪定量化关键技术研究 ［D］. 长春 ： 中国科学院长春光学精密机械与物理研究所 ， 2018 . 空间大气/太阳光谱仪定量化关键技术研究 DOU CH H . Research on the Key Quantitative Technology of Atmospheric/Solar Remote Sensing Instrument ［D］. Changchun ： Changchun Institute of Optics ， Fine Mechanics and Physics ， Chinese Academy of Sciences ， 2018 . （in Chinese） Research on the Key Quantitative Technology of Atmospheric/Solar Remote Sensing Instrument 王金元 ， 刘志惠 . 一种小型圆口径相控阵天线的指向误差分析 ［J］. 现代雷达 ， 2020 ， 42 （ 7 ）： 65 - 69 . 一种小型圆口径相控阵天线的指向误差分析 WANG J Y ， LIU ZH H . Analysis of beam-steering error of phased array antenna with small & circular shape ［J］. Modern Radar ， 2020 ， 42 （ 7 ）： 65 - 69 . （in Chinese） Analysis of beam-steering error of phased array antenna with small & circular shape 胡晓炜 ， 朱庆生 . 机架模型下指向误差数据的筛选与分析 ［J］. 天文研究与技术 ， 2021 ， 18 （ 1 ）： 77 - 86 . 机架模型下指向误差数据的筛选与分析 HU X W ， ZHU Q SH . Screening and analysis of pointing error data under frame model ［J］. Astronomical Research & Technology ， 2021 ， 18 （ 1 ）： 77 - 86 . （in Chinese） Screening and analysis of pointing error data under frame model 康跃然 ， 肖本龙 ， 傅亦源 ， 等 . 三轴转台系统姿态变化下指向误差建模与分析 ［J］. 制造业自动化 ， 2018 ， 40 （ 12 ）： 135 - 139 . 三轴转台系统姿态变化下指向误差建模与分析 KANG Y R ， XIAO B L ， FU Y Y ， et al . . Modeling and analyzing pointing error of the three-axis turntable system under posture change ［J］. Manufacturing Automation ， 2018 ， 40 （ 12 ）： 135 - 139 . （in Chinese） . Modeling and analyzing pointing error of the three-axis turntable system under posture change 吴伟平 ， 马庆军 ， 王淑荣 . 大气临边观测中视轴临边指向精度的在轨补偿 ［J］. 光学 精密工程 ， 2019 ， 27 （ 3 ）： 569 - 576 . 大气临边观测中视轴临边指向精度的在轨补偿 WU W P ， MA Q J ， WANG SH R . On-orbit correction method for limb pointing error of ozone detector ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 3 ）： 569 - 576 . （in Chinese） On-orbit correction method for limb pointing error of ozone detector LIAO T T . Modeling and analysis of laser shaft alignment using 4×4 homogeneous coordinate transformation matrix ［J］. Measurement ， 2009 ， 42 （ 1 ）： 157 - 163 . Modeling and analysis of laser shaft alignment using 4×4 homogeneous coordinate transformation matrix 胡春晖 ， 颜昌翔 . 基于高斯光学齐次坐标变换的光机装调 ［J］. 光学 精密工程 ， 2012 ， 20 （ 11 ）： 2353 - 2359 . 基于高斯光学齐次坐标变换的光机装调 HU CH H ， YAN CH X . Optical-mechanical assembly based on Gaussian optical homogeneous coordinate transformation ［J］. Opt. Precision Eng. ， 2012 ， 20 （ 11 ）： 2353 - 2359 . （in Chinese） Optical-mechanical assembly based on Gaussian optical homogeneous coordinate transformation LOU Y M ， HU J M . Coordinate transformation of three-dimensional image in integral imaging system illuminated by a point light source ［J］. Optik ， 2021 ， 226 ： 165885 . Coordinate transformation of three-dimensional image in integral imaging system illuminated by a point light source KULESZA J A ， Jr SOLOMON C J ， Jr KIEDROWSKI B C . Discrete ordinates analysis of the forced-flight variance reduction technique in Monte Carlo neutral particle transport simulations ［J］. Journal of Computational Physics ， 2021 ， 429 ： 109997 . Discrete ordinates analysis of the forced-flight variance reduction technique in Monte Carlo neutral particle transport simulations 10.37188/OPE.20212903.0474 V447.1 V448.222 A 1004-924X（2021）03-0474-10 郭旭,胡春晖,颜昌翔等.基于蒙特卡罗法的星载太阳辐照度光谱仪对日指向误差分析[J].光学精密工程,2021,29(03):474-483. 郭旭,胡春晖,颜昌翔等.基于蒙特卡罗法的星载太阳辐照度光谱仪对日指向误差分析[J].光学精密工程,2021,29(03):474-483. GUO Xu,HU Chun-hui,YAN Chang-xiang,et al.Analysis of sun pointing error of spaceborne solar spectroradiometer based on Monte Carlo method[J].Optics and Precision Engineering,2021,29(03):474-483. GUO Xu,HU Chun-hui,YAN Chang-xiang,et al.Analysis of sun pointing error of spaceborne solar spectroradiometer based on Monte Carlo method[J].Optics and Precision Engineering,2021,29(03):474-483. 国家自然科学基金资助项目（No.61805235） 2020-10-12 2020-11-11 2021-03-15 2021-04-08 光学精密工程 3 29