机械优化设计复合形算法程序设计报告

第第页机械优化设计复合形算法程序设计报告一、设计要求：

基于复合形算法思想，运用复合形法编写C语言程序，完成一次迭代。

已知条件：

1、目标函数：y=25/*1*2

2、初始值：

*1=[2.5,0.65],*2=[2.4,0.6],*3=[2.6,0.7],*4=[2.7,0.65]

反射点迭代求解：完成一次迭代并输出

二、方法原理

(1)

复合形算法是求解约束优化问题的一种重要的径直解法。它的基本思路是在可行域内构造一个具有K个顶点的初始复合形。对该复合形各顶点的目标函数值进行比较，找到目标函数值最大的顶点，然后按肯定的法那么求出目标函数值有所下降的可行的新点，并用此点代替最坏点，构成新的复合形，复合形的外形每转变一次，就向最优点移动一步，直至迫近最优点。

3

三、程序清单：

#includestdaf*.h

#includestdio.h

voidmain()

{

double*[4][2],y[4],*c[2]={0,0},m,n=1.3,g[2],k;inti,j,r;

for(i=0;i4;i++)

{

printf(input*%d:\n,i+1);

scanf(%lf,%lf,*[i][0],*[i][1]);

y[i]=25/(*[i][0]**[i][1]**[i][1]**[i][1]);printf(y%d=%3.2f\n,i+1,y[i]);

}

for(i=0;i3;i++)

{for(j=i+1;j4;j++)

if(y[i]y[j])

{m=y[i],k=y[j],y[j]=m,y[i]=k;

for(r=0;r2;r++)

{m=*[i][r],k=*[j][r],*[j][r]=m,*[i][r]=k;}}

}

printf(排序后,按函数值的优劣从小到大：);for(i=0;i4;i++)

printf(\ny%d=%3.2f

*%d=[%3.2f,%3.2f],i+1,y[i],i+1,*[i][0],*[i][1]);for(j=0;j2;j++)

for(i=0;i3;i++)

*c[j]=*c[j]+*[i][j]/3;

printf(\n中心*c=

[%3.2f,%3.2f],*c[0],*c[1]);

i=1;

for(j=0;j2;j++)

g[j]=*c[j]+n*(*c[j]-*[3][j]);

while(g[0]*g[1]*g[1]*g[1]0.6||g[0]*g[1]2.5){n=n*0.7;

printf(\n第%d次迭代后:*r=[%3.2f,%3.2f],i,g[0],g[1]);

m=25/(g[0]*g[1]*g[1]*g[1]);i++;

for(j=0;j2;j++)

g[j]=*c[j]+n*(*c[j]-*[3][j]);}m=25/(g[0]*g[1]*g[1]*g[1]);

printf(第%d次迭代Fr=%3.2fFh,计算完成,i,m);

}

四、试验结果〔要求附上程序运行结果截图〕

一、设计要求：

基于复合形算法思想，运用复合形法编写C语言程序，完成一次迭代。

已知条件：

1、目标函数：y=25/*1*2

2、初始值：

*1=[2.5,0.65],*2=[2.4,0.6],*3=[2.6,0.7],*4=[2.7,0.65]

反射点迭代求解：完成一次迭代并输出

二、方法原理

(1)

复合形算法是求解约束优化问题的一种重要的径直解法。它的基本思路是在可行域内构造一个具有K个顶点的初始复合形。对该复合形各顶点的目标函数值进行比较，找到目标函数值最大的顶点，然后按肯定的法