化工数值计算-2探究.ppt

行分析,这样才能获得合乎实际的合理性数据。 2.2.5　弦截法 牛顿迭代法收敛速度快,但它要求计算f(x)的值。在科学与工程计算中,常会碰到f(x)不易计算或算式复杂而不便计算的情况。下面介绍一种既具有牛顿法恒收敛和收敛速度快的优点,又不用求导数f(x)的弦截法。 2.2.5.1　方法概述 弦截法的基本思想与牛顿法相似,即将非线性函数f(x)线性化后求解。两者的差别在于弦截法实现函数线性化的手段采用的是两点间的弦线,而不是某点的切线。 设方程式f(x)= 0的一个近似根为x0,另一个近似根为x1。若用过点[x0,f(x0)]和[x1,f(x1)]间的弦线代替曲线y=f(x)(见图2-18),得直线方程为 ? 令式(2-54)中f(x)= 0可得 ? 将上式中的x作为方程f(x)=0的又一新的近似根,从而有如下迭代格式: ? 即弦截法的迭代计算公式。显而易见,上式即为用差商 ? 代替牛顿公式(2-44)中的导数所得的结果。 弦截法的几何意义就是通过点Pk-1和Pk的弦割线与x轴的交点xk+1不断逼近真实根x*的过程,如图2-18所示。 与牛顿法只需给出一个初值不同,弦截法需要给出两个迭代初值x0 和x1。如果与逐步扫描法结合起来,则最后搜索的区间的两个端点值常可作为初值x0和x1。 弦截法虽比牛顿法收敛速度稍慢,但在每次迭代中只需计算一次函数值,又不必求函数的导数