###############Goldstein方法#########################
建立一個minGS.m的文件,如下:

function [x,minf] = minGS(f,XMAX,sigma1,sigma2,alpha,eps) %目標函數：f； %搜索最大值：XMAX； %可接受系數 1：sigma1； %可接受系數 2：sigma2； %增大探索點系數：alpha %精度：eps； %目標函數取最小值時的自變量值：x； %目標函數的最小值：minf format long; if nargin == 5 eps = 1.0e-6; endif sigma1<=0 || sigma1>1 %輸入參數檢測disp('sigma1 參數不對！'); x = NaN; minf = NaN; return; else if sigma2 <= sigma1 disp('sigma2 參數不對！'); x = NaN; minf = NaN; return; else if alpha <= 1 disp('alpha 參數不對！'); x = NaN; minf = NaN; return; end end end df = diff(f); %一階導數 f0 = subs(f, findsym(f),0); df0 = subs(df, findsym(df),0); a = 0; b = XMAX; k = 0; t = (a+b)/2; %初始探測點while 1 ft = subs(f,findsym(f),t); %探測點的函數值f1 = f0 + sigma1*t*df0; %可接受函數值上限if ft <= f1 f2 = f0 + sigma2*t*df0; %可接受函數值下限if ft >= f2 x = t; break; else a = t; %更新極值點所在區間的左端點if b < XMAX t = (a+b)/2; %更新探測點else t = alpha*t; end end else b = t; %更新極值點所在區間的右端點t = (a+b)/2; end k = k+1; end minf = subs(f,findsym(f),x); format short;

運行方法如下:
>> syms t;
>>> f = t^2-2*t+7;
>>> x = minGS(f, 8, 0.4,0.75,2)
然后運行即可.

####################################Wolfe-Powell 法 ########################################
新建一個minWP.m,代碼如下:

function [x,minf] = minWP(f,XMAX,sigma1,sigma2,alpha,eps) %目標函數：f； %搜索最大值：XMAX； %可接受系數 1：sigma1； %可接受系數 2：sigma2； %增大探索點系數：alpha %精度：eps； %目標函數取最小值時的自變量值：x； %目標函數的最小值：minf format long; if nargin == 5 eps = 1.0e-6; end if sigma1<=0 || sigma1>1 %輸入參數檢測disp('sigma1 參數不對！'); x = NaN; minf = NaN; return; else if sigma2 <= sigma1 disp('sigma2 參數不對！'); x = NaN; minf = NaN; return; else if alpha <= 1disp('alpha 參數不對！');x = NaN;minf = NaN;return;endend enddf = diff(f); f0 = subs(f, findsym(f),0); df0 = subs(df, findsym(df),0); a = 0; b = XMAX; k = 0; t = (a+b)/2; %初始探測點 while 1 ft = subs(f,findsym(f),t); %探測點的函數值f1 = f0 + sigma1*t*df0; %可接受函數值上限if ft <= f1 f2 = sigma2*t*df0; %可接受函數值下限if ft >= f2 x = t; break; else a = t; %更新極值點所在區間的左端點if b < XMAX t = (a+b)/2; %更新探測點else t = alpha*t; endendelseb = t; %更新極值點所在區間的右端點t = (a+b)/2; end k = k+1; end minf = subs(f,findsym(f),x); format short;

運行方法如下:
syms=t;
f=t^2-2*t+7;
x1=minWP(f,8,0.1,0.65,2)

以上代碼來自《精通MATLAB最優化計算》(第二版)-龔純 王正林
這本書是下面這個鏈接的博主在評論區把百度網盤中的資源中分享給我的,在此表示感謝.
https://blog.csdn.net/u012366767/article/details/81563555

總結

以上是生活随笔為你收集整理的Armijo-Goldstein和wolfe-power的matlab代码实现(转)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。