数值分析5.docVIP

成 绩 11信计2012-2013（二） 《数值分析》 实验报告 设计题目　 Newton迭代法求根 所在班级　 11精算 指导教师　 杨扬 学号姓名　 20110402111 王兰花 　 20110402112 李源 　 20110402113 孙可可 　 20110402114 曹炜炜 　 20110402115 曹云 徐州工程学院数学与物理科学学院 Newton迭代法求根 一、二、实验 2）无求根公式的五次方程 3）超越方程 要求: （1）先用图象法求初始近似,再用Newton迭代法加工,计算结果达到5位有效数字. （2）设计一个程序找出Newton迭代的收敛域 三、实验过程 3.1算法思想 输入 ； 输出 方程在附近的根或失败信息； 步1 ； 步2 对执行步3——步5； 步3 若则，退出循环；否则：； 步4 ； 步5 若则，退出循环； 步6 若则输出，否则若则输出“迭代M次失败”，否则输出“奇异”； 3.2程序求解及结果 x1= 3; x0 = 0; while abs(x1 - x0) 1e-5 f = x1^3-2*x1-5; g = 3*x1^2 - 2; x0 = x1; x1 = x1 - f/g; end x1 x=0:0.001:10; plot(x,x.^3-2.*x-5) grid on x=-2:0.01:2; plot(x,x.^5-4.*x-2) grid on x=-4*pi:0.1:4*pi; ezplot(tan(x),x) hold on ezplot(x,x) ezplot(0,x) grid on solve(x-tan(x)=0,x) ans =0. 四、、附录（源程序清单）function x=Newton(f,df,x0,tol,M) dv=0; for k=1:M if feval(df,x0)==0 dv=2;break else x1=x0-feval(f,x0)/feval(df,x0); end e=abs(x1-x0);x0=x1; if e=tol dv=1;break end end if dv==1 x=x1; elseif dv==0 x=[迭代 ,num2str(M),次失败]; else x=奇异 ; end function y=f(x) y=x^3-2*x-5; function y=df(x) y=3*x^2-2; x0=2.2;tol=0.00005;M=100; x=Newton(f,df,x0,tol,M) x =2.0946 2) function y=f(x) y=x^5-4*x-2; function y=df(x) y=5*x^4-4; x0=0.;tol=0.00005;M=100; x0=0;tol=0.00005;M=100; x=Newton(f,df,x0,tol,M) x =-0.5085 3) function y=f(x) y=x-tan(x); function y=df(x) y=1-1/1+x^2; x0=0.1;tol=0.00005;M=100; x=Newton(f,df,x0,tol,M) x =迭代100次失败 1