第二次实验报告(非线性方程求根).doc_第1页
第二次实验报告(非线性方程求根).doc_第2页
第二次实验报告(非线性方程求根).doc_第3页
第二次实验报告(非线性方程求根).doc_第4页
全文预览已结束

下载本文档

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

文档简介

数值计算方法Matlab实验报告 数 学 实 验 报 告实验五学院:数学与统计学院 班级:数学与统计学院(2)班 姓名:石紫雲学号:201470010234非线性方程的数值解法实验1实验目的1)进一步熟练掌握求解非线性方程的牛顿迭代法和弦截法。2)根据牛顿迭代法和弦截法的原理,编写程序求解非线性方程,提高编程解决问题的能力。2 实验内容(1)用牛顿法和割线法求下列方程的根x2-ex=0;x*ex-1=0;(2)用牛顿迭代法求平方根,计算3实验原理(1)牛顿迭代公式:双点弦法公式:(2)令,再用牛顿法求根。4实验步骤1) 根据牛顿迭代法,双点弦法的算法编写相应的求根函数;2) 用牛顿迭代法和双点弦法分别对方程进行求解;5 程序设计牛顿迭代法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(Newton method failed) break end if(abs(x1-x0)delta | abs(fc1(x1)=delta) c=x1; x1=cutnext(x0,x1);x0=c; 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=1*exp(x)+x*exp(x);end3)原函数function fc1=fc1(x)fc1=x2-5;end导函数function fc2=fc2(x)fc2=2*x;end6实验结果及分析方程x2-ex=0x*ex-1=0牛顿迭代法结果牛顿法迭代次数弦截法结果弦截法迭代次数注:牛顿迭代法由于设置delta=1e-6,

温馨提示

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

评论

0/150

提交评论