計算方法考題B07(答案)

精品文檔西安交通大學 研究生考題考試 考查 成績課程 計算方法（B） 考試日期 2007-12-30學院姓名學號1、（4分）浮點數(shù)系F(10,5,2,3)中共有1080001個數(shù)（包括0），實數(shù)3.1415926??和100在該數(shù)系中的浮點化數(shù)fl()=0.31416E1，3fl(100)=0.33333E2，在浮點數(shù)系F(10,5,2,3)中計算3fl()fl(100)0.36475E2；32、（4分）按數(shù)值積分的復化梯形公式計算得：T40.94451,T80.94568。由此可估計誤差：I(f)T81(T8T4)0.00039；33、（4分）設函數(shù)p(x)x3qx2rxs,若p[0,1,t]1,那么p[1,2,t]3；100010004、（4分）矩陣A1/2100，則A11/2100，且A125，1/30101/3010121/40011/4001A3/2，范數(shù)意義下的條件數(shù)Cond(A)=9/4；5、（4分）用四等分法計算f(x)x24x2的極小點，若以[0，8]區(qū)間為初始搜索區(qū)間，那么第一步刪去部分區(qū)間后保留的搜索區(qū)間為：[0,4]；第二步保留的區(qū)間為：.[1,3]；6、（6分）已知如下分段函數(shù)為三次樣條，試求系數(shù)A,B,C：。1歡迎下載精品文檔3Axx12S(x)22xBx21x31x0222xCx2x30x1則A=1/2，B=3/2，C=3/2；7、（8分）求滿足下述插值條件的插值多項式p(x)x-1012y(x)-2-112y(x1)xiyi[,][,,][,,,][,,,,]12011解：建立差商表110011211/22211/23/45/12p(x)2(x1)1(x1)x25(x1)x2(x1)2128、（8分）已知yf(x)的數(shù)據(jù)表如下：xi-2-1012yi00.20.50.81求一次式p(x)axb，使得p(x)為f(x)的最小二乘一次近似；121111111解：G10,GT2101,212502.5GTG10,GTY02.6。2歡迎下載精品文檔GTGb2.5a2.6b0.5,a0.26p(x) 0.26x 0.5。3歡迎下載精品文檔9、（10分）將如下線性方程組的系數(shù)矩陣 A作Crout分解，即分解為矩陣乘積LU形式（L下三角、U單位上三角矩陣），并求解該線性方程組。x12x22x37x13x24x3x4102x27x37x44x12x211x314x4212207134110l211解：l31002774l4111211142

122071213277409145l32 2l42 0

122071220711213l43312132352352x19145011111220112201112111121A213502315/30103111091110、（10分）線性方程組：x1 2x2 2x3 1x1 x2 x3 32x1 2x2 x3 5考察用Jacobi迭代和Gauss-Seidel迭代解此方程組的收斂性；解：12202222A111,J101det(IJ)113(J)022122022又，det(IGS)detIDE1FdetDE1[(DE)F]0等價于解：22det(DE)F13424(2)2(GS)2；22結(jié)論：Jacobi 迭代收斂，Gauss-Seidel迭代不收斂。。4歡迎下載精品文檔11、（10分）方程x3x2110在1.5鄰近有根x，首先討論迭代xk1(xk1)2的收斂性；其次對此迭代實施改善，若不收斂，使改善后的迭代收斂； 若收斂，使改善后的迭代收斂加速； ( 2 1.414, 3 1.732, 5 2.236)解：(x)(x1)1/2,'(x)1/2(x1)3/2，取區(qū)間1.4,1.6，顯然x1.4,1.6。31此時， '(x)'(1.6)(1.61)22

1 1 52 0.6 3

1.076 1, x [1.4,1.6]，因此(x)x(x)(x)'()xxxxx[1.4,1.6]1所以，迭代xk1(xk1)2不收斂；改善：取'(1.5)1(1.51)2

32

