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

1.计算方法实际计算时，由于受计算机字长限制而导致的误差称为舍入误差。2.x*=1.1021是经过四舍五入得到的近似数，有5位有效数字，相对误差限为0.5*10-4。3.利用二分法求方程1-x-sinx=0在[0,1]内的根要二分15次。（0.5*10-4）4.写出用Newton法建立求b的迭代公式xk+1=(xk2+b)/2xk。5.使用矩阵分解法求解线性方程组时，平方根法适用于系数矩阵为对称正定矩阵的方程组，追赶法适用于系数矩阵为三对角阵的方程组。6.设线性方程组Ax=b，为7.07511072121xxxx，则||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)=)21(12)10(2323xcxbxxxxx，是以0,1,2为节点的三次样条函数，则b=-2，c=3。二、解答题（24分，每题6分）1.看书上或课件定义2.对于方程组试构造一收敛的高斯-赛德尔迭代格式，并说明收敛理由。-1422-310521x1x2x3598得分解：将方程组变换为：系数矩阵为严格对角占优阵，则方程组存在收敛的高斯-赛德尔迭代格式。把方程组等价变形为：10921031513453211412583512521xxxxxxxxx收敛的高斯-赛德尔迭代格式为：10921031513453211412583512521)1()1()1()()1()1()()()1(kkkkkkkkkxxxxxxxxx3.以线性拟合为例简述最小二乘原理。答：设近似函数为y=a+bx，R=niiiybxa12)(。根据极值理论，要使R达到最小，必有：0)(21niiiybxaaR，0)(21iniiixybxabR由方程组可以解出a，b的值，从而得到拟和曲线的表达式。4.确定下列求积公式的常数a，使其代数精度尽量高，并判定其具有的代数精度。hhffahhffhdxxf02)](')0('[)]()0([2)(解：当f(x)=1时：hhhdx0]11[21521-1422-310x1x2x3859当f(x)=x时：]0[2]11[]0[221202hhahhhhxdxh当f(x)=x2时：]20[]0[23120232hahhhhdxxh，解得：a=1/12当f(x)=x3时：]30[121]0[241220343hhhhhdxxh当f(x)=x4时：]40[121]0[251320454hhhhhdxxh说明所求求积公式具有三次代数精度。三、证明题（16分，每题8分）1.若f(x)=(x-x0)(x-x1)...(x-xn)，xi互异，证明当k=n+1时f[x0,x1,...,xk]=1。证明：由差商性质：kikiiiiiiixxxxxxxxxf0110))...()()...(()(xk]x1,...,f[x0,当k=n时))...()()...(()(110kiiiiiiixxxxxxxxxf=0))...()()...(())...()()()...((110110kiiiiiikiiiiiiiixxxxxxxxxxxxxxxxxx当k=n+1时101110))...()()...(()(xk]x1,...,f[x0,niniiiiiixxxxxxxxxif=0+1))...(())...((101101nnnnnnxxxxxxxx2.证明对于牛顿-科特斯求积公式的科特斯系数有10)(nknkC。证明：由牛顿-科特斯求积公式：nkbanknkbadxxnfxfCabdxxf0)1()()()1()()()()(！得分设f(x)=1则bandxxnf)()1()()1(！=0。所以：nkknkbaxfCababdx0)()()(1，即：1)(0)(nkknkxfC四、计算题（26分）1.（10分）给出sinx在[0.4,0.7]的数值表x0.40.50.60.7sinx0.389420.479430.564640.64422如果使用二次插值求sin0.63891的近似值，问如何选取结点，才使其近似值的误差较小？并求该近似值，小数点后保留5位数字。（注意：拉格朗日插值与牛顿插值两种方法任选，若采用牛顿插值，构造出差商表）解：应选三个节点，使截断误差|R2(x)|=|f(3)()||(x-x0)(x-x1)(x-x2)|尽量小。故最靠近0.63891的三个节点一定满足要求。显然，取{0.5,0.6,0.7}。（1）采用拉格朗日插值：L2(x)=220.644)6.07.0)(5.07.0()6.0)(5.0(640.564)7.06.0)(5.06.0()7.0)(5.0(430.479)7.05.0)(6.05.0()7.0)(6.0(xxxxxx=220.64402.0)6.0)(5.0(640.56401.0)7.0)(5.0(430.47902.0)7.0)(6.0(xxxxxx所以：sin0.63891L2(0.63891)得分=220.64402.0)6.00.63891)(5.00.63891(640.56401.0)7.00.63891)(5.00.63891(430.47902.0)7.00.63891)(6.00.63891(=220.64402.00.038910.13891640.56401.0)06109.0(0.13891430.47902.0)06109.0(0.03891=220.64402.0054000.640.56401.000849.0430.47902.000238.0=17394.047938.005705.0=0.59627（2）采用牛顿差值：xiyi一阶差商二阶差商0.50.479430.60.564640.85210.70.644220.7958-0.2815N2(x)=0.47943+0.8521(x-0.5)-0.2815(x-0.5)(x-0.6)所以sin0.63891N2(0.63891)=0.47943+0.8521*(0.63891-0.5)-0.2815*(0.63891-0.5)*(0.63891-0.6)=0.47943+0.8521*0.13891-0.2815*0.13891*0.03891=0.47943+0.11837-0.00152=0.596282.（8分）设max|f’’(x)|=81,x[2,8]，用复化梯形公式计算82)(dxxf的近似值时，为使截断误差的绝对值不超过51021，至少应将[2,8]分为多少等份？解：用复化梯形公式，截断误差：Rn(f)=)(''112)28()(''12232fnfhab因为max|f’’(x)|=81所以|Rn(f)|=5231021811126n解得n=671所以至少分为671等份。3.（8分）用欧拉预报-校正法求初值问题0)0(1'2yyy在x=0.3，0.6处的数值解，步长h=0.3，小数点后保留5位数字。解：由预报-校正公式有：])(11[2)1(20121201nnnnnnnyyhyyyhyyh=0.3,n=0,1,2,...利用上述公式，及y(0)=0得：3135.0]3.0101[*23.003.0)01(*3.002101yyy(0.3)y1=0.313869026.0]64298.013135.01[*23.03135.064298.0)3135.01(*3.03135.0221202yyy(0.6)y2=0.69026一.填空k1．已知=3.1415926…若其近似值的绝对误差限为0.5×10-5,则该近似值是什么？2、对于充分接近90度的x,为不损失有效数字，应对公式1-sin(x)做何变化？3、对于不动点迭代Xk+1=Φ(Xk),若在不动点x*满足Φ‘(x*)≠0，则该迭代格式是几阶收敛的4、牛顿迭代法的特点是什么？对于单根，它是几阶收敛的？5、关于线形方程组系数矩阵的条件数a、反映绝对误差放大倍数b、反映相对误差放大倍数c、条件数越大，方程组越呈“良”态6、写出两种非线形方程的解法7、追赶法适合解系数矩阵为的方程组8、设xi(i=0，1,2,3，4)为互异结点，li(x)为对应的插值基函数则：403)(ixlixi=402)()24(ixlixixi=9、什么是三次样条插值函数？，写出三个要点10、A=1a,当a=，A可做ＬＬT分解，12其中Ｌ的元素满足Ｌ＝11、向量X=（x1,x2,x3）T,则|x1+2x2|+|x1+x3|是不是一种向量范数？二．解答：1、当A有扰动δA和b有扰动δb时，如何用矩阵A的条件数去估计方程组的相对误差||δx||/||x||?2．写出gauss列主元的算法描述三、解方程组已知方程组Ax=b,其中A=12b=10.312(1)写出解此方程组的Jacobi迭代公式，讨论用Jacobi迭代解此方程组的收敛性(2)写出解此方程组的Gauss-Seidel迭代公式，讨论用Gauss-Seidel迭代解此方程组的收敛性四．求形如y=aebx(a,b为常数，且a0)的经验公式，使它能和下表数据相拟合：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)h[x0,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迭代收敛性定理如下：设f(x)在区间[a,b]上二阶导数存在，且对于x∈[a,b]满足：①f(a)f(b)0;②f′(x)≠0;③f〞(x)存在且不变号；④取x0∈(a,b),使f〞(x0)f(x0)0（2分）则Newton迭代法收敛于f(x)=0在[a,b]上的唯一根。（2）说明该定理每个条件的作用①保证有根②保证单根③保证凸凹性不变④-保证收敛（2分