數(shù)值計(jì)算答案-石瑞民

1、習(xí)題一1、取3.14,3.15, 22,啰作為n的近似值,求各自的絕對(duì)誤差,相 7113對(duì)誤差和有效數(shù)字的位數(shù).解：x1 =3.14121 一 1 J3n -Xi <- X10 =一父10 22所以,木有三位有效數(shù)字絕對(duì)誤差：e = n -3.14,相對(duì)誤差：er-3.14冗絕對(duì)誤差限：£ JM10： 相對(duì)誤差限：6r=,M10H=1M10/ 22 36X2 =3.15r： -3.15 =0.00840174 = 0.84074 10 <0.5 10'=0.5 101立所以,X2有兩位有效數(shù)字絕對(duì)誤差：e = n -3.15,相對(duì)誤差：er =-315絕對(duì)誤差限：

2、-2,101相對(duì)誤差限：皆?1.22X2=Y22M0.5 10/ =0.5 10147所以,X3有三位有效數(shù)字一=0.0012645 =0.12645M10_22絕對(duì)誤差：e=n-22,相對(duì)誤差：er7二絕對(duì)誤差限：=工M104,相對(duì)誤差限：露=工父1026X1 =355113<0.5 104 =0.5 101355_6-=0.00000032 =0.32 10113所以,X4有七位有效數(shù)字355JI 絕對(duì)誤差：e=n-竺5,相對(duì)誤差：a=-113113二絕對(duì)誤差限：,區(qū)相對(duì)誤差限:=610-63、以下各數(shù)都是對(duì)準(zhǔn)確數(shù)四舍五入后得到的近似數(shù),試分別指出它 們的絕對(duì)誤差限和相對(duì)誤差限,有效

3、數(shù)字的位數(shù).Xi =0.0315,X2 =0.3015,X3 =31.50,X4 =5000解：X1 =0.0315m=-1*X 一 X11 工- 1021j _3=102X2X3所以,n=3, X1有三位有效數(shù)字絕對(duì)誤差限： 仁工h.、 相對(duì)誤差：名r 2= 0.3015m=0* 1_4 10 _4X -X2 <-10 =父10 22所以,n=4, X1有四位有效數(shù)字絕對(duì)誤差限：名二1".'相對(duì)誤差："2=31.50m=2=7產(chǎn)110X44、* 1_2 12 _4X -X3 <-x10 =_ 父10 22所以,n=4, X1有四位有效數(shù)字絕對(duì)誤差限：君

4、=1"0：相對(duì)誤差：力2=5000m=4* - 10 14-4x -x4 <-10 =一父1022所以,n=4, X1有四位有效數(shù)字絕對(duì)誤差限：/M100=0.5,2相對(duì)誤差：；r =工 10" 1 =10, =10,2a2 5計(jì)算而的近似值,使其相對(duì)誤差不超過(guò)U優(yōu)7100.1%解：設(shè)取n位有效數(shù)字,由定理1.1知,%=工父102a由向=10父0.3162,所以,a1=3由題意,應(yīng)使父10"*<0.1%,即10父10n <6父106所以,n=4,即?而的近似值取4位有效數(shù)字近似值x =3.1626、在機(jī)器 數(shù)系下F(10,8,L,U)中取 三個(gè)數(shù)

5、x = 0.23371258父10, y =0.33678429 父10 Q1 Q同理；(x2)10；(x3)10 2 , z =0.33677811 父 102 , 試按(x + y) + z 和 x + (y + z)兩 種算法計(jì)算x + y+z的值,并將結(jié)果與精確結(jié)果比較.(x y) z_4_22二 (0.23371258 100.33678429 10 ) -0.33677811 10 一 _ _2 _ _ 2、 _2解. =(0.00000023371258 100.33678429 10 ) -0.33677811 10= 0.33678452371258 102 -0.33677

6、811 102= 0.33678452 102 -0.33677811 10223= 0.00000641 102 = 0.64100000 10x (y z) _4 2 _2= 0.23371258 10(0.33678429 10 -0.33677811 10 )43= 0.23371258 100.61800000000 10= 0.023371258 10 4 0.61800000000 10-0.641371258 10,一一 一 一一 一2-0.64137126 10x y z _4 _2 _2= 0.23371258 100.33678429 102 -0.33677811 10

