1、現(xiàn)代設(shè)計(jì)理論與方法實(shí)驗(yàn)一、實(shí)驗(yàn)名稱：復(fù)合形法參考程序上機(jī)實(shí)驗(yàn)二、實(shí)驗(yàn)?zāi)康模?:掌握復(fù)合形法優(yōu)化問(wèn)題的尋優(yōu)策略2：能讀懂程序并應(yīng)用程序計(jì)算一些目標(biāo)函數(shù)的最優(yōu)解三 實(shí)驗(yàn)流程：1、題目：求如下約束優(yōu)化問(wèn)題的最優(yōu)解F ( X)= (Xi 3)2(X2 4)2gi(x) xi ogs(x)x2 og3(X) 2.5 xi X2Og4(X) 5 xi X2 0已知：N=2, xi 0,6 , X2 0,8 ,取 k=4, E1 103o2 編寫黃金分割法的c語(yǔ)言程序代碼；3、利用visual c程序運(yùn)行C程序代碼；include ”include ”include “”#define E1#define

2、ep#define n 2#define k 4 double af;int i, j；double X0n,XXn,Xkn,FFk;double an,bn;double rm=;double F(double Cn)(double F;F=pow(C0-3,2)+pow(C1 -4,2);return F;int cons(double Dn)(if(D0=0)&(D1 =0)&(D0=6)&(D1 =0)&(5-D0-D1 )=0) return1；elsereturn 0;void bou()(a0=0;b0=6;a1=0;b1=8;double r()(double r1,r2,r3

3、,rr; r1=pow(2,35);r2=pow(2,36);r3=pow(2,37);rm=5*rm;if(rm=r3)rm=rm-r3;if(rm=r2)rm=rm-r2;if(rm=r1 )rm=rm-r1;rr=rm/r1;return rr;void produce(double An,double Bn)(int jj;double S;si: for(i=0;in;i+)S=r();XXi=Ai+S*(Bi-Ai);)if(cons(XX)=0)goto si;for(i=0;in;i+)X0i=XXi;)for(j=1;jk;j+)(for(i=0;in;i+)S=r()；Xj

4、i=Ai+S*(Bi-Ai);)for(j=1;jk;j+)(for(i=0;in;i+)X0i=0;for(jj=1 ;jjj+1 ；jj+)(xoi+=xoji；XOi=(1/j)*(XOi);)if(cons(X0)=0)goto si;)for(i=0;in;i+)XXi=Xji; while(cons(XX)=0)for(i=0;in;i+)xoi=xoi+*(xoi-xoi)；xxi=xoi；main()(double EE,Xcn,Xhn,Xgn,X1n,Xrn,Xsn,w; int I, Ip, Ip1: bou(); s111:produce(a,b); s222:for(j

5、=0;jk;j+)for(i=0;in;i+) xxi=X0i；FFj=F(XX);)for(l=0;lk-1 ;l+)for(lp=0;lpk-1 ;lp+)(lp1=lp+1;if(FFlpFFlp1)(w=FFlp;FFlp=FFlp1 ;FFlp1 =w; for(i=0;in;i+) XX卜 Xlpi;Xlp卜 Xlp1i;Xlp1卜 XXi;)for(i=0;ivn;i+)(Xhi=X0i;Xgi=X1i;X1 i=Xk-1i;for(i=0;ivn;i+)(Xsi=0;for(j=0;jk;j+)Xsi+=X0i；)Xsi=1/(k+*Xsi;EE=0;for(j=0;jvk;j

6、+)(EE+=pow(FF0-F(Xs),2);EE=pow(1/(k+*EE),;if(EE=E1)(goto s333;for(i=0;ivn;i+)(Xci=0; for(j=1 ;jk;j+)(Xci+=X0i；)Xci=1/*Xci;if(cons(Xc)=1)af=；ss:for(i=0;i=F(Xh)(if(af=ep)(for(i=0;in;i+)Xhi=Xgi;) af=;goto ss;else af=1/*af;goto ss;else (for(i=0;ivn;i+) XOi=Xri;goto s222;)elseaf=1/*af;goto ss;)elsefor(i=

7、0;in;i+)(if(X1iXci) ai=X1i;bi=Xci;elseai=Xci;bi=X1i;goto s111;s333:pri ntf(HF(Xmi n)=%fnn, F(X1);for(i=0;i n;i+)prin tf(un The X%d is %f.”, i, X1i);prin tf(Hn 姓名：丁文浩 學(xué)號(hào)：1”);四、實(shí)驗(yàn)結(jié)果:c， *C:Ducu*ent s and Set tingsAd*inisttaloiDebugCpp1eiw* F 0g(X) = 2 - x2 - Xi 0取：Xi -5 6x2 -5 8k=4 ?= 10 3實(shí)驗(yàn)結(jié)果:岡 4 ( ris10iw11. wuw/y 0.99558.A. 9941 14. Prt? x xkr?y 11 cout iiuie*F-0.008959c *C : DocTcn-t s and Sct-tAdTnxwTrat orXDebngkCppi . ctcH 1 2 3 K X X M1009 4 x g 60x ontfrquoo oKuv 6E9066 ST PX 叫丄9乙6666