非线性方程和非性方程组的解法.PPTVIP

《实用数值计算方法》*4.8.1《实用数值计算方法》*4.8.2算法思想小结1.迭代法，而不是解析法2.线性化，极值化方法3.降维法4.下降法5.插值逼近法6.组合法7.梯度方向共轭方向改变度量《实用数值计算方法》*第四章习题《实用数值计算方法》*第四章*《实用数值计算方法》*4.6.3变尺度法VariableMatrixMethods算法框图：图4.29变尺度法算法框图《实用数值计算方法》*4.6.3《实用数值计算方法》*4.6.3《实用数值计算方法》*4.6.3《实用数值计算方法》*4.6.3《实用数值计算方法》*4.6.3表4.5各种方法比较《实用数值计算方法》*4.7直接法（Simplex,Powell）大量的目标函数是很复杂的，有时连解析式都没有，因而它的导数?f(x)很难求，有时甚至不存在。4.7.1单纯形法SimplexMethodNelder--Mead(1965)提出这种简单的方法。它不需要求导数（梯度）对变元不多的情况是有效的。程序简单。《实用数值计算方法》*4.7.1单纯形的思想是在n维空间的(n+1)个点(它们构成单纯形)上引进函数值比较。丢弃最坏的点并代之以新点。它们仍然构成单纯形。以此逐步逼近极小点。《实用数值计算方法》*4.7.1图4.30单纯形法中的反射《实用数值计算方法》*4.7.1图4.31单纯形法中的延伸《实用数值计算方法》*4.7.1《实用数值计算方法》*4.7.1图4.32单纯形法中的收缩《实用数值计算方法》*4.7.1e)缩小边长图4.33单纯形法中的缩小边长《实用数值计算方法》*4.7.1单纯形法（Simplex）框图：解x*?x0图4.34单纯形法计算框图《实用数值计算方法》*以上的迭代过程直到满足精度为止。精度：则x0作为所求的近似解。4.7.2Powelll方法Powelll方法是一种不依赖于目标函数梯度的直接搜索法。它逐步构造共轭方向并作为搜索方向，因此Powell方法也是一种共轭方向法。它的基本过程如下：《实用数值计算方法》*4.7.2图4.35Powell搜索路径表4.6Powell方法解题过程5.02.5《实用数值计算方法》*4.7.2《实用数值计算方法》*4.7.2Powell方法过程图示：图4.36Powell方法计算过程图示《实用数值计算方法》*4.7.2?循环上面(1)--(3)，直至P0点函数值不再减小为止。?当循环k次(k?n)以后，un与它前面的k-1个向量un-k+1,?,un-1共轭。因此对于二次函数，理论上只要循环n次即可求得极小值。即具有二次收敛性。事实上，因为P0和Pn是沿相同方向un求得的极小值，所以Pn?P0与un方向共轭。图4.37共轭方向《实用数值计算方法》*4.7.2图4.38Powell方法计算过程示意《实用数值计算方法》*4.7.2表4.7Powell方法第一次循环计算结果《实用数值计算方法》*图4.39单纯形法求一维极值示意图（1）4.7.2《实用数值计算方法》*图4.40单纯形法求一维极值示意图（2）4.7.2《实用数值计算方法》*4.7.2?但是，实际计算中对二次函数也不能保证n步内达到极小值点。因为每一循环都用Pn--P0“挤掉”u1,所以新的向量系ui(I=1,…,n)有可能线性相关，例如，某一循环中，如果?1?０则这样，u2,u3,…,Pn--P0是线性相关的。当发生这种情况时，以后的搜索就在n维的子空间中进行。最后的解就不正确。解决的办法是Pn--P0不是挤掉u1。而是挤掉ur,而?r?０。一般取最大下降方向（theDirectionoftheLarge