第14卷第1期2011年1月
STUDIES1NC()LLEGE
高等数学研究
MATHEMATICS
V01.14,No.1
Jan.,2011
关于定积分近似计算中矩形法的误差估计
郑立飞,解小莉,王洁
(西北农林科技大学应用数学系,陕西杨凌712loo)
摘要利用拉格朗日中值定理给出复合矩形法的误差估计,并指出复合矩形法只有在等分的次数很大的时
候才能比较精确的估计所求定积分的值.最后.给出复合矩形法、复合梯形法及复合抛物线法的误差估计实例.
关键词
复合矩形法;分部积分法;拉格朗日中值定理;误差估计
013;()24
中图分类号文献标识码
A
文章编号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],(复合)梯形法的误差为
6
“形盛=掣,
高等数学研究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
7
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
L
Finny。MauriceDweir,FrankRGiordano.
M;攀≈2.59808.
Thomas’calculus:tenth
edition[M].影印.北京:高等
Z
教育出版社,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
中图分类号文献标识码
A
文章编号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],(复合)梯形法的误差为
6
“形盛=掣,
高等数学研究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
7
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
L
Finny。MauriceDweir,FrankRGiordano.
M;攀≈2.59808.
Thomas’calculus:tenth
edition[M].影印.北京:高等
Z
教育出版社,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