复合形法C语言程序.docx

复合形法C语言程序
作者:
日期:
现代设计理论与方法实验
一、实验名称复合形法参考程序上机实验
二、实验目的：1:掌握复合形法优化问题的寻优策略
2:能读懂程序并应用程序计算一些目标函数的最优
解
三、实验流程：
1、题目：求如下约束优化问题的最优解
F ( X ) = (% 3)2(X2 4)2
. gi(X) Xi 0
g2(X) X2 0
g3(X) X1 X2 0
g4(X) 5 X1 X2 0 3
已知:N=2, x1 0,6,x2 0,8，取 k=4, Ei 10
2、编写黄金分割法的C语言程序代码；
3、利用v 1 sual c程序运行C程序代码；
i n c lude " m ath. h"
#incl u de ""
#i n cl u d e ""
d e fine E 1 0. 00 1
#defi n e ep
#define n 2
#d e fin e k 4
double af;
int i , j ；
double X0[n],XX[n],X[k][n],FF [ k]; do ub le a[n ] ,b [ n];
double r m=2 6 5 7 ;
d oubl e F (doub 1 e C[n])
d o ubl e F;
F=pow(C[ 0 ] 3, 3,2)+pow( C[1] - 4 , 2);
?re t urn F ;
}
int cons(doub 1 e D [ n ])
{
?if((D⑼ >=0) &&(D [1]>=0 )&&(D[0] <=6 )&&( D[ 1 ]< = 8)&&((-D ⑼+D[1 ] )>=0)&& ((5-D[0] -D[1])>=0))
?? r e t urn 1;
e lse
return 0;
}
void bou ()
{
?a[0]=0;b[0] = 6 ;
?a[ 1 ] = 0;b [1] =8;
}
do u b l e r ()
{
? d ouble r 1 , r 2, r 3 ,rr;
r1=pow(2,35); r 2 = pow(2,36); r 3 = po w (2,3 7 );r m = 5* r m;
if(r m > =r3){rm= r m — r3 ;}
if(r m >=r2) {rm= r m — r2; }
i f(rm>= r 1 ) { r m = rm-rl ;}
?rr=r m /r1;
?r e t urn rr;
}
v o id p r od u ce (dou b le A[n],d o ub 1 e B[n])
{
int j j ; double S ;
si: f o r(i= 0 ;i< n ; i++)
{
S = r ();
? XX [i]=A[i]+S*(B[i]-A[ i ]);
}
i f ( con s (XX)==0)
{goto sl;}
f o r(i=0;i< n ;i++)
{
X ：0] [i]=XX[i];
?}
for( j =1 ； j v k;j++)
for(i=0;i<n; i ++)
? {
? S =r();
X [j] [i]= A[i]+S* (B [i] -A[ i ]);
}
}
for(j= 1 ； j< k ;j++)
{
?for (i= 0 ; i <n;i++)
?{
?X0[i]=0;
for(jj=1；jj<j+ 1 ；jj + +)
? {
?? X0[i] + =X[jj][i];
?}
?X0[i]=(1/j)* (X0[i]);
}
?i f (cons (X0)==0)
{
g oto sl;
}
?for(i= 0 ;i< n ;i+ + )
?{X X[i] =X[j][i];}
w hi 1 e(cons(XX = =0)
{
? for ( i=0; i < n ;i++)
?{
?X[j][i]=X0[ i ]+ 0 . 5* (X [j][i]-X0[i ]);
?? XX [ i]=X[ j] [ i ];
??}
?}
}
}
m a i n ()
{
?double EE,Xc [n],Xh[n],Xg[n] , X1[n] ,Xr [n],Xs[n],w;
i n t l, lp, lp1;
bou();
s111:prod u ce ( a , b);