wolfe算法確定步長的最速下降法

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)數(shù)學(xué)與計(jì)算科學(xué)學(xué)院實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)項(xiàng)目名稱wolfe算法確定步長的最速下降法所屬課程名稱 最優(yōu)化方法 實(shí) 驗(yàn) 類 型 算法編程 實(shí) 驗(yàn) 日 期 2015.11.6 班 級(jí) 學(xué) 號(hào) 姓 名 成 績 一、實(shí)驗(yàn)概述：【實(shí)驗(yàn)?zāi)康摹?、使用非精確線性搜索wolfe算法確定步長。2、再使用最速下降法解決無約束最優(yōu)化問題?！緦?shí)驗(yàn)原理】1、最速下降法的基本思想是：從當(dāng)前點(diǎn)出發(fā)，取函數(shù)在點(diǎn)出下降最快的方向作為我們的搜索方向。由的Taylor展開式知道略去t的高階無窮小項(xiàng)不計(jì)，可見取時(shí)，

2、函數(shù)值下降的最多。2、Wolfe線性搜索準(zhǔn)則為：給定常數(shù)，取使得3、解無約束問題的最速下降法計(jì)算步驟：Step1：選取初始點(diǎn)，給定終止誤差Step2：計(jì)算Step3：Step4：利用wolfe線性搜索求得步長，令【實(shí)驗(yàn)環(huán)境】Windows7Visual C+ 6.0二、實(shí)驗(yàn)內(nèi)容：【實(shí)驗(yàn)方案】1、本次試驗(yàn)的題目為：用wolfe線性搜索下的最速下降法求解：2、根據(jù)最速下降法的求解步驟，在求解步長時(shí)，用wolfe線性搜索確定，運(yùn)用C語言編寫程序，在計(jì)算機(jī)上實(shí)現(xiàn)?！緦?shí)驗(yàn)過程】（實(shí)驗(yàn)步驟、記錄、數(shù)據(jù)、分析）（一）由題目可得， 。選取初值點(diǎn)則，因?yàn)閨d|0.1,所以根據(jù)wolfe搜索法，則檢驗(yàn)其是否滿足，

3、其中，若滿足，則，檢驗(yàn)2、若不滿足（1），（2），令 中使(1)成立的最大者，令。若滿足(2)則停止，并取得步長，3、否則，令，使(1)成立的最大者，令i=i+1,繼續(xù)2.4、根據(jù)算法編寫程序，編譯運(yùn)行，得出結(jié)果。【實(shí)驗(yàn)結(jié)論】（結(jié)果）【實(shí)驗(yàn)小結(jié)】（收獲體會(huì)）本次試驗(yàn)，加深了最速下降法的掌握。在求解步長時(shí)，除了可以運(yùn)用wolfe線性搜索求解，還可以運(yùn)用精確線性搜索、Aimijo搜索、wolf搜索等方法。本次試驗(yàn)在運(yùn)用wolfe搜索求解步長時(shí)，編程遇到了不少困難。在今后的學(xué)習(xí)中，要加強(qiáng)對(duì)算法的理解和算法實(shí)現(xiàn)。三、指導(dǎo)教師評(píng)語及成績：評(píng) 語評(píng)語等級(jí)優(yōu)良中及格不及格1.實(shí)驗(yàn)報(bào)告按時(shí)完成,字跡清楚,文字

4、敘述流暢,邏輯性強(qiáng)2.實(shí)驗(yàn)方案設(shè)計(jì)合理3.實(shí)驗(yàn)過程（實(shí)驗(yàn)步驟詳細(xì),記錄完整,數(shù)據(jù)合理,分析透徹）4實(shí)驗(yàn)結(jié)論正確. 成 績： 指導(dǎo)教師簽名： 批閱日期：附錄1：源 程 序#include#include#define e 0.1#define p 0.2#define q 0.5#define r 1#define s 0.5#define t 0.5float f(float x)float a;a=2*x0*x0+x1*x1;return a;float g1(float x)float g;g=4*x0;return g;float g2(float x)float g;g=2*x1;re

5、turn g;int Ok1(float x1,float x0,float a)float b,c;b=f(x1);c=f(x0)-p*a*(g1(x0)*g1(x0)+g2(x0)*g2(x0);if(b=c) return 1;else return 0;float qiubuchang(float x0,float d)int a0=1,a1=r,a,r1,r2;float x12;x10=x00+a*d0;x11=x01+a*d1;if(Ok1(x1,x0,a) & (Ok2(x1,x0,a) a=a0;elsedoa1*=s;x10=x00+a1*d0;x11=x01+a1*d1;

6、while(Ok1(x1,x0,a1);a0=a1;while(!Ok2(x1,x0,a0)r1=a0/r;dor2=r1-a0; a1=a0+r2;x10=x00+a1*d0;x11=x01+a1*d1;r2*=t;while(Ok1(x1,x0,a1);a0=a1;a=a0;return a;void main()int i=1;float a;float x02,x12,g2,d2;printf(輸入初始值：);scanf(%f,%f,&x00,&x01);printf(%f,%f,x00,x01);while(sqrt(double)(g1(x0)*(g1(x0)+(g2(x0)*(g2(x0)=e)printf(第%d次迭代：n,i+);a=qiubuchang(x0,d);printf(