数值分析实验报告.docxVIP

数值分析实验报告 （第二章） 实验题目： 分别用二分法、牛顿迭代法、割线法、史蒂芬森迭代法求方程 f 的根x=1，观察不同初始值下的收敛性，并给出结论。 问题分析： 题目有以下几点要求： 不同的迭代法计算根，并比较收敛性。 选定不同的初始值，比较收敛性。 实验原理： 各个迭代法简述 二分法：取有根区间[a,b]的重点x0，确定新的有根区间[a1,b1]的区间长度仅为[a,b] 牛顿迭代法：不动点迭代法的一种特例，具有局部二次收敛的特性。迭代格式为 x 割线法：是牛顿法的改进，具有超线性收敛的特性，收敛阶为1.618. 迭代格式为 x 史蒂芬森迭代法：采用不动点迭代进行预估校正。至少是平方收敛的。迭代格式为 y z x 这里φx可采用牛顿迭代法的迭代函数 实验内容： 写出该问题的fx 代码如下： 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 分别写出各个迭代法的迭代函数 代码如下： 二分法： function y=dichotomie(a,b,e) i=2; m(1)=a; while abs(a-b)e t=(a+b)/2; s1=f(a); s2=f(b); s3=f(t); if