复合形法C语言程序.docx
- 文档编号:3051268
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:10
- 大小:104.54KB
复合形法C语言程序.docx
《复合形法C语言程序.docx》由会员分享,可在线阅读,更多相关《复合形法C语言程序.docx(10页珍藏版)》请在冰豆网上搜索。
复合形法C语言程序
现代设计理论与方法实验
一、实验名称:
复合形法参考程序上机实验
二、实验目的:
1:
掌握复合形法优化问题的寻优策略
2:
能读懂程序并应用程序计算一些目标函数的最优解
三、实验流程:
1、题目:
求如下约束优化问题的最优解
(
)=
S.t.
已知:
N=2,
取k=4,
。
2、编写黄金分割法的C语言程序代码;
3、利用visualc程序运行C程序代码;
#include"math.h"
#include"stdio.h"
#include"stdlib.h"
#defineE10.001
#defineep0.00001
#definen2
#definek4
doubleaf;
inti,j;
doubleX0[n],XX[n],X[k][n],FF[k];
doublea[n],b[n];
doublerm=2657863.0;
doubleF(doubleC[n])
{
doubleF;
F=pow(C[0]-3,2)+pow(C[1]-4,2);
returnF;
}
intcons(doubleD[n])
{
if((D[0]>=0)&&(D[1]>=0)&&(D[0]<=6)&&(D[1]<=8)&&((2.5-D[0]+D[1])>=0)&&((5-D[0]-D[1])>=0))
return1;
else
return0;
}
voidbou()
{
a[0]=0;b[0]=6;
a[1]=0;b[1]=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;
returnrr;
}
voidproduce(doubleA[n],doubleB[n])
{
intjj;doubleS;
sl:
for(i=0;i { S=r(); XX[i]=A[i]+S*(B[i]-A[i]); } if(cons(XX)==0) {gotosl;} for(i=0;i { X[0][i]=XX[i]; } for(j=1;j { for(i=0;i { S=r(); X[j][i]=A[i]+S*(B[i]-A[i]); } } for(j=1;j { for(i=0;i { X0[i]=0; for(jj=1;jj { X0[i]+=X[jj][i]; } X0[i]=(1/j)*(X0[i]); } if(cons(X0)==0) { gotosl; } for(i=0;i {XX[i]=X[j][i];} while(cons(XX)==0) { for(i=0;i { X[j][i]=X0[i]+0.5*(X[j][i]-X0[i]); XX[i]=X[j][i]; } } } } main() { doubleEE,Xc[n],Xh[n],Xg[n],X1[n],Xr[n],Xs[n],w; intl,lp,lp1; bou(); s111: produce(a,b); s222: for(j=0;j { for(i=0;i { XX[i]=X[j][i]; } FF[j]=F(XX); } for(l=0;l { for(lp=0;lp { lp1=lp+1; if(FF[lp] { w=FF[lp];FF[lp]=FF[lp1];FF[lp1]=w; for(i=0;i { XX[i]=X[lp][i];X[lp][i]=X[lp1][i];X[lp1][i]=XX[i]; } } } } for(i=0;i { Xh[i]=X[0][i];Xg[i]=X[1][i];X1[i]=X[k-1][i]; } for(i=0;i { Xs[i]=0; for(j=0;j { Xs[i]+=X[j][i]; } Xs[i]=1/(k+0.0)*Xs[i]; } EE=0; for(j=0;j { EE+=pow((FF[j]-F(Xs)),2); } EE=pow((1/(k+0.0)*EE),0.5); if(EE<=E1) { gotos333; } for(i=0;i { Xc[i]=0; for(j=1;j { Xc[i]+=X[j][i]; } Xc[i]=1/(k-1.0)*Xc[i]; } if(cons(Xc)==1) { af=1.3; ss: for(i=0;i { Xr[i]=Xc[i]+af*(Xc[i]-Xh[i]); } if(cons(Xr)==1) { if(F(Xr)>=F(Xh)) { if(af<=ep) { for(i=0;i { Xh[i]=Xg[i]; } af=1.3;gotoss; } else {af=1/2.0*af;gotoss;} } else { for(i=0;i { X[0][i]=Xr[i]; } gotos222; } } else {af=1/2.0*af;gotoss;} } else { for(i=0;i { if(X1[i] {a[i]=X1[i];b[i]=Xc[i];} else {a[i]=Xc[i];b[i]=X1[i];} } gotos111; } s333: printf("F(Xmin)=%f\n",F(X1)); for(i=0;i { printf("\nTheX%dis%f.",i,X1[i]); } printf("\n姓名: 丁文浩学号: 201007001"); } 四、实验结果: 五、应用举例: 复合形法程序测试题 1、 取: k=4 实验结果: 2、 实验结果: 3、 实验结果:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 复合 语言 程序