二重数值积分算法的研究与实现.docVIP

 本科毕业论文 题目： 二重数值积分算法的研究与实现 学院： 数学与计算机科学学院 班级： 数学与应用数学2004级四班 姓名： 　　 指导教师： 职称： 副教授 完成日期： 2008 年 5 月 30 日  二重数值积分算法的研究与实现 摘 要： 二重积分与定积分有着密切的联系，它是一元函数定积分在二元函数的推广，本文通过已有的单积分的Newton-Cotes公式和 Gauss型求积公式，讨论二元数值积分的Newton-Cotes公式和 Gauss型求积计算方法.推而广之，该两种二重数值积分方法同样适用于多元数值积分的计算．本文并给出了MATLAB实现方法，减轻了简单积分问题的繁琐计算量. 关键词：Newton-Cotes公式；Gauss型求积公式；MATLAB；二重数值积分 目 录 1 单重数值积分的构造 1 1.1 Newton-Cotes公式计算法 1 1.2 Gauss型求积公式计算法 3 2 二重数值积分公式的构造 3 2.1 二元Newton-Cotes积分方法 3 2.1.1二元Newton-cotes型积分方法的算法实现 6 2.1.2算法分析 6 2.2 比较 6 3 MATLAB的实现方法 7 计算定积分的值是工程技术计算中常常遇见的一个问题．关于单积分的求积公式以及误差估计式，研究是较为深入的，而二重积分与定积分有着密切的联系，它是一元函数定积分在二元函数的推广．计算二重积分的基本方法是化为二次定积分来计算．本文将对单积分中Newton-Cotes公式进行推广，用来计算二重积分以及二重以上的数值积分，解决函数的原函数无法计算出来的问题，同理也可以把 Gauss型求积公式进行推广． 二重积分在科学计算中起着重要作用，但因二重积分计算过程中涉及到多元之间的隐函数关系式，所以它是一门很难弄懂学透的数学课程，遇到较难一点的计算题，就不知道如何动手．如果我们利用MATLAB函数作为辅助工具，这种束手无策的问题就容易多．MATLAB是一门科学计算软件，具有输入简捷，运算高效，内容丰富，扩展性强等特点．目前以成为大学数学教育和科学研究中必不可少的工具，它大大提高了课程教学，解题作业和分析研究的效率．在二重积分学中应用MATLAB函数，可以帮助直接计算二重积分． 1 单重数值积分的构造 对区间上定义的黎曼可积函数构造定积分的近似计算公式，若直接计算定积分有困难，可以利用插值公式，先构造一个多项式逼近被积函数，并以在区间上的积分代替在区间上的积分即 ，而多项式的积分则是容易计算的，下面就根据这一想法构造计算定积分的多种近似计算公式． 1.1 Newton-Cotes公式计算法 把区间 进行n等分，其分点为： 过这n+1个节点，构造一个n 次多项式 ，其中，用代替被积函数则有：= == ， 其中， 使用Newton-Cotes公式的关键是计算系数，用变量替换于是而 这样 引进记号，则， 这时是不依赖与函数和区间的常数，可以事先计算出来叫Newton-Cotes系数,下面给出n取值1-6的Newton-Cotes系数表Ⅰ 表Ⅰ n 1 2 3 4 5 6 在n=1时 为梯形求积公式 在n=2时为Simpson公式（抛物线求积公式） 1.2 Gauss型求积公式计算法 对于一元的数值积分，在节点数目固定为n的条件下，适当地选择节点位置和相应的系数，通过积分插值公式，可得如下的插值型求积公式 ． ⑴ 这里是实轴上的有限或无限区间；是一个确定的函数，称为权函数；称为求积公式的系数，它只依赖于权函数和积分区间，而与求积函数无关；称为求积公式的基点． 若基点取作上带权函数的正交多项式的零点，则求积公式⑴称为Gauss型求积公式，它的基点称为Gauss型基点． 众所周知，n个基点的Gauss 型求积公式具有2n-1次代数精确度．Gauss 型求积公式因具有基点少精度高，数值稳定等显著特点而经常被采用． 2 二重数值积分公式的构造 将上述单重数值积分公式推广到二重积分的情形：考虑二重积分不妨设积分区域是上的连续可微函数，在上充分光滑（连续）则． 2.1 二元Newton-Cotes积分方法