冉启军090704010070数值分析数本.docVIP

贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数值分析 班级：数本(2)班 实验日期：2011年11月2 学 号：090704010070 姓名：冉啟军 指导教师： 闭海 实验成绩： 一、实验名称 二、实验目的及要求三、实验环境四、实验内容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 法; 3) Romberg 方法. 五、算法描述及实验步骤……n执行 x=x+h;Tn=Tn+2f(x); 步3 Tn=h/2*Tn 步4 输出Tn 结束 (2) 复合Simpson法 功能 用复合Simpson公式公式（3.14）……n执行 x=x+h;Sn=Sn+4f(x);x=x+h/2;Sn=Sn+2f(x); 步3 S=h/6*Sn 步4 输出Sn 结束 (3)Romberg积分法 输入 被积函数f(x)，积分区间端点a,b,允许误差 输出 Romberg积分值 步1 步2 反复执行步3→步9. 步3 步4 反复执行步5→步6. 步5 步6 若x≥b,则退出本层循环. 步7 执行 步8 执行 步9 若e≤且k≥5,则退出循环. 步10 步11 输出 2、实验二算法描述 （1）变步长梯形法 输入 被积函数f(x),积分区间端点a,b和允许误差. 输出 复合梯形积分值 步1 步2 步3 反复执行步4→步10 步4 步5 反复执行步6→步7 步6 步7 若x≥b，则退出本层循环 步8 步9 步10 若e≤,则退出循环 步11 步12 输出 （2）变步长Simpson·法 步1 定义被积函数f(x) 步2 输入积分区间端点a,b和允许误差 步3 开始计时：tic 步4 调用基于变步长Simpson积分函数[S,n]=quad(‘f’,a,b,tol) 步5 结束计时：t2=toc (3)Romberg方法 输入 被积函数f(x),积分区间端点a,b和允许误差 输出 Romberg积分值 步骤 调用Romberg积分程序 六、调试过程及实验结果t=[0 10 20 30 40 50 60 70 80]; a=[30.00 31.63 33.44 35.47 37.75 40.33 42.39 46.69 50.67]; h=10;v0=0; v50=v0+(h/2)*(a(1)+2*(a(2)+a(3)+a(4)+a(5))+a(6)) v80=v0+(h/2)*(a(1)+2*(a(2)+a(3)+a(4)+a(5)+a(6)+a(7)+a(8))+a(9)) p1 v50 = 1.7345e+003 v80 = 3.0803e+003 复合Simpson法： t=[0 10 20 30 40 50 60 70 80]; a=[30.00 31.63 33.44 35.47 37.75 40.33 42.39 46.69 50.67]; h=20;v0=0; v80=v0+(h/6)*(a(1)+4*(a(2)+a(4)+a(6)+a(8))+2*(a(3)+a(5)+a(7))+a(9)) x=[0 10 20 30 40 50];y=[30.00 31.63 33.44 35.47 37.75 40.33]; p=polyfit(x,y,2);poly2sym(p);x=[5 15 25 35 45 ];q=polyval(p,x); a=[30.00 q(1) 31.63 q(2) 33.44 q(3) 35.47 q(4) 37.75 q(5) 40.33]; h=10; v50=v0+(h/6)*(a(1)+4*(a(2)+a(4)+a(6)+a(8)+a(10))+2*(a(3)+a(5)+a(7)+a(9))+a(11)) p2 v80 = 3.0810e+003 v50 = 1.7336e+003 Romberg积分法： function R2n=Romberg(f,a,b,tol) h=b-a;T1=(h/2)*(feval(f,a)+feval(f,b));S1=0;C1=0;R1=0;k=0