欢迎来到数值积分(Numerical Integration)的世界!
在纯数学(Pure Mathematics)的学习中,你已经掌握了如何精确地积分许多函数。但如果遇到无法用标准规则积分的函数,例如 \( e^{-x^2} \),该怎么办呢?
这时候,数值积分就能派上用场了!与其寻找精确的代数答案,我们利用巧妙的图形来估算曲线下的面积。想象一下,就像计算圆形地毯的面积时,用许多小矩形砖块把它铺满——你使用的砖块越多,测量结果就越准确。在本章中,我们将探讨三种主要方法:中点规则(Midpoint Rule)、梯形规则(Trapezium Rule)和辛普森规则(Simpson’s Rule)。
1. 核心概念:条带(Strips)与宽度
为了估算曲线在 \(x = a\) 到 \(x = b\) 之间的面积,我们将该区域划分为 \(n\) 个条带。
每个条带的宽度称为 \(h\)。
条带宽度的公式为:
\(h = \frac{b - a}{n}\)
别担心这看起来有很多字母!只要记住:(终点减去起点)除以段数。
快速重温:关键术语
纵坐标(Ordinates):指曲线上各点的 \(y\) 值。
区间(Interval):指从 \(a\) 到 \(b\) 的总距离。
凹性(Concavity):指曲线是向上弯曲(像微笑一样)还是向下弯曲(像皱眉一样)。
2. 中点规则(\(M_n\))
中点规则使用矩形来估算面积。不过,我们不是使用条带起点或终点的高度,而是使用每个条带正中间的高度。
公式:
\(M_n = h(y_{\frac{1}{2}} + y_{\frac{3}{2}} + ... + y_{n-\frac{1}{2}})\)
注:此处的 \(y\) 值是每个条带中点处的高度。
类比:想象你在建造一个楼梯来适配曲线屋顶。如果你测量每个阶梯中间的高度,阶梯“突出”曲线的部分大约会抵消掉曲线下方的“空隙”!
3. 梯形规则(\(T_n\))
我们不再使用平顶的矩形,而是使用梯形(顶部倾斜的形状)。我们用直线将曲线上的点连接起来。
公式:
\(T_n = \frac{1}{2}h(y_0 + y_n + 2(y_1 + y_2 + ... + y_{n-1}))\)
记忆小撇步:将第一个和最后一个高度相加,再加上中间所有高度的两倍。最后,将整个结果乘以条带宽度的一半。
常见错误:学生常会忘记,对于 \(n\) 个条带,会有 \(n+1\) 个高度(\(y\) 值)。计算 \(y\) 值时一定要小心点算!
凹性与误差
梯形规则并不完美,其准确性取决于曲线的形状:
• 如果曲线是向上凹(Concave upwards)(呈 U 型),梯形规则会高估面积,因为直线位于曲线之上。
• 如果曲线是向下凹(Concave downwards)(呈倒 U 型),梯形规则会低估面积。
中点规则的误差方向通常与梯形规则相反!
4. 辛普森规则(\(S_{2n}\))
这是数值积分的“专业”版本。它不使用直线,而是使用抛物线(Parabolas)来拟合条带的顶部,这使得它准确得多!
公式:
\(S_{2n} = \frac{1}{3}h(y_0 + y_{2n} + 4(y_1 + y_3 + ... + y_{2n-1}) + 2(y_2 + y_4 + ... + y_{2n-2}))\)
等等,为什么是 \(2n\)? 在辛普森规则中,我们将条带成对分组,因此总需要偶数个条带。
记忆辅助:“首项 + 末项 + 4 倍奇数项 + 2 倍偶数项”。最后乘以 \(h/3\)。
你知道吗?尽管它是以 Thomas Simpson 命名,但早在几百年前,数学家就已经在使用这种方法了,包括开普勒(Kepler)在计算酒桶体积时也用过它!
5. 各种方法之间的关系
MEI 课程中最酷的部分之一就是了解这些方法是如何联系在一起的。你实际上可以利用简单的估算值来计算出更复杂的估算值!
“神奇”的关联:
1. \(T_{2n} = \frac{1}{2}(M_n + T_n\))
这意味着如果你知道 \(n\) 个条带的中点和梯形估算值,那么条带数加倍后的梯形估算值就是它们的平均值!
2. \(S_{2n} = \frac{1}{3}(2M_n + T_n\))
3. \(S_{2n} = \frac{1}{3}(4T_{2n} - T_n\))
快速重温:准确度阶数
• 中点规则与梯形规则是二阶(second-order)方法。如果你将条带宽度(\(h\))减半,误差大约会除以 4。
• 辛普森规则是四阶(fourth-order)方法。它强大得多!如果你将 \(h\) 减半,误差大约会除以 16。
6. 电子表格(Spreadsheets)的使用
在考试中,你可能会看到电子表格的截图。你不需要成为电脑专家,但必须知道公式看起来是怎样的。
• 单元格(Cell)如 B4 指的是 B 列、第 4 行的值。
• 要找出下一个 \(x\) 值,公式通常是 =A4 + $D$2(其中 D2 是条带宽度 \(h\))。
• $ 符号非常重要——它能“锁定”该单元格,这样当你向下拖曳公式时,它不会改变参考位置。
总结:关键要点
• 数值积分用于我们无法进行代数积分时寻找面积。
• h 是一个条带的宽度:\(h = \frac{b-a}{n}\)。
• 凹性能告诉你估算值是偏高还是偏低。
• 辛普森规则最准确,使用二次(抛物线)曲线拟合。
• 计算器:如果函数涉及三角函数,请务必确认你的计算器处于弧度(Radians)模式!
• 步骤:1. 找出 \(h\)。2. 列出 \(x\) 值。3. 计算 \(y\) 值。4. 将它们代入公式手册提供的公式中。
别担心,公式看起来吓人吗?考试时你随时都可以查阅公式手册。你的主要任务是识别正确的 \(y\) 值,并细心地代入公式中。