数值计算笔记之数值积分(一)

最新推荐文章于 2024-07-22 19:30:03 发布

MinJinFan

最新推荐文章于 2024-07-22 19:30:03 发布

阅读量1.1w

点赞数 7

分类专栏：数值计算 C++

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/qq_37264323/article/details/104370486

版权

本文详细介绍了数值积分的概念、方法及其应用。内容包括积分思想，如中矩形、梯形和辛普森公式，讨论了求积公式的余项和代数精度，介绍了牛顿-柯特斯公式和复化求积法，特别是复化梯形和辛普森公式。此外，还提供了C++和MATLAB的代码实现示例。

摘要由CSDN通过智能技术生成

目录

一、数值积分的积分思想

1、中矩形公式

2、梯形公式

3、辛普森公式

二、求积公式的余项和代数精度

三、插值型求积公式

四、牛顿--柯特斯公式 (N-C公式)

五、复化求积法

1、复化梯形公式

2、复化辛普森公式（要求 n 为偶数）

六、代码实现

0、引言

在高数中，可以根据 $\int_{a}^{b}f(x)dx = F(b)-F(a)$ ，求得积分。但是如果 f(x) 存在以下两种情况：

是离散函数
无法求得原函数

那么对于这种情况，就可以利用数值积分。

所谓数值积分，指利用被积函数 f(x) 在有限个点上的函数值来计算积分近似值的一种方法。

一、数值积分的积分思想

积分中值定理： $\int_{a}^{b}f(x)dx = (b-a)f(\xi ),\xi \in [a,b]$

一重积分可以理解为求面积，积分中值定理可以看成：在 $x\in [a,b]$ 上，函数 f(x) 与轴的面积等于以 (b-a) 为长， $\xi \in [a,b]$ 中一点函数值 $f(\xi )$ 为宽的矩形的面积。

1、中矩形公式

$\int_{a}^{b}f(x)dx \approx (b-a)\cdot f(\frac{a+b}{2})$

2、梯形公式

$\int_{a}^{b}f(x)dx\approx (b-a)\cdot \frac{f(a)+f(b)}{2}$

3、辛普森公式

$\int_{a}^{b}f(x)dx\approx (b-a)\cdot \frac{1}{6}[f(a)+4f(\frac{a+b}{2})+f(b)]$

一般情况下，可用： $f(\xi )\approx \sum_{k=0}^{n}w_{k}\cdot f(x_{k})$

$\Rightarrow \int_{a}^{b}f(x)dx = (b-a)\cdot f(\xi )\approx (b-a)\cdot \sum_{k=0}^{n}w_{k}\cdot f(x_{k}) = \sum_{k=0}^{n}A_{k}f(x_{k})$ ，其中 $A_{k}=(b-1)w_{k}$ 。

上式称为数值积分的基本公式， $x_{k}$ 称为求积节点， $A_{k}$ 称为求积系数。

二、求积公式的余项和代数精度

1、余项

$R\begin{bmatrix} f \end{bmatrix}=\int_{a}^{b}f(x)dx-\sum_{k=0}^{n}A_{k}\cdot f(x_{k})$

2、代数精度

若数值积分公式对任意的 $\leq m$ 次的代数多项式都准确成立，而对于 $x^{m+1}$ 却不能成立，则称该数值积分公式的代数精度为。

举个例子理解一下：求梯形公式的代数精度。

解： $\int_{a}^{b}f(x)dx\approx (b-a)\cdot \frac{f(a)+f(b)}{2}$

当，
当， $left = \frac{1}{2}(b^{2}-a^{2})= right$
当

最低0.47元/天解锁文章

关注

7
点赞
踩
65

收藏

觉得还不错? 一键收藏
0
评论
数值计算笔记之数值积分(一)

0、引言在高数中，可以根据，求得积分。但是如果存在以下两种情况：是离散函数无法求得原函数那么对于这种情况，就可以利用数值积分。所谓数值积分，指利用被积函数在有限个点上的函数值来计算积分近似值的一种方法。一、数值积分的积分思想积分中值定理：一重积分可以理解为求面积，积分中值定理可以看成：在上，函数与轴的面积等于以为长，中一点函数值为宽的矩形的面积。1、中矩形...
复制链接

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。