《數(shù)值計算方法》試題集和答案解析

1、WORD格式整理數(shù)值計算方法復(fù)習試題、填空題:一4-101|A =A = 14-11、 了 一1 4，則A的LU分解為 一答案:一4-10115/4-1-4/15 1_-56/15一 1A= -1/4 0專業(yè)知識分享2、已知f(1)=1.0，f(2)=12 f=3,則用辛普生(辛卜生)公式計算求得31 f (x)dx :用三點式求得f (1) ：答案：2.367, 0.25 3、f=T f(2) =2, f(3)=1 ,則過這三點的二次插值多項式中x2的系數(shù)為拉格朗日插值多項式為12(x-1)(x-2)1L2(x)(x -2)(x -3) -2(x -1)(x -3) 一答案：-1 ,24、近

2、似值x =0.231關(guān)于真值乂 = 0.229有(2 )位有效數(shù)字;5、設(shè)f(x)可微，求方程x=f(x)的牛頓迭代格式是()xn - f(xn )1- f (xn)6、對他廣二“十1，差商 f0,1,2,3=( 1 ), f0,123,4=( 0 )；7、計算方法主要研究( 截斷)誤差和( 舍入)誤差;8、用二分法求非線性方程 f (x)=0在區(qū)間(a, b)內(nèi)的根時，二分n次后的誤差限為b -a廣)9、求解一階常微分方程初值問題y = f (x,y) , y(X0)=y0的改進的歐拉公式為hyn 1 = ynf(xn,yn)f(xn；1，yn 1)(2)；10、已知f(1) =2, f (

3、2) =3, f(4) =5.9 ,則二次 Newton插值多項式中x2系數(shù)為(0.15 )11、兩點式高斯型求積公式11.1 r ,3-1 - 3 1廣汽虺=Jf(x)dxW口而)+ f(而")，代數(shù)精12、解線性方程組Ax=b的高斯順序消元法滿足的充要條件為(A的各階順序主子式均不為零)13、為了使計算4(x-1)26(X 1)3的乘除法次數(shù)盡量地少，應(yīng)將該表1x -1,為了減少舍入誤差，應(yīng)將表達式y(tǒng) =10 (3 (4 -6t)t)t,t達式改寫為_2&001 -41999 改寫為V2001 +1999314、用二分法求方程f(x)=x x-1 =0在區(qū)間0,1內(nèi)的根，

4、進行一步后根的所在區(qū)間為0.5、1 , 進行兩步后根的所在區(qū)間為 0.5 、 0.75。xdx15、計算積分0.5，取4位有效數(shù)字。用梯形公式計算求得的近似值為0.4268.用辛卜生公式計算求得的近似值為0.4309 、梯形公式的代數(shù)精度為1 .辛卜生公式的代數(shù)精度為 J；3x1 +5x2=1卜產(chǎn)=(1-5x2k)/316、求解方程組&.2x1 +4x2 =°的高斯塞德爾迭代格式為_&產(chǎn)=一x1(k，20 一該迭1代格式的迭代矩陣的譜半徑P(M)=12 _017、設(shè) f(0)=0, f(1)=16, f (2)=46，則 l1(x)= l1(x) = x(x2)_ f

5、(x)的二次牛頓插值多項式為由=16x 7x（x-1）18、求積公式af(x)dx r Akf(Xk)k=0的代數(shù)精度以(高斯型）求積公式為最高，具有（2n+1 ）次代數(shù)精度5、.f (x)dx19、已知f (1)=1, f (3)=5, f (5)=-3,用辛普生求積公式求,1=( 12 )。20、設(shè) f (1)=1 , f(2)=2 , f (3)=0,用三點式求 f'飛 2.5 )。21、如果用二分法求方程 x3 + x -4 = 0在區(qū)間1,2內(nèi)的根精確到三位小數(shù)，需對分(10)次。x30MxM 1S(x) =132(x -1)3 a(x -1)2 b(x -1) c 1 &l

6、t; x < 322、已知、2是三次樣條函數(shù)，則a=( 3),b=( 3),c=(1)。n:一 l k (x)=«(1),n* (x4 x23)lk(x) =4k=0( xy = f (x, y)24、解初值問題、y(x。)= y。的123、 l0(x),l1(x),，ln(x)是以整數(shù)點 x0,x1； ,xn為節(jié)點的Lagrange插值基函數(shù)，則 nxklj(xk) =xy( xj ), 當 n至2 時+ x2+3)。yn°L =yn +hf (xn,yn)|yn+ = yn +!f (xn, Yn )十“*口 書，y；*)i進歐拉法L2是2階方法。25、區(qū)間a,b

7、 I上的三次樣條插值函數(shù)S（x）在"b1上具有直到 2 階的連續(xù)導(dǎo)數(shù)。26、改變函數(shù)f （X） =，x+1 -豉 （x»1 ）的形式，使計算結(jié)果較精確3位小數(shù)，則需要對分1027、若用二分法求方程f（x）=0在區(qū)間1,2內(nèi)的根，要求精確到第次。28、S（x）=J 32x3,0MxM1ax2 bx c,1 E x 三 2 1 ，一 一 ,1 x 2是3次樣條函數(shù)，則a= 3 , b= -3 , c= 1。610 ,利用余項公式估計，至少用4771 exdx29、若用復(fù)化梯形公式計算°,要求誤差不超過個求積節(jié)點。x1 +1.6x2 =1 0 4x+x-2 30 、 與

8、出求解方程組 ”的 Gauss-Seidel 迭代公式'x!k* ) = 1 1.6x2k)遂尸)= 2+0.4xSk=0,1二迭代矩陣為0<0-1.6、一 0.641,此迭代法是否收斂 收斂用4)A =:.32、設(shè)矩陣33、一42若 f(x) =3x4276_a=lu，則 U =+ 2x +1,則差商 f2,4,8,16,30=334、1 .2 .f (x)dx f( -1) 8 f( 0) f (1)數(shù)值積分公式 T9的代數(shù)精度為235、線性方程組2的最小二乘解為36、設(shè)矩陣一3211145分解為A = LU ,則U11103212 1、單項選擇題:1、Jacobi迭代法解方

9、程組Ax = b的必要條件是（C ）A . A的各階順序主子式不為零B . P(A)<1Ca. #0,i =1,2,nD|心12 2-3A = 0 512、設(shè)L° 0 一71 則 P(A)為(C ).A . 2 B.5C. 7 D . 33、三點的高斯求積公式的代數(shù)精度為（B ）A . 2B. 5 C.3 D . 44、求解線Tt方程組Ax力的LU分解法中,A須滿足的條件是（B ）31、設(shè) 4 3,，則 1AL =色A.對稱陣 B.正定矩陣C.任意陣 D .各階順序主子式均不為零5、舍入誤差是(A ) 產(chǎn)生的誤差。A.只取有限位數(shù)B.模型準確值與用數(shù)值方法求得的準確值C.觀察與

10、測量D.數(shù)學模型準確值與實際值6、3.141580是冗的有(B )位有效數(shù)字的近似值。A.6B. 5 C . 4 D .77、用1 + x近似表示ex所產(chǎn)生的誤差是(C ) 誤差。A.模型 B .觀測C.截斷 D .舍入8、解線性方程組的主元素消去法中選擇主元的目的是(A ) oA.控制舍入誤差B .減小方法誤差C.防止計算時溢出D .簡化計算_x9、用1 + 3近似表示羽+x所產(chǎn)生的誤差是(D )誤差。A.舍入 B .觀測 C .模型 D.截斷10、-324. 7500是舍入得到的近似值，它有(C )位有效數(shù)字。A . 5 B . 6C. 7 D . 811、設(shè)f (-1)=1, f (0)

11、=3, f (2)=4,則拋物插值多項式中x2的系數(shù)為(A ) 0A.-0. 5 B .0.5 C . 2 D . -212、三點的高斯型求積公式的代數(shù)精度為(C )。A . 3 B . 4C. 5 D . 213、( D )的3位有效數(shù)字是0.236 X02。(A) 0.0023549 M03 (B) 2354.82 >10-2 (C) 235.418(D) 235.54 >10-114、用簡單迭代法求方程f(x)=0的實根，把方程f(x)=0表示成x=P(x),則f(x)=0的根是(B) o(A) y=中(x)與x軸交點的橫坐標(B) y=x與丫=邛僅)交點的橫坐標(C) y=

12、x與x軸的交點的橫坐標(D) y=x 與y=P(x)的交點13x - x2 - 4x 3 1一一 x1 +2x2 -9x3 = 015、用列主元消去法解線性方程組l-4x1 -3x2+x3=-第1次消元，選擇主元為(A )。(A) -4(B) 3 (C) 4(D)-916、拉格朗日插值多項式的余項是(B ),牛頓插值多項式的余項是(C )。(A) f(x,x0,x1,x2,xn)(x x1)(x x2) - (x xn 1)(x xn),(B)Rn(x) =f(x) - Pn(x) =f(n 1)()(n 1)!(C)f(x,x0,x1,x2,xn)(x x0)(x x1)(x x2)(x x

13、n 1)(x xn),(D)Rn(X) =f(x) Pn(X) =f(n - 1)!- -n i(X)17、等距二點求導(dǎo)公式f(x1) X A )。f(X1)-f(X。)(A)Xi -X0(B)f(X1)-f(X0)(C)f(x0) f(X1)(D)Xo - XiX0 - Xif(X1)- f(x。)Xi X018、用牛頓切線法解方程f(x)=0 ,選初始值x0滿足(A),則它的解數(shù)列xnn=0,1,2,一定收斂到方程f(x)=0的根。(A)f(x0)f (x) 0(B)f(xo)f(x) 0(C) f(Xo)f (x) ；0(D) f (Xo) f (x) < 019、為求方程x3x2

