复化梯形公式和复化辛普森公式的精度比较

实验四、复化梯形公式和复化Simpson公式的精度比较 （2学时）

一、实验目的与要求

1、熟悉复化Simpson公式和复化梯形公式的构造原理； 2、熟悉并掌握二者的余项表达式；

3、分别求出准确值，复化梯形的近似值，复化Simpson的近似值，并比较后两者的精度；

4、从余项表达式，即误差曲线，来观察二者的精度，看哪个更接近于准确值。 二、实验内容： 对于函数f(x)表格如下：

0 1 1/8 0.9973978 1/4 0.96158 3/8 0.9767267 1/2 0.9588510 5/8 0.9361556 3/4 0.9088516 7/8 0.8771925 1 0.8414709 1sinxsinx，试利用下表计算积分Idx。

0xx注：分别利用复化梯形公式和复化Simpson公式计算，比较哪个精度更好。 其中：积分的准确值I0.9460831。 三、实验步骤

1、 熟悉理论知识，并编写相应的程序；

2、 上机操作，从误差图形上观察误差，并与准确值相比较，看哪个精度更好； 3、 得出结论，并整理实验报告。 四、实验注意事项

1、复化梯形公式，程序主体部分： for n=2:10 T(n)=0.5*T(n-1) for i=1:2^(n-2)

T(n)=T(n)+(sin((2*i-1)/2^(n-1))/((2*i-1)/2^(n-1)))/2^(n-1)； end

1文档来源为:从网络收集整理.word版本可编辑.

文档收集于互联网，已重新整理排版.word版本可编辑.欢迎下载支持.

end

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

五．实验内容

复化梯形公式和复化辛普森公式的引入 复化梯形公式：

hTn[f(xkf(xk1)]；

k02n1复化辛普森公式：

hSn[f(xk4f(x1)f(xk1)]；

kk062n1根据题意和复化梯形公式、复化辛普森公式的原理编辑程序求解代码如下： Matlab代码 clc

s=quad('sin(x)./x',0,1) p1=zeros(10,1);

1文档来源为:从网络收集整理.word版本可编辑.

文档收集于互联网，已重新整理排版.word版本可编辑.欢迎下载支持.

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(i))*(y(i)+y(i+1)); end

for 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; end p1; p2;

s1=s+p1(4) s2=s+p2(4) format long

for k=1:length(p1) p1(k)=abs(p1(k)); p2(k)=abs(p2(k)); end p1 p2

plot(6:1:15,p1,'-r')

1文档来源为:从网络收集整理.word版本可编辑.

文档收集于互联网，已重新整理排版.word版本可编辑.欢迎下载支持.

hold on

plot(6:1:15,10000*(p2),'-c') hold off

部分程序结果输出： s =

0.0076534 s1 =

0.3582701 s2 =

0.5384947 结果分析

1根据结果输出可知：积分Isin(x)dx的准确值为：I= 0.0076534； x0通过复化梯形公式和复化辛普森公式得到的积分值为： s1 =0.3582701: s2 =0.5384947; 相对误差为：

1S1I1004.15104；

100ISI100221.62108；

100I显然，从相对误差可知通过辛普森公式得到的结果误差小精度高。

由于以上的算法只算了结点个数为9的情况，只能横向比较两公式的精确程度，而不能分别比较两公式随节点个数变化精度的变化，故而将以上程序重新编（以上程序为最终程序）可得出两公式随节点个数变化精度的变化情况所取得节点个数为从6到15，共计10种情况对应的误差值如下表：

节点数 T 6 0.001 7 8 9 10 0.0007 0.00051 0.00039 0.00031 S*10000 0.98762 0.47766 0.25913 0.15308 0.09666 节点数 11 12 13 14 15 1文档来源为:从网络收集整理.word版本可编辑.

文档收集于互联网，已重新整理排版.word版本可编辑.欢迎下载支持.

T 0.00025 0.00021 0.00017 0.00015 0.00013 S*10000 0.041 0.04491 0.03257 0.02444 0.012 （表1）

注：由于辛普森公式的精度较高，所得的误差值较小不宜比较，故而将辛普森公式计算出的误差值乘上10000得到以上表1的结果，其相应的曲线图如下（图1）。

（图1:两误差曲线比较）

备注：红色，青色分别比奥斯曲线复化梯形公式和复化辛普森公式的误差值曲线

从曲线图可知复化梯形公式和复化辛普森公式的随着节点个数的增加误差值越小即其精度逐渐增大，且复化辛普森公式的精度远高于复化梯形公式的精度。

1文档来源为:从网络收集整理.word版本可编辑.