数值积分 近似计算=f(x) 插值型积分公式 81插值型求积公式 思利用插值多项式P(x)≈对积分易算 在b上取a≤xx1<,A≤b,做/的n次插值多 项式L(x)=∑f(xk即得到 =0 误差R门 mk4由节点决定,1/()-((x 与/(0无关,-(15mxx)
数值积分 近似计算 = b a I f (x)dx 8.1 插值型求积公式 思 路 利用插值多项式 Pn (x) f (x) 则积分易算。 在[a, b]上取 a x0 < x1 <…< xn b，做 f 的 n 次插值多 项式 ，即得到 = = n k Ln x f xk l k x 0 ( ) ( ) ( ) = b a b a k n k f (x)dx f (xk ) l (x)dx 0 Ak − − = b a j k x x x x Ak dx k j j ( ) ( ) 由 决定， 与 无关。 节点 f (x) 插值型积分公式 = + = − + = = − = = − b a n k k x n b a n b a n b a n k k k x x dx n f f x L x dx R x dx f x dx A f x R f 0 ( 1) 0 ( ) ( 1)! ( ) [ ( ) ( )] ( ) ( ) ( ) [ ] 误差
2复化求积公式 如果积分区间比较大,直接地使用上述求积公式 精度难以保证 高次插值有 Runge现象,故采用分段低次插值 →分段低次合成的 Newton- Cotes复合求积公式。 (1)等分求积区间,比如取步长h b-a 分|a,b为n等分 分点为x=x+Mk=0,1,2,,n (2)在区间[x,xk1上使用以上求积公式求得 (3)取和值Ⅰ=∑/作为整个区间上的积分近似值
8.2 复化求积公式 如果积分区间比较大，直接地使用上述求积公式， 精度难以保证。 高次插值有Runge 现象，故采用分段低次插值 分段低次合成的 Newton-Cotes 复合求积公式。 （1）等分求积区间，比如取步长 ，分[a, b]为n等分， n b a h − = 分点为 xk = x0 + kh k = 0, 1, 2,…, n （2）在区间 [xk , xk+1]上使用以上求积公式求得Ik （3）取和值 ，作为整个区间上的积分近似值。 − = = 1 0 n k k I I
b-a 复化梯形公式〃n’x=a+kh(k=0,…,n) 在每个xx上用梯形公式:。8888。 ∫ f(x)≈,=∫(xk1)+f(xk),k=1,…,n (xk∑,[(x-)+1(x)=a)+2x)+/(b)=7 k=1 ∑∫(5) R=∑ h h ∫"(5k) (b-a) k=1 积分中值定理* 12 12 h2 (b-a)f"(4),∈(a,b)
➢ 复化梯形公式： , x a k h (k 0, ... ,n) n b a h k = + = − = 在每个 [xk−1 , xk ] 上用梯形公式： f x f x k n x x f x dx k k x x k k k k [ ( ) ( )], 1, ... , 2 ( ) 1 1 1 + = − − − − = + + − = 1 1 ( ) 2 ( ) ( ) 2 n k k f a f x f b h = − + b a n k k k f x f x h f x dx 1 1 [ ( ) ( )] 2 ( ) = Tn ( ) ( ), ( , ) 12 ( ) ( ) 12 ( )] 12 [ ] [ 2 1 2 1 3 b a f a b h n f b a h f h R f n k n k k k = − − = − = − − = = /*积分中值定理*/
复化 Simpson公式:h=b-a n Mk=a+kh (k=0, ∫(x)b≈If(xk)+4f(xk+4)+∫(xk+1) k k k+1 ∫(x)t≈么a)+(xn)+2(x,)+r(=Sn k=0 =-b-a/h) fn(5) 180(2 主为方便编程,可采用另一记法:令n=2n为偶数 这时n=b=a=B,x1=a+kM,有 Sn=,Ufa)+4∑f(x)+2∑f(x)+∫(b) odd k even k
➢ 复化 Simpson 公式： , x a k h (k 0, ... ,n) n b a h k = + = − = [ ( ) 4 ( ) ( )] 6 ( ) 1 2 1 1 + + + + + k k k x x f x f x f x h f x dx k k xk 2 1 k+ x xk+1 4 4 4 4 4 [ ( ) 4 ( ) 2 ( ) ( )] 6 ( ) 1 0 1 0 1 2 1 − = − = + + + + + n k n k k k b a f a f x f x f b h f x dx = Sn ( ) 180 2 [ ] (4) 4 f b a h R f − = − 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 h x a k h ，有 n b a h k = = + − = , 2 [ ( ) 4 ( ) 2 ( ) ( )] 3 + + + = odd k even k n k k f a f x f x f b h S
例8.1:利用数据表 3/8 12 5/8 3/4 7/8 f( 3.938463.764703.506853.200002876402.460002.26549 计算积分 2 1+x 这个问题有明显的答案r=4arcx=x=31415926 取n=8用复化梯形公式 T=×f(0)+2+21|+2f|。+2 8 2 +2 8/12/3 4/+2/ 十 8 =3.13899
例 8．1：利用数据表 xk 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (xk ) 4 3.93846 3.76470 3.50685 3.20000 2.87640 2.46000 2.26549 2 计算积分 dx x I + = 1 0 2 * 1 4 这个问题有明显的答案 I * = 4arctgx | 1 0 = = 3.1415926 取n = 8用复化梯形公式 ( ) 3.13899 1 8 7 2 4 3 2 8 5 2 2 1 2 8 3 2 4 1 2 8 1 (0) 2 2 1 8 1 8 = + + + + + + + = + f f f f T f f f f f