复合形法C语言程序

1、现代设计理论与方法实验一、实验名称：复合形法参考程序上机实验二、实验目的：1:掌握复合形法优化问题的寻优策略2:能读懂程序并应用程序计算一些目标函数的最优解三、实验流程：1、题目：求如下约束优化问题的最优解F （ X ）= （xi -3）2 （X2 -4）2S.t.gi（X） =Xi 一0Ig2（X） = 为 一 0g3（X）=2.5-x X2-0.-IJ,g4（X） =5-x -X2 -0已知：N=2, Xi 亡 10,6曝2W0,8，取 k=4, E11032、编写黄金分割法的C语言程序代码；3、利用visualc程序运行C程序代码；尸 I *#includemath.hI#include

2、stdio.h,/ J 11 y#includestdlib.h#defineE10.001#defineep0.00001#definen2 #definek4 doubleaf;inti,j;doubleX0n,XXn,Xkn,FFk;doublean,bn;doublerm=2657863.0;doubleF(doubleCn)doubleF;F=pow(C0-3,2)+pow(C1-4,2);returnF;intcons(doubleDn)if(D0=0)&(D1=0)&(D0=6)&(D1=0)&( (5-D0-D1)=0)return1;elsereturn。;voidbou()

3、a0=0;b0=6;a1=0;b1=8;doubler()doubler1,r2,r3,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; .-I J Jreturnrr; jf-.voidproduce(doubleAn,doubleBn) I *Iintjj;doubleS;*1 丁一_ Jsl:for(i=0;in;i+)S=r()；XXi=Ai+S*(Bi-Ai);if(cons(XX)=0) gotosl；fo

4、r(i=0;in;i+)(X0i=XXi;)for(j=1;jk;j+)(for(i=0;in;i+)(S=r();Xji=Ai+S*(Bi-Ai);)for(j=1;jk;j+)(for(i=0;in;i+)I X I (X0i=0;for(jj=1;jjj+1;jj+)(X0i+=Xjji;)X0i=(1/j)*(X0i);)if(cons(X0)=0)gotosl;for(i=0;in;i+)XXi=X皿;while(cons(XX)=0)for(i=0;i.二/i r -Ji A 、j I 、1 1main()doubleEE,Xcn,Xhn,Xgn,X1n,Xrn,Xsn,w;int

5、l,lp,lp1;bou()； s111:produce(a,b);s222:for(j=0;jk;j+)(for(i=0;in;i+)(XXi=Xji;)FFj=F(XX);)for(l=0;lk-1;l+) I(for(lp=0;lpk-1;lp+) .- I J,(lp1=lp+1;if(FFlpFFlp1)(Iw=FFlp;FFlp=FFlp1;FFlp1=w;for(i=0;in;i+)(XXi=Xlpi;Xlpi=Xlp1i;Xlp1i=XXi;)for(i=0;in;i+)(Xhi=X0i;Xgi=X1i;X1i=Xk-1i;)for(i=0;in;i+)(Xsi=0;for(j

6、=0;jk;j+)(Xsi+=Xji;)Xsi=1/(k+0.0)*Xsi;)EE=0; J LI Afor(j=0;jk;j+)r(x z y - - -EE+=pow(FFj-F(Xs),2);)EE=pow(1/(k+0.0)*EE),0.5);if(EE=E1)(gotos333;(Xci=0;for(j=1;jk;j+)(Xci+=Xji;)Xci=1/(k-1.0)*Xci;)if(cons(Xc)=1)(af=1.3;ss:for(i=0;i=F(Xh)(if(af=ep)for(i=0;in;i+)Xhi=Xgi;af=1.3;gotoss;elseaf=1/2.0*af;go

7、toss;elsefor(i=0;in;i+)X0i=Xri;gotos222;I、 |/J y-_ J ielseaf=1/2.0*af;gotoss;elseif(X1iXci) ai=X1i;bi=Xci;elseai=Xci;bi=X1i;gotoslll;s333:printf(F(Xmin)=%fn,F(X1);for(i=0;i 0 g(X) = 2 - x2 - X1 0取：【k k 】k=4:1。I ，实验结果：2、minF(X)(1-x3)2 + 10=1OO(X2 - X1)2 + (1 - Xi)2 + 90(X4 -+(xl_ I)2 + (x4一 I)2 + 19.8(x2 - 1)(x4 - 1)s.t. -10 三 xi 1 10,(i = 1,2,3,4)取：x1 -10,10, x2-10,10, x3 -10,10x4 -10,10, k = 6, ； =10*实验结果：口：1: Docu*ent s and Se1: ngsActB.iis ra!: orDebugCpp 1. exeThe Xtl is 1 .皿8239 口The XI is 1007979 =The X2 Is 0.99S58G.The K3 is 0,?94114HPj*es：s： any key to continue223、 minF(X)=Xj + x2 -