南京邮电大学-数值计算实践报告.pdfVIP

数值计算实践 I、方程求根 一、实验目的 熟悉和掌握 Newton 法,割线法,抛物线法的方法思路，并能够在 matlab 上编 程实现 二、问题描述 (1).给定一个三次方程,分别用 Newton 法,割线法,抛物线法求解. 方程的构造方法: (a)根:方程的根为学号的后三位乘以倒数第二位加 1 再除以 1000. 假设你的学号为 则根为 141*(4+1)/1000=0.564 (b)方程:以你的学号的后三位数分别作为方程的三次项,二次项,一次项的系数,根 据所给的根以及三个系数确定常数项. 例如: 3 2 你的学号是 则你的方程是 x +4x +x+a =0 的形式. 0 方程的根为 0.564,因此有 3 2 0.564 +4*0.564 +0.564+a0=0,于是 a0=-2.015790144 3 2 你的方程为 x +4x +x-2.015790144=0. 2 (2)假设方程是 sinx+4x +x+a0=0 的形式（三个系数分别是学号中的数字），重新 解决类似的问题 (3)构造一个五次方程完成上面的工作. 四次方程的构造:将三次多项式再乘以(x-p*)2 得到对应的五次多项式(p*为已经确 定的方程的根,显然,得到的五次方程有重根). (4)将（2 ）中的方程同样乘以(x-p*)得到一个新的方程来求解 注:(1)Newton 法取 0.5 为初值,割线法以 0,1 为初值,抛物线法以 0,0.5,1 为初值, (2)计算精度尽量地取高. 终止准则：根据| p n  p n1 |  来终止 (3)可供研究的问题：  （一） 的取值不同对收敛速度有多大的影响 （二）将注（1）中的初值该为其它的初值，对收敛性以及收敛速度有无影响 （三）能否求出方程的所有的根 （4 ）实验报告的撰写 实验报告包含的内容：（一）实验目的（二）问题描述（三）算法介绍（包 括基本原理）（四）程序（五）计算结果（六）结果分析（七）心得体会 三、算法介绍 在本问题中，我们用到了 newton 法，割线法，抛物线法。 1.Newton 法迭代格式为： f (x ) k x  x  k 1 k f (x ) k 当初值与真解足够靠近，newton 迭代法收敛，对于单根，newton 收敛速度 很快，对于重根，收敛较慢。 2.割线法：为了回避导数值的计算，使用上的差商代替，得到割线法迭代 公式： f (x )(x  x ) k k k 1 x  x  k 1 k f (x )  f (x ) k k 1 割线法的收敛阶虽然低于 newton 法，但迭代以此只需计算一次函数值，不 需计算其导数，所以效率高，实际问题中经常应用。 3.抛物线法：可以通过三点做一条抛物线，产生迭代序列的方法称为抛物 线法。其迭代公式为： p (x)  f [x ]  f [x , x ](x  x )  f [x , x , x ](x  x )(x  x ) k k k 1 k k k 1 k 2 k k 1 其中 f [xk , xk 1],f [xk , xk 1, xk 2 ] 是一阶均差和二阶均差。 收敛速度比割线法更接近于 newton 法。 对于本问题的解决就以上述理论为依据。终止准则为：| xn  xn1 |  本题中所有精