數(shù)值分析課后習題答案

1、第一章題12給定節(jié)點X。=-1,X1=1,X2=3,X3=4,試分別對下列函數(shù)導出拉格朗日插值余項：3(1)(1)f(x)=4x-3x2f(x)=x42X3解(1)f(4)(x)=O,f(4)()由拉格朗日插值余項得f(x)-p(x)=FT(x-Xo)(X-X1)(X-X2)(X-X3)=0f(4)(x)=4!由拉格朗日插值余項得4!f(x)-p(x)=(x-Xo)(x-x1)(xXz)(x-X3)4!=(x1)(x-1)(x-3)(x-4)p(x),插題15證明：對于f(x)以X。，X1為節(jié)點的一次插值多項式2(X1X。)“f(X)-p(x)<maxf(X)值誤差8x。坐坐1f()證由

2、拉格朗日插值余項得f(X)-p(X)=.(X-X0)(X-X1),其中X。'-wX1,f(X)p(x)f()2!maxf<X。名感(x)(X-x0)(X-X1)2七(X1X。)”<maxf(x)8x。a*題22采用下列方法構造滿足條件p(0)=p。)=。，p(1)=p'(1)=1的插值多項式p(x)：(1) (1)用待定系數(shù)法；(2) (2)利用承襲性，先考察插值條件p(。)=p'(0)=。，p(1)=1的插值多項式p(x).23解(1)有四個插值條件，故設p(x)=a。+a1X+a2X+a3X,a0=0a0+a1+a2+a3=1a1=02p(x)=2+2a

3、2x+3a3X,代入得方程組ai+2a2+3a3=1a。=0ai=0a2=2解之，得a3=-123p(x)=2x-x-(2)先求滿足插值條件p(0)=p=0,p(i)=i的插值多項式p(x),由0為二重零點，可設p(x)=ax2,代入p(1)=1,得a=1,j.p(x)=x2；再求滿足插值條件p(0)=p'(0)=0,p(1)=p，(1)=1的插值多項式p(x),可設p(x)=x2+bx2(x1)”p'(x)=2x+2bx(x-1)+bx2,代入p'(1)=1,得b=1,2223p(x)=x-x(x-1)=2x-x332x+x0<x<1S(x)=題33設分段

4、多項式、2x3+bx2+cx_11Wx2是以0,1,2為節(jié)點的三次樣條函數(shù)，試確定系數(shù)b,c的值.解由S(1)=2得2+b+c1=2,b+c=1;-23x2x0：:x：二1S(x)=2(6x+2bx+c1cx<2,由s'(1)=5得6+2b+c=5,j.2b+c=-1;聯(lián)立兩方程，得b=-2,c=3,6x20：:x：:1S(x)：且此時l12x+2bex",S?)=8=S*1),S(x)是以0,1,2為節(jié)點的三次樣條函數(shù).2x+4y=113x-5y=3x+2y=6題35用最小二乘法解下列超定方程組：、2x+y=7.解記殘差的平方和為2_2_2_2f(x,y)=(2x4y

5、-11),(3x5y3)(x,2y-6)(2xy-7)830x二273113y二91ffn二0:x二f36x_6y_102=0=0令，得6x+92y_96=0,解之得題37用最小二乘法求形如y=a+bx2的多項式，使與下列數(shù)據(jù)相擬合：x1925313844y19.032.349.073.397.82解擬合曲線中的基函數(shù)為九=1,中0(x)=x,Y。90)(90,*)YaJ(fW。)'其法方程組為©鼻)(中0$)八b廣9色"其中(：0)=5,(k：；)=(丁，；0)=5327,(；*)=7277699,(f,1)=271.4532a=0.9726547285b=0.0

6、52(f,*)=369321.5,解N得I5696,J.y=0.9726+0.05x.第二章題3確定下列求積公式中的待定參數(shù)，使其代數(shù)精度盡量地高,并指明求積公式所具有的代數(shù)精度：(2)1f(x)dx0：A0f(-)-A1f(-)-A2f(2)從結(jié)論“在機械求積公式中，代數(shù)精度最高的是插值型的求積公式”出發(fā)，1A0=il()(x)dx-011、，13、3(一)(一-)34244,(x1)(x-3-)414dx一一1113、3(一)(-)32424,-01131(x)(x-)12422dx二一A1=I1(x)dx=0-01A2=il2(x)dx=0031313(,)()444211123()-f

7、()f()43234,(x-)(x-)214dx12f(x)dx:.f.03當f(x)=x3時,、，3,，八rf(x)dx二xdx左邊二00右邊二2111-f()一一f()3f(-)4113-3233(一)4左邊=右邊，當f(x)=x4時,f(x)dx=左邊二004.xdx二一211123右邊=3f(；)-；f(2)3f(4)11-()32占(4)37左邊？右邊，所以該求積公式的代數(shù)精度為題8已知數(shù)據(jù)表3.x1.11.31.5xe3.00423.66934.48171.5x-1.1試分別用辛甫生法與復化梯形法計算積分edx辛甫生法1.51.5-1.1exdx3.004243.66934.481

