第五章 非线性方程的数值解法.pptVIP

第五章 非线性方程的数值解法 §5.1 二分法 §5.2 迭代法 §5.3 牛顿法 §5.4 弦截法 引 言 ※非线性方程： 指次数不低于二的代数方程和超越方程。 ※代数方程： 16世纪：找到三次、四次方程的求根公式； 19世纪：证明了5次以上的一般代数方程无公式解。 ※解题步骤： §5.1 二分法 §5.2 迭代法 §5.3 牛顿法 §5.4 弦截法 * * ※工程中常有此类方程，例如：在风道设计计算 中，通常采用柯氏公式： ——5.1 №1 确定初始区间： 一般题中给出或画出函数图形或由实际意义得到。 №2 通过二分法求得一个比较粗糙的近似值。 №3 利用其它解法求得精度较高的近似值。 一、二分法： a b x1 x2 a b 或 x* 图5－1 二分法的几何意义 ※几何意义： ［例］ ［解］ (2.09375,2.109375) f(2.109375)0 6 (2.09375,2.125) f(2.09375)0 5 (2.0625,2.125) f(2.0625)0 4 (2,2.125) f(2.125)0 3 (2,2.25) f(2.25)0 2 (2,2.5) f(2.5)0 1 (2,3) f(3)=160 f(2)=-10 有根区间 函数值符号 n 二、算法与流程图： 1.算法： 2.流程图：（略） 3.C－程序： # include math.h # include stdio.h # define A 2 # define B 3 float f(x) float x; { float y; y=x*x*x-2*x-5; return(y); } main() { float a=A,b=B,x; x=(a+b)/2; if(f(x)=1e-8) printf(The root is %f\n,x); while(fabs(f(x))1e-5) { if(f(a)*f(x)0) b=x; else a=x; x=(a+b)/2; } printf(x=%f\n,x); } 三、误差估计定理：（证明略） 如前例误差限为： ♀优点：程序简单，对函数性质要求低，只要连续就 可以了，收敛速度不算太低。 ♀缺点：不能求偶数重根，也不能求复根和虚根。 ♀ 迭代法是一种逐次逼近法，可用来求解代数方程、超越方程、微分方程、方程组以及求矩阵的特征值等。存在收敛性及收敛快慢的问题。 ♀ 已知非线性型方程的一个近似根后，可用迭代法使这个根逐步精确化，一直到满足精度为止。 一、迭代公式： 1.改写： ——5.2 ［如］ 2.迭代： 3.迭代序列的收敛性与几何解释： ——5.3 x y y = x x* y=g(x) x0 p0 x1 p1 ? x y y = x x* y=g(x) x0 p0 x1 p1 ? 图5－2 迭代法的几何意义 二、收敛条件： ★定理5.1： 证明： 证毕。 ※ 为了控制迭代次数，我们给出如下定理。 ——5.4 ★定理5.2： ——5.5 ※ 上机时，可利用(5.5)结束运算。 ※ 上述定理条件可减弱为 ♀优点：逻辑结构简单 （证略）。 ［例］ ［解］ 0.0111 0.57117 5 -0.00065 0.56691 10 -0.1964 0.56006 4 0.00115 0.56756 9 0.34446 0.57970 3 -0.00203 0.56641 8 -0.61929 0.54524 2 0.00358 0.56884 7 0.10653 0.60653 1 -0.00631 0.56486 6 0.5 0 n n ［例］ 这两种过程都收敛吗？ ［解］ ♀将非线性函数线性化； 在单根附近收敛快，且可计算代数方程的复根。 一、公式及误差分析： 1.公式： ——5.6 这就是牛顿(Newton)法的迭代公式。 ♀优点：收敛速度快，可计算复数根。 ♀缺点：需要计算导数值。 2.几何意义： 图5－3 牛顿法的几何意义 x y x* x0 3.收敛性： 由于 ——5.7 根据上节的结果，则可知牛顿法的迭代公式收敛的一个 充分条件是 ——5.8 ——5.9 利用微分中值定理知 具有线性收敛性。 代入（5.9）得 对(5.7)求导知 于是有 由于误差的这一特点而称牛顿法具有平方收敛性。 ［例］ ［解］ 0.56714 0.56716 0.57102 0.5 3 2 1 0 n ☆ 收敛速度远高于前一种迭代格式 二、算法与流程图： 1.流程图：（略） 2.C程序： #include math.h #include stdio.h #def