复合梯形公式、复合辛普森公式.doc

宁夏师范学院数学与计算机科学学院 《》实验报告 实验序号：　　　　　　　实验项目名称： 学　　号 姓　　名 专业、班 实验地点 指导教师 时间 一、实验目的及要求 二、实验设备（环境）及要求 1、环境要求： 硬件：一般要求486以上的处理器、16MB以上内存、足够的的硬盘可用空间 (随安装组件的多少而定)； 软件：。 实验内容. 1.复合梯形公式： 源程序： function y=fx(x) y=sin(x)./x; function T_n=fht(a,b,n) h=(b-a)/n; for k=0:n x(k+1)=a+k*h; if x(k+1)==0 x(k+1)=10^(-10); end end T_1=h/2*(fx(x(1))+fx(x(n+1))); for i=2:n F(i)=h*fx(x(i)); end T_2=sum(F); T_n=T_1+T_2; 复合辛普森公式： 源程序： function y=f(x) y=sin(x)./x; function S_n=S_P_S(a,b,n) h=(b-a)/n; for k=0:n x(k+1)=a+k*h; x_k(k+1)=x(k+1)+1/2*h; if (x(k+1)==0)|(x_k(k+1)==0) x(k+1)=10^(-10); x_k(k+1)=10^(-10); end end S_1=h/6*(f(x(1))+f(x(n+1))); for i=2:n F_1(i)=h/3*f(x(i)); end for j=1:n F_2(j)=2*h/3*f(x_k(j)); end S_2=sum(F_1)+sum(F_2); S_n=S_1+S_2; 四．实验结果与数据处理 T_1=fht(0,1,2) T_1 = 0.9398 vpa(T_1,7) ans = .9397933 T_2=fht(0,1,4) T_2 = 0.9445 vpa(T_2,7) ans = .9445135 T_3=fht(0,1,8) T_3 = 0.9457 vpa(T_3,7) ans = .9456909 以此类推，得到如下计算结果： k 1 2 3 4 5 Tn 0.9397933 0.9445135 0.9456909 0.9459850 0.9460596 k 6 7 8 9 10 Tn 0.9460769 0.9460815 0.9460827 0.9460830 0.9460381 清单： S_1=S_P_S(0,1,2) S_1 = 0.9495 vpa(S_1,7) ans = .9495483 S_2=S_P_S(0,1,4) S_2 = 0.9465 vpa(S_2,7) ans = .9465170 S_3=S_P_S(0,1,8) S_3 = 0.9461 vpa(S_3,7) ans = .9461373 S_4=S_P_S(0,1,16) S_4 = 0.9461 vpa(S_4,7) ans = .9460899 S_5=S_P_S(0,1,32) S_5 = 0.9461 vpa(S_5,7) ans = .9460839 五、分析与讨论 六、教师评语 成绩