大連理工優(yōu)化方法增廣拉格朗日方法MATLAB程序.doc
《大連理工優(yōu)化方法增廣拉格朗日方法MATLAB程序.doc》由會員分享,可在線閱讀,更多相關(guān)《大連理工優(yōu)化方法增廣拉格朗日方法MATLAB程序.doc(5頁珍藏版)》請在裝配圖網(wǎng)上搜索。
上機(jī)大作業(yè)II定義目標(biāo)函數(shù)funfunction f=fun(x)x1=x(1);x2=x(2);f=4*x1-x22-12;定義目標(biāo)函數(shù)梯度函數(shù)dfunfunction f=dfun(x)x2=x(2);f=4;-2*x2;定義等式約束函數(shù)hffunction qua=hf(x)qua=25-x(1)2-x(2)2;定義等式約束函數(shù)梯度函數(shù)dhffunction qua=dhf(x)qua=-2*x(1);-2*x(2);定義不等式約束函數(shù)gfunfunction inq=gfun(x)inq=10*x(1)-x(1)2+10*x(2)-x(2)2-34;定義不等式約束梯度數(shù)dgffunction inq=dgf(x)inq=10-2*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*sigma*s1;s2=0.0;for i=1:m s3=max(0.0, lambda(i) - sigma*gi(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(:,i);endfor i=1:m dpsi=dpsi+(sigma*gi(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=1e-4; k=0;n=length(x0); Bk=eye(n);while(kmaxk) gk=feval(dmpsi,x0,fun,hf,gfun,dfun,dhf,dgf,mu,lambda,sigma); if(norm(gk)epsilon1) break; end dk=-Bkgk; m=0; mk=0; while(m20) newf=feval(mpsi,x0+rhom*dk,fun,hf,gfun,dfun,dhf,dgf,mu,lambda,sigma); oldf=feval(mpsi,x0,fun,hf,gfun,dfun,dhf,dgf,mu,lambda,sigma); if(newf0) Bk=Bk-(Bk*sk)*sk*Bk)/(sk*Bk*sk)+(yk*yk)/(yk*sk); end k=k+1; x0=x;endval=feval(mpsi,x0,fun,hf,gfun,dfun,dhf,dgf,mu,lambda,sigma); 定義增廣拉格朗日乘子法函數(shù)multphrfunction answer=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*ones(m,1);btak=10;btaold=10;while(btakepsilon&kepsilon if(k=2&btak theta*btaold) sigma=eta*sigma; end for i=1:l mu(i)=mu(i)-sigma*he(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運(yùn)行求解 x0=0;0x0 = 0 0 multphr(fun,hf,gfun,dfun,dhf,dgf,x0)x = 1.00128148956437 4.89871784708758f = -31.9923105871169mu = 1.01559644571312lambda = 0.754451167977228k = 4- 1.請仔細(xì)閱讀文檔,確保文檔完整性,對于不預(yù)覽、不比對內(nèi)容而直接下載帶來的問題本站不予受理。
- 2.下載的文檔,不會出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入(下載+內(nèi)容+預(yù)覽)歸上傳者、原創(chuàng)作者;如果您是本文檔原作者,請點(diǎn)此認(rèn)領(lǐng)!既往收益都?xì)w您。
下載文檔到電腦,查找使用更方便
9.9 積分
下載 |
- 配套講稿:
如PPT文件的首頁顯示word圖標(biāo),表示該P(yáng)PT已包含配套word講稿。雙擊word圖標(biāo)可打開word文檔。
- 特殊限制:
部分文檔作品中含有的國旗、國徽等圖片,僅作為作品整體效果示例展示,禁止商用。設(shè)計者僅對作品中獨(dú)創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞:
- 大連理工 優(yōu)化 方法 增廣 拉格朗日 MATLAB 程序
鏈接地址:http://ioszen.com/p-6663865.html