7、2_2_2_2= 0.00000023371258 102 0.33678429 102 -0.33677811 102= 0.00000641371258 102= 0.64137126 10 ,所以,x+(y+z)比(x + y)+z 精確,且 x+(y+z)與 x + y+z 相同；因此,在做三個(gè)以上的數(shù)相加時(shí),需要考慮相加的兩個(gè)同號(hào)數(shù)的階數(shù) 盡量接近.8、對(duì)于有效數(shù)x=-3.105, x2 =0.001, x3 =0.100,估計(jì)以下算式的相對(duì)誤差限. y1=x+x2+x3, y1=x1x2x3, 丫3=迎x3解：x二-3.105, m=1;*1_311 4x -x1 E1父10 =父

8、1022z 、13所以；(xj = 2 10i 一 3e( Xi) W Mio1 一 3e(X2) W 父 I02e(X3) <i xi0、e Xier(Xi X2 X3)=er(Xi)er(Xi)erX3e(Xi)Xie(X2)X2e(X3)X31io-、i2 或；r(XI)io學(xué)-3.io252 3i_3-io< 或即(x2 )o.ooi-io,2-或；r(X3)o.iooX2X3e xi r e X2 廣 e X3XiX2X3X1X2X3i_3 io2 ier(yi) = e.(Xi +X2 +X3) W o.49975M io.er(y2) =0(XiX2X3) =e(Xi

9、X2) e(X3) =e(Xi) e(X2) e)所以,er(y2) <o.5o5i6er ( y3 ) = er () = er (X2 ) - er ( X3 )X3所以,4也)| <o.5o5綜合得： /(y1)=o.49975xio3 ,即(y2) =o.5o5i6 , "(y3) =.5.59、試改變以下表達(dá)式,使其結(jié)果比較精確其中x «i表示x充分接近0, |x Ai表示X充分大.(1)In xi _ In x2,xi : x2(2)i -xif i x(3)(4)i -cosX(5)i一一 cot X , X答案:(1) In 4；(3)X2(4)

10、法一：用i-c0sx十得出結(jié)果為:法二:i - cosx i - cosxsin x i - cosx xi 一 cosxsin xsin xsin x sin xsin xi cosx(x > o)或=3Xsin xx =tan212、試給出一種計(jì)算積分In =e.jxnexdx近似值的穩(wěn)定性遞推算法解：顯然,In>0,n=1,2,當(dāng) n=1 時(shí),得,I1 = fxexJLdx = 1 0e當(dāng)nA 2時(shí),由分部積分可得： 1 n -1n =jxe dx=1-n1n二,n=2,3,另夕卜,還有：In = fxnexJ1dx< fxndx = - 00 n 1由遞推關(guān)系In=1

11、-nIn-1 ,可得計(jì)算積分序列In的兩種算法: In =1 -nIn,n=2,3 -1-1 InL- n = 2,3,n下面比較兩種算法的穩(wěn)定性 假設(shè)的一個(gè)近似值In,那么實(shí)際算得的In的近似值為I n = 1 - n I n 所以,In -In = (-n)( I n- I n)I n - I n =nIn一In由此可以看出In的誤差放大n倍傳到了 In,誤差傳播速度逐 步放大由In計(jì)算In In,=n = N,N -1,1 n 假設(shè)In的一個(gè)近似值是In,那么實(shí)際計(jì)算的In的近似值為 1 - InI n 4 =n 所以,1n二- 1n 4(I n - I n)n1I n一 I n-I n

