(最新整理)大連理工優(yōu)化方法增廣拉格朗日方法MATLAB程序

1、(完整)大連理工優(yōu)化方法 增廣拉格朗日方法matlab程序(完整)大連理工優(yōu)化方法 增廣拉格朗日方法matlab程序 編輯整理：尊敬的讀者朋友們：這里是精品文檔編輯中心，本文檔內(nèi)容是由我和我的同事精心編輯整理后發(fā)布的，發(fā)布之前我們對文中內(nèi)容進行仔細校對，但是難免會有疏漏的地方，但是任然希望（(完整)大連理工優(yōu)化方法 增廣拉格朗日方法matlab程序）的內(nèi)容能夠給您的工作和學習帶來便利。同時也真誠的希望收到您的建議和反饋，這將是我們進步的源泉，前進的動力。本文可編輯可修改，如果覺得對您有幫助請收藏以便隨時查閱，最后祝您生活愉快 業(yè)績進步，以下為(完整)大連理工優(yōu)化方法 增廣拉格朗日方法matla

2、b程序的全部內(nèi)容。上機大作業(yè)ii定義目標函數(shù)funfunction f=fun（x)x1=x(1)；x2=x(2)；f=4x1-x22-12;定義目標函數(shù)梯度函數(shù)dfunfunction f=dfun（x)x2=x(2）;f=4;2*x2;定義等式約束函數(shù)hffunction qua=hf(x)qua=25x（1)2-x（2)2;定義等式約束函數(shù)梯度函數(shù)dhffunction qua=dhf（x)qua=2x（1);-2x（2）；定義不等式約束函數(shù)gfunfunction inq=gfun(x)inq=10*x(1)-x(1)2+10*x(2)-x（2)234;定義不等式約束梯度數(shù)dgffun

3、ction inq=dgf(x）inq=102*x（1）;10-2*x（2）;定義增廣拉格朗日函數(shù)mpsifunction psi=mpsi(x,fun，hf,gfun，dfun,dhf，dgf，mu，lambda，sigma)f=feval(fun，x)；he=feval(hf,x）；gi=feval(gfun,x）；l=length（he）;m=length（gi）；psi=f;s1=0;for i=1:l psi=psi-he（i）mu（i)； s1=s1+he（i)2;endpsi=psi+0.5*sigmas1;s2=0。0;for i=1：m s3=max(0。0， lambda(

4、i） sigmagi（i）； s2=s2+s32-lambda（i)2;endpsi=psi+s2/（2。0*sigma）;定義增廣拉格朗日函數(shù)梯度函數(shù)dmpsifunction dpsi=dmpsi(x,fun,hf，gfun,dfun,dhf,dgf，mu，lambda,sigma）dpsi=feval（dfun，x）;he=feval（hf,x);gi=feval(gfun,x）；dhe=feval（dhf，x）；dgi=feval(dgf,x）；l=length(he)；m=length（gi）;for i=1：l dpsi=dpsi+(sigma*he(i)-mu(i）)dhe（：,

5、i）;endfor i=1：m dpsi=dpsi+（sigmagi(i)-lambda(i)*dgi(：,i);end定義bfgs法函數(shù)函數(shù)bfgsfunction x，val，k=bfgs(mpsi，dmpsi,x0，fun，hf,gfun，dfun,dhf，dgf,mu,lambda，sigma）maxk=1000；rho=0.5;sigma1=0.4;epsilon1=1e4； k=0；n=length（x0)； bk=eye(n)；while(kmaxk) gk=feval(dmpsi，x0,fun，hf,gfun,dfun,dhf，dgf,mu,lambda,sigma）; if(

6、norm(gk）epsilon1) break; end dk=bkgk; m=0; mk=0; while（m20) newf=feval（mpsi，x0+rhomdk,fun,hf，gfun，dfun,dhf,dgf,mu，lambda,sigma）； oldf=feval（mpsi，x0,fun，hf,gfun，dfun,dhf,dgf，mu,lambda,sigma）; if（newfepsilon&k theta*btaold） sigma=eta*sigma； end for i=1：l mu(i)=mu（i)-sigmahe(i）； end for i=1：m lambda(i)=max（0.0，lambda（i)sigma*gi(i）); end end k=k+1； btaold=btak; x0=x;endf=feval（fun,x）；xfmulambdak運行求解 x0=0;0x0 = 0 0 multphr（fun,hf，gfun,dfun,dhf，dgf，x0）x = 1.0