数值积分的基本思想
 
由积分中值定理可知，在积分区间 
 

   内存在一点 
 
 

   ，成立
 
 

 

  
 
 

  
 
 

 

   式的几何意义即为：底为 
 
 

   而高为 
 
 

   的矩形的面积恰等于所求曲边梯形的面积 
 
 

   。因此，要想求出 
 
 

   式左端积分，我们只需要知道三个值： 
 
 

   即可，这里 
 
 

   是显然的，问题在于 
 
 

   的具体位置一般是不清楚的，从而 
 
 

   未知。我们暂且将 
 
 

   称为区间 
 
 

   上的平均高度，我们的目标就是寻求一种求出平均高度 
 
 

   的算法，这样我们的问题也就解决了。
 
 

 
首先我们想到的是用区间两端点的“高度” 
 

   与 
 
 

   的算术平均值作为平均高度 
 
 

   的近似值，从而导出下述求积公式
 
 

 

  
 
 

  
 
 

 

   式即为我们我们熟知的
 
 梯形公式。写出MATLAB函数为
 

  
function
 
 
而若用区间中点 
 

   的“高度” 
 
 

   近似取代平均高度 
 
 

   ，则可导出
 
 中矩形公式（简称
 矩形公式）
 

 

  
 
 

  
 
 
其MATLAB函数为
 

  
function
 
 
总结：更一般地，我们可以在区间 
 

   上适当选取某些节点 
 
 

   ，然后用 
 
 

   的加权平均得到平均高度 
 
 

   的近似值，这样构造出来的求积公式具有下列通式
 
 

 

  
 
 

  
 
 
其中 
 

   称为
 
 求积节点， 
 

   称为
 
 求积系数，也称作伴随节点 
 

   的
 
 权。
 
这种数值积分方法通常称为机械求积，其优势在于将积分求值问题归结为被积函数值的计算，很适合在计算机上使用。
 
代数精度
 
由于数值积分方法是近似方法，为保证精度，自然希望我们的求积公式能够对“尽可能多”的函数准确成立，这便提出了代数精度的概念。
 

  如果某个求积公式对于次数不超过 
  

    的多项式均能够准确成立，但对于 
  
  

    次多项式就不准确成立，则称该求积公式具有
  
  
  

    次代数精度
  
  （或
  代数精确度）。
 
 
前面提到的梯形公式 
 

   和矩形公式 
 
 

   均具有一次代数精度。
 
 
求积公式的收敛性与稳定性
 

  在求积公式 
  

    中，若 
  
  

   
  
  

    其中 
  
  

    ，则称求积公式 
  
  

    是
  
  收敛的。
  
对任给 
  

    ，若 
  
  

    ，只要 
  
  

    就有 
  
  

   
  
  

    成立，其中 
  
  

    为计算 
  
  

    时产生误差 
  
  

    后实际得到的值，即 
  
  

    ，则称求积公式 
  
  

    是
  
  稳定的。
 
 
定理 
 

  
 
 若求积公式 
 

  
 
 中系数 
 

  
 
 ，则此求积公式是稳定的。
 
牛顿-柯特斯(Newton-Cotes)公式
 
将积分区间 
 

   划分为 
 
 

   等份，此时步长为 
 
 

   ，选取等距节点 
 
 

   构造出的插值型求积公式
 
 

 

  
 
 

  
 
 
称为牛顿-柯特斯(Newton-Cotes)公式，其中 
 

   称为
 
 柯特斯系数。可以通过下式确定
 

 

  
 
 

  
 
 
当 
   

     时， 
   
   

     这时的求积公式便是之前的梯形公式 
   
   

     。
   
当 
   

     时， 
   
   

    
   
   

     此时对应的求积公式便是
   
辛普森(Simpson)公式： 
   

    
   
   

     写成MATLAB函数为
   
 

  
function
 
 
当 
   

     时的牛顿-柯特斯公式则特别称为
   
柯特斯公式，其形式为 
   

    
   
   

     这里 
   
   

     。
   
 
其MATLAB代码为
 

  
function
 
 
复合求积公式
 
由于牛顿-柯特斯公式在 
 

   时不具有稳定性，故不可能再通过提高阶的方法来提高求积精度。
 
 复合求积法便是通过把积分区间分成若干个子区间（通常是等分），再在每个子区间上使用低阶求积公式的方法来提高精度的。
 
其实细心的同学可以发现，前面我所展示的MATLAB代码使用的便是对每个求积公式使用复合法完成的。常用的有复合梯形公式和复合辛普森公式。
 
龙贝格(Romberg)求积公式
 
前面介绍的复合求积方法可提高求积精度，如若精度仍不够，则可通过将步长逐次减半的方式来提高精度。如对复合梯形公式可导出其递推公式
 

 

  
 
 

  
 
 
其中 
 

   表示在 
 
 

   基础上步长 
 
 

   减半后的积分值。
 
 
