第10章 数值分析_05(带matlab).pptVIP

* * 【案例10.3】用线性插值求 （x* = 10.723805） 解：设 ，取x0 = 100，x1 = 121 则 y0 = 10 y1 = 11 从而 Matlab(实现) yi = interp1([100 121],[10 11],115) * 10.2.2逼近 用简单的函数p(x)近似地代替函数f (x)，是计算数学中最基本的概念和方法之一。近似代替又称为逼近，函数f (x)称为被逼近的函数，p (x)称为逼近函数，两者之差 称为逼近的误差或余项。 * 1、一致逼近 以函数f (x)和p (x)的最大误差 作为度量误差f (x) - p (x)的“大小”的标准，在这种意义下的函数逼近称为一致逼近或均匀逼近，讲得更具体一点，也即对于任意给定的一个小正数? 0，如果存在函数p (x)，使不等式 成立，则称该函数p (x)在区间[a, b]上一致逼近或均匀逼近于函数f (x)。 * 2、平方逼近： 如果采用 作为度量误差 的“大小”的标准，在这种意义下的函数逼近称为平方逼近或均方逼近。这种方法要比一致逼近的相应问题简单得多。 * 10.3非线性方程求根 10.3.1 二分法 二分法是求方程f(x)=0根的近似值的常用方法。其基本思想是利用有根区间的判别方法确定方程根的区间[a,b] ,将有根区间平分为二；再利用有根区间的判别方法判断那一个区间是有根区间；重复上述步骤，直到小区间端点差的绝对值小于等于精度要求的数值，则用将上一区间的分半值作为方程的根的近似值。 * 根据精度要求可以事先计算出需执行步骤数n。 * 二分法的计算步骤如下： （1）计算区间端点的函数值f(a) , f(b)（不妨设f(a)＜0,f(b)＞0）；确定初始有根区间[a,b]。 （2）二分有根区间[a,b]，并计算 （3）判断：若 ，则方程的根为 ； 若 ，则有根区间为 若 ，则有根区间为 （4）如果│b-a│ε（ε为误差限），则方程的根为 ；否则转向步骤2，继续二分有根区间[a1,b1]，并计算中点值，继续有根区间的判断，直到满足精度要求为止，即│bn-an│ε。 n确定二分次数： * 【案例10.4】 用区间二分法求方程 f(x)=x3-5x+3在某区间内实根的近似值（精确到0.001） 解 ∵f(1.8)=-0.168＜0, f(1.9)=0.359＞0 ∴f(x)在区间[1.8,1.9]内有一个根。 由公式 取n=6， 计算结果列表如下： * 则方程在区间[1.8，1.9]内所求近似值为x*≈ x = 1.8328125 a=fzero(‘x^3-5*x+3’,[1:2]) fplot(‘[x^3-5*x+3,0]’,[1,2]) grid * 10.3.2 切线法（牛顿法） 1.切线法的基本思想：假设方程f(x)=0在区间[a,b]上有唯一根x*,过曲线y= f(x)上的一点（x0,f(x0)），作曲线的切线，用此切线与x轴的交点的横坐标x1作为方程的根x*的新的近似值, 再过点（x1,f(x1)），作曲线的切线，则又得到新的近似值,按此方法进行迭代计算，直到满足精度要求为止。 切线法（牛顿法）的迭代公式为 * 牛顿法有明显的几何意义。方程f (x) = 0的根x*在几何上表示曲线y = f (x)与x轴的交点。当求得x*的近似值xk以后，过曲线y= f (x) 上对应点（xk, f (xk)）作f (x) 的切线，其切线方程为 求此切线方程和x轴的交点，即得x*的新的近似值xk+1必须满足方程 这就是牛顿法的迭代公式 的计算结果。 由图可知，只要初值取得合适，点列{xk}就会很快收敛于x*。 * 2.切线法（牛顿法）的收敛性： 定理10.1. 设f(x) 在[a,b]上存在二阶连续导数，且满足条件 （1） f(a)f(b)0； （2） f’(x)在[a,b]上不等于零； （3） f”(x)在[a,b]上不变号； 则对任意初值x0∈[a,b] ， 只要满足 f(x0)f”(x0)=0 。 则由切线法迭代公式得到的近似根序列 {xn}平方收敛于方程f(x)=0 在区间[a,b]的唯一根x*。 牛顿迭代法的主要优点是收敛速度较快，它是平方收敛的。  * 3.切线法（牛顿法）的计算步骤 先判断有根区间[a,b]，然后选择初始值x0（一般地，若f”(x)0 ，则选择区间的右端点；若 f”(x)0，则选择区间的左端点），再建立迭代