定积分的近似计算99260.pptVIP

实验四 定积分的近似计算 我们已经学习了定积分的基本概念和定积分 的计算方法，那里所谓的计算方法，是基于 原函数的牛顿 - 莱布尼兹公式。但在许多实 际问题中遇到的定积分，被积函数往往不用 算式给出，而通过图形或表格给出；或虽然 可用一个算式给出，但是要计算它的原函数 却很困难，甚至于原函数可能是非初等函数。 本实验的目的，就是为了解决这些问题，介 绍定积分的“数值积分”，即定积分的近似 计算。 所谓定积分的近似计算，就是找到一个适当的 计算公式，利用被积函数在积分区间上若干个 点处的函数值，来计算定积分的近似值，并作 出误差估计。我们知道，定积分 在几 何上表示 曲线，直线 及 x 轴所围成的曲边梯形的面积。定积分近似计 算的思想，就是将积分区间分割成许多小区间， 然后在小区间上近似计算小曲边梯形的面积， 最后将小曲边梯形的面积求和，就得到了定积 分的近似值。 ? b a dx x f ) ( ) ( x f y ? b x a x ? ? , 1 、 观察黎曼和式的收敛性 由定积分的定义知道，定积分就是黎曼和式 的极限，因此可以用黎曼和式来近似 计算定积分。为计算方便，这里特殊的，将积分 区间等分为 段，并以小区间中点处的函数值作 近似，于是黎曼和式为： i n i i x f ? ? ? 1 ) ( ? n ? ? ? ? ? ? ? n k n a b k a f n a b 1 ) ) 5 . 0 ) 1 (( ( 因而 。 ? ? ? ? ? ? ? ? ? b a n k n a b k a f n a b dx x f 1 ) ) 5 . 0 ) 1 (( ( ) ( 例 1 计算 的黎曼和。 dx x ? 3 2 ln 1 解：输入命令如下： f x_ : 1 Log x ; a 2; b 3; n 200; s NSum f a k 1 0.5 b a n b a n , k, 1, n 述命令是将区间 [2, 3] 等分为 200 段，运行求 得黎曼和为： 1.11842 。 1 、 梯形法 大家可以看出，用上述方法进行的近似计算，其 实是对小曲边梯形的面积用矩形面积来近似，上 面取的特殊的黎曼和又称为中点积分公式。如果 不用矩形而改用梯形来近似，就可以得到定积分 的一个较好的近似方法 —— 梯形积分法。具体方 法如下： 将区间 用 等分为 n 个小区间， 小区间的长度为 。设 ，则每个小梯形的面积为 ，从而得到梯形法的公式为： ] , [ b a b x x x a n ? ? , , , 1 0 ? n a b ? ) ( ) ( n a b i a f x f y i i ? ? ? ? ) , , 1 , 0 ( n i ? ? n a b y y i i ? ? ? ? 2 1 n a b y y y y y y dx x f n n b a ? ? ? ? ? ? ? ? ? ? )] ( 2 1 ) ( 2 1 ) ( 2 1 [ ) ( 1 2 1 1 0 ? ] ) ( 2 1 [ 1 2 1 0 ? ? ? ? ? ? ? ? n n y y y y y n a b ? ] ) ( 2 ) ( ) ( [ 1 ? ? ? ? ? ? ? ? n i n a b i a f b f a f n a b 下面来估计梯形法的误差。第 个小曲边梯形的面 积为 ，做变换 ，则 当在 区间 上连续时，利用分部积分法可 以证明： i ? ? ? ? i i x x i dx x f A 1 ) ( t n a b x x i ? ? ? ? 1 dt n a b t n a b x f dx x f A i x x i i i ? ? ? ? ? ? ? ? ? ? ? 1 0 1 ) ( ) ( 1 ) ( x f ] ,