实验三实验报告2.docVIP

贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数 值 分 析 班级： 数本（1）班 实验日期： 2013年1月9日 学 号：110701010039 姓名： 欧 习 昌 指导教师： 杨 一 都 老 师 实验成绩： 一、实验名称 二、实验目的及要求三、实验环境四、实验内容0 10 20 30 40 50 60 70 80 a(m/s) 30.00 31.63 33.44 35.47 37.75 40.33 42.39 46.69 50.67 要求:分别用复化梯形法,复化Simpson 法和 Romberg公式计算. 题2 给定积分 和 ,分别用下列方法计算积分值要求准确到 ,并比较分析计算时间. 变步长梯形法; 变步长 Simpson 法 Romberg 方法 算法描述及实验步骤[30.00,31.63,33.44,35.47,37.75,40.33,42.39,46.69,50.67] 步1：输入h 步2：T=(h/2)*(a(1)+a(6)+2*(a(2)+a(3)+a(4)+a(5))) 步3：T=(h/2)*(a(1)+2*(a(2)+a(3)+a(4)+a(5)+a(6)+a(7))+a(8)) (2)根据复合simpson公式：Sn= 根据已知输入数据a=[30.00,31.63,33.44,35.47,37.75,40.33,42.39,46.69,50.67] 先用Lagrange插值得出t=5,15,25,35,45,55,65,75处的值 步1：h=;Sn=f(a)-f(b);x=a. 步2：对k=1,2,...,n执行 x=x+h/2;Sn=Sn+4f(x);x=x+h/2;Sn=Sn+2f(x). 步3：Sn=Sn 步4：输出Sn. (3)根据Romberg积分公式Rn= 2.(1). 求积分 ，允许误差为ε。 被积函数f（x）,a,b,ε. 复合梯形积分值T2n. h=b-a. T1=h/2(f(a)+f(b)). 反复执行步4—步10。 S=0;x=a+h/2. 反复执行步6—步7. S=S+f（x）;x=x+h. 若x≥b,则退出本层循环。 T2=T1/2+h/2S e=|T2-T1|;h=h/2;T1=T2. 若e≤ε则退出循环。 T2n=T2. 输出T2n. (2). 求积分 ，允许误差为ε。 被积函数f（x）,a,b,ε. 复合梯形积分值T2n. h=b-a. T1=h/2(f(a)+f(b));S1=0;k=0. 反复执行步3—步9. S=0;x=a+h/2. 反复执行步5—步6. S=S+f（x）;x=x+h. 若x≥b,则退出本层循环。 T2=T1/2+h/2S;S2=4/3T2-1/3T1. e=|S2-S1|;h=h/2;T1=T2;S1=S2;k=k+1 若e≤ε则退出循环。 S2n=S2. 输出S2n. (3).求积分 ，允许误差为ε。 被积函数f（x）,积分区间端点a,b,ε. 复合梯形积分值R2n. h=b-a. T1=h/2(f(a)+f(b));S1=0;C1=0;R1=0;k=0. 反复执行步3—步9. S=0;x=a+h/2. 反复执行步5—步6. S=S+f（x）;x=x+h. 若x≥b,则退出本层循环. 执行 T2=T1/2+h/2S；S2=4/3T2-1/3T1;C2=16/15S2-1/15S1;R2=64/63C2-1/63C1. 执行 e=|R2-R1|;h=h/2;T1=T2.T1=T2;S1=S2;C1=C2;R1=R2;k=k+1. 若e≤ε且k≥5则退出循环。 R2n=R2. 输出R2n. 六、调试过程及实验结果 复合Simpson法 题2 变步长梯形法 总结 八、附录（源程序清单）function cy=Lagrange(x,y,n,cx) m=length(cx);cy=zeros(1,m); for k=1:n+1 t=ones(1,m); for j=1:n+1 if j~=k t=t.*(cx-x(j))./(x(k)-x(j)); end end cy=cy+y(k).*t; end x=0:10:80; y=[30.00,31.63,33.44,35.47,37.75,40.33,42.39,46.69,50.67]; cx=5:10:75; cy=Lagrange(x,y,5,cx