常州大学数值分析12-13研究生试卷A参考答案.docVIP

一．（1）已知函数用秦九昭法计； （2）秦九昭法 解：（1） ………… 5 分 秦九昭法………… 5 分 A）防止大数“吃”小数； B）避免除数绝对值远远小于被除数绝对值的除法；C）避免相近数相减；D）避免使用不稳定的算法；E）注意简化计算步骤，减少运算次数；………… 5 分 二． （1）以分量形式写出解此线性方程组的Jacobi迭代格式和Gauss-Seidel迭代格式； （2）求和； （3）判断Gauss-Seidel迭代格式的敛散性。 解：（1）Jacobi迭代 ， Gauss-Seidel迭代 ， ………… 5 分 （2），；………… 5 分 （3）因为方程组系数矩阵严格对角占优，所以Gauss-Seidel迭代格式收敛。………… 5 分 三． 方程上存在唯一实根； （2）叙述牛顿法求方程根的方法思想； （3）以初值要求 。 解：（1）因为， 并且由于，所以方程在区间上存在唯一实根。 (2) 牛顿法求方程根的方法思想：求根的Newton方法是迭代方法, 对于给定的第k次初始近似解, Newton迭代法的思想是用曲线在点的切线与x轴的交点的横坐标作为曲线与x轴交点，即方程根的第k+1次近似。相应的迭代公式如下: , 当时可以认为为满足精度的近似解。………… 7 分 (3) 以初值X = 1.000000000000000 0.914155281832543 0.910017665783406 0.910007572548888 ………… 3 分 四．（1）求，使得数值求积公式具有尽可能高的代数精度，并出其代数精度计算。 解：（1）代数精度为3，其中。………… 5 分 （2） ………… 5 分 五．插值插值（2）对下述列表函数： 12.0 14.0 15.0 16.0 y 27.0 22.0 13.0 34.0 写出差商表；并写出逼近上述列表函数的三次Newton插值多项式，Newton插值方法通过构造如下形式的插值多项式 保证当新增加节点时，新的插值多项式满足 ， 其中系数为由确定的阶差商，仅与有关。这样就可以新增节点后的新的插值多项式能够在已有的插值多项式的基础上用较少的工作量得到。 Lagrange插值 ，其中 Newton插值插值 12.0000 27.0000 14.0000 22.0000 -2.5000 15.0000 13.0000 -4.6667 -2.1667 16.0000 34.0000 1.7500 2.1250 4.2917 ………… 4 分 ………… 3 分 六． 解：（1）L = 1.0000 0 0 0.3333 1.0000 0 0.6667 -1.0000 1.0000 U = 3.0000 -1.0000 1.0000 0 -1.6667 3.6667 0 0 -2.0000 ………… 5 分 （2）方程组的解为 x=[1, -2, -1]’. 七．，完成利用改进的Euler方法解下述初值问题的Matlab程序, 使得输出结果yy(n+1)为下述微分方程初值问题的解函数在处函数值的近似值。 function ex( ) a=0; b=0.3; y0= 2 ; n= 3 ; [x,yy]=improved_euler_method(@fxy,a,b,y0,n) return function z=fxy(x,y) z= 2*x+3 ; return function [x,y]=improved_euler_method(f,a,b,y0,n) h=(b-a)/n; y=zeros(1,n+1); x=a:h:b; y(1)=y0; for j=1:n temp= f(x(j),y(j)) ; y(j+1)=y(j)+h*temp; y(j+1)=y(j)+h*(temp+f(x(j+1),y(j+1)))/2; end return (2) 叙