数值分析—第2章 非线性方程数值解法课件.ppt

第2章 非线性方程的数值解法 Numerical Solutions of Nonlinear Equations;求 f (x) = 0 的根;f(x)=(x-x*)mg(x);引入 逐步搜索法;例 求方程;2.1 对分区间法（二分法）;这时有三种情况：;有下述误差估计式;误差 分析：;二分法的计算步骤:;例2-1 求方程 f(x)=x3-x-1=0在区间（1.0, 1.5）内的一个实根，要求准确到小数点后的第二位。;n;①简单; ② 对f (x) 要求不高(只要连续即可)，且误差估计容易 .;2.2 迭代法;f (x) = 0;分别按以上三种迭代公式，并取x0=1.5进行迭代计算，结果如下：;K;迭代法的几何意义:;x;k;证明：① ?(x) 在[a, b]上存在不动点？;④;例2-4 试用不动点定理分析例2-2中同一方程的，3种不同迭代函数?1 , ?2 , ?3是否满足收敛条件.;（2）对于;若存在 的某个邻域 ;k;定义2-2 设迭代过程xk+1=?(xk)收敛于方程x=?(x)的根x*，如果存在常数p≥1和非零常数C，使得;定理2-3 如果x*是?(x) 的不动点，??(x)在x*的邻域连续，且 ?? (x)≠0 ，;定理2-4 如果x*是?(x) 的不动点，对于整数p1 ，迭代函数?(x)及其p阶导数在x*的邻域上连续，且满足;因此对迭代误差，令k→∞时有;钨监肝撰母题久款伍祟茫马骇肆须昏丙滞峙走及苔挤雨敝材着哄敌拉陷诲数值分析—第2章 非线性方程数值解法课件数值分析—第2章 非线性方程数值解法课件;2.3 加速收敛迭代法;则序列{xk}满足;2.3.2 Steffensen(斯蒂芬森)迭代法;? Aitken 加速：;定理2-6 （1）若x*为? 的不动点， ?? 在x*为领域连续且满足?? ( x*)≠1，则x*是Ψ(x)的不动点；反之，若x*为Ψ(x)的不动点，则x* 也是?(x) 的不动点。 （2）若x*是? 的不动点， ???? (x)在x*领域连续，且??( x*)≠1，则 Steffensen迭代法至少2阶收敛. ;例2-6 求解方程 在x0 =1.5附近的根。;2.4 Newton迭代法;　　设xk是根x*的某个近似值，过曲线y=f(x)上横坐标为xk的点Pk引切线，并将该切线与x轴交点的横坐标xk+1作为x*的新的近似值. ;解 令f(x)=xex-1，则f?(x)=ex(1+x) ，;定理2-7 设f?C2[a, b], 若x*为f(x)=0在[a, b]上的根，且f?(x*)?0，则牛顿迭代法是二阶收敛的，且;由 Taylor 展开：;由定理2-3;? 问题;Q2: 当x* 是 f (x)=0的m重根, 是否平方收敛？;Q3: 如何加速重根的收敛？;2.4.2 牛顿下山法;　　例如，用牛顿法求解方程 x3-x-1=0，在x=1.5附近的一个根x*. ;　　为了防止迭代发散，我们对迭代过程再附加一项要求，即具有单调性.;　　若用此法解方程x3-x-1=0 ，当x0=0.6时，求得x1=17.9，它不满足条件| f (x1) | | f (x0) | ，通过λ逐次减半进行试算,当λ=1/32时可求得x1=1.140625. 有f(x0)=-1.384, f(x1)=-0.656643, 显然|f(x1)||f(x0)|. 计算x2,x3,?时，均能使条件| f (xk+1) | | f (xk) 成立. 得到x4=1.32472即为x*的近似.;2.5 正割法（弦截法）;　　定理6 假设f(x)在根x*的邻域内△: |x-x*|≤δ具有二阶连续导数，且对任意x?△有f?(x)≠0，所取的初值x0, x1?△，那么当邻域△充分小时，弦截法(5.2)将按阶;　　因为上式用到前两点xk-1和xk的值，故此方法又称为双点割线法.;　　例题 用牛顿迭代法和割线法求方程 f(x)=x4+2x2–x–3=0， 在区间(1, 1.5)内之根(误差为10-9).;例 用快速弦截法求方程xex-1=0的根. 设方程的两个初始近似根为x0=0.5 , x1=0.6.;例 用快速弦截法求方程xex-1=0的根. 设方程的两个初始近似根为x0=0.5 , x1=0.6.; 单变量非线性方程求解：[x,fval]=fzero(f,x0,tol) 该函数采用迭代法计算函数f(x)的一个零点，迭代初值为x0，当两次迭代结果小于tol时停止迭代过程。tol的缺省值是eps。 注意: 在调用函数fzero 之前，要使用m文件建立自己要计算的