... 389alpha=alpha0;g=grad(f,x);s=‐g;forn=1:nalpha=alpha0;fx1=feval(f,x+alpha*2*s);forn1=1:nmax1fx2=fx1;fx1=feval(f,x+alpha*s);if(fx0>fx1+tolfun)&(fx1<fx2‐tolfun)%fx0>fx1<fx2den=4*fx1‐2*fx0‐2*fx2;num=den‐fx0+fx2;alpha=alpha*num/den;x=x+alpha*s;fx=feval(f,x);break;elseifn1==nmax1/2alpha=‐alpha0;fx1=feval(f,x+alpha*2*s);elsealpha=alpha/2;endendx0=x;fx0=fx;ifn<ng1=grad(f,x,h);ifKC<=1s=‐g1+(g1‐g)*g1ʹ/(g*gʹ+1e‐5)*s;elses=‐g1+g1*g1ʹ/(g*gʹ+1e‐5)*s;endg=g1;endifn1>=nmax1warning=warning+1;%kgtimduockichthuoctoiuu 370CHƯƠNG 8: TỐI ƯU HOÁ §1.KHÁINIỆMCHUNGVỀTỐI ƯU HOÁ Tối ưu hoálàthuậtngữthườngđượcdùngđểcựctiểuhoáhaycựcđạihoámộthàm.Thôngthườngtachỉcầntìmcựctiểu ... 0.618033989 2Chúýlàmỗilầnthugọnkhoảngchứađiểmcựctiểuthìkhoảng(a,b)giảmtỉlệvớir. Điều nàylàmsốlầntínhlớnhơnphươngphápchiađôi.Tuynhiênphươngpháptỉlệvàngchỉđòihỏitínhgiátrịhàmmộtlầntrongkhiphươngphápchiađôicầntínhgiátrịhàm2lần.Sốlầntínhxácđịnhbằng:−=εnbar ... Trongchươngnàychúngtasẽlầnlượtxétcácthuậttoántìmcựctiểukhôngràngbuộcvàcóràngbuộc.§2.PHƯƠNGPHÁPTIẾTDIỆNVÀNG Taxétbàitoántìmcựctiểucủahàmmộtbiếnf(x).Điểmcựctiểuđượcxácđịnhtheo điều kiệndf/dx=0.Docóthểcónhiềuđiểmcựctiểunêntaphảivâyđiểmcựctiểu(xácđịnhlâncậnchứađiểmcựctiểu)trước.Thủthuậtvâyđiểmcựctiểukháđơngiản:chođiểmđầux0vàtínhgiátrịcủahàmđangđi...