发明名称 一种基于包络对偶算法体系的对地观测卫星多学科优化系统
摘要 一种基于包络对偶算法体系的对地观测卫星多学科优化系统,它涉及包络对偶算法体系的对地观测卫星多学科优化系统的5个分层:数学模型层L1,优化模型层L2,算法计算层L3,数据层L4和应用层L5;它们之间是信息传递关系,依次为L1‑L2‑L3‑L4‑L5;本发明是在对地观测卫星数学模型基础上,确定优化目标‑对地观测分辨率和相关约束,并构建卫星的优化模型,然后用包络对偶算法对卫星模型进行优化计算,如果优化结果满足收敛准则,则得到卫星最优的对地观测分辨率。它也可以实现对卫星设计目标和参数变量的快速计算,和传统的卫星设计‑试验‑更改设计‑试验的流程相比,提升了效率。
申请公布号 CN103761379B 申请公布日期 2016.08.17
申请号 CN201410016402.0 申请日期 2014.01.14
申请人 北京航空航天大学 发明人 王晓慧;徐岩
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种基于包络对偶算法体系的对地观测卫星多学科优化系统,其特征在于:它涉及包络对偶算法体系的对地观测卫星多学科优化系统的5个分层:数学模型层L1,优化模型层L2,算法计算层L3,数据层L4和应用层L5;它们之间是信息传递关系,依次为L1‑L2‑L3‑L4‑L5;所述数学模型层L1由任务模块L11,有效载荷模块L12,电源模块L13,发射模块L14四个模块组成;任务模块L11、有效载荷模块L12、电源模块L13和发射模块L14是相互并列的关系,通过卫星轨道高度和卫星相机孔径耦合;该任务模块L11,其结构是由下列经验公式构成:<maths num="0001"><math><![CDATA[<mrow><mi>D</mi><mi>L</mi><mi>D</mi><mo>=</mo><mfrac><mi>T</mi><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>&times;</mo><mn>2</mn><mi>&lambda;</mi><mo>&times;</mo><mfrac><mi>&pi;</mi><mn>180</mn></mfrac><mo>&times;</mo><mfrac><mn>1</mn><mn>60</mn></mfrac></mrow>]]></math><img file="FDA0000973699040000011.GIF" wi="454" he="103" /></maths><maths num="0002"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mn>2</mn><mi>&pi;</mi><msqrt><mfrac><msup><mrow><mo>(</mo><mn>6378</mn><mo>+</mo><mi>H</mi><mo>)</mo></mrow><mn>3</mn></msup><mn>39860</mn></mfrac></msqrt></mrow>]]></math><img file="FDA0000973699040000012.GIF" wi="373" he="134" /></maths><maths num="0003"><math><![CDATA[<mrow><mi>&lambda;</mi><mo>=</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mn>6378</mn><mrow><mn>6378</mn><mo>+</mo><mi>H</mi></mrow></mfrac><mo>&times;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mn>5</mn><mo>)</mo></mrow><mo>-</mo><mn>5</mn></mrow>]]></math><img file="FDA0000973699040000013.GIF" wi="574" he="111" /></maths>式中:DLD是下行链接持续时间,取值范围5分钟到15分钟;H是轨道高度,取值范围500公里‑750公里;T是卫星运行周期,是中间变量;λ是卫星覆盖区对应的地球中心角,是中间变量;卫星任务要考虑卫星与地面接收站的可见性及数据传输能力、卫星连续执行两次任务之间的调整校准时间,任务需求是下行链接持续时间即Down‑link data transmission;该有效载荷模块L12,其结构是由下列经验公式构成:<maths num="0004"><math><![CDATA[<mrow><mi>R</mi><mi>T</mi><mo>=</mo><mfrac><mrow><mn>6378</mn><mo>&times;</mo><mn>7.27</mn><mo>&times;</mo><mi>T</mi><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>5</mn></mrow></msup></mrow><mrow><mi>s</mi><mi>w</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000973699040000014.GIF" wi="458" he="103" /></maths><maths num="0005"><math><![CDATA[<mrow><mi>s</mi><mi>w</mi><mo>=</mo><mn>2</mn><mi>H</mi><mo>&times;</mo><mi>arc</mi><mi> </mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>&lsqb;</mo><mn>639</mn><mo>&times;</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi> </mi><mi>t</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mn>10</mn><mrow><mn>2000</mn><mo>&times;</mo><mn>8</mn><mo>&times;</mo><mi>D</mi></mrow></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000973699040000015.GIF" wi="813" he="117" /></maths>式中:RT是回访时间,取值范围90天‑150天;sw即swath width是行迹宽度,是中间变量;T是卫星运行周期,是中间变量;H是轨道高度,取值范围500公里‑750公里;D是相机孔径,取值范围50毫米‑150毫米;对地观测卫星有效载荷作用是地面图像的捕捉,所有的设计都是围绕有效载荷性能最大化,考虑回访时间即revisit time的影响作用;该电源模块L13,其结构是由下列经验公式构成:A<sub>sa</sub>=P<sub>sa</sub>/181.8<maths num="0006"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>s</mi><mi>a</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>a</mi><mi>v</mi></mrow></msub><mo>&times;</mo><mrow><mo>(</mo><mfrac><msub><mi>T</mi><mi>e</mi></msub><mn>0.6</mn></mfrac><mo>+</mo><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>e</mi></msub></mrow><mn>0.8</mn></mfrac><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>e</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000973699040000021.GIF" wi="574" he="118" /></maths><maths num="0007"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>a</mi><mi>v</mi></mrow></msub><mo>=</mo><mn>21</mn><msup><mi>e</mi><mrow><mo>(</mo><mn>0.006</mn><msub><mi>M</mi><mrow><mi>s</mi><mi>a</mi></mrow></msub><mo>)</mo></mrow></msup></mrow>]]></math><img file="FDA0000973699040000022.GIF" wi="350" he="85" /></maths><maths num="0008"><math><![CDATA[<mrow><msub><mi>T</mi><mi>e</mi></msub><mo>=</mo><mi>T</mi><mo>&times;</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mn>6378</mn><mrow><mn>6378</mn><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000973699040000023.GIF" wi="461" he="111" /></maths>式中:A<sub>sa</sub>是太阳帆板面积,取值小于1.2平方米;P<sub>sa</sub>是太阳能电池阵输出功率,是中间变量;T<sub>e</sub>是卫星地影时长,是中间变量;M<sub>sa</sub>是整星质量,取值小于200千克;T是卫星运行周期,是中间变量;H是轨道高度,取值范围500公里‑750公里;太阳帆板和蓄电池室是卫星电源设计的重要元素,为卫星在太空的运行提供动力;该发射模块L14,其结构是由下列经验公式构成:M<sub>sa</sub>+H‑720≤0式中:M<sub>sa</sub>是整星质量,取值小于200千克;H是轨道高度,取值范围500公里‑750公里;卫星发射与火箭组成的运载工具和用于地面支持的基础结构有关,在卫星设计过程中只考虑卫星质量与轨道高度的相互关系;所述优化模型层L2由两个模块组成,优化目标模块L21和约束模块L22;由数学模型层L1提供的经验公式,构建卫星的优化模型;该优化目标模块L21是对地观测卫星优化模型建立的关键模块,整个优化系统围绕优化目标进行;其中优化目标‑对地观测分辨率,指可以识别的最小地面距离和最小目标物的大小;该优化目标R由下列经验公式构成:<maths num="0009"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mrow><mn>10</mn><mo>&times;</mo><mi>H</mi></mrow><mrow><mn>8</mn><mo>&times;</mo><mi>D</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000973699040000024.GIF" wi="197" he="95" /></maths>式中:H是轨道高度,取值范围500公里‑750公里;D是相机孔径,取值范围50毫米‑150毫米;该约束模块L22作用是提供给优化过程限定条件;该限定条件包括整星质量取值小于200千克;太阳帆板面积,取值小于1.2平方米;回访时间,取值范围90天‑150天;下行链接持续时间,取值范围5分钟到15分钟;所述算法计算层L3核心是包络对偶算法计算模块,根据优化模型层L2提供的优化模型,利用包络对偶算法,对对地观测分辨率进行计算求解;该包络对偶算法,指针对地观测卫星设计优化问题,提出的一种求解方法,其含义表述如下:优化问题典型的数学模型为:{求目标函数f(X)的最小值,使得f(X)满足约束函数g<sub>i</sub>(X)≤0,(i=1,2,...,n)}根据该问题的数学意义,变换为如下的等价形式:{求目标函数f(X)的最小值,使得f(X)满足约束函数g<sub>max</sub>(X)≤0}式中,g<sub>max</sub>(X)等于g<sub>i</sub>(X)的最大值;基于经典的数学理论,进一步得到优化问题的数学模型:{求对偶函数d(λ)的最大值,对偶函数d(λ)={f(X)+λ·g<sub>max</sub>(X)}的最小值.}用一个光滑的包络函数E(x)替代约束函数g<sub>max</sub>(X),得到包络对偶问题数学模型:{求对偶函数d(λ)的最大值,对偶函数d(λ)={f(X)+λ·E(X)}的最小值.}该f(X)是目标函数,g<sub>max</sub>(X)是约束函数,d(λ)是对偶函数,E(x)是包络函数,λ是对偶变量;所述数据层L4的核心是数据分析模块,利用收敛准则,分析算法计算层L3给出的计算结果,符合收敛准则,则此结果为最终的优化结果,反之从优化模型层L2开始重新计算;该收敛准则是一种数学算法定义;所述应用层L5由软件组成,为用户提供一组卫星设计参数;该卫星设计参数,是指卫星设计制造过程中,所需的设计值。
地址 100191 北京市海淀区学院路37号