8、7)=1.477541.16復化梯形法1.50.2exdx93.00421.1223.66934.4817)=1.48245_14題17用三點高斯公式求下列積分值7T=7dx-、一一x=一(tI)先做變量代換，設')4.7127dx1-(t1)1+x=45+x94I=3.141068第三章用歐拉方法求解初值問題y=ax+by(0)=0:(1)試導出近似解yn的顯式表達式;解(1)其顯示的Euler格式為：Vn=yn4hf(XnLVn)=Vnh(aXnJ'b)故yn='n2h(aXnN'b)yi=y0h(ax。-b)將上組式子左右累加，得yn=y。-ah(x0,x

9、。/,x。)nhb=ah(0-h2h-(n2)h(n1)h)nhb2=ahn(n-1)/2nhb題10選取參數(shù)p、q,使下列差分格式具有二階精度:yn1=ynhKiKi=f(Xnph,ynqhKJ解將K1在點(XnTn)處作一次泰勒展開，得-2K1=f(Xnph,ynqhK1)二f(4,yn)phfx(%,y0)qhK1fy(%,y0)O(h)二f(Xn，Yn)'Phfx(Xn，yn)_2_2qhf(Xn，yn),phfx(Xn，yn),qhK1fy(Xn，yn),O(h)fy(Xn，yn),O(h)2=f(Xn,yn)phfx(Xn,yn)qhf(Xn,yn)fy(Xn,yn)O(h

10、)._._._._.一.2代入,得yn1-ynhf(Xn,yn),phfx(Xn,yn),qhf(Xn,yn)fy(Xn,yn),O(h).2.2._3yn1=ynhf(Xn，yn)phfxX，丫0)qhf(Xn,yn)fy(Xn，yn)O(h)2h3y(Xn1)=y(Xnh)=y(Xn)hy(Xn)y(Xn)O(h)2h23二y(Xn)hf(Xn,y(Xn)fx(Xn,y(Xn)f(Xn,y(Xn)fy(Xn,y(Xn)O(h)1q=2時,1考慮其局部截斷誤差，設yn=y(Xn),比較上兩式，當p=Z,-3y(Xn1)7n1=O(h)第四章1xcosx題2證明方程-2有且僅有一實根；試確定這

11、樣的區(qū)間a，b,1使迭代過程Xk,=28sXk對一切xoWa,b均收斂.1解設"'：'/0'、，則f(x)在區(qū)間(q*)上連續(xù),1 1二二1二二f(0)=_cos0=_：-0f()=-cos=02 2,22222,Tl所以“*)在0,/上至少有一根；1二f(x)=1sinx.00,-又2,所以f(x)單調(diào)遞增，故f(x)在2上僅有一根.11迭代過程*0s"，其迭代函數(shù)為9(x)=2C0s*,111.-x0,0<9(x)=-cosx_.9(x)0,2,222,2；9(x)-1sinx2,19(x)<-<12冗由壓縮映像原理知一,。0,

12、了，1x.1ccosx.kIk2兀0,一汪這里取a,b為區(qū)間2，也可取a,b為區(qū)間(3*)等.2x,=4cosx題5考祭求解方程123x+2cosx=0的迭代法3(1) (1)證明它對于任意初值x0均收斂；(2) 證明它具有線性收斂性；29(x)=4cosx(1)迭代函數(shù)為3VxW(-,F),9(x)w(g,")；-sinx<<12由壓縮映像原理知以0,xk'=4'3cosxk均收斂;kxk_x-sinx1,二0*3(否則，若sinx*=0,貝Ijxmm,Z三_2不滿足方程)，所以迭代九十二43cosxk具有線性收斂速度;題7求方程x3x2代過程在區(qū)間1,

13、3,1,6_1=0在x0=1,5附近的一個根，證明下列兩種迭上均收斂：解(1)1-1卜一2xk；(2)31-x2其迭代函數(shù)為VxW1,3,1,6,.g(x)-1,3,1,6改寫方程為改寫方程為1g(x)=1'x1.3<1,3906：-1,11-1x,相應的迭代公式為,相應的迭代公式為迭代公式為xk”11.32：1,5917：1,62g(x)二一二又x3,-0.9103=21.33-2<一3x2_-=-0,48831,63，g(x)<0.9103由大范圍收斂定理知一x01,3,1,6xk1=114xk均收斂;x3-x-1=0：=x=1迭代公式為xk.131-x2其迭代函

14、數(shù)為g(x)+X2VxW1,3,1,6,g(x)1,3,1,61,3<1,390831-1,32<31-x2W11+"2江5269<1,6，g(x)一又332x2x0_3"33.(1x)2x33(2x)2=*2_'0,4912g(x)<0,4912<1,由大范圍收斂定理知-x1,3,1,6Xk+M+Xi均收斂.題5分別用雅可比迭代與高斯-塞德爾迭代求解下列方程組:x15x2i3x35x1-2x22x1x2-5x3二-11(k1)x1(k)(k)二2-5x23x3(k1)x2k1)(2)其雅可比迭代格式為5-2-x12112,=X55,(k)1(k)1+21十一5(k)x3(k)*2取初始向量（0）x,迭代發(fā)散;(k1)x1_(k)(k)=2-5x23x3(k1)x2-2其高斯-塞德爾迭代格式為(k1)x