• 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
    • / 15
    • 下载费用:30 金币  

    重庆时时彩计划软件APP: 一种对空间合作目标跟踪成像的卫星姿态控制方法.pdf

    关 键 词:
    一种 空间 合作 目标 跟踪 成像 卫星 姿态 控制 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610498033.2

    申请日:

    2016.06.29

    公开号:

    CN106125745A

    公开日:

    2016.11.16

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G05D 1/08申请日:20160629|||公开
    IPC分类号: G05D1/08; G05D1/12 主分类号: G05D1/08
    申请人: 中国人民解放军国防科学技术大学
    发明人: 项军华; 张学阳; 李泰博; 曾国强; 张琦; 崔凯凯; 税海涛; 吴国福; 李志军; 袁福; 高玉东; 连一君; 韩大鹏
    地址: 410073 湖南省长沙市开福区德雅路109号
    优先权:
    专利代理机构: 北京中济纬天专利代理有限公司 11429 代理人: 陈立新
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610498033.2

    授权公告号:

    106125745B||||||

    法律状态公告日:

    2017.05.31|||2016.12.14|||2016.11.16

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明公开了一种对空间合作目标跟踪成像的卫星姿态控制方法,首先确定成像时间窗口需满足的条件,并由给定的卫星轨道状态和目标轨道信息计算期望姿态四元数和角速度,进而计算姿态跟踪误差四元数与实际角速度与期望角速度之间的误差量,然后建立姿态动力学模型,最后选取控制律参数,计算飞轮力矩控制量。该方法基于姿态的四元数表示,避免了大角度姿态运动时使用欧拉角可能出现的奇异问题。且该方法算法简单,运算量小,易于工程实现,具有较高的指向控制精度和较好的指向稳定度。

    权利要求书

    1.一种对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,包括以下步骤:
    S1、给定视频卫星与空间合作目标的轨道状态信息;
    S2、确定跟踪成像时间窗口的条件;
    S3、由给定的视频卫星与空间合作目标的轨道状态信息计算期望姿态四元数和角速
    度;
    S4、计算姿态跟踪误差四元数与实际角速度与期望角速度之间的误差量,然后建立姿
    态动力学模型;
    S5、选取控制律参数,计算飞轮力矩控制量。
    2.根据权利要求1所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S1中,给定的视频卫星与空间合作目标的轨道状态信息如下:
    地球惯性系下,视频卫星为S,空间合作目标为T,地心为O,令rst=
    rt-rs,rs=||rs||,rt=||rt||,rst=||rst||。
    3.根据权利要求2所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S2中,对于搭载可见光视频载荷的视频卫星,实现对空间合作目标的跟踪成像,需要满
    足以下三个条件:
    (1)视频卫星和空间合作目标可见;
    (2)距离要求;
    设视频卫星上的星载相机的焦距为f,像元大小为d,目标尺寸为L,像距为v,为使成像
    有意义,空间合作目标在像平面上的投影应超过一个像素,得到距离要求
    <mrow> <msub> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> </msub> <mo>&lt;</mo> <mfrac> <mrow> <mi>L</mi> <mi>f</mi> </mrow> <mi>d</mi> </mfrac> <mo>+</mo> <mi>f</mi> <mover> <mo>=</mo> <mo>&CenterDot;</mo> </mover> <mfrac> <mrow> <mi>L</mi> <mi>f</mi> </mrow> <mi>d</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    (3)光线要求
    为使空间合作目标良好成像,需要空间合作目标迎光且不处于地球的阴影区。
    4.根据权利要求3所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S2中,视频卫星和空间合作目标可见即地心O到视频卫星与空间合作目标连线ST的距
    离大于地球半径,则
    <mrow> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>r</mi> <mi>s</mi> </msub> <mo>&times;</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <msub> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> </msub> </mfrac> <mo>&gt;</mo> <mi>R</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    其中R为地球半径。
    5.根据权利要求3所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S2中,
    在地球惯性系下,认为太阳光是来自太阳-地心方向的平行光,设日地方向单位矢量为
    sun,则空间合作目标迎光即为在惯性系下,
    rst·sun>0 (3)
    空间合作目标不处于地球的阴影区需满足
    <mrow> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mfrac> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <mi>s</mi> <mi>u</mi> <mi>n</mi> </mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> </mfrac> <mo>&gt;</mo> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mi>R</mi> <msub> <mi>r</mi> <mi>t</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mi>&kappa;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    其中κ是太阳光下地球的阴影区域形成圆锥的半顶角,为0.264°。
    6.根据权利要求3所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S3包括以下步骤:
    S3.1计算期望姿态四元数
    定义视频卫星轨道坐标系为Oo-XoYoZo,它以视频卫星的质心为坐标系原点,OoXo轴沿着
    视频卫星的速度方向,OoZo轴从地心指向视频卫星质心;定义视频卫星本体坐标系为Ob-
    XbYbZb,它以视频卫星的质心为坐标系原点,三个坐标轴方向分别沿着视频卫星本体惯量主
    轴三个方向;定义参考空间合作目标坐标系为Ot-XtYtZt,它以视频卫星的质心为坐标系原
    点,以视频卫星的质心与目标点T的连线并且指向目标点的方向为坐标系的OtZt轴,以OtZt
    轴与视频卫星轨道坐标系OoZo的叉乘为参考空间合作目标坐标系的OtYt轴;
    假设星载相机固连在视频卫星上,星载相机的光轴与视频卫星本体坐标系的Z轴重合,
    于是跟踪成像就是控制视频卫星本体坐标系跟踪参考空间合作目标坐标系的过程,即视频
    卫星本体坐标系相对于地球惯性系的姿态四元数qb跟踪参考空间合作目标坐标系相对于
    地球惯性系的姿态四元数qt的过程;
    已知视频卫星的轨道倾角、升交点赤经和升交点赤经起算的轨道幅角分别为i,Ω,u,
    那么视频卫星轨道坐标系相对于地球惯性坐标系的方向余弦阵Roi可以表示为
    <mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>o</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>R</mi> <mi>Y</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>&Omega;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi> </mi> <mi>u</mi> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>cos</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi> </mi> <mi>u</mi> <mo>+</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>cos</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>cos</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi> </mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi> </mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <mi>i</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>u</mi> <mo>+</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>sin</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>u</mi> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>sin</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi> </mi> <mi>i</mi> <mi> </mi> <mi>sin</mi> <mi> </mi> <mi>u</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
    其中RX(θ),RY(θ),RZ(θ)分别表示绕X轴,Y轴,Z轴旋转角度θ的坐标变换矩阵;
    参考空间合作目标坐标系的Z轴在地球惯性系下的单位矢量为
    <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <mfrac> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>o</mi> </msubsup> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>i</mi> </mrow> <mi>o</mi> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
    其中的分母表示长度,分子是黑体,表示矢量,二者相除得到单位矢量,且:
    Rio=RoiT (7)
    参考空间合作目标坐标系的Y轴在地球惯性系下的单位矢量为
    <mrow> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
    其中是视频卫星轨道坐标系Z轴在地球惯性系下的单位矢量:
    <mrow> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
    参考空间合作目标坐标系的X轴在地球惯性系下的单位矢量为
    <mrow> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
    于是得到参考空间合作目标坐标系相对于地球惯性系的方向余弦阵为Rti
    <mrow> <msub> <mi>R</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mo>&lsqb;</mo> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&rsqb;</mo> <mo>=</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mrow> <mn>3</mn> <mo>&times;</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
    其中rij表示矩阵元素,矩阵是3*3矩阵,i和j分别表示矩阵中的行和列,共9个元素;
    式(11)就是跟踪目标成像的期望姿态矩阵;
    令期望姿态四元数qt=(q1:3,q4)T,其中q1:3=(q1,q2,q3)T,则
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>11</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>1</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>22</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>2</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>33</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>3</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>4</mn> </msub> <mi>q</mi> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
    对上式取模,即可得到期望姿态四元数qt;
    S3.2计算期望姿态角速度
    由姿态运动学方程得到期望姿态角速度:
    <mrow> <msub> <mi>&omega;</mi> <mi>t</mi> </msub> <mo>=</mo> <mn>2</mn> <msup> <mi>&Xi;</mi> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <msub> <mover> <mi>q</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mn>4</mn> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>3</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
    7.根据权利要求6所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S4包括以下步骤:
    S4.1计算姿态跟踪误差四元数
    令视频卫星本体坐标系相对于地球惯性系的姿态四元数其中qb1:3=
    (qb1,qb2,qb3)T,则姿态跟踪误差四元数
    <mrow> <mi>&delta;</mi> <mi>q</mi> <mo>=</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&CircleTimes;</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>=</mo> <mo>&lsqb;</mo> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>&rsqb;</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>4</mn> </mrow> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>4</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
    S4.2计算实际角速度与期望角速度之间的误差量
    ωe=ωb-ωt (17)
    其中,ωb为卫星当前角速度,ωt为期望角速度;
    S4.3得到卫星姿态动力学方程
    <mrow> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>u</mi> <mo>+</mo> <msub> <mi>T</mi> <mi>d</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
    其中J为卫星转动惯量,ωb是视频卫星本体坐标系的角速度,h为飞轮的角动量,u为控
    制力矩,Td为外部干扰力矩,且
    <mrow> <mi>u</mi> <mo>=</mo> <mo>-</mo> <mover> <mi>h</mi> <mo>&CenterDot;</mo> </mover> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
    8.根据权利要求7所述的对空间合作目标跟踪成像的卫星姿态控制方法,其特征在于,
    步骤S5中选取控制参数,计算飞轮力矩控制量u
    <mrow> <mi>u</mi> <mo>=</mo> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <msub> <mi>&delta;q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> <msub> <mi>&omega;</mi> <mi>e</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
    其中Kp和Kd为常值正定矩阵。

    说明书

    一种对空间合作目标跟踪成像的卫星姿态控制方法

    技术领域

    本发明涉及一种航天航空领域的卫星姿态控制方法,它为视频卫星对空间合作目
    标跟踪成像提供一种姿态控制方法,属于自动控制技术领域。

    背景技术

    视频卫星是一种采用视频成像、视频数据实时传输、人在回路交互式操作工作方
    式的新型天基信息获取类微小卫星。与传统卫星相比,视频卫星可提供实时视频图像,与静
    止单幅图像相比增加了时域信息,可获取目标的动态过程信息,能够探测到动态事件的发
    生,并可以基于视频图像中的序列图像进行图像重构获得更高分辨率的图像,为抗灾救灾、
    战时监控、计划决策提供第一手资料。

    空间合作目标是指已知运动轨道信息的空间飞行器,利用视频卫星对空间合作目
    标跟踪成像,相当于在太空架起了摄像机,不受气象条件的限制,且避免了大气对目标信号
    的影响,实时性更好,可以为卫星健康监测、故障诊断提供有力支持,未来甚至能直播太空
    的重要事件。

    已有文献对视频卫星对地面固定目标的凝视成像姿态控制方法进行了较为深入
    的研究,但未有涉及对空间合作目标的跟踪成像姿态控制。

    发明内容

    针对视频卫星对空间合作目标的跟踪成像问题,本发明提供一种对空间合作目标
    跟踪成像的卫星姿态控制方法。该方法首先给出了成像的时间窗口需满足的条件,以及期
    望姿态四元数与期望角速度的求解方法,建立了其姿态运动的数学模型;以此模型为受控
    对象,设计了PD控制器。本发明所提出的跟踪成像姿态控制器结构框图如图1所示。实践表
    明,由该方法控制的闭环系统能够稳定指向空间合作目标,且具有良好的指向精度和稳定
    度。

    一种对空间合作目标跟踪成像的卫星姿态控制方法,首先确定成像时间窗口需满
    足的条件,并由给定的卫星轨道状态和目标轨道信息计算期望姿态四元数和角速度,进而
    计算姿态跟踪误差四元数与实际角速度与期望角速度之间的误差量,然后建立姿态动力学
    模型,最后选取控制律参数,计算飞轮力矩控制量。实际应用中,卫星姿态四元数和角速度
    由姿态确定系统测量得到,将由该方法计算得到的控制量传输至执行机构即可实现对空间
    合作目标跟踪成像姿态控制功能。

    一种对空间合作目标跟踪成像的卫星姿态控制方法,包括以下步骤:

    步骤一:给定视频卫星与空间合作目标的轨道状态信息

    地球惯性系下,视频卫星为S,空间合作目标为T,地心为O,令rst
    =rt-rs,rs=||rs||,rt=||rt||,rst=||rst||。

    步骤二:确定跟踪成像时间窗口的条件

    对于搭载可见光视频载荷的视频卫星,实现对空间合作目标的跟踪成像,需要满
    足三个条件:

    (1)卫星和目标可见;

    视频卫星和空间合作目标可见,即地心O到视频卫星与空间合作目标连线ST的距
    离大于地球半径,则

    <mrow> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>r</mi> <mi>s</mi> </msub> <mo>&times;</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <msub> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> </msub> </mfrac> <mo>&gt;</mo> <mi>R</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>

    其中R为地球半径。

    (2)距离要求;

    设视频卫星上的星载相机的焦距为f,像元大小为d,目标尺寸为L,像距为v,为使
    成像有意义,空间合作目标在像平面上的投影应超过一个像素,得到距离要求

    <mrow> <msub> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> </msub> <mo>&lt;</mo> <mfrac> <mrow> <mi>L</mi> <mi>f</mi> </mrow> <mi>d</mi> </mfrac> <mo>+</mo> <mi>f</mi> <mover> <mo>=</mo> <mo>&CenterDot;</mo> </mover> <mfrac> <mrow> <mi>L</mi> <mi>f</mi> </mrow> <mi>d</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

    (3)光线要求。

    为使空间合作目标良好成像,需要空间合作目标迎光且不处于地球的阴影区。在
    地球惯性系下,可以认为太阳光是来自太阳-地心方向的平行光,设日地方向单位矢量为
    sun,则空间合作目标迎光即为在惯性系下,

    rst·sun>0 (3)

    空间合作目标不处于地球的阴影区域需满足

    <mrow> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mfrac> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <mi>s</mi> <mi>u</mi> <mi>n</mi> </mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> </mfrac> <mo>&gt;</mo> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mi>R</mi> <msub> <mi>r</mi> <mi>t</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mi>&kappa;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

    其中κ是太阳光下地球的阴影区域形成圆锥的半顶角,为0.264°。

    步骤三:期望姿态四元数与期望角速度计算

    1)计算期望姿态四元数

    定义视频卫星轨道坐标系为Oo-XoYoZo,它以视频卫星的质心为坐标系原点,OoXo轴
    沿着视频卫星的速度方向,OoZo轴从地心指向视频卫星质心;定义视频卫星本体坐标系为
    Ob-XbYbZb,它以视频卫星的质心为坐标系原点,三个坐标轴方向分别沿着视频卫星本体惯
    量主轴三个方向;定义参考空间合作目标坐标系为Ot-XtYtZt,它以视频卫星的质心为坐标
    系原点,以视频卫星的质心与目标点T的连线并且指向目标点的方向为坐标系的OtZt轴,以
    OtZt轴与视频卫星轨道坐标系OoZo的叉乘为参考空间合作目标坐标系的OtYt轴。以上坐标系
    都为右手系。

    假设星载相机固连在视频卫星上,星载相机的光轴与视频卫星本体坐标系的Z轴
    重合,于是跟踪成像就是控制视频卫星本体坐标系跟踪参考空间合作目标坐标系的过程,
    即视频卫星本体坐标系相对于地球惯性系的姿态四元数qb跟踪参考空间合作目标坐标系
    相对于地球惯性系的姿态四元数qt的过程。

    已知视频卫星的轨道倾角、升交点赤经和升交点赤经起算的轨道幅角分别为i,
    Ω,u,那么视频卫星轨道坐标系相对于地球惯性坐标系的方向余弦阵Roi可以表示为

    <mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>o</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>R</mi> <mi>Y</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>&Omega;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>u</mi> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>u</mi> <mo>+</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>i</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>u</mi> <mo>+</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>u</mi> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

    其中RX(θ),RY(θ),RZ(θ)分别表示绕X轴,Y轴,Z轴旋转角度θ的坐标变换矩阵。

    参考空间合作目标坐标系的Z轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <mfrac> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>o</mi> </msubsup> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>i</mi> </mrow> <mi>o</mi> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

    其中的分母表示长度,分子是黑体,表示矢量,二者相除得到单位矢量,且:

    Rio=RoiT (7)

    参考空间合作目标坐标系的Y轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

    其中是视频卫星轨道坐标系Z轴在地球惯性系下的单位矢量:

    <mrow> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>

    参考空间合作目标坐标系的X轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>

    于是得到参考空间合作目标坐标系相对于地球惯性系的方向余弦阵为Rti

    <mrow> <msub> <mi>R</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mo>&lsqb;</mo> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&rsqb;</mo> <mo>=</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mrow> <mn>3</mn> <mo>&times;</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>

    其中rij表示矩阵元素,矩阵是3*3矩阵,i和j分别表示矩阵中的行和列,共9个元
    素,也是本领域常用表示方法

    这就是跟踪目标成像的期望姿态矩阵。

    令期望姿态四元数qt=(q1:3,q4)T,其中q1:3=(q1,q2,q3)T,则

    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>11</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>1</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>22</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>2</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>33</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>3</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>4</mn> </msub> <mi>q</mi> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>

    对上式取模,即可得到期望姿态四元数qt。

    2)计算期望姿态角速度

    由姿态运动学方程得到期望姿态角速度:

    <mrow> <msub> <mi>&omega;</mi> <mi>t</mi> </msub> <mo>=</mo> <mn>2</mn> <msup> <mi>&Xi;</mi> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <msub> <mover> <mi>q</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>

    其中

    <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mn>4</mn> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>3</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>

    步骤四:建立姿态动力学模型

    1)计算姿态跟踪误差四元数

    令视频卫星本体坐标系相对于地球惯性系的姿态四元数qb=(qb1:3,qb4)T,其中
    qb1:3=(qb1,qb2,qb3)T,则姿态跟踪误差四元数

    <mrow> <mi>&delta;</mi> <mi>q</mi> <mo>=</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&CircleTimes;</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>=</mo> <mo>&lsqb;</mo> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>&rsqb;</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>

    其中

    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>4</mn> </mrow> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>4</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>

    2)计算实际角速度与期望角速度之间的误差量

    ωe=ωb-ωt (17)

    其中,ωb为卫星当前角速度,ωt为期望角速度。

    3)得到卫星姿态动力学方程

    <mrow> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>u</mi> <mo>+</mo> <msub> <mi>T</mi> <mi>d</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>

    其中J为卫星转动惯量,ωb是视频卫星本体坐标系的角速度,h为飞轮的角动量,u
    为控制力矩,Td为外部干扰力矩,且

    <mrow> <mi>u</mi> <mo>=</mo> <mo>-</mo> <mover> <mi>h</mi> <mo>&CenterDot;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>

    步骤五:控制律设计:选取控制参数(Kp和Kd),计算飞轮力矩控制量u

    <mrow> <mi>u</mi> <mo>=</mo> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <msub> <mi>&delta;q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> <msub> <mi>&omega;</mi> <mi>e</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>

    其中Kp和Kd为常值正定矩阵。

    本发明的有益效果是:

    1)该方法基于姿态的四元数表示,避免了大角度姿态运动时使用欧拉角可能出现
    的奇异问题。

    2)该方法算法简单,运算量小,易于工程实现。

    3)该方法具有较高的指向控制精度和较好的指向稳定度。

    控制工程师在应用过程中可以指定空间合作目标,并将由该方法得到的控制量传
    输至执行机构实现跟踪成像姿态控制功能。

    附图说明

    图1为本发明所述卫星跟踪成像姿态控制系统结构

    图2为本发明所述卫星跟踪成像控制方法步骤

    图3为本发明所述欧拉角控制误差变化曲线

    图4为本发明所述角速度控制误差变化曲线

    图5为本发明所述飞轮控制力矩变化曲线。

    具体实施方式

    以下将结合具体实施例和说明书附图对本发明做进一步详细说明。

    一种对空间合作目标跟踪成像的卫星姿态控制方法,其具体步骤如下:

    步骤一:给定视频卫星和空间合作目标的轨道根数

    视频卫星将采用某型号视频卫星设计参数,目标则选为天宫一号设计参数,仿真
    起始时间为20160501UTC,它们的轨道根数如表1所示。

    表1某型号视频卫星与天宫一号的轨道根数


    某型号视频卫星
    天宫一号
    半长轴(km)
    6821.235
    6767.416
    偏心率
    0.000689
    0.001972
    轨道倾角(°)
    97.314
    42.794
    升交点赤经(°)
    201.542
    168.994
    近地点幅角(°)
    345.682
    51.774
    平近点角(°)
    349.531
    289.998

    由轨道根数可以解算出卫星与空间合作目标的运动状态。

    步骤二:确定成像时间窗口

    某型号视频卫星搭载的高清相机焦距为1000mm,像元尺寸8.33μm,天宫一号最大
    直径3.35m,代入式(2)得到对天宫一号的可成像距离小于402km。

    利用时间窗口的三个条件,在一周内可以得到的跟踪成像时间窗口如表2所示。

    表2跟踪成像时间窗口

    访问次数
    起始时间
    终止时间
    持续时间(s)
    1
    1May 2016 00:59:20.721
    1May 2016 01:00:54.408
    93.688
    2
    4May 2016 05:50:49.191
    4May 2016 05:52:08.380
    79.189
    3
    7May 2016 09:08:31.097
    7May 2016 09:08:48.445
    17.348
    4
    7May 2016 10:41:17.614
    7May 2016 10:42:16.530
    58.915

    仿真时间选为1May 2016 00:57:20到1May 2016 01:02:55,即在第一个时间窗口
    前后各加120s。

    步骤三:计算期望姿态四元数与期望角速度

    1)计算期望姿态四元数

    已知视频卫星的轨道倾角、升交点赤经和升交点赤经起算的轨道幅角分别为i,
    Ω,u,那么视频卫星轨道坐标系相对于地球惯性坐标系的方向余弦阵Roi可以表示为

    <mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>o</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>R</mi> <mi>Y</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>X</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>Z</mi> </msub> <mrow> <mo>(</mo> <mi>&Omega;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>u</mi> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>u</mi> <mo>+</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mi>i</mi> <mi>cos</mi> <mi>u</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mi>&Omega;</mi> <mi>sin</mi> <mi>i</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>i</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>u</mi> <mo>+</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>u</mi> <mo>-</mo> <mi>cos</mi> <mi>&Omega;</mi> <mi>cos</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mi>i</mi> <mi>sin</mi> <mi>u</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>

    其中RX(θ),RY(θ),RZ(θ)分别表示绕X轴,Y轴,Z轴旋转角度θ的坐标变换矩阵。

    参考空间合作目标坐标系的Z轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <mfrac> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>o</mi> </msubsup> <msubsup> <mi>r</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>o</mi> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>

    其中

    Rio=RoiT (23)

    参考空间合作目标坐标系的Y轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>

    其中是视频卫星轨道坐标系Z轴在地球惯性系下的单位矢量:

    <mrow> <msubsup> <mi>Z</mi> <mi>o</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>o</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>

    参考空间合作目标坐标系的X轴在地球惯性系下的单位矢量为

    <mrow> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>

    于是得到参考空间合作目标坐标系相对于地球惯性系的方向余弦阵为Rti

    <mrow> <msub> <mi>R</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mo>&lsqb;</mo> <msubsup> <mi>X</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Y</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>,</mo> <msubsup> <mi>Z</mi> <mi>t</mi> <mi>i</mi> </msubsup> <mo>&rsqb;</mo> <mo>=</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mrow> <mn>3</mn> <mo>&times;</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>

    这就是跟踪目标成像的期望姿态矩阵。

    令期望姿态四元数qt=(q1:3,q4)T,其中q1:3=(q1,q2,q3)T,则

    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>11</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>1</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>22</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>2</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>13</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>31</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mn>33</mn> </msub> <mo>-</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>3</mn> </msub> <mi>q</mi> <mo>,</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>31</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>13</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mn>12</mn> </msub> <mo>-</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>trR</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mn>4</mn> <msub> <mi>q</mi> <mn>4</mn> </msub> <mi>q</mi> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>

    对上式取模,即可得到期望姿态四元数qt。

    2)计算期望角速度

    由姿态运动学方程得到期望姿态角速度:

    <mrow> <msub> <mi>&omega;</mi> <mi>t</mi> </msub> <mo>=</mo> <mn>2</mn> <msup> <mi>&Xi;</mi> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <msub> <mover> <mi>q</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>

    其中

    <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mn>4</mn> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>3</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>

    步骤四:计算飞轮控制力矩u

    1)计算姿态跟踪误差四元数

    <mrow> <mi>&delta;</mi> <mi>q</mi> <mo>=</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&CircleTimes;</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>=</mo> <mo>&lsqb;</mo> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>&rsqb;</mo> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>

    其中

    <mrow> <mi>&Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>4</mn> </mrow> </msub> <msub> <mi>I</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>&lsqb;</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>&times;</mo> <mo>&rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mrow> </mtd> <mtd> <msub> <mi>q</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>

    <mrow> <msubsup> <mi>q</mi> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>q</mi> <mn>4</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>32</mn> <mo>)</mo> </mrow> </mrow>

    2)计算实际角速度与期望角速度之间的误差量

    ωe=ωb-ωt (33)

    其中,ωb为卫星当前角速度,ωt为期望角速度。

    3)得到闭环控制系统

    将式(20)代入式(18),忽略外部干扰力矩Td,得到闭环控制系统方程

    <mrow> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>e</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <msub> <mi>&delta;q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> <msub> <mi>&omega;</mi> <mi>e</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>

    某型号视频卫星转动惯量为控制器的参
    数选为Kp=diag(0.5,0.8,0.4),Kd=diag(5,8,4),状态初值如表3所示。

    表3状态变量初值

    变量
    数值
    变量
    数值
    q1
    0.064496
    ωx
    3.5806348009×10-5rad/s
    q2
    -0.608277
    ωy
    -0.0010860287rad/s
    q3
    0.186158
    ωz
    -0.0002351377rad/s
    q4
    0.768886

    4)飞轮控制力矩为:

    <mrow> <mi>u</mi> <mo>=</mo> <mi>J</mi> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <msub> <mi>&delta;q</mi> <mrow> <mn>1</mn> <mo>:</mo> <mn>3</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> <msub> <mi>&omega;</mi> <mi>e</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>

    实施例中的视频卫星对空间合作目标跟踪成像姿态控制结果如图3~图5所示。为
    了更加直观,通过将四元数转换成欧拉角,图3给出了欧拉角控制误差曲线,图4给出了角速
    度控制误差曲线,由图3与图4可得:在60s内可使视频卫星光轴从初始状态指向目标,在交
    会时控制误差变大,但仍然有欧拉角控制精度<0.3°,角速度控制精度<0.03°/s,满足成像
    要求。图5给出了飞轮控制力矩变化曲线,已知飞轮的最大控制力矩为20mN·m,由图5可得:
    所提出的控制方法不会造成飞轮过载。

    以上包含了本发明优选实施例的说明,这是为了详细说明本发明的技术特征,并
    不是想要将发明内容限制在实施例所描述的具体形式中,依据本发明内容主旨进行的其他
    修改和变型也受本专利?;?。本发明内容的主旨是由权利要求书所界定,而非由实施例的
    具体描述所界定。

       内容来自专利网重庆时时彩单双窍门 www.4mum.com.cn转载请标明出处

    关于本文
    本文标题:一种对空间合作目标跟踪成像的卫星姿态控制方法.pdf
    链接地址://www.4mum.com.cn/p-5842706.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • .江苏十一选五开奖号码 6场半全场随机 爱玩棋牌斗牛游戏 广东十一选五 7星彩开奖结果查询1 美容保健赚钱吗 以太坊公司如何赚钱 大乐透预测139期历史同期 14场胜负 重庆幸运农场几点开始 体彩新疆11选5开奖结果 甘肃快3开奖 英雄杀卡牌 湖北30选5最新开奖结果查询 山西快乐十分玩法 云南快乐十分50