共軛梯度法,機械優(yōu)化設計

1、實驗報告實驗課程名稱機械優(yōu)化設計實驗項目名稱共軻梯度法年級專業(yè)學生姓名學號實驗時間：2012年11月2日學生所在學院：專業(yè)：班級:姓名學號實驗組實驗時間指導教師成績實驗項目名稱共軻梯度法求函數(shù)極小值1、實驗目的及要求：實驗目的：掌握機械優(yōu)化設計方法并能夠理論聯(lián)系實際地加以應用，任務是將課程所學的知識應用于實踐，通過實際編寫調(diào)試及運行程序加深理論知識的掌握并提高解決優(yōu)化問題的能力。根據(jù)實驗指導書的要求應能夠獨立的編寫優(yōu)化程序并在計算機上運行，學會判斷結(jié)果及程序的正確性，學會建立機械優(yōu)化設計的數(shù)學模型，合理選用優(yōu)化方法，獨立的解決機械優(yōu)化設計的實際問題。實驗(或算法)原理：共軻梯度法是共軻方向法中

2、的一種，該方法中每一個共軻向量都是依賴與迭代點處的負梯度而構(gòu)造出來。它通過梯度來尋找極小點。先通過一維搜索確定搜索區(qū)間，然后再通過共軻梯度法運用c語言編程求解。實驗硬件及軟件平臺：VS2010實驗步驟：1 .確定所需求解的函數(shù)y=pow(x0+t*p0,2)+25*pow(x1+t*p1,2)2 .確定搜索區(qū)間3 .畫出程序框圖4 .用c語言在vs2010上寫出源代碼5 .運行程序6 .檢驗試驗結(jié)果，分析結(jié)果實驗內(nèi)容（包括實驗具體內(nèi)容、算法分析、源代碼等等）本實驗通過c語百編程,運用共軻梯度法求解函數(shù)y極小值；程序框圖源程序#include#include#defineN10#defineep

3、spow(10,-6)doublef(doublex口，doublep口，doublet)doubles;s=pow(x0+t*p0,2)+25*pow(x1+t*p1,2);returns;voidsb(double*a,double*b,doublex口，doublep口)doublet0,t1,t,h,alpha,f0,f1;intk=0;t0=2.5;/*初始值*/h=1;/*初始步長*/alpha=2;/*加步系數(shù)*/f0=f(x,p,t0);t1=t0+h;f1=f(x,p,t1);while(1)if(f1f0)h=alpha*h;t=t0;t0=t1;f0=f1;k+;else

4、if(k=0)h=-h;t=t1;else*a=tt1?t:t1;break;t1=t0+h;f1=f(x,p,t1);doublehjfg(doublex口，doublep口)doublebeta,t1,t2,t;doublef1,f2;doublea=0,b=0;double*c,*d;c=&a,d=&b;,x0,x1,p0,p1);pause);sb(c,d,x,p);printf(nx1=%lf,x2=%lf,p1=%lf,p2=%lfprintf(na,b=%lf,%lf,a,b);system(beta=(sqrt(5)-1.0)/2;t2=a+beta*(b-a);f2=f(x,

5、p,t2);t1=a+b-t2;f1=f(x,p,t1);while(1)if(fabs(t1-t2)eps)break;elseif(f1f2)t=(t1+t2)/2;b=t2;t2=t1;f2=f1;t1=a+b-t2;f1=f(x,p,t1);elsea=t1;t1=t2;f1=f2;t2=a+beta*(b-a);f2=f(x,p,t2);t=(t1+t2)/2;returnt;voidgtd()doublexN,gN,pN,t=0,f0,mod1=0,mod2=0,nanda=0;inti,k,n;printf(請輸入函數(shù)的元數(shù)值n=);scanf(%d,&n);printf(n請輸

6、入初始值n);for(i=0;ieps)p0=-g0;p1=-g1;k=0;while(1)t=hjfg(x,p);printf(np1=%lf,p2=%lf,t=%lf,p0,p1,t);x0=x0+t*p0;x1=x1+t*p1;g0=2*x0;g1=50*x1;/*printf(nx1=%lf,x2=%lf,g1=%lf,g2=%lf,x0,x1,g0,g1);*/mod2=sqrt(pow(g0,2)+pow(g1,2);if(mod2=eps)break;elseif(k+1=n)g0=2*x0;g1=50*x1;p0=-g0;p1=-g1;k=0;elsenanda=pow(mod2,2)/pow(mod1,2);printf(nnanda=%lf,mod=%lf,nanda,mod2);P0=-g0+nanda*p0;P1=-g1+nanda*p1;mod1=mod2;k+;printf(n);printf(n最優(yōu)解a為ax1=%lf,x2=%lf,x0,x1);printf(n最終