大连理工优化方法 增广拉格朗日方法MATLAB程序

《大连理工优化方法 增广拉格朗日方法MATLAB程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、上机大作业II定义目标函数funfunctionf=fun(x)x1=x(1);x2=x(2);f=4*x1-x2^2-12;定义目标函数梯度函数dfunfunctionf=dfun(x)x2=x(2);f=[4;-2*x2];定义等式约束函数hffunctionqua=hf(x)qua=25-x(1)^2-x(2)^2;定义等式约束函数梯度函数dhffunctionqua=dhf(x)qua=[-2*x(1);-2*x(2)];定义不等式约束函数gfunfunctioninq=gfun(x)inq=10*x(1)-x(1)^2+10*

2、x(2)-x(2)^2-34;定义不等式约束梯度数dgffunctioninq=dgf(x)inq=[10-2*x(1);10-2*x(2)];定义增广拉格朗日函数mpsifunctionpsi=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;fori=1:lpsi=psi-he(i)*mu(i);s1=s1+he(i)^2

3、;endpsi=psi+0.5*sigma*s1;s2=0.0;fori=1:ms3=max(0.0,lambda(i)-sigma*gi(i));s2=s2+s3^2-lambda(i)^2;endpsi=psi+s2/(2.0*sigma);定义增广拉格朗日函数梯度函数dmpsifunctiondpsi=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)

4、;dgi=feval(dgf,x);l=length(he);m=length(gi);fori=1:ldpsi=dpsi+(sigma*he(i)-mu(i))*dhe(:,i);endfori=1:mdpsi=dpsi+(sigma*gi(i)-lambda(i))*dgi(:,i);end定义BFGS法函数函数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;e

5、psilon1=1e-4;k=0;n=length(x0);Bk=eye(n);while(k

6、un,dfun,dhf,dgf,mu,lambda,sigma);if(newf0)Bk=Bk-((Bk*sk)*sk'*Bk)/(sk'*Bk*sk)+(yk*yk')/(yk'*sk);endk=k+1;x0=x;endval=feval(mpsi,

7、x0,fun,hf,gfun,dfun,dhf,dgf,mu,lambda,sigma);定义增广拉格朗日乘子法函数multphrfunctionanswer=multphr(fun,hf,gfun,dfun,dhf,dgf,x0)maxk=5000;sigma=2.0;eta=2.0;theta=0.8;k=0;ink=0;epsilon=1e-4;x=x0;he=feval(hf,x);gi=feval(gfun,x);l=length(he);m=length(gi);mu=0.1*ones(l,1);lambda=0.1*one

8、s(m,1);btak=10;btaold=10;while(btak>epsilon&&k