12、 - I nn由此可以看出In的誤差將縮小n倍傳到了 In,誤差傳播速度逐 步衰減.1綜上可看出,計(jì)算積分I n = eI xnexdx的一種穩(wěn)定性算法為1 - InInJL =n n = N,N -1,N -2 ,1.n習(xí)題二1、利用二分法求方程x3. =ln(2cosx)x2-4s-7 = 03, 4內(nèi)的根,精確到1.)口1 1.,即誤差不超過(guò)2.解：令 f(x) =x3 -2x2 -4x-7f =1.<., f(40=9.,說(shuō)明在3, 4內(nèi)有根, 利用二分法計(jì)算步驟得出 x10 =3.632324219 , x11 = 3.63218359381b11 -a11 =x11 -x10

13、 =0.4882181x10<父10 滿足精度要求 2所以,x:xn =3.6321 ,共用二分法迭代11次.2、證實(shí)1-x-sinx=0在0,1內(nèi)有一個(gè)根,使用二分法求誤差不大于1M10工的根.2證實(shí)：令 f (x) = 1 - x - sin xf(0) =1 >0; f (1) =-sin1 <0 ,所以,f(0) f (1) <0由零點(diǎn)定理知,f(x)在0,1內(nèi)有一根根據(jù)計(jì)算得出：x* % x15 =0.98283 ,此時(shí)共迭代15次.4、將一元非線性方程2cosx-ex=0寫成收斂的迭代公式,并求其在x0 =0.5附近的根,精確到10,解：令 f (x) =

14、2cosx-ex令f (x) =0,得到兩種迭代格式x(p *ej 匕=arccosx(x):-e,不滿足收斂定理.2) (x)= -2s1nx 二4口乂 2cosx叼(Xo)=率2(0.5) =0.008727 <1 ,滿足收斂定理由方程寫出收斂的迭代公式為xk 1 = ln(2cosxk)取初值為x° =0.5 ,得出近似根為：x*之x2 =0.693074175、為方程x3 x2.1=0在x0 =1.5附近的一個(gè)根,設(shè)方程改寫為以下等價(jià)形式,并建立相應(yīng)的迭代公式：(1)x =1 + 4 ,迭代公式xk由=1 +' ；xxk(2) x3=x2+1,迭代公式 乂3乂2

15、+1)1/3(3)迭代公式xk 1 =(xk -1)解：(1)利用局部收斂定理判斷收斂性,判斷初值x0 =1.5附近的局部收斂(2)局部收斂(3)不滿足局部收斂條件但由于汽,(x)|.|叼(x),所以%(x)比%(x)收斂的慢取第二種迭代格式x.1 =(xj 1)1/3取初值x0 =1.5 ,迭代9次得x* 5tx9 = 1.4667、用牛頓法求解x3-3x-1=0在初始值% =2臨近的一個(gè)正根,要求人書xk <10 O由牛頓迭代法知：Xk 1 = Xk -3,f (Xk)3Xk 1 .,、 一 2f (Xk)3(Xk -1)迭代結(jié)果為:Xk1.888891.8794593987滿足了精

16、度要求,X* ：. X3 =1.879398、用牛頓法解方程1-C=0,導(dǎo)出計(jì)算C的倒數(shù)而不用除法的一種X簡(jiǎn)單迭代公式,用此公式求0.324的倒數(shù),設(shè)初始值x0=3,要求計(jì)算 結(jié)果有5位有效數(shù)字.解：f(x)= 1-0.325 Xf (x)=-3,由牛頓迭代公式Xk Xk-fx"X'f (Xk)迭代結(jié)果為：k0123Xk33.083.086413.08486420滿足精度要求* X : x3 = 3.0864所以,0.324的倒數(shù)為3.086411、用快速弦截法求方程X3-3x-1=0在X0=2附近的實(shí)根,(取X1 = 1.9, 要求精度到10).解：f (x) = x3 -

17、3x -1 ,迭代結(jié)果：Xk2滿足精度要求*x <- X4 = 1.879391.91.88101.87949411601.8793912、分別用以下方式求方程4cosx=eX在x°=±附近的根,4有效數(shù)字要求有三位(1)用牛頓法,取Xo =-4(2)用弦截法,取X.二三X1二三42(3)用快速弦截法,取X.二三%;三42解：求出的解分別為：X1 =0.905 X2 =0.905 X3 =0.905習(xí)題三1、用高斯消元法解以下方程組2X1 -x2 +3x3 = 1(1 )-4x1 +2x2 +5x3 = 4x1 2x2= 7工 11X1 -3X2 -2X3 =3