定理 
 

  
 
  设 
 

   ，则有
 
 
 

 

  
 
 
 

   其中
 
  
 

   ，
 
 系数 
 

   与 
 
 
 

   无关。
 
 
 
从 
 

   式可以看出， 
 
 

   是 
 
 

   阶，若用 
 
 

   代替 
 
 

   ，有
 
 

  
 
 

  
 
 
再用 
 

   乘 
 
 

   式再减去 
 
 

   式后再除以 
 
 

   得
 
 

 

  
 
 

  
 
 
这里 
 

   是与 
 
 

   无关的系数。
 
 
从 
 

   式可看出，用近似积分值 
 
 

   ，其误差阶为 
 
 

   ，显然误差阶是提高了。类似这种将计算 
 
 

   的近似值的误差阶由 
 
 

   提高到 
 
 

   的方法称为
 
 外推算法，也称为
 理查德森(Richardson)外推算法。
 

  这是“数值分析”中一个重要的技巧，只要真值与近似值的误差能够表示成 
  

    的幂级数，如 
  
  

    式所示，都可以使用外推算法，提高精度。
  
 
 
龙贝格(Romberg)算法
 

 

   次外推加速为
 
 

 

  
 
 

  
 
 
余项为
 

 

  
 
 

  
 
 
此方法常称为理查德森外推加速方法。
 
设用 
 

   表示二分 
 
 

   次后求得的梯形值，且以 
 
 

   表示序列 
 
 

   的 
 
 

   次加速值，则由递推公式 
 
 

   可得
 
 

 

  
 
 

  
 
 
上次则称为龙贝格求积算法，计算过程如下：
 
取 
   

     ，求 
   
   

     。令 
   
   

     ( 
   
   

     记区间 
   
   

     的二分次数)。
   
求梯形值 
   

     ，即按递推公式 
   
   

     计算 
   
   

     。
   
求加速值。
若 
   

     (预先给定的精度)，则终止计算，并且 
   
   

     ；否则令 
   
   

     转 
   
   

     继续计算。
   
 
MATLAB代码为
 

  
function
 
 
自适应积分方法
 
复合求积方法通常适用于被积函数变化不太大的积分，如果在求积区间中被积函数变化很大，有的部分函数值变化剧烈，另一部分却变化平缓。这时如果统一将区间等分再用复合求积公式计算积分将会导致计算量很大，我们想实现的是在满足误差要求的前提下，对变化剧烈部分将区间细分，而平缓部分则可使用大步长，也即针对被积函数在区间上不同情形采用不同的步长，使得再满足精度前提下积分计算工作量尽可能小。其算法技巧是在不同区间上预测被积函数变化的剧烈程度确定相应步长。这就是自适应积分方法的基本思想。
 
下面是一个基于复合辛普森法的自适应积分算法的MATLAB代码：
 

  
function
 
 
下面我们计算积分 
 

   ，可输入下列语句调用上述MATLAB函数进行计算
 
 

  
clear

 
 
结果为
 

  
I 
 
 
高斯求积公式
 
下面研究带权积分 
 

   （ 
 
 

   为权函数）的求积公式
 
 

 

  
 
 

  
 
 
下面看定义
 

  如果求积公式 
  

    具有 
  
  

    次代数精度，则称其节点 
  
  

    为
  
  高斯点，相应公式 
  

    称为
  
  高斯型求积公式。
 
 
这里主要包括四种高斯型求积公式，即
 
高斯-勒让德求积公式
高斯-切比雪夫求积公式
高斯-拉盖尔求积公式
高斯-埃尔米特求积公式
 
多重积分
 
这里主要说一下多重积分计算的思路，归结为一句话，就是：多重积分化累次积分，再由里到外使用数值积分公式进行求积计算。

 
对于积分区间[a, b]，若
 

  
 
 
则成立
 

  
 
 
辛普森公式
 
辛普森公式可看作是改良的梯形公式。梯形公式是以直线逼近实际曲线，而辛普森公式则以二次曲线(即抛物线)逼近。
 

  
  
以二次曲线逼近实际曲线
 
 
根据辛普森公式可得各积分小区间内的积分值如下式所示：
 

  
 
 
其中
 

  
 
 
于是，全积分区间[a, b]的积分值为
 

  
 
 
为便于计算机计算，上式可写成
 

  
 
 
其中
 

  
 
 
C语言程序实现辛普森公式计算积分
 
计算积分
 

  
 
 
代码
 
//simpson.c
 
#include 
 
#include 
 
/*函数定义*/
 
double f(double x)
 
{
 
 return (1.0-x)*exp(-1*x);
 
}
 
/*辛普森公式计算积分*/
 
double simpson(
 
 double a, double b, int n)
 
