牛顿法和割线法.doc_第1页
牛顿法和割线法.doc_第2页
牛顿法和割线法.doc_第3页
牛顿法和割线法.doc_第4页
牛顿法和割线法.doc_第5页
已阅读5页，还剩6页未读， 继续免费阅读

下载本文档

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

文档简介

作业十（第五章）：1. 在区间(0,1.5)上分别用二分法、牛顿法和割线法编程求下面的函数的零点，精度要求10-10。二分法function X=bisection(fx,xa,xb,n,delta)% 二分法解方程% fx是由方程转化的关于x的函数，有fx=0。% xa 解区间上限% xb 解区间下限%解区间人为判断输入% n 最多循环步数，防止死循环。%delta 为允许误差x=xa;fa=eval(fx);x=xb;fb=eval(fx);for i=1:n xc=(xa+xb)/2;x=xc;fc=eval(fx); X=i,xc,fc; if fc*fa0 xb=xc; else xa=xc; end if (xb-xa)delta,break,endEnd二分法结果：迭代34次，xc=0.5149牛顿法function X=newton(fx,e,x0,m)x=x0;k=0;F=eval(fx);if abs(F)=e X=x F; disp(X); return endwhile k=m x=x0;g=eval(diff(fx); x1=x0-F/g; x=x1;F=eval(fx);k=k+1; if abs(F)m fprintf(牛顿法迭代M次没有找到方程的根) return end x0=x1;endfprintf(n%s%.4ft%s%d,X=,X,k=,k) %输出结果牛顿法结果：迭代5次 结果0.5149割线法：function X=gx9(fx,x0,x1,m,e)x=x0;f0=eval(fx);x=x1;f1=eval(fx);if abs(f0)=e X=x0,f0;endfor k=2:m if abs(f0)abs(f1) b=x0;x0=x1;x1=b; b=f0;f0=f1;f1=b; end t=(x1-x0)/(f1-f0);x0=x1;f0=f1; x1=x1-t*f1; x=x1;f1=eval(fx); if abs(f1)0 y0=y1;continue end x0=x-k;x1=x; X=gx9(fx,x0,x1,m,e);%割线法 if xb X=x;disp(X); return end y0=y1;end%function X=gx9(fx,x0,x1,m,e)x=x0;f0=eval(fx);x=x1;f1=eval(fx);if abs(f0)=e X=x0,f0;endfor k=2:m if abs(f0)abs(f1) b=x0;x0=x1;x1=b; b=f0;f0=f1;f1=b; end t=(x1-x0)/(f1-f0);x0=x1;f0=f1; x1=x1-t*f1;
人人文库> 全部分类> 应用文书 > 产品手册

温馨提示

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

最新文档

评论

0/150

提交评论