18、71; -23x1 + 11x2 +x3 =0x1 +2x2 +2x3 = -1(2)解：(1)等價(jià)的三角形方程組為4x1 2x2 5x3 = 4« -2x2 +0.5x3 = -1 ,回代求解為 *721一二 X3 = 一、84(2)等價(jià)的三角形方程組為x1 = 9x2 = -1X3 - -6-23x1 11x2 x3 =05747 X22323193X3 = -1 ,回代求解為57223X3 :5741Xi -193106X2193223 x31932、解：L二,1020-0 111* A =2 0-11C 0111 000、0 111 U =2 010 ,U-10 0 1、5)

19、3、用 LU作LU分解.j,1 020、0 1110 0-510 006.5)5768方程組7 10579108610 Yx1、1、9x217x315>1x4 /工解:,16/5L =7/5L10010-1/2103/520 0、5107-2/500910-4/5-3-5-17/201/10-1/5,X-1/2(3/10-12-534、用列主元的三角分解法求解L方程組-12-2-T解：A=3-147、2-3-20100、/3L = 2/31/31-7/50 ,U = 00；<0-1-7/30x12 x2 - 2x3 - - 13x1 - x2 4x3 = 72x1 - 3x2 -

20、2x3 = 04-14/3 ,Y =-4 j7-14/3 , X =-221<1/2>5、 用追趕2100-12-10A= 0-1 2-100-12100-1角方程組Ax = b-12T0000J2-1-1/23/2 -1解：L =1/2-2/31-3/4-4/5 14/3-15/4 -16/51/31/4J/55/6、2/31/21/3J/66.用改進(jìn)的Cholesky分解法解方程組4x1 -2x2 -4x3 =10,-2x1 +17x2 +10x3 =3-4x1 +10x2 +9x3 = -77、用改進(jìn)的cholesky分解法解方程組1-13-1X2X38-44-1解：U =I

21、011/4-3/425/450/1178/25-6/11456/25-1100、4 -2 -4、10、2、解：L =-1/210,U =0 168,丫 =8,X =11-11/2 1 ,d 0 1J1/-1I 1J1-1047X18、設(shè)x=(1,-2,3)T,求岡乂2和網(wǎng)厘.xn解：x1bi 1xmaxx =31 19、設(shè) A = 2 22 40-3,求A|,網(wǎng)2和同g1 ,解：11Ali =8, I A0G = 10, |A2 =v'P(ATAT) =7.14171 110、設(shè) A =22<5 45X= 3 ,計(jì)算 |乂8,|Alc及 |AxL,<2 J并比較Ax二和也,

22、網(wǎng)8的大小.解：也=3,網(wǎng)L=1.,|Ax|g=9,1-22 Yx1 11-12、11、給定方程 -111 X2 = 02 -21 AX3 ； <10 >(1)寫出Jacobi和Gauss-Seide迭代格式；(2)證實(shí)Jacobi迭代法收斂而Gauss-Seide迭代法發(fā)散；(3)給定x(0) =(0,0,0)T ,用迭代法求出該方程的解,精確到|x(k+)-x(k)L<10解：x1 =(1) Jacobi 迭代公式x2 = x1風(fēng)=2x12x2-2x3 -122x210Gauss-Seide迭代公式4(k 1)(k)(k)x1 2x2- 3x3-12(k 1)(k)(k)

23、-x2 2x2-12(k 1)(k)(k) ooxn 8x2- 6 x3- 38(3)用 Jacobi迭代得,X* 定X(4) =(12,46,58)T5x1 一 x2 - x3 "_413、 rx1 +10x2 -x3 -x4 =12,考察 Jacobi迭代格式和 Gauss-Seidel - x 一 x2 5x3 -x4 = 8-x1 一x2 -x3z 10x4 =34迭代格式的收斂性.14、方程組Ax=b,其中1 a a'A= 4a 1 0 , x,bw R3<a 0 1利用迭代收斂的充分必要條件確定使 迭代法均收斂的a的取值范圍.0- a解:Jacobi迭代矩陣

24、為Bj = -4a 0、a 0當(dāng) P(Bj )<1 得,a < 5Jacobi 迭代法和 Gauss-Seidel-a00Gauss-Seide迭代矩陣為:當(dāng)儀Bs)<1得,a兩5430 Yx1 "15、設(shè)方程組3 4 -1 ix22 一14 人x3/Id-a-aBj =04a24a2922aa24 '30分別用Gauss-Seidel迭代法和24,w=1.25的SOR法求解此方程,準(zhǔn)確到4位有效數(shù)字(取x(0) =(l,l,l)T) 解：Gauss-Seidel迭代法共迭代17次,此時(shí)近似解為*(17)-Tx x( ) =(3.000,4.000,5.00

