数值分析第四版课后习题答案.pdf
《数值分析第四版课后习题答案.pdf》由会员分享,可在线阅读,更多相关《数值分析第四版课后习题答案.pdf(96页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第一章1、设x 0,x的相对误差为3,求I nx的误差。解 设x*0为x的近似值,则有相对误差为;(x)=6,绝对误差为*(x)=廉”,从而 I nx 的误差为 (I nx)=|(l nx*“(x)=二 加=S,相对误差为;(I n x)=皿*=-AI nx I nx2、设x的相对误差为2除 求x 的相对误差。解 设x*为x的近似值,则有相对误差为;(%)=2%,绝对误差为*(x)=2%卜,从而 x 的误差为 *(l nx)=,(x*)=2 ,相对误差为 ;(l nx)=;:|;j)=2。3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是儿位有效数字:x;=
2、1.1 0 21 ,x;=0.0 31,x;=38 5.6,x;=5 6.430,x;=7 x 1.0。解x;=1.1 0 21有5位有效数字;x;=0.0 0 31有2位有效数字;x;=38 5.6有4位有效数字;x;=5 6.430有5位有效数字;x;=7 x 1.0有2位有效数字。4、利用公式(3.3)求下列各近似值的误差限,其中x:,x;,无;,x;均为第3题所给的数。/1 *辛 辛(1 )匹 +%+%;e*(x:+X;+*;)=萼 卜(X;)=()+(X;)+(X;)解 .=-x l 0-4+-X 1 0-3 4-l x l O _3=1.0 5 x l 0-32 2 2(2)x;x
3、;x;*/*、d f/*、/*、/*、/*、/*、/*、,*、e(七4 乙)二Z 丁 (%)=(X2X3)(演)+(x/3)(Z)+(x/2)()k=7 (0.031x385.6)-xl0-4+(1.1O21X385.6)-X1O-3+(1.1O21XO.O31)-X1O-3;2 2 2=0.59768 X10-3+21 2.48488 xl0-3+0.01708255 x 10-3=213.09964255 xlO-J=0.21309964255(3)%;/x;。e*(%;/%;)=k=(x:)=r (芯)+?)(%;)乙(乙)解 八,x L l k+O 3 l,xLl。56.430 2(5
4、6.430)2 2*4 6 1 x_LxI。、=0 88654x107(56.430)2 25、计算球体积要使相对误差限为现,问度量半径R允许的相对误差是多少?4*(1%(/?*)3)解 由 1%=;(兀(R*)3)=-可知,3 加*)3*(g(R*)3)=l%x g%(R*)3=:%(R*)3 *(/?*)=4乃(R*)2X*(R*),从而 (R)=1%X-R,故 ;(R*)=*f)=1%X =-o3 R 3 3006、设 公=2 8,按递推公式匕=匕1-高(=1,2,-一)计算到匕0 0,若取V783 27.982(五位有效数字,)试问计算匕皿将有多大误差?解 令匕表示,的近似值,e*(%
5、)=匕一匕,则e*(y0)=0,并且由匕=匕 一 看 X27.982,/=一 击*历 可 知,Yn-Yn=;_I_yn_1_-l_x(27.982-V 783),即e*(y)=e*(%T)-焉 X(27.982-7783)=e*(Yn_2)一 高 x(27.982-A/783)=,从W/(yiOo)=*()-(27.982-7783)=7783-27.982,而-27.9 8 2 “)=工-打,则由卜。=及 与J”=1。-Tv =14 1 1 可知,*(%)=彳 Xi。-工-%=1 0(-h1),即口=10九-1-1 2*(%)=10*(y,T)=10*(y。),从而*(乃0)=101*(汽)
6、=10隈3 1 0-2=91 0 8,因此计算过程不稳定。12、计算/=(&-1)6,取0=1.4,利用下列公式计算,哪一个得到的结果最好?(V2 +1)6(3-2扬 3-L,9 9-7 07 2 o(3+2 V2)3解因为 *(/)=;x l()T ,所以对于 =(V2 +1)6e*(/;)=M|e*(L4)=x,xl OT=6.5 4X1(T 1*1()-2,有一位有效数字I 1 I (1.4 +1)7 2 2对于=G-2/儿(/2)=|/2y(1.4)=6(3-2 x 1.4)2x 1 x 10-1=0.12 x 10-1 1 x l 0-1,没有有效数字;对于八=(3+2扬 3e*()
7、=L|e*(L4)=-Tx-x l O-l=2.6 5 X10-3 -x l0-2,有一位有效数3 T I (3+2 x 14)4 2 2字;对于=99 7 0匹,e*()=e(L4)=7 0 x g x l()T =35X1()T g x l()i ,没有有效数字。13、/(x)=ln(x-A/X2-1).求/(30)的值。若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式i m x-T 7=Tx-ina+TT%)计算,求对数时误差有多大?解 因为,3()2 一 1=次 旃=2 9.9 8 33(六位有效数字),f*(x)=-X10-4,所以2e*(/i)=(/:)*e*(x)=
8、-1t-x-x l O(3 0-V302-1)2-1-x-x l O-4=0.2 9 9 4 x 10-230-2 9.9 8 33 2e*(32)=|(/2 )e*(x)=-x +x-x l O-42-1-x 1x l()T =0.8 336 x l()F30+2 9.9 8 33 214、试用消元法解方程组X+i n10 r =i n10 2,假定只有三位数计算,问结果是否X1+工2 =2可靠?解 精 确 解 为X =1 A1 0 1 A10 _ 7T ,2=二。当 使 用 三 位 数 运 算 时,得到10-1 2 10,-1xx=l,x2=1 ,结果可靠。15、已知三角形面积,=,出!1
9、。,其中c为弧度,0 c s i n c|A a|dXk 21 a2+s i n c c o s c|A c|,;b s i n c +1一 s i n e2+-abcosc2所以1 absincA c+Nbb+A ct a n c n(x-Xj)证明 由 =:.j=可得求证。n-l=V,T(Xo,X D X,L i n(x-X j)j=o2、当x =l,-l,2时,/(x)=0-3,4,求/(x)的二次插值多项式。(x-x)(x-x2)(x-x0)(x-x,)(X-Xo)(X-X|)L,(x)=y0-!-=-+y,-=-+y,-(x0-x,)(x0-X2)(X|Xo)(X|一/)(x2-x0
10、)(x2-X1)r*-i (x+l)(x 2)(x l)(x 2)(x l)(x +1)解=0 x -+(-3)x -+4 x -(1+1)(1-2)(-1-1)(-1-2)(2-1)(2 +1)=-(%2-3x +2)+(x2-1)-x +x 2 3 6 2 33、给出f M=In x的数值表用线性插值及二次插值计算In 0.5 4的近似值。X0.40.50.60.70.8In x-0.9 16 2 9 1-0.6 9 314 7-0.5 108 2 6-0.35 7 7 6 5-0.2 2 314 4 解 若取 为=05,%,=0.6 ,则%=/(/)=/(。5)=-0.6 9 314 7
11、,y,=/(%,)=/(0.6)=-0.5 108 2 6,则L*.(/x)、=yx x.x xn,x 0.6 /s ee/x _ 0.50-+y.-=-0.6 9 314 7 x-0.5 108 2 6 x-1 Xo-/-1 x,-x0 0.5-0.6 0.6-0.5 ,=6.9 314 7(x-0.6)-5.108 2 6(x -0.5)=1.8 2 32 lx -1.6 04 7 5 2从而 L,(0.5 4)=1.8 2 32 lx 0.5 4-1.6 04 7 5 2 =0.9 8 4 5 334-1.6 04 7 5 2 =-0.6 2 02 18 6。若取3 =0.4 ,x,=0
12、.5 ,x2=0.6 ,则 氏=f(x0)=/(0.4)=-0.9 16 2 9 1,%=/(/)=/(0.5)=-0.6 9 314 7,y2=/(x2)=/(0.6)=-0.5 108 2 6 ,则,、(j-x.X x-x J (x-x0)(x-x2)(%-%0)(%-%)L,(x)=y()+y2(x0-X)(x0-x2)(x,-%0)(%,-x2)(x2-x0)(x2-x j=-0.9 16 2 9 lx (x-5Xx-6)+(-0,6 9 314 7)x(0.4-0.5)(0.4-0.6)(0.5-0.4)(0.5-0.6)+(-0,5 108 2 6)x *一 04)(x-0-5)(
13、0.6-0.4)(0.6-0.5)=-4 5.8 14 5 5 x(x2-1.lx +0.3)+6 9.314 7 x (x2-x +0.2 4)-2 5.5 4 13(r -0.9 x +0.2)=-2.04 115 x 2 +4.06 8 4 7 5 x -2.2 17 09 7从而&(0-5 4)=-2.04 115 x 0.5 4 2 +4.06 8 4 7 5 x 0.5 4 -2.2 17 09 7=-0.5 9 5 19 9 34 +2.19 6 9 7 6 5-2.2 17 09 7 =-0.6 15 319 8 44、给出c o s x,0 Wx 9 0的函数表,步长力=1
14、=(1/6 0),若函数具有5位有效数字,研究用线性插值求c o s x近似值时的总误差界。解 设插值节点为/X X =/+/?,对应的C O S X值为先,必,函数表值为了0,不 则由题意可知,瓦焉|()一 歹0)+(y i 2!x0.X x-x0C O S J/、/、/_ x x-x./_ x X-XG 7 X=一 一-(x-x0)(x-x,)+(y0-y0)-+(%一%)-,穴(x(),x J2%一七%,-x0,从而|R(X)K 1 c o s J(X-x0)(x -X|)|+瓦一 y0|A _A +一 y jx -x0 (X-Xg)(x_ X)+x l O XX-X,1 5 x-x。-
15、+-X 1 0 x-人 一 匹 21 h2-2 4+-X 1 0-52=-x+-X10-5=-x 6.9 4 x l 0-5+-X 1 0-5=3.4 7x 1 O-52 14 4 00 2225、设占=%()+劭,k=1,2,3,求 m a x,2(x)|x0-x-x2m a x|/2(x)|.VQA*(X-Xo)(X-X|)(X-X3)m a x-!-:%(x2-x0)(x2-Xj)(x2-/)解=m a x(x x0)(x x0 h)(x x0 3 h)(2/2)A(-/2)=y m a x (x-x0)(x-x0-=(x -x0-3/i)|f(x)=(x-x0X x-xQ-h)(x-x
16、Q-3A)令,则=x (3XQ+4/?)x +(3x j +8/?+3 h )x (x(:+4/z x j +3 h )fx)=3x*2-2(3X(,+4h)x+(3x;+8xoh+3 h2),从而极值点可能为6,又因为_(3x0+4/i)V 7/?_ 4 V 7-=X(fl3 34 -y l 4 V 7 1 -1 5 V 7 1 /z-3/(x0 H-h)=-h x-h x-h=(14 V 7 20)/i,3 3 3 3 274 +V 7,、4 +7 7,1 +7 7,7 7 5,1 一 不,3f(xn H-h)=-h x-h x-h=-(20+14 J7)/7,3 3 3 3 27显然/G
17、 o +4,h)=(y-x)*,则左侧是/()=(y-x)k的n阶拉格朗日多项式,令 y =x ,即得求证。7 设/(x)e目 且/=/(b)=0,求证 m a x|/(x)|(。一a)2 ma x|/(x)|。_ 2(3x0+4/)74(3x0+4/?)2-12(3x +8xoh+3 h2)解 见补充题3,其中取/(幻=/3)=0即得。8、在-4W x4上给出/(x)=的等距节点函数表,若用二次插值求e 的近似值,要使截断误差不超过1C T6,问使用函数表的步长h应取多少?解 由题意可知,设x使用节点x 0=x/z,再,/=玉+力进行二次插值,则此(x)=(%-X|)(x -x2)插值余项为
18、/,e s=-x-(xt-72)(X-X)X(X +/),Jw (x0,x2)6令 f M =(-1 _ h)(x-x)x-(x1+h)=x3-3 xx2+(3x;-h2)x +x1-h2),则/V)=3x2-6 X X+(3x:-/z2),从而 f(x)的极值点为 x =xi h,故I /.z J V 3,n V 3 V 3.,2V 3 3 而m a x /(x)=-/?.(14-)h-(1-)/2=-A,而e0 3 3 3 9|/?2(x)|m a x|/(x)|/z3=h3,要使其不超过 10,则有6 “0文 秘2 6 9 27,+2+(T“(J4、x+i+(T)(a4、+(T(口4、小
19、 +(T4、)。=2n+2-4 x 2n+1+6 x 2-4x 2-+2-2=16 x 2-32x 2T +24X2T 一 8x 2一+2 证明:Z x=4 v 一八打。j=OJ LI.?一!证明 X x=Z(加一划)=-)()。j=o j=o14、若 f(x)=a0+a,x +-+an_,xn +%x 有 n 个不同实根x”/,怎,证明 x:JO,0 k n-2马 广(乙)=后,k=n-证明 由题意可设 f(x)=an(x-)(x-x2)(x-xn)=a“什(x-x j ,故(=if,(xj)-anY(xj-%,),再由差商的性质1和 3 可知:/=1沏=3=,工 区,,招=富?,从而得证。
20、“(为7,)明 明(i=M j15、证明n 阶均差有下列性质:1)若 F(x)=t/(x),则“XtpXi,,x“=c/Xo,X,七;2)若 F(x)=/(x)+g(x),则 Fx(),X1,x,=/xo,x”,x,+gxo,x n x,。e F(x)cf(x:)尸 0,芯,x,=-=-网 x f j(X j-巧)f=o/=0 证明 1),由V f(Xj)=CL-=cfx(),xl,-,xn川口-七)i=0Hi(尸(Xj)0 与)2*4+|)2/4!,(4,4+1),并由此求出分段三次埃尔米特插值的误差限。Q 解 见 P30与P 33,误差限为。)+hmax园。27 118、X X X X X
21、 X X X X X.19、求一个次数不 高 于 4 次的多项式P(x),使它满足P(O)=P(O)=O ,F(l)=PQ)=1,P=1 角 星 设 P(x)=24%4+aix+的/+/,贝U P (x)=4 tz4x3+3a3x2+2a2x +at,再由尸(0)=(0)=0,p(l)=p/(l)=l,P(2)=l可得:0=P(0)=劭0=a八00=/(0)=%0=%1 =尸=a4+a3+a2+a1+0 解得 8 时,e.(x)在以上一致收敛到/(x)。sup/(x)+i n f /(x)解 令处(x)=3当 产巴一,i =l,2,3,*2 1、设/*)=1/(1 +/),在 5Wx 5上取”
22、=1 0,按等距节点求分段线性插值函数(x),计算各节点中点处的,(x)与/(x)的值,并估计误差。解 由题意可知,h=,从而当xw x*,x*+J时,/,%(/x)、=fr,1 X X,.+|1 X-Xkklk+fk+L+l =-7-L+-r-l +k-xk-xk+l l +(k+l)-xk+l-xko11=-(X -4+1)*-(x-)h(+k2)/i(i+a+i)22 2、求x)=/在 上 的 分 段 线 性 插 值 函 数/%(x),并估计误差。解 设 将 a,划 分 为 长 度 为h的小区间a =x 0 W X|W W x.=8 ,则当x e x,4+J,%=0,1,2,一1 时,,
23、/一,2 元-4+1 2 XXkIh(%)=f Jk+fk+Jk+=Xk+Xk+XXk+i KfXk+(X 4)一 Xk(X Xk+)_ X(X;+1-年)+4+1 X;Xk+Xk _/.r r Y一8 4+十 4 J 4+1 43-4从而误差为此(X)=萼(f )(x-)=(f-力2故艮(幻|二|(%一九)。一匕+)K彳。2 3、求/。)=/在 鼠”上的分段埃尔米特插值,并估计误差。解 设 将 划 分 为 长 度 为h的小区间。=与 玉工工%=。,则当xe 民,Z+J,k=0,1,2,,-1 时,/力(x)=fk%+fk+l%+fkBk+fk+lBk+1(X)可Ud(f)+4认二:卜一心r(
24、4)(Q从而误差为 R2(x)=-(x-xA.)2(x-xk+)2=(x-xk)2(x-xk+)2,故网=|(X7*)2(X-X*+1)2归 而。2 4、给定数据表如下:Xj0.2 50.300.390.4 50.5 3yj0.5 0000.5 4 7 7 0.6 2 4 5 0.6 7 08 0.7 2 8 0试求三次样条函数S(x),并满足条件:1)5/(0.2 5)=1.0000,5z(0.5 3)=0.6 8 6 8 ;2)S(0.2 5)=S (0.5 3)=0。解 由%=0.30-0.2 5 =0.05 ,4=0.3 9 0.30=0.09,h2=0.4 5-0.39 =0.06,
25、h./?3=0.53-0.45=0.0 8,及(8.1 0)式儿=2 7 hJ,T,+/?,J.hhM,(J=1,-1)可知,40.090.05+0.09 4 =0.062/?1 +h 0.09+0.06 54=0.084h2+h3 0.06+0.08 74=%30.0560+h 0.05+0.095-,214 20.094+0.09+0.06 5力2出=0.063h2+/z3 0.06+0.08 7由(8.1 1)式勺+(j=l,一1)可知,g,=3(A1/X0,XI+/1/X1,X2)=39/(x j-/(x )5/(4)-/(花)14%1-%014 x2-Xjc/9 0.5477-0.5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 第四 课后 习题 答案
限制150内