中国石油大学计算方法模拟试题

1、1. 计算方法实际计算时，由于受计算机字长限制而导致的误差称为舍入误差C2. x *=1.1021是经过四舍五入得到的近似数，有 _5_位有效数字，相对误差限为0.5*10 -4 c3. 利用二分法求方程1-x-sinx=0 在0,1内的根要二分 15 次。(0.5*10-4)4. 写出用Newton法建立求b的迭代公式 k+1=(xk2+b)/2x。5. 使用矩阵分解法求解线性方程组时，平方根法适用于系数矩阵为对称正定矩阵的方程组，追赶法适用于系数矩阵为三对角阵的方程组C6. 设线性方程组 Ax=b,为 7x1 10x2 1，则 |A| 2= 14.933 ,Cond (A)为5x1 7x2

2、0.7289，若右端向量有扰动b=(0.01,-0.01) T，则解的相对误差限为 2.89。7. 求解数值积分的Simpson公式的代数精度为：3 ，若将积分区间n等分，步长为h,贝U复化Simpson公式的截断误差为h的几阶无穷小，即O(h ? 4 )8. 应用龙贝格求积公式求积分，其整个计算过程的特点是：将积分区间逐次分半，并将每一公式先后两次的计算结果按一定线性组合构成新的精度较高近似值。9. 常微分方程初值问题的数值解法分为单步和多步，显式和隐式，下列方法属于哪一类？龙格-库塔法：单步、显式，阿当姆斯内插公式:多步、隐式10.若 s(x)=x3 x2(0 x 1)2x3bx2cx 1