25、0)SOR法w=1.25時(shí),迭代11次,此時(shí)的近似解為x* ： x(11) =(3.000,4.000,5.000)T16、用SOR方法解方程組(分別取松弛因子 w=1.03, w=1, w=1.1)J_ 4xi -x2=1« xi +4x2-3x3 =4 精確解 x* =(1/2,1,1/2),要求當(dāng) |x* - x(k)| < 5 10 -x2 +4x3 = -3時(shí),終止迭代,并且對(duì)每一個(gè)w值確定迭代次數(shù).解：當(dāng) W=1.03 時(shí),迭代 5 次,x* 鄭 x(5) =(0.5,1,-0.5)T當(dāng) W=1 時(shí),迭代 6 次,x* «x(6) =(0.5,1,-0.5

26、)T當(dāng) W=1.1 時(shí),迭代 6 次,x* %x(6) =(0.5,1,-05),習(xí)題四1、設(shè)xo =0, xi =1 ,與出f(x)=e.的一i次插值多項(xiàng)式Li(x),并估計(jì)插值 誤差.解：Li(x) = yo *-也(x - xo) = 1 J -1)x xi -xoeI Ri(x)|EM|(x x0)(x xi),其中 M = max f "(x) =12%.咤一,、一 1 ,、,、 J1 Ri (x) F (x Xo)(x Xi) 一二282、給定函數(shù)表Xi-0.10.30.71.1f (Xi)0.9950.9950.7650.454選用適宜的三次插值多項(xiàng)式來(lái)近似計(jì)算f (0

27、.2)和f (0.8)解：、求 f (0.2),選用插值節(jié)點(diǎn)為 X0 =-0.1 , xi =0.3, X2 =0.7,用 lagrange插值多項(xiàng)式為：,(x -Xi)(X -X2)(x -Xo)(X -X2)(x -Xo)(X- Xi)L2 (x):yo yi y2(Xo -Xi)(Xo -X2)(Xi -Xo)(Xi -X2)(X2 - Xo)(X2 Xi)解得 f(-0.i) :. L2(-0.i) =0.979、求 f(0.8),選用插值節(jié)點(diǎn) x0=0.3, xi=0.7, X2=i.i,L2(X) n(X - xi )(X - x2 )(X0 - Xi )( X0 - X2 )y.

28、(x - x0)( X - x2)( X - x0 )( X - Xi)-yi y2(Xi - X0)( Xi - X2 )(X2 - X0 )( X2 - Xi )解得：f (0.8) ：. L2(0.8) =0.69754、給定數(shù)據(jù)(f(x) =4)X 2.02.i2.22.4f(xj |i.i42i4 i.449i38 i.48320 i.549i7 (i)試用線性插值計(jì)算f(2.3)的近似值,并估計(jì)誤差.一 (2)試用二次Newton插值多項(xiàng)式計(jì)算f (2.i5)的近似值,并估計(jì)誤差 解：(i)取 X0 =2.2 , Xi =2.4Li(x)= yyi -y0 (x -x0) -0.3

29、2995X 0.7573iXi -X0f(2.3) : Li(2.3) =i.5i6i95|Ri(x)|<M(x-x0)(x-xi) , M = max f"(x) =0.0766 2x0<?<?i0 0766| Ri(2.3)|< 00266 (2.3 -2.2)(x -2.4) =0.000383i(2)寫出二次Newton插值差商表Xif(Xi)一階差商二階差商2.0i.i42i42.ii.449i380.349242.2i.483200.34062-0.043iN2(x) =i.4i42i4 0.34924(x -2) -0.043i(x - 2)(x

30、 - 2.i) f(2.i5) ： N2(2.i5) = i.4663 &(2.i5) M 0.000004i435、給出函數(shù)值X0i234y0i646880試求各階差商,并寫出Newton插值多項(xiàng)式和差值余項(xiàng) 解：f (6)()6!(x -0)(x-1)(x-2)(x -4)(x -5)xiy一階差商二階差商三階差商四階差商001161624630738821-3-5/240-88-109/3-25/2-7/6N4(x) =16x -7x(x -1) -5/2x(x -1)(x -2) -7/6x(x - 1)(x - 2)(x -4)R(x)= f (x) -N4(x)= fxo,

