Powell法求解极小化问题.docVIP

数学与计算科学学院 实 验 报 告 实验项目名称 Powell法求解极小化问题 所属课程名称 最优化方法 实 验 类 型 算法编程 实 验 日 期 2015.11.20 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 1、掌握Powell方法求解； 2、了解Powell的二次终止性； 【实验原理】 1、基本思想：Powell方法把整个计算过程分成若干个阶段，每一阶段由n+1次一维搜索组成，在算法的每一阶段中，先依次沿着已知的n个方向搜索得一个最好点，然后沿着本阶段的初点与该最好点连线方向进行搜索，求得这一阶段的最好点，再用最后搜索的方向取代前n个方向之一开始下一阶段。 2、二次终止性：设，A为n阶对称正定矩阵。任意给定方向从出发，沿方向d搜索，得极小点；从出发，沿方向d搜索，得极小点；则与d关于A共轭。 当极小化正定二次函数时，如果每轮迭代中前n个方向均线性无关，那么Powell方法至多经过n轮迭代达到极小点。 3、计算步骤： Step1：给定初始点，n个线性无关方向：， 误差； Step2：让，从出发，依次沿方向进行搜索，得点，再从出发，沿方向作一维搜索得到。 Step3：若，则停止计算，得，否则令 Step4：让，返回步2 【实验环境】 Windows7 Visual C++ 6.0 二、实验内容： 【实验方案】 1、本次试验的题目为： 用Powell方法求解下列问题： 2、根据Powell方法的求解步骤，结合本题，运用C语言编写程序，在计算机上实现。 【实验过程】（实验步骤、记录、数据、分析） 1、由题目可得，其中海瑟矩阵是正定矩阵，并且是一组G共轭的非零方向，所以本题具有二次终止性。 2、通过定义函数，将一阶导输入，用一维数组存储海瑟矩阵： float g1(float x[]){ float g=4*x[0]+2*x[1]-2; return g; } float g2(float x[]){ float g=2*x[1]+2*x[0]; return g; } 3、根据一维线性搜索，编写函数，求步长a： float qiubuchang(float x[],float g[],float d[],float G[]){ g[0]=g1(x); g[1]=g2(x); float a,b,c; b=g[0]*d[0]+g[1]*d[1]; c=d[0]*(d[0]*G[0]+d[1]*G[2])+d[1]*(d[0]*G[1]+d[1]*G[3]); a=-b/c; return a; } 4、在主函数中调用以上函数，根据Powell法的步骤，经过两轮迭代，求出结果。 【实验结论（结果）小结结论#includemath.h #includestdio.h float g1(float x[]){ float g=4*x[0]+2*x[1]-2; return g; } float g2(float x[]){ float g=2*x[1]+2*x[0]; return g; } float qiubuchang(float x[],float g[],float d[],float G[]){ g[0]=g1(x); g[1]=g2(x); float a,b,c; b=g[0]*d[0]+g[1]*d[1]; c=d[0]*(d[0]*G[0]+d[1]*G[2])+d[1]*(d[0]*G[1]+d[1]*G[3]); a=-b/c; return a; } void main(){ int i=1,j=1,k=1; float a; float x10[2],x11[2],x12[2],x13[2],x20[2],g[2],d11[2],d12[2],d13[2],d21[2],d22[2],G[4]={4,2,2,2}; printf(输入初始值：); scanf(%f,%f,x10[0],x10[1]); printf(输入搜索方向1：); scanf(%f,%f,d11[0],d11[1]); printf(输入搜索方向2：); scanf(%f,%f,d12[0],d12[1]); printf(第%d轮迭代：\n,i++); a=qiubuchang(x10,g,d11,G); printf(a=%f\n,a); x11[0]=x10[0]+a*d11[0