2022年复化梯形公式和复化辛普森公式的精度比较 .pdf

《2022年复化梯形公式和复化辛普森公式的精度比较 .pdf》由会员分享，可在线阅读，更多相关《2022年复化梯形公式和复化辛普森公式的精度比较 .pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、实验四、复化梯形公式和复化Simpson 公式的精度比较（2 学时）一、实验目的与要求1、熟悉复化 Simpson公式和复化梯形公式的构造原理；2、熟悉并掌握二者的余项表达式；3、分别求出准确值，复化梯形的近似值，复化Simpson 的近似值，并比较后两者的精度；4、从余项表达式，即误差曲线，来观察二者的精度，看哪个更接近于准确值。二、实验内容：对于函数sin( )xf xx，试利用下表计算积分10sin xIdxx。表格如下：x0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 ( )f x1 0.9973978 0.9896158 0.9767267 0.9588510 0.9

2、361556 0.9088516 0.8771925 0.8414709 注：分别利用复化梯形公式和复化Simpson公式计算，比较哪个精度更好。其中：积分的准确值0.9460831I。三、实验步骤1、熟悉理论知识，并编写相应的程序；2、上机操作，从误差图形上观察误差， 并与准确值相比较， 看哪个精度更好；3、得出结论，并整理实验报告。四、实验注意事项精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 6 页1、复化梯形公式，程序主体部分：for n=2:10 T(n)=0.5*T(n-1) for i=1:2(n-2) T(n)=T(n)

3、+(sin(2*i-1)/2(n-1)/(2*i-1)/2(n-1)/2(n-1)；endend 2、复化 Simpson公式，程序主体部分：for i=1:10 n=2.i x=0:1/n:1 f=sin(x)./x f(1)=1 s=0 for j=1:n/2 s=s+f(2*j) end t=0 for j=1:(n/2-1) t=t+f(2*j-1) end S(i)=1/3/n*(f(1)+4*s+2*t+f(n+1) end五实验内容复化梯形公式和复化辛普森公式的引入精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 6 页复化

4、梯形公式：110()2nnkkkhTf xf x；复化辛普森公式：11102(4 ()()6nnkkkkhSf xf xf x；根据题意和复化梯形公式、复化辛普森公式的原理编辑程序求解代码如下：Matlab 代码clc s=quad(sin(x)./x,0,1) p1=zeros(10,1); p2=zeros(10,1); for k=6:15 s1=0; s2=0; x=linspace(0,1,k); y=sin(x)./x; z=(1/(2*(k-1):(1/(k-1):1; sz=sin(z)./z; y(1)=1; for i=1:(k-1) s1=s1+0.5*(x(i+1)-x

5、(i)*(y(i)+y(i+1); endfor j=1:(k-1) s2=s2+(1/6)*(x(j+1)-x(j)*(y(j)+y(j+1)+4*sz(j); end p1(k-5)=s1-s; p2(k-5)=s2-s; endp1; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 6 页p2; s1=s+p1(4) s2=s+p2(4) format longfor k=1:length(p1) p1(k)=abs(p1(k); p2(k)=abs(p2(k); endp1 p2 plot(6:1:15,p1,-r) hold

6、onplot(6:1:15,10000*(p2),-c) hold off部分程序结果输出：s = 0.946083070076534 s1 = 0.945690863582701 s2 = 0.946083085384947 结果分析根据结果输出可知： 积分10sin( )xIdxx的准确值为： I= 0.946083070076534；精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 6 页通过复化梯形公式和复化辛普森公式得到的积分值为：s1 =0.945690863582701: s2 =0.946083085384947; 相对误

7、差为：4111004.15 10100SII；8221001.62 10100SII；显然，从相对误差可知通过辛普森公式得到的结果误差小精度高。由于以上的算法只算了结点个数为9的情况，只能横向比较两公式的精确程度，而不能分别比较两公式随节点个数变化精度的变化，故而将以上程序重新编 （以上程序为最终程序） 可得出两公式随节点个数变化精度的变化情况所取得节点个数为从 6到15 ，共计 10种情况对应的误差值如下表：节点数6 7 8 9 10 T 0.001 0.0007 0.00051 0.00039 0.00031 S*10000 0.98762 0.47766 0.25913 0.15308

8、0.09666 节点数11 12 13 14 15 T 0.00025 0.00021 0.00017 0.00015 0.00013 S*10000 0.06441 0.04491 0.03257 0.02444 0.01892 （表1）注：由于辛普森公式的精度较高，所得的误差值较小不宜比较，故而将辛普森公式计算出的误差值乘上10000得到以上表 1 的结果，其相应的曲线图如下（图 1） 。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 6 页（图1: 两误差曲线比较）备注：红色，青色分别比奥斯曲线复化梯形公式和复化辛普森公式的误差值曲线从曲线图可知复化梯形公式和复化辛普森公式的随着节点个数的增加误差值越小即其精度逐渐增大，且复化辛普森公式的精度远高于复化梯形公式的精度。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 6 页