3.1.非线性方程求解.pdf

程序设计与问题求解 ——非线性方程求解 非线性方程求解 • 非线性方程是许多学科用来表达客观规律的数学工具； • 非线性方程求根是个难题 ，一般无解析方法，仅有数值解法； • 某些时候，还需要求非线性方程组的根 ，这就更加困难。 本节介绍只有一个未知数的单个方程求根问题 非线性方程常见类型 一、代数多项式方程 a a x  a xn 0 0 1 n 例如： 2 x 5x 1 0 二、超越方程 凡是左侧不是代数多项式的方程，都称为超越方程 例如：2x 5x 2 0 非线性方程求解 定理 若函数 f (x ) 在 [a,b] 上连续（即f C [a ,b ]），且 f (a )f (b ) 0 ，那么方程 f (x ) 0 在区间 (a,b) 内 至少有一个根 非线性方程求解——二分法 二分法思想: 通过判断函数值的符号 ，将有根区间逐步二等分, 直到区间长度缩短到容许误差范围之内 非线性方程求解——二分法  求 f(x)=0 的根，二分法核心代码如下： while ( fabs(a-b)=0.000001 ) { c=(a+b)/2; // 计算中点 fa=f(a); // 计算f(a) fmid=f(c); //计算f(c) if ( fmid==0 ) // c 就是根，退出 break; else if (fa*fmid 0) // f(a)和f(c)同号，移动a a=c; else // f(a)和f(c)异号 ，移动b b=c; } printf(方程的解是：%6.6f,(a+b)/2); 非线性方程求解——牛顿法 假设求 f (x )=0 的一个近似根 ，且f (x )  0 以及 f (x ) 不变号 非线性方程求解——牛顿法  算法推导 首先给出经过 (x , f (x )) 切线方程 ： k k y f (x k )  f (x k )(x x k ) 该直线过x 轴上的点 (xk 1,0) ，于是 0 f (x k ) f