实验报告二 一元非线性方程解法...docVIP

浙江大学城市学院实验报告 课程名称 科学计算 实验项目名称 一元非线性方程的解法 实验成绩 指导老师（签名 ） 日期 2013/10/10 一. 实验目的和要求 用Matlab软件掌握求解非线性方程的二分法、迭代法和牛顿法，并对结果作初步分析； 通过实例练习用非线性方程求解实际问题。 二. 实验内容和原理 分析应用题2-1，2-2，2-3，2-4，2-5要求将问题的分析过程、Matlab源程序和运行结果和结果的解释、算法的分析写在实验报告上。 程序注释 对用二分法、迭代法和牛顿法求一元非线性方程数值解的Matlab程序的每一句添上注释。 二分法： [x,n]=BisectionMethod(a,b,epsi,Nmax) 迭代法： [x,n]=IterationMethod(x0,epsi,Nmax) Newton法： [x,n]=NewtonMethod(x0,epsi,Nmax) 分析应用题 用二分法求方程和根的近似值，使误差不超过，输出每次二分之后解的近似结果以及二分的次数，其中的根只需求最接近原点的那个根。 分析应用题 已知方程有以下三种迭代格式，分析三种迭代格式的收敛性，求出迭代精度为的数值结果，并比较迭代序列的收敛速度。 1） 2） 3） 分析应用题 用下列方法求方程的近似根，要求误差不超过，并比较计算量。 1）在区间[0, 1]上用二分法； 2）取初值，并用迭代过程； 3）取初值用牛顿法求解。 分析应用题 以定期存储为基础的储蓄帐户的累积值可由定期年金方程确定 在这个方程中，是帐户中的数额，是定期存储的数额，是个存储期间的每期利率。一个工程师想在20年后退休时储蓄帐户上的数额达到750000美元，而为了达到这个目标，他每个月能存1500美元。为实现他的储值目标，最小利率应是多少？假定利息是月复利的。 三. 操作方法与实验步骤（包括实验数据记录和处理） 分析应用题 对用二分法、迭代法和牛顿法求一元非线性方程数值解的Matlab程序的每一句添上注释。 二分法： [x,n]=BisectionMethod(a,b,epsi,Nmax) Function[c,err,yc]=erfen(f,a,b,delta) Ya=feval(f,a); Yb=feval(f,b); If ya*yb0,break,end Maxl=1+round((log(b-a)-log(delta))/log(2)); For k=1:max1 C=(a+b)/2; Yc=feval(f,c); ?If yc = 0 a=c; b=c; elseif yb*yc0 b=c; yb=yc; else a= c; ya=yb; end if b-adetal,break,end end c=(a+b)/2; err= abs(b-a); yc = feval(f,c); 输入：a,b值及精确控制量ε； If f(a)*f(b)0 then 返回第一步，重新输入a,b值else转第三步； While|a-b|ε 时做 x=1/2（a+b）,计算f（x）；if f（x）=0 then 输出x，停机。 if f（a）*f（x）0 then [a,b]替换为[a,x] else [a,b]替换为[x,b] endwhile; 输出x=1/2(a+b). 迭代法： [x,n]=IterationMethod(x0,epsi,Nmax) 设f(x)=0,改写为 X=ψ（x）Xn+1=ψ（xn）（n=1,2,3……） 若产生的序列{xn}收敛于x*，即x*即为方程的根。 Newton法： [x,n]=NewtonMethod(x0,epsi,Nmax) Function[p0,err,k,y]=Newton(f,df,p0,delta,epsilon,max1) For k=1:mak1 P1=p0-feval(f,p0)/feval(df,p0); Err=abs(p1-p0); Relerr =2*err/(abs(p1)+delta); Pp(k)=p1; P0=p1; Y=feval(f,p0); If(errdelta)|(relerrdelta)|(abs(y)epsilon),break,end End Pp （1） 输入：x0，ε； （2） f0=f（x0），f（1）=f’(x0);