数值分析与算法 变步长梯形求积法计算定积分_第1页
数值分析与算法 变步长梯形求积法计算定积分_第2页
数值分析与算法 变步长梯形求积法计算定积分_第3页
数值分析与算法 变步长梯形求积法计算定积分_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数值分析与算法变步长梯形求积法计算定积分数值分析与算法变步长梯形求积法计算定积分数值分析与算法变步长梯形求积法计算定积分数值分析与算法变步长梯形求积法计算定积分编制仅供参考审核批准生效日期地址：电话：传真:邮编:变步长梯形求积法计算定积分1．原理：变步长求积法的主要思想是利用若干小梯形的面积代替原方程的积分，当精度达不到要求时，可以通过增加点数对已有的区间再次划分，达到所需精度时即可；其中由于新的式子中有原来n点中的部分项，故可以省略一些计算，符合了计算机计算存储的思想。主要公式：T2n=Tn/2+(h/2)*Σf(xk+;2．C++语言实现方式：通过每次的Tn值和新增的函数值点计算T2n，再通过判断|Tn-T2n|的大小来判断是否达到精度要求。3．源程序如下：#include""#include""doublef(doublex)//预先输入的待积分函数{ doubles; s=log(x*x); return(s);}doubleffts(doublea,doubleb,doubleeps){ intn,k; doublefa,fb,h,t1,p,s,x,t; fa=f(a); fb=f(b); n=1; h=b-a; t1=h*(fa+fb)/2; p=eps+1; while(p>=eps) { s=0; for(k=0;k<=n-1;k++) { x=a+(k+*h; s=s+f(x); } t=t1/2+h*s/2; p=fabs(t1-t); cout<<"步长n为:"<<n<<"时的"<<"Tn="<<t1<<'\t'<<"T2n="<<t<<'\t'<<"误差变化:"<<p<<endl; t1=t; n=n*2; h=h/2; } return(t);}voidmain(){ doubleresult,a,b,eps; cout<<"需要求解的积分式为f(x)=log(x^2)"<<endl; cout<<"输入边界值a="<<'\t'; cin>>a; cout<<"输入边界值b="<<'\t'; cin>>b; cout<<"输入误差限"<<'\t'; cin>>eps; result=ffts(a,b,eps); cout<<"经过变步长梯形求积法得方程结果为:"<<result<<endl;}4．程序运行结果：根据程序提示依次输入积分上限和积分下限，然后输入误差限；本程序需要预先在程序中输入需要积分方程的表达式。程序运行的结果如下图所示：5．实验小结：通过本次实验，我学会了变步长梯形求积法的实现方式。其本质在于通过增加函数值点来获得精度的提高。而定积分的实质就是面积的计算，这样就可以

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论