 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
数值分析第五版答案
下载积分：0
内容提示：数值分析第五版答案
文档格式：PDF|
浏览次数：124|
上传日期： 12:33:53|
文档星级：
全文阅读已结束，此文档免费下载
下载此文档
该用户还上传了这些文档
数值分析第五版答案
关注微信公众号 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
数值分析第五版_李庆扬__课后习题答案[内部资料]
下载积分：999
内容提示：数值分析第五版_李庆扬__课后习题答案[内部资料]
文档格式：PDF|
浏览次数：242|
上传日期： 04:07:44|
文档星级：
全文阅读已结束，如果下载本文需要使用
 999 积分
下载此文档
该用户还上传了这些文档
数值分析第五版_李庆扬__课后习题答案[内部资料]
关注微信公众号 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
数值分析-李红-华中科技大学出版社
下载积分：0
内容提示：数值分析-李红-华中科技大学出版社
文档格式：PDF|
浏览次数：228|
上传日期： 23:19:39|
文档星级：
全文阅读已结束，此文档不支持下载
该用户还上传了这些文档
数值分析-李红-华中科技大学出版社
关注微信公众号您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
哈工大数值分析上机报告.doc 66页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
需要金币：100 &&
哈工大数值分析上机报告
你可能关注的文档：
··········
··········
哈工大数值分析上机报告
哈工大-数值分析上机实验报告
实验报告一
非线性方程求解
摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。
前言：（目的和意义）
掌握二分法与Newton法的基本原理和应用。
数学原理：
对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。
对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)&0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)&0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式
xk?1?xk?f(xk) 'f(xk)
产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为
xk?1?xk?rf(xk) f'(xk)
其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。
程序设计：
本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下
正在加载中，请稍后...查看: 2049|回复: 6
什么是计算数学
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
才可以下载或查看，没有帐号？
现代的科学技术发展十分迅速，他们有一个共同的特点，就是都有大量的数据问题。&BR&&BR&比如，发射一颗探测宇宙奥秘的卫星，从卫星世纪开始到发射、回收为止，科学家和工程技术人员、工人就要对卫星的总体、部件进行全面的设计和生产，要对选用的火箭进行设计和生产，这里面就有许许多多的数据要进行准确的计算。发射和回收的时候，又有关于发射角度、轨道、遥控、回收下落角度等等需要进行精确的计算。&BR&&BR&有如，在高能加速器里进行高能物理试验，研究具有很高能量的基本粒子的性质、它们之间的相互作用和转化规律，这里面也有大量的数据计算问题。&BR&&BR&计算问题可以数是现代社会各个领域普遍存在的共同问题，工业、农业、交通运输、医疗卫生、文化教育等等，那一行那一业都有许多数据需要计算，通过数据分析，以便掌握事物发展的规律。&BR&&BR&研究计算问题的解决方法和有关数学理论问题的一门学科就叫做计算数学。&BR&&BR&计算数学属于应用数学的范畴，它主要研究有关的数学和逻辑问题怎样由计算机加以有效解决。&BR&&BR&计算数学的内容&BR&&BR&计算数学也叫做数值计算方法或数值分析。主要内容包括代数方程、线性代数方程组、微分方程的数值解法，函数的数值逼近问题，矩阵特征值的求法，最优化计算问题，概率统计计算问题等等，还包括解的存在性、唯一性、收敛性和误差分析等理论问题。&BR&&BR&我们知道五次及五次以上的代数方程不存在求根公式，因此，要求出五次以上的高次代数方程的解，一般只能求它的近似解，求近似解的方法就是数值分析的方法。对于一般的超越方程，如对数方程、三角方程等等也只能采用数值分析的办法。怎样找出比较简洁、误差比较小、花费时间比较少的计算方法是数值分析的主要课题。&BR&&BR&在求解方程的办法中，常用的办法之一是迭代法，也叫做逐次逼近法。迭代法的计算是比较简单的，是比较容易进行的。迭代法还可以用来求解线性方程组的解。求方程组的近似解也要选择适当的迭代公式，使得收敛速度快，近似误差小。&BR&&BR&在线性代数方程组的解法中，常用的有塞德尔迭代法、共轭斜量法、超松弛迭代法等等。此外，一些比较古老的普通消去法，如高斯法、追赶法等等，在利用计算机的条件下也可以得到广泛的应用。&BR&&BR&在计算方法中，数值逼近也是常用的基本方法。数值逼近也叫近似代替，就是用简单的函数去代替比较复杂的函数，或者代替不能用解析表达式表示的函数。数值逼近的基本方法是插值法。初等数学里的三角函数表，对数表中的修正值，就是根据插值法制成的。&BR&&BR&在遇到求微分和积分的时候，如何利用简单的函数去近似代替所给的函数，以便容易求到和求积分，也是计算方法的一个主要内容。微分方程的数值解法也是近似解法。常微分方程的数值解法由欧拉法、预测校正法等。偏微分方程的初值问题或边值问题，目前常用的是有限差分法、有限元素法等。&BR&&BR&有限差分法的基本思想是用离散的、只含有限个未知数的差分方程去代替连续变量的微分方程和定解条件。求出差分方程的解法作为求偏微分方程的近似解。&BR&&BR&有限元素法是近代才发展起来的，它是以变分原理和剖分差值作为基础的方法。在解决椭圆形方程边值问题上得到了广泛的应用。穆恰，有许多人正在研究用有限元素法来解双曲形和抛物形的方程。&BR&&BR&计算数学的内容十分丰富，它在科学技术中正发挥着越来越大的作用。&BR&laixi:博士论坛
[ 本帖最后由 咕噜噜 于
11:43 编辑 ]
发信人: likemath (libai), 信区: Mathematics &BR&标&&题: 侃侃计算数学 （数值优化） &BR&发信站: 北大未名站 (日11:36:15 星期二), 转信 &BR&&&&BR&&&&BR&侃侃计算数学 （数值优化） &BR&&&&BR&谈到数值优化，不能不提的是单纯形算法。这被誉为20世纪最受人欢迎的算法为人们带来 &BR&&&&BR&巨大的经济效益。不过有趣的是，这个最好的算法在算法复杂度理论里面却解释不通， &BR&因为它不是多项式算法。 &BR&&&&BR&数值优化以求解有约束或无约束条件下函数最值为目标。我想数值优化里面最 &BR&令人头疼的是如何判断你找到的不是极值而是最值。因为二者的区分 &BR&似乎只能从函数值上得到，其它的信息包括各种导数似乎都没有什么区别 &BR&&&&BR&。但是，实际中的很多问题都有大量的极值点，如果挨个寻找根本不可能。 &BR&对付这个问题，现在最有效的武器 &BR&&&&BR&应该是随机算法包括遗传算法等等。但是，其庞大的计算量有时也让人望而却步。 &BR&优化里面另外一个困难的问题是整数优化，凡是涉及的整数的问题总是令人头疼的 &BR&，因为限制太为严格。直到 &BR&今天，人们连线性方程组的整数解都没有完全解决，何况在此基础上考虑整数规划等等。 &BR&&&&BR&其它的诸如不可微优化、非线性规划等等发展到今天似乎很难有什么突破，也局限于在 &BR&理论上推导满足一些条件的算法，但实际中有几个问题能满足这些条件（我的愚见，未必 &BR&正确）。 &BR&&&&BR&现在，与计算机组合优化密切相关的计算复杂度理论异军突起，新千年7个悬赏问题之一 &BR&就是与之相关的P是否等于NP.我想，结合图论组合优化计算机等学科，这一方面的发展是 &BR&很有空间的。 &BR&
说的对，数值优化最怕掉到局部最优的陷阱里去了
回复：（jacqueline）说的对，数值优化最怕掉到局部最...
&DIV class=quote&&B&以下是引用&I&jacqueline&/I&在 8:39:24的发言：&/B&&BR&说的对，数值优化最怕掉到局部最优的陷阱里去了&/DIV&
&br&但是要做到全局最优又挺困难的，大家在这方面有没有什么心得可以谈谈吗？
回复：（milly2005）回复：（jacqueline）说的对，数值...
&DIV class=quote&&B&以下是引用&I&milly2005&/I&在 15:59:25的发言：&/B&&BR&&BR&&BR&但是要做到全局最优又挺困难的，大家在这方面有没有什么心得可以谈谈吗？&/DIV&
&br&转一篇文章大家看看，《论全局最优》&BR&&BR&一、 从爬山谈起&BR&每个人都会迷路。但是也许很少有人会对迷路这种现象思考甚多。原因很简单，一个人对周围的环境不熟悉，自然就可能迷路。而且，越是生疏的地方，迷路的可能就越大。如果是自己非常熟悉的地盘，也许闭着眼睛也不会迷路。&BR&但其实迷路只是一种现象，却非本质。对人们来说，迷路不是偶然的，而是必然的。是否迷路，迷路的程度多少，完全取决于当事人对周围环境的信息的了解程度。比如，从一个地方A到另外一个地方B，其中有N条岔道，其中只有一条连接A和B；假设当局者对从A到B一无所知，那么，他要花多少时间才能到达目的地呢？答案是，当局者也许一辈子也找不到从A去B的路----如果他不在探索的时候记录信息的话。我说的意思是，如果当局者在寻找这条唯一的道路的时候，不能记录他以前的探索结果，不能记录这N条路的特征，那么，当局这就可能永远也找不到正确的从A到B的道路。&BR&这是一件令人伤心的事实。但毕竟，很幸运的是，我们是人，毕竟能够发挥主观能动性，在寻找这条道路的时候可以记录信息，从而在搜索N次以后，终于可以找到从A到B的路。但是在另外一些情况中，我们可能连这样的乐观都没有。比如，假设有人让某人在一座山脉中寻找一个海拔最低的山谷，但是这个山脉有多少山谷，怎样到达这些山谷，我们根本不知道。在这样的情况下，假设该人能够记录他搜索的历史信息，，那么他要花时间才能够断言他找到的是海拔最低的山谷呢？答案是，无穷时间。我的意思是，该人一辈子都不能够断言他找到的是最低的山谷。&BR&相似的情形在社会生活中很多领域出现。一件困惑人们千百年之久的事情就是，怎样才能找到最适合自己的对象？是否真的存在自己的那根肋骨？是否真的有对自己最好的人？我们可以悲观的说，人们是不可能找到最适合自己的对象的，也不存在自己的那根肋骨，不存在对自己最好的人。所有的这些我们都不知道，因为我们不能够一一了解存在的人类个体，不能亲身体验。这样的一个事实的结果是这样的一个忠告：惜取眼前人。正如苏格拉底让他的学生们所做的，在摘麦穗的时候，自己觉得最大了的就摘。别老是想着摘最大的，否则可能结果恰恰相反。&BR&这些问题就是全局最优的问题。这些问题涉及到社会生活的方方面面。因此，对这种问题作一个探讨是有必要的，也可能在某些方面能够改进人们的某些观念。比如，一个让喜欢美女的人们惋惜不已的问题是，人类为什么要一夫一妻，为什么不能够多取几个老婆，这样，至少也能让我们多亲芳泽，不至于私下里经常流着口水？我们将会看到，解释是哲学层面的，但是诠释确是非常合理。&BR&&BR&二、 局部最优和全局最优&BR&既然我们要考察全局最优的问题，那么一个必然要设计的概念就是局部最优。什么是全局最优，什么是局部最优呢？这样的解释可能比较容易理解：假设当事人在一定的范围内寻找一个解决方案，使得当事人期望的目标达到最大或者最小。如果一个方案使得当事人的目标得以实现，就称该方案为全局最优方案，而目标就称为全局最优的目标。假设一个方案虽然不能让当事人的目标达到全局最优，但是在一个小范围内，已经没有其他的方案比这个方案更好了，那么就成该方案为局部最优。用数学语言可以如下描述：&BR&设f(x)是当事人的目标函数，F是当事人的选择范围。我们希望找到使得f(x)最小的x。如果有一个x^*使得f(x^*)&=f(x),对任意的x属于F都成立，就称x^*使全局最优的。不然，如果存在一个集合E包含在F中，x^*属于E，使得f(x^*)&=f(x),对任意的x属于E都成立，就称x^*是局部最优的。&BR&求解全局最优的问题就称为全局最优化问题。由以上阐述可以看到，一个全局极小问题可以描述如下&BR&global min f(x) &BR&s.t. x 属于 F&BR&其中F={x|g(x)&=0}，g(x)是一个向量函数。&BR&假设f(x),g(x)都是曲线，那么全局最优解就是让f(x)在g(x)&=0限定的区域F内达到最低点的点。&BR&&BR&三、 全局最优的解法&BR&要考察如何求解全局最优问题，先从几何观点来审视这个问题是有好处的。假如在一个斜坡上放一个球，如果摩擦力适当小，那么球将会往更低的地方跑。从物理上我们知道，在一个斜角为v的斜坡上，物体的下滑力，如果忽略摩擦力，当为mgsin v。设斜坡斜率为K，则有&BR&(sin v)^2+(cos v)^2=1, sin v/ cost v= K&BR&由此可得&BR&sin v =K/(1+K^2)^{1/2}&BR&根据牛顿定律，可得&BR&mds/dt =mg sin t=m K/(1+K^2)^{1/2}=&&BR&ds/dt=K/(1+K^2)^{1/2}&BR&其中s是沿着斜坡运动的速度，所以ds/dt就是沿着斜坡运动的加速度。&BR&我们来看物体在平行地面（即可设为x轴）运动的加速度是多少。根据矢量分解的平行四边形法则，该加速度应该为&BR&dx/dt=ds/dt * cos v=K/(1+K^2)^{1/2}* 1/(1+K^2)^{1/2}=K/(1+K^2)&BR&现在考虑球放在曲线f(x)上在点x处的加速度。由于在该点的梯度为\nabla f(x),故可类似推得：&BR&dx(t)/dt=-\nabla f(x)/(1+||\nabla f(x)||^2)&BR&其中||.||表示\nabla f(x)这个向量的范数。如向量a的范数为||a||=(a^Ta)^{1/2}。&BR&注意到-\nabla f(x)是一个向量，而1+||\nabla f(x)||^2是一个数，所以方程&BR&dx(t)/dt=-\nabla f(x)/(1+||\nabla f(x)||^2)就表示沿着\nabla f(x)下降的运动方程。更进一步，我们可以略掉1+||\nabla f(x)||^2，而直接考虑：&BR&dx(t)/dt=-\nabla f(x) （A）&BR&对t进行离散，并用差分离散该方程，就可以得到诸如最速下降法（或者称为下山法），牛顿法之类的局部最优的求解方法，如，一阶离散，可得：&BR&x_{k+1}-x_k=-\nabla f(x_k)*(t_{k+1}-t_k)&BR&这就是最速下降法。&BR&如对方程（A）先对t求导数，即&BR&d^2x(t)/dt^2=-\nabla^2 f(x)*d x(t)/dt=-\nabla^2f(x)*\nabla f(x) (B)&BR&由此可得牛顿法。&BR&注意到当\nabla=0的时候，上述方程(A),(B)都处于稳定状态了，这时候球打到最低点，这时候就不能下降了。&BR&由上面可以看到，要求的一个局部最优还是很容易的。只要求系统(A),(B)的稳定点就可以了。&BR&但是求解全局最优就没有那么幸福了。因为这个时候我们根本就无法判断一个点是否全局最优。即使我们能够找到很多局部最优，我们也不能断言我们已经找到全局最优----除非函数f(x)只有一个局部最优。这是一个悲剧。&BR&一种很直接的想法是通过局部最优来求解全局最优。这又分为几种方法。一种十多取几个初始点，求得多个局部最优，然后比较这些局部最优，取其最好的。这种方法就叫做随机算法。这种算法很直接，缺点也很明显---我们无法判断我们是否达到了全局最优。因为很多局部最优只是一个简单的罗列，而没有任何实质性的改进。&BR&另外一种想法几何直观也很明显。这种方法希望在达到一个局部最优以后，能够跳出该局部最优，然后继续搜索，达到另外一个更好的局部最优点。这种方法包括填充函数法和隧道函数法。这样看起来，如果这种想法真正的能实现，那么，就能够实现全局最优了。因为假设f(x)在F的局部最优有限，通过该方法，一直找下去，总是能够达到全局最优。但是很不幸，这些方法是无法实现真正的跳出已经求得的局部极小的。填充函数法通过构造一个新的辅助函数，把局部极小变成极大，然后对该点进行随机扰动（或者是确定性的扰动），再利用经过扰动以后的点为初始点进行搜索，寻找该辅助函数的局部极小点，该局部极小点是原来函数f(x)的更好的局部极小点或者至少和原来一样好。问题就是出在对原局部极小点进行扰动上。该扰动是没有目的的，因而很盲目。而且，该扰动的方案有无穷多种，每种方案可能导致不同的结果，而在实际计算中又无法全部实现这些方案。如对x是一维变量，对局部极小点的扰动就有两种方案，对x是n维变量的情况，对局部极小点的扰动有无穷多种方案。而且，利用扰动以后的点求得的辅助函数的局部最优可能根本就使原来的局部最优点，或者在一些计算以后回到该局部最优点，从而形成循环。这样就意味着我们无法判断利用辅助函数搜索到的点是全局最优点。相似的情况发生在隧道函数上。隧道函数是通过一个辅助函数到达另外一个点上，该点在另外一个更好的局部最优点的吸收域中。它的本质方法是水平集法。隧道函数的缺点就在于构造隧道函数求另外一个点的时候，同样存在无穷多个搜索方向，因而也无法保证最终求得的局部最优为全局最优。&BR&积分水平集方法是另外一个求解全局最优化的方法。该方法的本质是这样的：如果找到了一个可行点x^*属于F，如果能够构造水平集E={x|f(x)&=f(x^*),x属于F}。在该区域中搜索，求得新的局部最优，一次迭代，就能够求得全局最优。构造水平集是一件困难的工作。为此，研究者们采用随机投点的方法来近似构造该水平集，并用f(x)&=f(x^*)该评判准则来剔除大于f(x^*)的点。这样，近似构造水平集就取决于随机所投的点。为了获得在F上全局最优，就需要投无穷多个点。随着投点的越来越多，计算量也越来越大。因此，从理论上来说，积分水平集方法能够求得全局最优，但是从计算角度来看，需要耗费非常多的计算量。&BR&有朋友可能会问，既然如此，为什么不对&BR&min f(x)&BR&s.t. f(x)&f(x^*)-eps,x属于F&BR&其中eps是一个很小的正数，这样一个问题来求解呢？这样不是在水平集E上求解么？但是，根据笔者用计算软件计算的经验，这样的问题始终还是停留在原来的局部最优点x^*，没有什么改进。背后的机理如何，现在还不清楚。&BR&最直观的方法是分割区域F来求解全局最优。这样就有区域分割算法、单纯分割算法等等方法。显然，区域分割方法为了求得全局最优，需要大量的计算量。&BR&启发式算法，包括遗传算法、神经网络算法、蚁群算法等等，都不能保证求得全局最优。而且计算量相当巨大。&BR&就目前来说，除了对某些特殊的全局最优问题能够有效求解，对于一般的问题，目前还没有有效的手段来求解。&BR&
回复：（TNC）什么是计算数学
四、 全局最优与正性问题&BR&全局最优问题和正性问题是等价的。所谓正性问题是指如下问题：&BR&[正性问题]：&BR&判断是否对任意的x满足g(x)&=0，都有f(x)&=0成立。显然，如果正性问题能够判定，那么因为全局最优问题等价于,对于全局最优的t=f(x^*)，其中x^*使全局最优解，有&BR&t-f(x)&=0,对任意的x满足g(x)&=0&BR&所以全局最优问题与正性问题等价。&BR&已经被Richardson证明，对于一般的函数类f(x)和g(x)，该正性问题是不可&BR&解的。所谓不可解，就是无法作出判定。因此，对函数f(x),g(x)作出限定是有必要的。&BR&如果把f(x),g(x)限定在多项式函数，那么很多结果已经在代数几何中得到。Hilbert提出的23个著名问题中的Hilbert’s17th problem就是判定是否对任意的多项式f(x)&=0,对于任意的x属于R^n，f(x)都是可以表示成有限多个有理多项式的平方的和。&BR&Hilbert’s 17th problem已经由Artin在1927年的出否定的解决。他证明，对于在实数空间R^n上恒非负的多项式f(x)，它可以表示成为有限个有理多项式的平方和，但是可以举出反例不可以表示为多项式的平方和。&BR&但是对以更广泛的问题f(x)&=0,\forall x:g(x)&=0这样的判定问题，其中g(x)是多项式函数，现在仍然有很多问题需要研究。如果g(x)十多项式向量函数，那么区域F={x|g(x)&0}就称为半代数集。因此一般的正性问题可以叙述如下：&BR&[正性问题]：&BR&判定，是否对任意的满足g(x)&=0的x，都可以推出f(x)&=eps?&BR&目前对该正性问题已经有不少结果。如Hilbert’s Positivstellensitz.，Schugmzen’s Positivstellensitz,Putinar’s Positivstellensitz,等等结果。但是已经证明，虽然正性问题能够判定，但是计算量却是关于变量x的维数n以及1/eps的指数。显然，这样的计算量是非常巨大的。&BR&用代数几何中的这些结果到最优化中，最近有若干研究。有两个人比较有代表性，他们是Parrilo和Lassarre。但是可以看到，用这样的方法求解多项式全局最优化，只能限于一些低维问题。&BR&和正性问题相关的有半无限最优化问题以及半无限最优化问题的特例robust最优化问题。这些问题都归结成正性问题。半无限最优化问题形式如下：&BR&min c^Tx&BR&s.t. f(x,s)&=0,\forall s属于U&BR&其中U是一个确定性的集合。显然，求解这类问题的关键在于能否判定对于任意的s属于U，都有f(x,s)&=0。或者说，满足该问题的x是否可以解析地描述出来。一般来说，这样的问题需要指数时间的计算量。&BR&&BR&五、 全局最优化、组合最优化与计算复杂性&BR&组合最优化在实际中有广泛的应用。大部分实际的最优化问题都归结成组合最优化问题。组合最优化考虑如下问题&BR&min f(x)&BR&s.t. x属于F&BR&其中F是一个离散点集。&BR&组合最优化只是全局最优化的特例。自从Cook1972年提出NP理论，衡量计算复杂性就有了标准。已经证明了相当大一类组合最优化问题是NP困难的问题和NP完全的问题。而这些问题有等价于一些全局最优化问题。因此，计算复杂性的概念就被引入了全局最优化。&BR&所谓P问题，就是可以被关于问题本身的参数，如维数，约束个数等的多项式时间内求解的问题。NP问题就是对于一个给定的点，能多项式时间内判定它是否给定问题的解的问题。NP包含P是以显然的事实。但是P是否也包含NP，就是一个非常困难的问题。目前这个问题被列为全世界7大数学难题之首。&BR&有一类NP问题，它们之间相互等价，求解其中一个问题就求解了全部问题。大部分组合组优化问题属于此类。这类问题称为NP完全问题类。单个问题就称为NP完全问题。一般相信，这类问题不存在多项式时间算法。&BR&由于全局最优化包含了组合最优化，因此也可以用计算复杂性的概念到全局最优化上面。能多项式时间求解的问题是容易的，不能多项式时间内求解的问题是困难的。目前这已经成了判断算法好坏的一个标准。&BR&&BR&六、 全局最优化与积分&BR&我们看到，求解全局最优化关键是利用全局信息。全局信息，就是能够描述函数整体性质的信息。导数、梯度、Hessian矩阵等等，都是局部信息。利用全局信息的方法，一个是水平集，一个是积分。利用水平集，如积分水平集方法，就是其中的代表。水平集的概念我们上面已经介绍过了，这里就不多说。这种方法最终归结到计算近似积分上。另外的一个方法，是直接计算积分。目前的一个结果是，如果给定一个点，判断它是否一个全局最优点，只要判断一个积分函数是否存在零点就够了。详细地说，给定x^*，考虑&BR&global min f(x)&BR&s.t. x属于[0，1]^n&BR&其中[0,1]^n表示n维空间上的一个箱子。&BR&通过构造辅助函数&BR&g(t)=int_{x属于[0,1]^n}e^{t(f(x^*)-f(x))}dx&BR&那么，可以证明g(t)是凸函数，并且，如果x^*不是全局最优，那么g(t)存在一个0点解t^*，而且当t&t^*时，g(t)是增函数。而如果x^*使全局最优，那么对g(t)的近似0点t^*，当t&t^*时，g(t)是减函数。利用该方法，就可以判断x^*是否全局最优。这也同时建立了全局最优与积分的关系。而且，也可以构造算法，来求解f(x)在[0,1]^n的全局最优。&BR&利用上述方法的时候要计算高维积分。而这仍然是一个困难问题。求高维积分一般采用的是概率方法。但这样就无法得到一个确定性的误差，这样，就不知道在计算积分的时候损失了多少。从这个意义上说来，我们无法采用上述方法来求解全局最优。但是，有一种基于数论的一致分布的方法来求高维积分，而且误差是确定性的。这样，就可以利用它来计算近似全局最优解。&BR&目前，探索积分与全局最优、利用积分来求解全局最优的方法和理论还很缺乏，有必要也有潜力大力发展。&BR&&BR&七、 从哲学观点看全局最优&BR&从能量守恒原理来看，物质运动都趋于使维持物质的能量最小这个方向发展。如肥皂泡在收缩过程中形成的曲面要使曲面面积最小，从而使支撑该曲面的能量最小；又如原子为了自身的稳定，总是释放能量，让电子轨道向低能级跃迁；其他，无论是生物学、进化论，还是天体演化，等等，都表明了这一原理的合理性，即使不敢说正确性。而耗散理论的观点也在这方面说明了这一点。&BR&物质运动既然趋向低能和稳定结构，那么，一般来说，不存在普遍的完美的解决方案，使得即可以解决问题，而解决的手段又漂亮的无懈可击。因为如果这样，就意味着解决方案处于高能状态，而这是不稳定的。只要条件稍微改变，旧有可能引致该解决方案发生质的变化，由完美变成不完美。完美的存在是高能耗的，不稳定的，其存在当且仅当在很严格的条件满足。如就最优化来说，人们目前能自豪地宣称能够找到全局最优的问题，也就是自己实现能判定出来的局部最优唯一的问题类----甚至对这样的问题，如果要用计算复杂性理论来度量的话，也还有一些问题不能够多项式时间求解。人类能作用的范围太小了。&BR&既然完美的活的需要高代价，那么寻找全局最优有必要向两个方向发展。第一个方向是保证能够求得全局最优，另外一个是能尽可能快地求得相对好的解。目前的研究状态就是这样的。前者导致了启发式算法、随机算法的热烈研究，而后者则导致了近似算法的研究。显然，我们也希望，在进行寻找全局最优的时候，不必要用各种准则来限制研究方法。“黑猫白猫，抓到老鼠的是好猫”。&BR&有时候我们常常通过外界映照，来推知自身所处的位置。如通过观察其他动物来推知人类本身所处于世界中的位置，通过别人的行为推知自己行为的合理性。对世界观察，并用数学同构的观点来认知世界，也许不失为拓广知识的一种方法。.但是这些对于无知觉的算法来说，都统统失效。智能算法，如能记录搜索过程，自适应调整搜索方向的算法，也许值得研究。&BR&但无论如何，我们处于一种尴尬的境地，一方面，我们需要全局最优，一方面，我们对大部分问题无能为力。也许这恰好是人类几千年来挣扎着发展的历史。&BR&&BR&八、 对生活的建议&BR&对别人的生活提出忠告是一件很冒险的事情。因为任何人不见得比别人好得更多---如果单纯从人在世界所处的位置来看。所以这里只是提出建议。&BR&在惆怅现在找对象还是将来找对象，现在找了，万一将来遇到更好的怎么办的朋友，建议赶快找一个对自己好，自己也相对喜欢的身边人。好没有穷尽，任任何时候也不能判断自己的选择是否最优，既然如此，不如赶快动手，别丢失了你现在最大的那颗麦穗。&BR&对勇于立志，要成大事业但是不知道干那个好，感觉什么都能做的朋友，不如从现在做起，从眼前做起，把当前要做的事情做好，同时也要智能地观察和调整自己的路。也许这样，才是这个人生的最优解。&BR&对研究感到失望，诚如若干年前我陷入的疑惑的朋友，不如先判断那种方法能有一些改进，然后就埋头进去苦干。人类的发展需要前赴后继的努力，能给大厦添砖加瓦就相当不错了。毕竟人能力有大小，做好本分也是值得尊重的。&BR&对决策影响到很多人、但是手中的权力无监督的朋友，请遇事情反复思量。一种决策很难达到最优，永远不要武断地认为自己已经找到最佳的解决途径了。不妨多不耻下问，也许至少能够改善决策，提高效果。 &BR&
Powered by}