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

释放双眼，带上耳机，听听看~！

目录

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

\\[ 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\\)。

理论表明n个点的Gauss公式代数精度为\\(2n-1\\)，其选择如下表，(这里仅仅举1-4个点情况，实际使用的时候一般2点或者3点的精度已经完全够了)更多积分点可参考 gauss表.

gauss点个数 \\(n\\)

gauss 点 \\(x_i\\)

权重 \\(A_i\\)

精度

1

\\(x_1\\)=0

\\(A_1