3、(1，是以0,1,2为节点的三次样条函数，则x 2)得分b= -2，c= 324分，每题6分)1. 看书上或课件定义2.对于方程组142 fx1ry52-310x2二9521 、x3.8试构造一收敛的高斯-赛德尔迭代格式，并说明收敛理由解：将方程组变换为:f-,521x1f8-142x2二52-310 M3 t 9系数矩阵为严格对角占优阵，则方程组存在收敛的高斯-赛德尔迭代格式。把方程组等价变形为:小，必有:(a bxiyi)0，x1 2x2 1x3 8555115x2x1x3 -424139x3x1x251010收敛的高斯-赛德尔迭代格式为:x1(k 1Zx2(k)1 x3(k)8555x2

4、(k 1丄屮11x3(k)5424x3(k 11x1(k 1-x2(k1)9510103.以线性拟合为例简述最小二乘原理。答：设近似函数为y=a+bx，nR= (abXiyi)2。根据极值理论，要使R达到最i 1nn2 (a bx yJXii 1由方程组可以解出a，b的值，从而得到拟和曲线的表达式4.确定下列求积公式的常数a,使其代数精度尽量高，并判定其具有的代数精度解：当f(x)=1时:h0 f(x)dxh1dx h02f(h) ahf(O)f(h)hxdx0!h22h ah211h0 h2h 2 ,x dx!h3h0h2 ah202h，解得:a=1/12032h 3x dx丄h4h0h3丄

5、 h203h204212h 4.x dx!h5h0h4 h204h305212当 f(x)=x 时：当 f(x)=x2 时:当 f(x)=x 3 时:当 f(x)=x 4 时:说明所求求积公式具有三次代数精度得分三、证明题（16分，每题8分）fx0, x1,., xkf(Xi)Xi i)( XiXi i).(XiXk)k=1 o1.若 f(x)=(x-x o)(x-x 1)(x-x n)，Xi 互异，证明当 k=n+1 时 fx o,xi,.,x证明：由差商性质：当k=n时f (Xi)(XiXo).(XiXi i)(Xi Xi i).(Xi Xk)(XjXo).(xXii)(xXi)(x Xi

6、i).(xXk)(XiXo).(XiXi i)(XiXi i).(XiXk)当k=n+i时fx0, xi,., xkf(xi)Xi i)( XiXn)Xi i).(XiXn i)(Xn iX).(Xn i(Xn iX).(Xn iXn)n2. 证明对于牛顿-科特斯求积公式的科特斯系数有Ckn）ik 0证明：由牛顿-科特斯求积公式:f(x)dx(bna)Ckn)f(Xk)k 0b f (n1)() a (n 1)!(x)dx所以:b1dx b aan(n)Ckf(Xk)得分k 0四、计算题(26 分)设 f(x)=1 则(x)dx=O(n 1)!n(b a)Ckn)f(xk)，即:k 01. (

7、10分)给出sinx在0.4,0.7的数值表x0.40.50.60.7si nx0.389 420.479 430.564 640.644 22如果使用二次插值求sin0.63891的近似值，问如何选取结点，才使其近似值的误差较小？并求该近似值，小数点后保留 5位数字。(注意：拉格朗日插值与牛顿插值两种方法任选，若采用牛顿插值，构造出差商表)解：应选三个节点，使截断误差|R2(x)|v=|f (3)( )|(x-x0)(x-x1)(x-x2)| 尽量小故最靠近0.63891的三个节点一定满足要求显然，取0.5,0.6,0.7(1)采用拉格朗日插值:(x .6)(x 50.47943 (x .5

8、)(x。刀0.56464l_2(X)=(0.5 0.6)(0.5 0.7)(0.6 0.5)(0.6 0.7)0.64422(x 0.5)(x 0.6)(0.7 0.5)(0.7 0.6)=(x 0.6)(x 0.7)(x 0.5)(x0.7)0.479 430.564 640.02 0.01(x 0.5)( x 0.6)0.644 220.02所以：sin0.63891L 2(0.63891)(0.638910.6)(0.638910.7)047943 638910.5)(0.638910.7) 0 564 640.02 . 0.01 .(0.638910.5)(0.638910.6)0.6

9、44 220.020.479 4313891( .06109)0.564 640.010.03891 ( 0.06109) =0020.13891 0.038910.020.644 22.00238 0.479 43 0490.02 0.010.564 640.005400.020.644 22=0.05705 0.47938 0.17394=0.59627(2)采用牛顿差值:xiyi一阶差商二阶差商0.50.479 430.60.564 640.85210.70.644 220.7958-0.2815N(x) = 0.479 43 +0.8521(x-0.5) - 0.2815(x-0.5)

10、(x-0.6)所以 sin0.638912(0.63891)=0.479 43 + 0.8521*(0.63891-0.5) - 0.2815*(0.63891-0.5)*( 0.63891-0.6)=0.479 43 + 0.8521*0.13891 - 0.2815*0.13891*0.03891=0.479 43 + 0.11837 - 0.00152=0.596281 82. (8分)设max|f(|v= - ,x 2,8，用复化梯形公式计算2 f(x)dx的近似值8时，为使截断误差的绝对值不超过 -10 5，至少应将2,8分为多少等份？ 解：用复化梯形公式，截断误差：3Rn(f) =

11、b a h2f() 也2i1 f()1212 n21因为 max|f(x)|=-所以| Rn|=671所以至少分为671等份。3. (8分)用欧拉预报-校正法求初值问题步长h=0.3，小数点后保留5位数字。y 1 y在x=0.3，0.6处的数值解,y(0) 0解：由预报-校正公式有:0yn 1yn 12yn h(1 yn)h 20 2h=0.3, n=0,1,2,yn尹 y； 1(y0J2利用上述公式，及y(0)=0 得:y0 0.3* (1 0)0.30 32y10*101 0.32 0.31352y(0.3)y1=0.3138y；0.3135 0.3* (1 0.31352)0.64298

12、0.322y10.3135*1 0.31352 1 0.6429820.690262y(0.6)y2=0.69026一.填空k1已知 =3.1415926若其近似值的绝对误差限为 0.5 X 10-5,则该近似值是什么？ 2、对于充分接近90度的x,为不损失有效数字，应对公式1- sin(x)做何变化?3、 对于不动点迭代Xk+1 = O (X k),若在不动点x*满足 (x*)工0,则该迭代格式是几阶收敛的4、 牛顿迭代法的特点是什么？ 对于单根，它是几阶收敛的？ 5、关于线形方程组系数矩阵的条件数a、反映绝对误差放大倍数b、反映相对误差放大倍数c、条件数越大，方程组越呈“良”态6、写出两种

13、非线形方程的解法 7、追赶法适合解系数矩阵为的方程组8设xi (i=0, 1,2,3 , 4)为互异结点，li(x)为对应的插值基函数4贝U:xi3li (x) = i 04(xi2 4xi2)li(x) = 9、什么是三次样条插值函数？，写出三个要点10、A= 1 a ,当 a=，A 可做 LL T分解，VJ ，1 2r、其中L的元素满足 L=11、 向量 X= (x1,x2,x3 ) T ,贝U | x1+2x2|+| x1+x3|是不是一种向量范数?二.解答：1、 当A有扰动S A和b有扰动S b时，如何用矩阵A的条件数去估计方程组的相对误差| S x| / |x|?2.写出gauss列

14、主元的算法描述三、解方程组已知方程组Ax=b ,其中A=1 2 b=1尸10 )的经验公式，使它能和下表数据 相拟合：xi1.001.251.501.752.00yi5.105.796.537.458.46已知对数表x5.105.796.537.458.46lnx1.631.761.882.012.12五、已知函数表:x1246y0311231构造差商表，写出Newton插值多项式2、写出Laglanre插值多项式3、写出该插值多项式的余项六、设 f (x) =g(x)h(x)证明:f x0 , x1 = g(x0) hx0 , x1 + g x0 ,x1 h(x1)七、用最小二乘法解矛盾方程组2x + 3y = 6x + y = 22x + y = 42. 补充Newton迭代的大范围收敛性定理，并完成所给问题(8分)(1) Newton迭代收敛性定理如下：设f(x)在区间a, b上二阶导数存在，且对于 x a, b满足:则Newton迭代法收敛于f(x)=0在a, b上的唯一根。(2) 说明该定理每个条件的作用(3) 图示Newton迭代法的几何意义(4) 推导用Newton迭代法求正数a的平方根的迭代格式2.补充Newton迭代的大范围收敛性定理，并完成所给问题(8分)(1) Newton 迭代收敛性定理如