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

1.计算方法实际计算时，由于受计算机字长限制而导致的误差称为 舍入误差 。 2.x*=1.1021是经过四舍五入得到的近似数，有 5 位有效数字，相对误差限为 0.5*10-4 。 3.利用二分法求方程1-x-sinx=0在[0,1]内的根要二分 15 次。（??0.5*10-4） 4.写出用Newton法建立求的迭代公式 xk+1=(xk2+b)/2xk 。 5.使用矩阵分解法求解线性方程组时，平方根法适用于 系数矩阵为对称正定矩阵的方程组 ，追赶法适用于 系数矩阵为三对角阵的方程组 。 6.设线性方程组Ax=b，为，则||A||2= 14.933 ,Cond?(A)为 289 ，若右端向量有扰动?b=(0.01,-0.01)T，则解的相对误差限为 2.89 。 7.求解数值积分的Simpson公式的代数精度为： 3 ，若将积分区间n等分，步长为h, 则复化Simpson公式的截断误差为h的几阶无穷小,即 O(h ？ 4 ) 8.应用龙贝格求积公式求积分，其整个计算过程的特点是：将积分区间逐次分半，并将每一公式先后两次的计算结果按一定线性组合构成新的精度较高近似值。 9.常微分方程初值问题的数值解法分为单步和多步，显式和隐式，下列方法属于哪一类？龙格-库塔法： 单步、显式 ，阿当姆斯内插公式： 多步、隐式 。 10.若s(x)=，是以0,1,2为节点的三次样条函数，则b= -2 ，c= 3 。 得分 得分 二、解答题（24分，每题6分） 1.看书上或课件定义 2.-1 4 22 -3 10 -1 4 2 2 -3 10 5 2 1 x1 x2 x3 5 9 8 试构造一收敛的高斯-赛德尔迭代格式，并说明收敛理由。 5 2 1-1 4 2 5 2 1 -1 4 2 2 -3 10 x1 x2 x3 8 5 9 系数矩阵为严格对角占优阵，则方程组存在收敛的高斯-赛德尔迭代格式。把方程组等价变形为： 收敛的高斯-赛德尔迭代格式为： 3.以线性拟合为例简述最小二乘原理。 答：设近似函数为y=a+bx，R=。根据极值理论，要使R达到最小，必有：， 由方程组可以解出a，b的值，从而得到拟和曲线的表达式。 4.确定下列求积公式的常数a，使其代数精度尽量高，并判定其具有的代数精度。 解：当f(x)=1时： 当f(x)=x时： 当f(x)=x2时：，解得：a=1/12 当f(x)=x3时： 当f(x)=x4时： 说明所求求积公式具有三次代数精度。 得分三、证明题（16分，每题8分） 得分 1.若f(x)=(x-x0)(x-x1)...(x-xn)，xi互异，证明当k=n+1时f[x0,x1,...,xk]=1。 证明：由差商性质： 当k=n时= 当k=n+1时 = 0 + 2.证明对于牛顿-科特斯求积公式的科特斯系数有。 证明：由牛顿-科特斯求积公式： 设f(x)=1则=0。 所以：，即： 得分 得分 四、计算题（26分） 1.（10分）给出sinx在[0.4,0.7]的数值表 x 0.4 0.5 0.6 0.7 sinx 0.389 42 0.479 43 0.564 64 0.644 22 如果使用二次插值求sin0.63891的近似值，问如何选取结点，才使其近似值的误差较小？并求该近似值，小数点后保留5位数字。（注意：拉格朗日插值与牛顿插值两种方法任选，若采用牛顿插值，构造出差商表） 解：应选三个节点，使截断误差|R2(x)|=|f(3)(?)||(x-x0)(x-x1)(x-x2)|尽量小。 故最靠近0.63891的三个节点一定满足要求。 显然，取{0.5,0.6,0.7}。 （1）采用拉格朗日插值： L2(x) = = 所以：sin0.63891 L2(0.63891) = = = = =0.59627 （2）采用牛顿差值： xi yi 一阶差商 二阶差商 0.5 0.479 43 0.6 0.564 64 0.8521 0.7 0.644 22 0.7958 -0.2815 N2(x) = 0.479 43 + 0.8521(x-0.5) - 0.2815(x-0.5)(x-0.6) 所以sin0.63891 N2(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