计算方法第二章方程求根上机报告.doc_第1页
计算方法第二章方程求根上机报告.doc_第2页
计算方法第二章方程求根上机报告.doc_第3页
计算方法第二章方程求根上机报告.doc_第4页
全文预览已结束

下载本文档

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

文档简介

计算方法与实习实验报告 实验报告名称 班级： 学号： 姓名： 成绩：1实验目的1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点。2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法。2 实验内容用牛顿法和割线法求下列方程的根x2-ex=0;x*ex-1=0;lgx+x-2=0;3实验步骤1) 根据二分法和牛顿迭代法，割线法的算法编写相应的求根函数；2) 将题中所给参数带入二分法函数，确定大致区间；3) 用牛顿迭代法和割线法分别对方程进行求解；3 程序设计牛顿迭代法x0=1.0;N=100;k=0;eps=5e-6;delta=1e-6;while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if kN disp(Newmethod failed) break end if(abs(x1-x0)delta | abs(fc1(x1)=delta) c=x1; x1=cutnext(x0,x1); x0=c; %x0 x1x1 x2 x0 x1 end k=k+1; if kN disp(Cutline method failed) break; end if(abs(x1-x0)delta | abs(fc1(x1)delta) break; end endfprintf(%.10fn,x1); function y=cutnext(a,b)y=b-fc(b)/(fc(b)-fc(a)*(b-a);1） 原函数function fc1=fc1(x)fc1=x2-exp(x);end导函数function fc2=fc2(x)fc2=2*x-exp(x);end2）原函数function fc1=fc1(x)fc1=x*exp(x)-1;end导数function fc2=fc2(x)fc2=(x+1)*exp(x);end3）原函数function fc1=fc1(x)fc1=log10(x)+x-2;end导函数function fc2=fc2(x)fc2=1/x/log(10)+1;end 4实验结果及分析1）牛顿法结果-0.7034722378割线法结果-0.70346742252）牛顿法结果0.5671435302割线法结果0.56714329043）牛顿法结果1.7553985566割线法结果1.7555794993牛顿迭代法由于设置delta=1e-6，所以算出的误差e1.0*10-5;割线法由于设置delta=5e-8，所以误差e1.0*10-7;5总结编程时由于将迭代的代码x0=x1放在 if(abs(x1-x0)delta | abs(fc1(x1)delta) break; end之前导致程序没有执行就跳出，通过Debug发现了问题，将x0=x
人人文库> 全部分类> 应用文书 > 作业报告

温馨提示

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

最新文档

评论

0/150

提交评论