现代设计一维优化方法——黄金分割法c程序源代码.docVIP

现代设计一维优化方法——黄金分割法c程序 该程序是现代设计优化方法，黄金分割法的c语言源代码 该程序无错误，已经调试，后面附有实例，仅供参考 希望对你有所帮助………… #include stdio.h #include math.h #include conio.h #define e 0.001 #define tt 0.01 double function(double x) { double y; y=8*pow(x,3)-2*pow(x,2)-7*x+3; return(y); } void searching(double a[3],double f[3]) {double t=tt,a1,f1,ia; int i; a[0]=0; f[0]=function(a[0]); for(i=0;;i++) {a[1]=a[0]+t; f[1]=function(a[1]); if(f[1]f[0]) break; if(fabs(f[1]-f[0])=e) { t=-t;a[0]=a[1];f[0]=f[1]; } else{ if(ia==1) return; t=t/2;ia=1; } } for(i=0;;i++) { a[2]=a[1]+t; f[2]=function(a[2]); if(f[2]f[1]) break; t=2*t; a[0]=a[1];f[0]=f[1]; a[1]=a[2];f[1]=f[2]; } if(a[0]a[2]) { a1=a[0];f1=f[0]; a[0]=a[2];f[0]=f[2]; a[2]=a1;f[2]=f1; } return; } double gold(double *ff) { double a1[3],f1[3],a[4],f[4]; double aa; int i; searching(a1,f1); a[0]=a1[0];f[0]=f1[0]; a[3]=a1[2];f[3]=f1[2]; a[1]=a[0]+0.382*(a[3]-a[0]); a[2]=a[0]+0.618*(a[3]-a[0]); f[1]=function(a[1]); f[2]=function(a[2]); for(i=0;;i++) { if(f[1]=f[2]) {a[0]=a[1];f[0]=f[1]; a[1]=a[2];f[1]=f[2]; a[2]=a[0]+0.618*(a[3]-a[0]); f[2]=function(a[2]); } else{ a[3]=a[2];f[3]=f[2]; a[2]=a[1];f[2]=f[1]; a[1]=a[0]+0.382*(a[3]-a[0]); f[1]=function(a[1]); } if(a[3]-a[0]e) { aa=(a[1]+a[2])/2;*ff=function(aa); break; } } return(aa); } void main() {double xx, ff; xx=gold(ff); printf(\nThe Optimal Design Result Is:\n); printf(\n\tx*=%f\n\tf*=%f,xx,ff); getch(); } 实例： f(x)=8*x^3-2*x^2-7*x+3 取：a[0]=0,e=0.001,tt=0.01 运行以上程序，最优解： x*=0.629867 f*=-0.203425