《数学规划 (最速下降法,c语言编程).doc》由会员分享，可在线阅读，更多相关《数学规划 (最速下降法,c语言编程).doc（9页珍藏版）》请在金锄头文库上搜索。
1、数 学 规 划 课 程 设 计题目：用最速下降法求解无约束非线性规划问题姓名： 学号： 成绩： 2011年6月用最速下降法求解无约束非线性规划问题摘要：无约束非线性规划问题是一类重要的数学规划问题。文主要研究了用最速下降法也就是梯度法对无约束非线性规划问题进行求解。对于一个无约束非线性规划利用最速下降法求解，首先需要确定其优化方向，此优化方向应该选择为f在当前点处的负梯度方向，利用一维搜索法找出沿此方向上的最小值及其对应点，此后将该点作为新的出发点重复上述过程，直到达到允许的误差为止。本文最后利用c+语言编程得到满足允许误差内的最优解。本文主要对一个无约束非线性规划问题的实例，首先利用上述迭代
2、的方法，计算出各迭代点的函数值，梯度及其模。然后应用c+语言编程，得到精确地最优解，需迭代六次才使得，得到的最优解为，。关键词：最速下降法 无约束非线性规划 最优解一、问题重述用最速下降法求解无约束非线性规划问题：，设初始点取为,迭代到满足允许误差=0.01为止的精确解。二、问题分析2.1 无约束非线性问题的最优条件该问题是一个无约束非线性规划问题，利用最少下降法求解该问题，无约束非线性规划问题的最优解所要满足的必要条件和充分条件是我们设计算法的依据，为此有以下几个定理。定理1 设f：在点处可微，若存在，使，则向量p是f在点处的下降方向。定理2设f：在点处可微，若是无约束问题的局部最优解，则有
3、数学分析中我们已经知道，使的点x为函数f的驻点或平稳点。函数f的一个驻点可以是极小值点；也可以是极大值点；甚至也可能既不是极小值点也不是极大值点，因此称它为函数f的鞍点，以上定理告诉我们，是无约束问题的局部最优解的必要条件是：是其目标函数f的驻点。定理3 设f：在点处的Hesse矩阵存在，若，并且正定，则是无约束非线性问题的严格局部最优解。一般而言，无约束非线性问题的目标函数的驻点不一定是无约束非线性问题的最优解，但对于其目标函数是凸函数的无约束凸规划，下面定理证明了它的目标函数的驻点就是它的整体最优解。定理4设f：，f是上的可微凸函数。若有，则是无约束问题的整体最优解。2.2最速下降法的基本
4、思想最速下降法又称为梯度法，是1847年由著名数学家Cauchy给出的，他是解析法中最古老的一种，其他解析方法或是他的变形，或是他的启发得到的，因此它是最优化方法的基础。设无约束非线性规划问题中的目标函数f：在点处可微。最速下降法的基本思想是：从当前点出发，取函数在点处下降最快的方向作为我们收索方向，由的Taylor展示知 ，略去的高阶无穷小项不计，可见取时，函数值下降的最多 ，于是，我们可以够造出最速下降法的迭代步骤。2.3无约束非线性规划问题的迭代步骤解无约束非线性规划问题的最速下降法计算步骤第1步 选取初始点，给定终止误差 0，令k=0；第2步 计算，若，停止迭代，输出，否则进行第3步；
5、第3步 取；第4步 进行一维搜索，求，使得，令，k=k+1。转第2步。由以上计算步骤可知，最速下降法迭代终止时，求得的是目标函数驻点的一个近似点。 三、问题求解3.1对原无约束非线性规划迭代首先进行第一次迭代 则 令，即=0，解得：所以，此时，又因为 则进行第二次迭代则 令，代入即可解得：所以此时，又因为 则进行第三次迭代则 令，代入即可解得：所以此时，又因为 则进行第四次迭代则 令，代入即可解得：所以此时，又因为 以上仍然没有达到要求，即还需继续迭代，直到满足为止。3.2对原无约束非线性规划进行c+语言编程求解就这样无限的迭代下去，直到为止，为此，我们可以用c+语言编程得到，其算法设计如下图
6、（图3-1）停取，k:=0计算是否求令k:=k+1图（3-1）利用c+语言编程（源代码如下）：#include#includedouble lambda(double x2,double p2,double a2)double lam1,lam2;lam1=4*(pow(a0,3)*x0*x0+pow(a1,3)*x1*x1);lam2=-4*(pow(a0*x0,2)+pow(a1*x1,2);double s;s=-lam2/(2*lam1);return s;void main()double lamb,x2,a2,p2,g2,e=0.01,y; int i=0; x0=4.0; x1=
7、4.0;cout输入函数的系数：a0,a1:endl;for(i=0;iai; p0=2*a0*x0; p1=2*a1*x1; g0=-p0; g1=-p1;i=0; /开始迭代将次数赋值为0coute&i=200) lamb=lambda(x,g,a); x0=x0+lamb*g0; x1=x1+lamb*g1; p0=2*a0*x0; p1=2*a1*x1; g0=-p0; g1=-p1; i+; /cout(di %d ci mo=%f x1=%ftx2=%ftbuchang a=%fn,+i,sqrt(g0*g0+g1*g1),x0,x1,a); cout迭代次数为i p的模sqrt(
8、g0*g0+g1*g1)endlx的值x0 x1endlendl;y=a0*x0*x0+a1*x1*x1; coutendl分别输出x1,x2：x0 x1endl及极小值y： yendl;3.3结果分析运行即可得到如下结果：由以上结果可以得出：要想使得，则需要迭代6次，此时，。参考文献1范玉妹，徐尔.数学规划及其应用 M.北京：冶金工业出版社，2009.9.2倪勤. 最优化方法与程序设计 M.北京：科学出版社，2009.6.3孙文瑜，徐成贤，朱德通. 最优化方法 M.北京：高等教育出版社，2004.8.4徐宝文. C+程序设计语言 M.北京：机械工业出版社，2009.8.5冉崇善. C+程序设计语言设计教程 M.北京：机械工业出版社，2009.12.6张建勋，纪纲. C程序设计语言设计教程 M.北京：清华大学出版社，2008.2.1
