機(jī)械優(yōu)化設(shè)計(jì)復(fù)合形算法程序設(shè)計(jì)報(bào)告

第第頁機(jī)械優(yōu)化設(shè)計(jì)復(fù)合形算法程序設(shè)計(jì)報(bào)告一、設(shè)計(jì)要求：

基于復(fù)合形算法思想，運(yùn)用復(fù)合形法編寫C語言程序，完成一次迭代。

已知條件：

1、目標(biāo)函數(shù)：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]

反射點(diǎn)迭代求解：完成一次迭代并輸出

二、方法原理

(1)

復(fù)合形算法是求解約束優(yōu)化問題的一種重要的徑直解法。它的基本思路是在可行域內(nèi)構(gòu)造一個(gè)具有K個(gè)頂點(diǎn)的初始復(fù)合形。對(duì)該復(fù)合形各頂點(diǎn)的目標(biāo)函數(shù)值進(jìn)行比較，找到目標(biāo)函數(shù)值最大的頂點(diǎn)，然后按肯定的法那么求出目標(biāo)函數(shù)值有所下降的可行的新點(diǎn)，并用此點(diǎn)代替最壞點(diǎn)，構(gòu)成新的復(fù)合形，復(fù)合形的外形每轉(zhuǎn)變一次，就向最優(yōu)點(diǎn)移動(dòng)一步，直至迫近最優(yōu)點(diǎn)。

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(排序后,按函數(shù)值的優(yōu)劣從小到大：);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,計(jì)算完成,i,m);

}

四、試驗(yàn)結(jié)果〔要求附上程序運(yùn)行結(jié)果截圖〕

一、設(shè)計(jì)要求：

基于復(fù)合形算法思想，運(yùn)用復(fù)合形法編寫C語言程序，完成一次迭代。

已知條件：

1、目標(biāo)函數(shù)：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]

反射點(diǎn)迭代求解：完成一次迭代并輸出

二、方法原理

(1)

復(fù)合形算法是求解約束優(yōu)化問題的一種重要的徑直解法。它的基本思路是在可行域內(nèi)構(gòu)造一個(gè)具有K個(gè)頂點(diǎn)的初始復(fù)合形。對(duì)該復(fù)合形各頂點(diǎn)的目標(biāo)函數(shù)值進(jìn)行比較，找到目標(biāo)函數(shù)值最大的頂點(diǎn)，然后按肯定的法