收藏 分享(赏)

无约束优化搜索方法鲍威尔法(齐童).doc

上传人:weiwoduzun 文档编号:5632644 上传时间:2019-03-10 格式:DOC 页数:8 大小:125KB
下载 相关 举报
无约束优化搜索方法鲍威尔法(齐童).doc_第1页
第1页 / 共8页
无约束优化搜索方法鲍威尔法(齐童).doc_第2页
第2页 / 共8页
无约束优化搜索方法鲍威尔法(齐童).doc_第3页
第3页 / 共8页
无约束优化搜索方法鲍威尔法(齐童).doc_第4页
第4页 / 共8页
无约束优化搜索方法鲍威尔法(齐童).doc_第5页
第5页 / 共8页
点击查看更多>>
资源描述

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 ,int n,double a ,double b )int i;double *x3,h,f1,f2,f3;for (i=0;i=f1)h= -h0;for (i=0;if2)for(i=0;ieps);for(i=

3、0;idlt)dlt=df;m=j;sdx=0.;for (i=0;in;i+)sdx=sdx+fabs(xi-(*(xx1+i);if(sdxeps)free(ss);free(s);for (i=0;i4;i+)free(xxi);return(f);for (i=0;in;i+)*(xx2+i)=xi;f2=f;for (i=0;in;i+)*(xx3+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)

4、;if(f3f1)|(qd)if(f2=f3)for (i=0;in;i+)*(xx0+i)=*(xx2+i);elsefor (i=0;in;i+)*(xx0+i)=*(xx3+i);elsefor (i=0;in;i+)*(ss+(i+1)*(n+1)=xi-(*(xx1+i);*(s+i)=*(ss+(i+1)*(n+1);f=oneoptim(xx0,s,h0,epsg,n,x);for(i=0;in;i+)*(xx0+i)=xi;for (j=m+1;j=n;j+)for (i=0;in;i+)*(ss+i*(n+1)+j-1)=*(ss+i*(n+1)+j);void main()

5、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输出结果:shuchuzuiyoujiex1=0.333246,x2=0.333246,f=-0.333333算例 2ff=3*x0*x0+6*x1*x1+2*x0+(x0+5)*x1输出结果:shuchuzuiyoujiex1=-0.394357,x2=-0.394357,f=-1.253521感想本次上机让我感触颇多,让我学会了使用鲍威尔法进行无约束优化问题的分析,感谢老师对我的指导。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 企业管理 > 管理学资料

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报