31、x1,x2,x3,x4,xw5(x)6、給定數(shù)據(jù)表x 0.1250.250.3750.5000.6250.750f(x)| 0.796180.773340.743710.704130.656320.60228試用三次牛頓差分插值公式計(jì)算 f (0.158)和f (0.636).解：、求 f (0.158),取 x0 =0.125, x1= 0.25, x2 = 0.375, x3 = 0.500 , h=0.125差分表為xif (xi)一階差分二階差分三階差分0.1250.796180.250.77334-0.022840.3750.74371-0.02963-0.006790.50.704

32、13-0.03958-0.00995-0.00316由公式riwJykL疝由牛頓插值公式有f(0.158) : N3(0.158) =0.79061、求 f (0.636 ),取 x0 =0375 , x1 =0.500, X2 = 0.625, X3 = 0.750 , h=0.125xif(xi)一階差分二階差分三階差分0.3750.743710.50.70413-0.039580.6250.65632-0.04781-0.008230.750.60228-0.05404-0.006230.002求解得 f 0.636) : N3(0.636) =0.651799、給出sinx在0,pi的

33、等距節(jié)點(diǎn)函數(shù)表,用線性插值計(jì)算sinx的近似值,使其截?cái)嗾`差為1M10,問該函數(shù)表的步長(zhǎng)h應(yīng)取多少才能滿足 2要求解：設(shè)插值節(jié)點(diǎn)為xj =ih , (i=0,1h) ,h = ±n.h2. h2由 Rn(x) < f (x) =m282F(x)=sinx, f"(x)=-sinx,所以 f"(x)W1,即 m2 =1所以h三0.02步長(zhǎng)h應(yīng)取為0.02才能滿足要求.14、實(shí)驗(yàn)數(shù)據(jù)如下xi1925313844y19.032.349.073.397.8用最小二乘法求形如y=a+bx2的經(jīng)驗(yàn)公式,并計(jì)算均方差解：設(shè)擬合多項(xiàng)式為y=a+bx2,那么正規(guī)方程組為Ya！

34、08 S3S4 人b J j2jS0S1S2S1S2S31575327、即:1575327532719233119233107277699 人b ,2 271.49776.1©69321.5,a =0.968b =0.05所以,經(jīng)驗(yàn)公式為：y =0.968 0.05x2均方誤差為0.00301915、觀測(cè)物體的直線運(yùn)動(dòng),得出以下數(shù)據(jù)時(shí)間t(s)00.91.93.03.95.0距 離S(m)010305080110求運(yùn)動(dòng)方程.解：設(shè)擬合多項(xiàng)式為y=a+bx + cx2,那么正規(guī)方程組為S2S2aT0S3S30S4bT253.63218.9079510302328010780533.2 ,614.7即： 14.753.6363.63 218.907a=-0.5834, b=11.0814, c=2.2488所以擬合多項(xiàng)式為 y =-0.5834+11.0814x + 2.2488x2.習(xí)題五1、分別用梯形公式和辛普森公式計(jì)算以下積分,并比較結(jié)果.(1) f；dx (n=8)0 4 x2解：用復(fù)合梯形公式 h Jn =8,