1 12 1.4142 0.51(x)x(x1)21.414x，(x)2.4141迭代xk1(xk),x01.5必收斂。12、（14分）試求如下數(shù)值積分公式的結(jié)點xi(i0,1,,n)及求積系數(shù)Ai(i0,1,,n)，使公式具有最高代數(shù)精度，并求其誤差：1n2f(x)dxAif(xi)（n0,1)x1i0解：方法一具有最高代數(shù)精度的求積公式必是 Gauss型求積公式，對 (x) x2,積分區(qū)間(1,1)，先求相應的正交多項式。內(nèi)積：1x2f(x)g(x)dx(f,g)10(x)1(x0,0)1(x)x0x(0,0)0,0)12dx2(x0,0)1x2xdx0,1(x)x(x10013(x1,1)]1(x)(1,1)2(x)(x1)1(x)10(x)[x0(x)(1,1)(0,0)112xxdx2(x1,1)x2xxxdx0,(1,1)x115。5歡迎下載精品文檔10,132(x)x2355n0，結(jié)點：x012dx212f(x)dx20，A0xxf(0)1313代數(shù)精度為1E(f)rf''()令f(x)x2有r2!1202x2x2dx135因此r1，由此得誤差公式：E(f)1f"()55n1，結(jié)點：x03/5,x13/53/5dx求積系數(shù)A0x2l0(x)dxx2x1111123/53A1x2l1(x)dxx2x3/5dx1111123/5311[f(3/5)f(3/5)]x2f(x)dx13代數(shù)精度為3E(f)rf(4)()令f(x)x4有r4!1x2x4dx1[(3)2(3)2]81355725811f(4)()因此r4!，由此得誤差公式：E(f)725525525方法二： 待定系數(shù)法1x2f(x)dxn01取f(x)1f(x)x

A0f(x0)2A030 A0x0解得：A02,x001x2f(x)dx2f(0)313將f(x)x2代入，等式不再成立。因此代數(shù)精度為1.令f(x)x2有r2!x2x2dx20211351，由此得誤差公式：1因此rE(f)f"()551n1x2f(x)dxA0f(x0)A1f(x1)1。6歡迎下載精品文檔f(x)123A0A1f(x)x0A0x0A1x1f(x)x225A0x02A1x12f(x)x3A0x03A1x130解得：A0A11x03,x13；公式：3,5511f(3)f(3)，將f(x)x4代入，等式不再成立。x2f(x)dx5513因此代數(shù)精度為3.誤差：由于12f(x)dx135)f(35)E(f)rf(4)()1xf(3令f(x)x4，有r4!22(9)268，因此r1，由此得誤差公式：17325725175525E(f)f(4)()52513、（14分）解常微分方程初值問題 y f(t,y), a x b, y(a) y0的一個算法（Milne公式）有如下形式：yi 1 yi3h(0fi 1fi 1 2fi 2)1）確定系數(shù) , 0,1, 2，使算法具有盡可能高的精度，并給出局部截斷誤差；2）請將所得公式與以下公式結(jié)合，組成“預估 -修正-校正-改進”公式：yi1yih[fi25fi119fi9f(ti1,yi1)],E(ti,h)19h5y(5)(i)24720解：yi1y(ti1),Taylor展開y(ti1)y(tih)yihyi1h2yi1h3yi1h4yi(4)1h5yi(5)O(h6)2!3!4!5!yi3y(ti3h)yi(3h)yi1(3h)2yi1(3h)3yi2!3!1(3h)4yi(4)1(3h)5yi(5)O(h6)4!5!h0fi1h0yih1fi1h1y(tih)h1yih21yi1h31yi1h41yi(4)1h51yi(5)O(h6)23!4!。7歡迎下載精品文檔h2fi2h2y(ti2h)h2yi2h22yi2h32yi8h42yi(4)16h52yi(5)O(h6)3!4!11 39292

01201221112222

183yi1yi34h[2fi2fi12fi]43383E(ti,h) y(ti1) yi1h4yi(4)4!1h5yi(5)5!1h5yi(5)5!14h5yi(5)45

1[1(3h)4yi(4)4![1(3h)5yi(5)5!h5yi(5)5!O(h6)

1h4(4)yi(4)88h42yi(4)]3!333!1h5(4)yi(5)816h52yi(5)]O(h6)4!334!O(h6)112h5yi(5)O(h6)35!預估：pi1yi34h[2fi2fi12fi]3校正：yi1yih[fi25fi119fi9f(ti1,pi1)]24對應局部截斷誤差公式：y(ti1)pi114h5yi(5)O(h6),y(ti1)yi119h5y(5)()45720兩式相減yi1pi1243h5y(5)()h5y(5)()720(yi1pi1)72024314720pi1)224pi1)y(ti1)pi1(yi1(yi145243243y(ti1)yi119720(yi1pi1)19(yi1