数值分析第五版_李庆扬_王能超_易大义主编课后习题答案.pdf
《数值分析第五版_李庆扬_王能超_易大义主编课后习题答案.pdf》由会员分享,可在线阅读,更多相关《数值分析第五版_李庆扬_王能超_易大义主编课后习题答案.pdf(65页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第 一 章 绪 论1.设x O,x的 相 对 误 差 为 求Inx的误差。p*Y *一 丫解:近似值父的相对误差为5=e;=-%*%*而Inx 的误差为 e(lnx*)=lnx*-lnx进而有e(lnx*)x b2.设x的相对误差为2%,求/的相对误差。解:设y(x)=x ,则函数的条件数为c =i卫 卫I*/(%)又,:/(x)=nx,C=1 A 1=nn又:r(x*)Cp-(x*)且4(x*)为2J(x*)a0.023.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:X:=1.1021,x;=0.031,x;=385.6,x;=56.430
2、,x;=7 x 1.0.解:x;=1.1021是五位有效数字;x;=0.031是二位有效数字;x;=385.6是四位有效数字;x;=56.430是五位有效数字;x;=7x1.0.是二位有效数字。4.利用公式(2.3)求下列各近似值的误差限:x:+x;+x;,(2)x:x;x;,(3)x2lx4.其中x:,x;,x;,x;均为第3题所给的数。解:(X:)=;x l(r4(x;)=gx l()-3(X;)=gx l()T(X;)=gx l O-3(x;)=gx l()T(1)(X:+X;+X:)=(X;)+(X;)+(X:)=-x l O-4+-x l O-*3+-x l O-35 6.4 3 0
3、 x 5 6.4 3 0=1 0一55计算球体积要使相对误差限为1,问度量半径R时允许的相对误差限是多少?解:球体体积为丫=一万/?33则何种函数的条件数为 R V R 4球 、C,=-=:-=3V ,叱3,(/?*)=3 j(R*)2 2 2=1.0 5 x 1 0-3(2)(X:X;X;)=忖司 (X;)+|x;X;|(X:)+忖 X;,(X;)=|1.1 0 2 1 x 0.0 3 1|x l x l 0-1+|0.0 3 1 x 3 8 5.6|x l x l 0_4+|1.1 0 2 1 x 3 8 5.6|x i x l 0-3 0.2 1 5(3)(X;/X;)同 (X;)+|x
4、 (X;)i ,i ,0.0 3 1 x-x l 0-3+5 6.4 3 0 x-x l O 32又,*)=1故度量半径R时允许的相对误差限为j(R*)=;x l。0.336.设=2 8,按递推公式匕=工1 -V783(n=l,2,.)计算到九加。若取J旃a 27.982(5位有效数字),试问计算地将有多大误差?解:.%=%+屈Y oo=%-Y gg K)g-V783100L9 8%9 7 -J78310()-V7830 100依次代入后,有K0G=%IOOX击即几0=乂-,若取27.982,./0G=%-27.982(X M)=(/)+(27.982)=1xl0-3.y1 二1Go的误差限为
5、5x10-3。7.求 方 程 炉-56x+l=0的两个根,使它至少具有4位有效数字(J旃=27.982)。解:X2-56X+1 =0,故方程的根应为4 2 =28 土 7783故 玉=28+7 2 8 +27.982=55.982.王具有5位有效数字X,=28-V783=x-=-0.01786328+V783 28+27.982 55.982%具有5位有效数字阴+1 18.当N充分大时,怎 样 求:-dx2卜 +X2rN+1解 Jv-2 dx=ar ctan(N +1)-ar ctan N设 a =ar ctan(N +1),=ar ctan N。则 tan a =N+1,tan B =N.J
6、w+iN11 +x2dx-a-3-ar ctan(tan(a 一 夕)ar ctantan cr -tan /?1 +tan ez tan (3ar ctanN+l-N1 +(N +1)N=ar ctanM+N+I9.正方形的边长大约为了 1 0 0 cm,应怎样测量才能使其面积误差不超过1 cm2?解:正方形的面积函数为A(X)=2(A*)=2 A*(尤*).当 x*=1 0 0 时,若 (4*)4 1,则(X*)LX1 0-22故测量中边长误差限不超过0.0 0 5 cm 时,才能使其面积误差不超过l cJ1 0.设S=;g/,假定g 是准确的,而对t的测量有0.1 秒的误差,证明当t增加
7、时S的绝对误差增加,而相对误差却减少。解:5 =洌1 ,/()(S*)=gt2 (f*)当f*增加时,S*的绝对误差增加?(5*)=芾=gr g(r*)*)2_ 2 /*)当f*增加时,“*)保持不变,则 S*的相对误差减少。1 1.序列 为 满足递推关系先=1 0 先_1 一1 (n=l,2,.),若治=及=1.4 1 (三位有效数字),计算到凹。时误差有多大?这个计算过程稳定吗?解:0=血;1.4 11 ,,(%*)=于 1 0-2又,.”=1 0 几-1(%*)=1。(%*)又:y2=1 0%T(%*)=1。(/*).,.(%*)=1。2 (%*)(w*)=1。1 (%*)=10l(,x
8、 lx l0-22=-x l 082计算到M o 时误差为:X 10 8,这个计算过程不稳定。12.计算/=(血 一 1)6,取 血“1.4,利用下列等式计算,哪一个得到的结果最好?-r1 ,(3-2V2)3,9 9-7 0 0。(V2 +1)6(3+2 V 2)3解:设 y=(x 1卜,若=&,x*=1.4,则式x*)=g x l()T。若通过一计算y 值,则(V2+1)6*1()=6 x;-(x +1)(丁)6 */*、=j 7 y s(x )(X+i y=2.5 3yb(x*)若通过(3-2 夜)3计算y 值,则8(y*)=|-3 x 2 x (3-2x*)21 8(x*)6 */*j
9、*ya 工)3-2x=30 y*(x )若通过计算y 值,则(3+23(*)=3 x-(x*)(3+2x y=6x-r z-y*(x*)(3+2.r)7=1.0 34 5),*(尤*)通过二 二计算后得到的结果最好。(3+2 V 2)313./(x)=l n(x Jx 2-l),求。(30)的值。若开平方用6位函数表,问求对数时误差有多大?若改用另一等价公式。l n(x _ Jx 2_ l)=Tn(x +Jx 2_ l)计算,求对数时误差有多大?解 /(x)=l n(x-7x2-l),.-./(30)=l n(30-/899)设 =屈 =/(30)则“*=29.9833(*)=;x l 0-4
10、故底卜一岛而)=焉而)。3x 10-3若改用等价公式n(x-yl x2-1)=-l n(x +/x2-1)贝 i J/(30)=-l n(30 +屈此时,(),*)一二 J (*)30 +=-E(M)5 9.9833a 8x 10-7第二章插值法1.当x =l,-1,2时,/(x)=0,-3,4,求/(x)的二次插值多项式。解:=1,4 =1,%2=2,/(%)=0,/(%)=-3,/()=4;/(%)=_ l(x +i)(x_ 2)(%一%)(%一 工2)2/,(x)=(X 虫一玉)=-(x-l)(x -2)(七一 Xo)(X|一)6(X-Xo)(X-X|)1L(x)=-!=(x l)(x
11、+1)(乙一公)(一玉)3则二次拉格朗日插值多项式为2k=O=-3Z0(X)+4/2(X)1 4=-(x-l)(x-2)+-(x-l)(x +l)2.给出/(x)=l n x的数值表X0.40.50.60.70.8Inx-0.916291-0.693147-0.510826-0.356675-0.223144用线性插值及二次插值计算In 0.5 4 的近似值。解:由表格知,XQ 0.4,%)0.5,x,=0.6,X3 0.7,0.8;/(x0)=-0.916 291J (x j =-0.6 9314 7/(x2)=0.5 10 826,/=0.35 6 6 75/(x4)=-0.22314 4
12、若采用线性插值法计算In 0.5 4 即/(0.5 4),则 0.5 0.5 4 0.6=-=-10(x-0.6)当一/2(x)=-10(x-0.5)2-占“X)=/(X|K(x)+/()4 (x)=6.9314 7(%-0.6)-5.10 826(x -0.5)L,(0.5 4)=-0.6 20 2186 -0.6 20 219若采用二次插值法计算In 0.5 4 时,,o(x)4(x)(X-X)。-/)(X0-X,)(X0-J C2)(x-x0)(x-x2)(七一%)(占一)(X_ Xo)(X_ X|)(x2-x0)(x2-x,)=5 0(x-0.5)(x-0.6)-10 0(x-0.4)
13、(x-0.6)=5 0(x-0.4)(x-0.5)4(%)=/(x()为(%)+/*1R (x)+/(4 M(%)=-5 0 X 0.916 29 l(x -0.5)(x -0.6)+6 9.314 7(x -0.4)(x -0.6)-0.5 10 826 x 5 0(x -0.4)(x -0.5)4(0.5 4)=-0.6 15 31984 -0.6 15 3203.给全c o s x,0 Y x M 9(T 的函数表,步长 =1 =(1/6 0),若函数表具有5 位有效数字,研究用线性插值求c o s x近似值时的总误差界。解:求解C OSX近似值时,误差可以分为两个部分,一方面,X 是近
14、似值,具 有 5位有效数字,在此后的计算过程中产生定的误差传播;另一方面,利用插值法求函数C OSX 的近似值时.,采用的线性插值法插值余项不为0,也会有一定的误差。因此,总误差界的计算应综合以上两方面的因素。当0 K x 90 时,令/(x)=c o s x取/川=白磊畸=品令若=x04-i/i,z =0,1,.,5 4 0 07T则5 =5=90“当x e x*,x j 时,线性插值多项式为4(x)=/(x j x_%+/(%)-*一心4一次+1 4+i x*插值余项为R(x)=|c o s x -A,(x)|=;/)(x -)(x x+1)又.在建立函数表时,表中数据具有5位有效数字,且
15、c o s x e 0,1,故计算中有误差传播过程。.(/*&-)小炉R,(x)=(f xk)-+(/*(1)(/*(X*)(Xf L+Xf)Xk Xk+Xk+Xk=(/*(4 )(小 _ X+X _ X。h=(/*(“)总误差界为R=/?|(x)+7?2()=g(cosj)(x xj(x 4+1)+(/6)|x(x-xt)(xt+|-x)+(f xk);力)2+(/*1)=1.06xl0-8+-x l0-52=0.50106x10 54.设为互异节点,求证:(1)三3 (k=0,1,/!);J=0JJ(2)(Xj-x)kIj(x)=0(k=0,1,);;=o证明(1)令/(.)=%若插值节点
16、为Xj,j=0,,则函数/(x)的n次插值多项式为L.(x)=。j=o插值余项为此=-,仆)=需%又,:k 4.尸)(0.R.(x)=0Xjlj(x)=xk(女=0,篦);y=o f (xL)六0=(c W(r)J%(x)y=o z=o=力。5产 之 叽(x)/=0 j=0又04注 由上题结论可知Z x;/j(x)X1j=0原 式=c;(-x)y1=0=(x x)”=0得证。5 设/(x)G C2 a,h且/(a)=f(b)=0,求证:max|/(x)|(Z?-a)2max|/ff(x)|.axb 8 axb解:令4=a,X1=b,以此为插值节点,则线性插值多项式为4(x)=/(%)X-X|+
17、/(%)xQ-x x-xQ“、x-b”,、x-a=;()-+/0)a-b x-a又/(a)=/(b)=0.-.Ll(x)=0插值余项为 R(x)=/(x)-L(x)=;r(x)(x-x()(x X)/(X)=g/(x)(x-Xo)(x X 1)又|(x-X o)(x-X|)|w;(x-X o)+a -x)=-(X I-x0)2=;(j)2max|/(x)|-a)2 max|/ff(x)|.axb 11 8 幼6.在-4 V x4上给出/(x)=e,的等距节点函数表,若用二次插值求e*的近似值,要使截断误差不超过10M,问使用函数表的步长h应取多少?解:若插值节点为王_”王和毛+|,则分段二次插
18、值多项式的插值余项为&(x)=不/e)(x )(x%)(x -X,.+1)火2 )|,(x -X-)(x -X,)(x -X,.+1)ma x|/w(x)|设步长为h,即=x h,xM=xi+h:.Rx)-e4-h3 e4h1-1 6 3G 2 7若截断误差不超过1 0 ,则|2W|i o-6.-.eV 1 02 7/.h 0.006 5.7.若”=2 ,求 心 及b”,解:根据向前差分算子和中心差分算子的定义进行求解。K =2 ”=(-炉 先=(2-D4Z,=2 1 _ 1色.=(E2 _ E2)4%I=(E)4(-D4X,=E-2 yn yn-2rn-28.如 果/(x)是m次多项式,记V
19、(x)=/(x +/z)/(x),证 明/(x)的k阶差分/(x)(OWA W?)是m k次多项式,并且AM+i/OOnO(/为正整数)。解:函数/(x)的2”展式为f)=x)+r 3 仆诉+)+岛/吸加“其中 Je(x,x+/?)又/(x)是次数为机的多项式 尸)=。Af(x)=/(X +/i)-y(x)=f M h+!_T(x)2 +-+f(m)(x)A2 m l.g(x)为2 1阶多项式A2/(x)=A(AfW)乂/口)为 机-2阶多项式依此过程递推,得K/(x)是机-左次多项式A(x)是常数.当/为正整数时,!(x)=O9.证明(&)=力公W+取+/证明M f kgk)=f k+igk
20、+f kgk=f k+SM-f k8k+f k8k+f k8k=SM(A+I-f j +f k(Sk+i 8k)=g z/+f i A g k=.%+gnM得证,一!10.证明=f 8-f g o-E g k+Mk=0k=0证明:由上题结论可知力 坳=(4)-gnMk=0n-1=Z(/)-g +i 颂)k=0=(力gJ-g*+Mk=0 k=0(f kgk)=f k+l gk+l-f kgk;$(&)k=0=(/g -/o g o)+(&g 2 一/g)+(力g 一 l g-l)=f n8n f()8o,一!E f Z k=f“g n-f o g o-E g k+N kk=0J c=O得证。i
21、i .证 明 工 1 力=,一%J=o证明 A2y.=(5-旬)j=0 j=0=(%)+(%-A X)+(”一 AX.-i)=/“一%得证。1 2 .若/(x)=al)+alx+-+al l_lxl+a x 有”个不同实根%,马,x”,白 x*(0,0k n-2;证明:,占/(Xj)。4 =-1证明:、f(x)有个不同实根国,尤2,X且 f(x)=a0+alx+-+a“_ i x i +anxn./(%)=an(x-x,)(x-x2)-(x-x)令 七(x)=(x 一 再)(x 一 )(x x“)n k kYK n Y则白 r(Xj)而 式。)=。_ 苫2)(一曰一。一苍,)+。一芯)。-%3
22、)一。一乙)+(%-%)(X-X2)-(X-Xn_ 1)”Xj)=(Xj-X|)(Xj _ 4)(Xj -Xj T)(Xj -Xj+J(Xj -x“)令 g(x)=x*,n xkg X,X2,=则 g X,X2,x“=fW Xk J又。下 入=g X,X2,X,1月/(Xj)a“9 x:Q,Ok (Xj -x)-c(t _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _)=0(Xj Xo)(Xj -Xy_ 1)(x.-xj+i)(勺-x)=(/%,王,x“得证。./(x)=/(x)+g(x)
23、y _ _ _ _ _ _ _ _ _ _ _ _ _ _ _/(x,)+g(W)_ _ _ _ _ _ _ _ _ _ _ _ _ _ _占(Xj 尤o )一(勺-Xj _ 1 )(Xj Xj+I)(勺 x“)y_ 0_)白(Xj Xo)(Xj -Xj _|)(Xj _ X*)(Xj -x)+y-2-)尸)(x,-x0)-(xy-Xj T)(Xj _X,+I)(Xj _x“)=fx0,-,xn+gx0,-,xn.得证。1 4./()=/+/+3+1,求产2 ,2 一,2 7及尸2 ,2、2 8。f t?:v/(x)=x7+x4+3 x +l若 Xj =2 ,i =0,1,-,8心)则/1%,
24、石,x“=n!乎aI了%,玉,闻=勺=01 5.证明两点三次埃尔米特插值余项是R3M =/(4,OX-X,)2(X-X,+1)2/4!,G ,王小)解:若XWX,XJ,且插值多项式满足条件/)=尤。既/)=广 )插值余项为R(x)=/(x)-3(x)由插值条件可知R(4)=R(k J =0且 RU)=R(4M)=0R(x)可写成R(x)=g(x)(x xj2(x x*+I)2其中g(x)是关于X的待定函数,现把x看成 XA,X/+J上的一个固定点,作函数阿=%)-%-g(x)-xk)t-xk+l)2根据余项性质,有夕心)=0,夕(+)=09(x)=/(x)-%(x)-g(x)(x-xJ 2%=
25、/(X)3(X)-R(X)=0(p(t)=(t)-一 g(x)2(f-x,)(Z-x,+1)2+2(Z-xt+1)(/-x,)2(pxk)=0由罗尔定理可知,存在Je(x,x)和Je(x,x+),使”&)=0 2)=0即夕(X)在X,X/+J上有四个互异零点。根据罗尔定理,夕 。)在(f)的两个零点间至少有一个零点,故夕(f)在(,xM)内至少有三个互异零点,依此类推,“4)在(X*,七华)内至少有一个零点。记为 qe(x,4+1)使=-4)c)_34)4!g(x)=0又 =0f(4).g(x)=/,会 人 加)其中J依赖于XR(x)/(94!(x-XjJ-Cx xt+1)2分段三次埃尔米特插
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 第五 李庆扬 王能超 大义 主编 课后 习题 答案
限制150内