研一计算方法实验报告.doc

数值实验1.1 迭代法的设计和运行 实验目的 迭代法是解各种方程(含方程组)的基本方法, 它常通过构造一个定常的迭代格式，重复计算而产生一个收敛的解的序列（标量序列，向量序列等），逐步逼近问题的真解。对同一个问题，常可设计不同的迭代格式，这些格式的计算效果可能相差很大，对初始值的选择的要求不同，收敛速度也不同。 实验题 求方程x3+4x2-10=0在[1,1.5]内的根, 设计下列3个迭代格式： k=0,1,… (1)xk+1=10+xk-4xk2-xk3 (2)xk+1=1210-xk3（原迭代格式有误，根据方程修改得之。） (3)xk+1=10xk+4 ? 实验步骤和要求 1.选取若干个可行的初始值, 分别用3 种迭代格式进行计算, 记录每步的计算结果, 并在坐标纸上画出, 分析其变化趋势,判别其收敛性； 2.观察和分析迭代法结果关于初始值的依赖关系； 3.你还能设计出新的迭代格式吗？ 4.用Steffen加速法尝试下。 实验内容 本次实验采用Matlab编程。 1．因求方程在区间[1,1.5]内的根，不妨选取初始值分别为1,1.25,1.5，并对以上三个迭代格式分别进行迭代运算，序列前后项误差小于12×10-4即认为该迭代格式下序列收敛。迭代序列计算结果和分析如下： （1）初始值X0=1时，各迭代格式的迭代序列图如图1所示： 图1 初始值X0=1时各迭代格式下的迭代序列图 初始值X0=1时，各迭代格式的迭代计算结果及分析如下： 迭代格式（1），发散。迭代次数n=7。第七个数达到10^205数量级，导致前面的数约等于0。 1.0e+205 * Columns 1 through 8 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -2.1086 NaN Columns 9 through 16 NaN NaN NaN NaN NaN NaN NaN NaN Columns 17 through 20 NaN NaN NaN NaN 迭代格式（2），收敛。迭代次数n=15 Columns 1 through 8 1.0000 1.5000 1.2870 1.4025 1.3455 1.3752 1.3601 1.3678 Columns 9 through 16 1.3639 1.3659 1.3649 1.3654 1.3651 1.3653 1.3652 1.3652 Columns 17 through 20 1.3652 1.3652 1.3652 1.3652 迭代格式（3），收敛。迭代次数n=7 Columns 1 through 8 1.0000 1.4142 1.3590 1.3660 1.3651 1.3652 1.3652 1.3652 Columns 9 through 16 1.3652 1.3652 1.3652 1.3652 1.3652 1.3652 1.3652 1.3652 Columns 17 through 20 1.3652 1.3652 1.3652 1.3652 小结：从以上计算结果以及迭代序列图可以很明显的发现在初始值X0=1时迭代格式（1）发散，迭代格式（2）和（3）均收敛，且（3）比（2）收敛速度快。 （2）初始值X0=1.25时，各迭代格式的迭代序列图如图2所示： 图2 初始值X0=1.25时各迭代格式下的迭代序列图 迭代格式（1），发散。迭代次数n=7。第七个数达到10^138数量级，导致前面的数约等于0。 1.0e+138 * Columns 1 through 8 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -2.0514 Inf Columns 9 through 16 NaN NaN NaN NaN NaN NaN NaN NaN Columns 17 through 20 NaN NaN NaN NaN 迭代格式（2），收敛。迭代次数n=13 Columns 1 through 8 1.2500 1.4184 1.3367 1.3795