计算方法考试试题及答案.docVIP

2006 级《计算方法》试题 B 卷答案 一、 填空，每题 4 分，共 34 分 1） a 的绝对误差界为 1 10 2 ， a 的相对误差界为 1 10 4 ； 2 4 2 2）法方程组为： 2 8 a 1 2 2 ； b 1 8 24 2 1 2 2 0 , 0 , 0 , 1 , , , 0 , sin x 1 , 1 , sin x x sin xdx 1 , 1 1 0 2 8 3 2 3）设 A 17 ， cond A 17 × 17=289 ； 7 10 7 10 5 7 5 7 1 1 4）应改写为 x 16 x 0 x 8 x 1 16 x 17 x 18 x 14 x 13 x 1 5）均差 f [0,1, 2] 0 ， l 0 x 2 x 1 x 2 ； 6）此数值求积公式的代数精度为： 3 ； 7）求解 u u n t n 1 e 1 h ； u t e 1 的隐式 Euler ： u n 1 1 h 8）用二分法进行一步后根所在区间为 : [ 1, 2]。 9) LL T 分解为： 1 0 1 2 ； 2 1 0 1 10) [0,1] 上以 ( x ) ln 1 权函数的正交多项式 0 x 1 ， 1 x x 1 。 x 4 11) ； xk 1 1 xk e xk 0, 1, 2, xk ex k 1 , k 1 2 2 1 12) 正交矩阵 H 2 1 2 ： 3 1 2 2 二、计算题 1．（ 15 分）解： 已知， 2 1, 1 0, 0 1, 20, 1 1 , 0 3 。 2 2 c 0 1 1 0, c1 2 1 3 0, c2 1 1 3 ，故此为二步一阶方法。 2 2 2 2 0 2 2 局部误差主项为： 3 h 2u tn O h 3 。 2 2 又 1 ，满足根条件，故此差分格式收敛。 又考虑模型问题 u u 则，有特征多项式： ( ) h ( ) 2 h 1 3 h ，其中 h h 2 0 2 由判别式可知 1 的充要条件是： h 3 h 2 ，而 h 3h 自然成立， 2 2 2 2 则由 3 h 2 得出 h 4 , 0 。 2 3 由于 0 h h 40 h 4 ，故 h 的取值范围是： 0 h 1 。﹟ 3 30 2．（15 分）解： m 1 2 m x m 0,1, 2, ，则 x dx , x 2 1 x 3dx 0 ， x 4 dx 2 ， 3 x 5dx 0 。 0 d x 2 ， 1 2 1 1 1 1 1 3 1 1 5 1 2 0 1 3 4 4 4 1 1 2 x 2 ，令 2 x 0 即得， x 2 2x 0 x 15 25 2 x 0 5 5 3 5 2 0 x 2 5 得 Gauss点：x0 ,1 3 1 2 3 3 f x 1 ，x ，令 A 0 f 。取 x f x dx 5 A1 f 5 1 5 即得到方程组： 2 A 0 A1 ， 0 A1 A0 ，解之，得 A0 A1 1 ，从而得到 3 3 2 1 3 3 4 1 f f ，又取 f x x x f x dx 3 5 5 1 2 4 4 2 1 1 5 5 1 5 25 x 4 x 4 dx 1 9 3 7 7 3 7 147 故所得到的数值求积公式是具有 m=3 次代数精度 Gauss求积公式。 2 2 1 3 3 2 3 1 2 5 e e e 5 ﹟ 5 1 3 3 3．（ 10 分）解： ATA 1 0 1 1 1 1 1 0 0 0 1 ， 1 det ( I - A H A ) 1 1 2 1 2 0 1 1 1 则 A H A 的特征值为 1 2, 2 0 1 2 , 2 0 , 所以 Σ 2 。 下面求对应的标准正交的特征向量 (正规直交 )，即 1 1 1 ; 1 , v 2 1 1 p1 1 , v1 1 p1 1 2 ; 2 1 1 1 1 即 V v1 v2 2 2 , 因 rank(A)=1，故有 V1 2 。 计算得 1 1 1 2 2 2 1 U 1 AV 1Σ 1 1 1 2 1 1 = 0 2 == u1， 0 1 0 2 得约化的奇异值分解 A U 1ΣV 1H 1 1 1 1 1 = 2 0 0 0 2 2 0 计算 u2, 使其与 U1 构成 R2 的一组标准正交基，可取 U 2 =u2= , 1 则 U U 1 U 2 是酉阵 , 故矩阵 A 的奇异值分解（满的奇异值分解）为 1 1 1 1 1 0 2 0 2 2 ﹟ 0 0 0 1 0 0 1 1 2 2 k 1 k x1 1 a x 2 4