14、1=0在區(qū)間1.3,1.6內(nèi)的一個根，把方程改寫成下列形式，并建立相應(yīng)的迭代公式，迭代公式不收斂的是 (A )21, -1x =，迭代公式：Xk4=/(A)X -1Xk -1(B)X=1十'，迭代公式：Xk書=1 + X12Xk(C)x3=1+x2，迭代公式：X"2x1/3二(1Xk)(D)X3 -1=X2，迭代公式：Xk.二12Xk2,XkXk 17= f(x, y)20、求解初值問題 J(X =y。歐拉法的局部截斷誤差是();改進歐拉法的局部截斷誤差是()；四階龍格庫塔法的局部截斷誤差是(A )(A)O(h2)(B)O(h3)(C)O(h4)(D)O(h5)21、解方程組

15、 Ax = b的簡單迭代格式(1)二1, (2) M：1Bx + g收斂的充要條件是(P(A)a1,P(B)>122、在牛頓-柯特斯求積公式:bf(x)dx :an(b-a卜 Ci(n)f (Xi)C(n)用中，當系數(shù)Ci是負值時，公式的穩(wěn)定性不能保證，所以實際應(yīng)用中，當()時的牛頓-柯特斯求積公式不使用。(1) n >8, n >7, n >10,(4)n >6,23、有下列數(shù)表X00.511.522.5f(x)-2-1.75-10.2524.25所確定的插值多項式的次數(shù)是()。(1)二次；(2)三次；(3)四次；(4)五次24、若用二階中點公式y(tǒng)n 1 = y

16、nf(Xn , yn )求解初值問題y -2y, y(0)- 1 ,試問為保證該公式名對穩(wěn)定，步長h的取值范圍為()。 0 <h -1, (2)0Mh M1, ( 3) 0 ：h：1, (4)0M h ； 125、取6%1.732計算x =(6-1),下列方法中哪種最好?(A)28-16 石；(4-2百)2.26、已知3xS( x) =32(x -1)a(x -2) b_16_16_C) (4 + 2®)2 ；(D)(出 + 1)4。0M x <22-x - 4是三次樣條函數(shù)，則a，b的值為(A)6 , 6;(B)6, 8;(C)8, 6;(D)8, 8。27、由下列數(shù)表

17、進行 Newton插值，所確定的插值多項式的最高次數(shù)是()x11.522.533.5f(Xi)-10.52.55.08.011.5(A) 5；(B)4；(C)3；( D) 2。bf f(x)dxA f(x1) + A2 f(x2)+ A3 f( x3)28、形如a的高斯(Gauss)型求積公式的代數(shù)精度為(A)9；(B)7；( C) 5；(D)3O29、計算熱的Newton迭代格式為() xk3xk 3Xk i - Xk 1- -2x22x.(A) xk ； ( B)k ； (C)xk ； (D)30、用二分法求方程x3 +4x2 -10=0在區(qū)間1,2 內(nèi)的實根，要求誤差限為 次數(shù)至少為（）

18、（A）10 ;（B）12;（C）8;（D）9。13；=10231、經(jīng)典的四階龍格一庫塔公式的局部截斷誤差為（）_4_2_5(A)。；(B)0(h);( C)O(h);(D)O(h 3)932、設(shè)li（x）是以xk =k（k =0,1川為節(jié)點的 Lagrange(A) x；(B) k ；(C) i ；33、5個節(jié)點白牛頓-柯特斯求積公式，至少具有 （）插值基函數(shù)，則(D) 1。次代數(shù)精度、kli(k) =k =0(A)5；34、已知(A)6, 6;(B)4;(C)6;(D)3。x30 _ x _2S(x)=3.2(x -1) *a(x-2)+ b 2 < x<4是三次樣條函數(shù)，則(B

19、)6, 8;(C)8, 6;(D)8, 8。a,b的值為（）x0=2不收斂的是（）335、已知方程x -2x-5=0在x = 2附近有根，下列迭代格式中在2 x； 5(D) Xk +一 3 x2-2。x01234f(x)1243-5(Q3xk 1 = xkxk-5；36、由下列數(shù)據(jù)確定的唯一插值多項式的次數(shù)為（）(A) 4 ;(B)2;(C)1;(D)3。37、5個節(jié)點的Gauss型求積公式的最高代數(shù)精度為()(A)8 ;(B)9 ;(C)10;(D)11。、是非題（認為正確的在后面的括弧中打 &否則打父）1、已知觀察值（Xi，yi）（i =0"，m）,用最小二乘法求n次擬合

20、多項式Pn（x）時,Pn（x）的次數(shù) n可以任意取。（）2x2、用1- 2近似表示cosx產(chǎn)生舍入誤差。（）（X - X0 ）（ X - X2 ）()311、-2 5 3.112 5 一. ，,一.5、矩陣A=L 2 "具有嚴格對角占優(yōu)。()四、計算題：4x1 2x2 x3 =11x1 4x2 2x3 = 181、用高斯-塞德爾方法解方程組12" +x2 +5x3 =22,取x=(0，0，0)T ,迭代四次(要 求按五位有效數(shù)字計算)。答案：迭代格式丫 x1= -(11- 2x2k)-4x3k)e= l(18.xr).4-2x3k)x3k+)(222乂尸)-x")

