发明名称 针对单星任务规划问题的数学模型的构建方法及求解方法
摘要 本发明提供一种针对单星任务规划问题的数学模型的构建方法及求解方法,针对单星多圈的成像卫星任务规划问题,对问题进行了合理简化与假设,并构建数学模型,介绍了软件求解框架和9大分功能模块,其中排序模块采用了启发式的排序策略,活动安排模块创新性的采用了时间窗消减算法,较好的解决了时间窗的冲突问题,提升了时间窗的利用率,使得可行解接近最优。
申请公布号 CN104050324A 申请公布日期 2014.09.17
申请号 CN201410282628.5 申请日期 2014.06.23
申请人 中国人民解放军国防科学技术大学 发明人 邢立宁;张雪婷;义余江;刘嵩;杨振宇;刘晓路;王沛;朱剑冰;郭坚;汪路元
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 一种针对单星任务规划问题的数学模型的构建方法,其特征在于,包括以下几个步骤:步骤S101:设定以下输入条件:(1)每次进行规划的任务总数量为N,回传窗口的数目为N<sub>BW</sub>;(2)每个任务T<sub>i</sub>的优先级为p<sub>i</sub>,需要的观测时间为Vt<sub>i</sub>,回传时间为Bt<sub>i</sub>,任务的容量为C<sub>i</sub>,其中:i=1,2,…N;(3)每个任务观测窗口VW<sub>i</sub>的时间窗为[VS<sub>i</sub>,VE<sub>i</sub>],该任务只能在其观测时间窗内安排观测活动,其中:i=1,2,…N;(4)回传窗口BW<sub>j</sub>的时间窗为[BS<sub>j</sub>,BE<sub>j</sub>],只能在该回传窗口的时间窗内安排回传活动,其中:j=1,2,…N<sub>BW</sub>;(5)卫星所能承载的最大载荷量C<sub>max</sub>,t时刻卫星的承载容量为C<sub>t</sub>;步骤S102:设定以下输出条件:(1)决策变量:<img file="FDA0000525186430000011.GIF" wi="666" he="148" /><img file="FDA0000525186430000012.GIF" wi="725" he="152" />单个任务的收益值为分段函数,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>pf</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0.6</mn><msub><mi>&omega;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>&le;</mo><msub><mi>t</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>t</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0.8</mn><msub><mi>&omega;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>&le;</mo><msub><mi>t</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>t</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mn>1.0</mn><msub><mi>&omega;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mn>2</mn></msub><mo>&le;</mo><msub><mi>t</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>t</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><mn>0.8</mn><msub><mi>&omega;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mn>3</mn></msub><mo>&le;</mo><msub><mi>t</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>t</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><mn>0.6</mn><msub><mi>&omega;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mn>4</mn></msub><mo>&le;</mo><msub><mi>t</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>t</mi><mn>5</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000525186430000013.GIF" wi="455" he="362" /></maths>其中[t<sub>0</sub>,t<sub>5</sub>]为任务收益函数pf<sub>i</sub>的定义域,t<sub>0</sub>=VS<sub>i</sub>,t<sub>5</sub>=VE<sub>i</sub>,即每个观测时间窗被均匀划分为5段;(2)N个任务的排序结果:以N×4的矩阵形式输出,其中第一列表示任务的ID序号T<sub>i</sub>,第二列表示该任务实际开始观测时刻RVS<sub>i</sub>,第三列表示该任务实际开始回传时刻RBS<sub>i</sub>,第四列表示该任务的收益pf<sub>i</sub>;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><msub><mi>RVS</mi><mn>1</mn></msub></mtd><mtd><msub><mi>RBS</mi><mn>1</mn></msub></mtd><mtd><msub><mi>pf</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><msub><mi>RVS</mi><mn>2</mn></msub></mtd><mtd><msub><mi>RBS</mi><mn>2</mn></msub></mtd><mtd><msub><mi>pf</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>N</mi></msub></mtd><mtd><msub><mi>RVS</mi><mi>N</mi></msub></mtd><mtd><msub><mi>RBS</mi><mi>N</mi></msub></mtd><mtd><msub><mi>pf</mi><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000525186430000021.GIF" wi="527" he="290" /></maths>步骤S103:将该单星多圈的成像卫星任务规划问题的优化目标设置为在保证成像质量的前提下,实现以尽可能少的资源消耗实现最大的任务收益,即以完成观测和回传的任务总收益作为衡量求解结果的度量指标,目标函数为:<img file="FDA0000525186430000022.GIF" wi="461" he="128" />步骤S104:设定以下约束条件:(1)在安排观测任务时,任务的总容量大小不能超过卫星的最大载荷量:C<sub>t</sub>≤C<sub>max</sub>;(2)每个任务只能在其观测时间窗内安排观测活动;设任务T<sub>i</sub>实际开始观测的时刻为RVS<sub>i</sub>:VS<sub>i</sub>≤RVS<sub>i</sub>≤RVS<sub>i</sub>+Vt<sub>i</sub>≤VE<sub>i</sub>,i=1,2,...,N;(3)所有任务只能在回传时间窗内安排回传活动;设任务T<sub>i</sub>实际开始回传的时刻为RBS<sub>i</sub>:BS<sub>j</sub>≤RBS<sub>i</sub>≤RBS<sub>i</sub>+Bt<sub>i</sub>≤BE<sub>j</sub>,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>,</mo><mo>&Exists;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>N</mi><mi>BW</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000525186430000023.GIF" wi="617" he="74" /></maths>(4)每个任务只有在安排了观测活动之后才能安排回传活动,且每个任务最多只能安排一次观测和一次回传:y<sub>i</sub>≤x<sub>i</sub>≤1,i=1,2,...,N;RBV<sub>i</sub>+Vt<sub>i</sub>≤RBS<sub>i</sub>,i=1,2,...,N。
地址 410073 湖南省长沙市开福区砚瓦池正街47号