數(shù)值分析習(xí)題答案

本文格式為Word版，下載可任意編輯——數(shù)值分析習(xí)題答案第一章緒論

3．以下各數(shù)都是經(jīng)過四舍五入得到的近似數(shù)，即誤差限不超過最終一位的半個(gè)單位，試指出它們是幾位有效數(shù)字：x1?1.1021,x2?0.031,x3?385.6,x4?56.430,x5?7?1.0.解：x1?1.1021是五位有效數(shù)字；

*x2?0.031是二位有效數(shù)字；*x3?385.6是四位有效數(shù)字；*x4?56.430是五位有效數(shù)字；*x5?7?1.0.是二位有效數(shù)字。

******4．利用公式(2.3)求以下各近似值的誤差限：(1)x1?x2?x4,(2)x1x2x3,(3)x2/x4.其中x1,x2,x3,x4均為第3題所給的數(shù)。解：

************121*?(x2)??10?321*?(x3)??10?1

21*?(x4)??10?321*?(x5)??10?12?(x1*)??10?4***(1)?(x1?x2?x4)***??(x1)??(x2)??(x4)111?4?3?3??10??10??10222?1.05?10?3***(2)?(x1x2x3)*********?x1x2?(x3)?x2x3?(x1)?x1x3?(x2)111?1.1021?0.031??10?1?0.031?385.6??10?4?1.1021?385.6??10?3222?0.215

**(3)?(x2/x4)?****x2?(x4)?x4?(x2)*x42110.031??10?3?56.430??10?322?56.430?56.430?10?56．設(shè)Y0?28，按遞推公式Y(jié)n?Yn?1?

1783（n=1,2,…）100計(jì)算到Y(jié)100。若取783?27.982（5位有效數(shù)字），試問計(jì)算Y100將有多大誤差？解：?Yn?Yn?1?1783100?Y100?Y99?17831001Y99?Y98?7831001Y98?Y97?78310017831001783100……

Y1?Y0?依次代入后，有Y100?Y0?100?即Y100?Y0?783，

若取783?27.982,?Y100?Y0?27.982

1*??(Y100)??(Y0)??(27.982)??10?3

21?Y100的誤差限為?10?3。

29．正方形的邊長大約為了100cm，應(yīng)怎樣測(cè)量才能使其面積誤差不超過1cm？解：正方形的面積函數(shù)為A(x)?x

22??(A*)?2A*??(x*).

當(dāng)x*?100時(shí)，若?(A*)?1,則?(x*)?1?10?22故測(cè)量中邊長誤差限不超過0.005cm時(shí)，才能使其面積誤差不超過1cm10．設(shè)S?212gt，假定g是確鑿的，而對(duì)t的測(cè)量有?0.1秒的誤差，證明當(dāng)t增加時(shí)S的212gt,t?022絕對(duì)誤差增加，而相對(duì)誤差卻減少。解：?S???(S*))?gt??(t*當(dāng)t*增加時(shí)，S*的絕對(duì)誤差增加

?r(S*)??(S*)S*

gt2??(t*)?1*2g(t)2?(t*)?2*t當(dāng)t*增加時(shí)，?(t*)保持不變，則S*的相對(duì)誤差減少。11．序列?yn?滿足遞推關(guān)系yn?10yn?1?1(n=1,2,…),若y0?2?1.41（三位有效數(shù)字），計(jì)算到y(tǒng)10時(shí)誤差有多大？這個(gè)計(jì)算過程穩(wěn)定嗎？

2?1.41

解：?y0?1??(y0*)??10?2

2又?yn?10yn?1?1?y1?10y0?1??(y1*)?10?(y0*)又?y2?10y1?1??(y2*)?10?(y1*)

??(y2*)?102?(y0*)

0??(y10*)?110?y(0*)?10?101?1?022

18??102計(jì)算到y(tǒng)10時(shí)誤差為

其次章插值法

1．當(dāng)x?1,?1,2時(shí)，f(x)?0,?3,4,求f(x)的二次插值多項(xiàng)式。解：

1?108，這個(gè)計(jì)算過程不穩(wěn)定。2x0?1,x1??1,x2?2,f(x0)?0,f(x1)??3,f(x2)?4;l0(x)?l1(x)?l2(x)?(x?x1)(x?x2)1??(x?1)(x?2)

(x0?x1)(x0?x2)2(x?x0)(x?x2)1?(x?1)(x?2)(x1?x0)(x1?x2)6(x?x0)(x?x1)1?(x?1)(x?1)(x2?x0)(x2?x1)3則二次拉格朗日插值多項(xiàng)式為

L2(x)??yklk(x)

k?02??3l0(x)?4l2(x)??(x?1)(x?2)?124(x?1)(x?1)3?5237x?x?6234．設(shè)為互異節(jié)點(diǎn)，求證：（2）證明

?(xj?0nj?x)klj(x)?0(k?0,1,?n,)(2)?(xj?x)klj(x)j?0n??(?Ckjxij(?x)k?i)lj(x)

j?0ni?0iknn??C(?x)(?xijlj(x))k?ii?0j?0n又?0?i?n由上題結(jié)論可知

?xl(x)?x

kjjij?0n?原式??Cki(?x)k?ixii?0n?(x?x)k?0

?得證。

5設(shè)f(x)?C2?a,b?且f(a)?f(b)?0,求證：

1maxf(x)?(b?a)2maxf??(x).a?x?ba?x?b8解：令x0?a,x1?b，以此為插值節(jié)點(diǎn)，則線性插值多項(xiàng)式為

L1(x)?f(x0)=?f(a)x?x1x?x0?f(x1)

x0?x1x?x0x?bx?a?f(b)a?bx?a

又?f(a)?f(b)?0?L1(x)?0插值余項(xiàng)為R(x)