21、kx1(k)x2k)(k) x3000012.75003.81252.537520.209383.17893.680530.240432.59973.183940.504202.48203.70191 一一 一 一一 1 一 1f(x)dx ; Af(-1) f Bf(-) f()2、求A、B使求積公式22的代數(shù)精度盡量I高，并求其代數(shù)精度；利用此公式求2 1 dxM x (保留四位小數(shù))。2答案：f(x)=1，x，x是精確成立，即2A 2B =22A 1B =223求積公式為i18.1.1Lf(x)dx=9f(f + f個匕葉當f(x) = x3時，公式顯然精確成立；當f(x) = x4時,

22、1右=3 。所以代數(shù)精度為3。2 1GdXt=2xJ3 11_出- t 39 -1 3 1 39-1/2 3 1 2 3970.692861403、已知xi1345f(xi)2654分別用拉格朗日插值法和牛頓插值法求f(x)的三次插值多項式P3(x),并求f(2)的近似值(保留四位小數(shù))L3(x) =2(x-3)(x-4)(x-5) 6(x-1)(x-4)(x-5)(1-3)(1-4)(1-5)(3-1)(3-4)(3-5).5(x-1)(x-3)(x-5) , 4(x-1)(x-3)(x-4) (4 -1)(4 -3)(4 -5)(5-1)(5 -3)(5 -4)xiyi一階均差二階均差三階

23、均差1236245-1-1差商表為54-10141P3(x) = N3(x) =2 2(x -1) -(x -1)(x -3) - 一(x -1)(x - 3)(x - 4)4f(2) ： P3(2) =5.54、取步長h=02,用預(yù)估-校正法解常微分方程初值問題y'= 2x +3y1y(0)=1(0<x<1)yn0)i=yn0.2(2xn3yn)答案:解:yni=yn0.1(2xn3yn)(2xni3yn0)i)yn 1 =0.52xn 1.78yn0.04n012345xn00.20.40.60.81.0y n11.825.879610.713719.422435.02

24、795、已知xi-2-1012f(xi)42135求f(x)的二次擬合曲線P2(x),并求f '(0)的近似值答案:解:ixiyi2 xi3 xi4 xixi yi2xi yi0-244-816-8161-121-11-22201000003131113342548161020015100343415a0 10a2 =1510ai =3正規(guī)方程組為2（制=”工小7101410ao +34a2 =4110311a。-, a1- -,a2-71014/、311P2（x）x1073f （0） : P2（0）=-106、已知sinx區(qū)間0.4 , 0.8的函數(shù)表xi0.40.50.60.70.

25、8yi0.389420.479430.564640.644220.71736如用二次插值求sin。63891的近似值，如何選擇節(jié)點才能使誤差最?。坎⑶笤摻拼鸢福航猓?應(yīng)選三個節(jié)點，使誤差M3|R2（X）|-J3| '3（X）|盡量小，即應(yīng)使 僧3（X）I盡量小，最靠近插值點的三個節(jié)點滿足上述要求。即取節(jié)點 0.5，0.6,0.7最好，實際計算結(jié)果Sin0.63891 : 0.596274且sin 0.63891 -0.596274 1 .< (0.63891 -0.5)(0.63891 -9-0.6)(0.63891 -0.7) 3!1-4< 0.55032 10-7、構(gòu)

26、造求解方程ex +10x-2=0的根的迭代格式Xn”中n =0,12，討論其收斂 性，并將根求出來，|Xn省一。|<10工。答案：解：令 f(x)=ex 10X-2, f(0) = -2；0, f (1)=10 e 0且f (x)=eX+10A0對VXW(-8,+叼，故f(x)=0在1)內(nèi)有唯一實根.將方程f (x) =0變形為X=110(2-eX)(x) =/(2 -eX)10 ,| ;： (x) | =10e ：二 110故迭代格式Xn 1 = (2。eXn10收斂。取X0 =0.5,計算結(jié)果列表如下:n0123xn0.50.035 127 8720.096 424 7850.089

27、 877 325n4567xn0.090 595 9930.090 517 3400.090 525 9500.090 525 008 4 一6* 4 一 4且滿足 |x7x6 產(chǎn)0.000 000 95<10所以 x % 0.090 525008x1 +2x2 +3x3 = 142xi 5X2 2x3 =188、利用矩陣的LU分解法解方程組I 3X1 +X2 + 5X3=20。WORD格式整理答案:解:A = LU一1213 -5 1 j|1 213 -4 -24令 Ly = b 得 y = (14,10,72)t9、對方程組3xi - 2x2 10x3 =15 * 10x1 4x2

28、X3 52xi +10x2 4x3 =8(1)試建立一種收斂的Seidel迭代公式，說明理由;(2)取初值 x(°)=(0,0,0)T ,利用(1)中建立的迭代公式求解，要求| x(k 1) - x(k) |二；10 3解：調(diào)整方程組的位置，使系數(shù)矩陣嚴格對角占優(yōu)10x1 - 4x2 - X3 - 5« 2x1 +10x2 -4x3 = 83x1 +2x2 +10x3 = 15故對應(yīng)的高斯塞德爾迭代法收斂.迭代格式為x"=x2k 1) =x3k J101(-2x101(-3x10(k 1) 14x2k)，x3k)，5)4x3k) 8)1(k1)-2x2k1)15)

29、取x=(0Q0)T,經(jīng)7步迭代可得:x* x(7) = (0.999 991 459, 0.999 950 326, 1.000 010)T10、已知下列實驗數(shù)據(jù)xi1.361.952.16f(xi)16.84417.37818.435試按最小二乘原理求一次多項式擬合以上數(shù)據(jù)1 x解：當 0<x<1 時，f "(x) =ex,則 f (x) =e,且 2 dx有一位整數(shù).專業(yè)知識分享WORD格式整理專業(yè)知識分享要求近似值有5位有效數(shù)字，只須誤差R(n)(f) <(b -a)312n2”(“,只要Rin) (ex):二 £_ 12n2即可，解得所以n=68,

30、因此至少需將0,1 68一1511、用列主元素消元法求解方程組解:一15：2-11-43 -12111r1r27一51R1(n)14(f) <-x10-2* 10-4= 67.30877等份。-1-4-4-1iirxii -4ix2-121Jl_X3jL11 J3 -121 41112 一15r 2r3 一5回代得12、取節(jié)點X0解:=0, X115135-4135251515513-1285795-12795513 1-'4-12 1r2-13515152579585X3 = -1，X2 =6, X1 =3o= 0.5,X2 =1，求函數(shù) f(X)=今(x),并估計誤差。在區(qū)間

31、0,1e4 (x- j (X-0)(X-1)(0-0.5)(0-1)(0.5-0)(0.5-1)上的二次插值多項式(x-0)(x.0.5)e (1 -0)(1 -0.5)= 2(x -0.5)(x-1)-4e-0.5x(x-1) 2e,x(x0.5)f (x) = e, f (x) = -e ",M 3 = max | f (x) |= 1又X0,1故截斷誤差/1但他近 e(x)1x(xs5)(xF13、用歐拉方法求,、 X J2 ,y(x)=0e dt在點x：0.5，1。1.5,2.0處的近似值。X 2解：y（X）=10e出等價于F 1y2 _x e eiy(0) =0( X &g

32、t; 0)記 f(x,y) =e"，取h=0.5, x0= 0,X1 -0.5, x2 =1.0, x3 =1.5, X4 =2.0則由歐拉公式y(tǒng)n書二ynhf (Xn ,yn)J0 =0n = 0,1,2,3可得y(0.5) ： y1 = 0.5,y(1.0) =、2 0.88940y(1.5) ： y3 =1.07334, y(2.0) = y4 ： 1.12604X14、給定方程 f(x)=(x-1)e -仁01）分析該方程存在幾個根;2）用迭代法求出這些根，精確到5位有效數(shù)字;3）說明所用的迭代格式是收斂的。X解：1）將方程（X-1）e 7=0（1）改寫為x 7 = e&qu

33、ot;x作函數(shù) f1(x)=X，f2(x)=e"的圖形(略)知(2)有唯一根x三(1,2)。2)將方程(2)改寫為x = 1 e，(k =0,12 )構(gòu)造迭代格式xk 中=1 +e"k=xo =1.5計算結(jié)果列表如下:k123456789xk1.223131.294311.274091.279691.278121.278561.278441.278471.278463) 平(x)=1+e-x中(x)=-e«當 xw1,2時，P(x)wW(2)W(1)u1,2,且| ： (x) |Me,：二 1所以迭代格式 xk41 =*(xk) (k=0,1,21為任意 xo71

34、,2均收斂。15、用牛頓(切線)法求)3的近似值。取xo=1.7,計算三次，保留五位小數(shù)。2解：也是f(x) =x 一3 = 0的正根,f'(x)=2x,牛頓迭代公式為x2 -3xn 1 = xn - 2xn,)Xn -1xn22Xn(n =0,12 )取xo=1.7,列表如下:n123xn1.732351.732051.7320516、已知f (-1)=2 , f (1)=3 , f (2)=-4 ,求拉格朗日插值多項式L2(x)及f (1 , 5)的近似值，取五位小數(shù)。(x-1)(x-2)(x1)(x-2)(x1)(x7)解:L2(x)=2 3 - 4 (-1 -1)(-1 -2)

35、(11)(1 -2)(21)(2-1)2 34= -(x-1)(x-2) -(x 1)(x-2)-(x 1)(x-1)3 23.1f(1.5) ： L2(1.5) = 0.04167241 exdx 17、n=3,用復(fù)合梯形公式求Jo的近似值(取四位小數(shù))，并求誤差估計。1 x1 -0 01 32 31e dx T3 =e 2(e e ) e ： 1.7342解：02 3f(x) =ex, f "(x) =ex, 0WxE1時，|f"(x)|Me|R|=|ex-T3 歸e12 32e=0.025 0.05108至少有兩位有效數(shù)字。18、用Gauss-Seidel迭代法求解線

36、性方程組取x(0)=(0,0,0) T,列表計算三次，保留三位小數(shù)。解：Gauss-Seidel迭代格式為：婢劃=1(-x3k)+5)3.尸:-夕一尸-x3k)-1)x3k+)(k+)+x2k+)-8)、43011 -31系數(shù)矩陣1T41嚴格對角占優(yōu)，故Gauss-Seidel迭代收斂.取x(0)=(0,0,0) T,列表計算如下：kx1(k)x2k)x3k)11.6670.889-2.19522.3980.867-2.38332.4610.359-2.52619、用預(yù)估一校正法求解 )=1 (0<x<1) , h=0。2,取兩位小數(shù)。解：預(yù)估一校正公式為1,、yn 1 = yn2

