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

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

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

3、以0,1,2为节点的三次样条函数，则x 2)b=-2， c=3。得分24分，每题6分)1. 看书上或课件定义2.对于方程组-142x152-310x2二9-521 _.x3丿8试构造一收敛的高斯-赛德尔迭代格式，并说明收敛理由 解：将方程组变换为：521r rx1r、8-142x2二5-2-310 、x3.9丿系数矩阵为严格对角占优阵，则方程组存在收敛的高斯-赛德尔迭代格式。把方程组等价变形为:218x1 x2 - x3555c 1 o 5x2 x1 x3 -424139x3 -x1 x251010收敛的高斯-赛德尔迭代格式为:x1(k1)Zx2(k)1 x3(k)8555x2(k1)1(k

4、1)x11 x3(k)5424x3(k1)切1)?x2(k1)510_9103.以线性拟合为例简述最小二乘原理n答：设近似函数为y=a+bx ，R= (a bXj y2。根据极值理论，要使R达到最i 1小，必有:(a bXiyi)o，nyJXi2 (a bxii 1由方程组可以解出a，b的值，从而得到拟和曲线的表达式4.确定下列求积公式的常数a,使其代数精度尽量高，并判定其具有的代数精度h0 f(x)dx护（0）f(h)ah2 f(0)f(h)解：当f(x)=1时:当 f(x)=x 时：当 f(x)=x 2 时:当 f(x)=x 3 时:当 f(x)=x 4 时:1dxhxdx0h尹1!h2-

5、0 h ah211-0 h2 2h213hX dxh-0032h3 .1.4h sx dxh-0042h1h sx4dxh5-00522h2 ah20 2h，解得:3122h h 0 3h 12h4丄 h20 4h312a=1/12说明所求求积公式具有三次代数精度得分三、证明题（16分，每题8分）1. 若 f(x)=(x-x O)(X-X 1)(X-X n)，xi 互异，证明当 k= n+1时 fX0,X1,.,xk=1 o证明：由差商性质:fx0, x1,., xkkf(X)i 0 (XiX).(XiXi 1 )(XiXi1).(XiXk)当k=n 时f (Xi)(Xi Xo)(XiXi 1

6、)(Xi Xi 1)(XiXk)(xX。)(XiXi1)(XjXj)(Xj Xj 1)(Xi Xk)(XiX。)(Xj Xi (XjXi J(Xj Xk)当k=n+1 时fx0, x1,., xkf(xi)Xi 1 )(XiXi 1).(XiXn 1)(xn 1X).(Xn 1 Xn)2. 证明对于牛顿-科特斯求积公式的科特斯系数有Cj 1k 0证明：由牛顿-科特斯求积公式:f(x)dx(bna)Ckn)f(xk)k 0b f(n 1)() a (n 1)!(x)dx设 f(x)=i 则 bUa (x)dx=0 a (n 1)!所以：bn1dx b a (b a) ckn)f(Xk)，即:ak

7、 0nckn)f(Xk)1得分k 0四、计算题(26 分)1.( 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

8、。(1)采用拉格朗日插值：L2(X)=g脇 0-43(x 0.5)(x 0.7)(0.6 0.5)(0.6 0.7)0.56464& 0恥 .6)0.64422(0.7 0.5)(0.7 0.6)(x 0.6)(x0.7)0020.479 43(x 0.5)(x0.7)0010.564 64(x 0.5)(x 0.6)0.020.644 22所以：sin0.63891L2(0.63891)(0.638910.6)(0.638910.7)0.020.479 43(0.638910.5)(0.638910.7)0.010.564 64(0.638910.5)(0.638910.6)0.020.64

9、4 220.03891( 0.06109)= 0.02=0.13891 0.038910.020.479 430.644 22严)0.564 64000空 0.479 43 000竺!0.02 0.010.564 640.644 220.02=0.05705 0.47938 0.17394=0.59627(2 )采用牛顿差值:xiyi一阶差商二阶差商0.50.479 430.60.564 640.85210.70.644 220.7958-0.2815N2(x) = 0.479 43 + 0.8521(x-0.5) - 0.2815(x-0.5)(x-0.6)所以 sin0.638912(0.

10、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 (x)|= 8 ,x 2,8，用复化梯形公式计算2f(x)dx的近似值时，为使截断误差的绝对值不超过1 10 5，至少应将2,8分为多少等份?解：用复化梯形公式，截断误差:2h f()(82)31212nf()1 因为 max|f， (x

11、)|=-863112 n2所以 | Rn(f) |=671所以至少分为671等份。3.(8分)用欧拉预报-校正法求初值问题y 1y(0)2y在x=0.3，0.6处的数值解,0步长h=0.3，小数点后保留5位数字。解：由预报-校正公式有：y01yn h(1 yj)h 21 ( 0 )2h=0.3,nn，1,2,yn 1yn 尹 yn 1 (ym)利用上述公式，及y(0)=0 得:y；0 0.3* (1 0)0.3y10 0*1 0 1 0.320.31352y(0.3) y1 =0.3138y 0.3135 0.3* (1 0.31352)0.64298y10.3135 字* 1 0.31352

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

13、 (i=0，1,2,3，4)为互异结点，li(x)为对应的插值基函数44贝U:xi3li(x)= (xi2 4xi 2)li(x) = i 0i 09、什么是三次样条插值函数？，写出三个要点r10、A= (1 a丿 ，当a= ，A可做LL T分解,其中L的元素满足 L =11、向量 X= (x1,x2,x3 ) T ,贝U | x1+2x2|+| x1+x3|是不是一种向量范数?二.解答：1、 当A有扰动SA和b有扰动Sb时，如何用矩阵A的条件数去估计方程组的相对误差|冈| / |x|?2 写出gauss列主元的算法描述三、解方程组已知方程组Ax=b ,其中A=1 2b=10 )的经验公式，使

14、它能和下表数据相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) 图示Newt on迭代法的几何意义(4) 推导用Newt on迭代法求正数a的平方根的迭代格式2 补充Newt on迭代的大范围收敛性定理，并完成所给问题(8分)(1) Newton 迭代收敛性定理如下： 设f(x)