高精度数值积分公式的构造及其应用 数学专业(设计) 学位论文
- 格式:doc
- 大小:2.75 MB
- 文档页数:43
一种构建高精度求积公式的新策略郑华盛;徐伟【摘要】In this paper, a new general tactics of improving and constructing high-order accurate quadrature formula is presented to use low-order quadrature formula as a elementary building block by applying its expression of remainder term and concepts of algebraic accuracy. Similarly, finite number improving can be carried out. Finally, several examples of common low-order quadrature for- mulae are given to show higher validation of this method.%以低阶求积公式为基本模块,基于它的余项表达式及代数精度概念,提出了一种改进和构造高精度求积公式的普适性新策略。
该方法可实现求积公式的有限次改进。
最后,应用于几个常用低阶求积公式,以验证本文方法的有效性。
【期刊名称】《江西科学》【年(卷),期】2012(030)005【总页数】4页(P559-561,602)【关键词】求积公式;代数精度;余项【作者】郑华盛;徐伟【作者单位】南昌航空大学数学与信息科学学院,江西南昌330063;南昌航空大学数学与信息科学学院,江西南昌330063【正文语种】中文【中图分类】O241.4科学与工程实际中的许多问题最后都归结为定积分的计算问题。
由于被积函数过于复杂或原函数不易求得等原因,实际计算中大多得不到定积分的精确解,此时,就必须借助数值求积公式计算定积分的近似值,因而研究高精度数值求积公式具有很重要的意义。
数值积分的理论及其应用研究数值积分的多种问题及其在现代工程中的广泛应用张冲聪(西安文理学院数学系陕西西安 710065)摘要:数值积分的多种问题及其在现代工程中的广泛应用的探讨是计算数学的一个重要课题,数值积分是数学上的重要课题之一,是数值分析中的重要内容之一,也是数学的研究重点。
并在实际问题及应用中有着广泛的应用。
常用于科学与工程的计算中,如涉及到积分方程,工程计算,计算机图形学,金融数学等应用科学领域都有着相当重要的应用,所以研究数值积分问题有很重要的意义。
研究方法有插值法和抽样插值法等。
当然大家都知道计算积分可以借助原函数和查找积分表,但是,用这些方法只能解决很狭隘的一类积分,而且在计算的过程中,肯定会产生误差,我们要想法子使得误差尽可能的小。
因此,数值积分的公式应满足:计算简单,误差小,代数精度高等。
首先,我们通过构造函数并运用罗必达法则探讨数值积分中的矩形公式,梯形公式吧,抛物线公式,高斯公式的渐进性质。
结果表明,当积分区间的长度趋于零时,不但可以确定求积公式余项中的中介点的位置,还可以得到与之相应的修正公式,而且通过数值试验还能发现经过修正后的求积公式具有较高的代数精确度,我们可以通过构造函数运用分部积分的方法得到矩形公式,梯形公式吧,抛物线公式,高斯公式的推广并结合一些例子。
关键字:数值积分,矩形公式,梯形公式,抛物线公式,高斯公式近些年来,有关数值积分的研究已经成为一个很活跃的研究领域,所以研究数值积分有很重要的意义。
设f是闭区间[]ba,上某一给定的可积函数,现在要计算定积分⎰dtxf)(,我们可以借助原函数,或借助函数逼近的方法来计算,对于不熟悉的我们也可以借助参考积分表。
但都有一定的局限性,由于许多函数的无定积分无法用简单的函数表达出来,如一些离散点上的函数。
在微积分理论中,我们知道了牛顿—莱布尼茨(N ewtou_Leibniz)公式⎰-=)()()(aFbFdxxf其中f(x)在闭区间[]ba,上连续,F(x)是被积函数f(x)的某一个原函数,但是对于很多实际问题都无能为力。
(2012 届)本科毕业论文(设计)题目:柯西积分定理与柯西积分公式的由来及其应用学院:教师教育学院专业:数学与应用数学(师范)班级:数学082学号:姓名:指导教师:完成日期:教务处制诚信声明我声明,所呈交的论文(设计)是本人在老师指导下进行的研究工作及取得的研究成果。
据我查证,除了文中特别加以标注和致谢的地方外,论文(设计)中不包含其他人已经发表或撰写过的研究成果,也不包含为获得______或其他教育机构的学位或证书而使用过的材料。
我承诺,论文(设计)中的所有内容均真实、可信。
论文(设计)作者签名:签名日期:年月日授权声明学校有权保留送交论文(设计)的原件,允许论文(设计)被查阅和借阅,学校可以公布论文(设计)的全部或部分内容,可以影印、缩印或其他复制手段保存论文(设计),学校必须严格按照授权对论文(设计)进行处理,不得超越授权对论文(设计)进行任意处置。
论文(设计)作者签名:签名日期:年月日柯西积分定理与柯西积分公式的由来及其应用王莉莉(嘉兴学院数学与信息工程学院)摘要:复变函数是综合性大学或师院类院校理工专业的必修课,是实变函数微积分的推广和发展.其中柯西积分定理和柯西积分公式是复变函数理论的基础,是研究复变函数理论的关键.它的核心内容是柯西积分定理,即解析函数沿围线的积分值为零.本文研究了柯西积分定理和柯西积分公式的相关概念、证明、推广及在代数基本定理证明、实积分计算中的应用,论述了柯西积分定理与复变函数的积分有着密切的联系,利用柯西积分定理很容易导出著名的柯西积分公式,还对留数定理作了简要介绍,利用留数定理可以分别得到复变函数中的柯西积分定理、柯西积分公式和高阶导数公式.关键词:复变函数;柯西积分定理;柯西积分公式;留数定理Cauchy Integral Theorem and Cauchy Integral Formulas ofthe Origin and its ApplicationWanglili(College of Mathematics and Information Engineering , Jiaxing University)Abstract:Complex-variable function is a comprehensive university or institute of technology of normal colleges and universities professional required courses. It is real veriables function of the promotion and development of calculus. One cauchy integral theorem and cauchy integral formula is a complex function theory foundation. It is the key of sresearching complex function theory. One of its important contents is the cauchy integral theorem, which says that the integral along a contour of an analytic function is zero. This paper studies the cauchy integral theorem and cauchy integral formula related concepts and prove, promotion and in algebra fundamental theorem of integral proof, in the calculation of the application. It discusses the closely related of the cauchy integral theorem and cauchy complex functions. The famous cauchy integral formula can follows easily from the cauchy integral theorem. Also residue theorem are briefly introduced. Use of residue theorem can get the complex functions respectively cauchy integral theorem and cauchy integral formulas and high derivatives formula.Key words:Complex-variable function;cauchy integral theorem;cauchy integral formula;residue theorem目录1 绪论 (1)1.1 研究背景 (1)1.1.1 复变函数概况 (1)1.1.2 复积分的定义 (2)1.1.3 柯西积分定理的引入 (3)1.2 本文的研究工作 (4)1.3 本文的未来工作 (4)2 柯西积分定理 (5)2.1 柯西积分定理 (5)2.2 柯西积分定理的证明 (5)2.3 柯西积分定理的推广 (6)2.4 柯西积分定理的应用 (9)3 柯西积分公式 (12)3.1 柯西积分公式 (12)3.2 柯西积分公式的证明 (12)3.3 柯西积分公式的推广 (13)3.4 柯西积分公式的应用 (14)4 复变函数积分之间的关系 (18)4.1 柯西积分定理与柯西积分公式的关系 (18)4.2 复变函数积分与留数定理的关系 (19)参考文献 (22)1 绪论1.1 研究背景在18世纪后半叶到19世纪初,开始了复函数的偏导数与积分性质的探索.复分析真正作为现代分析的一个研究领域是在19世纪建立起来的,主要奠基人:柯西、黎曼和魏尔斯特拉斯.柯西建立了复变函数的微分和积分理论.1814年、1825年的论文《关于积分限为虚数的定积分的报告》建立了柯西积分定理,1826年提出留数概念,1831年获得柯西积分公式,1846年发现积分与路径无关定理[3].1.1.1 复变函数概况复数的概念起源于求方程的根,在二次、三次代数方程的求根中就出现了负数开平方的情况.在很长时间里,人们对这类数不能理解.但随着数学的发展,这类数的重要性就日益显a ,其中i是虚数单位.以复数作为自变量的函数就叫做复现出来.复数的一般形式是:bi变函数,而与之相关的理论就是复变函数论.解析函数是复变函数中一类具有解析性质的函数,复变函数论主要就是研究复数域上的解析函数,因此通常也称复变函数论为解析函数论.复变函数论的全面发展是在十九世纪,就像微积分的直接扩展统治了十八世纪的数学那样,复变函数这个新的分支统治了十九世纪的数学.当时的数学家公认复变函数论是最丰饶的数学分支,并且称为这个世纪的数学享受,也有人称赞它是抽象科学中最和谐的理论之一.复变函数中的许多概念、理论和方法是实变函数在复数领域内的推广和发展,因而它们之间有着许多的相似之处.但是,复变函数又有与实变函数不同之点,它是数学分析在研究领域的扩展.在我们学习中,要勤于思考,善于比较,既要注意共同点,更要弄清不同点.这样,才能抓住本质,融会贯通.复变函数论主要包括单值解析函数理论、黎曼曲面理论、几何函数论、留数理论、广义解析函数等方面的内容.如果当函数的变量取某一定值的时候,函数就有一个唯一确定的值,那么这个函数解就叫做单值解析函数,多项式就是这样的函数.复变函数研究多值函数,黎曼曲面理论是研究多值函数的主要工具.由许多层面安放在一起而构成的一种曲面叫做黎曼曲面.利用这种曲面,可以使多值函数的单值枝和支点概念在几何上有非常直观的表示和说明.对于某一个多值函数,如果能作出它的黎曼曲面,那么,函数在黎曼曲面上就变成单值函数.黎曼曲面理论是复变函数论和几何间的一座桥梁,能够使我们把比较深奥的函数的解析性质和几何性质联系起来.近来,关于黎曼曲面的研究还对另一门数学分支拓扑学产生了比较大的影响,逐渐地趋向于讨论它的拓扑性质.复变函数论中用几何方法来说明、解决问题的内容,一般叫做几何函数论,复变函数可以通过共形映像理论为它的性质提供几何说明.导数处处不是零的解析函数所实现的映像都都是共形映像,共形映像也叫做保角变换.共形映像在流体力学、空气动力学、弹性理论、静电场理论等方面都得到了广泛的应用.留数理论是复变函数论中一个重要的理论.留数也叫做残数,它的定义比较复杂.应用留数理论对于复变函数积分的计算比起线积分计算方便.计算实变函数定积分,可以化为复变函数沿闭回路曲线的积分后,再用留数基本定理化为被积分函数在闭合回路曲线内部孤立奇点上求留数的计算,当奇点是极点的时候,计算更加简洁.把单值解析函数的一些条件适当地改变和补充,以满足实际研究工作的需要,这种经过改变的解析函数叫做广义解析函数.广义解析函数所代表的几何图形的变化叫做拟保角变换。
编号学士学位论文牛顿-柯特斯求积公式的代数精度研究学生姓名:买买提克热木·亚森学号:20080101025系部:数学系专业:信息与计算科学年级: 2008-5班指导教师:阿米娜·沙比尔完成日期: 2013年4月20日本文主要对求数值积分公式的代数精度进行探讨。
首先描述了数值积分的矩形法,梯形法,插值求积公式等求积方法的基本思路和代数精度概念,进行了余项估计。
然后重点讨论对牛顿—柯特斯求积公式当15n 的情形及其代数精度.最后用数值例题验证牛顿—柯特斯求积公式的代数精度的重要性。
关键词:数值积分;梯形公式;代数精度;牛顿-柯特斯公式;2目 录摘要 (1)引言 (2)1基本概念 (4)1.1代数精度的概念 (4)1.3 插值型的求积公式 (7)1.4 求积公式的余项 (8)1.5 求积公式的收敛性和稳定性 (9)2. 牛顿—柯特斯公式 (10)2.1 柯特斯系数与辛普森公式 (10)2.2偶阶求积公式的代数度 (18)2.3 辛普森公式的余项 (19)参考文献 (23)致谢 (24)2引言1问题的提出实际问题当中常常需要计算积分,有些数值方法,如微分方程和积分方程的求解,也都和积分计算相联系.依据人们所熟知的微积分基本定理,对于积分()b af x dx I =⎰,只要找到被积函数)(x f 的原函数)(x F ,便有下列牛顿—莱布尼茨(Newton —Leibniz )公式()()()ba f x dx Fb F a =-⎰但实际使用这种求积方法往往有困难,因为大量的被积函数,比如xx )sin()0(≠x ,2x e -等,其原函数不能用初等函数表达,故不能用上述公式计算,即使能求得原函数的积分有时计算也十分困难.例如对于被积函数611)(xx f +=,其原函数 c x x x x x x x x F ++-+++-+=1313ln 341)1arctan(61arctan 31)(22 计算)(),(b F a F 仍然很困难.另外,当)(x f 是有测量或数值计算给出的一张数据表时,牛顿—莱布尼茨公式也不能直接运用.因此很有必要研究计算方便,计算量少,精度高而且稳定的数值积分方法.2我的想法积分中值定理告诉我们,在积分区间[]b a ,内存在一点ξ,成立()()()b a f x dx b a f ζ=-⎰就是说,低为a b -而高为)(ζf 的矩形面积恰等于所求曲边梯形的面积I ,如图4.1 ,问题在于点ξ的具体位置一般是不知道的,因而难以准确算出)(ζf 的值,我们将)(ζf 称为区间[]b a ,上的平均高度。
目录1.数值积分的历史 (3)1.1数值积分的起源 (3)1.2数值积分的经典方法 (3)1.3著名数学家—辛普森 (3)2.常用的数值积分方法 (4)2.1插值型求积公式 (4)2.2Newton-Cotes公式 (5)2.2.1梯形公式 (5)2.2.2 Simpson 公式 (5)2.2.3柯特斯公式 (6)2.2.4 牛顿—柯特斯公式 (6)2.3复合求积公式 (7)2.3.1复合梯形公式 (8)2.3.2 复合的辛普森公式 (8)2.3.3 复合的柯特斯公式 (9)2.4逐次分半技术与龙贝格公式 (10)2.4.1梯形公式的递推化 (10)2.4.2龙贝格公式 (10)2.5高斯型求积公式 (12)2.5.1高斯—勒让德求积公式 (14)2.5.2高斯—切比雪夫求积公式 (14)G求积公式 (15)2.5.3auss-Laguerre2.6奇异积分的数值计算 (15)2.6.1反常积分的计算 (15)2.6.2无穷区间积分的计算 (16)2.7振荡函数的积分 (17)2.7.1分部积分公式 (18)2.7.2 Filon法 (19)3.数值积分在Matlab中的应用 (20)参考文献 (26)1.数值积分的历史1.1数值积分的起源数值积分是求定积分的近似值的数值方法,即用被积函数的有限个抽样值的离散或加权平均近似值代替定积分的值。
求某函数的定积分时,在多数情况下,被积函数的原函数很难用初等函数表达出来。
另外,许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数,对这类函数的定积分,也不能用不定积分的方法求解。
由于以上原因,数值积分的理论与方法一直是计算数学研究的基本课题。
对微积分学做出杰出贡献的数学大师,如牛顿,欧拉,高斯等人也在数值积分这个领域做出了各自的贡献,并奠定了它的理论基础。
1.2数值积分的经典方法构造数值积分最通常的方法是用积分区间上的n次插值多项式代替被积函数,由此导出的求积公式为插值型求积公式。
二重数值积分公式的高精度对偶公式及其应
用
1二重数值积分公式
二重数值积分公式是一种用于计算复杂数学函数的结果的数学方法。
它可以更加准确和准确地计算函数值,是计算复杂数学函数值的一个非常有用的方法。
2高精度对偶公式
高精度对偶公式是二重数值积分公式的一种重要选择。
它使用比较复杂的步骤来计算最终结果,但是其总体效率可以比普通方法大大提高。
这种方法可以使计算非常快速,甚至在高精度应用中也可以达到最佳效果。
3应用
高精度对偶公式的应用非常广泛,它可以应用于多种不同的数学和物理问题,其中包括惯性导航、频率响应、热传导和声学问题等等。
此外,它还可以用来解决一维或多维区域的分析问题,使结果的精度和精确度得到更大的提高。
这种技术也被广泛用于科学研究和工程领域,对分析复杂数学函数具有重要的作用。
高精度数值求积公式的构造
龙爱芳;胡军浩
【期刊名称】《河南师范大学学报:自然科学版》
【年(卷),期】2012(40)5
【摘要】构造了一类求积公式,此类求积公式只需计算节点上的函数值,避免计算导数值,它比复化梯形公式的计算量小,但收敛阶却大大的提高了.利用这类积分公式进行计算可以得到十分精确的结果.
【总页数】4页(P31-34)
【关键词】数值积分;导数;插值
【作者】龙爱芳;胡军浩
【作者单位】中南民族大学数学与统计学学院
【正文语种】中文
【中图分类】O241.4
【相关文献】
1.一个高精度数值求积公式的重构及其渐近性 [J], 管林挺;郑华盛
2.二维弱奇异积分高精度数值求积公式的构造 [J], 曾光;黄晋;雷莉;宁德圣
3.带重结点的三角求积公式的迭代构造 [J], 蔡晖;金国祥
4.一种构建高精度求积公式的新策略 [J], 郑华盛;徐伟
5.含Cauchy核奇异积分高精度求积公式 [J], 陈传希;金国祥
因版权原因,仅展示原文概要,查看原文内容请购买。
数值积分方法在工程计算中的应用在工程计算领域,准确求解各种积分问题是至关重要的。
数值积分方法作为一种有效的工具,为解决复杂的积分计算提供了可靠的途径。
它在诸如力学、电学、热学等多个工程学科中都有着广泛而重要的应用。
数值积分方法的基本思想是通过对被积函数在给定区间内进行离散化处理,用一系列简单的函数值之和来近似积分值。
常见的数值积分方法包括梯形法则、辛普森法则、高斯求积公式等。
梯形法则是一种简单直观的数值积分方法。
它将积分区间划分为若干个小梯形,通过计算这些小梯形的面积之和来近似积分值。
虽然梯形法则的精度相对较低,但在计算简单、对精度要求不高的情况下,它仍然能够提供较为快速的估算结果。
辛普森法则则在梯形法则的基础上进行了改进。
它将积分区间划分为若干个小段,然后用二次曲线来拟合被积函数,从而提高了积分的精度。
辛普森法则在处理一些光滑的函数时,能够给出较为准确的结果。
高斯求积公式则是一种具有较高精度的数值积分方法。
它通过选择适当的积分节点和权重系数,使得积分的精度得到显著提高。
然而,高斯求积公式的计算相对复杂,需要根据具体的积分区间和精度要求进行选择和计算。
在工程实际中,数值积分方法有着众多的应用场景。
例如,在结构力学中,计算梁的弯曲变形、结构的内力分布等问题时,常常需要对各种复杂的函数进行积分。
通过数值积分方法,可以快速有效地得到近似解,为结构的设计和分析提供重要依据。
在电气工程中,求解电路中的电流、电压、功率等参数时,也经常会遇到积分运算。
特别是在处理非线性电路或时变电路时,数值积分方法能够帮助工程师准确地计算出相关的电气量,从而保证电路的正常运行和设计的可靠性。
在热学工程中,计算热传递过程中的热量分布、温度变化等问题时,数值积分同样发挥着重要作用。
例如,在研究物体的冷却过程或热交换器的性能时,需要对热传导方程进行积分求解,数值积分方法为解决这类问题提供了有效的手段。
此外,在流体力学中,计算流体的流量、压力分布等参数,以及在控制系统的分析和设计中,计算系统的响应、稳定性等指标时,数值积分方法都不可或缺。
数值积分方法与应用数值积分方法是一种数值计算技术,用于计算函数在给定区间上的定积分。
在实际应用中,我们经常会遇到无法通过解析方法求解的定积分,这时候就可以借助数值积分方法来进行近似计算。
本文将介绍数值积分的基本原理、常用方法以及在实际问题中的应用。
一、基本原理在介绍数值积分方法之前,我们先来回顾一下定积分的几何意义。
对于函数f(x),在区间[a, b]上的定积分∫[a, b]f(x)dx表示函数f(x)在区间[a, b]上与x轴之间的面积。
当函数f(x)是非常复杂的时候,我们往往无法通过解析方法求解定积分,这时候就需要借助数值积分方法进行近似计算。
数值积分方法的基本原理是将积分区间分割成若干个小区间,然后在每个小区间上选取一个节点进行函数值的采样,最后通过对这些采样值的加权和来近似表示定积分的值。
常用的数值积分方法包括Newton-Cotes公式、Gauss求积法等。
二、常用方法1. Newton-Cotes公式Newton-Cotes公式是最简单的数值积分方法,其基本思想是将积分区间均匀分割成若干个小区间,然后在每个小区间上取若干个节点进行函数值的采样。
最常见的Newton-Cotes公式为梯形公式和Simpson 公式。
梯形公式是将积分区间[a, b]分割成n等分,然后在相邻两个节点上计算函数值,最后通过梯形面积的加权和来近似表示定积分的值。
Simpson公式是将积分区间[a, b]分割成2n等分,然后在每个子区间的两个端点和中点上计算函数值,最后通过三次多项式的插值来近似表示定积分的值。
2. Gauss求积法Gauss求积法是通过选取一定的节点和权重来提高数值积分方法的精度。
其基本思想是在给定区间上选取一些特定的节点和权重,然后通过这些节点和权重的组合来构造一个更高阶的数值积分公式。
Gauss求积法的优点是可以通过适当选择节点和权重来提高数值积分的精度,适用于高阶多项式的数值积分。
三、应用案例数值积分方法在科学计算、工程建模等领域有着广泛的应用。
毕业设计(论文)设计(论文)题目:数值积分算法与MATLAB实现摘要在求一些函数的定积分时,由于原函数十分复杂难以求出或用初等函数表达,导致积分很难精确求出,只能设法求其近似值,因此能够直接借助牛顿-莱布尼兹公式计算定积分的情形是不多的。
数值积分就是解决此类问题的一种行之有效的方法。
积分的数值计算是数值分析的一个重要分支;因此,探讨近似计算的数值积分方法是有着明显的实际意义的。
本文从数值积分问题的产生出发,详细介绍了一些数值积分的重要方法。
本文较详细地介绍了牛顿-科特斯求积公式,以及为了提高积分计算精度的高精度数值积分公式,即龙贝格求积公式和高斯-勒让德求积公式。
除了研究这些数值积分算法的理论外,本文还将这些数值积分算法在计算机上通过MATLAB软件编程实现,并通过实例用各种求积公式进行运算,分析比较了各种求积公式的计算误差。
【关键词】数值积分牛顿-科特斯求积公式高精度求积公式MATLAB软件ABSTRACTWhen the solution of the definite integral of some function values,because the original function is very complex and difficult to find the elementary function expression, the integral is difficult to accurately calculate, only managed to find the approximate value, and the case is small that allows to direct interface with the Newton - Leibniz formula to calculate the definite integral. Numerical integration is an effective method to solve such problems. The numerical integration is an important branch of numerical analysis; therefore, exploring the approximate calculation of the numerical integration method has obvious practical significance. This article departure from the numerical integration problem, described in detail some important numerical integration methods.This paper has introduced detail the Newton - Coates quadrature formula, and in order to improve the calculation accuracy of numerical integration formulas, More precise formulas have Romberg quadrature formulas and the Gauss - Legendre quadrature formula. In addition to the study of these numerical integration algorithm theory, the article also involve what these numerical integration algorithm be programmed by matlab software on the computer, and an example is calculated with a variety of quadrature formulas, finally analysis and comparison to various quadrature formulas calculation error.【Key words】Numerical integration Newton-Cotes quadratureformula High-precisionquadrature formula Matlab software目录前言 ..................................................................第一章牛顿-科特斯求积公式..............................................第一节数值求积公式的构造...........................................第二节复化求积公式.................................................第三节本章小结.....................................................第二章高精度数值积分算法...............................................第一节梯形法的递推.................................................第二节龙贝格求积公式...............................................第三节高斯求积公式.................................................第四节高斯-勒让德求积公式..........................................第五节复化两点高斯-勒让德求积公式 ..................................第六节本章小结.....................................................第三章各种求积公式的MATLAB编程实现与应用 ...........................第一节几个低次牛顿-科特斯求积公式的MATLAB实现...................第二节复化求积公式的MATLAB实现..................................第三节龙贝格求积公式的MATLAB实现................................第三节高斯-勒让德求积公式的MATLAB实现...........................第五节各种求积算法的分析比较 .......................................第六节本章小结.....................................................结论 ..................................................................致谢 ..................................................................参考文献 ................................................................附录 ..................................................................一、英文原文......................................................二、英文翻译......................................................前言对于定积分,在求某函数的定积分时,在一定条件下,虽然有牛顿-莱布里茨公式可以计算定积分的值,但在很多情况下的原函数不易求出或非常复杂。
基于Hermite插值的高精度数值积分公式龙爱芳;胡军浩【摘要】构造Hermite插值多项式,得到插值型求积公式.分析积分中值定理中间点的渐近性,得到具有更高精度的数值求积公式.对数值积分公式中的导数进行处理,最终得到不用计算导数值,只需计算节点处函数值的高精度数值求积公式.【期刊名称】《华侨大学学报(自然科学版)》【年(卷),期】2013(034)003【总页数】4页(P349-352)【关键词】Hermite插值;中间点;渐近性;高精度【作者】龙爱芳;胡军浩【作者单位】中南民族大学数学与统计学学院,湖北武汉430074;中南民族大学数学与统计学学院,湖北武汉430074【正文语种】中文【中图分类】O241随着科学的进步,计算机技术的发展,很多领域涉及到定积分的计算,因此研究高精度的数值积分公式是有实际意义的.数值积分常见的有梯形公式和Simpson公式,它们的计算虽无需提供导数值,但代数精度不高.梯形公式有一次代数精度,Simpson公式有3次代数精度[1-6].文献[7-8]虽给出一个高精度的数值求积公式,但必须提供求积节点的一阶导数值.文献[9]给出了Newton-cotes求积公式的渐近性,虽可大大提高数值求积公式的代数精度,但同样必须提供n+1阶导数值,文献[10]的求积公式没有承袭性.本文从Hermite插值多项式出发,构造了具有误差量级为O(h5),且不需要计算导数值,只需要提供求积节点函数值的高精度数值求积公式.1 高精度数值求积公式1.1 Hermite插值多项式的构造构造满足插值条件f(xk)=p(xk),f(xk+1)=p(xk+1),f′(xk)=p′(xk),f′(xk+1)=p′(xk+1)次数不超过3的Hermite插值多项式的p (x).即上式中:αk(x),αk+1(x),βk(x),βk+1(x)为插值基函数,它们分别满足αk(x),αk+1(x),βk(x),βk+1(x)的表达式分别为其中:h=xk+1-xk.构造的插值多项式余项表达式为式(1)中:ξ在xk与xk+1之间,并且与x有关.1.2 数值求积公式的构造对式(1)两边求积分,应用广义积分中值定理,可得到积分中值定理.定理1 设函数f(x)在区间[x,xk+1]有4阶连续导函数,则成立,η在xk与xk+1之间.由式(2)得到带有一阶导数的数值求积分公式,即由式(2)可知,数值求积公式(3)具有3次代数精度.在式(2)中,令xk=a,xk+1=x,则有为了提高数值求积公式的代数精度,分析定理1中间点η的渐近性,则可得到定理2.定理2 设f(x)充分光滑,且存在正整数α及常数,则式(4)的中间点η的渐近性质为证明令应用3次洛必达法可得应用式(4)可得应用式(5),(6)可得记若f(5)(a)≠0,由定理2知α=1,λ=1/2,可修正数值积分公式(3).定理3 设f(x)充分光滑,且f(5)(xk)≠0,则数值求积公式至少有5次代数精度.证明显然,当f(x)=1,x,x2,x3 时,式(7)准确成立.当f(x)=x4 时,左边右边=即有左边=右边,故当f(x)=x4时,式(7)准确成立.当f(x)=x5时,左边右边左边=右边,故当f(x)=x5时,式(7)准确成立.综上所述,求积公式(7)至少有5次代数精度.由式(7),可得到具有5次代数精度的第2个数值求积公式,记为应用复化求积得为了避免求导数,对上式的数值求积公式进行修正,应用公式得到只需计算求积节点函数值,无须提供求积节点导数值的两个数值求积公式,分别为公式(9),(10)的误差量级分别为O(h4),O(h5).2 数值试验例1 计算积分,精确值为1-cos 1=0.459 697 694 131 86….应用梯形公式、求积公式(3),(8)计算,计算结果如表1所示.应用复化梯形公式、求积公式(9),(10)计算,计算结果如表2所示.表1 梯形公式、求积公式(3),(8)的计算结果Tab.1 Numerical experiment of trapeziod formula and formula(3),(8)0.420 735 492 403 0.038 962 20求积公式(3) 0.459 043 633 581 0.000 654 06求积公式(8)方法计算结果绝对误差梯形公式0.459 709 502 385 0.000 018 08表2 复化梯形公式,求积公式(9),(10)的计算结果Tab.2 Numerical experiment of compound trapeziod formula and formula(9),(10)05复化梯形公式1.533 3×10-3 3.381 5×10-4 9.577 4×10方法绝对误差h=0.2 h=0.1 h=0.-5求积公式(9)1.355 1×10-5 1.033 6×10-6 7.027 2×10-8求积公式(10)1.191 2×10-6 6.908 2×10-8 4.091 7×10-9从计算结果看,公式(3)与公式(8)比梯形公式的精度高很多,但必须提供一阶导数或四阶导数;而公式(9)和公式(10),却不用计算导数,计算的节点函数值的个数与复化梯形公式一样.因此,计算量与复化梯形公式相当.复化梯形公式的误差量级为O(h2),而公式(9)的误差量级为O(h4),公式(10)的误差量级为O(h5).因此,公式(9)和公式(10)是非常有效的,无须计算导数的两个数值积分公式.参考文献:[1]李庆杨.数值分析[M].武汉:华中科技大学出版社,1986:83-106. [2]李庆杨.数值计算原理[M].北京:清华大学出版社,2000:251-282. [3]李毅夫.Simpson公式余项中间点渐进性定理及Simpson公式的改进[J].贵州师范大学学报:自然科学版,2007,25(4):67-69.[4]程海来.关于Simpson公式的两点注记[J].数学的实践与认识,2007,27(21):91-93.[5]杜跃鹏,肖泽昌.改进Simpson公式及误差分析[J].高师理科学刊,2008,28(4):27-29.[6]刘彬清.一类高斯求积公式的极限性质[J].工程数学学报,2003,20(4):137-139.[7]吴新元.一个高精度数值积分公式[J].计算物理,1988,5(4):473-477.[8]郑华盛.高精度数值积分公式的构造及其应用[J].数学的实践与认识,2007,37(15):141-148.[9]刘彬清.Newton-cotes数值求积公式的渐近性[J].上海大学学报:自然科学版,2002,8(6):503-506.[10]刘彬清.一类高斯求积公式的极限性质[J].工程数学学报,2003,20(4):137-139.。
数值积分公式的代数精度好的,以下是为您生成的关于“数值积分公式的代数精度”的文章:在数学的广袤世界里,数值积分公式就像是一把神奇的钥匙,能帮助我们打开解决各种复杂问题的大门。
而这其中,代数精度可是一个相当重要的概念。
咱先来说说啥是数值积分公式。
想象一下,你面前有一个形状怪异的曲线,要计算它下面的面积,这可不好直接算吧?这时候数值积分公式就派上用场啦!它能通过一些巧妙的方法,给咱一个大概的面积值。
那代数精度又是啥呢?简单说,就是一个数值积分公式能精确计算出多少种多项式函数的积分值。
比如说,如果一个公式对于一次多项式能精确计算积分,那它的代数精度至少是 1;要是能精确计算二次多项式的积分,代数精度至少就是 2 啦。
我记得有一次给学生们讲这个的时候,有个学生一脸懵地问我:“老师,这东西到底有啥用啊?”我笑着跟他们说:“就好比你要去买水果,你不知道每种水果的价格,但是你可以通过一些估算的方法,大致算出你要花多少钱。
数值积分公式就是帮我们做这种估算的,而代数精度就是衡量这个估算准不准的一个标准。
”咱们来看看一些常见的数值积分公式,比如梯形公式。
它的形式很简单,就是把曲线下的面积近似看成梯形的面积。
那它的代数精度是多少呢?经过一番推导和计算,我们会发现梯形公式的代数精度是1 。
这意味着对于像直线这样的一次多项式函数,它算出来的积分值是比较准的。
再说说辛普森公式。
这个公式就更厉害一点啦,它把曲线下的面积分成了一个抛物线形的部分来近似计算。
它的代数精度能达到 3 呢,对于二次、三次多项式函数都能给出比较精确的积分结果。
在实际应用中,我们要根据具体的问题来选择合适的数值积分公式。
如果问题比较简单,对精度要求不高,梯形公式可能就够用了;要是问题复杂,精度要求高,那就得用上像辛普森公式这样代数精度更高的公式。
就像上次我帮一个工程师朋友计算一个零件的不规则形状的面积。
他一开始用了个简单的方法,结果误差太大,差点导致生产出来的零件不合格。
重庆科技学院毕业设计(论文)题目高精度数值积分公式的构造及其应用学院专业班级指导教师职称讲师评阅教师职称年月日注意事项1.设计(论文)的内容包括:1)封面(按教务处制定的标准封面格式制作)2)原创性声明3)中文摘要(300字左右)、关键词4)外文摘要、关键词5)目次页(附件不统一编入)6)论文主体部分:引言(或绪论)、正文、结论7)参考文献8)致谢9)附录(对论文支持必要时)2.论文字数要求:理工类设计(论文)正文字数不少于1万字(不包括图纸、程序清单等),文科类论文正文字数不少于1.2万字。
3.附件包括:任务书、开题报告、外文译文、译文原文(复印件)。
4.文字、图表要求:1)文字通顺,语言流畅,书写字迹工整,打印字体及大小符合要求,无错别字,不准请他人代写2)工程设计类题目的图纸,要求部分用尺规绘制,部分用计算机绘制,所有图纸应符合国家技术标准规范。
图表整洁,布局合理,文字注释必须使用工程字书写,不准用徒手画3)毕业论文须用A4单面打印,论文50页以上的双面打印4)图表应绘制于无格子的页面上5)软件工程类课题应有程序清单,并提供电子文档5.装订顺序1)设计(论文)2)附件:按照任务书、开题报告、外文译文、译文原文(复印件)次序装订3)其它学生毕业设计(论文)原创性声明本人以信誉声明:所呈交的毕业设计(论文)是在导师的指导下进行的设计(研究)工作及取得的成果,设计(论文)中引用他(她)人的文献、数据、图件、资料均已明确标注出,论文中的结论和结果为本人独立完成,不包含他人成果及为获得重庆科技学院或其它教育机构的学位或证书而使用其材料。
与我一同工作的同志对本设计(研究)所做的任何贡献均已在论文中作了明确的说明并表示了谢意。
毕业设计(论文)作者(签字):年月日摘要求解函数()f x 在[,]a b 区间上的定积分()ba f x dx ⎰时,如果被积函数()f x 在[,]ab 区间上原函数很难用初等函数表达,我们就不能够借助牛顿-莱布尼兹公式来计算此定积分()ba f x dx ⎰。
另外,许多实际问题中的被积函数()f x 往往是列表函数或其他形式的非连续函数,显然其原函数没有意义,所以对这类函数的积分,也不能用经典的不定积分方法求解。
因此,数值积分的理论与方法一直是计算数学研究的基本课题。
本文首先总结了数值积分的基本思想和几类常用的数值积分方法,并且给出了数值积分稳定的一般性条件。
接着,我们基于文献[10]提出了一个改进的三点高斯公式,通过理论分析,此公式具有7次代数精度。
最后,利用一个数值算例验证了我们提出公式相比文献[10]的两点高斯公式无论是在代数精度还是在数值精度方面都有较大提高。
并将该公式应用到2010年数学建模A 题,取得了比较好的结果。
关键词: 数值积分方法 三点高斯公式 代数精度ABSTRACTIf the primitive function of integrand f(x) cannot be expressed by the elementary function in definite integral computation. We cannot calculate the definite integral()baf x dxby using Newton-Leibniz formula. In real world, integrand is often list function or other forms of discontinuous function in many practical problems. For this kind of function of definite integral ,its primary function is obviously unmeaningful; it cannot use the indefinite integral method to solve. Therefore, the theory and method of numerical integration is always the computational mathematics basic topic.This paper first summarizes the basic ideas of numerical integration and some common numerical integration formates, in addition, represents the general stability condition of numerical integration. Next, we propose an improvement two-point gaussian formula based on the literature [10]. It has seven-time algebraic precision essentially. Finally, the experimental results are represented, which indicate that our numerical format is superior in algebraic precision and numerical precision compared to the three-point Gaussian formula in [10] and some of the other classical numerical format. and the formula is applied to the 2010 mathematical modeling problem A, achieved good results.Keywords : Numerical integral method ;Three-point gauss formula ;Algebra precision目录摘要 (I)ABSTRACT (II)1绪论 (1)2数值积分的计算方法 (2)2.1 数值积分的基本思想与评价指标 (2)2.2 几种常用数值积分方法 (5)2.2.1 插值型求积公式 (5)2.2.2 Newton-Cotes公式 (6)2.2.3 复合求积公式 (8)2.2.4 逐次分半技术与Romberg公式 (9)2.2.5 Gauss型求积公式 (10)2.2.6 Gauss-Legendre求积公式 (13)2.2.7 Gauss-Chebyshev求积公式 (14)3改进三点Gauss公式 (16)3.1改进三点Gauss公式的介绍 (16)3.2 数值算例 (17)3.3 2010年数学建模A题求解 (18)4 结束语 (28)参考文献 (29)致谢 (30)附录 (31)1绪论数值积分是求定积分的近似值的数值方法。
即用被积函数的有限个抽样值的加权平均近似值代替定积分的值。
数值积分是计算方法或数值分析课程中非常重要的教学内容,数值积分方法也是解决实际计算问题的重要方法。
我们知道计算定积分是采用牛顿—莱布尼兹(Newton-Leibniz)公式:()()()baf x dx F b F a =-⎰但由于其适用范围有限,不能普遍适用,因此有其局限性。
所有说牛顿—莱布尼兹公式不是万能的,而数值积分公式却具备这种良好的性质。
只要数值积分构造得当,就能很好的计算出某个定积分的近似值。
以下罗列出牛顿—莱布尼兹公式不适用的三种情况:(1) ()x f 的解析式没有给出,只给出了()x f 的一些离散点。
(2) ()x f 的原函数不能用初等函数表示,如:()dx e dx x x x ⎰⎰-10102sin ,。
(3)原函数表达式相当复杂,计算十分不便。
由于以上原因,数值积分的理论与方法一直是计算数学研究的基本课题。
对微积分学作出杰出贡献的数学大师,如牛顿、欧拉、高斯等人也在数值积分这个领域作出了各自的贡献,并奠定了它的理论基础。
数值积分还是微分方程数值解法的重要依据。
许多重要公式都可以用数值积分方程导出。
因此,数值积分的理论与方法还是其他学科的理论依据。
由于数值积分是求解定积分近似值的数值方法,所以它的意义在于能够求出定积分的近似值。
而定积分又具有广泛应用,它几乎是所有课程的公共基础课。
以下列举了它的一部分应用:1、 计算图形面积、曲线弧长、立体图像体积。
2、 在力学中的应用,计算力做的功、位移、能量等作用。
3、 在电学中的应用,计算场强、电势差、电压等作用。
4、利用定积分证明不等式。
由于定积分的应用广泛,作用巨大,而高精度数值积分是计算定积分近似值的良好数值方法。
由此可以看出构造高精度数值积分公式是十分必要的。
由于高精度数值积分具有计算结果准确、代数精度高、使用方便、稳定性好等优点。
因此,探讨高精度数值积分的构造及其应用具有明显的实际意义。
2数值积分的计算方法2.1 数值积分的基本思想与评价指标[1]对于数值积分法的思想来源于定积分的定义,即1()lim ()nbk k ak f x dx f x x λ→==∆∑⎰其中m a x{}k kx λ=∆,一般的提法是:用()f x 在点01n a x x x b =<<⋅⋅⋅<=处的函数值()(0,1,,)k f x k n =的线性组合 00110()()()()()nn n n k k k I f A f x A f x A f x A f x ==++⋅⋅⋅⋅⋅⋅+=∑作为积分()()baI f f x dx =⎰的近似值,即()()()nbk k ak I f f x dx A f x ==≈∑⎰ (2.1)并称此为数值求积公式,也称为机械求积公式。
形如:()()()n R f I f I f =-为求积公式(2.1)的余项或误差,k x 及(0,1,2,,)k A k n =⋅⋅⋅分别称为求积公式(2.1)的求积节点及求积系数,这里求积系数(0,1,2,,)k A k n =⋅⋅⋅只与积分区间[,]a b 有关,而与()f x 无关。
为保证机械求积公式的精度, 自然希望它对尽可能多的简单函数是准确成立的,如果要求它对一切不超过n 次多项式都准确成立, 而对1n +次多项式不一定准确成立。
则得到关于系数k A 的1n +阶线性方程组:0122001111001121n n nn n nn nn nA A A b a b a A x A x A x b a A x A x A xn +++++=-⎧⎪-⎪+++=⎪⎨⎪⎪-+++=⎪+⎩由于系数行列式为范德蒙德行列式,其值不为零,则解k A 是唯一确定的。