关于定积分近似计算中矩形法的误差估计

第14卷第1期2011年1月

STUDIES1NC()LLEGE

高等数学研究

MATHEMATICS

V01.14,No.1

Jan.,2011

关于定积分近似计算中矩形法的误差估计

郑立飞,解小莉,王洁

(西北农林科技大学应用数学系,陕西杨凌712loo)

摘要利用拉格朗日中值定理给出复合矩形法的误差估计,并指出复合矩形法只有在等分的次数很大的时

候才能比较精确的估计所求定积分的值.最后.给出复合矩形法、复合梯形法及复合抛物线法的误差估计实例.

关键词

复合矩形法;分部积分法;拉格朗日中值定理;误差估计

013;()24

中图分类号文献标识码

文章编号1008—1399(2011)0l一0005—02

误差估计是数值计算中对算法好坏衡量的一个基本标准.文[1]虽给出了定积分的几个近似计算公式,包括矩形法(复合左矩形法和复合右矩形法)、梯形法和抛物线法,但并没有给出这几个算法的误差估计.

梯形法和抛物线法的误差估计可参考文[2—4],而复合左矩形法和复合右矩形法的近似计算的误差估计却在多数情形下仅以练习题的形式出现.本文在这里借用高等数学的方法给出复合左矩形法近似计算的误差估计,以供教师和学生参考.复合右矩形法的近似计算估计可以同理给出.

为了叙述的方便,首先给出定积分左矩形法计算定积分值的算法公式,该公式来源于文献[1].

设厂(z)在■,6]上连续,在区间■,6]中插入分点

口一zo<z1<…<z,1<z。=6,

将其划分成竹等分,并记

厂(z卜1)=y-l,

凼此,

