发明名称 一种飞行器系统余度配置优化方法
摘要 本发明提供了一种飞行器系统余度配置优化方法,首先确定机载设备设计目标及约束条件,并固化为数学模型或仿真模型,模型输入参数需要包括不同机载设备或其内部部件的余度参数;然后把每个设计目标模型作为NSGA‑II算法的适应度函数,每个设计约束条件作为NSGA‑II算法的输入,余度参数作为决策变量;之后对决策变量进行编码,判断约束条件是否满足,满足则计算适应度函数,对种群进行选择、交叉、变异并进行非支配排序;最后进行拥挤度计算,直到达到终止条件,此时获得最优余度配置向量集合,通过用户依据实际需要进行选取,最终获得最优余度配置方案。
申请公布号 CN106599411A 申请公布日期 2017.04.26
申请号 CN201611075903.1 申请日期 2016.11.30
申请人 中国航空工业集团公司沈阳飞机设计研究所 发明人 费思邈;孙国强;宁彬;王斌
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京航信高科知识产权代理事务所(普通合伙) 11526 代理人 周良玉
主权项 一种飞行器系统余度配置优化方法,所述优化方法基于NSGA‑II算法,其特征在于,包括以下步骤:步骤一,确定机载设备的设计目标及约束条件,将所述设计目标固化为数学模型或仿真模型,所述数学模型或仿真模型的输入参数包括不同机载设备的余度参数或该机载设备内部组成单元的余度参数;步骤二,将每个所述数学模型或仿真模型均作为NSGA‑II算法的适应度函数,将每个所述约束条件作为NSGA‑II算法的输入,将所述余度参数作为决策变量;步骤三,对决策变量进行编码,然后判断决策变量的约束条件是否满足,所属约束条件包括第一约束条件,第一约束条件为余度参数是整数,若满足约束条件则计算每个适应度函数;其中适应度函数即优化目标,有时候优化目标可能有多个,也就是多个适应度函数,每个适应度函数都需要计算一遍;步骤四,对个体进行选择、交叉、变异,并进行非支配排序;选择的标准是选择通过适应度函数算出的优化目标的最大值;交叉过程中的交叉算子采用拉普拉斯算子,通过公式(1)对交叉过程进行控制,<maths num="0001"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>a</mi><mo>-</mo><mi>b</mi><mi> </mi><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>a</mi><mo>+</mo><mi>b</mi><mi> </mi><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&gt;</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001166077810000011.GIF" wi="1366" he="289" /></maths> 公式(1)中β<sub>i</sub>为随机变量,a为位置参数,b为尺度参数,μ<sub>i</sub>为服从0~1均匀分布的随机数,将β<sub>i</sub>代入公式(2)中,通过公式(2)对个体进行交叉,<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>y</mi><mi>i</mi><mn>1</mn></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>i</mi><mn>1</mn></msubsup><mo>+</mo><msub><mi>&beta;</mi><mi>i</mi></msub><mo>|</mo><msubsup><mi>x</mi><mi>i</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mi>i</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>&beta;</mi><mi>i</mi></msub><mo>|</mo><msubsup><mi>x</mi><mi>i</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup><mo>|</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001166077810000021.GIF" wi="1247" he="397" /></maths> 公式(2)中x是父代,y是子代,xi是父代决策变量,yi是子代决策变量,xi和yi的上标1和上标2代表从父代个体种群中任意选择的两个中的第一个个体和第二个个体;变异过程中的变异算子采用幂函数算子,通过公式(3)进行变异,<maths num="0003"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>s</mi><mrow><mo>(</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>-</mo><msup><mi>x</mi><mi>l</mi></msup><mo>)</mo></mrow><mo>,</mo><mi>t</mi><mo>&lt;</mo><mi>r</mi></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>+</mo><mi>s</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>u</mi></msup><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>,</mo><mi>t</mi><mo>&GreaterEqual;</mo><mi>r</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001166077810000022.GIF" wi="1116" he="322" /></maths>公式(3)中x<sup>u</sup>是决策变量的上界,x<sup>l</sup>是决策变量的下界,r是服从0~1均匀分布的随机变量,t为决策变量均值与极大极小值的关系,并且<img file="FDA0001166077810000023.GIF" wi="323" he="215" />s为利用幂函数生成的随机变量s=(s<sub>1</sub>)<sup>p</sup>,其中p为变异强度,且p为整数,s1为调节变异强度的随机数;步骤五,若未达到终止条件,则对种群进行拥挤度计算,选择拥挤度满足设定标准的个体,然后返回步骤三,若达到终止条件则得到最优余度配置集合; 步骤六,通过F值评判标准对最优余度配置集合进行评判并选出F值最高的余度配置方案,F值评判标准的公式如公式(4)所示,<maths num="0004"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfrac><mrow><munderover><mo>&Pi;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>G</mi><mi>i</mi></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>G</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001166077810000031.GIF" wi="735" he="344" /></maths> 公式(4)中F值为余度值,Gi为第i个优化目标经过正则化后的值,n为优化目标的总数。
地址 110035 辽宁省沈阳市皇姑区塔湾街40号