估计积分 I=∫ₐ^(b)(sin x)/(x)dx 的值,其中 0lt alt b.
解
如果利用 |sin x| ≤slant |x|,则只能得到
| I | ≤slant b - a. (11.17)这个估计在 b-a 较大时当然很差.
利用积分第一中值定理,由于因子 1 / x 不变号,就得到
| I | = | sin ξ ∫_ a ^ b frac d xx | ≤slant ∫_ a ^ b frac d xx = ln frac ba = ln (1 + frac b - aa) < frac ba - 1. (11.18)如果 a > 1 ,则这个估计比 (11.17) 要好.但是对于 b - a 很大的情况仍然不好
利用积分第二中值定理和因子 1 / x 单调非负,就有
| I | = | frac 1a ∫_ a ^ ξ sin x d x | ≤slant frac 1a | cos ξ - cos a | ≤slant frac 2a. (11.19)如果 a > 1 ,则对于 b - a 很大的情况这个估计明显比前两个要好.观察该被积函数的特性(其图像见图 4.2(a)),当积分区间较大时,必须将 sin x 所起的正负抵消的作用考虑进去,而不能如前两个估计那样只利用 |sin x| ≤slant |x| 和 |sin x| ≤slant 1.这就是估计 (11.19) 优于它们的理由所在.
估计 (11.19) 与区间右端无关.下面一个估计则对任何 0 ≤slant a < b 都成立:
| ∫_ a ^ b frac sin xx d x | < 3. (11.20)这里不妨设 0 ≤slant a < 1 < b,否则下面的估计更为简单.这时将积分拆开,并利用 (11.17) 和 (11.19) 就得到
| I | ≤slant | ∫_ a ^ 1 frac sin xx d x | + | ∫_ 1 ^ b frac sin xx d x | < 1 + 2 = 3. (11.21)当然这还是一个很粗的估计.与区间无关的最优估计问题留作为本章第二组参考题 5.又若 a, b 是给定的具体数值,则还需要利用被积函数在 [a, b] 上的具体特性才能作出较好的估计.
下面就是一个例子 (即是 [27] 中的 2328 题).
估计定积分 I = ∫₁00π²00π(sin x)/(x)dx 的值.
解 1
将积分区间按 π 的整倍数拆开,就不难证明 I > 0 (细节从略). 又利用上面已有的估计式 (11.19),就得到 [27] 中的答案:
I = frac θ5 0 π, 0 < θ ≤slant 1.解 2
利用被积函数的分母大于 100π ,而如下分部积分后会变得更大,就有
r I = (- frac cos xx - frac sin xx ^ 2) g | _ 1 0 0 π ^ 2 0 0 π - 2 ∫_ 1 0 0 π ^ 2 0 0 π frac sin xx ^ 3 d x = frac 12 0 0 π - 2 ∫_ 1 0 0 π ^ 2 0 0 π frac sin xx ^ 3 d x.对于右边的积分用积分第二中值定理估计,就有
| I - frac 12 0 0 π | ≤slant frac 21 0 ^ 6 π^ 3 | ∫_ 1 0 0 π ^ ξ sin x d x | ≤slant frac 41 0 ^ 6 π^ 3 ≈ 0. 1 2 9 × 1 0 ^ - 6.也就是说积分 I ≈ (1)/(200π) ≈ 0.00159.
用 Mathematica 计算得到 I ≈ 0.00159149.实际上,如果对上面的最后一个积分再用分部积分和第二中值定理,则可以得到更为精确的近似值.
各种不等式在估计中都可能有用.其中 Cauchy 不等式与 Schwarz 不等式更是常用的工具.例如,椭圆 (x²)/(a²) + (y²)/(b²) = 1 (a, b > 0) 的周长为
s = 4 ∫_ 0 ^ π / 2 sqrt a ^ 2 sin^ 2 t + b ^ 2 cos^ 2 t d t,其中被积函数的原函数不是初等函数,因此不可能用 Newton-Leibniz 公式来计算 (可参看 [50] 中的 138–142 页). 下面我们分别用 Schwarz 不等式和 Cauchy 不等式来估计这个积分的上界和下界.它们的平均值即是在某些数学手册中关于椭圆周长的近似公式之一.
证明: π (a + b)≤slant s≤slant π √(2a² + 2b²).
证
首先,不难由 Schwarz 不等式得到上界估计:
l s = 4 ∫_ 0 ^ π / 2 sqrt a ^ 2 sin^ 2 t + b ^ 2 cos^ 2 t d t ≤slant 4 [ ∫_ 0 ^ π / 2 (a ^ 2 sin^ 2 t + b ^ 2 cos^ 2 t) d t ] ^ frac 12 (∫_ 0 ^ π / 2 d t) ^ frac 12 = 4 [ frac π4 (a ^ 2 + b ^ 2) ] ^ frac 12 (frac π2) ^ frac 12 = π sqrt 2 a ^ 2 + 2 b ^ 2.然后 (反方向) 用 Cauchy 不等式求出被积函数的下界:
r l & sqrt a ^ 2 sin^ 2 t + b ^ 2 cos^ 2 t = sqrt a ^ 2 sin^ 2 t + b ^ 2 cos^ 2 t · sqrt sin^ 2 t + cos^ 2 t & ≥slant a sin t · sin t + b cos t · cos t & = a sin^ 2 t + b cos^ 2 t,对两边积分再乘 4 就得到积分的下界估计:
s ≥slant 4 ∫_ 0 ^ π / 2 (a sin^ 2 t + b cos^ 2 t) d t = π (a + b).在数学分析教科书中对于积分的近似计算一般是介绍三种方法,即梯形公式、矩形公式和抛物线公式 (也称为 Simpson (辛普森) 公式). 这些公式以及更深入的数值积分方法都是以下面的定理为基础的.
设函数 f ∈ C^((2m+2))[a, b] , h = (b - a)/n , xᵢ = a + ih , i = 0, 1, …, n ,则
l frac b - an ∑_ i = 1 ^ n frac 12 [ f (x _ i - 1) + f (x _ i) ] - ∫_ a ^ b f (x) d x = ∑_ k = 1 ^ m frac B _ 2 k(2 k) ! h ^ 2 k [ f ^ (2 k - 1) (b) - f ^ (2 k - 1) (a) ] + frac B _ 2 m + 2(2 m + 2) ! h ^ 2 m + 2 f ^ (2 m + 2) (ξ) (b - a), (11.22)其中 ξ∈[a,b] , B₂k(k=1,2,…,m+1) 是 Bernoulli 数 (见 7.2.3 小节), 其中前三个是:
B _ 2 = frac 16, B _ 4 = - frac 13 0, B _ 6 = frac 14 2.容易看出,上述公式的左边就是对于区间 [a, b] 的 n 等距分划下的梯形公式与积分之差,因此公式给出了梯形公式的误差表达式.通过组合就可以得到矩形公式和抛物线公式的误差估计.
公式 (11.22) 的证明并不困难,主要是用分部积分法.下面给出的几个例题主要是介绍方法.将例题中所得的结果用于等距分划的每个子区间,并加以合并就可以得到 m = 0, 1 时的 Euler-Maclaurin 公式.它们提供了梯形公式和矩形公式的误差估计.
设 f ∈ C²[0, h],则存在 ξ ∈ [0, h],使成立
∫_ 0 ^ h f (x) d x = frac h2 [ f (0) + f (h) ] - frac 11 2 f ^ prime prime (ξ) h ^ 3. (11.23)证
如下用两次分部积分得到
r l ∫_ 0 ^ h f (x) d x & = ∫_ 0 ^ h f (x) (x - frac h2) ^ prime d x & = f (x) (x - frac h2) | _ 0 ^ h - ∫_ 0 ^ h f ^ prime (x) (x - frac h2) d x & = frac h2 [ f (0) + f (h) ] - frac 12 ∫_ 0 ^ h f ^ prime (x) [ x (x - h) ] ^ prime d x & = frac h2 [ f (0) + f (h) ] + frac 12 ∫_ 0 ^ h f ^ prime prime (x) [ x (x - h) ] d x. (11.24)由于 x(x - h) 不变号,对右边的积分用第一中值定理就得到所要的结果:
frac 12 ∫_ 0 ^ h f ^ prime prime (x) [ x (x - h) ] d x = frac 12 f ^ prime prime (ξ) ∫_ 0 ^ h x (x - h) d x = - frac 11 2 f ^ prime prime (ξ) h ^ 3.如果引进变上限积分 F(x) = ∫₀ˣ f(t)dt, 0 ≤slant x ≤slant h ,就可以看出本题与第七章第一组参考题 9(3) 相同。两者在条件和结论上的差异不是本质的,只要应用第十章第一组参考题 11 就可以解决。此外,这也说明在积分学中的许多问题用微分学也是可以解决的。
设 f ∈ C⁴[0, h] ,证明:存在 ξ ∈ [0, h] ,使成立
∫_ 0 ^ h f (x) d x = frac h2 [ f (0) + f (h) ] - frac h ^ 21 2 [ f ^ prime (h) - f ^ prime (0) ] + frac 17 2 0 f ^ (4) (ξ) h ^ 5. (11.25)证
从上一例题中推导得到的等式 (11.24) 继续做下去:
r l & ∫_ 0 ^ h f (x) d x - frac h2 [ f (0) + f (h) ] & = frac 12 ∫_ 0 ^ h f ^ prime prime (x) [ x (x - h) ] d x & = - frac h ^ 21 2 ∫_ 0 ^ h f ^ prime prime (x) d x + frac 12 ∫_ 0 ^ h f ^ prime prime (x) (x ^ 2 - h x + frac h ^ 26) d x & = - frac h ^ 21 2 [ f ^ prime (h) - f ^ prime (0) ] + frac 12 ∫_ 0 ^ h f ^ prime prime (x) (frac x ^ 33 - frac h x ^ 22 + frac h ^ 2 x6) ^ prime d x & = - frac h ^ 21 2 [ f ^ prime (h) - f ^ prime (0) ] - frac 12 ∫_ 0 ^ h f ^ prime prime prime (x) (frac x ^ 33 - frac h x ^ 22 + frac h ^ 2 x6) d x & = - frac h ^ 21 2 [ f ^ prime (h) - f ^ prime (0) ] - frac 12 ∫_ 0 ^ h f ^ prime prime prime (x) (frac x ^ 41 2 - frac h x ^ 36 + frac h ^ 2 x ^ 21 2) ^ prime d x & = - frac h ^ 21 2 [ f ^ prime (h) - f ^ prime (0) ] + frac 12 4 ∫_ 0 ^ h f ^ (4) (x) [ x ^ 2 (x - h) ^ 2 ] d x,在最后一个积分中,利用因子 x²(h - x)² 不变号,再用积分第一中值定理就可得到所要的等式.square
以上结果对于梯形公式和矩形公式的误差估计已经够用.为了对抛物线公式作出误差估计,还需要 m = 2 时的 Euler-Maclaurin 公式,其证明方法与上面完全一样,读者可自己完成.此外,在 [14, 17] 等教科书中均对抛物线公式采用微分学方法作出误差估计.下面只是对于抛物线方法中的基本公式作一点介绍.
若 p(x) 是不超过 3 次的多项式,则有
∫_ a ^ b p (x) d x = frac 16 [ p (a) + 4 p (frac 12 (a + b)) + p (b) ] (b - a). (11.26)证
令 q(t) = f(a + t(b - a)) ,就可以将要证明的公式变为等价的
∫_ 0 ^ 1 q (t) d t = frac 16 [ q (0) + 4 q (frac 12) + q (1) ].然后利用积分为线性运算,分别用 q(t) = 1, t, t², t³ 代入验算即可.
这个公式在初等数学的体积计算中有万能公式的美名:只要将一个立体的顶截面、中截面和底截面的面积分别乘 1:4:1 并相加,然后除以 6,再乘高度即可.容易验证它对于球、圆锥和圆台等形体的体积都给出了准确的答案.
11.3.3 练习题
Section titled “11.3.3 练习题”- 证明:
(1) ∫₀^(√(2π))sin x²dx>0; (2) (1)/(203√(2))<∫₀¹fracx¹9sqrt[3]1+x⁶dx<(1)/(20); (3) ∫₀^(π/2)x((sin nx)/(sin x))⁴dx<fracn²π²4; (4) 0<(π)/(2)-∫₀^(π/2)(sin x)/(x)dx<fracπ³144; (5) 0.005<∫₀¹00frace⁻xx+100dx<0.01; (6) 2/9π²<∫_(π/6)^(π/2)(2x)/(sin x)dx<1/3π².
- 设 f 在 [0, a] (a > 0) 上有可积的导函数,证明:
- 设 f 在 [0,1] 上有可积的导函数,证明:
- 设函数 f 在 [a, b] 上可微,|f'(x)| ≤slant M,且 ∫ₐ^(b) f(x) dx = 0.对于函数 F(x) = ∫ₐˣ f(t) dt ,
(1) 证明: |F(x)| ≤slant (M(b - a)²)/(8) ;
(2) 在增加条件 f(a) = f(b) = 0 时证明: |F(x)| ≤slant (M(b - a)²)/(16) .
- 证明: 对每个正整数 n,成立
- 设 f 在 [a, b] 上可微,f(a) = f(b) = 0 , |f'(x)| ≤slant M,证明:
- 设 f 在 [a, b] 上二阶可微,f(a) = f(b) = 0 , |f''(x)| ≤slant M,证明:
- (矩形公式) 设 f ∈ C²[a, b],证明: 存在 ξ ∈ (a, b),使成立
- 设 f 于 [-1, 1] 上可微,且有 a ∈ (0, 1),使得 ∫_(-a)^(a) f(x) dx = 0,证明:
- 设 f 于 [0,1] 上可微,|f'(x)| ≤slant M,证明:
- 设 f 在区间 [0,1] 上可微,且 f' ∈ R[0,1],对正整数 n 定义
证明: limₙ→ ∞nAₙ = 1/2 (f(0) - f(1)).
- 设 f 在区间 [0,1] 上二阶可微,且 f'' ∈ R[0,1],对正整数 n 定义
证明: limₙ→ ∞n² Bₙ = (1)/(24) [f'(1) - f'(0)].