发明名称 一种求解全局优化问题的正交逐次逼近方法
摘要 本发明提出了一种基于正交试验设计和变尺度邻域搜索的正交逐次逼近方法,包括以下步骤:设置初始化参数,由所求问题维度和离散水平数选择合适的正交表;由初始点x<sup>0</sup>出发在可行域内开展正交试验,采用惩罚函数法评估计算各试验方案;从试验方案中选择新的迭代点,若x<sup>1</sup>优于x<sup>0</sup>,则令x<sup>0</sup>=x<sup>1</sup>,同时扩大搜索步长增强全局搜索,否则缩小搜索空间加强局部寻优;重复上述步骤,反复迭代,逐次逼近全局最优解,直至满足收敛条件。本发明提出一种求解全局优化问题的正交逐次逼近方法,具有原理简单、计算参数少、收敛速度快等优点,可快速获取全局优化问题的最优解或其近似解。
申请公布号 CN103984675A 申请公布日期 2014.08.13
申请号 CN201410191150.5 申请日期 2014.05.06
申请人 大连理工大学 发明人 程春田;冯仲恺;廖胜利;牛文静;武新宇;李刚;申建建
分类号 G06F17/14(2006.01)I;G06F17/30(2006.01)I 主分类号 G06F17/14(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 赵连明;梅洪玉
主权项 1.一种求解全局优化问题的正交逐次逼近方法,其特征包括如下步骤,步骤1:设定初始点<img file="FDA0000500995790000011.GIF" wi="432" he="78" />初始搜索步长<img file="FDA0000500995790000012.GIF" wi="433" he="78" />终止精度ξ,邻域离散水平数t,扩张因子a<sub>1</sub>(a<sub>1</sub>≥1)、压缩因子a<sub>2</sub>(0<a<sub>2</sub><1)及最大迭代次数<img file="FDA00005009957900000110.GIF" wi="73" he="67" />置计数器I=1;步骤2:由离散水平数t及维度n选择相应正交表;步骤3:根据正交表构造以x<sup>0</sup>为中心的正交方案,并检查各因素水平是否越限,若越限则取为边界值;构造公式如下:<img file="FDA0000500995790000013.GIF" wi="443" he="148" /><maths num="0001"><![CDATA[<math><mrow><mmultiscripts><mi>x</mi><mi>i</mi><none/><mprescripts/><none/><mi>r</mi></mmultiscripts><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mi>i</mi></msub></mtd><mtd><mi>if</mi><mrow><mo>(</mo><mmultiscripts><mi>x</mi><mi>i</mi><none/><mprescripts/><none/><mi>r</mi></mmultiscripts><mo><</mo><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mi>i</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mmultiscripts><mi>x</mi><mi>i</mi><none/><mprescripts/><none/><mi>r</mi></mmultiscripts></mtd><mtd><mi>if</mi><mrow><mo>(</mo><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mi>i</mi></msub><mo>&le;</mo><mmultiscripts><mi>x</mi><mi>i</mi><none/><mprescripts/><none/><mi>r</mi></mmultiscripts><mo>&le;</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></mtd><mtd><mi>if</mi><mrow><mo>(</mo><mmultiscripts><mi>x</mi><mi>i</mi><none/><mprescripts/><none/><mi>r</mi></mmultiscripts><mo>></mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中<img file="FDA0000500995790000015.GIF" wi="68" he="69" />为因素i的r水平相应取值;<img file="FDA00005009957900000111.GIF" wi="91" he="69" />表示不大于x的最大整数;<img file="FDA0000500995790000016.GIF" wi="104" he="68" />分别为因素i的上限与下限;步骤4:利用惩罚函数法计算各试验方案指标值,并从中选择较优方案x<sup>1</sup>,然后采用极差分析法统计获取较优解x',若F(x')<F(x<sup>1</sup>),则令x<sup>1</sup>=x';<maths num="0002"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>M</mi><mi>i</mi></msub><mo>&times;</mo><mi>max</mi><mo>{</mo><mn>0</mn><mo>,</mo><msubsup><mi>g</mi><mi>i</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>}</mo><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>C</mi><mi>j</mi></msub><mo>&times;</mo><msubsup><mi>h</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中M<sub>i</sub>为第i项不等式约束惩罚系数;C<sub>j</sub>为第j项等式约束惩罚系数;步骤5:若F(x<sup>1</sup>)≤F(x<sup>0</sup>),则令x<sup>0</sup>=x<sup>1</sup>,同时增大步长,h<sub>0</sub>=a<sub>1</sub>×h<sub>0</sub>,转至步骤3;否则缩小步长h<sub>0</sub>=a<sub>2</sub>×h<sub>0</sub>,转至步骤6;步骤6:I=I+1;步骤7:若<img file="FDA0000500995790000018.GIF" wi="116" he="65" />或<img file="FDA0000500995790000019.GIF" wi="462" he="95" />则转至步骤8;否则转步骤3;步骤8:停止计算,输出最优解x<sup>*</sup>=x<sup>0</sup>。
地址 116024 辽宁省大连市甘井子区凌工路2号