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

第五章 非线性方程的数值解法 §1 引言 §1 引言 (1) 画图法 例1 方程f(x)=x3-x-1=0 确定其有根区间 解：用搜索的方法，不难发现 f(0)0 f(2)0 在区间（0，2）内至少有一个实根 设从x=0出发,取h=0.5为步长向右进行根的 搜索,列表如下 §2 二分法 §2 二分法 将含方程根的区间平分为两个小区间，然后判断根在哪个小区间，舍去无根的区间，而把有根的区间再一分为二，再判断根属于哪个更小的区间，如此周而复始 ，直到求出满足精度要求的近似根。 §3 迭代法 §3 迭代法 二.迭代求解过程 求解过程分以下二步： (1) 建立迭代公式。 由公式f(x)=0出发将其分解为等价形式 x= φ(x) 式中φ(x)叫做方程的迭代函数 例如：f(x)=x3+2x2-4=0 可以分解为 x=x+f(x)=x3+2x2+x-4; x=x-(x3+2x2-4)/(3x2-4x) (２) 进行迭代计算。由初值x0出发，按迭代函数进行计算 xn+1=φ(xn) （n=0,1,2,…）称为迭代公式 x0,x1,x2,…　称为迭代序列，其值相应地称为根的零次、一次、二次……近似值 序列的计算过程称为迭代过程 (1)将方程改写为 x=(1+x)1/3 取x=1.5逐次迭代得 x0=1.5, x1=1.35721, x2=1.33086, x3=1.32588, x4=1.32494, x5=1.32476, x6=1.32473, x7=1.32472, x8=1.32472 最后取稳定至小数后５位的迭代值x8 =1.32472为方程的近似根 (2)如果将方程改写为 x=x3-1， 并建立迭代公式xk+1=xk3-1，则得到的结果是 x0=1.5, x1=2.375, x2=12.39, …… 迭代值愈来愈大，这时迭代序列的极限不存在，迭代法失效。 我们把迭代序列有极限存在的迭代过程称为收敛的迭代过程；否则就是发散的。 五.迭代序列的误差估计 可用来预估迭代次数， 迭代序列的误差估计 上述公式仅当L1时才成立，否则L≈1可能出现假收敛的情况 八、迭代法加速（埃特金方法）（1） （Steffense法） 例 ：用埃特肯法解x=e-x 解：取x0=0.5,按照埃特肯法公式得 x0= 0.5 y0= e-0.5 = 0.60653,z0= e-0.60653 = 0.54524 2．局部收敛性：若x* 是方程x= ?(x)的根?’(x)在x*的邻近连续，且 | ?’(x*) | 1，则迭代公式 局部收敛。 1、优点：牛顿迭代法具有平方收敛的速度，所以在迭代过程中只要迭代几次就会得到很精确的解。这是牛顿迭代法比简单迭代法优越的地方。 2、缺点：选定的初值要接近方程的解，否则有可能得不到收敛的结果。再者，牛顿迭代法计算量比较大，因每次迭代除计算函数值外还要计算导数值。 小结 （一）Newton迭代法的基本思想 （二）Newton迭代法的构造 （三）几何意义 （四）Newton迭代法的局部收敛性 （五） Newton迭代法的计算流程 （六）重根时牛顿迭代法的改进 作业 1、教材第159页第5,7题 §5 弦截法(正割法 )和抛物线法 5.1 弦截法(正割法 ) 5.2、抛物线法* 抛物线法计算公式 0-1004 0-1003 0-1002 -0-1001 -0-010 x[k]-x[k-1] d(f(xk)) f(xk) x[k] k please input initial x0=1.0 演示 k xk f(xk) 0 3.000000000000 -1.472366552741015 1 -1.156000022150 3.213549548141341 2 0.189437548249 0.726814111584673 3 0.714043350055 0.075722015564973 4 0.782542444682 0.001129739710533 5 0.783