数值积分与数值微分习题课.docxVIP

数值积分与 数值微分 习题课 一、已知x0 1,x1 1,x2 3，给出以这3个点为求积节 4 2 4 点在0.1上的插值型求积公式 解：过这3个点的插值多项式基函数为 l02 x x1 x x1 x0 l12 x x0 x x0 x1 l22 x x0 x x0 x2 1  x2x0x2 x2x1x2 x1x2x1 Ak0lk2xdx,k0,1,2 A 1 x x1 x x2dx 0 0 x x x x 2 0 1 0  1 x 2 011 42  3 x 4 3 4  2 dx 3 x  1  3 x A1 1 x x x x x1 0 x1 2dx 0 x0 x2 A 1 x x0 x x1dx 2 0 x x x x 0 2 2 1  4 011 24 1 x 4 031 44  4 3 4 1 2 1 2  1 dx 3 2 dx 3 故所求的插值型求积公式为 fxdx 2f 1 1f 1 2f 3 1 0 3 4 3 2 3 4 二、确定求积公式 1 15f0.68f05f 0.6 fxdx 1 9 的代数精度，它是Gauss公式吗？ 证明：求积公式中系数与节点全部给定，直接查验 依次取f x1,x,x2,x3,x4,x5，有 2 1 1 51 81 5 1 1dx 1 9 0 1 1 5 0.6 8 05 0.6 xdx 1 9 2 3 0 2 5 0  1 x2dx 1 1 x3dx 1 1 x4dx 1 1 x5dx 1  5 5 5 1 5 9  0.6 2 02 5 2 8 0.6 0.6 3 03 5 3 8 0.6 0.6 4 04 5 4 8 0.6 0.6 5 05 5 5 8 0.6 此题已经达到2n-1=5。故它是Gauss公式。 三、试应用复合梯形公式计算积分 21 Idx 12x 要求误差不超过103，并把计算结果与正确值比较。 解：复合梯形公式的余项为 Rf,T b f(x)dxT bah2f() n a n 12 ba n1 Tn f(a)f(b)2 f(xk) 2n k1 xakh,h ba,k0,1,2,L,n k n 此题fx 此题余项为 Rf,Tn 要使Rf,Tn =b-a 2-1 n 0.1 h  1 ，fx 13,M2maxfx1 2x x x1,2 21h2f() h2 maxf(x) h2 12 12x1,2 12 h2 ，得h 0.109545，取h0.1 103 12 于是有 IT10 1 1 1 1 1 ... 1 2 10 2 2 2 1.1 2 1.2 0.346886 4 2 1.9 查验：I T10 1ln2 T10 3103 2 四、证明若函数fxC1a,b，则其上的一阶差商函数是连续函数，并借助此结果用Newtong插值余项证明梯形求积公式的余项为 Rf fxdxbafafb ba 3 f b 1 a 2 12 证明：不妨设一阶差商函数为 fx,a，x0 a,b，有 lim fx0 h,a f x0 h f a lim x0 h a h0 h 0 f x0 f h fa lim x0 h a h 0 fx0 fa f h fx0 fa lim x0 ha x0 ha x0 fx0,a h 0 a 由x0的随意性，可知一阶差商函数是连续函数。 由插值特点，显然有 R1 b fxL1xdx b f fxN1xdx a a 线性插值的Newton余项公式为 fxN1xfx,a,bxaxb 故有 R1f b fx,a,bxaxbdx a 由 lim f x h,a,b fx h,a f a,b lim xh b h0 h 0 lim f xh,a f a,b f x,a f a,b h 0 xb x b fx,a,b 可知fx,a,b是变量x在a,b上的连续函数，而函数 axb在a,b上可积，不变号，根据积分中值 定理，存在a,b，使 b afxN1xdxf,a,b 由差商性质，存在 a,b，使f b xN1x f f dx a 2 b 3 a f 12  b a xbdx x a ,a,b f 。所以 2 b x bdx xa a 结论得证。 五、导出中矩形公式 b xdx b a b f af 2 a 的余项。 解：将f x 在x a 2 b处进行泰勒展开 fx fab fab xab 1f xab2 2 2 2 2 2 a,b。对上式两边在a,b 上积分，有 b b ab b ab ab 1b ab 2 f dx f x x dx fxdx 2 2 dx f 2 a a a 2 2a ab b ab ab 1 b x ab 2 baf 2 f 2 x dx f 2 dx a 2 2a 中矩形公式的余项 R