If厂(z)dz一厂(以)(6一口)l=m化)叫口)]dzl-眇㈣(z刊dzl≤MIJ.:cz一口,d辱l=等c6一n,2.

定理2

设,(工)在[口,6]上连续,在[口,6]上的

至多有限个点之外有连续导数,且存在正常数M使

I厂(z)I≤M.

对区间[口,6]进行咒等分划分:

口一zo<zl<…<z,1<z。一6,

并记

弛一一舶,^=z。一珀=字,

为b么,

则复合左矩形公式为

r6厂(z)d工≈蔓亍旦(y。+y,+…+j,r,).

定理1

lf厂(z)出一,l宴弦。I≤篆(6一口)2.

J_

f=l

厶,‘

证明利用定理l可得。

设厂(z)在[口,阳上连续,在[口,6]上的

至多有限个点之外有连续导数,且存在正常数M使

忆m)出一Ⅳ(珀)f≤丢胁;.

由此即知

I厂(z)I≤M,

那么,

lIl,(z)出一厂((£)(6一n)l≤÷M(6一n)2.f6,(z)出一厂((£)(6一n)l≤丢M(6一n)z.

证明

据拉格朗日中值定理,任给z∈[n,6],

№汕一^喜弛一I≤

客l罡。化,出叫c珀,I≤

号胁2

存在e∈(n,6),使得

,(z)一厂(口)=/(拿)(z一口).

收稿日期:2009—12—07;修改日期:2010—08—3】.

基金项目:西北农林科技大学教学改革项目资助(JY0902104).作者简介。郑立飞(1973一),男,陕西临潼人,硕士,讲师,主要从事生

物数学研究.Email:zhenglifei@126.com.

解小莉(1973一),男,陕西临潼人,硕士,讲师,主要从事生物数学研究.Email;小englifei@126.com.

2扣(6~)2.

对于复合右矩形法,可以同样证明,其误差估计与复合左矩形法相同.

从定理2可知,矩形法在使用的时候,要求有足够多的分点,即咒要足够得大,否则精度难以保证.

根据文[5],(复合)梯形法的误差为

“形盛=掣,

高等数学研究2011年1月

蕊形法2赤×2・59808×1=o・129904,

而(复合)抛物线法的误差为

如物线法2志M(6一盘)2・

这里的M定义同定理2.由此可知,复合梯形法要优“铂缱击:垡掣:o.0360844.妇形击2—丽一2

确形击:掣:咖64592,

u・∞43了z’

魄铂缱击2—1矿2

u・u弘u苍44・

于复合矩形法,而复合抛物线法优于复合梯形法.因这个误差上界有些大.此在实际中,计算定积分近似值时经常使用复合抛下面利用文[4]中给出的复合梯形法和复合抛

物线法.但是复合矩形法作为定积分近似计算的一物线法的误差公式

种朴素思想,对理解其它方法具有启发作用,因此重

£梯形岳==——■i—主一n口o’

要性是不言而喻的.

£…=%笋Mo,

文[5]给出的误差公式只需用到一阶导数,计算e-…=%菩Mo,

量较少,也方便对三种算法的误差进行比较.而文[4]其中Mo为l/,(z)l的上界,计算上述例子中复合

给出的误差公式需要用到二阶导数,计算量较大,不梯形法和复合抛物线法的误差,分别可得

方便对三种算法的误差进行比较.但是,文[4]给出的误差上限小于文[5]给出的误差上限,因此取文[4]e梯形法:彳黑:o.67×10一2,6梯形法2西丽2

u.b

x1u‘’

的方法来计算误差上限更好,而对复合矩形法误差的

计算可以利用本文给出的误差公式.

.£抛铀线法=志×10叫≈o.4×10一.

以下通过实例来说明如何比较复合矩形法、复这里,

合梯形法与复合抛物线法的误差.

.M0=8,疗一10.

例1

试给出用复合矩形法、复合梯形法和复

仍然可以看出:复合矩形法计算的误差较大,而复合合抛物线法计算定积分fr与dz的误差.

抛物线法的计算误差最小.

参考文献

不妨令

・,(z)2南,

[1]同济大学数学系.高等数学:上册[M].北京:高等教育

出版社,2007:229—231.

[2]华中理工大学数学系.计算方法[M].北京:高等教育出

那么,

I,,(训=l志l,

版社,1999:93-94.

[3]关治,陆金甫.数值分析基础[M].北京:高等教育出版

社,1998;183・188.

可以求得上式的最大值为

[4]Ross

Finny。MauriceDweir,FrankRGiordano.

M;攀≈2.59808.

Thomas’calculus:tenth

edition[M].影印.北京:高等

教育出版社,2004:376—378.

将积分区间均等分为10份,利用文[5]给出的[5]刘征,丁桂艳.关于定积分近似计算的误差估计[J].鞍

误差公式,上述三种算法进行计算后的误差分别为:

山科技大学学报.2003.26(4):313.317.

ErrorEstimationofApproximatingDefiniteIntegrals

byReCtan2uIarRule

ZHENGLi—fei,

XlEXia扩li,

WANGJie

(AppliedMathematicaIDepartment.Northwest

A&FUniversity.Yangling712lOO.PRC)

Abstract:

ThispapermainlyusedLagrangeMean

Value

Theorem

to

formulate

anerror

estimationofapproximatingdefiniteintegralsbyrectangularrule.ItindicatesthattheCompositeRectangularRuleyieldsbetterresultsifandonlyiftheintervalispartitionedintosufficientmanysman

intervals.

Anexample

is

used

to

compare

the

errors

generated,respectively,

by

the

compositerectangular,trapezoidalandparabo王icrules.

Keywords:

Composite

Rectangular

Rule,integrationbyparts,Lagrange

MeanValue

Theorem。errorestimation

关于定积分近似计算中矩形法的误差估计

作者:作者单位:刊名:英文刊名:年,卷(期):

郑立飞, 解小莉, 王洁, ZHENG Li-fei, XIE Xiao-li, WANG Jie西北农林科技大学应用数学系,陕西,杨凌,712100高等数学研究

STUDIES IN COLLEGE MATHEMATICS2011,14(1)

参考文献(5条)

1. 刘征;丁桂艳 关于定积分近似计算的误差估计 2003(04)

2. Ross L Finny;Maurice D Weir;Frank R Giordano Thomas' Calculus 20043. 关治;陆金甫 数值分析基础 19984. 华中理工大学数学系 计算方法 19995. 同济大学数学系 高等数学 2007

本文链接:http://d.g.wanfangdata.com.cn/Periodical_gdsxyj201101003.aspx

第14卷第1期2011年1月

STUDIES1NC()LLEGE

高等数学研究

MATHEMATICS

V01.14,No.1

Jan.,2011

关于定积分近似计算中矩形法的误差估计

郑立飞,解小莉,王洁

(西北农林科技大学应用数学系,陕西杨凌712loo)

摘要利用拉格朗日中值定理给出复合矩形法的误差估计,并指出复合矩形法只有在等分的次数很大的时

候才能比较精确的估计所求定积分的值.最后.给出复合矩形法、复合梯形法及复合抛物线法的误差估计实例.

关键词

复合矩形法;分部积分法;拉格朗日中值定理;误差估计

013;()24

中图分类号文献标识码

文章编号1008—1399(2011)0l一0005—02

误差估计是数值计算中对算法好坏衡量的一个基本标准.文[1]虽给出了定积分的几个近似计算公式,包括矩形法(复合左矩形法和复合右矩形法)、梯形法和抛物线法,但并没有给出这几个算法的误差估计.

梯形法和抛物线法的误差估计可参考文[2—4],而复合左矩形法和复合右矩形法的近似计算的误差估计却在多数情形下仅以练习题的形式出现.本文在这里借用高等数学的方法给出复合左矩形法近似计算的误差估计,以供教师和学生参考.复合右矩形法的近似计算估计可以同理给出.

为了叙述的方便,首先给出定积分左矩形法计算定积分值的算法公式,该公式来源于文献[1].

设厂(z)在■,6]上连续,在区间■,6]中插入分点

口一zo<z1<…<z,1<z。=6,

将其划分成竹等分,并记

厂(z卜1)=y-l,

凼此,

If厂(z)dz一厂(以)(6一口)l=m化)叫口)]dzl-眇㈣(z刊dzl≤MIJ.:cz一口,d辱l=等c6一n,2.

定理2

设,(工)在[口,6]上连续,在[口,6]上的

至多有限个点之外有连续导数,且存在正常数M使

I厂(z)I≤M.

对区间[口,6]进行咒等分划分:

口一zo<zl<…<z,1<z。一6,

并记

弛一一舶,^=z。一珀=字,

为b么,

则复合左矩形公式为

r6厂(z)d工≈蔓亍旦(y。+y,+…+j,r,).

定理1

lf厂(z)出一,l宴弦。I≤篆(6一口)2.

J_

f=l

厶,‘

证明利用定理l可得。

设厂(z)在[口,阳上连续,在[口,6]上的

至多有限个点之外有连续导数,且存在正常数M使

忆m)出一Ⅳ(珀)f≤丢胁;.

由此即知

I厂(z)I≤M,

那么,

lIl,(z)出一厂((£)(6一n)l≤÷M(6一n)2.f6,(z)出一厂((£)(6一n)l≤丢M(6一n)z.

证明

据拉格朗日中值定理,任给z∈[n,6],

№汕一^喜弛一I≤

客l罡。化,出叫c珀,I≤

号胁2

存在e∈(n,6),使得

,(z)一厂(口)=/(拿)(z一口).

收稿日期:2009—12—07;修改日期:2010—08—3】.

基金项目:西北农林科技大学教学改革项目资助(JY0902104).作者简介。郑立飞(1973一),男,陕西临潼人,硕士,讲师,主要从事生

物数学研究.Email:zhenglifei@126.com.

解小莉(1973一),男,陕西临潼人,硕士,讲师,主要从事生物数学研究.Email;小englifei@126.com.

2扣(6~)2.

对于复合右矩形法,可以同样证明,其误差估计与复合左矩形法相同.

从定理2可知,矩形法在使用的时候,要求有足够多的分点,即咒要足够得大,否则精度难以保证.

根据文[5],(复合)梯形法的误差为

“形盛=掣,

高等数学研究2011年1月

蕊形法2赤×2・59808×1=o・129904,

而(复合)抛物线法的误差为

如物线法2志M(6一盘)2・

这里的M定义同定理2.由此可知,复合梯形法要优“铂缱击:垡掣:o.0360844.妇形击2—丽一2

确形击:掣:咖64592,

u・∞43了z’

魄铂缱击2—1矿2

u・u弘u苍44・

于复合矩形法,而复合抛物线法优于复合梯形法.因这个误差上界有些大.此在实际中,计算定积分近似值时经常使用复合抛下面利用文[4]中给出的复合梯形法和复合抛

物线法.但是复合矩形法作为定积分近似计算的一物线法的误差公式

种朴素思想,对理解其它方法具有启发作用,因此重

£梯形岳==——■i—主一n口o’

要性是不言而喻的.

£…=%笋Mo,

文[5]给出的误差公式只需用到一阶导数,计算e-…=%菩Mo,

量较少,也方便对三种算法的误差进行比较.而文[4]其中Mo为l/,(z)l的上界,计算上述例子中复合

给出的误差公式需要用到二阶导数,计算量较大,不梯形法和复合抛物线法的误差,分别可得

方便对三种算法的误差进行比较.但是,文[4]给出的误差上限小于文[5]给出的误差上限,因此取文[4]e梯形法:彳黑:o.67×10一2,6梯形法2西丽2

u.b

x1u‘’

的方法来计算误差上限更好,而对复合矩形法误差的

计算可以利用本文给出的误差公式.

.£抛铀线法=志×10叫≈o.4×10一.

以下通过实例来说明如何比较复合矩形法、复这里,

合梯形法与复合抛物线法的误差.

.M0=8,疗一10.

例1

试给出用复合矩形法、复合梯形法和复

仍然可以看出:复合矩形法计算的误差较大,而复合合抛物线法计算定积分fr与dz的误差.

抛物线法的计算误差最小.

参考文献

不妨令

・,(z)2南,

[1]同济大学数学系.高等数学:上册[M].北京:高等教育

出版社,2007:229—231.

[2]华中理工大学数学系.计算方法[M].北京:高等教育出

那么,

I,,(训=l志l,

版社,1999:93-94.

[3]关治,陆金甫.数值分析基础[M].北京:高等教育出版

社,1998;183・188.

可以求得上式的最大值为

[4]Ross

Finny。MauriceDweir,FrankRGiordano.

M;攀≈2.59808.

Thomas’calculus:tenth

edition[M].影印.北京:高等

教育出版社,2004:376—378.

将积分区间均等分为10份,利用文[5]给出的[5]刘征,丁桂艳.关于定积分近似计算的误差估计[J].鞍

误差公式,上述三种算法进行计算后的误差分别为:

山科技大学学报.2003.26(4):313.317.

ErrorEstimationofApproximatingDefiniteIntegrals

byReCtan2uIarRule

ZHENGLi—fei,

XlEXia扩li,

WANGJie

(AppliedMathematicaIDepartment.Northwest

A&FUniversity.Yangling712lOO.PRC)

Abstract:

ThispapermainlyusedLagrangeMean

Value

Theorem

to

formulate

anerror

estimationofapproximatingdefiniteintegralsbyrectangularrule.ItindicatesthattheCompositeRectangularRuleyieldsbetterresultsifandonlyiftheintervalispartitionedintosufficientmanysman

intervals.

Anexample

is

used

to

compare

the

errors

generated,respectively,

by

the

compositerectangular,trapezoidalandparabo王icrules.

Keywords:

Composite

Rectangular

Rule,integrationbyparts,Lagrange

MeanValue

Theorem。errorestimation

关于定积分近似计算中矩形法的误差估计

作者:作者单位:刊名:英文刊名:年,卷(期):

郑立飞, 解小莉, 王洁, ZHENG Li-fei, XIE Xiao-li, WANG Jie西北农林科技大学应用数学系,陕西,杨凌,712100高等数学研究

STUDIES IN COLLEGE MATHEMATICS2011,14(1)

参考文献(5条)

1. 刘征;丁桂艳 关于定积分近似计算的误差估计 2003(04)

2. Ross L Finny;Maurice D Weir;Frank R Giordano Thomas' Calculus 20043. 关治;陆金甫 数值分析基础 19984. 华中理工大学数学系 计算方法 19995. 同济大学数学系 高等数学 2007

本文链接:http://d.g.wanfangdata.com.cn/Periodical_gdsxyj201101003.aspx


相关内容

  • 一阶微分方程解的存在定理
  • 第三章 一阶微分方程解的存在定理 [教学目标] 1. 理解解的存在唯一性定理的条件.结论及证明思路,掌握逐次逼近法,熟练近似解的误差估计式. 2. 了解解的延拓定理及延拓条件. 3. 理解解对初值的连续性.可微性定理的条件和结论. [教学重难点] 解的存在唯一性定理的证明,解对初值的连续性.可微性定 ...

  • 计算曲边梯形面积
  • 78 上海中学数学2014年第1-2期 计算曲边梯形面积 362261 福建省晋江市养正中学刘华湘 计算曲边梯形面积的过程蕴涵着定积分的基本思想方法,笔者从现实背景导人问题,从几何直观到数列求和,从代数推理以及算法编程求和两方面求 快捷便利的瞬间,即可体现"无限". (一) 教学 ...

  • 蒙特·卡罗方法
  • 蒙特·卡罗方法 百科名片 蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法.是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法.与它对应的是确 ...

  • 高斯求积型公式及其程序开发
  • 2011年1月19日 摘要 在实际中,求非代数函数的积分往往要求精度很高,因为高斯型求积公式据有最高代数精度且高斯型求积公式是收敛和稳定的,同时它可以使更多的函数准确成立,所以研究高斯型求积公式及其程序开发是很必要的. 目的是总结分析高斯型求积公式,在掌握其基本思想的基础上,深入学习几种常见的高斯型 ...

  • 数值分析中牛顿迭代法的引入方法探讨
  • 第25卷第5期2010年lO月 天中学刊 JournalofTianzhong .,01.25No.5 oct.20lO 数值分析中牛顿迭代法的引入方法探讨 王霞,张启虎 (郑州轻工业学院数学与信息科学系,河南郑州450002) 摘要:数值分析中牛顿迭代法是求解非线性方程的基本方法.与一般教材上牛顿 ...

  • 浅谈几个著名的大数定律及应用
  • 浅谈几个著名的大数定律及应用 李 蕊 (青海大学成教学院,青海 西宁,810001) 摘 要 大数定律以严格的数学形式表达了随机现象最根本的性质--平均结果的稳定性,是随机现象统计规律性的具体表现,本文介绍了几种常用的大数定律,并给出一些简单应用. 关键词 大数定律 随机变量 数学期望 概率 1 引 ...

  • 有限差分法
  • 有限差分法 有限差分法 finite difference method 微分方程和积分微分方程数值解的方法.基本思想是把连续的定解区域用有限个离散点构成的网格来代替, 这些离散点称作网格的节点:把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似:把原方程和定解条件中的微商用差商来近 ...

  • 电磁场数值分析基础
  • 教材与参考书 电磁场数值分析基础 邵 维 办公室:科研楼704 [email protected]  教材: <电磁场数值分析基础>讲义  参考书: (1)<数值分析>,钟尔杰等,高等教育出版社, 2004. (2) <数值分析> (第四版),李庆扬等, ...

  • 基于容积卡尔曼滤波的卫星姿态估计
  • 第34卷第2期2013年2月 宇航学报 JournalofAstronautics V01.34February No.22013 基于容积卡尔曼滤波的卫星姿态估计 魏喜庆,宋申民 (哈尔滨工业大学控制理论与制导技术研究中心,哈尔滨150001) 摘要:为了获得更好的估计精度和滤波稳定性,提出了一种 ...