<p><font size="4">使用罚函数方法求这个模型的解时,需将上面的模型重新写为如下的无约束最小化形式: <br/><br/>min p(X,M)=c1x1+c2x2+c3x3+Mi[min(0,gi(x1,x2,x3,yl))]2 </font></p><p><font size="4">其中,g1?x1?x2?x3?yl =x1-yf_min<br/><br/>g2?x1?x2?x3?yl =yf_max-x1<br/><br/>g3?x1?x2?x3?yl =x2-sf_min<br/><br/>g4?x1?x2?x3?yl =sf_max-x2<br/><br/>g5?x1?x2?x3?yl =x3-mei_min<br/><br/>g6?x1?x2?x3?yl =mei_max-x3 <br/><br/>g7?x1?x2?x3?yl =yl-yl_min<br/><br/>g8?x1?x2?x3?yl =yl_max-yl</font></p><p><font size="4">Mi(i=1?2?…,8)是罚函数系数。</font></p><p><font size="4">优化模型?1 的求解步骤为:</font></p><p><font size="4">(1)取Mi?i=1?2?…,8 初始值为1000,允许误差为ε?k=1;</font></p><p><font size="4">(2)求无约束极值问题优化模型?2 的最优解;</font></p><p><font size="4">(3)若对某一个j1≤j≤8,有-gj(x1,x2,x3,yl≥ε,</font></p><p><font size="4">则:?</font></p><p><font size="4">Mk+1j=10×Mkj,令k=k+1,转第2步,否则停止迭代。</font></p> |