數(shù)值分析試卷及其答案4

1、1、 (本題 5分 取 99的 6位有效數(shù)字 94987. 9, 問以下這種算法有幾位有效數(shù)字05013. 094987. 9109910=-解:令99=x , 94987. 9*=x則5*1021 (-=x x x e (2分由于( 10(*x e x e -故5*1021 ( 10(-x e x e另一方面05013. 0991-故在這里 2-=m , 由 51-=+-n m 有 4=n . (3分 即算式至少有 4位有效數(shù)字 .2、(本題 6分用列主元 Gauss 消去法解線性方程組 .-=+=+-=-+3344531213332321321x x x x x x x x解:-=3344

2、747470453312334131134533123344533121311312r r r B -1649164900334045331247474703340453312 (4分故等價方程組為:=-=+=+-16491649334453312332321x x x x x x (1分同代得13-=x , 02=x , 43=x (1分3、(本題 6分已知 -=6134A ,求 1A , A, 2A .解:96314max 1=+-+=A (1分76134max =+-+=A(1分-=-=4518181761346314A A T18 45(17(451818172=-=-=-A A E

3、T即0441622=+-(3分 解得 2311+=, 2312-=, 231 (+=A A T 231 (2+=A A AT (1分4、(本題 7分給定線性方程組-=-71420328112315321x x x (1 試分別寫出 Jacobi 迭代格式和 Gauss-Seidel 迭代格式;(2 分析 Gauss-Seidel 迭代格式的收斂性 . 解:(1 Jacobi迭代格式為:+-=-=-+=+20/ 327( 1/( 81(15/ 234( (3 (1 1(3(3 (1 1(2(3 (2 1(1k k k k k k k k k x x x x x x x x x (2 分Gauss

4、-Seidel 迭代格式:+-=-=-+=+20/ 327( 1/( 81(15/ 234( 1(3 1(1 1(3(3 1(1 1(2(3 (2 1(1k k k k k k k k k x x x x x x x x x (2分(2Gauss-Seidel迭代格式的迭代矩陣 G 的特征方程為02032823=-0 48418300(2=+-解得127. 0600483004 418(418126. 1600483004 418(418023221-=>-+=則 1 (2>=G故 Gauss-Seidel 迭代格式發(fā)散 . (3分 5、(本題 8分用下列方法求 013 (3=-=

5、x x x f 在 20=x 附近的根,根的準(zhǔn) 確值 87938524. 1*=x ,要求計算結(jié)果準(zhǔn)確到四位有效數(shù)字 . (1用牛頓法;(2用弦截法,取 20=x , 9. 11=x 解:(1 33 (2-='x x f 牛頓法的迭代公式為=-+=-=+2 1(3123313023231x x x x x x x x k k k k k k k計算得888889. 11=x , 879452. 12=x3*21021-<-xx故 879. 1*=x (4分 (2弦截法的迭代公式為=-+=-=-+-+9. 1, 21(31( ( (102212111x x x x x x x x

6、f x f x x x f x x k a k k k k k k k k k k k 計算得881094. 12=x 879411. 12=x 3*21021-<-xx故 879. 1*=x (4分 6、(本題 8分給定數(shù)據(jù)如下 (1 寫出 (x f 的 3次 Lagrange 插值多項式 (3x L (2 寫出 (x f 的 3次 Newton 插值多項式 (3x N 解:(1由題設(shè)條件有00=x1 (0=x f21=x 3 (1-=x f 32=x 4 (2-=x f53=x2 (3=x f由于 n 次 Lagrange 插值多項式的基函數(shù)為.(.( .(.(1101110n k

7、k k k k k n k k k x x x x x x x x x x x x x x x x x x l -=+-+-故三次 Lagrange 插值多項式的基函數(shù)為5(3(2(301 50(30(20( 5(3(2( ( (3020103210-=-=-=x x x x x x x x x x x x x x x x x x x l5(3(61 52(32(02( 5(3(0( ( (3121013201-=-=-=x x x x x x x x x x x x x x x x x x x l5(2(61 53(23(03( 5(2(0( ( (3212023102-=-=-=x x x

8、 x x x x x x x x x x x x x x x x l3(2(301 35(25(05( 3(2(0( ( (2313032103-=-=-=x x x x x x x x x x x x x x x x x x x l (3分故所求三次 Lagrange 插值多項式3(2(151 5(2(32 5(3(21 5(3(2(3013(2(30125(2( 61( 4( 5(3(63( 5(3(2(301(1 ( ( ( ( ( ( ( ( (332211003-+-+-=-+-+-+-=+=x x x x x x x x x x x x x x x x x x x x x x x

9、x x l x f x l x f x l x f x l x f x L (1分(2由題中所給數(shù)據(jù),構(gòu)造下列差商表 (3分 由于(, , , (, , (, ( (21032101021001003x x x x x x x x x x f x x x x x x x f x x x x f x f x N -+-+-+=故所求三次 Newton 插值多項式3(2(51 2(31213(2(0(51 2(0(31 0(21 (3-+-+-=-+-+-=x x x x x x x x x x x x x N (1分7、(本題 8分設(shè) xa x f -=1 (,且 n x x x a ,., ,

10、 , 10互不相同,證明=-=kj jk xa x x x f 010(1,., , n k , . . . , 1, 0=并寫出 (x f 的 n 次 Newton 插值多項式 . 證:用數(shù)學(xué)歸納法來證明 當(dāng) 1=k 時(1(111( (, 01011001010110110x a x a x a x a x a x a x x x x x a x a x x x f x f x x f -=-+-=-=-=即當(dāng) 1=k 時公式成立 . (2分假設(shè)當(dāng) n m k <=時等式成立 即=-=mi im xa x x x f 010(1,., , +=+-=1110(1,., , m i i

11、 m x a x x x f那么當(dāng) 1+=m k 時+=+=+=+=+-=-=-=-=1101001010011012101110(1( ( (1 (1 (11,., , ,., , 1 ,., , (m i im i i m m mi i m i i m m m m m xa x a x a x a x x x a x a x x x x x f x x x f x x x x x f即公式對 1+=m k 亦成立有歸納法原則知原等式對任意 n k 均成立 (4分 我們以 n x x x ,., , 10為插值節(jié)點來求 n 次 Newton 插值多項式 因為.(,., , . (, , (,

12、 ( (110101021001003-+-+-+=n n x x x x x x x x x f x x x x x x x f x x x x f x f x N故所求插值多項式為(1(103x xa x N k nk ki i-=-=其中.( (101k k x x x x x x x -=- (4分8、 (本題 5分求滿足條件 的艾爾米特差值多項式 .解:令 10=x , 21=x ,代入艾爾米特差值多項式121112101102000201003 ( ( (1(21 ( ( (1 (21 (y x l x x y x l x x x x y x l x x y x l x x x x

13、 x H '-+-+'-+-= (2分這里 1010 (x x x x x l -=, 0101 (x x x x x l -=,得5982 (233+-+-=x x x x H(3分9、(本題 6分求函數(shù) x y arctan =在 0, 1上的一次最佳平方逼近多項式 . 解:設(shè) 1 (0=x , x x = (1, , 10span =,所求函數(shù)為x a a x *1*0*1 (+=,則214arctan , (2ln 214arctan , (31 , (21 , (11 , (11101211110100-=-=dx x x y dx x y dx x dx x dx(

14、3分由正規(guī)方程組-=+-=+21431212ln 21421*1*0*1*a a a a (1分 解得042909. 022ln 23*0-=a791831. 02ln 3623*1+-=ax x 791831. 0042909. 0 (*1+= (2分 10、 (本題 9分 運用梯形公式、 辛普森公式、 柯特斯公式分別計算積分 dx e x 10,并估計各種方法的誤差(要求小數(shù)點后至少保留 5位 . 解:運用梯形公式:8591409. 121101=+e e dx e x(2分其誤差2265235. 012101(1213=-=e e f R (1分運用辛普森公式:7188612. 1461

15、12101=+e e e dx e x(2分其誤差00094385. 0288028801=-=e ef R (1分運用柯特斯公式:718282688. 17321232790114321411=+e e e e e dx e x(2分其誤差000001404. 0494524945 (266=-=e a b a b f R (1分11、 (本題 6分已知 (x f 的函數(shù)值如下: 用復(fù)合梯形公式和復(fù)合辛普森公式求 dx x f 6. 28. 1 (的近似值 .解:用復(fù)合梯形公式,小區(qū)間數(shù) 4=n ,步長2. 0 8. 16. 2(41=-=h則99. 40. 10 0. 80. 64. 4(

16、21. 322. 04=+=T (3分復(fù)合辛普森公式,小區(qū)間數(shù) 2=n ,步長4. 0 8. 16. 2(21=-=h則98. 40. 10 0. 84. 4(40. 621. 364. 02=+=S (3分12、 (本題 8分用 3, 2=n 高斯 -勒讓德公式計算積分 dx x ex31sin .解:由于高斯求積公式為=-nk k kx f Adx x f 011( (其中 k x 是 (1x P n +的零點 首先將積分區(qū)間轉(zhuǎn)化為 1, 1-令 2+=t x 則 3, 1x 時 1, 1-t (1分 而dt t edx x e I t x-+=112312sin(sin (2分令 2si

17、n( (2+=+t e t g t2=n 時7745967. 01-=t 02=t 7745967. 03=t 5555556. 01=A 8888889. 02=A 5555556. 03=A (2分3=n 時8611363. 01-=t 3399810. 02-=t 3399810. 03=t 8611363. 04=t 3478548. 01=A 6521452. 02=A 6521452. 03=A 3478548. 04=A(2分I A1 g (t1 + A2 g (t 2 + A3 g (t 3 + A4 g (t 4 = 10.95014051 （1 分） y' = x

18、+ y 13、 （本題 6 分） 用改進歐拉法求解 (0 x 1 ，h = 0.2 ， 取兩位小數(shù)。 y ( 0 = 1 1 y n+1 = y n + (k1 + k 2 2 改進歐拉法格式為 k1 = hf ( x n , y n k = hf ( x + h, y + k n n 1 2 解 n = 0,1,2, L （2 分） 其中 f ( x, y = x + y, y 0 = 1, h = 0.2, n = 0,1,2,3,4 代入上式得： n xn yn 1 0.2 1.24 2 0.4 1.58 3 0.6 2.04 4 0.8 2.64 5 1.0 3.42 （4 分） 14

19、、 （本題 6 分）寫出用四階經(jīng)典的龍格庫塔方法求解下列初值問題的計算公 式： y' = x + y 1 ,0 < x <1 y ( 0 = 1 3y y' = 2 1+ x ,0 < x <1 y ( 0 = 1 解：令 h = 0.2 k1 = k 2 = 1 k 3 = k 4 = f ( xn , y n = xn + y n h h h h f ( x n + , y n + k1 = x n + + y n + k1 = 1.1( x n + y n + 0.1 2 2 2 2 h h h h f ( x n + , y n + k 2 = x n + + y n + k 2 = 1.11( x n + y n + 0.11 2 2 2 2 f ( x n + h, y n + hk 3 = x n + h + y n + hk 3 = 1.222( x n + y n + 0.222 h (k1 + 2k 2 + 2k 3 + k 4 = 0.2214 x n + 1.2214 y n + 0.0214 6 y n +1 = y n + （3 分） k1 = 3 y n (1 + x n k = 3( y + 0.1k (1 + x + 0.1 n 1 n 2 2