数值分析试题A091.docxVIP

中国石油大学（北京）2017--2017 学年第一学期 研究生期末考试试题A (闭卷考试) 课程名称：数值分析 所有试题答案写在答题纸上，答案写在试卷上无效 所有试题答案写在答题纸上，答案写在试卷上无效 题号 题号 一 二 三 四 五 六 七 八 总分 得分 注：计算题取小数点后四位 一、填空题（共 30 分，每空 3 分） 1、已知x k (k ? 0,1, , n) 是互异节点，l k ?x?是对应节点的Lagrange 插值基函数， P(x) 是任意一个首项系数为 1 的n ? 1次多项式，则 P(x) ? ?n k ?0 ??x3 ? x2 , 0 ? x ? 1 2、设分段多项式 S (x) ? ? ??2x3 ? bx2 ? cx ?1, 1 ? x ? 2 P(x )l (x) = 。 k k 是以0,1,2 为节点的三次样条函数，则b ? ， c ? 。 3、如果 A 是正交矩阵，则Cond ( A) = 。 2 4、用 x = 3.141 作为? 的近似值，则 x 有 位有效数字，其绝对误差限为 。 5、数值积分公式? 3 0 f (x)dx ? 3 ? f (1)? f (2)?是否为插值型求积公式： ，其代数 2 精度为 。 6、下列matlab程序中s2计算的是 ， 并指明s1与s2的区别为 。 其中： aex ? a ?10x ; a, x ? R 。 t=0; s2=1e14; for i=1:1e6 temp= 1/(1e3+i); t=t+temp; s2=s2+temp; end s1= t+1e14; 二、（8 分）已知函数表 x x y 0 1 1 0 2 1 y? 1 试利用重节点Newton 差商构造满足插值条件P(0) ? 1, P(1)? 0, P (1)? 1, P(2) ? 1, 的三次多项式 P(x) 。（要求构造出差商表） 三、（8 分）已知向量 x ? (2,0,2,1) T ，试构造 Householder 变换阵，使 Hx ? (0,0, k,0) T ， 其中k ? R 。 1 ? ? 四、（12 分）已知勒让德（Legendre）正交多项式P 0 ? 1,P 1 ? x, P ? 2 2 3x2 ?1 ，试利用勒 让德正交多项式在二次多项式类 H ? span ?1, x2?中求一个多项式 S ?x?，使其成为 f ?x?? ex在??1，1?上的最佳平方逼近函数，并计算出平方误差。五、（10 分）写出求解线性代数方程组 ?? x ? 2 x ? ?1 2 ? 2 x ? 5 3 ? x ? 3 x ? ? 1 ?? 2 x1 ? 7 x2 ? 2 1 3 的 Gauss-Seidel 迭代格式，并分析此格式的敛散性。 六、（12 分）用追赶法求解三对角方程组。（要求写出 LU 分解的具体计算过程） ? 2 ?1 0 0 ?? x ? ? 1 ? ? ?? 1 ? ? ? ? ?1 2 ?1 0 ?? x2 ? ? ? 0 ? ? 0 ?1 2 ?1?? x ? ? 0 ? ? ?? 3 ? ? ? 4? 0 0 ?1 2 4 ?? x ? ? 0 ? 2 ? 2 ? 2 ?七、（12 分）给出计算 x 2 ? 2 ? 2 ? 的迭代格式，讨论迭代格式的收敛性, 八、（8 分）求解常微分方程初值问题 ? y? ? f (x, y) ? ?y(x ) ? y ? 0 0 的改进欧拉公式 y ? y ? h [ f (x , y ) ? f (x , y ? hf (x , y ))] 是几阶方法？ n?1 n 2 n n n?1 n n n 其中h ? x ? x 为常数， n ? 0,1, 。 n?1 n 中国石油大学（北京）2017--2017 学年第一学期 研究生期末考试试题标准答案A (闭卷考试) 课程名称：数值分析 题号 题号 一 二 三 四 五 六 七 八 总分 得分 一、（30 分） 1、?n i?0 (x ? x ) ； 2、b ? ?2, c ? 3 ； 3、1； i 4、3， 1 ?10?2 ； 5 、 是， 1； 2 1 1 1 6、计算S ? 1014 ? ? 的值， 103 ?1 103 ? 2 103 ?106 1 1 1 其中 s2 ? (((1014 ? ) ? ) ? ) ， 103 ?1 103 ? 2 103 ?106 1 1 1 s1 ? 1014 ? ( ? ) ，s2 会产生大数吃小数的问题。 103 ?1 103 ? 2 103 ?106 二、（8 分）构造差商表： i i 0 1 1 0 -1 1 0 1 2 2 1 1 0 -1 所以插值多项式 P(x)