复合形法C语言程序

1、现代设计理论与方法实验一、实验名称：复合形法参考程序上机实验二、实验目的： 1：掌握复合形法优化问题的寻优策略2：能读懂程序并应用程序计算一些目标函数的最优解三、实验流程：1、题目：求如下约束优化问题的最优解F ( X)=(x1 3)2 (x2 4)2S.t.g1(X) x1 0g2(X) x2 0g3(X) 2.5 x 1 x2 0g4(X) 5 x 1 x2 0已知： N=2, Xi 0,6 , x20,8 ,取 k=4, E1 1032 、 编写黄金分割法的 C 语言程序代码；3、 利用 visual c 程序运行 C 程序代码；#include "math.h"#i

2、nclude "stdio.h"#include "stdlib.h"#define E1 0.001#define ep 0.00001#define n 2#define k 4 double af;int i, j;double X0n,XXn,Xkn,FFk;double an,bn;double rm=2657863.0;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>

3、=0)&&(D0<=6)&&(D1<=8)&&(2.5-D0+D1)>=0)&&(5-D0-D1)>=0)return 1;elsereturn 0;void bou()a0=0;b0=6;a1=0;b1=8;double r()double r1,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=r

4、m/r1;return rr;void produce(double An,double Bn)int jj;double S;sl: for(i=0;i<n;i+)S=r();XXi=Ai+S*(Bi-Ai);if(cons(XX)=0)goto sl;for(i=0;i<n;i+)X0i=XXi;for(j=1;j<k;j+)for(i=0;i<n;i+)S=r();Xji=Ai+S*(Bi-Ai); for(j=1;j<k;j+)for(i=0;i<n;i+)X0i=0;for(jj=1;jj<j+1;jj+)X0i+=Xjji;X0i=(1/j

5、)*(X0i); if(cons(X0)=0)goto sl;for(i=0;i<n;i+)XXi=Xji; while(cons(XX)=0)for(i=0;i<n;i+)Xji=X0i+0.5*(Xji-X0i); XXi=Xji; main()double EE,Xcn,Xhn,Xgn,X1n,Xrn,Xsn,w; int l, lp, lp1; bou(); s111:produce(a,b); s222:for(j=0;j<k;j+) for(i=0;i<n;i+) XXi=Xji;FFj=F(XX); for(l=0;l<k-1;l+) for(lp=

6、0;lp<k-1;lp+) lp1=lp+1; if(FFlp<FFlp1) w=FFlp;FFlp=FFlp1;FFlp1=w; for(i=0;i<n;i+) XXi=Xlpi;Xlpi=Xlp1i;Xlp1i=XXi; for(i=0;i<n;i+)Xhi=X0i;Xgi=X1i;X1i=Xk-1i; for(i=0;i<n;i+)Xsi=0; for(j=0;j<k;j+) Xsi+=Xji;Xsi=1/(k+0.0)*Xsi;EE=0;for(j=0;j<k;j+)EE+=pow(FFj-F(Xs),2); EE=pow(1/(k+0.0)*

7、EE),0.5); if(EE<=E1)goto s333; for(i=0;i<n;i+)Xci=0;for(j=1;j<k;j+)Xci+=Xji;Xci=1/(k-1.0)*Xci; if(cons(Xc)=1)af=1.3;ss:for(i=0;i<n;i+)Xri=Xci+af*(Xci-Xhi);if(cons(Xr)=1) if(F(Xr)>=F(Xh) if(af<=ep)for(i=0;i<n;i+) Xhi=Xgi; af=1.3;goto ss; else else af=1/2.0*af;goto ss;for(i=0;i<

8、;n;i+)X0i=Xri;goto s222;elseaf=1/2.0*af;goto ss;elsefor(i=0;i<n;i+)if(X1i<Xci)ai=X1i;bi=Xci; elseai=Xci;bi=X1i;goto s111;s333:pri ntf("F(Xmi n)=%fn", F(X1);for(i=0;i< n ;i+)prin tf("n The X%d is %f.", i, X1i);prin tf("n 姓名：丁文浩 学号：201007001");四、实验结果:五、应用举例:复合形法程

9、序测试题1、mi nF (X) = (x 1- 2)2 + (x2- 1)2g(X) = x 2 - X12 > 0g(X) = 2 - x 2 - X1 > 0取：刘？【-5 6】X2?【-5 8k=4 ?= 10-3、) 2 X U (D +UE + LO 袋 XLLMX TOVTO、二 T ；二D co=O、7>7c' "C: Doca®«irt fl and SettiitcaXAclfciMxsr at ©r Debu Cppl _ «ze* F<Xnin；-0.60959i Im XH i s 1 .HMHZJ9 If to XI X3 丄 *.r>» X