龙贝格算法的matlab实现知识讲解_第1页
龙贝格算法的matlab实现知识讲解_第2页
龙贝格算法的matlab实现知识讲解_第3页
全文预览已结束

下载本文档

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

文档简介

1、龙 贝 格 算 法 的 ma t l a b实现作业三龙贝格算法的matlab实现程序流程图：输出“所求次数不够或不可程序源代码：文件 f.m function fx = f(x) if x = 0fx = 1;elsefx = sin(x) / x; end end文件 longbeige.mclc clear all; format longa=input('请输入你要求得积分的下限：')； b=input('请输入你要求得积分的上限:'); e=in put('请输入你要求得积分的结束精度:'); k=i nput('请输入你要求得

2、积分的最大次数:'); fx=(x)sin(x)/x;lbg(f,a,b,k,e)文件 lbg.mfunction lbg(fx,a,b,k,e) T=zeros(k,k);T(1,1)=(b-a)*(1+fx(b)/2;for i=1:k m=0;for j=1:2A(i-1) m=m+fx(a+(2*j-1)*(b-a)/(29);endT(i+1,1)=0.5*T(i,1)+(b-a)*m/2Ai;for n=1:iT(i+1,n+1)=(4An*T(i+1,n)-T(i,n)/(4An-1); endif abs(T(i+1,i+1)-T(i,i)<=e & i&

3、gt;=4 break;elseJend endfor i=1:kif T(i,1)=0 j=i； break;elseendendif j=kerror('所求次数不够或不可积')else endT=T(1:j-1,1:j-1) disp('所求的积分值为：') disp(T(j-1,1)鼻 C&Timand Window曰。 gdit De hug DeskTop Window Help” New to MATLAB?讥咕tzh this亡6 ice Mmosi or read Getting Started.済输入邯辱审獰和甘的下龌沖诲辎入冊耍未律积対的上眼.】倩辎丸惭宴玉骨机分區结束轴虞：J.请辅扎怖姜朮潯视廿菸最光宓埶:"00. 920 北沾241)39肝0L9397933e48OSr?0.914513521665300QL 9490903582701L9i6PSEC20e34»6所求的稅芬徨为乂0.0. 945D36&3j5517'JI0,945DS3J10E884:0.P150?3OS5；8491700C.&1603300i0e36_ld.&4SOB3Q&935Q917C.&160B30703E13790000

人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

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

最新文档

评论

0/150

提交评论