收敛性与收敛速度的比较实验目的:通过用不同迭代法解同一非线性方程,比较各种方法的收敛性与收敛速度。实验内容:求解非线性方程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方法直接求根收敛速度较快,迭代次数较少。