09—10数值分析(研)试题及答案-东北大学.docxVIP

班 级 学 号 姓 名  ? 东北大学研究生院考试试卷 一（ 1－8） 一（ 9-10 ? 总分 ） 二 三 四 五 六 ? 2009 —2010 学年第 1 学期 ? ○ 课程名称： 数值分析 ? ? 一、解答下列各题：（每题 5 分，共 50 分） 5．求满足条件 f (0) 0, f (1) 1, f ( 2) 0, f (1) 0 的三次插值多项式 H 3 ( x) 的表 ? ? 1．设近似值 x 321.235 近似 x* 具有 5 位有效数字，求 x 的相对误差限。 达式。 密 解 设 H 3 (x) x( x 2)( ax b) ，则 (a b) 1， a 0 。 ? 解 0.5 10 2 0.000015564 ( 0.15564 10 4 0.0015564 %) r ? 321.235 于是， H 3 ( x) x(x 2) 。 ? 2 1 3 x1 5 ? b n n ○ 2．用 LU 分解法求方程组 4 5 1 x2 3 的解。 6. 设求积公式 Ak f ( xk ) 是插值型求积公式，求 Ak . f ( x)dx ? 6 3 2 x3 1 a k 0 k 0 ? 2 1 3 2 1 3 1 0 0 2 1 3 n n b b n b ? 解 Ak l k ( x) dx l k ( x)dx 1dx b a 。 ? 解 由于 4 5 1 2 3 5，所以A 2 1 0 0 3 5 k 0 k 0 a a k 0 a 封 6 3 2 3 0 7 30100 7 或，由于公式对 f ( x) 1 精确成立，所以 n b a 。 ? T T Ak 1dx b ? 解 Ly=b 得： y (5, 7, 14) ，再解 Ux=y 得 x ( 1,1,2) 。 k 0 a ? 3．解线性方程组的迭代格式 x (k Mx ( k) g , k 0,1,2,...是否收敛，为什么？ 7．求区间 [ － 1， 1] 上权函数为 (x) 2 的二次正交多项式 P2 ( x) 。 ? 1) x ○ 2 1 0 解 P0 ( x) 1， P1 ( x) x (x,1) x ， ? 其中 M 0 4 1 。 1 (1,1) ? 0 2 3 ? (x 2 ,1) (x 2 , x) 3 P2 ( x) x 2 1 x 2 解 不收敛。因为 2 是 M的特征值，所以 (M ) 2 1 . (1,1) x 5 ? ( x, x) ? 4．求简单迭代法 xk xk 1 , (k 0,1,2,...) 的收敛阶。 ? 1 8．设 f ( x) 5x3 x 2 3 ，求差商 f [ 0,1], f [ 7,6,3,5], f [ 3,1,2,6,4] 。 ? 2 xk ? 解 设 lim x k 得 1 ，即 2 。又由于 ( x) x 1 ，所以， 解 f [0,1] f (1) f (0) 4 ， f [7,6,3,5] 5 ， f [3,1,2,6,4] 0 。 ? k 2 2 x 1 0 ? ( ) 1 1 0 ， ( ) 2 0 ，所以，迭代法收敛阶为 2。 ? 2 2 2 2 1 ? ? ? ? ○ ? ? ? ? 密 ? ? ? ? ○ ? ? ? ? 封 ? ? ? ? ○ ? ? ? ? ? ? ? ? ? ? ?  9．给定离散数据 xi －1 0 1 2 yi 3 1 2 4 试求形如 y a bx 2 的拟合曲线。 解 基函数为 0 ( x) 1 ， 1 (x) x2 ，于是 0 (1, 1, 1, 1)T , 1 (1, 0, 1, 4)T , f (3, 1, 2, 4)T ， 正则方程组为： 4a 6b 10 3/ 2, b 2 / 3 6a 18b ，解之得： a 21 所以，拟合曲线为： y 3 2 x 。 2 3 10．求解初值问题 y yex 1 x 2 的改进 Euler 方法是否收敛？为 y(1) 2 什么？ 解 因为 f ( x, y) yex 关于变量 y 满足 Lipschitz 条件，故收敛。 2x1 x2 x3 1 0 二、（ 11 分）用 Jacobi 法解线性方程组 x1 3x2 x3 2 ，取 x( 0) 0 ， 2x1 x2 4x3 3 0 若使 x( k) x* 10 3 ，问应迭代多少步？ 1 0 1/ 2 1/ 2 5 . 解 由于 Jacobi 迭代矩阵为 B 1/ 3 0 1/ 3 , B 1 1/ 2 1/ 4 0 6 迭代一步得： x(1) (1/2, 2/ 3, 3/4)T ，若使 x (k) x* 10 3 ，则有： 1  k ln