进退法和黄金分割法.docVIP

进退法函数 1）程序头部 #includestdio.h #define N 2 /*定义N为常数2*/ double x[N]; double a,b; /*定义x[N]、a,b为全程变量*/ 2）目标函数 （例3-1） double fx(double x[]) { double f; f=x[0]*x[0]+x[1]*x[1]-8*x[0]-12*x[1]+52; return(f); } 3）调试主函数 （调试进退法专用） void main() { void jtf(double x0[],double s[],double t0); double t0=1; double x0[N]={2,2}; double s[N]={0.707,0.707}; jtf(x0,s,t0); printf(%f,%f\n,a,b); } void jtf(double x0[],double s[],double t0) /*jtf为空函数，结果a、b由全程变量返回*/ { doublet,a1,a2,f1,f2; double fx(double x[]); int i; t=t0; a1=0; a2=t; f1=fx(x0); /*调用目标函数*/ for(i=0;i=N-1;i++) x[i]=x0[i]+a2*s[i]; f2=fx(x); /*调用目标函数*/ if (f2f1) /*选择结构*/ { /*条件成立语句部分—左分支*/ do { /*条件循环语句部分*/ t=2*t; a2=a2+t; f1=f2; for(i=0;i=N-1;i++) x[i]=x0[i]+a2*s[i]; f2=fx(x); if(f1f2) a1=a2-t; } while (f1f2); /*条件成立则重复执行循环体，与直到F2F1相反*/ } else /*条件不成立部分—右分支*/ t = -t do {/*右分支条件循环语句部分*/ } while (f2f1); /*条件成立则重复执行循环体，与直到F1F2相反*/ } a=a1; b=a2; } /*函数定义结束*/ 黄金分割法 1）程序头部 #includestdio.h #includemath.h /*数学函数库*/ #define N 2 double x[N]; double a,b,fmin,alpha; /*最优值、最优步长*/ double epslx=0.0001; /*收敛精度*/ double epslf=0.0001; 2）目标函数 （例3-1） double fx(double x[]) { double f; f=x[0]*x[0]+x[1]*x[1]-8*x[0]-12*x[1]+52; return(f); } 3）主函数（修改） void main() { void jtf(double x0[],double s[],double t0); void hjf(double x0[],double s[]); double t0=1; double x0[N]={2,2}; double s[N]={0.707,0.707}; jtf(x0,s,t0); printf(%f,%f\n,a,b); hjf(x0,s); printf(%f,%f\n,x[0],x[1]); printf(%f,%f\n,alpha,fmin); } void hjf(double x0[],double s[]) { double x1[N],x2[N]; double f1,f2,a1,a2,ef; double Q=0.618; int j=0,i; do { ...... do { if(f1f2) ...... 左分支 else ...... 右分支 j=j+1; if(j50) break; } while(fabs((f2-f1)/f1)epslf); if(fabs((a2-a1)/a1)epslx) ...... } while(fabs((a2-a1)/a1)epslx); if(f1f2) { alpha=a1; for(i=0;i=N-1;i++) x[i]=x0[i]+a1*s[i]; fmin=f1; } else { alpha=a2; for(i=0;i=N-1;i++) x[i]=x0[i]+a2*s[i]; fmin=f2; } } /*函数结束*/