(用牛顿运动定律解决问题3.docVIP

实验七 非线性方程求根 一、实验目标 1. 掌握常用的非线性方程求根算法（二分法、不动点迭代法与Newton法）及加速技术（Aitken加速与Steffsen加速）. 2. 会编写计算机程序实现给定迭代函数的迭代算法及其加速；掌握迭代算法的精度控制方法. 二、实验问题 求代数方程 的实根. 三、实验要求 1．方程有一个实根：. 将方程以下面六种不同方式等价地改写，构造迭代格式，计算 : (a) , (b) , (c) , (d) (e) , (f) . 2. 对每一种迭代格式，编制一个程序进行运算，观察每种格式的敛散情况；用事后误差估计来控制迭代次数，并且输出迭代的次数；观察不同初值的结果. 3. 从理论上分析各种格式的收敛性及收敛阶. 4. 将收敛较慢的一种格式分别用Atken 方法及 Steffsen 方法加速，通过输出结果了解加速效果. 5. 将一种不收敛的方法用Steffsen 方法加速得到收敛的迭代. 附录一： 《数值分析》实验报告（模板） 【实验课题】 用牛顿迭代法求非线性方程根 【实验目标】 明确实验目标 1. 掌握常用的非线性方程求根算法（二分法、不动点迭代法与Newton法）及加速技术（Aitken加速与Steffsen加速）. 2. 会编写计算机程序实现给定迭代函数的迭代算法及其加速；掌握迭代算法的精度控制方法. 3探索不同方式改写方程的收敛程度 【理论概述与算法描述】 牛顿法 设已知方程f(x)=0有近似根 xk,将函数f(x)在点xk展开，有 f(x)=f(xk)+f’(xk)(x-xk), 于是方程可表示为 f(xk)+f’(xk)(x-xk)=0, 这是个线性方程，记其根为x(k+1), 则x(k+1)=xk-f(xk)/f’(xk),这就是牛顿迭代法求根. 埃特金加速收敛方法 设是根的某个近似值，用迭代一次得， 而由微分中值定理，有 其中介于和之间。 假设改变不大，近似地取某个近似值L，则有 若将校正值再迭代一次，又得由于 将它与前面的式子联立，消去未知的L，有 由此推知 ， 记 称为埃特金加速方法。 斯特芬森迭代法 将埃特金加速技巧与不动点迭代结合，则可得到如下的迭代法 即为斯特芬森迭代法 【实验问题】 1.求代数方程 的实根. 2．方程有一个实根：. 将方程以下面六种不同方式等价地改写，构造迭代格式，计算 : (a) , (b) , (c) , (d) (e) , (f) . 3. 对每一种迭代格式，编制一个程序进行运算，观察每种格式的敛散情况；用事后误差估计来控制迭代次数，并且输出迭代的次数；观察不同初值的结果. 4. 从理论上分析各种格式的收敛性及收敛阶. 5. 将收敛较慢的一种格式分别用Atken 方法及 Steffsen 方法加速，通过输出结果了解加速效果. 6. 将一种不收敛的方法用Steffsen 方法加速得到收敛的迭代. 【实验过程与结果】 用matlab编程计算代数方程的根 分别编写6个迭代法编程，对结果进行分析 【结果分析、讨论与结论】 迭代公式1： x1 = 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 2.0000 1.5000 迭代公式2： x2 = 1.0e+142 * 0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -1.4947 -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf 迭代公式3： x3