数值分析实验三.docVIP

贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数值分析 班级： 12级信本 实验日期： 2014年 11月5日 学 号： 120703010031 姓名： 李盼 指导教师： 实验成绩： 一、实验名称 二、实验目的及要求三、实验环境四、实验内容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 法. 五、算法描述及实验步骤。 复合Simpson 输入：被积函数，积分区间[a,b]和n, 输出：复合Simpson积分值。 步1 :; ;. 步2：对k=1,2,…,n执行 步3：。 步4：输出。 Romberg公式 根据已知数据对其进行多项式拟合得出p(x);输入被积函数，积分区间端点a,b,允许误差输出Romberg积分值。 步1 :; ;. 步2：对k=1,2,…,n执行 步3：。 步4：输出。 六、调试过程及实验结果 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)) v50 = 1.7345e+003 v80=v0+(h/2)*(a(1)+2*(a(2)+a(3)+a(4)+a(5)+a(6)+a(7)+a(8))+a(9)) v80 = 3.0803e+003 复合Simpson 法： 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)) v80 = 3.0810e+003 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)) 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;? while?1 S=0; x=a+h/2; while?1 S=S+feval(f,x); x=x+h; if?x=b bread end end T2=T1/2+(h/2)*S; S2=(4/3)*T2-(1/3)*T1; C2=(16/15)*S2-(1/15)*S1; R2=(64/63)*C2-(1/63)*C1; e=abs(R2-R1); h=h/2; T1=T2;S1=S2;C1=C2;R1=R2;k=k+1;??if?e=tolk=5?????break????endendR2n=R2; function?y=f(x)y=(5930286613325939*x^3)/295147905179352825856?-?(6590784751299553*x^2)/9223372036854775808?+?(3387527922312091*x)/18014398509481984?+?74033/2475?x=[0:10:80]x?=?0?10 20 30 40?50?60?70?80?y=[30.00,31.63,33.44,35.47,37.