发明名称 一种圆参考轨道下给定边界的卫星初始相对状态确定方法
摘要 本发明提供一种圆参考轨道下给定边界的卫星初始相对状态确定方法,给出了分别适用于周期性相对运动和绕飞相对运动的优化模型,然后采用一定的优化算法求解所建立的优化模型,得到给定边界的初始相对状态数值解,并通过仿真验证了该卫星初始相对状态确定方法的有效性和可行性,因此,能够快速精确地确定符合给定边界的卫星初始相对状态,同时,为具有给定边界的卫星编队或集群任务设计与分析奠定了理论基础。
申请公布号 CN104200030B 申请公布日期 2017.03.29
申请号 CN201410450128.8 申请日期 2014.09.05
申请人 清华大学 发明人 张育林;王兆魁;党朝辉;张学阳;安梅岩;蒋超;刘红卫;侯振东;张斌斌;程子龙
分类号 G06F17/50(2006.01)I;G06N3/12(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京市盛峰律师事务所 11337 代理人 席小东
主权项 一种圆参考轨道下给定边界的卫星初始相对状态确定方法,其特征在于,所述圆参考轨道是指参考星所在轨道的偏心率为零的轨道;包括以下步骤:S1,建立用于确定周期性相对运动条件下给定边界的初始相对状态的第一类优化模型,该第一类优化模型的目标函数为:<maths num="0001"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>X</mi></munder><mi>J</mi><mo>=</mo><mo>|</mo><mo>|</mo><msubsup><mi>d</mi><mi>min</mi><mrow><mi>r</mi><mi>e</mi><mi>q</mi></mrow></msubsup><mo>-</mo><msub><mi>d</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mo>|</mo><mo>|</mo><msubsup><mi>d</mi><mi>max</mi><mrow><mi>r</mi><mi>e</mi><mi>q</mi></mrow></msubsup><mo>-</mo><msub><mi>d</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0001165389090000011.GIF" wi="630" he="79" /></maths>该第一类优化模型的约束条件为:s.t.<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mi>n</mi></mfrac><mo>,</mo><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>2</mn><mi>n</mi></mfrac><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mo>+</mo><mn>3</mn><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>c</mi><mn>4</mn></msub><mo>=</mo><msub><mi>y</mi><mn>0</mn></msub><mo>-</mo><mfrac><mn>2</mn><mi>n</mi></mfrac><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>c</mi><mn>5</mn></msub><mo>=</mo><mfrac><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mi>n</mi></mfrac><mo>,</mo><msub><mi>c</mi><mn>6</mn></msub><mo>=</mo><msub><mi>z</mi><mn>0</mn></msub><mo>,</mo><msub><mi>c</mi><mn>3</mn></msub><mo>=</mo><mn>2</mn><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mfrac><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mi>n</mi></mfrac><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>g</mi><mn>1</mn></msub><mo>=</mo><mn>4</mn><msubsup><mi>c</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>4</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>6</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><msub><mi>c</mi><mn>1</mn></msub><msub><mi>c</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>g</mi><mn>2</mn></msub><mo>=</mo><mn>4</mn><mrow><mo>(</mo><mo>-</mo><mn>3</mn><msub><mi>c</mi><mn>1</mn></msub><msub><mi>c</mi><mn>2</mn></msub><mo>+</mo><mn>2</mn><msub><mi>c</mi><mn>2</mn></msub><msub><mi>c</mi><mn>4</mn></msub><mo>-</mo><msub><mi>c</mi><mn>5</mn></msub><msub><mi>c</mi><mn>6</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>g</mi><mn>3</mn></msub><mo>=</mo><mn>2</mn><mrow><mo>(</mo><mo>-</mo><mn>2</mn><msubsup><mi>c</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><mn>7</mn><msubsup><mi>c</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>4</mn><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><msubsup><mi>c</mi><mn>5</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>c</mi><mn>6</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>g</mi><mn>4</mn></msub><mo>=</mo><mn>4</mn><mrow><mo>(</mo><mn>3</mn><msub><mi>c</mi><mn>1</mn></msub><msub><mi>c</mi><mn>2</mn></msub><mo>+</mo><mn>2</mn><msub><mi>c</mi><mn>2</mn></msub><msub><mi>c</mi><mn>4</mn></msub><mo>+</mo><msub><mi>c</mi><mn>5</mn></msub><msub><mi>c</mi><mn>6</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>g</mi><mn>5</mn></msub><mo>=</mo><mn>4</mn><msubsup><mi>c</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>4</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>c</mi><mn>6</mn><mn>2</mn></msubsup><mo>+</mo><mn>4</mn><msub><mi>c</mi><mn>1</mn></msub><msub><mi>c</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>S</mi><mo>=</mo><mo>{</mo><msub><mi>s</mi><mi>i</mi></msub><mo>&Element;</mo><mo>|</mo><msub><mi>g</mi><mn>2</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>4</mn></msubsup><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>g</mi><mn>3</mn></msub><mo>-</mo><mn>2</mn><msub><mi>g</mi><mn>1</mn></msub><mo>)</mo></mrow><msubsup><mi>s</mi><mi>i</mi><mn>3</mn></msubsup><mo>+</mo><mn>3</mn><mrow><mo>(</mo><msub><mi>g</mi><mn>4</mn></msub><mo>-</mo><msub><mi>g</mi><mn>2</mn></msub><mo>)</mo></mrow><msubsup><mi>s</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>g</mi><mn>5</mn></msub><mo>-</mo><msub><mi>g</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>s</mi><mi>i</mi></msub><mo>-</mo><msub><mi>g</mi><mn>4</mn></msub><mo>=</mo><mn>0</mn><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mi>max</mi></msub><mo>=</mo><munder><mi>max</mi><mrow><msub><mi>s</mi><mi>i</mi></msub><mo>&Element;</mo><mi>S</mi></mrow></munder><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msubsup><mi>s</mi><mi>i</mi><mn>2</mn></msubsup></mrow></mfrac><msqrt><mrow><msub><mi>g</mi><mn>1</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>4</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>3</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>3</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>4</mn></msub><msub><mi>s</mi><mi>i</mi></msub><mo>+</mo><msub><mi>g</mi><mn>5</mn></msub></mrow></msqrt></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mi>min</mi></msub><mo>=</mo><munder><mi>min</mi><mrow><msub><mi>s</mi><mi>i</mi></msub><mo>&Element;</mo><mi>S</mi></mrow></munder><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msubsup><mi>s</mi><mi>i</mi><mn>2</mn></msubsup></mrow></mfrac><msqrt><mrow><msub><mi>g</mi><mn>1</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>4</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>3</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>3</mn></msub><msubsup><mi>s</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>g</mi><mn>4</mn></msub><msub><mi>s</mi><mi>i</mi></msub><mo>+</mo><msub><mi>g</mi><mn>5</mn></msub></mrow></msqrt></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001165389090000012.GIF" wi="1382" he="815" /></maths>其中,J为目标函数;<img file="FDA0001165389090000013.GIF" wi="72" he="55" />—期望的星间距离上界;<img file="FDA0001165389090000014.GIF" wi="70" he="55" />—期望的星间距离下界;d<sub>max</sub>—实际的星间距离上界;d<sub>min</sub>—实际的星间距离下界;<img file="FDA0001165389090000015.GIF" wi="402" he="61" />—待求解的跟踪星初始相对状态向量;x<sub>0</sub>—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对位置在x轴上的分量;y<sub>0</sub>—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对位置在y轴上的分量;z<sub>0</sub>—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对位置在z轴上的分量;<img file="FDA0001165389090000016.GIF" wi="41" he="53" />—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对速度在x轴上的分量;<img file="FDA0001165389090000021.GIF" wi="44" he="47" />—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对速度在y轴上的分量;<img file="FDA0001165389090000022.GIF" wi="36" he="47" />—跟踪星在参考星当地轨道坐标系下,相对参考星的初始相对速度在z轴上的分量;n—参考星的轨道角速度,<img file="FDA0001165389090000023.GIF" wi="214" he="87" />其中,μ为地球引力常数;a为参考星的轨道半长轴,是已知量;r<sub>0</sub>—参考星的矢径,由于参考星运行在圆轨道上,则:r<sub>0</sub>=a;S2,以<img file="FDA0001165389090000024.GIF" wi="403" he="62" />作为设计变量,以<img file="FDA0001165389090000025.GIF" wi="236" he="61" />n、r<sub>0</sub>作为已知输入量,按照一定的优化算法对各个设计变量进行优化,最终得到符合期望星间距离上界和期望星间距离下界的卫星初始相对状态。
地址 100084 北京市海淀区清华园1号