网站大量收购闲置独家精品文档,联系QQ:2885784924

计算方法-数值积分.ppt

  1. 1、本文档共21页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多

数值积分小组成员:引言在数学分析中,当f(x)在区间[a,b]上连续且具有原函数F(x)时,我们往往采用Newton-Leibniz公式来求积分:然而,随着学习的深化,发现牛顿-莱布尼兹公式存在的很大的局限性Newton-Leibniz公式的局限性对大多数f(x)而言,找原函数困难,即使存在原函数也不能用初等函数表示01原函数表达式过于复杂02被积函数由表格给出,没有解析形式,也无法使用Newton-Leibniz公式来求积分03为了避免上述积分过程中存在的问题,我们可以采用数值积分的方法来求解,这样就避免了原函数的求解过程,同时对于由测量或计算得到的数据表表示的f(x)也可以求解主要有五种方法,对应五种计算公式:梯形法中值法辛普森积分法高斯积分法高斯积分法-三点公式数值积分梯形公式x0x0+原理介绍:用一个梯形来近似代替每个子区间的面积,如下图所示,用图中阴影部分梯形面积替代曲边梯形的面积,积分过程如下图所示:梯形公式而对于整个区间[x0,x1],可通过累加求和得到,其中区间[x0,x1]被分为n等份,每个区间长度为,因此区间[x0,x1]上的积分可通过下述式子得到梯形公式算法for(inti=0;i20;i++){ area1=0; x=PI/intervals[i]; for(intj=0;jintervals[i];j++) { area1=area1+0.5*(f(j*x)+f(j*x+x))*x; }doublee=(area1-2)/area2; area2=area1-2; coutsetiosflags(ios::left)setw(6)intervals[i]setw(7)intervals[i]+1setw(20)area1-2setw(6)eendl;}辛普森积分法原理介绍:把区间[x0,x1]分为2n等分,n个区间,在长度为的区间上进行泰勒展开,可得区间[x0,x0+]上的积分形式如下所示:x0x0+x0+辛普森积分法辛普森求积公式的几何意义是用一条过三点的抛物线(如上图中三点)近似代替被积函数的曲线,从而用一个二次抛物线所围成的容易计算的曲边梯形面积(图中阴影部分)来近似代替原来的曲边梯形的面积.x0x0+x0+通过对n个区间按上述公式累加,可得区间[x0,x1]上的积分形式为注意:因为该公式是把区间[x0,x1]划分为2n等份,n个区间,因此在用循环求积分时,要注意其上下限0102辛普森积分法辛普森算法代码for(inti=0;i20;i++){area1=0;x=PI/intervals[i];for(intj=0;jintervals[i];j++){ area1=area1+x*(f(j*x)+4*f(j*x+0.5*x)+f(j*x+x))/6;}doublee=(area1-2)/area2;area2=area1-2;coutsetiosflags(ios::left)setprecision(12)setw(6)intervals[i]setw(7)intervals[i]*2+1setw(20)area1-2setw(6)eendl;}高斯积分通过待定系数法及泰勒展开找到两个相对精确的评估点高斯积分for(i=0;i10;i++){s=0;x=(x1-x2)/intervals[i];for(j=0;jintervals[i];j++)s+=function(x1+(0.5-sqrt(3.0)/6)*x+j*x)+function(x1+(0.5+sqrt(3.0)/6)*x+j*x);s=s*x/2;error1=s-2;Error_Ratio=error1/error2;error2=error1;coutsetiosflags(ios::left)setw(6)intervals[i]setw(7)intervals[i]*2

文档评论(0)

SYWL2019 + 关注
官方认证
内容提供者

权威、专业、丰富

认证主体四川尚阅网络信息科技有限公司
IP属地四川
统一社会信用代码/组织机构代码
91510100MA6716HC2Y

1亿VIP精品文档

相关文档