最优化理论与算法.docx

最优化理论与算法最优化理论与算法笔记 在老师的指导下，我学习了最优化理论与算法这门课程。最优化理论与算法是一个重要的数学分支，它所研究的问题是讨论在众多方案中什么样的方案最优以及怎样找出最优方案。 由于生产和科学研究突飞猛进的发展，特别是计算机的广泛应用，使最优化问题的研究不仅成为了一种迫切的需要，而且有了求解的有力工具，因此迅速发展起来形成一个新的学科。至今已出现了线性规划、整数规划、非线性规划、几何规划、动态规划、随机规划、网络流等许多分支。 整个学习安排如下，首先介绍线性与非线性规划问题，凸集和凸函数等基本知识及线性规划的基本性质；然后再这个基础上学习各种算法，包括单纯形法、两阶段法、大M法、最速下降法、牛顿法、共轭梯度法等，以及各种算法相关的定理和结论；最后了解各种算法的实际应用。 主要学习的基础知识： 1、一般线性规划问题的标准形式 minnåcxjj=1ijjnjs.tåaxj=1=bi,i=1,.,m,xj³0,j=1,.,n.学会引入松弛变量将一般问题化为标准问题；同时掌握基本可行解的存在问题，通过学习容易发现线性规划问题的求解，可归结为求最优基本可行解的问题。 2、熟练掌握单纯形法、两阶段法和大M法的概念及其计算步骤。 单纯形法是一种是用方便、行之有效的重要算法，它已成为线性规划的中心内容。其计算步骤如下： 1）解BxB=b,求得xB=B-1b=b，令xN=0,计算目标函数值f=cBxB； 2)求单纯形乘子w，解wB=cB ，得到w=cBB-1； 3)解Byk=pk，若yk£0，即yk的每个分量均非正数，则停止计算，问 题不存在有限最优解，否则，进行步骤； 4)确定下标r，使步。 两阶段法：第一阶段是用单纯形法消去人工变量，即把人工变量都变换成非基变量，求出原来问题的一个基本可行解；第二阶段是从得到的基本可行解出发，用单纯形法求线性规划的最优解。 大M法：在约束中增加人工变量xa，同时修改目标函数，加上罚项MeTxa，其中M是很大的正数，这样，在极小化目标函数的过程中，由于M的存在，将迫使人工变量离基。 3、掌握最速下降法的概念及其算法，并且能够讨论最速下降算法的收敛性。掌握牛顿法，能够熟练运用牛顿迭代公式：x(k+1)brb=minryrk>0，得到新的基矩阵B，返回第一 yrkyrk=x(k)-Ñ2f(x(k)(x-x(k) ，掌握共轭梯度法及其相关结论，以及其收敛性的讨论，掌握最小二乘法及其基本步骤。 最速下降法：迭代公式为x(k+1)=x(k)-lkd(k)。 (1)n计算步骤：1）给定点xÎR，允许误差e>0,臵k=1； 2）计算搜索方向d(k)=-Ñf(x(k)； (k)£e，则停止计算，否则，从x(k)出发，沿d(k)进行一维搜索，3）若d(k)(k)(k)(k)求lk，使f(x+lkd)=minf(x+ld)； l³04）令x (k+1)=x(k)-lkd(k)，臵k:=k+1，转步骤。 共轭梯度法：共轭梯度法是把共轭性与最速下降法相结合，利用已知点处的梯度构造一组供各方向，并沿这组方向进行搜索，求出目标函数的极小点，根据共轭方向的基本性质，这种方法具有二次终止性。 具体求解方法：首先，任意给定一个初始点x(1)，计算出目标函数f(x)在这点的梯度，若g1=0，则停止计算；否则，令d(1)=-Ñf(x(1)=-g1，沿方向d(1)搜索，得到点x(2)。计算在x(2)处的梯度，若g2¹0，则利用-g2和d(1)构造第二个搜索方向d(2)，再沿d(2)搜索；若已知点x(k)的搜索方向d(k)，得到)x(k+1)=xk(+lkdk(，再从x(k+1)出发，沿方向d(k+1)搜索。 结论：由共轭梯度法产生的方向d，d步必达到极小值。注：初始方向比为d(1)(1)(2)，d(m)是A共轭的，经有限=-Ñf(x(1)。 牛顿修正法：对于原始牛顿法和阻尼牛顿法都有共同缺点，一是可能出现2Hesse矩阵奇异的情形，因此不能确定后继点；二是即使 Ñf(x)非奇异，也未必正定，因而牛顿方向不一定是下降方向，这就可能导致算法失效。 牛顿修正法的策略，记d(k)=x-x(k)，得到Ñ2f(x(k)d(k)=-Ñf(x(k)，阻(k)尼牛顿法所用搜索方向是上述方程的解：d=-Ñ2f(x(k)-1Ñf(x(k)，这里假设2(k)-12(k)逆矩阵Ñf(x)存在。解决Hesse矩阵Ñf(x)非正定问题的基本思想是，2(k)(k)(k)修正Ñf(x)，构造一个对称矩阵Gk，得到方程Gkd=-Ñf(x)，解此方程，得到在点x(k)处的下降方向d(k)=-Gk-1Ñf(x(k)，再沿此方向作一维搜索。 2Tf(x)=(Ax-b)(Ax-b)，åii=1m最小二乘法：线性最小二乘问题 F(x)=求F(x)平稳点，令ÑF(x)=2ATAx-2ATb=0，则F(x)的平稳点满足ATAx=ATb，所以F(x)的平稳点为x=(ATA)-1ATb，由于F(x)为凸函数，所以x比为全局最小点。 非线性最小二乘 基本思想是通过解一系列线性最小二乘问题求非线性问题的解。基本步骤如下： 1）给定初点x(1)，允许误差e>0，臵k=1； éf1(x(k)ùê(k)úf(x)ú2）计算函数值fi(x(k)(i=1,2,.,m)，得到向量fk=ê2，再计算一ê.úêú(k)êëfm(x)úû¶fi(x(k)阶偏导数aij=,¶xii=1,.,m;j=1,.,n,得到矩阵Ak=(aij)m´n。 T(k)f，求得Gauss-Newton方向d(k)。 3）解方程组AkTAkd=-Ak4）从x(k)出发，沿d(k)作一维搜索，求步长lk，使得 F(x(k)+lkd(k)=minF(x(k)+ld(k). l令x(k+1)=x(k)+lkd(k). (k+1)-x(k)£e，则停止计算，得解x=x(k+1)，否则臵k:=k+1，返 5）若x回步骤二。 在整个学习过程中，我们学习的侧重点都是以应用为主，尤其是对各种相应的算法。而这个内容的学习除了可以解决很多实际问题外，对学习数值计算这门课程的帮助也是非常大的，其中所讨论的各种算法，对数值计算中求解线性方程组的解也能适用，这就将两门课程有机的结合在一起。 当然在学习过程中，有很多问题也是我没人理解透彻的，需要今后继续努力改进的。比如像各种算法的收敛性的讨论，以及在计算机中的编程实现等都将是我以后长期的课题。