用牛顿法求非线性方程.docVIP

实验 一、实验目标 1算法 2. 会编写计算机程序实现给定迭代函数的迭代算法及其加速；掌握迭代算法精度控制 二、实验问题 方程三、实验要求 有个实根. 将方程以下面六种不同，计算 , (b) , (c) , (d) (e) , (f) . 2. 对每一种迭代格式，编制一个程序进行运算，每种格式的敛散情况；用事后误差估计来控制迭代次数，并且出迭代的次数； 附录一： 《数值分析》实验报告（模板） 【实验课题】 【实验目标】 明确实验目标【理论概述与算法描述】 ’(xk)(x-xk), 于是方程可表示为 f(xk)+f’(xk)(x-xk)=0, 这是个线性方程，记其根为x(k+1), 则x(k+1)=xk-f(xk)/f’(xk),这就是牛顿迭代法求根. 埃特金加速收敛方法 设是根的某个近似值，用迭代一次得， 而由微分中值定理，有 其中介于和之间。 假设改变不大，近似地取某个近似值L，则有 若将校正值再迭代一次，又得由于 将它与前面的式子联立，消去未知的L，有 由此推知 ， 记 称为埃特金加速方法。 斯特芬森迭代法 将埃特金加速技巧与不动点迭代结合，则可得到如下的迭代法 即为斯特芬森迭代法 【实验问题】 方程有个实根. 将方程以下面六种不同，计算 , (b) , (c) , (d) (e) , (f) . 3. 对每一种迭代格式，编制一个程序进行运算，每种格式的敛散情况；用事后误差估计来控制迭代次数，并且出迭代的次数；【实验过程与结果】 【结果分析、讨论与结论】 从上述的运算结果可以看出，迭代公式1、2、4不收敛，3虽然收敛，但与其他迭代法的结果差异太大，对5和6分别用埃特金加速和斯特芬森迭代得到结果如下： 对于5埃特金加速结果： B = 2.0000 2.2804 2.2791 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 0 斯特芬森迭代结果： x = 2.0000 2.1547 2.2792 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 2.2790 对于6埃特金加速结果： B = 2.0000 2.2878 2.2790 2.2790 2.2790 2.2790 斯特芬森迭代结果： x = 2.0000 2.1544 2.2838 2.2790 　　　2.2790 　　　 　　　从以上结果可以看出，埃特金加速方法和斯特芬森迭代法确实可以加快收敛速度，且在此题的情况下，两种方法的加速效果差不多，但埃特金加速方法较斯特芬森迭代法来说更为简单易理解，运算步骤也少一些，因此对于此题，我们可以选用埃特金加速方法。 【附程序】 function x=Aitken(A); n=length(A); x=zeros(n,1); t=0; x(1)=A(1); for i=1:n-2 x(i+1)=A(i)-((A(i+1)-A(i))^2)/(A(i)-2*A(i+1)+A(i+2)); end function x=Steffsen(A,B) n=length(B); x=zeros(n,1); x(1)=B(1); for i=2:n x(i)=A(i)-((B(i-1)-A(i))^2)/(B(i)-2*B(i-1)+A(i)); end ④ %构造迭代算法x=(3*x+5)/