玩转 matlab 之一维 gauss 数值积分公式及matlab源代码

最新推荐文章于 2024-10-03 09:13:22 发布

孙开心2020

最新推荐文章于 2024-10-03 09:13:22 发布

阅读量1.6w

点赞数 21

分类专栏： matlab gauss积分文章标签： gauss积分数值积分 matlab 有限元数值分析

本文链接：https://blog.csdn.net/sunzhenwei1990/article/details/90084751

版权

本文介绍了如何在MATLAB中应用一维Gauss数值积分公式进行数值积分，并通过两个函数的实验对比了Gauss积分与MATLAB内置integral函数的结果，展示了随着Gauss点数增加，积分精度提升的过程。实验表明，对于多数问题，2点或3点的Gauss公式已具备高精度。

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

在数值分析中，尤其是 有限元刚度矩阵、质量矩阵等的计算中，必然要求如下定积分:

I=\int_a^b f(x)dx

学好 gauss积分也是学好 有限元的重要基础，学过高等数学的都知道，手动积分能把人搞死（微笑脸），而且有些函数还不存在原函数，使用原始的手动算出原函数几乎是不现实的。因此非常有必要学习数值积分，简单讲就是近似计算，只要这个近似值精确度高和稳定性好就行。Gauss积分公式就是这么一个非常好用的工具。本文介绍高斯积分公式的使用以及简单的数值算例。

标准区间

先考虑特殊情况，对于一般区间呢？待会会处理这个问题。
$I=\int_{-1}^1 f(x)dx$
不加证明的直接给出gauss公式如下：详情参阅任何一本数值分析书都有详细的证明过程：
$I=\int_{-1}^1 f(x)dx=\Sigma_{i=1}^n A_if(x_i)$
其中 $A_i$ 称作权， $x_i$ 称作 gauss 点。
下面的问题就是如何选择 $n,A_i,x_i$ 。