昆工工研上机题.doc

昆明理工大学工科研究生 《数值分析》上机实验课题 理学院 信息与计算科学 教研室 任课教师： 结合课程教学，配备适当的上机实验以便加深课堂教学的实践性，同时通过实验可以加强学生对数学模型的总体分析，算法选取，程序结构，上机调试和结果分析等环节的训练,为使实验环节有成效,需要写出课题计算的实验报告,并以此作为《数值分析》课程期末成绩评定的一部分。 实验报告内容要求： 课题名称 班级、姓名、学号 三、目的和意义 方法的理论意义和实用价值。 四、计算公式 五、结构程序设计 六、结果讨论和分析 如初值对结果的影响；不同方法的比较；该方法的特点和改进；整个实验过程中（包括程序编写，上机调试等）出现的问题及其处理等广泛的问题，以此扩大知识面和对实验环节的认识。 实验课题数量要求： 　　请从以下12个课题中选择1个不同的课题，题中带﹡号的要求为选做。 课题一 迭代格式的比较 问题提出 设方程f(x)=x- 3x –1=0 有三个实根 x=1.8793 , x=-0.34727 ,x=-1.53209现采用下面六种不同计算格式，求 f(x)=0的根 x 或x x = x = 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，松弛因子的选取，对计算结果的影响。 课题四 矩阵求逆与矩阵行列式 一、问题提出 应有列主元Gauss-Jordan消去法求满秩矩阵的逆矩阵，并计算A的行列式det的值，如下列矩阵之一： 二、要求 分析列主元Gauss-Jordan消去法的计算公式； 确定选主元，换行，计算行（约化非主元行）和交换列序等四个子程序； ﹡应用结构程序设计编出计算n阶非奇异方阵的通用程序﹡； ﹡计算考核题 结果： ﹡ 三、目的和意义 通过该课题的实验，主要掌握Gauss-Jordan消去法求非奇异矩阵的逆阵的程序设计方法； 体会Gauss-Jordan消去法求解线性方程组的关键步骤； 提高科学计算和编程的能力。 课题五 函数插值方法 一、问题提出 对于给定的一元函数 的n+1个节点值 。试用L