清華大學(xué)楊頂輝數(shù)值分析第6次作業(yè)

1、19.令Tn*(x) Tn(2x 1),x 0,1，試證T；(x)是在0,1上帶權(quán)(x)的正交多項(xiàng)式，并求 T0(x),T1*(x),T2*(x),T；(x).證明:(x)Tn*(x)T；(x)dxT=rTn(2x 1)Tm(2x 1)dx0.X x令t 2x 1,那么(x)Tn*(x)T；(x)dx1 1 11 Tn(t)Tm(t)dt2 1 t 1 (t 1)2:2 ( 2 )1由切比雪夫多項(xiàng)式Tn(x)關(guān)于權(quán)函數(shù)。1 2的正交性，可得V1 x20 m n1 * * 1 1 1(x)Tn(x)Tm(x)dx2Tn(t)Tm(t)dt=； m n 001 1 t22m n 0所以T；(x)是

2、在0,1上帶權(quán)(x)的正交多項(xiàng)式Jx X2T0(x) T°(2x 1) 1T1* (x) T1(2x 1) 2x 1T2 (x)E(2x 1)2(2 x1)21 8x2 8xT；(x)T3(2x 1)4(2x1)33(2x 1) 32x348x218x 114.實(shí)驗(yàn)數(shù)據(jù)如下:Xi1925313844yi用最小二乘法求形如y a bx2的經(jīng)驗(yàn)公式，并求均方誤差解：法方程為(1,1) (1,x2) a(1,y)(x2,1) (x2,x2) b(x2,y)55327 a271.45327 7277699 b369321.5解得a 0.972579b 0.050035擬合公式為 y 0.97

3、2579 0.050035X2均方誤差422yi a bx20.015023i 021給出f (x) In x的函數(shù)表如下:xln x用拉格朗日插值求ln 0.54的近似值并估計(jì)誤差計(jì)算取n 1及n 2解：n 1 時(shí)，取 x00.5, x10.6由拉格朗日插值定理有1x 06x 05L.(x)f (xj)lj(x)0.6931470.510826j 00.5 0.60.6 0.51.82321 0x 1.604752所以 ln0.54 L1(0.54)0.620219誤差為 ln0.54 ( 0.620219)= 0.004032n2 時(shí)，取 x° 0.4,捲 0.5, X2 0.6

4、由拉格朗日插值定理有2L2(x)f(Xj)lj(x)j o(x0.5)(x 0.6)(x0.4)(x 0.6)(x0.4)(x 0.5)0.9162910.6931470.510826 (0.40.5)(0.40.6)(0.50.4)(0.50.6)(0.60.4)(0.60.4)22.041150x4.068475x 2.217097所以 InO.54 L2(0.54)0.615320誤差為InO.54 ( 0.615320)8.662994 1023.建立三次樣條插值函數(shù)s(x)，并求f(O)的近似值s(O)，這里已給函數(shù)表。Xif(Xi)邊界條件 s''( 0.3) s&

5、#39;'(O.3) O1 1解：由剖分節(jié)點(diǎn)可知hO h1 h2 °2 12 2，12 -d1 6 f x0, x1, x2 6.4215, d2 6 f x1, x2, x3 6.4215得到方程組M1M26.42156.4215解得 M1 M 22.5686注意到Mo M3 0得到三次樣條插值函數(shù)s(x)2.14O5X3 1.92645X2 1.O64215X 0.000632,x O.3, O.1s(x) 1.2843X2 x O.OO2773,x O.1,O.12.14O5X3 1.92645X2 0.935785x 0.000632,x O.1,O.329.確定以下

6、求積公式的待定參數(shù)，使其代數(shù)精度盡量高，并指出所構(gòu)造出的求積公式所具有的代數(shù)精度:hBf (0)Cf(h)1 f(x)dx Af ( h)h解:取 f(x) 1,x,x2 ，有2hh1dx A B Chhxdx Ah Chh2h33hx2dxhAh2 Ch2解以上方程,ih，B-h3求積公式為hf (x)dxhf(h)4hf(0)3-hf (h)3取 f(x)3 x，x3dxh0,Af(h)Bf(0) Cf(h)取 f (x)4 xh2,5 1 , zB 01 ,25x4dxh5 -h(h)4hh4-h5h5333h所以x4dxAf(h)Bf(0)Cf(h)hh0h因此構(gòu)造的求積公式代數(shù)精度為

7、30證明求積公式f(x)dx -5f( ,|)8f(0)5f(：)的代數(shù)精度為5證明：取f(x) xi1那么對(duì)任意的i 0，有 f(x)dx11xidx11 i +1,1 lx 1|1_1i1|-1(1)i1i 10, i為奇數(shù),i為偶數(shù)i 1i為奇數(shù)時(shí)，f (x)為奇函數(shù)95f(，5)8f(0)0時(shí)5f (95匕 3)=15f(8f(0)5f(3)3)=訥5f5=22時(shí),扣(3)8f(0)5f(315)=954時(shí)5f(93)8f(0)5f(3)T56時(shí)5f (98f(0)5f(1= 9-5f(1x0dx1=l1x2dx1從而有E(xi)0,0i 5,E(x6)所以求積公式的代數(shù)精度為、 、1

8、 133.求 X1,X2,A1,A,使公式=f (x) 0 Vx解：構(gòu)造區(qū)間0,1上關(guān)于權(quán)函數(shù)由P2(x)與俐X正交，可得0° 丄 P,x)dx0 vxax b)dx92527125251755卒2551x4dx1276125= 251x6dx15f=03)Af(X!)A2f(X2)為高斯型求積公式1-的二次正交多項(xiàng)式P2(x) X2 ax b00 丄xP2(x)dx0 VxXx(x0 -、Xax b)dx-a 2b53聯(lián)立以上方程,可得a§ b 27353_35求P2(X)的零點(diǎn)即為高斯求積公式的節(jié)點(diǎn)所以 F2(x) x2x1 1(3 2.675)X2 丄(3 2、675)