37、(k1k2)&1 =hf (Xn,yn)k2 =hf (xn h,yn *1)20、n = 0,1,2,解:2 .q，= span1, x AT解方程組1252AT AC131211 382=19.0 32.3 49.0 73.3其中解得:21、（ATA= 339l33913529603t 173.6 ATy =179980.70.9255577C =0.0501025所以 a =0.9255577,15分)用n=8的復(fù)化梯形公式(或復(fù)化Simpsonb =0.0501025公式）計算dx一e時，試用余項估計其誤其中 f(x,y)=x+y, y0=1 , h=0.2, n=0，1,2,

38、3,4,代入上式得:n12345xn0.20.40.60.81.0yn1.241.582.042.643.42（8分）用最小二乘法求形如 y = a *bx的經(jīng)驗公式擬合以下數(shù)據(jù):x19253038yi19.032.349.073.3差。用n =8的復(fù)化梯形公式（或復(fù)化Simpson公式）計算出該積分的近似值。解:R"f =-b12ah2f P)111e00.00130212 82768h7T(8)=)f(a) 2% f(xQ f(b)2 k 11=11 2 (0.8824969 0.7788008 0.606530660.5352614 0.47236655 0.41686207)

39、 0.36787947= 0.632943422、（15分）方程x3 -x-1 =0在x= 1.5附近有根，把方程寫成三種不同的等價形式（1） x = 3/x+1_ 3 1 x = . 1 ： -xn 1對應(yīng)迭彳t格式xn41 （2） x x對應(yīng)迭代格式113xn ; (3) x = x3-1 對應(yīng)迭代格式“斗=xn -1。判斷迭代格式在x0 =1.5的收斂性，選一種收斂格式計算 x = 1.5附近的根, 精確到小數(shù)點后第三位。解：（1）1f:(x) (x 1) 336（1到=0.18<1,故收斂;(2)(3)（x）二2x2 ：1 +1x,自（到=0.17<<故收斂;x）=3

40、x2,仔。5）二5、1,故發(fā)散。選擇(1”x°=1.5x=1.3572x2=1.3309x3=1.3259x，=1.3249x5 =1.32476 x6 =1.3247223、（8分）已知方程組 AX = f，其中一431 一 241A = 3 4-1 f = 30-14 J 24（1） 列出Jacobi迭代法和Gauss-Seidel迭代法的分量形式。（2） 求出Jacobi迭代矩陣的譜半徑。解：Jacobi迭代法:婿川=1(24-3x2。)4聞)=：(30-3x1(k) +x3k)x尸=1(-24 + x2k)4k=0,1,2,3:Gauss-Seidel 迭代法:x，*) =(

41、24-3x2k)4x2k *)=；(30-3x1(k +)+x3k)x3k 卻(_24 + x2由4. k = 0,1,2,31Bj - -D (L U)=-340340 % °,：(Bj.10(或)=0.79056943 = _y+1d dx24、1、(15分)取步長h=0.1,求解初值問題I y(0)=1用改進的歐拉法求y(0.1)的值；用經(jīng) 典的四階龍格一庫塔法求y(0.1)的值。y2 二 yhf (Xn, yn) =0.9yn 0.1h(0)yn 1 = yn f (Xn, yn) f (Xn 1, yn 1 ) = 0.905丫口 0.095解：改進的歐拉法：2所以 y(O

42、。=y1 =1 ； 經(jīng)典的四階龍格一庫塔法：hyn 書=Yn +M匕 +2k2 +2k3 +k4 6k1 = f (Xn, yn),一 , h, h,、k2 = f (Xn +)n + k1)22hhk3 = f (Xn 十萬，Yd 十 ”2)、 k4 = f (Xn + h, yn + hk3)k1 = k2 = k3 = k4 = 0 所以 y(0.1) = Y1 = 1。25、數(shù)值積分公式形如1Xf(X)dX : S(x) =Af (0) Bf(1) Cf (0) Df (1)A B D0試確定參數(shù)A B,C, D使公式代數(shù)精度盡1f(x) C401 ,， R(x) = °xf

43、(x)dxS(x)量高；(2)設(shè)f(X) C 0,1,推導(dǎo)余項公式0，并估計誤差。a 3 r 7 r 1c解：將f(X) =1,X,x2,x3分布代入公式得:A = , B = , B = , D =20202030H3(x) = f(Xi) -構(gòu)造Hermite插值多項式H 3(x)滿足F3(為) = f (xi)i =0,1其中X0 = 0,X1 = 1 f ()32R(X) = 0xf (x) -S(x)dx = 04X (x -1) dx1一、一，、 f (4)( ) 2, 八 2xH3(x)dx=S(x) f(x) -H3(X) = X (X -1)則有：03,4!f(4)()4!1

44、 32 f ()x (x -1) dx =(-04! 60f(4)()144026、用二步法Yn 1 = ： 0Yn . ： 1Ynq h寸(Xn,Yn) (1 -？)f(Xni,Yn7)：y'= f (x, y)=并求局求解常微分方程的初值問題Ly(X0) =Y0時，如何選擇參數(shù)"0,a1出使方法階數(shù)盡可能高,部截斷誤差主項，此時該方法是幾階的 解:h2h3Rn,h = y(xn 1) - yn 1 = y(xn) hy (%) 可 y (%) 可 y (%) 2!3!h2h3-二 0y(xn) - ： 1(y(xn) -hy (xn)-y (xn)y (xn)2!3!比曠

45、（4） （1-口）（丫（4）-八丫乂）23五 y (xn) -3 y(4)(xn)=(1 - ： 0 - ： 1)y(xn) h(1 -1 ： 1)y (xn)2 1-1.3 1h2(-V 1-)y (xn) h3(-226)y (Xn) O(h4)所以21 T-0 - -1 = 0、工1 = 0-三 1 -1-02“0=1=y =03 =-223h3主項：12y (Xn)該方法是二階的。27、(10 分)已知數(shù)值積分公式為:h0 f(x)dx :h .2-f(0) f(h)-h2f (0) -f (h)2，試確定積分公式中的參數(shù)，使其代數(shù)精確度盡量高，并指出其代數(shù)精確度的次數(shù)。解：f（x）

46、=1顯然精確成立；f(x) = x 時,hxdxh2 h2=h0 h h21 -122f(x)2二x時,h 2x dxf(x)3x時,h 3x dxh 22h30 h2h20 -2h2 h =22f(x)所以,4x時,)4hh 54x dx =)5h 3122= U0 h3 h20 -3h2212.h 4123 h0 h h 0 - 4h ：2126 .其代數(shù)精確度為3。1 , xk 1 = -(xk228、（8分）已知求'；a（a >0）的迭代公式為:a)x00 k = 0,1,2證明：又一切k =1,2， 從而迭代過程收斂。, xk之石，且序列&,是單調(diào)遞減的,Xka

47、 = Jak =0,1,2 -xk1 a 1xk 1 =-(xk -)-2證明：xkxk2故對一切k =1,2,，xk之7a。2 j 馬至 1（1 1）=1；.，又xk2 xk 2所以xk41 -xk,即序列xk是單調(diào)遞減有下界，從而迭代過程收斂。33f（x）dx 二f f（2）29、（9分）數(shù)值求積公式L02是否為插值型求積公式？為什么？其代數(shù)精 度是多少?解：是。因為f（x） 在基點1、2處的插值多項式為x - 2 P(x)= 1 2x -1f二1f330p(x)dx =a1f(1) f (2)其代數(shù)精度為30、（6分）寫出求方程4x=cosx）+1在區(qū)間0,1的根的收斂的迭代公式，并證明其收斂性。xn 1(6分)二xn=1 1 cos xn 14, n=0,1,2, L 1. ， 1儼(x ) = /sin (x )<1：14對任意白初值x0 0,1 ,迭代公式都收斂。31、（12分）以100,121,144為插值節(jié)點，用插值法計算,