2008级计算机专业《数值分析》实验课题.docVIP

东北大学本科 《数值分析》实验课题 计算数学研究室 数学实验室 结合课程教学，配备适当的上机实验以便加深课堂教学的实践性，同时通过实验可以加强对数学模型的总体分析，算法选取，程序结构，上机调试和结果分析等环节的训练,为使实验更为有成效,需要写出课题计算的实验报告,以此可作为《数值分析》课程成绩评定的参考。 实验报告内容要求： 程序设计流程图 源程序代码 三、程序运行结果 四、程序运行结果分析 课题一 迭代格式的比较 问题提出 设方程f(x)=x- 3x –1=0 有三个实根 x=1.8793 , x=-0.34727 ,x=-1.53209现采用下面六种不同计算格式，求 f(x)=0的根 x 或x x = x = x = 二、要求 1、编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； 2、用事后误差估计来控制迭代次数，并且打印出迭代的次数； 3、初始值的选取对迭代收敛有何影响； 4、分析迭代收敛和发散的原因。 三、目的和意义 1、通过实验进一步了解方程求根的算法； 2、认识选择计算格式的重要性； 3、掌握迭代算法和精度控制； 4、明确迭代收敛性与初值选取的关系。 课题二 线性方程组的直接算法 一、问题提出 给出下列几个不同类型的线性方程组，请用适当算法计算其解。 1、设线性方程组 = x= ( 1, -1, 0, 1, 2, 0, 3, 1, -1, 2 ) 2、设对称正定阵系数阵线方程组 = x = ( 1, -1, 0, 2, 1, -1, 0, 2 ) 3、三对角形线性方程组 = x= ( 2, 1, -3, 0, 1, -2, 3, 0, 1, -1 ) 二、要求 对上述三个方程组分别利用Gauss顺序消去法与Gauss列主元消去法；平方根法与改进平方根法；追赶法求解（选择其一）； 应用结构程序设计编出通用程序； 比较计算结果，分析数值解误差的原因； 尽可能利用相应模块输出系数矩阵的三角分解式。 三、目的和意义 1、通过该课题的实验，体会模块化结构程序设计方法的优点； 2、运用所学的计算方法，解决各类线性方程组的直接算法； 3、提高分析和解决问题的能力，做到学以致用； 4、通过三对角形线性方程组的解法，体会稀疏线性方程组解法的特点。 课题三 线性方程组的迭代法 一、问题提出 对课题二所列目的和意义的线性方程组，试分别选用Jacobi 迭代法，Gauss-Seidol迭代法和SOR方法计算其解。 二、要求 1、体会迭代法求解线性方程组，并能与消去法做以比较； 2、分别对不同精度要求，如由迭代次数体会该迭代法的收敛快慢； 3、对方程组2，3使用SOR方法时，选取松弛因子=0.8，0.9，1，1.1，1.2等，试看对算法收敛性的影响，并能找出你所选用的松弛因子的最佳者； 4、给出各种算法的设计程序和计算结果。 三、目的和意义 1、通过上机计算体会迭代法求解线性方程组的特点，并能和消去法比较； 2、运用所学的迭代法算法，解决各类线性方程组，编出算法程序； 3、体会上机计算时，终止步骤 或k （予给的迭代次数），对迭代法敛散性的意义； 4、体会初始解 x，松弛因子的选取，对计算结果的影响。 课题四 数值积分 一、问题提出 选用复合梯形公式，复合Simpson公式，Romberg算法，计算 I = I = I = I = 二、要求 编制数值积分算法的程序； 分别用两种算法计算同一个积分，并比较其结果； 分别取不同步长，试比较计算结果（如n = 10, 20等）； 给定精度要求，试用变步长算法，确定最佳步长。 三、目的和意义 深刻认识数值积分法的意义； 明确数值积分精度与步长的关系； 根据定积分的计算方法，可以考虑二重积分的计算问题。 1 数值分析实验课题 1