优化设计鲍威尔方向法实验程序0.docVIP

/*鲍威尔共轭方向法C语言程序*/ /* a[0]--初始区间的下界值； ad——一维迭代精度及收敛精度； tt——一维搜索的初始迭代步长； xk[n]——迭代初始点； */ # include stdio.h # includestdlib.h # includemath.h # includemalloc.h # includeconio.h # define n 2 # define tt 0.005 # define ad 0.0000001 float ia; float fny(float *x) { float f; f=10*pow((x[0]+x[1]-5),2)+pow((x[0]-x[1]),2); /*目标函数*/ return(f); } float *iterate(float *x,float a, float *s) { float *x1; int i; x1=(float *)malloc(n*sizeof(float)); for( i=0;in;i++) x1[i]=x[i]+a*s[i]; return(x1); } float func(float *x,float a,float *s) { float *x1; float f; x1=iterate(x,a,s); f=fny(x1); return(f); } void finding(float a[3],float f[3], float *xk, float *s) { float t=tt; float a1,f1; int i; a[0]=0; f[0]=func(xk,a[0],s); for( i=0;;i++) { a[1]=a[0]+t; f[1]=func(xk,a[1],s); if(f[1]f[0]) break; if(fabs(f[1]-f[0])=ad) { 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]=func(xk,a[2],s); 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; } /* second insert*/ float lagrange(float *xk,float *ft, float *s) { float a[3],f[3]; float b,c,d,aa; int i; finding(a,f,xk,s); for( i=0;;i++) { if(ia==1) {aa=a[1]; *ft=f[1]; break; } d=(pow(a[0],2)-pow(a[2],2))*(a[0]-a[1])-(pow(a[0],2)-pow(a[1],2))*(a[0]-a[2]); if(fabs(d)==0) break; c=((f[0]-f[2])*(a[0]-a[1])-(f[0]-f[1])*(a[0]-a[2]))/d; if(fabs(c)==0) break; b=((f[0]-f[1])-c*(pow(a[0],2)-pow(a[1],2)))/(a[0]-a[1]); aa=-b/(2*c); *ft=func(xk,aa,s); if(fabs(aa-a[1])=ad) { if(*ftf[1]) aa=a[1]; break; } if(aaa[1]) { if(*ftf[1]) { a[2]=aa; f[2]=*ft; } else if(*ftf[1]) {a[0]=a[1]; a[1]=aa; f[0]=f[1]; f[1]=*ft; } else if(*ft==f[1]) { a[2]=aa; a[0]=a[1]; f[2]=*ft; f[0]=f[1]; a[1]=(a[0]+a[2])/2; f[1]=func(xk, a[1],s); } } else { if(*ftf[1]) { a[0]=aa; f[0]=*ft; } else if (*ftf[1]) { a[2]=a[1]; a[1]=aa; f[2]=f[1]; f[1]=*ft; } else if (*ft==f[1]) { a[0]=aa; a[2]=a[1]; f[0]=*ft; f[2]