惩罚函数的乘子法.doc

《惩罚函数的乘子法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2013-2014（1）专业课程实践论文题目：惩罚函数的乘子法一、算法理论乘子法是Powell和Hestenes于1969年针对等式约束优化问题同时独立提出的一种优化算法,后于1973年经Rockfellar推广到求解不等式约束优化问题。其基本思想是从原问题的拉格朗日函数出发,再加上适当的罚函数,从而将原问题转化为求解一系列的无约束优化子问题。由于外罚函数法中的罚参数,因此增广目标函数变得“越来越病态”。增广目标函数的这种病态性质是外罚函数法的主要缺点,而这种缺陷在乘子法中由于引入拉格朗日函数及加上适当的罚函数而得以有效的克服。我们考虑同时带有等

2、式和不等式约束的优化问题的乘子法:其基本思想是把解等式约束优化问题的乘子法推广到不等式约束优化问题,即先引进辅助变量把不等式约束化为等式约束,然后再利用最优性条件消去辅助变量。为叙述的方便计,我们先考虑如下只带有不等式约束的最优化问题引进辅助变量，可以将上面的优化问题化为等价的等式约束优化问题：利用外发函数法求解，此时增广拉格朗日函数为为了消去辅助变量，可考虑关于变量的极小化，由一阶必要条件，令可得即故当时，有否则，由可推得。综合起来。，有既有（1）因此，当时，我们有而当时，有综合上述两种情形，将结果代回到中去得于是，将式（1）带入乘子迭代公式得

3、即回到一般约束优化问题，此时，增广拉格朗日函数为乘子迭代的公式为令则终止准则为二、算法框图三、算法程序importjava.awt.*;importjava.awt.event.*;importjavax.swing.*;importjava.text.NumberFormat;publicclassShuextendsFrameimplementsActionListener{Labellabelx12,labelx1,labelx22,labelx2,labelx1x2,labela,labelb,labelc,label11,label12

4、,labelnum,labelnumber;TextFieldtextx12,textx1,textx22,textx2,textc1,textx1x2,texta,textb,textc,textnum,textnumber;Buttonbutton;TextAreatextarea;//10doublex1,x12,x22,x2,x1x2,c1,a,b,c,num,number;doublea1x1,a1x2,a1c,a2x1,a2x2,a2c;intn=1;doubleans[],answ[];publicShu(){ans=newdoub

5、le[20];answ=newdouble[20];labelx12=newLabel("x1^2+");labelx1=newLabel("x1+");labelx22=newLabel("x2^2+");//20labelx2=newLabel("x2+");labelx1x2=newLabel("x1*x2+");labela=newLabel("x1+");labelb=newLabel("x2+");labelc=newLabel("=0");label11=newLabel("min");label12=newLabel("s.t."

6、);textx12=newTextField(3);textx1=newTextField(3);textx22=newTextField(3);//30textx2=newTextField(3);textx1x2=newTextField(3);textc1=newTextField(3);texta=newTextField(3);textb=newTextField(3);textc=newTextField(3);labelnumber=newLabel("δ:");labelnum=newLabel("μ:");textnum=new

7、TextField(3);textnumber=newTextField(3);button=newButton("enter");button.addActionListener(this);textarea=newTextArea(10,10);Boxbox1=Box.createHorizontalBox();//44box1.add(label11);box1.add(textx12);box1.add(labelx12);box1.add(textx1);box1.add(labelx1);box1.add(textx22);box1.

8、add(labelx22);box1.add(textx2);box1.add(labelx2);box1.add(textx1x2);