数值分析实验报告62338.pdf

数 值 分 析 实 验 报 告 6 2 3 3 8 用心整理的精品 word 文档，下载即可编辑！！ 数值分析实验报告 （第二章） 实验题目： 分别用二分法、牛顿迭代法、割线法、史蒂芬森迭代法求方程 的根 ，观察不同初始值下的收敛性，并给出结论。 问题分析： 题目有以下几点要求： 1. 不同的迭代法计算根，并比较收敛性。 2. 选定不同的初始值，比较收敛性。 实验原理： 各个迭代法简述 二分法：取有根区间 的重点 ，确定新的有根区间 的区间长度 仅为 区间长度的一版。对压缩了的有根区间 重复以上过程，又得到 新的有根区间 ，其区间长度为 的一半，如此反复，……，可得一系 列有根区间，区间收敛到一个点即为根。 牛顿迭代法：不动点迭代法的一种特例，具有局部二次收敛的特性。迭代 格式为 割线法：是牛顿法的改进，具有超线性收敛的特性，收敛阶为 1.618. 迭代 格式为 精心整理，用心做精品 2 用心整理的精品 word 文档，下载即可编辑！！ 史蒂芬森迭代法：采用不动点迭代进行预估校正。至少是平方收敛的。迭 代格式为 这里 可采用牛顿迭代法的迭代函数。 实验内容： 1. 写出该问题的 函数 代码如下： function py= f(x) syms k ; y=(k^2+1)*(k-1)^5; yy=diff(y,k); py(1)=subs(y,k,x); py(2)=subs(yy,k,x); end 2. 分别写出各个迭代法的迭代函数 if s1(1)*s3(1)=0 b=t; 代码如下： else 二分法： a=t; function y=dichotomie(a,b,e) end i=2; m(i)=t; m(1)=a; i=i+1; end while abs(a-b)e t=(a+b)/2; y=[t,i+1,m]; s1=f(a); end s2=f(b); 牛顿迭代法： s3=f(t); 精心整理，用心做精品 3 用心整理的精品 word 文档，下载即可编辑！！ function x1=x2; y=NewtonIterative(x,e) x2=t; i=2;