《数值计算方法》.docVIP

华北科技学院基础部综合性实验 实 验 报 告 课程名称 数值计算方法 实验学期 2013 至 2014 学年 第 1 学期 学生所在系部 建筑工程学院 年级 11级 专业班级 测绘B113班 学生姓名 赵静 学号 201105064329 任课教师 李 慧 实验成绩 《 数值计算方法 》课程综合性实验报告 开课实验室：数学应用实验室 年 月 日 实验题目 非线性方程的常微分方程初值问题的数值解—库塔法与阿当姆斯法等。 3）针对习题编制程序，并上机计算其所需要的结果。 4）通过对各种求解方法的计算实习，体会各种解法的功能、优缺点及适用场合，会选取适当的算法。 二、设备与环境 Matlab软件等。 三、实验内容及要求 1. 用二分法和牛顿法求解非线性方程的根，并就收敛性、收敛速度等进行比较分析； 2．常微分方程初值问题的数值解法有：欧拉法、后退欧拉法、梯形法、改进欧拉法、四阶龙格—库塔法与阿当姆斯法。选用至少两种方法解同一个题目，就其收敛性，局部阶段误差，稳定性等进行分析。 3. 题目来源：实验参考书的习题，和教材上的例题或习题。 实验结果及分析 解非线性方程的根 1.二分法 function[c,err,yc]=bisect(fun,a,b,delta) ya=fun(a) yb=fun(b) if yb==0 c=b; return end if ya*yb0 disp(‘（a,b）不是有根区间’); return end max1=1+round((log(b-a)-log(delta))/log(2)); for k=1:max1 c=(a+b)/2; yc=fun(c); if yc==0 a=c; b=c; return elseif yb*yc0 b=c; yb=yc; else a=c; ya=yc; end if(b-a)delta,return,end end k; c=(a+b)/2; err=abs(b-a); yc=fun(f,c); 结果 f=@(x)(x^3-x-1) f = @(x)(x^3-x-1) bisect(f,1.0,1.5,0.0005) ya = -1 yb = 0.8750 ans = 1.3247 2. 牛顿法 function[p1,err,k,y]=newton(fun,dfun,p0,delta,max1) y=f(p1) p0=fun(p0) for k=1:max1 p1=p0-fun(p0)/dfun(p0); err=abs(p1-p0); p0=p1 ; y=fun(p1) if(errdelta)|(y==0), break, end y=f1041(p1) end p1,err,k 结果 f=@(x)(x^3-3*x+2) f = @(x)(x^3-3*x+2) newton(‘f1041’,’df1041’,1.2,10^(-6),18) err = 8.1277e-007 k = 18 y= 1.9817e-007 ans = 1.0000 说明经过18次迭代能得到满足精度的值 常微分方程初值问题的数值解法 1.欧拉法 function[t,y]=ode_euler(f,a,b,ya,M) h=(b-a)/M; T=zeros(1,M+1); Y=zeros(1,M+1); T=a:h:b, Y(1)=ya; for k=1:M Y(k+1)=Y(k)+h*f(T(k),