首页 收敛性与收敛速度的比较

收敛性与收敛速度的比较

举报
开通vip

收敛性与收敛速度的比较收敛性与收敛速度的比较实验目的:通过用不同迭代法解同一非线性方程,比较各种方法的收敛性与收敛速度。实验内容:求解非线性方程X2-3x+2-ex二0的根,准确到10-6。实验要求:用你自己设计的一种线性收敛的迭代法求方程的根(利用算法4.1(简单迭代法)计算);用斯蒂芬森加速迭代(算法4.2)计算。输出迭代初值、各次迭代值及迭代次数用牛顿法(算法4.3)求方程的根,输出迭代初值、各次迭代值及迭代次数,并与(1)、(2)的结果比较。用MATLAB内部函数solve直接求出方程的所有根,并与(1)、(2)、(3)的结果进...

收敛性与收敛速度的比较
收敛性与收敛速度的比较实验目的:通过用不同迭代法解同一非线性方程,比较各种方法的收敛性与收敛速度。实验内容:求解非线性方程X2-3x+2-ex二0的根,准确到10-6。实验 要求 对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗 :用你自己设计的一种线性收敛的迭代法求方程的根(利用算法4.1(简单迭代法)计算);用斯蒂芬森加速迭代(算法4.2)计算。输出迭代初值、各次迭代值及迭代次数用牛顿法(算法4.3)求方程的根,输出迭代初值、各次迭代值及迭代次数,并与(1)、(2)的结果比较。用MATLAB内部函数solve直接求出方程的所有根,并与(1)、(2)、(3)的结果进行比较。实验程序:(1)简单迭代法:初值选用0.5formatlongc=10「(-6);x0=0.5;symsxfxfx=(x."2-exp(x)+2)/3;j=1;fx0=(x0."2-exp(x0)+2)/3while(abs((x0-fx0)/fx0)>c)x0=fx0;fx0=(x0."2-exp(x0)+2)/3j=j+1;endj结果:fx0=0.200426243099957fx0=0.272749065098375fx0=0.253607156584130fx0=0.258550376264936fx0=0.257265636335094fx0=0.257598985162190fx0=0.257512454514832fx0=0.257534913615251fx0=0.257529084167956fx0=0.257530597238330fx0=0.257530204510457fx0=0.257530306445639j=12j=12迭代次数:12(2)斯蒂芬森加速迭代法:formatlongc=10.人(-6);y=0.5;symsxfxgxhxfx=(x42-exp(x)+2)/3;gx=(fx42-exp(fx)+2)/3;hx=x-(fx-x).人2/(gx-2*fx+x);j=1;fy=(y.A2-exp(y)+2)/3;gy=(fy.A2-exp(fy)+2)/3;hy=y-(fy-y).A2/(gy-2*fy+y)while((abs(y-hy)/hy)>c)y=hy;fy=(y.A2-exp(y)+2)/3;gy=(fy.A2-exp(fy)+2)/3;hy=y-(fy-y).A2/(gy-2*fy+y)j=j+1;endj结果:hy=hy=0.257530317719808hy=0.257530285439861迭代次数:3(3)牛顿法formatlongc=10.人(-6);x0=0.5;symsxfxfx=x-(x42-3*x-exp(x))+2/(2*x-exp(x)-3);j=1;fx0=(x0.A2-exp(x0)+2)/3while((abs(x0-fx0)/fx0)>c)x0=fx0;fx0=x0-(x0.A2-3*x0-exp(x0)+2)/(2*x0-exp(x0)-3)j=j+1;end结果:fx0=0.200426243099957fx0=fx0=0.257530275750830fx0=0.257530285439861迭代次数:4(4)solve函数直接求根x=solve('xA2-3*x+2-exp(x)=0')x=0.25753028543986076045536730493724用solve方法直接求根收敛速度较快,迭代次数较少。
本文档为【收敛性与收敛速度的比较】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
is_769254
暂无简介~
格式:doc
大小:10KB
软件:Word
页数:5
分类:
上传时间:2019-11-18
浏览量:20