变尺度法课程设计.docVIP

变尺度法课程设计 班级：数学与应用数学07-2班 姓名：李彬彬 学号合肥工业大学 2011年4月3日 一：设计题目： 变尺度法求函数的极值计算机实现及其比较研究。 二：软件环境： Turbo C,Turbo C++,VC,VC++,Win-TC等能运行C语言的软件。 三：算法设计思路： 1：变尺度法的由来：尺度：指空间中对于点与点之间距离的一种度量，而对于正定矩阵A，若满足||x||j=?，称其为尺度矩阵，若f()为一般非线性函数，设其Hessen矩阵为正定阵，则拟牛顿条件为：（1）：（2）：；如果若以它为尺度矩阵，取最速下降方向为 , 在保持阻尼Newton法全局，快速收敛的同时，为了避免求逆矩阵，而通过边迭代边构造另外一个矩阵去近似，即令，由于尺度矩阵随迭代次数的变化而变化，故称此类算法为变尺度法 2：DFP法： a：基本思想：取初始矩阵（一般取为单位矩阵），用递推关系，满足拟牛顿条件确定，即第I次修正函数矩阵 b：DFP的一些性质：对于函数：具有二次终止性质，即，为Hessen逆阵；当A=E时产生共轭方向和共轭梯度。对与一般的函数：校正保持正定性，迭代下降的性质成立；具有超线性的收敛速度；且如果是凸函数，该方法具有总体收敛性。 四：程序流程图 五：演示运行结果 （1）当输入初值为5，7，精度为0时； （2）当输入的初值为2，12，精度为0时： 其他演示： 六：设计总结和体会。 总结：本程序的缺点，面向用户的功能比较差，无法让用户通过可执行文件直接输入要求的函数，而必须改动源程序，我很想弥补这个缺点，但由于水平和时间的有限，无法完成这部分工作，但我想这并不影响课程设计，因为这个课程设计的目的是实现常用的变尺度算法，而不是追求程序的完美性。 体会：通过这次《最优化理论》关于变尺度法的课程设计，强化了自己对C语言的掌握，同时更加深入的理解了变尺度法的算法原理，与计算过程。 附：算法源代码： 用变尺度法求函数+的最小值 #includestdio.h void dpfmethod(float a,float b,float m) { static float H[2][2]={1,0,0,1}; float x[2]; x[0]=a,x[1]=b; float c; c=4*(x[0]-5)*(x[0]-5)+(x[1]-6)*(x[1]-6); float df[2]; df[0]=8*(x[0]-5); df[1]=2*(x[1]-6); float hdf[2]; hdf[0]=(H[0][0]*df[0]+H[0][1]*df[1]); hdf[1]=(H[1][0]*df[0]+H[1][1]*df[1]); float namu=(-8*x[0]*hdf[0]+40*hdf[0]-2*hdf[1]*x[1]+12*hdf[1])/(8*hdf[0]*hdf[0]+2*hdf[1]*hdf[1]); float G[2]; G[0]=(-8*namu*hdf[0]); G[1]=(-2*namu*hdf[1]); float k=m; if(df[0]*df[0]+df[1]*df[1]=k){printf(最小值点是x1=%f，x2=%f\n,x[0],x[1]); printf(此时原函数的值为%f\n,c); getchar();} else { x[0]+=namu*hdf[0]; x[1]+=namu*hdf[1]; H[0][0]+=namu*hdf[0]*namu*hdf[0]/(G[0]*namu*hdf[0]+G[1]*namu*hdf[1])-((H[0][0]*G[0]+H[0][1]*G[1])*G[0]*H[0][0]+(H[0][0]*G[0]+H[0][1]*G[1])*G[1]*H[1][0])/(G[0]*(G[0]*H[0][0]+G[1]*H[1][0])+G[1]*(G[0]*H[0][1]+G[1]*H[1][1])); H[0][1]+=namu*hdf[0]*namu*hdf[1]/(G[0]*namu*hdf[0]+G[1]*namu*hdf[1])-((H[0][0]*G[0]+H[0][1]*G[1])*G[0]*H[0][1]+(H[0][0]*G[0]+H[0][1]*G[1])*G[1]*H[1][1])/(G[0]*(G[0]*H[0][0]+G[1]*H[1][0])+G[1]*(G[0]*H[0][1]+