{
 
 double h = (b-a)/n;
 
 double s0,s1=0.0,s2=0.0;
 
 double S;
 
 double x;
 
 int i;
 
 s0 = f(a) + f(b);
 
 for( i = 1; i < n; i += 2 )
 
 {
 
 x = a + h * i;
 
 s1 += f(x);
 
 }
 
 for( i = 2; i < n; i += 2 )
 
 {
 
 x = a + h * i;
 
 s2 += f(x);
 
 }
 
 S = h*(s0+4*s1+2*s2)/3;
 
 return S;
 
}
 
int main(void)
 
{
 
 //积分区间[a, b]
 
 double a = -1.0, b = 1.0;
 
 //区间分割数
 
 int N[5]={10,20,30,50,100};
 
 int i = 0;
 
 for( i = 0; i < 5; i++ )
 
 printf("[%3d等分]积分近似值:
 
 I = %.7lf

 

 

  
 
 
摘自各种资源，希望帮助到大家，欢迎关注我！
 
目录
 
导数公式
基本积分表
三角函数的有理式积分
初等函数
两个重要极限
三角函数公式
和差角公式
和差化积公式
倍角公式
半角公式
正玄余弦公式
莱布尼兹公式、中值定理、曲率
定积分应用相关公式
空间几何
平面方程
二次曲面
多元函数微分
多元复合函数及应用
隐函数求导
微分在几何上的应用
方向导数与梯度
多元函数极值
重积分及应用
柱面坐标和球面坐标
曲线积分
曲面积分
高斯公式
斯托克斯公式
常数项级数
审敛法
收敛判别
幂级数
三角级数
傅里叶级数
微分方程的相関概念
线性微分方程
 
导数公式
 

 

  
 
 

基本积分表
 

 

  
 
 
三角函数的有理式积分
 

 

  
 
 
初等函数
 

 

  
 
 
两个重要极限
 

 

  
 
 
三角函数公式
 

 

  
 
 
和差角公式
 

 

  
 
 
和差化积公式
 

 

  
 
 
倍角公式
 

 

  
 
 
半角公式
 

 

  
 
 
正玄余弦公式
 

 

  
 
 
莱布尼兹公式、中值定理、曲率
 

 

  
 
 

 

  
 
 
定积分应用相关公式
 

 

  
 
 
空间几何
 

 

  
 
 
平面方程
 

 

  
 
 
二次曲面
 

 

  
 
 
多元复合函数及应用
 

 

  
 
 
隐函数求导
 

 

  
 
 

 

  
 
 
微分在几何上的应用
 

 

  
 
 
方向导数与梯度
 

 

  
 
 
多元函数极值
 

 

  
 
 
重积分及应用
 

 

  
 
 
柱面坐标和球面坐标
 

 

  
 
 
曲线积分
 

 

  
 
 

 

  
 
 
曲面积分
 

 

  
 
 
高斯公式
 

 

  
 
 
斯托克斯公式
 

 

  
 
 
常数项级数
 

 

  
 
 
审敛法
 

 

  
 
 
收敛判别
 

 

  
 
 
幂级数
 

 

  
 
 

 

  
 
 
三角级数
 

 

  
 
 
傅里叶级数
 

 

  
 
 

 

  
 
 
线性微分方程
 

 

  
 
 

 

  
 
 

 

  
 

 
写到后面的时候忽然发现自己好多小点没有写，这里赶紧亡羊补牢一下：
 
pi：在matlab里代表圆周率
 
abs（a）对a求绝对值
 
mod（a,b）：a除b的余数，emmm C语言中一般用a%b
 
反三角函数：
 
反正弦：asin（）
 
反余弦：acos（）
 
反正切：atan（）
 

 
 上一节啰嗦了一下matlab的操作界面，这节我们正式开始学习一下怎么用matlab
 
基本运算（我们先假设几个常数a,b,c,d...）：
 
 加法：a+b
 
 减法：a-b
 
 乘法：分为点乘和叉乘，叉乘就是a*b，点乘是a.*b（matlab就是这么直接...）
 
 这里补充一点，数值运算点乘叉乘是一样的，矩阵才分点乘叉乘。
 
 除法：也是分为点除和叉除，即a/b和a./b
 
 乘方：a的b次方就是a^b，注意自然对数e的n次方是exp(n),其中n为常数
 
 对数：log(a）指以e为底，真数为a(也就是lna）；
 
 此外还有以2为底log2(a),以10为底log10(a）
 
 其它的对数需要用换底公式，比如：
 

 

   前面的是想求的，写在程序里的就是等号后面的式子
 
 
 开方：a的平方根就是sqrt(a)，当然用a^0.5也可以
 
其它：
 
赋值a=b就是把b的值赋给a
 
 特殊的一点：a=a+c就是把a的值加上c
 
 这节主要记录了matlab的一些基本运算，下一节简单讲讲matla的一些矩阵操作。
 
欢迎喜欢的朋友点赞收藏啊：）