数学建模与科学计算第三章.doc

第3章湘江流量估计模型——数值积分法 3.1湘江水流量估计的实际意义 水流量是水文特征的一个重要指标，而水文特征对于水资源的合理利用，防洪以及抗旱具有指导性的作用，因此，湘江水流量估计对于湘江流域的社会经济和人民生活具有重大的影响．现根据实际测量得到湘江某处河宽700m，其横截面不同位置某一时刻的水深如表3-1所示．若此刻湘江的流速为0.5m/s，试估计湘江此刻的流量． 表3-1 湘江某处横截面不同位置的水深数据 单位：m x 0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 h(x) 4.2 5.9 5.8 5.2 4.5 5.7 5 5.5 4.8 5.9 4.1 5.1 4.6 5.7, 4.7 要计算湘江水流量就需要知道其横截面面积，如果知道此处江的水深曲线函数，则其横截面面积为．但是在实际中是不可能精确得到的，那么怎样求出足够高精度的横截面面积的近似值． 3.1.1数值求积的必要性 在高等数学中，有牛顿-莱布尼兹（Newton-Leibniz）公式 , 其中是的一个原函数.此公式常被用来计算定积分．但是，在工程技术和科学研究中，常常遇到如下情况： (1) 的结构复杂，求原函数很困难； (2) 的原函数不能用初等函数表示； (3) 的精确表达式不知道，只给出了一张由实验提供的函数表． 对于这些情况，要计算积分的精确值都是十分困难的，这就要求建立积分的近似计算方法．此外，积分的近似计算又为其他一些数值计算，如微分方程数值解、积分方程数值解等，提供了必需的基础． 3.1.2 构造数值求积公式的基本方法 可以从不同的角度出发，通过各种途径来构造数值求积公式．但常用的一个方法是，利用插值多项式来构造数值求积公式．具体做法如下： 在积分区间上取一组点的次插值多项式 ， 其中为次Lagrange插值基函数．用近似代替被积函数，则得 （3.1.1） 若记 （3.1.2）则得数值求积公式 （3.1.3）形如式（3.1.3）的求积公式称为机械求积公式，其中称为求积节点，称为求积系数．若求积公式（3.1.3）中的求积系数是由式（3.1.2）确定的，则称该求积公式为插值型求积公式．本章主要讨论插值型求积公式． 3.1.3求积公式的余项 积分的真值与由某求积公式给出的近似值之差，称为该求积公式的余项，记作．例如，求积公式（3.1.3）的余项为 . 如果求积公式（3.1.3）是插值型的，则由式（3.1.3）知 ， 于是，由插值余项公式得 （3.1.4）其中． 3.1.4求积公式的代数精度 为了使一个求积公式能对更多的积分具有良好的实际计算意义，就应该要求它对尽可能多的被积函数都能准确地成立．在数值计算中，常用代数精度这个概念来描述这种准确的程度． 定义3.1 若求积公式  对任意不高于次的代数多项式都准确成立，而对于却不能准确成立，则称该公式的代数精度为． 例3.1证明梯形公式 （3.1.5）的代数精度． 证在几何上就是用梯形面积近似代替曲边梯形面积，如图3-1所示． 图3-1 梯形求积公式几何直观示意图 (1) 当时，在式（3.1.5）中： 左端=， 右端=， 左端=右端， 这表明求积公式（3.1.5）对是准确成立的. (2) 当时，在式（3.1.5）中： 左端=， 右端=， 左端=右端， 这表明求积公式（3.1.5）对也是准确成立的. 综上所述，容易看出求积公式（3.1.5）对函数和的任一线性组合（不高于一次的代数多项式）都准确成立，故公式（3.1.5）的代数精度至少等于1． (3) 当时，在式（3.1.5）中： 左端， 右端， 左端(右端（设）， 故由定义3.1可知，梯形公式（3.1.5）的代数精度． 显然，一个求积公式的代数精度越高，它就越能对更多的被积函数准确（或较准确）地成立，从而具有更好的实际计算意义．由插值型求积公式的余项（3.1.4）容易得出下面的定理. 定理3.1含有个节点的插值型求积公式（3.1.3）的代数精度至少为. 3.2 数值求积的牛顿-柯特斯方法 在3.1中，介绍了插值型求积公式及其构造方法．在实际应用时，考虑到计算的方便，常将积分区间等分，并取分点为求积节点．这样构造出来的插值型求积公式就称为牛顿-柯特斯（Newton-Cotes）公式．本节在介绍一般Newton-Cotes公式的基础上，介绍几个常用的Newton-Cotes公式以及这些公式在实际计算时的用法． 3.2.1 Newton-Cotes公式 若将积分区间作等分，取分点作为求积节点，并作变量替换，则插值型求积公式(3.1.3)的系数由式（3.1.2）可得