定积分的近似计算及MATLAB实现.pdfVIP

第29卷第4期 淮北煤炭师范学院学报(自然科学版) Vo1．29 No．4 2008年 12月 JournalofHuaibeiCoalIndust~TeachersCollege(NaturalScience) Dec．2008 定积分的近似计算及MATLAB实现 吴庆丰 (淮北煤炭师范学院数学科学学院，安徽 淮北 235000) 摘 要：文章考虑了定积分近似计算的矩形法、梯形法、抛物线法和一种在抛物线法基础上的定积分近似计算公式 并给出了每种方法的MATLAB程序． 关键词：定积分；近似计算；误差估计；MATLAB 中图分类号：O242．2 文献标识码：C 文章编号：1672—7177(2008)04—0086—03 定积分计算的基本公式是Newton—Leibniz公式，但对于较复杂的函数要求其精确值是比较困难的，而 且很多形式较简单的初等函数的原函数也不是初等函数，如』 ，fe ，fdx等等，其原函数都不是 初等函数，无法利用Newton—Leibniz公式计算．因此，研究定积分的近似计算是非常必要的，而且，可以通过 计算机编程来方便计算定积分的近似值． 1 矩形法 定理 11【 设函数 ／()在闭区间fa，6]上连续，则 ()在 [a，hi~可积，用一组分点 a=XoXl…‰= b把[a，b】n等分，则 』圳一 宝 一 )(0 )． (1) 例1用矩形法求l dx的近似值． 解按公式(1)有』)d一旦宝厂(一鱼)(0 1)，这里取am10000，￡=0．5． 下面给出矩形法 MATLAB求解程序： ~rmatlong sum：0；a=0；b=1；n=10000；t=0．5；d=(b—a)／n； f0ri=1：l：n sum=sum+d：l=sin(t：l=(2=I：a+(2$i一1)：I：d))．／(t：}：(2 a+(2$i一1) d)) end 下面的定理给出了矩形法的误差估计． 定理2t 设函数 ．厂()在闭区间fa，b】上连续，在开区间(a，b)内至多有限个点外处处有连续的导数， 且存在正数 使得1． 厂()IsM，用一组分点 n= o … =b把[a，b]n等分，则 一 耋( )1 其中常数÷是最佳的． 收稿 日期：2008—02—20 作者简介：吴庆丰(1979 一 )，男，湖北麻城人，硕：匕讲师，研究方向：优化理论及算法 第4期 吴庆丰：定积分的近似计算及MATLAB实现 87 2 梯形法 定理3 设函数 )在闭区间[a，b】上连续，在开区间(a，b)内至多有限个点外处处有连续的导数， 且存在正数 使得l-厂()l 用一组分点 a= o … =b把[a，b】n等分，则 )dx～ ㈦+fb(耋)]l 其中常数÷是最佳的． 梯形法公式：』() [)+6)+2塞 ] 现就上一题给出梯形法MATLAB求解程序： formatlong SUlfl=0；a=0；b=1；n=10000；d=(b—a)／n； fori=1： n—l sum=sum+d术sin(a+i水d)．／(a+i宰d) end sum=sum+0。5$d：I=(1+sin(b)．／b) 3 抛物线法 定理4怕 设 [a，6]一R是四次连续可微的函数，且 M= max I ()I；用一组分点 a=戈。 … =b把 [a，6】rt等分，则 一百b-a +2 4砉厂( )]I 其中常数 是最佳的． 抛物线法公式：』()出 [口)+6)+2 (+4耋(堕 )】 现就第一题给出抛物线法 MATLAB求解程序： formatlong SUm=0；a=0；b=1；n=10000；t=0．5；d=(b—a)／n；