下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、共扼梯度法程序源代碼#include#iiiclude#define N 10#define eps pow( 10.-6)double f(double x,double p,double t)double s;s=pow(x0+t*p 0 ,2)+2 5 *pow(x 1 +t*p 1 ,2);leturn s;/*以下是進(jìn)退法搜索區(qū)間源程序*/void sb(double *a,double *b.double x,double p)double tO,tl,t,h,alpha,ft),fl;mt k=0;t0=2.5;/* 初始值 */h=l;/*初始步長(zhǎng)*/alpha=2;/*加步系
2、數(shù)*/fO=f(x,ptO);tl=tO+h;fl=f(x,p,tl);wlule(l)(h=alpha*h; t=tO;tO=tl;fO=f 1;k+;)else(if(k=0)h=-h;t=tl;elsebreak;)tl=tO+h;fl=Rx,p,tl);/*以下是黃金分割法程序源代碼可double lijfg(double x.double p)double beta,tl,t2,t;doubledouble a=0,b=0;double *c,*d;c=&a,d=&b;sb(c,d,x,p);/*調(diào)用進(jìn)退法搜索區(qū)間*/pnntf(”nxl=%lf.x2=%lf.pl=%lf,p2=%
3、lT,x0,xl,p0,pl);printf(na,b=%lf,%lf”,a,b);beta=(sqrt(5)-l .0)/2;t2=a+beta*(b-a); f2=f(x.p,t2);tl =a+b-t2;fl =Rx,p,t 1);wlule(l)if(fabs(tl-t2)eps)break;elset=(tl+t2)/2;b=t2; t2=tl;f2=fl; tl=a+b-t2;elsea=tl; tl=t2;fl =2t2=a+beta*(b-a);f2=f(x,p、t2);)t=(tl+t2)/2;return t;/*以卜是共粒梯度法程序源代碼*/void gtd() doub
4、le xN.gN.pNj=0.fD,modl=0,mod2=0aianda=0; mt i,k.n;pnntf(”請(qǐng)輸入函數(shù)的元數(shù)值n=);scanf(”d”,&n);pnntf(”n請(qǐng)輸入初始值:n”);foi(i=O;ieps)pO=-gO;pl=-gl;k=O;while(l)(t=hjfg(x,p);/*調(diào)用黃金分割法求t的值*/pnntfCnpl=%lf;p2=%lft=%lfp0,pl,t);xO=xO+t*pO; xl=xl+t*pl;g0=2*x0;gl=5O*xl;/*printfl,uxl=%lf,x2=%lf,gl=%lf,g2=%lf,x0,xl,g0,gl);*/mo
5、d2=sqit(pow(g 0 ,2 )+pow(gl,2); /* 求梯度的長(zhǎng)度 */if(mod2=eps) break;else(if(k+l=n)(g0=2*x0; gl=50*xl;pO=-gO;pl=-gl;k=O;else(nanda=pow(mod2,2)/pow(mod 1,2);piiiitfi(n,iinanda=%lmod=%lfnandajiiod2);pO=-gO+nanda*pO;Pl=-gl+nanda*pl;modl=mod2;k+;i)pnntfCn”);)printfCNn 最優(yōu)解為 x 1=%lf x2=%lf x0,x 1 );pnntf(n 最終的函
6、數(shù)值為lf,,f(x,g,t);main()gtd();運(yùn)行結(jié)果如下:請(qǐng)輸入函數(shù)的元數(shù)值n=2請(qǐng)輸入初始值:22xl=2.000000,x2=2.000000.pI=4000000,p2=100.000000a,b=-4.500000,1.500000pl=4000000,p2=100.000000.t=0.020030nanda=0.001474jnod=3.842730 xl=L919879,x2=-0.003022,pl=3.845655,p2=0.003665a,b=-4.500000,1.500000pl=3.845655,p2=0.003665,t=0.499240 xl=0.000026,x2=-0.001192,pl=0.000052,p2=0.059610a,b=-4.500000,1.500000pl=0.000052,p2=0.059610.t=0.020000nanda=0.000000jnod=0.000050 xl=0.000025,x2=-0.000000.pl=0.000050,p2=0.000001a,b=-4.500000,1.500000pl=0.000050,p2=0.000001.t=0.495505 xl=0.000000,x2=0.000000,pl=0.000000,p2=-0.000023a,b=-4.500000,1.
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 試析工序質(zhì)量控制在工民建施工中的意義
- 新北師大二年級(jí)數(shù)學(xué)教學(xué)研究討論
- 提升寫作水平的書目
- 二年級(jí)元角分?jǐn)?shù)學(xué)練習(xí)題解析
- 蘇教版六年級(jí)語(yǔ)文復(fù)習(xí)教學(xué)設(shè)計(jì)
- 蘇教版高中物理教材練習(xí)目錄
- 高效學(xué)習(xí)公頃蘇教版課件幫你搞定
- 蘇教版高中新語(yǔ)文教材解析
- 蘇教版二年級(jí)上冊(cè)科學(xué)教學(xué)計(jì)劃教學(xué)改革策略
- 數(shù)學(xué)王國(guó)蘇教版小升初復(fù)習(xí)之旅
- 廣東省深圳市2024年八年級(jí)下學(xué)期語(yǔ)文期中試卷(附答案)
- 《奇山異水、多彩民俗-西南旅游區(qū)導(dǎo)學(xué)案-中國(guó)旅游地理》
- 2024年畜禽屠宰企業(yè)獸醫(yī)衛(wèi)生檢驗(yàn)人員考試試題
- 三病反歧視制度
- (2024年)面神經(jīng)炎課件完整版
- 義務(wù)教育2o22版書法課程標(biāo)準(zhǔn)
- 醫(yī)院常見(jiàn)不安全因素及防范課件
- 抗腫瘤藥物培訓(xùn)課件
- 施工現(xiàn)場(chǎng)事故案例分析與總結(jié)
- 高壓電工復(fù)審培訓(xùn)課件
- 辦公軟件常用操作培訓(xùn)課件
評(píng)論
0/150
提交評(píng)論