无约束优化搜索方法鲍威尔法

1、无约束优化搜索方法上机实践报告1、写出所选择的无约束优化搜索算法的基本过程、原理（可附流程图说明）鲍威尔（Powell ）法是直接利用函数值来构造共轭方向的一种方法 在鲍威尔基本算法中，每一轮迭代都用连结始点和终点所产生出的搜 索方向去替换原向量组中的第一个向量，而不管它的“好坏”，这是产生向量组线性相关的原因所在。在改进的算法中首先判断原向量组是否需要替换。如果需要替换，还要进一步判断原向量组中哪个向量最坏，然后再用新产生的向量替换 这个最坏的向量，以保证逐次生成共轭方向。2、程序计算结果分析：中间各步骤的结果分析及与理论计算结果分析对比。 算例：ff=x0*x0+2*x1*x1-4*x0-

2、2*x0*x1;(1 )鲍威尔改进算法 C+ 程序清单：#include "stdio.h"#include "stdlib.h"#include "math.h"double objf(double x)double ff;ff=x0*x0+2*x1*x1-4*x0-2*x0*x1;return(ff);void jtf(double x0 ,double h0,double s ,intn,double a ,double b )int i;double *x3,h,f1,f2,f3;for (i=0;i<3;i+)xi=(

3、double *)malloc (n*sizeof(double);h=h0;for(i=0;i<n;i+)*(x0+i)=x0i;f1=objf(x0);for(i=0;i<n;i+)*(x1+i)=*(x0+i)+h*si;f2=objf(x1);if(f2>=f1)h= -h0;for (i=0;i<n;i+)*(x2+i)=*(x0+i);f3=f1;for(i=0;i<n;i+)*(x0+i)= *(x1+i);*(x1+i)= *(x2+i);f1=f2;f2=f3;for(;)h=2. *h;for(i=0;i<n;i+)*(x2+i)=* (

4、x1+i) +h*si;f3= objf(x2);if(f2<f3)break;else for(i=0;i<n;i+)*(x0+i)= *(x1+i);*(x1+i)= *(x2+i);f1=f2;f2=f3;if(h<0. )for(i=0;i<n;i+)ai=*(x2+i);bi=*(x0+i);elsefor(i=0;i<n;i+)ai=*(x0+i);bi=*(x2+i);for(i=0;i<3;i+)free(xi);double gold(double a,double b,double eps,int n,double xx) int i;d

5、ouble f1,f2,*x2,ff,q,w; for(i=0;i<2;i+)xi=(double*)malloc (n*sizeof(double); for(i=0;i<n;i+)*(x0+i)=ai+0.618*(bi-ai); *(x1+i)=ai+0.382*(bi-ai);f1=objf(x0);f2=objf(x1);doif(f1>f2)for(i=0;i<n;i+)bi=*(x0+i);*(x0+i)=*(x1+i);f1=f2;*(x1+i)=ai+0.382*(bi-ai);f2=objf(x1);elsefor(i=0;i<n;i+)ai=

6、*(x1+i);*(x1+i)=*(x0+i);f2=f1;for(i=0;i<n;i+)*(x0+i)=ai+0.618*(bi-ai); f1=objf(x0);q=0;for(i=0;i<n;i+)q=q+(bi-ai)*(bi-ai);w=sqrt(q);while(w>eps);for(i=0;i<n;i+)xxi=0.5*(ai+bi);ff=objf(xx);free(xi);return(ff);double oneoptim(double x0,double s,double h0,double n,double x)double *a,*b,ff;a

7、=(double *)malloc(n*sizeof(double);b=(double *)malloc(n*sizeof(double);jtf(x0,h0,s,n,a,b);ff=gold(a,b,epsg,n,x);free(a);free(b);return(ff);double powell(double p,double h0,double eps,double n,double x)int i,j,m;double *xx4,*ss,*s;double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d;ss=(double *)malloc(n*(n+1)*siz

8、eof(double);epsg,intepsg,ints=(double *)malloc(n*sizeof(double);for (j=0;j<=n;j+)*(ss+i*(n+1)+j)=0;*(ss+i*(n+1)+i)=1;for (i=0;i<4;i+)xxi=(double *)malloc(n*sizeof(double);for (i=0;i<n;i+)*(xx0+i)=pi;for(;)for (i=0;i<n;i+)*(xx1+i)=*(xx0+i);xi=*(xx1+i);f0=f1=objf(x);dlt=-1;for (j=0;j<n;

9、j+)for (i=0;i<n;i+)*(xx0+i)=xi;*(s+i)=*(ss+i*(n+1)+j);f=oneoptim(xx0,s,h0,epsg,n,x);df=f0-f;if(df>dlt)dlt=df;m=j;sdx=0.;for (i=0;i<n;i+)sdx=sdx+fabs(xi-(*(xx1+i);if(sdx<eps)free(ss);free(s);for (i=0;i<4;i+)free(xxi);return(f);for (i=0;i<n;i+)*(xx2+i)=xi;f2=f;for (i=0;i<n;i+)*(xx

10、3+i)=2.*(*(xx2+i)-(*(xx1+i);xi=*(xx3+i);fx=objf(x);f3=fx;q=(f1-2*f2+f3)*(f1-f2-dlt)*(f1-f2-dlt);d=0.5*dlt*(f1-f3)*(f1-f3);if(f3<f1)|(q<d)if(f2<=f3)for (i=0;i<n;i+)*(xx0+i)=*(xx2+i);elsefor (i=0;i<n;i+)*(xx0+i)=*(xx3+i);elsefor (i=0;i<n;i+)*(ss+(i+1)*(n+1)=xi-(*(xx1+i);*(s+i)=*(ss+(

11、i+1)*(n+1);f=oneoptim(xx0,s,h0,epsg,n,x);for(i=0;i<n;i+)*(xx0+i)=xi;for (j=m+1;j<=n;j+)*(ss+i*(n+1)+j-1)=*(ss+i*(n+1)+j);void main()double p=1,1;double ff,x2,x1,x2,f;ff=powell(p,0.3,0.001,0.0001,2,x); printf("shuchuzuiyoujie ：n"); x1=x1;x2=x2;f=ff;printf("x1=%f,x2=%f,f=%fn",x1,x2,f);getchar();(2)运行结果为：算例 1ff=x0*x0+x1*x1+x0+x0*x1输出结果：shuchuzuiyoujie u|x1=