10数值分析A评分标准.docVIP

常州大学2010~2011学年第 2 学期硕士生考试试题评分标准 课程名称： 数值分析 试题类型：A （考生不得在此试题上填写答案，如填写一律无效） 本试卷考试时间为90分钟，满分100分。 （10分）当时, 试比较与在理论上与在数值计算上的差异？并叙述常见的防止误差的一些原则。 解：当时，两个表达式在理论上恒等， 但其数值计算结果不同，前者会出现相近数相减，失去有效数位，降低计算结果精度的问题；后者避免了相近数相减的问题，尽可能地保证了计算结果的精度。 ………… 5 分 防止误差的几个基本原则主要有： 防止大数“吃”小数； 避免除数绝对值远远小于被除数绝对值的除法； 避免相近数相减； 避免使用不稳定的算法； 5) 注意简化计算步骤，减少运算次数； ………… 5 分 （15分）叙述Newton插值法的方法思想，并根据列表函数 -1 0 1 2 2 -1 1 2 利用Newton插值方法求的插值逼近多项式。 解：Newton插值法的方法思想：对于给定节点，Newton插值方法通过构造如下形式的插值多项式 保证当新增加节点时，新的插值多项式满足 ， 其中系数为由确定的阶差商，仅与有关。这样就可以新增节点后的新的插值多项式能够在已有的插值多项式的基础上用较少的工作量得到。 ………… 7 分 Newton差商表： D= 2.0000 -1.0000 -3.0000 1.0000 2.0000 2.5000 2.0000 1.0000 -0.5000 -1.0000 ………… 5 分 ………… 3 分 （15分）写出变步长的梯形积分公式，并用复化辛普森（Simpson）公式计算 ，其中。 解：复化梯形求积公式之间的关系如下： （3-1） 其中。（3-1）式为变步长的梯形积分公式。 ………… 5 分 在变步长的梯形积分公式的基础上有 ； ………… 5 分 应用上述公式计算上述积分可得结果如下： 1.7271 =1.8006 1.7822 1.8001 1.7956 由此可得1.8001 ………… 5 分 （10分）设, ，试用LU分解法解线性方程组。 解： ………… 5 分 解可以转化为解，先解，再解可得所求解。具体如下：解可得；解可得。 ………… 5 分 （15分）对下述线性方程组 给出一个收敛的Gauss-seidel迭代格式，并用此收敛的Gauss-seidel迭代格式求此方程组的一个误差不超过的近似解。 解：改写上述方程组可得 ………… 5 分 新方程组的系数矩阵为严格对角占优矩阵，相应的Gauss-seidel迭代为 ………… 5 分 从初始近似解出发，可得满足精度要求的近似解为 ；迭代次数为 3。精确解为。 ………… 5 分 （15分）写出Newton迭代方法求非线性方程根的方法思想；并选用适当的方法求方程在附近的一个根，要求误差不超过。 解：求非线性方程根的Newton方法的方法思想： Newton法是一种迭代法, 在已知方程的第次近似解时, Newton法是利用曲线在处的切线来近似曲线, 把切线与x轴的交点作为与轴的交点的第次近似, 由此得到方程根的第次近似。即有 ………… 8 分 算法优点：对于方程的单根，方法具有二阶收敛速度；缺点：局部收敛性。 ………… 2 分 由Newton法解方程在附近的一个根可得迭代结果如下： -0.4239， -0.4263， 其中 -0.4263满足精度要求。 ………… 5 分 (10分)写出解常微分方程初值问题的四阶Runge－Kutta方法，并用四阶Runge－Kutta方法解方程 所确定的函数在处的近似值，取。 解：四阶龙格-库塔（Runge－Kutta）方法如下: ………… 5 分 利用Runge-kutta方法可得函数在的近似值为：y1= 2.0347，其中 k1= 2.0000 k2= 3.0125 k3= 3.4304 k4= 5.8074 y1= 2.0347。 ………… 5 分 （10分）已知列表函数 1 2 3 4 5