[专业课]09数值分析A评分标准.docVIP

常州大学2009~2010学年第 2 学期硕士生考试试题评分标准 课程名称： 数值分析 试题类型：A （考生不得在此试题上填写答案，如填写一律无效） 本试卷考试时间为120分钟，满分100分。 （10分）当充分大时, 试比较与在理论上与在数值计算上的差异？并叙述常见的防止误差的一些原则。 解：当充分大时，两个表达式在理论上恒等， 但其数值计算结果不同，前者会出现相近数相减，失去有效数位，降低计算结果精度的问题；后者避免了相近数相减的问题，尽可能地保证了计算结果的精度。 ………… 5 分 防止误差的几个基本原则主要有： 防止大数“吃”小数； 避免除数绝对值远远小于被除数绝对值的除法； 避免相近数相减； 避免使用不稳定的算法； 注意简化计算步骤，减少运算次数； ………… 5 分 （15分）已知列表函数 0 1 2 3 1 -1 3 2 利用Newton插值方法求的插值逼近多项式，利用插值多项式近似计算。 解： Newton差商表： D = 1.0000 -1.0000 -2.0000 3.0000 4.0000 3.0000 2.0000 -1.0000 -2.5000 -1.8333 ………… 5 分 ………… 5 分 。 ………… 5 分 （10分）已知列表函数 -1 0 1 2 0.21 0.94 1.7 2.3 若用函数拟合上述列表函数，试用最小二乘法确定参数的大小。 解：写出正规方程组 ………… 5 分 解上述正规方程组得 ………… 5 分 （15分）写出龙贝格（Romberg）方法的数值积分公式，并用龙贝格方法计算 ，要求误差不超过。 解：龙贝格（Romberg）方法计算定积分的数值积分公式如下： ， 其中。 ………… 7 分 利用上述公式计算可得 romberg_table = 1.6599 1.6301 1.6319 1.6375 1.6318 1.6332 ………… 8 分 （10分）设, ，试用高斯消去法或LU分解法解线性方程组。 解：利用LU分解法可得 ………… 5 分 分别解可得方程组的解为。 ………… 5 分 （15分）写出解线性方程组的Jacobi迭代法方法及Gauss-Seidel迭代方法的分量形式；对下述线性方程组 给出一个收敛的Gauss-seidel迭代格式，并说明收敛的理由。 解：设为方程组的第k次迭代解，则解线性方程组的Jacobi迭代法方法及Gauss-Seidel迭代方法的分量形式分别为 ， 。 ………… 8 分 将方程组改写为， 则新方程组的系数矩阵严格对角占优， 从而解上述方程组的一个收敛的Guass-Seidel迭代格式 ， 对于任意收敛。 ………… 7 分 （15分）写出Newton迭代方法求非线性方程根的方法思想；并选用适当的方法求方程在附近的一个根，要求误差不超过。 解：解非线性方程根的Newton方法是迭代方法, 对于给定的第k次初始近似解, Newton迭代法的思想是用曲线在点的切线与x轴的交点的横坐标作为曲线与x轴交点，即方程根的第k+1次近似。相应的迭代公式如下: , 当时可以认为为满足精度的近似解 ………… 8 分 本题取, 由上述迭代公式, 容易得到: x1 = 4.4936 x2 = 4.4934 The root x is 4.4934 ………… 7 分 (10分)用四阶Runge－Kutta方法解方程 取。 解：写出Runge－Kutta方法公式： ………… 5 分 利用y0=1及上述公式可得： k1= -1.0000 k2= -0.5117 k3= -0.8104 k4= -0.1894 y1= 0.8083 k1= -0.3988 k2= -0.1747 k3= -0.3239 k4= -0.0183 y2= 0.7376 所得列表函数为： x = 1.0000 1.3000 1.6000 y = 1.0000 0.8083 0.7376 ………… 5 分 5