牛顿迭代法解方程及根.docVIP

北 方 民 族 大 学 学生上机报告 课程名称 数值计算方法 上机名称 牛顿迭代法求方程的根 院系 数学与信息科学学院 专业 数学与应用数学 班级 12级应数2班 第 组 组长 周海龙 同组成员 周海龙 　 上机日期 2014.6 　　　　　　　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　的全部根，要求绝对误差限小于(先确定含根区间，然后构造迭代公式进行求解) 算法描述（流程图） 上机程序 function [k,x,wuca,yx]=newton(x0,tol) k=1; x1=x0-fun(x0)/fun1(x0); disp(x=) while abs(x1-x0)tol x0=x1; k=k+1; x1=x0-fun(x0)/fun1(x0); disp(x1) end k; x=x1; wuca=abs(x1-x0)/2; yx=fun(x); end %分程序1： function y1=fun(x) y1=x^3+3*x^2-1; end %分程序2： function y2=fun1(x) %函数fun(x)的导数 y2=3*x^2+6*x; end %调用函数输入数据 x0=1.0; xr=-2.5; xp=-0.5; tol=0.5*10^-8; [k,x,wuca,yx]=newton(x0,tol) [k,x,wuca,yx]=newton(xr,tol) [k,x,wuca,yx]=newton(xp,tol) 运行结果 x= 0.5486 0.5324 0.5321 0.5321 0.5321 k = 6 x = 0.5321 wuca = 5.9952e-015 yx = 2.2204e-016 x= -2.9009 -2.8797 -2.8794 -2.8794 -2.8794 k = 6 x = -2.8794 wuca = 2.6645e-015 yx = -3.5527e-015 x= -0.6528 -0.6527 -0.6527 k = 4 x = -0.6527 wuca = 1.0850e-009 yx = -1.1102e-016 组号 6 结束 输出迭代失败标志 K=N? ∣x1-xo∣e? X0-f(x0)/f’(x0)=x1 f’(x0)=0? 是 否 1=k 输入x0，e，N, 开始 输出x1 输出奇异标志 k+1=k x1=x0