定积分问题的数值求解及Matlab实现.pdfVIP

第28卷第5期 哈 尔滨 商 业 大 学 学报 (自然科学版) Vo1．28No．5 2012年 10月 JournalofHarbinUniversityofCommerce(NaturalSciencesEdition) Oct．2012 定积分问题的数值求解及 Matlab实现 张长耀 ，刘秀丽 (西藏农牧学院公共教学部 ，西藏 林芝860000) 摘 要：以一道定积分题 目为例，分别利用矩形法、梯形法和抛物线法等数值方法给出该定积分的近 似求解方法，借助Matlab编程，得到该定积分问题的近似解，并且给出该类型定积分的数值求解方 法 ． 关键词：定积分；Matlab；误差分析 中图分类号：0241．4 文献标识码：A 文章编号：1672—0946(2012)05—0620—03 Studyonnumericalsolutionon definiteintegralsand itsrealizationbyM atlab ZHANG Chang—yao，LIU Xiu—li (DepartmentofPublicTeaching，TibetAgricultureandAnimalHusbandryCollege，Linzhi860000，China) Abstract：Inthispaper，thenumericalsolutionondefiniteintegralswasintroducedbasedon anexampleusingthreemethodsrespectivelysuchasrectangularformula，trapezium formula andparabolaformula．AndthecorrespondingMatlabprogramsweregiven． Keywords：definiteintegral；Matlab；erroranalysis 实际问题当中常常需要计算积分．依据人们 求解定积分的常用数值方法有矩形法、梯形法 所 熟知 的微积分基本 定理，对 于积分 ，= 和抛物线法，下面我们借助Matlab编程分别实现． ) ，只要找到被积函数 )的原函数，()， 1 中点矩形法 便有下列牛顿 一菜布尼茨公式： rb 根据定积分的定义，对于定积分I_厂()dx，将 I )dx=F(b)一F(0)． Ja 积分区间[a，b]分成n个小区间，当积分区间被分 但实际使用这种求积方法往往有困难，因为大 割得很细时，即 很大时，积分和可以看作是定积 量的被积函数，诸如 ( ≠0)，e 等，其原函数 分的一个近似值，即 不能用初等函数表达 ，故不能用上述公式计算．下 Jnb厂() n ) 面我们用教学中遇到的一个定积分的题 目作为例 在Matlab程序中，我们将[a，b]n等分，对等分 子，来说明此类定积分题 目的求解． 区间 婆 例n求定积分I,／1+sin。 ． X0 = a l … Xi JO a++— — 0 …… = b 虽然这个题 目的被积函数是形式较简单的初 n 等函数，但它的原函数不是初等函数，无法利用牛