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

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

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

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

4、敘述流暢,邏輯性強(qiáng)2.實驗方案設(shè)計合理3.實驗過程（實驗步驟詳細(xì),記錄完整,數(shù)據(jù)合理,分析透徹）4實驗結(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(